Le portail francophone de la géomatique

Bonjour,

select st_area(st_intersection(a.geom, b.geom))
from parcelle a, ebc b
where st_intersects(a.geom, b.geom)

Creer des index spatiaux sur les tables avant.
Si volumetrie importante, ca peut prendre du temps !

(http://postgis.refractions.net/document … ction.html)

Nicolas

Certaines données doivent etre invalides:
select st_isvalidreason(geom) from parcelle where not st_isvalid(geom)
select st_isvalidreason(geom) from ebc  where not st_isvalid(geom)

Si ca se confirme, il faudra nettoyer les données (snapToGrid, st_buffer(geom, 0) par exemple.

Nicolas

En fait, il y a plusieurs facon de corriger des données qui ne sont pas valides au sens OGC. Le format Shapefile autorisant le stockage de géometries "bizarres", il arrive souvent apres chargement des shapefiles en base postgis que les polygones ne soient pas valides.

La methode snapToGrid réduit la précision des coordonnées en les "collant" sur une grille de taille donnée. Par exemple avec des données France métrique, il ne sert a rien de garder des chiffres apres la virgule. Ca n'a pas de sens au niveau geo (sauf si les données ont été numérisées avec une précision 10cm ou cm, mais c'est plutot rare ).

La methode buffer(geom, 0) arrive souvent a "nettoyer" des geometries invalides en reconstruisant le polygone et en virant les petites intersections. Le 0 veut dire: buffer de 0 m, soit la geometrie elle-meme.

Je vous recommande par ailleurs la lecture de la doc Postgis, qui contient deux chapitres d'introduction bien utiles pour comprendre le mécanisme de PostGIS.

Nicolas