Identifier le système de coordonnées du fichier de forme lorsqu'il est inconnu?

J'ai un fichier de forme mais son système de coordonnées est inconnu et il n'y a pas de fichier * .prj.

Comment puis-je l'identifier maintenant?

Il y a toujours la méthode "force brute":

1. Prenez un calque avec un système de coordonnées connu qui est supposé superposer avec votre calque inconnu.

2. Maintenant, faites des suppositions éclairées sur la projection de la couche inconnue. (UTM, Plate Carree, etc.). Projetez votre couche de système de coordonnées connue dans chaque projection jusqu'à ce que vous en trouviez une qui corresponde autant que possible à la couche inconnue.

Je fais la promotion des commentaires de Mapperz et de Brandon Copeland en ajoutant une réponse qui utilise leur technique.

1. Si vous avez ArcGIS Desktop, ajoutez des données de référence couvrant la même zone à ArcMap. Ces données de référence doivent avoir une définition valide du système de coordonnées (projection).
2. Ajoutez les données dans un système de coordonnées inconnu.
3. Recherchez les systèmes de coordonnées possibles utilisés dans ce domaine. Le registre EPSG est un lieu à vérifier .
4. Définissez le système de coordonnées du bloc de données en fonction des possibilités. Si les données avec des coordonnées connues sont alignées avec celles sans les coordonnées connues, vous avez probablement identifié le système de coordonnées.

Cette technique tire parti du fait qu'ArcMap peut convertir des données en un nouveau système de coordonnées en mémoire. Les données sans coordonnées connues ne peuvent pas être converties, elles sont simplement affichées. En projetant les données connues en mémoire, nous pouvons voir rapidement quelles sont ses coordonnées dans les différents systèmes de coordonnées possibles.

Comme cette question ne vieillit jamais, j'ai construit un site utilisant la méthode Brute Force . Si vous faites glisser un shp + shx zippé sur la carte, il le mappera dans tous les systèmes de coordonnées disponibles dans PostGIS. En supposant que vous sachiez à quoi "correct" ressemble, vous pouvez zoomer sur cette zone et cliquer sur le polygone pour obtenir le fichier .prj à partir de epsg.io.

Esri propose deux liens intéressants:

• Identifier un système de coordonnées projetées inconnu à l'aide d'ArcMap
• Identifier la référence spatiale, la projection ou le système de coordonnées de données

Lors de mon dernier emploi, j'avais une couche de données géologiques (polygone) nommée "FSU_Geol.shp". Mon patron me l'a donné et m'a demandé de découvrir un certain nombre de choses. Tout d’abord, le client lui a remis ce fichier de formes, et il n’existait pas de fichier .prj; il voulait donc que je le découvre. Il voulait également que je découvre quelles étaient les catégories de géologie. Je pourrais continuer encore et encore, mais passons à la chasse ... Je l'ai googlé, et je me suis retrouvé ICI .

L'acronyme utilisé dans la convention de dénomination signifiait "géologie de surface de l'ancienne Union soviétique" et la recherche sur Google m'a conduit "directement" à la source (USGS). Tout ce que je pouvais avoir besoin de savoir sur ce fichier de formes était au sommet du lien que je découvrais. Je ne dis pas que Google peut trouver n'importe quoi et tout , mais je venais de sortir de l'université et je prenais juste un "coup dans le noir", et regardez les retours que j'ai obtenus!

Dans un autre cas, une personne du bureau a téléchargé un certain nombre de fichiers de formes à l’aide d’une fonction de traitement par lots. Je ne me souviens pas du nom du progiciel pour le moment, mais les fichiers .prj étaient manquants. Je suis simplement allé dans les métadonnées d'ArcCatalog et j'ai effectivement trouvé ma réponse. Encore une fois, ne me citez pas, mais je pense qu'il y avait une URL vers la source.

Depuis ce jour, je "commence" avec une recherche Google s'il n'y a pas d'indice dans les métadonnées!

Demandez au producteur.

Si vous connaissez des SIR communs à votre région géographique, vous pouvez en essayer quelques-uns. Mais demander, c'est mieux.

J'aime vraiment cet outil Web: http://projfinder.com/ . Vérifiez votre fichier et recherchez des coordonnées. Zoomez à peu près au même endroit sur le globe et laissez-le deviner.

