Le portail francophone de la géomatique

emmamillefeuille · Wed 19 April 2023 19:10

Bonjour à tous !

Je découvre depuis peu le langage SQL et je me heurte aujourd'hui à un soucis :

Je cherche à faire l'intersection spatial entre une couche shp linéaire (des cours d'eau) et une couche polygone (des communes). Mon objectif est d'obtenir une nouvelle table/couche contenant le nom de chaque cours d'eau et la/les commune qu'il traverse. Les deux couches n'ont aucun id en commun.

Pour cela, j'aimerais utiliser une requête SQL car mes données sont dans une base (et je n'ai pas juste envie de faire une intersection "manuelle" via la boîte à outils car c'est un exercice que je vais réappliquer sur plusieurs jeux de données).

J'étais initialement partie sur quelque chose comme ceci :

Code:

Select nom_cours_eau, commune, CE.geometry
From cours_eau CE, communes
Where st_within(CE.geometry,communes.geometry)

mais ce code renvoie les cours d'eau en plusieurs exemplaires (je ne sais pas pourquoi, je me retrouve donc avec une table et plusieurs lignes identiques avec la même géométrie) et d'autre part, ne renvoie pas du tout les cours d'eau à l'intersection de plusieurs communes.

Dans l'idéal, j'aimerais bien obtenir une table du style :

cours d'eau // commune
CE 1 // com1
CE2 // com 1
CE2 // com2

Merci d'avance pour votre aide (en espérant avoir assez détaillé)

Miniopterine · Wed 19 April 2023 19:37

Bonjour,

https://gis.stackexchange.com/questions … -st-within
A priori le st_within ne renverra que les cours d'eau entièrement contenu dans une seule commune.

Essayez plutôt le St_intersects

Lsam · Thu 20 April 2023 00:05

Bonjour,

Si vous ne souhaitez qu'une ligne par cours d'eau, alors il faut utiliser une fonction d'agrégation sur vos communes, qui va vous permettre de les regrouper par cours d'eau.
Exemple :

Code:

SELECT ce.nom, group_concat(com.commune, ', ') communes, ce.geometry
FROM cours_eau ce, communes com
WHERE st_intersects(ce.geometry,com.geometry)
GROUP BY ce.nom, ce.geometry
;

La fonction group_concat est du SQL utilisé par SQLite. Elle permet de concaténer les valeurs d'un champ texte de plusieurs enregistrements. J'ai choisi comme séparateur ', ' mais ça peut être n'importe quelle autre chaîne de caractère(s).

Lsam · Thu 20 April 2023 09:27

J'avais peut-être mal compris votre demande. Si vous souhaitez n'avoir que les tronçons de cours d'eau coupant chaque commune :

Code:

SELECT ce.nom, com.commune, st_intersection(ce.geometry, com.geometry) geometry
FROM cours_eau ce, communes com
WHERE st_intersects(ce.geometry,com.geometry)
;

emmamillefeuille · Thu 20 April 2023 09:30

je vous remercie pour les réponses ! je vais tester tout ça. En effet je souhaite uniquement des tronçons mais je n'avais pas penser à grouper et concaténer, je pense que ce sera parfait

emmamillefeuille · Thu 20 April 2023 12:30

le 2e code fonctionnement parfaitement pour ce que je recherchais, merci !

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Wed 19 April 2023 19:10

QGIS BDD SQL: Intersection spatial entre lineaire et polygone

Code:

#2 Wed 19 April 2023 19:37

Re: QGIS BDD SQL: Intersection spatial entre lineaire et polygone

#3 Thu 20 April 2023 00:05

Re: QGIS BDD SQL: Intersection spatial entre lineaire et polygone

Code:

#4 Thu 20 April 2023 09:27

Re: QGIS BDD SQL: Intersection spatial entre lineaire et polygone

Code:

#5 Thu 20 April 2023 09:30

Re: QGIS BDD SQL: Intersection spatial entre lineaire et polygone

#6 Thu 20 April 2023 12:30

Re: QGIS BDD SQL: Intersection spatial entre lineaire et polygone

Pied de page des forums