Identifizieren des Koordinatensystems von Shapefile, wenn unbekannt?

Antworten:

49

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.

brenth
quelle
1
das ist auch meine methode. Ich benutze zwei verschiedene Software, um es zu "beweisen". Sie können diese Frage auch für einige grundlegende Anweisungen zu Projektionen verwenden. gis.stackexchange.com/questions/2769/…
Brad Nesom
6
Sie können Schritt 2 ein wenig vereinfachen, indem Sie einfach die Datenrahmenprojektion (ArcMap) ändern und die bekannte Ebene sofort neu projizieren.
Brandon Copeland
13
Beachten Sie, dass Sie die Daten nicht projizieren sollten. Sie sollten die Projektion definieren. Durch die Projektion der Daten werden die Koordinaten geändert. Durch das Definieren der Projektion wird diesen Daten lediglich eine Projektion zugewiesen, und der Software wird mitgeteilt, in welcher Projektion sich die Daten bereits befinden. Dies sind zwei verschiedene Tools in ArcMap: "Define Projection" (Projektion definieren).
Sean
Ich ziehe es vor, die Koordinatenwerte zu lesen. Normalerweise kann ich feststellen, ob es sich um Dezimalgrade, Meter oder Fuß handelt (schwerer zu unterscheiden), und dann zum besten System übergehen. Beginnen Sie mit WGS 84, wenn es dezimal erscheint, UTM, wenn Meter, Staatsebene, wenn Fuß. Ich hatte einmal eine, die niemand reparieren konnte (für Gravitationsanomalien). Es war ein US-amerikanisches Staatsflugzeug mit einer Längeneinheit von ...... Kilometern. Ich konnte an der Tatsache erkennen, dass es nah war, dass ich das richtige System hatte, aber der Zustand war winzig, ahh, Kilometer waren die Lineareinheiten.
Wenn Sie nicht wissen, nur GIS
33

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.

mkennedy
quelle
4
Randnotiz: In QGIS wird die Neuprojektion im Speicher als "on-the-fly" bezeichnet.
Underdunkel
3
In ArcGIS heißt es auch "on-the-fly"! Ich finde immer wieder Leute, die diesen Ausdruck nicht verstehen, also habe ich versucht, stattdessen 'in-memory' zu verwenden.
Mkennedy
Ist dieser Link noch gültig? Scheint ein bisschen verwickelter zu sein,
GISKid
1
@ GISKid Du hast recht. Ich glaube, die Website epsg-registry.org hat das Datum angegeben, an dem ich diese Antwort geschrieben habe. Ich werde aktualisieren.
Mkennedy
1
Eine weitere großartige Ressource: epsg.io
GISKid
29

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.

jmk
quelle
15
hervorragende Idee, aber ich denke, die Website ist nicht wirklich funktioniert
Fábio Dias
22

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

Simon
quelle
1
downloads2.esri.com/MappingCenter2007/arcGISResources/more/… ist ein kürzlich veröffentlichter Beitrag im Ressourcenzentrum (Mapping-Bereich), der den Vorgang durchläuft, ähnlich dem Vorschlag von
brenth
22

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!

Dano
quelle
18

Frage den Produzenten.

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

Underdunkel
quelle
Das einzige Problem ist, dass ich ihn nicht fragen kann ... Dies war auch der erste Gedanke, den ich hatte.
com
Dann musst du brutale Gewalt anwenden, wie bereits erwähnt.
Underdunkel
18

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.

johanvdw
quelle
Dies ist tatsächlich das beste Tool, das ich gefunden habe. Vermutungen in ArcMap sind auf lange Sicht nicht wirklich nützlich oder wenn Sie nicht wissen, wie sie aussehen sollen. Kann nicht mit Homebrew-Koordinatensystemen und dem Fehlen eines vollständig helfen, aber wer kann!
Martin
Schade, dass es nicht mehr funktioniert - weiß jemand, an wen er sich wenden kann?
Sventechie
16

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).

Sasa Ivetic
quelle
13

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!

Jason Baker
quelle
11

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.

mkennedy
quelle
Danke für den Tipp Melita. Der Code ist dort (* .frm). Die Lizenz ist gemeinfrei. PrjFinder ist eine Ergänzung zum Buch "Arbeiten mit Projektionen und Datumstransformationen in ArcGIS" (2006).
Matt Wilkie
4

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.

nmtoken
quelle
1
Dies beantwortet die gestellte Frage nicht.
Dan C
1
Wenn keine prj-Datei vorhanden ist, gibt hen ogrinfo keine Ergebnisse für die Koordinatensysteminformationen zurück.
Devdatta Tengshe
4

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.

  1. 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

  1. 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

Whyzar
quelle
Was soll ich tun, wenn ich eine Gitterebene habe, die nur ein unbekanntes Koordinatensystem anzeigt und der ich nichts zuweisen kann? Vielen Dank.
AndrewLebron
2

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.

ndr
quelle
0

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

Manfred Egger
quelle
0

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. ..

Natan
quelle
0

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'}
Sandip Nepal
quelle
Sind Sie sicher, dass dies ohne eine PRJ-Datei funktioniert? wie wird Fiona es wissen?
JR
Nein, ohne .prj-Datei wird es nicht funktionieren
Sandip Nepal