0

経度と緯度があり、予想される結果は、ポイントがどのマルチポリゴンにあるかに関係なく、マルチポリゴンの名前または ID を取得することです。

import geopandas as gpd

world = gpd.read_file('/Control_Areas.shp')
world.plot()

出力

   0    MULTIPOLYGON (((-9837042.000 6137048.000, -983...
   1    MULTIPOLYGON (((-11583146.000 5695095.000, -11...
   2    MULTIPOLYGON (((-8542840.287 4154568.013, -854...
   3    MULTIPOLYGON (((-10822667.912 2996855.452, -10...
   4    MULTIPOLYGON (((-13050304.061 3865631.027, -13.

以前の試み: フィオナ、シェイプリー、ジオパンダを試してみましたが、これを進めるのにひどく苦労しました。私が得た最も近いものは within と contains 関数ですが、私が苦労した作業の領域は、マルチポリゴンからポリゴンへの変換も成功させ、inside と contains の力を利用して目的の出力を取得することです。

シェープファイルはこちらからダウンロードしました。

4

1 に答える 1