Identifizieren des Koordinatensystems von Shapefile, wenn unbekannt?

Ich habe ein Shapefile, aber sein Koordinatensystem ist unbekannt und es gibt keine * .prj-Datei.

Wie kann ich es jetzt identifizieren?

Es gibt immer die "Brute Force" -Methode:

1. Nehmen Sie eine Ebene mit einem bekannten Koordinatensystem, die mit Ihrer unbekannten Ebene überlagert werden soll.

2. Machen Sie jetzt einige fundierte Vermutungen darüber, welche Projektion die unbekannte Ebene sein könnte. (UTM, Plate Carree usw.). Projizieren Sie Ihren bekannten Koordinatensystem-Layer in jede Projektion, bis Sie einen finden, der dem unbekannten Layer so gut wie möglich entspricht.

Ich fördere Kommentare von Mapperz und Brandon Copeland, indem ich eine Antwort hinzufüge, die ihre Technik verwendet.

1. Wenn Sie über ArcGIS Desktop verfügen, fügen Sie in ArcMap einige Referenzdaten hinzu, die denselben Bereich abdecken. Diese Referenzdaten müssen eine gültige Definition des Koordinatensystems (Projektion) haben.
2. Fügen Sie die Daten in einem unbekannten Koordinatensystem hinzu.
3. Erforschen Sie, welche möglichen Koordinatensysteme in diesem Bereich verwendet werden. Ein zu überprüfender Ort ist das EPSG-Register .
4. Stellen Sie das Koordinatensystem des Datenrahmens auf die Möglichkeiten ein. Wenn die Daten mit bekannten Koordinaten mit den Daten ohne bekannte Koordinaten übereinstimmen, haben Sie wahrscheinlich das Koordinatensystem identifiziert.

Diese Technik nutzt die Tatsache, dass ArcMap Daten in ein neues Koordinatensystem im Speicher konvertieren kann. Die Daten ohne bekannte Koordinaten können nicht konvertiert werden und werden nur angezeigt. Indem wir die bekannten Daten in den Speicher projizieren, können wir schnell sehen, welche Koordinaten sie in den verschiedenen möglichen Koordinatensystemen haben.

Da diese Frage nie alt wird, habe ich eine Site mit der Brute-Force-Methode erstellt . Wenn Sie ein gezipptes shp + shx auf die Karte ziehen, wird es in jedem in PostGIS verfügbaren Koordinatensystem abgebildet. Vorausgesetzt, Sie wissen, wie "richtig" aussieht, können Sie in diesen Bereich zoomen und auf das Polygon klicken, um die .prj-Datei von epsg.io abzurufen.

Es gibt zwei großartige Links von Esri, die diesbezüglich ins Detail gehen:

• Identifizieren Sie mit ArcMap ein unbekanntes projiziertes Koordinatensystem
• Identifizieren Sie den Raumbezug, die Projektion oder das Koordinatensystem der Daten

Bei meinem letzten Job hatte ich einen Layer mit Geologiedaten (Polygon) mit dem Namen "FSU_Geol.shp". Mein Chef gab es mir und bat mich, einige Dinge herauszufinden. Zuerst wurde ihm dieses Shapefile vom Kunden ausgehändigt, und es gab keine .prj-Datei. Deshalb wollte er, dass ich das herausfinde. Er wollte auch, dass ich herausfinde, was die Kategorien der Geologie sind. Ich könnte weitermachen, aber lassen Sie uns auf den Punkt kommen ... Ich habe gegoogelt und bin HIER gelandet .

Das Akronym in der Namenskonvention stand für "Oberflächengeologie der ehemaligen Sowjetunion", und die Google-Suche führte mich "direkt" zur Quelle (USGS). Alles, was ich jemals über dieses Shapefile wissen musste, war der oberste Link, den ich getroffen habe. Ich sage nicht, dass Google alles und jedes finden kann , aber ich war frisch von der Universität und habe nur einen "Schuss in die Dunkelheit" gemacht und mir das Feedback angesehen, das ich erhalten habe!

In einem anderen Fall hatte jemand im Büro mithilfe einer Batch-Funktion eine Reihe von Shapefiles heruntergeladen. Ich erinnere mich momentan nicht an den Namen des Softwarepakets, aber die .prj-Dateien fehlten. Ich habe einfach die Metadaten in ArcCatalog aufgerufen und dort meine Antwort gefunden. Auch hier zitiere mich nicht, aber ich denke, es gab eine URL zu der Quelle darin.

Seit diesem Tag "beginne" ich mit einer Google-Suche, wenn die Metadaten keine Hinweise enthalten!

Frage den Produzenten.

Wenn Sie ein allgemeines CRS für Ihre geografische Region kennen, können Sie einige davon ausprobieren. Aber fragen ist besser.

Ich liebe dieses Web-Tool wirklich: http://projfinder.com/ . Überprüfen Sie Ihre Datei und suchen Sie einige Koordinaten. Zoomen Sie auf ungefähr dieselbe Stelle auf dem Globus und lassen Sie es raten.

