[Geobd] PostgreSQL/Postgis 9.6

Message posté par : Léandre Béron ---------------------------------------- Bonjour, Je cherche à réaliser des statistiques de population à partir des données insee sur nos emprises de quartiers de notre ville. A ma disposition, 3 données : - Les quartiers de la ville - Les emprises iris - les points RIL avec le nombre de logements Je vais donc calculer le nombre de logement par quartier de la ville et ensuite multiplier par la valeur de population d'un logement en fonction du quartier iris. Pour cela, je calcule la valeur de population pour 1 logement, par quartier iris. Juste que là tout va bien. ----------------- Citation : SELECT i.nom_iris, sum(r.nombre_log) AS nb_log_iris FROM insee.ril_millesime r, insee.iris_pop17 i WHERE (st_intersects(st_transform(i.geom,3950), st_transform(r.geom,3950))) AND ((r.categorie) LIKE 'HABIT'::text) GROUP BY i.nom_iris ----------------- Ensuite, je dois connaître le nombre de logement au sein de mes quartiers de la ville pour ensuite multiplier par la valeur pour 1 logement de ma population. Toutefois, mes emprises de mes quartiers chevauchent plusieurs emprises des IRIS. Je dois donc éclater mes emprises de quartier de ville en fonction du chevauchement de mes emprises IRIS.. arrivez-vous à me suivre jusque là ? Cela doit me permettre de sortir : QUARTIER VILLE : A, nombre logement : X, iris correspondant : iris1 QUARTIER VILLE : A, nombre logement : Y, iris correspondant : iris 2 Ainsi de suite.... C'est là que je coince. Je n'arrive pas à éclater mes emprises pour obtenir ce rendu. L'outil de géotraitement UNION sous qgis fait très bien le boulot, mais sous Postgis, je n'arrive pas à reproduire ce traitement. le st_union de postgis ne fait pas du tout la même chose (sauf si j'ai mal compris sa description...) J'ai donc tenté la requête suivante : ----------------- Citation : SELECT s.oid, st_collectionextract(st_intersection(i.geom, s.geom), 3) AS inter, s."Nom", i.nom_iris--, i.pop FROM decoupage_administratif.comite_quartier s, insee.iris_pop17 i ORDER BY s."Nom", i.nom_iris ----------------- Cela me sort plus de 1508 objets, alors que le traitement QGIS en sort environ 200... Après coup, je pourrais lancer mon produit à l'aide des valeurs à l'échelle d'un logement de ma première requête. Quelle fonction géométrique permet de faire cela ? Cordialement, Léandre ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=341793#p341793 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net