Le portail francophone de la géomatique

Bonjour,

Ne faites pas de buffer, mais utilisez st_dwithin pour trouver les objets a moins de x m d'autres objets.
Pour faire un jointure sur la meme couche, listez-la 2 fois dans le from, avec des alias différents

select ...
from table1 t1 join table1 t2 on t1.id < t2.id and st_dwithin(t1.geom, t2.geom, 500)
...

Nicolas

Ca permet de ne tester qu'une seule fois un couple de 2 polygones.
Sinon, on teste A vs B et B vs A si on ne fait qu'une jointure spatiale.

Avec un exemple c'est plus parlant:
premiere requete: t1.id < t2.id: un seul test pour un couple de 2 id
deuxieme requete: t1.id <> t2.id: chaque couple d'id est testé 2 fois

with tmp as (
    select * from (values (1), (2), (3)) as v(id)
) select t1.id, t2.id
from tmp t1 join tmp t2 on t1.id < t2.id;
-- 1|2
-- 1|3
-- 2|3

-- vs un self join en excluant uniquement les memes objets (meme id)
with tmp as (
    select * from (values (1), (2), (3)) as v(id)
) select t1.id, t2.id
from tmp t1 join tmp t2 on t1.id <> t2.id;
-- 1|2
-- 1|3
-- 2|1
-- 2|3
-- 3|1
-- 3|2

Nicolas

Dernière modification par Nicolas Ribot (Fri 25 February 2022 12:41)