Ce n'est pas un outil (je n'en connais aucun qui vous permettrait de le faire), mais consultez la réponse de @ mkennedy à Comment convertir ce point en WKID 4326? . Elle explique comment elle est arrivée à la bonne référence spatiale. SpatialReference.org et la patience seront vos amis.

De plus, Esri fournit un guide sur la manière de deviner un système de coordonnées (bien que je préfère la méthode de mkennedy si vous en savez un peu plus sur les données).

La meilleure réponse que j'ai trouvée à cette question est non technique: découvrez d'où proviennent vos données. Les agences et organisations ont tendance à être cohérentes avec leur utilisation des projections. Sachez que cela vient de votre état DOT? Regardez le reste de leurs données et voyez ce que cela vous dit. Je ne sais pas d'où ça vient? Une supposition éclairée est tout aussi susceptible de vous envoyer dans la bonne voie.

Au moins, il est un peu plus facile de résoudre le problème avec la force brute!

Sérieusement obsolète à ce stade, mais Werner Flacke et Birgit Klaus ont publié Trouver Projection sur ArcScripts en 2007. Je ne pense pas que le code source existe, malheureusement. Il est basé sur VBA et ne peut donc être utilisé que dans ArcGIS Desktop version 9.2 et éventuellement 9.3. Il inclut deux fichiers de formes avec les zones d'intérêt du jeu de données de paramètres géodésiques EPSG qui pourraient être utilisés pour réduire les possibilités.

Blue Marble Geographic Calculator et Geographic Transformer ont des outils de récupération du système de coordonnées.

Essayez le programme ogrinfo fourni avec GDAL.

Voir Comment accéder aux métadonnées Shapefile en utilisant OGR?

Donc, ogrinfo pourrait ne pas vous donner les informations de projection sans un fichier .prj, mais cela reste un outil utile pour vous aider à explorer la liste des projections potentielles.

Par exemple:

Geometry: Polygon
Feature Count: 269 
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000) 

Cette réponse me suggère que le fichier de formes utilise un système de référence basé sur des mètres et non sur des degrés.

En supposant que vous sachiez d'où proviennent approximativement les données, vous disposez maintenant d'une liste plus courte de projections possibles.

Les autres meilleures pratiques pourraient être:

• Demandez aux personnes qui vous ont fourni les données au départ, même si elles ne savent pas exactement, elles devraient avoir une liste de projections possibles qui pourraient être examinées en premier lieu.
• Regardez s'il y a des métadonnées associées.

Il est probablement préférable de commencer par une position connue, plutôt que d’essayer de forcer de force.

Pour ajouter aux autres réponses dans ce post, j'ajouterais ce qui suit:

Comment: identifier un système de coordonnées projetées inconnu à l'aide d'ArcMap

Fournir une vérification spécifique sur NAD1927

Si les coordonnées sont exprimées en degrés décimaux, par exemple entre les longitude -180 et +180 et les latitudes -90 et +90, identifiez le système de coordonnées géographiques (datum) utilisé pour les données. Avant la version 9.2, ArcMap affectait par défaut GCS_Assumed_Geographic_1 aux données. Cela place les données sur le datum NAD_1927

en outre

Si les données se trouvent aux États-Unis et indiquent une étendue dans laquelle les coordonnées à gauche de la virgule sont 6, 7 ou 8 chiffres, les données sont probablement projetées dans les systèmes de coordonnées State Plane ou UTM.

FAQ: Notions de base sur la projection: ce que le professionnel des SIG doit savoir

Cette page fournit plusieurs vérifications sur la détermination des projections, ce qui est important pour vous aider à identifier un système de coordonnées inconnu.

Les concepts suivants sont fondamentaux pour comprendre l'utilisation des projections cartographiques dans ArcGIS. Veuillez noter cependant que le sujet des projections est extrêmement large et que cet article ne peut que toucher à quelques sujets importants.

1. Les systèmes de coordonnées, également appelés projections cartographiques, sont des désignations arbitraires de données spatiales. Leur objectif est de fournir une base commune pour la communication sur un lieu ou une zone particulière de la surface de la Terre. Le problème le plus important dans le traitement des projections cartographiques est de savoir quelle est la projection et d’avoir les informations correctes du système de coordonnées associées à un jeu de données.

