空間データセットのオーバーレイ

Pythonで扱う地理空間データ

Joris Van den Bossche

Open source software developer and teacher, GeoPandas maintainer

ポリゴンとの交差

countries.intersection(circle)

ポリゴンとの交差

countries.intersection(circle) の制約:

GeoSeries と単一ポリゴンのみを交差
属性情報が保持されない

2つのデータセットをオーバーレイ

countries.plot()

geologic_regions.plot()

2つのデータセットをオーバーレイ

geopandas.overlay(countries, geologic_regions, 
                  how='intersection')

オーバーレイ vs 交差

交差 (単一ポリゴン)

countries.intersection(geologic_region_A)

0                            ()
1    POLYGON ((-1.661 48.803...
2    POLYGON ((1.201 51.145,...
dtype: object

オーバーレイ手法

geopandas.overlay(countries, geologic_regions,
                  how='intersection')

     name geologic_region              geometry
1  France               C  POLYGON ((2.5 51....
2      UK               C  POLYGON ((0.7 52 ...
3  France               B  POLYGON ((-1.7 46...
..    ...             ...                   ...

練習してみましょう！

Pythonで扱う地理空間データ