python - geopandas 成对空间操作

Question

我需要在两个地理数据框之间进行成对包含操作。geo_listings包含点几何并geo_listings_buffer包含与在第一个数据帧上运行的结果相对应的多边形几何buffer。

该contains函数在执行元素操作后返回一个系列。我希望返回一个方阵。

到目前为止，我有这段代码可以工作，但速度很慢：

pair_contains = np.zeros((geo_listings.shape[0],geo_listings.shape[0]))
for i, row in geo_listings.iterrows():
    pair_contains[i,:] = geo_listings_buffer.contains(row['geometry'])

我有大约 8000 行，这段代码需要 8-9 分钟才能运行。

任何改进这将是有帮助的。

标签： pythonpandasnumpygeopandas