Mit den in "Statistics on Maps (Statistical GIS)" der japanischen Regierungsstatistik verwendeten GIS-Daten können dieselben Daten in verschiedenen Formaten heruntergeladen werden. Ich habe versucht herauszufinden, was anders ist.
shp ist ein altes Binärformat und GML ist ein kürzlich entwickeltes XML-Format. Verwenden Sie beim Lesen mit Python geopandas (intern fiona) shp. Es scheint, dass die von "fiona" intern verwendete "GDAL" die erforderlichen XML-Spezifikationen nicht vollständig unterstützt.
Danach war es aufgrund des Unterschieds in der Projektion, als ich mir die Definition von "proj4" ansah, so ([Erklärung] von "Geopandas" (http://geopandas.org/projections.html)). Es wurde in Kobe City herausgenommen, und das orthogonale Koordinatensystem der Ebene sollte je nach Gebiet etwas unterschiedlich sein.
| Name | proj4 string |
|---|---|
| Breiten- und Längengrad des geodätischen Systems Japans | +ellps=bessel +no_defs +proj=longlat |
| Breiten- und Längengrad des Weltvermessungssystems | +ellps=GRS80 +no_defs +proj=longlat |
| Orthogonales Koordinatensystem der Ebene des geodätischen Systems in Japan | +lon_0=134.3333333333333 +y_0=0 +lat_0=36 +ellps=bessel +no_defs +units=m +proj=tmerc +x_0=0 +k=0.9999 |
| Orthogonales Koordinatensystem der Ebene des geodätischen Weltsystems | +lon_0=134.3333333333333 +y_0=0 +lat_0=36 +ellps=GRS80 +no_defs +units=m +proj=tmerc +x_0=0 +k=0.9999 |
Das ebene orthogonale Koordinatensystem ist eine Zahl, die flach wie der Boden aussieht. Einheiten ist jetzt m. Bedeutet das, dass sich die elliptische Kugel zwischen dem japanischen geografischen System und dem geografischen Weltsystem im ebenen rechteckigen Koordinatensystem unterscheidet? Was als Japan-Messsystem bezeichnet wird, ist die Gefäßellipse, dh das ehemalige Japan-Messsystem. Da es heutzutage viele WGS84-Daten gibt, scheint es besser, GRS80 des "World Survey System" zu verwenden, das 2000 des Japan Survey System ist.
Wenn Sie verschiedene Berechnungen durchführen möchten, sollten Sie "Breite und Länge des Weltvermessungssystems" herunterladen, berechnen und bei der Erstellung einer endgültigen Zahl in ein ebenes rechteckiges Koordinatensystem konvertieren. Sie können es auf einmal mit to_crs () tun.
data.to_crs("+lon_0=134.3333333333333 +y_0=0 +lat_0=36 +ellps=GRS80 +no_defs +units=m +proj=tmerc +x_0=0 +k=0.9999").plot()
Muss ich beim Zuweisen von Werten über die Fläche im orthogonalen Koordinatensystem der Ebene berechnen?