Message posté par : Sylvain M.
----------------------------------------
Je vois mieux la problématique.
La solution proposée par Nicolas est la plus propre : il faut calculer la surface de la
zone d'intersection des bâtis avec le RFP, via quelque chose comme :
ST_Area(ST_Intersection(bati.geom,rfp.geom))
(à diviser par la surface d'origine du bâti, pour connaitre le pourcentage dans la
RFP).
Mais dans ce cas, ce qui marche souvent est de se baser non pas sur la géométrie des
bâtiments, mais uniquement sur leur centroide, via ST_Centroid(bati.geom), ou
ST_PointOnSurface(bati.geom) (et donc =>
ST_Intersects(ST_PointOnSurface(bati.geom),rfp.geom) )
En faisant l'intersection dessus, ça marche dans 99.9% des géométries, et c'est
plus rapide que la découpe des géométries intersectées (si le volume de données est
important par exemple).
A+
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=359881#p359881
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