Kein Tool (ich kenne kein Tool, mit dem Sie das tun könnten ), aber lesen Sie die Antwort von @ mkennedy auf Wie würde ich diesen Punkt in WKID 4326 konvertieren? . Sie erklärt, wie sie zum richtigen Raumbezug gekommen ist. SpatialReference.org und Geduld werden Ihre Freunde sein.

Außerdem bietet Esri eine Anleitung zum Erraten eines Koordinatensystems (obwohl ich die Methode von mkennedy bevorzuge, wenn Sie etwas mehr über die Daten wissen).

Die beste Antwort, die ich auf diese Frage gefunden habe, ist nicht technisch: Finden Sie heraus, woher Ihre Daten stammen. Agenturen und Organisationen neigen dazu, Projektionen zu verwenden. Wissen Sie, dass es von Ihrem Zustand PUNKT kam? Schauen Sie sich den Rest ihrer Daten an und sehen Sie, was sie Ihnen sagen. Weißt du nicht, woher es kam? Eine fundierte Vermutung wird Sie wahrscheinlich auf den richtigen Weg bringen.

Zumindest macht es das Problem mit brachialer Gewalt ein wenig mehr machbar!

Ernsthaft veraltet zu diesem Zeitpunkt, aber Werner Flacke und Birgit Klaus haben 2007 Find Projection auf ArcScripts veröffentlicht. Ich glaube, der Quellcode ist leider nicht vorhanden. Es ist VBA-basiert und daher nur in ArcGIS Desktop Version 9.2 und möglicherweise Version 9.3 verwendbar. Es enthält zwei Shapefiles mit den interessierenden Bereichen aus dem EPSG Geodetic Parameter Dataset, mit denen die Möglichkeiten eingegrenzt werden könnten.

Blue Marble Geographic Calculator und Geographic Transformer verfügen über Tools zur Wiederherstellung des Koordinatensystems.

Probieren Sie das im Lieferumfang von GDAL enthaltene Programm ogrinfo aus.

Siehe So greifen Sie mit OGR auf Shapefile-Metadaten zu

Ogrinfo gibt Ihnen möglicherweise keine Projektionsinformationen ohne eine .prj-Datei, ist aber dennoch ein nützliches Tool, mit dem Sie die Liste potenzieller Projektionen untersuchen können.

Zum Beispiel:

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

Diese Antwort legt nahe, dass das Shapefile ein auf Metern und nicht auf Grad basierendes Referenzsystem verwendet.

Vorausgesetzt, Sie wissen, woher die Daten stammen, haben Sie jetzt eine kürzere Liste möglicher Projektionen.

Andere bewährte Methoden könnten sein:

• Fragen Sie zunächst die Personen, die Ihnen die Daten zur Verfügung gestellt haben, auch wenn sie nicht genau wissen, ob sie eine Liste möglicher Projektionen haben sollten, die in erster Linie in Betracht gezogen werden könnten.
• Überprüfen Sie, ob Metadaten zugeordnet sind.

Es ist wahrscheinlich besser, von einer Position der Bekannten aus zu beginnen, als zu versuchen, Gewalt von Grund auf neu anzuwenden.

Um die anderen Antworten in diesem Beitrag zu ergänzen, würde ich Folgendes hinzufügen:

Gewusst wie: Identifizieren eines unbekannten projizierten Koordinatensystems mithilfe von ArcMap

Stellen Sie eine spezifische Prüfung für NAD1927 bereit

Wenn die Koordinaten in Dezimalgrad angegeben sind, z. B. zwischen -180 und +180 Längengrad und -90 und +90 Breitengrad, geben Sie das für die Daten verwendete geografische Koordinatensystem (Datum) an. Vor Version 9.2 weist ArcMap den Daten standardmäßig GCS_Assumed_Geographic_1 zu. Dies platziert die Daten auf dem NAD_1927-Datum

In Ergänzung

Wenn sich die Daten in den USA befinden und eine Ausdehnung aufweisen, in der die Koordinaten links von der Dezimalstelle aus 6, 7 oder 8 Ziffern bestehen, werden sie wahrscheinlich in das State Plane- oder UTM-Koordinatensystem projiziert.

FAQ: Grundlagen der Projektion: Was der GIS-Fachmann wissen muss

Diese Seite enthält mehrere Überprüfungen zur Bestimmung von Projektionen, die wichtig sind, um zu versuchen, ein unbekanntes Koordinatensystem zu identifizieren.

Die folgenden Konzepte sind für das Verständnis der Verwendung von Kartenprojektionen in ArcGIS von grundlegender Bedeutung. Bitte beachten Sie jedoch, dass das Thema Projektionen sehr weit gefasst ist und dieser Artikel nur einige wichtige Themen ansprechen kann.

