Le portail francophone de la géomatique

Bonjour,

Effectivement si vous avez ce genre de messages d'erreur il s'agit probablement d'un problème de topologie :

L'outil pour vérifier la validité des géométrie est native:checkvalidity la doc en français se trouve au lien suivant : https://docs.qgis.org/3.40/fr/docs/user … k-validity

Ensuite, une fois mes problèmes détectés, j'essaie de corriger manuellement les erreurs. S'il y a trop d'erreurs pour que la corrections soit faite manuellement, il y a l'outil native:fixgeometries qui peut corriger automatiquement ces erreurs (lien vers la doc : https://docs.qgis.org/3.40/fr/docs/user … ometries).


En fouillant sur le net, j'ai trouvé cette formation éditée par l'INRAP, qui a un chapitre interessant sur la validité des géométrie et la topologie : https://formationsig.gitlab.io/fiches-t … _geom.html

Bonne soirée.

Bonjour,
Il faudrait aussi vérifier qu'il y a bien des index spatiaux pour ces couches.
Votre table des communes métropolitaines est-elle dans la même base de données Geopackage ?

L'intersection se fait :
- avec les outils QGis ?
- ou par une requête SQL dans le gestionnaire DB (et si oui directement dans le fichier geopackage, ou dans les couches virtuelles ?)

Bonjour,
les 2 GPKG sont sur le reseau de l'entreprise, l'intersection se fait avec les outils QGIS. Je vais tenter de mettre le GPKG sur C:\temp et voir ce qui se passe.

Avez vous essayé d'intégrer vos données dans Postgis pour réaliser ce traitement, bien plus adapté qu'un gpkg ?

Bonjour, pas testé Postgis, j'ai testé à peu près tous les formats compatibles Qgis y compris les bases de données Google version spatiale. Spatialite en mode SQL par ex n'etait pas adapté.
Question si j'installe Postgis ce sera sur PC perso de travail. Les performances seront acceptables ?

@GlaDal => Le traitement Réparer les géométries est efficace et il permet de mener a bien les intersections mais qui demeurent tres longues.

Bonjour,
Etant un utilisateur d'outils SIG depuis 28 ans, je vous conseille d'effectuer les calculs d'intersection dans une base PostgreSQL/Postgis si vous préférez les outils libres.
Il faudra alors écrire un peu de code SQL pour réaliser l'opération d'intersection (ST_intersection) sur vos 2 couches indexées comme dans cet exemple :

On suppose que les 2 couches ont été préalablement importées dans la base Postgis (ici dans le schéma public)

CREATE INDEX ON public.layer1 USING GIST (geom);
CREATE INDEX ON public.communes USING GIST (geom);

DROP TABLE IF EXISTS public.layer1_par_commune;
CREATE TABLE public.layer1_par_commune AS
SELECT
    g1.*,
    c.nom AS nom_commune
   
    CASE
        WHEN ST_CoveredBy(g1.geom, c.geom)
            THEN g1.geom
        ELSE
            ST_Multi(
                ST_Intersection(g1.geom, c.geom)
                     )
           
    END AS geom

FROM public.layer1 AS g1
INNER JOIN public.communes AS c
ON ST_Intersects(g1.geom, c.geom);

Je n'avais pas suivi ce fil plus tôt, mais selon moi, le goulot d'étranglement se situe ici !
Travailler sur des données si volumineuses en SIG nécessite un accès très rapide aux données, de préférence donc un disque local SSD.
Si vous travaillez sur un lecteur réseau, il est normal que vos temps de traitements soient interminables !

Après, pour utiliser PostGis très régulièrement, c'est vrai que c'est encore plus performant que de travailler sur des fichiers (GPKG ou autres).
Mais si les fichiers ne sont pas locaux, c'est catastrophique !

Dernière modification par Sylvain M. (Wed 01 April 2026 16:32)