ポイントシェープファイルとポリゴンシェープファイルがあります。各ポリゴン内の最高点を見つけたいと思います。各ポリゴンに属するポイントを見つけるために交差を行いました。
import geopandas as gpd
from geopandas.tools import sjoin
point = gpd.GeoDataFrame.from_file(pointSHP)
print("POINT", point)
poly = gpd.GeoDataFrame.from_file(polygonSHP)
print("POLY", poly)
points_within_poly = gpd.sjoin(point, poly, how="inner", op='intersects')
print(points_within_poly.head(10))
ここで、各 index_right の最高点を選択したいと思います。ジオメトリ列の Z 値によるソートの問題だと思いますが、それを行うのに問題があります。geopandas を使用してジオメトリから Z 座標を抽出する方法がわかりません。最後に、空間結合を行い、Z 値を最も近い点 (別のシェープファイル) に設定したいと思います。
ありがとうございました