Le portail francophone de la géomatique

Sylvain230 · Tue 18 August 2026 20:33

Bonjour,

J'apprends depuis peu à effectuer des requêtes SQL sur QGIS, et je commence tout juste à croiser les données entre les couches.

Je me heurte à une requête que je n'arrive pas à trouver :

Je souhaite découper la couche "cours d'eau" aux limites d'un département. (Je connais la possibilité "Vecteur/ Outils de géotraitement / Couper") mais je cherche la requête SQL qui permet de faire cela. (Donc découper une couche polylignes par rapport à un polygone). je précise que ma couche département est filtrée sur le département ciblé.

J'ai testé plusieurs options, entre autre :

Code:

SELECT *
FROM "Cours","Departements_EPSG_3857"
WHERE "Departements_Français".'nom'='Ain'
AND st_intersects("Cours".'geometry',"Departements_EPSG_3857".'geometry’);

ou cette requête que j'ai adaptée à partir d'une solution trouvée sur le forum

Code:

SELECT "Cours.nom", group_concat(Departements_Français, ', ') "Departements_Français", "Cours".'geometry'
FROM "Cours", "Departements_Français"
WHERE "Departements_Français".'nom'='Ain'
AND st_intersects("Cours".'geometry',"Departements_Français".'geometry')
GROUP BY "Departements_Français".'nom', "Cours".'geometry';

Enfin, je précise que j'ai testé, avec succès! la requête suivante : (Là ce sont 2 couches polygones)

Code:

SELECT *
FROM "ContextePiscicole_FXX","Departements_Français"
WHERE "ContextePiscicole_FXX".'CdCtxPisci'='Salmonicole'
AND "Departements_Français".'nom'='Ain'
AND st_intersects("ContextePiscicole_FXX".'geometry',"Departements_Français".'geometry');

pour conserver le contexte piscicole souhaité dans le département ciblé.

Si vous avez la solution SQL pour découper ma couche cours d'eau aux limites du département, je suis preneur...
Merci !

Dernière modification par Sylvain230 (Tue 18 August 2026 21:20)

Sylvain230 · Sun 23 August 2026 12:29

Bonjour,
Je remets un message par rapport à mes "avancées"...
Si j'ai bien compris, St intersect ou intersection permettent de découper des polygones par rapport à des polygones.
La fonction qui permettrait de découper une couche de lignes par rapport à un polygone serait la fonction St_Split

J'ai donc essayé la requête suivante, sans succès...

Code:

SELECT *
FROM "Cours_Eau_EPSG4326_WGS84","Departements_Français"
WHERE "Departements_Français".'nom'='Ain'
AND ST_Split("Cours_Eau_EPSG4326_WGS84".'geometry',"Departements_Français".'geometry');

La requête s'exécute, j'ai bien les champs de ma couche "cours d'eau" qui apparaissent dans le résultat en dessous, mais sans aucunes valeurs dans ma table... Et ensuite impossible de charger la couche (QGIS "plante")

Une solution ?

Sebastien Wamster · Mon 24 August 2026 08:43

Salut Sylvain,

Le souci vient probablement de ST_Split, ce n'est pas la bonne fonction ici. Elle sert à découper une géométrie A avec une géométrie B, mais elle ne fonctionne pas comme ça dans un join classique, d'où le plantage.

Ce qu'il te faut, c'est ST_Intersection, qui renvoie directement la portion de ligne contenue dans le polygone :

Code:

SELECT 
    "Cours_Eau_EPSG4326_WGS84".nom,
    "Cours_Eau_EPSG4326_WGS84".longueur,
    CastToMultiLinestring(
        ST_Intersection("Cours_Eau_EPSG4326_WGS84".geometry, "Departements_Français".geometry)
    ) AS geom
FROM "Cours_Eau_EPSG4326_WGS84", "Departements_Français"
WHERE "Departements_Français".nom = 'Ain'
AND ST_Intersects("Cours_Eau_EPSG4326_WGS84".geometry, "Departements_Français".geometry);

Dernière modification par Sebastien Wamster (Tue 25 August 2026 10:42)

Sylvain230 · Mon 24 August 2026 20:18

Merci beaucoup Sébastien pour votre réponse, je suis en train de tester je vous tiens au courant.

Donc cela veut dire aussi qu'il faut écrire sous le "SELECT" une ligne pour chaque champ que l'on veut dans la requête si je comprends bien ?
Je précise que je passe par l'onglet "gestionnaire de bases de données" "couches virtuelles" et "requêtes SQL" pour écrire mes requêtes.

Les requêtes SQL sont (en tout cas sur mon ordi) extrêmement longues à aboutir qd il s'agit de croiser les couches entre elles (plusieurs heures)...
Je ne sais pas si c'est normal ou pas, et s'il y a un moyen de les accélérer ou pas... ?

La fonction CastToMultilinestring est la fonction SQL qui permet de couper une couche de polylignes aux limites d'un polygone c'est bien ça ?

Dernière modification par Sylvain230 (Mon 24 August 2026 20:25)

Sebastien Wamster · Tue 25 August 2026 10:49

Oui, il faut lister chaque champ voulu après le SELECT. En SQLite/Spatialite, il n'existe pas de syntaxe pour faire une exception sur juste une seule colonne.

Pour la lenteur, oui c'est probablement plus rapide avec un index spatial, il faut utiliser CreateSpatialIndex sur la couche, si c'est un GeoPackage ou du Spatialite car ça ne se fait pas nativement.

Pour CastToMultiLinestring il ne sert pas à découper mais à éviter le bug avec le mélange de points et lignes qui fait planter QGIS.

Sylvain M. · Hier 09:50

Si je peux me permettre, voici comment j'écrirais la requête de Sébastien proposée dans le message de lundi matin :

Code:

SELECT
    ce.nom,
    ce.longueur,
    CastToMultiLinestring(
        ST_Intersection(ce.geometry, df.geometry)
    ) AS geom
FROM
    "Cours_Eau_EPSG4326_WGS84" ce,
    "Departements_Français" df
WHERE
    df.nom = 'Ain'
AND
    ST_Intersects(ce.geometry, df.geometry);

C'est beaucoup plus lisible et maintenable en utilisant des alias.
Tout sur les Alias ici : https://sql.sh/cours/alias

Sinon, je confirme aussi que les index spatiaux sont fondamentaux pour ce genre de requêtes spatiales d'intersections. Les temps de traitements sont souvent divisés par 100 ou plus !

Sylvain230 · Hier 19:49

Bonjour,
C'est parfait ! Je m'en suis sorti en adaptant la requête que m'avait proposé Sébastien

Bon alors pour avoir fait quelques recherches (mais comme je disais je ne suis pas encore très habitué aux requêtes spatiales) l'erreur principale que j'avais faite était d'utiliser une couche "Cours d'eau" qui était en format GeoJson, format qui ne dispose pas d'index spatial et est donc incompatible avec les requêtes spatiales comme celles-ci.

J'ai donc converti ma couche en format Géopackage, bien mis mes couches dans les mêmes EPSG/SCR, adapté la requête que m'avait suggéré Sébastien, et ça a fonctionné !

Merci beaucoup Sébastien !
Et Sylvain j'ai noté aussi, effectivement ça à l'air plus simple, je testerai probablement cette solution aussi.
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

Les inscriptions sont ouvertes ! L'association GeoRezo sera présente

#1 Tue 18 August 2026 20:33

QGIS: Decouper couche "Cours d eau" sur un depart. avec requete SQL

Code:

Code:

Code:

#2 Sun 23 August 2026 12:29

Re: QGIS: Decouper couche "Cours d eau" sur un depart. avec requete SQL

Code:

#3 Mon 24 August 2026 08:43

Re: QGIS: Decouper couche "Cours d eau" sur un depart. avec requete SQL

Code:

#4 Mon 24 August 2026 20:18

Re: QGIS: Decouper couche "Cours d eau" sur un depart. avec requete SQL

#5 Tue 25 August 2026 10:49

Re: QGIS: Decouper couche "Cours d eau" sur un depart. avec requete SQL

#6 Hier 09:50

Re: QGIS: Decouper couche "Cours d eau" sur un depart. avec requete SQL

Code:

#7 Hier 19:49

Re: QGIS: Decouper couche "Cours d eau" sur un depart. avec requete SQL

Pied de page des forums