Pythonで地物のポリゴン(境界データ)を統合する
/
Pythonで地物のポリゴン(境界データ)を統合します.
# ファイル構成
project_dir
├── /data
│ ├── /R2_boundary/04_miyagi/04_miyagi.shp
│ └── (省略)
└── merge_bound.ipynb <- 実行用ノートブック
# ポリゴン(境界データ)の統合
ポリゴン(境界データ)を統合します.
import pandas as pd
import geopandas as gpd
import matplotlib.pyplot as plt
from shapely.ops import unary_union
# 境界データのshpファイル読込
pref_gdf = gpd.read_file("./data/R2_boundary/04_miyagi/04_miyagi.shp")
fig = plt.figure()
ax = fig.add_subplot()
# 境界データをmerge
bound = gpd.GeoSeries(unary_union(pref_gdf.geometry))
# 描画
bound.plot(ax=ax, color='none', edgecolor='b', linewidth=1)
plt.show()
# 出力
bound["pref"] = "miyagi"
# まとめ
Pythonで地物のポリゴン(境界データ)を統合しました.