有没有一种方法可以使用geopandas/shapely/fiona获得多边形的最高点？

Question

我有点形状文件和多边形形状文件。我想找出每个多边形内的最高点。我做了一个交集，以确定哪些点属于每个多边形：

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))

现在我想为每个索引选择最高点。我认为这是一个在几何列中按Z值排序的问题，但我在做这件事时遇到了问题。我不知道如何使用geopandas从几何体中提取Z坐标。最后，我想进行空间连接，并将Z值填充到最近的点（另一个shapefile）。

非常感谢。

Answer 1

我将重点讨论第一个问题，用它来提取z坐标

import geopandas as gpd
gdf = gpd.read_file("file.shp")
gdf['z'] = None
gdf['z'] = gdf.geometry.apply(lambda x: list(x.coords)[0][2])

或者，

z_values = [list(x.coords)[0][2] for x in gdf.geometry]
gdf['z'] = None
gdf['z'] = z_values

说明：对于每个几何体，从几何体中列出坐标，取坐标列表中的第一个[0]，取第三个[2]值，即几何体的z值（因为第一个[0]是x，第二个[1]是y值）。最近邻居的空间连接部分太长，无法写入。我建议就这个话题问一个新问题。