非常に単純なジオパンダの例を作成しました(以下を参照)。それは機能しますが、私が世界のカスタム部分を持つことができることが重要であることに気付きました。時にはドイツ、時にはベルリンだけ。 (また、geopandasファイルでポリゴンとして定義する領域ごとに、持っているデータを集約したいのですが、別の質問でこれを追加します。)
どうすれば異なる「ベースマップ」を取得できますか
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
視覚化のために?
# 3rd party modules
import pandas as pd
import geopandas as gpd
import shapely
# needs 'descartes'
import matplotlib.pyplot as plt
df = pd.DataFrame({'city': ['Berlin', 'Paris', 'Munich'],
'latitude': [52.518611111111, 48.856666666667, 48.137222222222],
'longitude': [13.408333333333, 2.3516666666667, 11.575555555556]})
gdf = gpd.GeoDataFrame(df.drop(['latitude', 'longitude'], axis=1),
crs={'init': 'epsg:4326'},
geometry=[shapely.geometry.Point(xy)
for xy in Zip(df.longitude, df.latitude)])
print(gdf)
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
base = world.plot(color='white', edgecolor='black')
gdf.plot(ax=base, marker='o', color='red', markersize=5)
plt.show()
geopandas.datasets.get_path(...)
documentation に書かれているように、実行する必要があります
_>>> geopandas.datasets.available
['naturalearth_lowres', 'naturalearth_cities', 'nybb']
_
どこ
「ドイツのシェイプファイル」を検索すると arcgis.comのURL が得られ、「BundesamtfürKartographie undGeodäsie」をソースとして使用しました。 _vg2500_geo84/vg2500_krs.shp
_を使用した結果は次のようになります。
ソース:
©BundesamtfürKartographie undGeodäsie、Frankfurt am Main、2011Vervielfältigung、Verbreitung undöffentlicheZugänglichmachung、auch auszugsweise、mit Quellenangabe gestattet。
base.set_aspect(1.4)
も設定しなければなりませんでした。値_1.4
_は試行錯誤によって見つかりました。
ベルリンのこのようなデータの別のソースは daten.berlin.de です
ジオパンダがシェープファイルを読み取るとき、それは列を持つジオパンダのデータフレームです
_['USE', 'RS', 'RS_ALT', 'GEN', 'SHAPE_LENG', 'SHAPE_AREA', 'geometry']
_
で:
USE=4
_すべての要素RS
は、16077または01003のような文字列ですRS_ALT
_は、160770000000または010030000000などの文字列ですGEN
は_'Saale-Holzland-Kreis'
_または_'Erlangen'
_のような文字列ですSHAPE_LENG
_は_202986.1998816
_や_248309.91235015
_のようなフロートですSHAPE_AREA
_は_1.91013141e+08
_や_1.47727769e+09
_のようなフロートですgeometry
は形状の良いジオメトリです-ほとんどがポリゴン