2017-03-09 3 views
4

Ich handle jetzt shp Dateien und ich stieß auf Probleme mit den Projektionen. Lassen Sie mich Ihnen meinen Code unten geben.Geopandas Projektionsprobleme

Import Pandas als pd Import geopandas als gpd aus geopandas importieren GeoSeries, GeoDataFrame import os

Aelly = gpd.read_file(r'C:\Users\Hyun Mo\Downloads\조인 (1)\after_join.shp', encoding = 'utf-8') 
base_map = gpd.read_file(r'C:\Users\Hyun Mo\Downloads\11000 (3)\TL_SCCO_SIG.shp', encoding = 'ANSI') 
Aelly_to_crs = Aelly.to_crs(base_map.crs) 
Aelly_to_crs.plot(ax=base_map.plot()) 

Und hier ist mein Datenaufbau

print(base_map.head()) 

data construction

print(Aelly.head()) 

data construction

Wenn ich print(base_map.crs), print(Aelly_to_crs.crs) ausgeführt habe, habe ich die Ergebnisse wie folgt erhalten.

results

Aelly_to_crs.plot(ax=base_map.plot()) 

data construction

Das obige Bild ist das Ergebnis der Ausführung Aelly.plot(ax=base_map.plot()) Und Sie können sehen, dass die beiden Bilder einander nicht entsprechen. Wie kann ich dieses Problem lösen?

----------- bearbeiten

Meine gewünschte Ausgabe Bild unten ist.

enter image description here

Hier sind meine Daten Links: http://blog.naver.com/khm2963/220929301892

unten Bilder sind Verfahren für meine flie Download

enter image description here enter image description here

+0

Ich verstehe das Problem nicht. Was ist die erwartete Ausgabe? Was könnte das Problem sein? Was hast du probiert?Die Leute werden wahrscheinlich Ihre Dateien nicht herunterladen, deshalb ist es besser, wenn Sie uns ein minimales Beispiel mit Daten geben, die wir kopieren/einfügen können. – cd98

+0

Oh ~ Entschuldigung ~ Ich bin spät dran. Und ich nehme an, dass Leute meine Daten von meinem Link bekommen. Aber es funktioniert nicht. Ich ändere meine Datenverbindung und füge weitere Informationen zu meiner Frage hinzu. Entschuldigung –

+1

Bitte erwägen Sie, ein paar (oder alle) Polygone auf http://geojson.io/ oder einem ähnlichen Dienst zu teilen, damit wir keine Dateien herunterladen müssen! – cd98

Antwort

0

Aus den Daten, die Sie gedruckt haben es aussieht Alles funktioniert wie es soll! Die Koordinaten zwischen den Shapefiles sind sehr unterschiedlich, aber crs ist das selbe, also macht die Handlung total Sinn.

GeoPandas kann Ihnen nicht sagen, ob die Daten und die Nummer in einer realen Welt sinnvoll sind. Du hast ihm zwei Shapefiles mit gut definierter Projektion (EPSG: 32652) und mit fest codierten Koordinaten gegeben und GeoPandas ist damit zufrieden.

Wenn Sie wissen, dass in der Realität beide Shapefiles denselben Bereich darstellen, dann müssen Sie erkennen, dass die Datenquelle irgendwie korrupt ist. Ich denke, dass eines der Shapefiles versehentlich eine andere crs-Definition als Metadaten erhalten hat (stelle es dir zum Beispiel als eine falsche Textkodierung vor).

Der einfachste Weg, das herauszufinden und zu korrigieren, ist die Verwendung von ArcGIS oder QGIS, wo Sie mit verschiedenen Projektionen spielen können, um herauszufinden, wie die richtige Projektion war. Dann können Sie das Shapefile mit neuen Projektionsmetadaten speichern und der Rest funktioniert sofort.

Verwandte Themen