1. Koordinatensysteme, auch Kartenprojektionen genannt, sind beliebige Bezeichnungen für räumliche Daten. Ihr Zweck ist es, eine gemeinsame Grundlage für die Kommunikation über einen bestimmten Ort oder ein bestimmtes Gebiet auf der Erdoberfläche zu schaffen. Das kritischste Problem beim Umgang mit Kartenprojektionen ist die Kenntnis der Projektion und die Zuordnung der richtigen Koordinatensysteminformationen zu einem Datensatz.

2. Als die ersten Kartenprojektionen erstellt wurden, wurde fälschlicherweise angenommen, dass die Erde flach war. Später wurde die Annahme revidiert und die Erde als perfekte Kugel angenommen. Im 18. Jahrhundert begannen die Menschen zu begreifen, dass die Erde nicht perfekt rund war. Dies war der Beginn des Konzepts des kartographischen Sphäroids.

3. Um Orte auf der Erdoberfläche genauer darzustellen, untersuchten die Kartenhersteller die Form der Erde (Geodäsie) und entwickelten das Konzept des Sphäroids. Dann wurden geografische Koordinatensysteme (GCS) entwickelt, die ein Datum, Maßeinheiten und einen Nullmeridian enthalten. Ein Datum verbindet einen Sphäroid mit einem bestimmten Teil der Erdoberfläche. Neuere Daten sind so konzipiert, dass sie gut auf die gesamte Erdoberfläche passen.

4. Die in Nordamerika am häufigsten verwendeten Bezugspunkte sind:

• NAD 1927 (North American Datum 1927) unter Verwendung des Sphäroids Clarke 1866
• NAD 1983 (North American Datum 1983) unter Verwendung des Sphäroids GRS 1980
• WGS 1984 (World Geodetic Survey 1984) unter Verwendung des Sphäroids WGS 1984

Neuere Sphäroide basieren auf Satellitenmessungen und sind genauer als die von Clarke im Jahr 1866 entwickelten. Die Begriffe "geografisches Koordinatensystem" und "Datum" werden austauschbar verwendet und ein Nullmeridian.

5. Die Datenkoordinaten ändern sich je nach Bezugspunkt und Sphäroid, auf denen diese Koordinaten basieren, auch wenn sie dieselbe Kartenprojektion und dieselben Parameter verwenden.

Die folgenden geografischen Koordinaten beziehen sich beispielsweise auf einen einzelnen Punkt in der Stadt Bellingham, Washington, und verwenden drei verschiedene Bezugspunkte:

Code: DATUM X-Koordinate Y-Koordinate NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.7438798534299

6. Ein Prinzip einer guten Datenverwaltung besteht darin, die Projektionsparameter von der Datenquelle zu erhalten, die die Daten bereitstellt. Machen Sie keine fundierten Vermutungen über die Projektion von Daten, da dies zu einer ungenauen GIS-Datenbank führen wird. Die notwendigen Parameter sind die folgenden:

• Projektion
• Maßeinheiten
• ZONE (für UTM)
• FIPS-Zone (für State Plane)
• Datum

Je nach Projektion können andere Parameter erforderlich sein. Beispielsweise erfordern Albers- und Lambert-Projektionen die folgenden Parameter:

• 1. Standardparallele in Grad, Minuten und Sekunden (DMS)
• 2. Standardparallele (DMS)
• Mittelmeridian (DMS)
• Breite des Projektionsursprungs (DMS)
• Falscher Osten und Maßeinheiten
• Falscher Norden und Maßeinheiten
• X-Verschiebung und Maßeinheiten
• Y-Verschiebung und Maßeinheiten

Ich hatte das gleiche Problem, als ich an meinem Projekt (DEMs) arbeitete. Eine der DEM-Koordinaten war unbekannt und kann nicht projiziert werden . Daher habe ich die Georeferenzierung mit dem Georeferenzierungs-Tool in ArcMap durchgeführt . Nach der Georeferenzierung erhalten Sie Ihre prj.file.

Schauen Sie sich diese Website an:

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

Ich bin der Entwickler dieses kostenlosen Tools. Vielleicht kann dieses Tool Ihnen helfen, die richtige Projektion Ihres Shapefiles zu finden und zu definieren. Es basiert auf der Idee von http://projfinder.com/ von Aaron Racicot.

Sie können diese speziellen Lösungen auch ausprobieren für:

Österreich (deutsch): https://www.data.gv.at/anwendungen/checkaustrianprojection/

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

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

Sie können zwischen geografischen und kartografischen Koordinaten unterscheiden, die ersten Breiten- und Längengrade und die zweiten Nord- und Ostkoordinaten werden in Metern angegeben. Wenn die Datei eine Form hat, können Sie die Koordinaten der Geometrie mit shapelib-1.3.0 ausgeben und die Werte dann anzeigen. ..

Du kannst es mit Fiona machen.

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

Sie erhalten eine Antwort als

    {'init': 'epsg:4326'}