Le constat : les calculs ne sont pas bons#

Date de publication initiale : 18 juillet 2024

Dans nos SIG, les opérations de superposition (overlay dans la langue de Shakespeare) telles que les intersections, les unions, les différences, etc. ainsi que l'accrochage utilisé par les dessinateurs, sont omniprésentes. Ces processus s'appuient sur des calculs similaires, simplifiés ici pour une meilleure compréhension dans cette présentation générale.

Cet article est la première partie de la série d'été sur la gestion de la géométrie dans les SIG.

Le dossier

Commenter cet article

Identification du problème#

Chargement des données#

Toutes les données utilisées sont disponibles sur mon GitHub, et pour simplifier la compréhension et la transposition de ces données dans différents SIG, j'utiliserai les formats WKB et WKT¹.

Revenons-en au fait, prenons un exemple avec une géométrie de type ligne, ici fermée, mais cela serait similaire pour un polygone (puisqu'un polygone est une ligne). Les géométries utilisées sont projetées dans le système de coordonnées EPSG:3946, projection de mon coin magnifique.

Exemple d'une géométrie de ligne au format WKB :

0102000000050000007997c6b68d3c3e4139eb62c260d55341ac9ea7316a3c3e41cbeb40e073d55341403e0bfbc33c3e41b3fc06f380d55341387a2a800c3d3e41f256b8176dd553417997c6b68d3c3e4139eb62c260d55341

À partir de cette ligne, je génère d'autres lignes qui s'accrochent aux points d'intersection/accrochage avec la ligne suivante :

010200000002000000ea9c6d2b873c3e41a03d941b7cd5534133db7796ce3c3e413fba569864d55341

Sur QGIS, on peut charger ces WKB - ainsi que les EWKB, EWKT, WKT - avec le très utile plugin QuickWKT.

Ce qui nous donne :

Intersections de ces lignes#

Nous avons désormais notre base pour étudier ce problème de précision/tolérance/intersection.

Justement, calculons l'intersection entre ces lignes. Pour cela, nous allons utiliser l'outil native:lineintersections de QGIS.

Nous obtenons deux points d'intersection. Visuellement, les résultats sont conformes à nos attentes, les points d'intersection se trouvant précisément sur les lignes.

Quand je dis « précisément », si l'on visualise sur une échelle totalement « absurde », on a toujours cette superposition.

Pour plus tard, on notera leurs WKB :

0101000000a899efc8c83c3e4175e5698166d55341

dont l'équivalent WKT est POINT(1981640.7849060092 5199258.022088398)

et

0101000000b5ebdd9e8f3c3e416bf8515379d55341

respectivement POINT(1981583.6205737416 5199333.301878075)

Création de lignes depuis ces intersections#

QGIS propose différentes options pour s'accrocher. On en utilisera deux, l'accrochage sur un sommet et l'accrochage aux intersections.

Accrochage sur les points d'intersection#

Tout d'abord, on va utiliser la fonctionnalité d'accrochage sur les intersections. Attention, ce n'est pas sur les points que l'on vient de générer, mais sur les intersections entre les géométries. On repère l'icône d'accrochage d'intersection avec une croix. L'accrochage sur un sommet avec un carré.

Dans la vidéo ci-après, je montre comment j'ai généré des lignes de part et d'autre de la ligne principale aux points d'intersection.

Accrochage sur les intersections#

Je répète l'opération, cette fois en me focalisant sur les points d'intersection. L'objectif est de garantir l'accrochage précis sur la ligne principale, indépendamment des variations des sommets adjacents.

Comparaison des points d'accrochage#

On fait confiance à QGIS et, visuellement, les lignes apparaissent bien accrochées à la ligne de base.

Comparons maintenant nos géométries textuelles. On a 8 lignes/segments dans les deux cas. Si l'accrochage se déroule correctement, on devrait avoir les mêmes points de départs.

line_intersection :

line_snap:

On remarque que nos points d'origine sont dans les deux cas :

• 1981583.62057374161668122 5199333.30187807511538267
• 1981640.78490600921213627 5199258.02208839822560549

respectivement en WKB :

• 0101000000b5ebdd9e8f3c3e416bf8515379d55341
• 0101000000a899efc8c83c3e4175e5698166d55341

Les lignes semblent donc avoir été correctement accrochées, comme le montrent nos représentations en WKT et WKB.

En effet, on retrouve bien les coordonnées d'origine :

b5ebdd9e8f3c3e416bf8515379d55341 et a899efc8c83c3e4175e5698166d55341 dans les points.

Sélection des lignes intersectant la base#

Par conséquent, nous avons des lignes qui sont sur le point d'intersection. Si l'on souhaite vérifier l'accrochage, on utilise le prédicat « intersecte » de QGIS. Pour vérifier cela, utilisons l'outil « sélection par localisation » :

Le prédicat intersects avec line_snap#

On regarde si line_snap est accroché sur base. Dans les exemples ci-après, même si je devrais n'utiliser que touches ou intersects, je vais tout cocher sauf disjoint.

Le résultat :

Aïe, seulement 2 sur les 4…

Le prédicat intersects avec line_intersection#

On fait de même pour line_intersection :

Et voici le résultat :

Re Aïe. Seulement 2 sur les 4, et du même côté. C'est au moins ça de cohérent dans notre incohérence…

Est-ce que ce n'est que pour ces deux points ?#

Pour aller plus loin, je réalise des tests aléatoires en essayant d'accrocher la ligne principale à différents endroits. C'est la couche test_line dans le GeoPackage.

Dans notre franchise OSGeo, QGIS est un de nos superhéros, tel Spider-Man. Notre petite toile d'araignée manque de glu, car notre accrochage n'est toujours pas bon :

Testons tous les 1 mètre sur la base#

Il y a quelques années, j'ai développé un outil pour créer des transects. Il génère des bandes, suivant un angle donné. Il est principalement utilisé pour des analyses sur des profils en travers, en génie civil, en écologie, etc. Ici, on va réaliser un transect avec l'outil… « transect », tous les mètres de part et d'autre de la ligne principale.

Pour commencer, nous allons densifier la géométrie :

Ensuite, on génère de cette couche les transects à gauche et à droite de la ligne (suivant le sens de la ligne) :

Notre dessin ressemble à :

Et maintenant, on va sélectionner les entités de transect_left et transect_right qui intersectent la base :

Et, non, ce n'est pas l'algo de transect qui est tout buggué.

Mais, comme on dit… "Caramba ! Encore raté !"

Et si je teste sur la géométrie densifiée ?

Et notre résultat ?

Cela fonctionne ! Mais pourquoi ?

Et bien, nous verrons cela plus tard.

2 : QGIS et GEOS

Auteur·ice#

Loïc Bartoletti#

Après un cursus en Histoire, je me suis orienté vers l'urbanisme sur l'aménagement des territoires.

J'ai travaillé pendant environ 10 ans dans une station touristique dans les Alpes, Megève, en tant qu'urbaniste puis responsable du bureau d'études et administrateur SIG.

Bidouilleur et partisan des solutions OpenSource, j'ai commencé à toucher à GRASS, puis QGIS et PostGIS. Au fil du temps j'ai contribué à ces logiciels, principalement pour migrer des outils DAO vers le SIG et je suis aujourd'hui commiter QGIS, PostGIS et FreeBSD où je m'occupe des paquets des outils OSGeo et plus si affinité.

Licence Beerware #

Ce contenu est sous licence Beerware (Révision 42).
Les médias d'illustration sont potentiellement soumis à d'autres conditions d'utilisation.

Commentaires