Message posté par : Nicolas Ribot ---------------------------------------- Bonjour, Pkoi calculer des centroids pour avoir la distance entre objets ? Ca sera faux dans certains cas (batiments en U, batiment avec cour, etc. (le centroid ne tombant pas dans le pg pour ces cas) Calculez directement la distance la plus proche, elle sera 0 si le batiment en touche un autre ! (utilisez l'operateur distance KNN <-> pour cela, et surtout pas st_distance(geom, geom), qui prendrait un temps fou. Sinon, pour regrouper des objets qui se touchent, plusieurs facons: st_union(geom) => tous les objets qui se touchent seront unis en un seul objet, ou plus rapide: st_clusterINtersecting: ca renvoie un tableau de geometryCollection contenant tous les objets qui se touchent, st_UnaryUnion pour avoir l'union de ces objets. Mais dans le cas que vous décrivez, ca me parait inutile. Nicolas ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=331583#p331583 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

Message posté par : Nicolas Ribot ---------------------------------------- Un exemple de SQL pour regrouper les objets qui se touchent (st_clusterIntersecting), en faire un polygone unique (st_unaryUnion), croiser cette couche avec elle même pour trouver le plus proche voisin (t.geom <-> t2.geom), croiser avec les batiments pour faire un lien entre chaque cluster de batiments et les batiments qui le composent (array_agg(id), fabriquer une geom qui est la plus petite distance entre deux clusters de batiments (st_shortestLine) : https://zupimages.net/viewer.php?id=20/14/095n.png ----------------- Code : with tmp as ( select st_clusterintersecting(geom) as arr from batiment ), tmp1 as ( select u.ordinality as idcluster, u as geom from tmp, unnest(tmp.arr) with ordinality as u ), tmp2 as ( select t.idcluster, t2.idcluster as idcluster2, t2.dist, st_unaryunion(t.geom) as geom, t2.geom as geom2 from tmp1 t cross join lateral ( select b.idcluster, b.geom <-> t.geom as dist, b.geom from tmp1 b where b.idcluster != t.idcluster order by b.geom <-> t.geom limit 1 ) as t2 ) select t.idcluster, t.idcluster2, dist, array_agg(b.id) as id_batiments, st_shortestline(t.geom, t.geom2) as geomline from tmp2 t join batiment b on st_contains(t.geom, st_pointonsurface(b.geom)) group by 1,2,3,5; idcluster idcluster2 dist id_batiments 1 3 49.8 {1} 2 3 97.2 {2} 3 1 49.8 {3,4,5,6,7,8} 4 2 176. {9} 5 6 78.3 {10,13} 6 7 41.7 {11} 7 6 41.7 {12} 8 1 99.1 {14} 9 8 100. {15} ----------------- Avec une vraie grosse table de batiments, il faut passer par des tables intermédiaires et créer des index spatiaux qui vont bien (geom, st_pointOnSurface(geom), ...) pour optimiser la requête. Nico ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=331590#p331590 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