Message posté par : Pascal PLUVINET
----------------------------------------
Bonjour,
En SQL sous PostGIS, on écrirait :
-----------------
Code :
create table mon_resultat as
select (st_dump(st_intersection(a.geom,st_union(b.geom)))).geom as geom
from data a, data b
where st_intersects(a.geom,b.geom)
group by a.geom
;
-----------------
Même si je maitrise R, je ne connais pas bien la librairie sf. Mais peut-être que ca donne
des pistes pour la rédaction sous R ou sous python.
De plus, si les polygones sont nombreux, il faudrait qu'ils soient indexés
préalablement pour une meilleur performance.
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=368158#p368158
Pour y répondre : geomatique(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