Le portail francophone de la géomatique

Bonjour,
Depuis la version 2.14 il est possible de faire du SQL (de spatialite) avec n'importe quelle couche.

https://www.qgis.org/fr/site/forusers/v … ual-layers

Bonjour à vous,

Je reviens sur cette question.
J'ai bien noté (même si pas encore pratiqué) que QGis 2.14 permettait de faire des requêtes SQL spatiales avec n'importe quelle couche.
Pour le géomaticien que je suis c'est une très bonne nouvelle :-)
Mais pour mes collègues utilisateurs de QGis, le SQL spatial n'est pas encore à leur portée.  :-/
Hors, une requête de type "les communes qui sont dans le périmètre administratif Truc" me semble être une opération basique qui devrait être à leur portée.
Je n'y suis pas parvenu avec des données pourtant cohérentes d'un point de vue topologique (issue du même millésime de GeoFla Communes).
La capture d'écran ci-jointe sera plus parlante.

La seule solution que j'ai trouvé (et utilisé), c'est de créer une couche des centroïdes des communes, pour faire la sélection sur une couche de points (intersection).
N'y aurait-il pas plus simple ? (dans ArcGis Basic, on a bien une sélection spatiale de type "les entités de la couche cible ont leur centroïde dans les entités de la couche source")

Merci à vous.

Sylvain M.

SOURCE :
[img]http://georezo.net/forum/attachment.php?item=6983&download=1[/img]

RÉSULTATS DE SÉLECTIONS :
[img]http://georezo.net/forum/attachment.php?item=6984&download=1[/img]

Dernière modification par Sylvain M. (Tue 03 May 2016 14:49)

Bonjour,

J'ai été confronté à ce problème car je voulais lister les communes qui se trouvaient à l'intérieur d'une CDC. Les résultats obtenus étaient erronés à cause de l'imprécision de mes géométries.

J'ai réussi à contourner le problème en utilisant la fonction suivante :

select * from public.t_commune JOIN public.t_epci ON st_within(ST_centroid(t_commune.geom) , t_epci.geom)

Cette requête SQL permet donc de retourner les données des deux tables et la liaison entre les deux tables se fait sur la fonction st_within(ST_PointOnSurface(t_commune.geom) , t_epci.geom) qui permet dans un premier temps de récupérer le centroïde de la commune(ST_centroid) et dans un second temps de vérifier la présence de ce centroïde à l'intérieur de l'EPCI (st_within)

Cordialement,
Sylvain THIBAULT

Salut,

On ne parle donc pas des mêmes outils, ce sont tous des extensions Core (le mien sous Ftools, le tien. On pourrait même y rajouter un troisième dans le menu Traitement (select by location). Il y a peut-être un ménage à faire de ce côté...
Essaie donc les deux autres (d'ailleurs avec une interface assez proche tous les trois) et tiens-nous informés des résultats.

Pour info, je viens d'installer QGis 2.14
- la requête en SQL Spatial [ST_centroid()] sur couches virtuelles marche nickel : pour moi c'est top.
- les différentes requêtes spatiales "presse-bouton" ne permettent toujours pas de faire cette sélection "basique"
Si quelqu'un veut tester, je peux lui envoyer mes données en MP, pour vérifier si c'est moi qui suis pas bien réveillé
A+

Sylvain M.

Salut Sylvain,
Tu peux toujours faire la demande mais je ne pense pas qu'elle couvrirait forcément le besoin. Peut-être dans ce cas-ci mais quid si une commune avait une forme telle que son centroide n'était pas en son intérieur (st_centroid vs st_pointonsurface) et que cette commune était en limite de territoire agrégé? Par exemple, le bout d'île au nord de La Manche? tu pourrais ne pas la sélectionner celle-là, s'il s'agit d'une commune en deux parties et que le centroïde tombe en pleine mer.

Perso, je préférerais une solide comparaison entre géométries polygonales (ou de base), plutôt qu'avec un centroïde.

J'ai donc aussi fait des tests avec geofla et une couche dep normands générée virtuellement à partir des communes avec

select CODE_DEP, NOM_DEP, CODE_REG, NOM_REG, ST_UNION(geometry) as geom
from (select * from COMMUNES where CODE_REG='NORMANDIE')
group by CODE_DEP, NOM_DEP, CODE_REG, NOM_REG

Et mes tests de sélection (avec l'outil du menu Recherche, une option à la fois) m'amènent, sans avoir consulté le code des fonctions toutefois, à ces règles:
- "intersectent" renvoie toutes les communes qui partagent une surface ou une limite avec un département (ici, même chose que toi)
- "touchent" concerne les communes qui partagent une limite avec un département sans partage de surface, donc les communes limitrophes de la Normandie
- "chevauchent" : une partie de la surface (mais pas toute) de la commune est contenue dans un département, soit aucune commune renvoyéevu que tout est topologiquement nickel
- "contenues" renvoie les communes qui ont toute leur surface à l'intérieur d'un département, soit bien TOUTES les communes du département ou de la région Normandie.

Je n'ai pas testé "Requête Spatiale" car extrêmement long (j'ignore si c'est ma config qui en est la cause, et je sais être patient mais là j'abandonne l'idée de reprendre les tests dessus- il n'y a pas photo entre les deux outils, d'autant que tu ne peux même pas annuler la requête). Quelqu'un pourrait me confirmer cet état afin que je signale aux devs, svp? Merci.

J'en viens donc à m'interroger sur la qualité topologique de tes données. Pourrais-tu utiliser le vérificateur topologique de QGIS pour voir les relations spatiales entre les deux couches? Et en relisant tes premiers messages, tu parles de superposer Geofla et bd-Topo qui n'ont pas la même précision et, si c'est ce que tu fais, cela renverra logiquement des résultats imprévisibles.

Bonjour,

Je me permet de relancer le sujet. Un peu comme Sylvain, j'aimerai sélectionner des communes qui sont à l'intérieur des entités de ma couche EPCI afin de renseigner automatiquement mes ma table commune. Sauf qu'il semble y avoir un léger décalage entre ma couche commune et ma couche EPCI.
La fonction within ne marche pas correctement dans mon cas. En lisant votre message, j'ai donc combiner les fonctions ST_Within et ST_Centroid sur pour une couche virtuelle.
Mais la couche virtuelle ne permet pas de faire des modifications.
Si quelqu'un a une solution, je suis preneur.

Pour info, j'ai un message presque similaire sur le forumsig qui pour l'instant n'a pas trouvé de réponse: http://www.forumsig.org/showthread.php/ … e-spatiale
Si ça change, je vous tiendrai au courant.

Merci d'avance!
Cordialement,
Elminster

Non, mais tu pourrais ensuite mettre à jour ta couche d'origine via une jointure attributaire à ta couche virtuelle ?
Mais, comme toi, je suis toujours aussi intéressé par un mode de jointure spatiale basée sur le centroïde de la couche source.
A+

Sylvain M.

Dernière modification par Sylvain M. (Tue 24 January 2017 17:36)

+1
Et sous ArcGis aussi d'ailleurs, de manière très simple : c'est la dernière option des méthodes de sélection spatiale ("selon l'emplacement") : "entités de la couche cible ont leur centroïde dans l'entité de la couche source"

Bonjour,
Utilisez l'algorithme de sélection par localisation ou son comparse "extraction par localisation", accessible depuis la boîte à outils du menu Traitement.