Geopandas
Pythonの新しいシリーズ本を依頼されたので勉強中なのだが,geopandasというのが面白い.いわゆるGISのためのPythonパッケージなのだが,データ解析の定番のpandasのようにデータフレームで地理的なデータを解析できる.
まずモジュールを読み込む.
import geopandas as gpd
国土交通省データベース
http://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-P20.html
の避難所データをもってくる.これはshapeファイルといってGISの基本データだが,geopandasだと直接読める.例として熊本の避難所データを読んでみる.
df = gpd.read_file("kumamoto/P20-12_43.shp")
データフレームと同じように条件でフィルタできる.たとえば,収容できる人数(列名は単なる番号で005が避難数)が500人以上の避難所をプロットすると,上の図が出てくる.
df[ df["P20_005"]>=500 ].plot()
全部で140箇所ある.内閣府のプッシュ方式だと500人以上の避難所に運んだと言われている.ただし,データを調べてみると,収容人数データはほとんどが-1でいわゆる欠損値だ.きちんとデータを整備するとともに,報告書で嘘を書かないようにすることが重要だ.
この記事が気に入ったらサポートをしてみませんか?