Comment trouver le systèmes de coordonnées d’une couche avec QGis 3.X?

Quand vous travaillez avec des données spatiales dans un logiciel SIG, vous rencontrez régulièrement le problème de fichiers qui n’ont pas de système de coordonnées clairement défini.

Lorsqu’un fichier manque de son système de coordonnées (SCR), il peut entraîner des couches mal alignées, des mesures incorrectes et une foule d’autres problèmes qui entravent la visualisation et l’analyse précises des données.

Dans cet article nous allons essayer de vous donner les conseils et outils disponibles en 2023 pour identifier les systèmes de coordonnées des fichiers avec des projections inconnues.

Avant de traiter le problème, jouons cartes sur table : vous pouvez y être pour partie ou en totalité responsable du problème. On va donc avoir deux chapitres dans cet article: le premier pour ceux qui n’ont pas l’habitude des systèmes de coordonnées, le deuxième pour ceux qui auront déjà fait par eux-mêmes les vérifications du premier chapitre et qui, donc, pourront aller directement au deuxième.

N’oubliez pas que les solutions suggérées ici ne résoudront jamais 100% des problèmes…

Vous êtes débutant en systèmes de coordonnées

Le problème peut se présenter à deux occasions:

  • Vous avez créé une nouvelle couche de données avec un traitement spatial de QGis et le résultat ne se superpose pas aux autres données.
  • Vous essayez de charger de nouvelles données et elles ne se superposent pas du tout avec les autres données du projet.

La couche est le résultat d’un traitement

Dans ce cas, votre souci vient du programme de traitement qui n’a pas su, ou pu, affecter le SCR de votre projet à sa sortie.

Pour connaître le SRC du projet, allez dans le menu Projet -> Propriétés, et dans l’onglet SCR

Le panneau inférieur vous indique le nom et les paramètres du SCR de votre projet, mais retenez surtout le code qui apparaît en bout de ligne dans la liste des systèmes récemment utilisés, dans cet exemple EPSG:2154.

Pour indiquer à QGis que votre couche utilise ce même système, cliquez avec le bouton droit sur la couche -> SCR de la couche

La première ligne, Aucun SCR, correspond à l’état de description du SCR connu par QGis. Pour lui affecter le SCR du projet, sélectionnez la ligne correspondante au SCR de la fenêtre précédente.

Si par contre la première ligne n’est pas Aucun SCR comme ici:

Cela veut dire que votre couche a été affectée d’un mauvais SCR par le programme de traitement. Ceci peut arriver quand votre projet et des couches qu’il contient ne sont pas dans le même système de coordonnées. QGis sait gérer ça, mais les traitements ont des origines diverses et peuvent bugger.

Regardez le SCR de la couche en entrée du traitement. S’il n’est pas le même que celui du projet, la solution consiste, généralement, à changer le SCR de la couche problématique par l’autre.

La couche provient de données externes que vous ne connaissez pas.

Vérifiez les métadonnées et les fichiers qui les accompagnent. Dans le plus grand nombre de cas, les fichiers sont au format shapefile. Un fichier shape se compose de plusieurs fichiers obligatoires et optionnels, à savoir:

  • .shp : le fichier principal contenant les données géométriques. (obligatoire)
  • .shx: Le fichier index pour un accès rapide aux données géométriques. (obligatoire)
  • .dbf: La table dBASE contenant les données d’attributs. (obligatoire)
  • .prj: Le fichier de projection contenant les informations du système de coordonnées (si disponible)

En général le problème de définition du SCR vient de l’absence du fichier .prj . Si ce fichier est manquant ou vide, vous devrez identifier le système de coordonnées vous-même. Si le fichier est présent et rempli, vous pouvez l’ouvrir avec un éditeur de texte pour afficher la description du système de coordonnées.

Voici un exemple de fichier .prj rempli:

S’il est absent ou vide, il faudra que vous le définissiez dans propriétés de la couche. Mais pour cela il faudra trouver quel est le système: recherchez les fichiers facultatifs, comme . txt ou . xml, qui pourraient contenir des métadonnées ou des renseignements supplémentaires sur les données; parfois, les fichiers shapefile sont accompagnés de fichiers de documentation, comme un README.txt ou un rapport PDF, qui fournissent des informations.

Si vous ne trouvez pas l’information dont vous avez besoin dans les métadonnées ou les fichiers connexes, songez à communiquer avec le fournisseur de données ou la source. Ils peuvent être en mesure de fournir les renseignements nécessaires sur le système de coordonnées.

Si vous trouvez la réponse, il faudra indiquer à QGis le SCR à utiliser. Cliquez sur la couche, menu SCR de la couche ->SCR de la couche :

Vous ouvrez la fenêtre de définition de SCR.

A l’aide du champ filtre, rentrez des parties du nom du SCR recherché. Vous aurez la liste des SCR contenant cette partie dans la liste Systèmes de Coordonnées de Référence Prédéfinis.

Si arrivé jusqu’ici vous n’avez pas résolu votre problème, il va vous falloir prendre votre courage à deux mains et passer au chapitre suivant.

