Le portail francophone de la géomatique

Il me semble que c'est une question qui irait mieux dans le forum MapInfo (je me trompe ?)


Pour davantage de lisibilité du code, je rajouterai quelques espaces.

Select Parcellaire_A.Commune, Parcellaire_A.Section, Parcellaire_A.Parcelle, 
Parcellaire_B.Parcelle, Parcellaire_B.Surface, 
proportionoverlap(Parcellaire_A.obj, Parcellaire_B.obj) 
from Parcellaire_A, Parcellaire_B 
where Parcellaire_A.obj intersects Parcellaire_B.obj

Et désolé, je connais mal la fonction proportionoverlap.

Dernière modification par ppluvinet (Tue 09 November 2010 08:31)

Bonjour geobrie,

si je comprends bien ta recherche, tu veux renseigner un champ à partir d'une valeur de référence, et au prorata de l'intersection entre deux objets.

Si c'est bien cela, voici un exemple fonctionnel avec la table commune de la bd_topo et une table contenant les perimtres d'inventaire et réglementaires environnementaux.

Je calcule la part de la population communale (valeur fixe popul) potentiellement concernée par la ZPS.

select a.code_insee, a.popul, 
round(((area2d(intersection(a.geometrie, b.geometrie))/area2d(a.geometrie))*100)::numeric,2) AS  "%_de la commune_en_en_zps", 
round(((area2d(intersection(a.geometrie, b.geometrie))/area2d(a.geometrie))*popul)::numeric,2) AS popul_concernee
FROM ign_bd_topo.commune a, perimetre b
WHERE type_perimetre = 'zps'AND st_intersects(a.geometrie, b.geometrie)

[edit]
Si on vaut la population non concernée, on selectionne

round(((1-area2d(intersection(a.geometrie, b.geometrie))/area2d(a.geometrie))*popul)::numeric,2) AS popul_non_concernee

[/edit]
J'espère que ça t'aidera.

Dernière modification par Mathieu BOSSAERT (Tue 23 November 2010 11:41)