一連のポイント (約 240 億ポイント、400 万行) に最も近いラインを見つけようとしています。ポイントは 1 つの GeoDataFrame に存在し、ラインは別の 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 ラインとポイントに適用していることに注意してください)。