我试图找到离一堆点最近的线(大约 240 亿点,400 万条线)。点存在于一个 GeoDataFrame 中,而线存在于另一个中。我试图遵循这个:https://github.com/geopandas/geopandas/issues/140,并做到了:
lines_sidx = lines_df['geom'].sindex
[list(lines_sidx.intersection((points.loc[i,'geom'].y, points.loc[i,'geom'].x))) for i in range(len(points))]
这只是返回一个空的列表列表。这里发生了什么?
(请注意,我将其应用于两个数据集中的前 100 行和点)。