Vous savez que le SCR des données est inconnu et vous cherchez à le découvrir.

En 2015 j’avais mis en ligne un document « Déterminer la projection inconnue d’une couche pour QGis 2.X » qui est toujours disponible en cliquant ici. Il n’est valable que pour les projections courantes de la France Métropolitaine et se base sur l’analyse des valeurs des XY des données inconnues. Heureusement que QGis et ses développeurs sont passés par là depuis, car QGis 3.x vous propose un outil de recherche de système inconnu, dans la boîte de traitements -> Outils généraux pour les vecteurs -> Trouver une projection.

L’outil de recherche de projection dans QGIS est conçu pour créer une liste restreinte de systèmes de référence de coordonnées (SCR) candidats pour une couche avec une projection inconnue, basée sur une zone cible spécifiée à l’aide d’une autre couche de référence, avec un SCR bien défini.

Par exemple, si vous savez ou pensez que vos données se trouvent quelque part autour d’Arcachon, vous chargez une couche de données comprenant les Landes :

  • Limites administratives (p. ex., frontières de pays, de Région, de Département)
  • Réseaux routiers
  • Imagerie satellite ou photographie aérienne (p. ex., Google Maps, OpenStreetMap)
  • etc…

Dans cet exemple nous avons une couche nbm avec un SCR inconnu, mais nous savons qu’elle est en relation avec le Bassin d’Arcachon.

L’outil Trouver une projection va prendre les données de nbm et itérer dans tous les SCR connus : Pour chaque SCR, il va l’affecter à la couche inconnue, et vérifier si les données ainsi renseignées se trouvent à l’intérieur ou à proximité de la zone cible définie par la couche de référence. Si c’est le cas, le traitement ajoute le SCR à la liste des candidats potentiels. Le traitement continue à parcourir tous les SCR jusqu’à ce que la liste complète des candidats ait été créée.

L’algorithme génère donc une liste restreinte de SCR candidats qui, s’ils étaient utilisés pour la couche inconnue, la placeraient près de la zone cible. A vous de tester manuellement chaque SCR candidat à l’aide de QGis, en définissant le SCR pour la couche et en le comparant à la couche de référence avec un SCR connu, en fonction de l’alignement et de la compatibilité avec la couche de référence.

Il est à noter que l’algorithme peut produire plusieurs candidats ou aucun candidat, selon la précision de l’entrée de la zone cible et l’étendue de la couche. Plus grande est l’étendue de la couche de référence, plus il y aura de SCR candidats. Il est également possible que le bon SCR ne figure pas parmi les candidats présélectionnés, il est donc crucial de vérifier les résultats et de combiner cette approche avec d’autres méthodes pour assurer une évaluation complète.

Pour poursuivre notre exemple, nous ouvrons l’outil Trouver une projection:

Vous pouvez définir la Zone cible pour la couche en cliquant sur l’icône à droite du champ. Dans ce cas la zone sera celle affichée sur la fenêtre cartographique de QGis. Vous pouvez aussi ouvrir la liste déroulante qui vous propose d’autres méthodes de définition. Cliquez sur Exécuter et patientez car QGis va projeter votre couche inconnue dans chacun des SCR qu’il connaît.

Une fois terminé, il affiche la liste restreinte dans l’onglet journal, et il aura ajouté une table dans le panneau des couches avec un enregistrement par SCR candidat.

Pour vous compliquer la tâche, et là QGis n’y est pour rien, le même système de coordonnées va avoir des codes différents selon l’organisme de référence. Ce sont les 4 lettres en début de chaque code. EPSG correspond à European Petroleum Survey Group, premier organisme à normaliser les systèmes de coordonnées de par le monde, IGNF correspond à l’IGN France, et ESRI à ESRI évidemment. Pour le même système chacun à un code différent.

Si vous regardez les définitions de tous les SCR de cet exemple:

EPSG:27573
EPSG:27583
IGNF:LAMB3C
ESRI:102583
EPSG:7403
EPSG:7413
EPSG:7422
ESRI:102587
IGNF:NTFLAMB3C
IGNF:NTFLAMB3C.IGN69

Ils correspondent tous au même système de projection Lambert 3 zone, utilisé jusqu’aux années 2000 pour le Sud de la France. Quel qu’il soit le SCR utilisé, le résultat sera le même:

Mais ce n’est pas toujours aussi simple. Si vous avez des données non projetées, en coordonnées géographiques Latitude-Longitude vous aurez le problème du choix de DATUM, en termes moins compliqués, de choix du centre de la terre de référence. Les mêmes coordonnées Lat/long, selon le DATUM utilisé peuvent se retrouver dans des positions distantes de 350 mètres…

Bref, si une fois suivies toutes les étapes précédentes vous n’avez pas réussi à résoudre votre problème, il vous reste l’option de me contacter à travers le formulaire de contact. Essayez de décrire au mieux votre problème, ça nous fera gagner du temps 😉.

Si cet article vous a intéressé et que vous pensez qu'il pourrait bénéficier à d'autres personnes, n'hésitez pas à le partager sur vos réseaux sociaux en utilisant les boutons ci-dessous. Votre partage est apprécié !

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *