Le portail francophone de la géomatique

Bonsoir,

L'utilisation de l'opérateur @ seul (et tous les autres aussi d'ailleurs &&, &<, etc) seront faux car ils utilisent les BBOX, donc il faut faire un deuxième passage. En réalité c'est ce que fait st_within() (ie ils utilisent les index donc les bbox puis fait un deuxième passage).

Par contre je ne sais pas quel opérateur est utilisé par st_within, peux être testé avec ST_DWithin et utiliser une distance de 0.000001 ?

Y.

Bonjour,

Outre le fait que c'est un résultat approché comme vous le soulignez, la requête est un peu bizarre dans sa forme

select A.insee_comm, A.geom @ B.geom as test

va renvoyer une géométrie et un booléen.

de même

where A.geom @ B.geom = true

Le "=true" est inutile.

Concernant les fonctions postgis je ne connais pas leur code, mais il existe deux algorithmes principaux pour savoir si un point est à l'intérieur d'un polygone. Soit ce que je nommerai l'algo de la ligne de balayage (1) soit l'algo topologique (2).

1) : On parcours le polygone de bas en haut avec une ligne parallèle à l'axe des X et on vérifie si le point appartient aux intervalles résultant de l'intersection de cette ligne et du polygone.(on l'utilise pour de l'image).

2) On parcours le polygone selon ses sommets (les anneaux intérieurs étant orientés dans le sens inverse de l'anneau extérieur) et on vérifie que le point concerné est toujours du même coté avec une fonction du type EstAGauche() : boolean .


Ces deux algos dépendent du nombres de sommets du polygone.

En revanche avec l'utilisation de l'opérateur @ seuls les points appartenant à au moins deux BBox ont a subir une deuxième passe.

Bonjour,
Etant donné que ST_WITHIN fait appel à la bbox, il me semble que ce n'est plus utile d'ajouter une comparaison de bbox avec

 AND territoire.geom && chl.geom

Pour ce qui est des performances, je ne connais pas le nombre d'entités de la table adm.territoire ni la complexité de leur géométrie mais peut-être que finalement, le temps de requête n'était pas si énorme (A part la sous-sélection, on a pratiquement la même requête que celle du début).
Que vaut la même requête sans le préfiltrage sur les territoires par ex?

WITH territoire AS (
            SELECT geom 
            FROM adm.territoire 
            )
        SELECT chl.gid, chl.insee_comm, chl.geom FROM adm.chl AS chl, territoire
        WHERE st_within(chl.geom, territoire.geom);

Bonjour,

Je reviens un peu tard sur ce post mais voici les résultats obtenus :
- avec la dernière requête de Nicolas Ribot : j'obtiens mes 7800 lignes en 19 secondes
- avec celle de Santanna : 17 secondes
- avec la requête de Santanna, mais en remplaçant le ST_Within par ST_Dwithin : 19 secondes

Du coup, bizarrement, ST_Within est plus rapide sur ce coup ?

Sinon, en prenant une couche "Territoire" avec des entités plus petites, comme les communes au lieu des régions par exemple (et donc des bbox à la fois plus petites et nombreuses), la requête est beaucoup plus rapide : moins de 2 secondes, que ce soit avec ST_Within ou ST_DWithin.
Pour le coup, ST_Dwithin est légèrement plus efficace mais ça se joue à un chouilla...

Merci en tout cas et bonne journée

Bonjour Jérémie,

Je viens de tester ta requête (car j'ai effectivement seulement 2 territoires mais assez grands) que j'ai dû modifier car st_within tout court me renvoyait une erreur ("Relate Operation called with a LWGEOMCOLLECTION type.  This is unsupported").
Je suis donc passé en st_dwithin.
Mais elle se termine en 38 secondes.