2. Lors de la conception des premières projections cartographiques, il était supposé, à tort, que la terre était plate. Plus tard, l’hypothèse fut révisée et la terre supposée être une sphère parfaite. Au 18ème siècle, les gens ont commencé à réaliser que la terre n'était pas parfaitement ronde. Ce fut le début du concept de sphéroïde cartographique.

3. Pour représenter plus précisément les emplacements sur la surface de la Terre, les cartographes ont étudié la forme de la Terre (géodésie) et créé le concept de sphéroïde. Des systèmes de coordonnées géographiques (GCS) ont ensuite été mis au point. Ils comprennent un datum, des unités de mesure et un premier méridien. Une donnée relie un sphéroïde à une partie particulière de la surface de la Terre. Les systèmes de référence récents sont conçus pour s'adapter à toute la surface de la Terre.

4. Les systèmes de référence les plus couramment utilisés en Amérique du Nord sont:

• NAD 1927 (datum nord-américain 1927) utilisant le sphéroïde Clarke 1866
• NAD 1983 (datum nord-américain 1983) utilisant le sphéroïde GRS 1980
• WGS 1984 (World Geodetic Survey 1984) utilisant le sphéroïde WGS 1984

Les sphéroïdes plus récents sont développés à partir de mesures satellitaires et sont plus précis que ceux développés par Clarke en 1866. Les termes «système de coordonnées géographiques» et «datum» sont utilisés de manière interchangeable. Toutefois, comme indiqué ci-dessus, un SCG comprend un datum, un sphéroïde, des unités de mesure. et un premier méridien.

5. Les coordonnées des données changent en fonction du datum et du sphéroïde sur lesquels elles sont basées, même s'ils utilisent la même projection cartographique et les mêmes paramètres.

Par exemple, les coordonnées géographiques ci-dessous correspondent à un seul point situé dans la ville de Bellingham, dans l’État de Washington, en utilisant 3 points de référence différents:

Code: DATUM Coordonnée X Coordonnée Y NAD_1927 - 1242.466903686523 48.7440490722656 NAD_1983 - 122.46818353793 48.7438798543649 WGS_1984 - 122.46818353737 48.74387985432999

6. Un principe de bonne gestion des données consiste à obtenir les paramètres de projection de la source de données fournissant les données. Ne faites pas de conjecture éclairée sur la projection des données, car une base de données SIG inexacte en résultera. Les paramètres nécessaires sont les suivants:

• Projection
• Unités de mesure
• ZONE (pour UTM)
• Zone FIPS (pour State Plane)
• Datum

D'autres paramètres peuvent être nécessaires en fonction de la projection. Par exemple, les projections Albers et Lambert nécessitent les paramètres suivants:

• 1er parallèle standard, en degrés, minutes et secondes (DMS)
• 2ème parallèle parallèle (DMS)
• Méridien central (DMS)
• Latitude de l'origine des projections (DMS)
• Fausse est et unités de mesure
• Faux nord et unités de mesure
• X-shift et unités de mesure
• Y-shift et unités de mesure

J'ai eu le même problème pendant que je travaillais sur mon projet (DEM). L'un des MNT était des coordonnées inconnues et ne pouvait pas être projeté. J'ai donc utilisé le "géoréférencement" à l'aide de l'outil Géoréférencement dans ArcMap . Après le géoréférencement, vous aurez votre fichier prj.file.

Regardez ce site:

http://www.egger-gis.at/shapefile-projectionfinder/

Je suis le développeur de cet outil gratuit. Peut-être que cet outil peut vous aider à trouver et à définir la projection correcte de votre fichier de formes. Il s’appuie sur l’idée de http://projfinder.com/ d’Aaron Racicot.

Vous pouvez également essayer ces solutions spéciales pour:

Autriche (allemand): https://www.data.gv.at/anwendungen/checkaustrianprojection/

Australie: https://maegger.github.io/map_australia.html

UTM - Zones: https://maegger.github.io/map_utm.html

Vous pouvez faire la distinction entre les coordonnées géographiques et cartographiques, les premières lat / lon sont en degré, les secondes nord / est en mètres. Si le fichier est en forme, vous pouvez vider les coordonnées de la géométrie avec shapelib-1.3.0 et voir les valeurs. ..

Vous pouvez le faire avec Fiona.

    import fiona as f
    a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
    print(a.crs)

Vous obtiendrez réponse comme

    {'init': 'epsg:4326'}