python - 如何使用 Python 中的 geopandas 测试 Point 是否在 Polygon/Multipolygon 中？

Question

我从网站上有来自美国的多边形数据
arcgis
我还有一个包含城市坐标的excel文件。我已将坐标转换为几何数据(点)。
现在我想测试点是否在美国。
两者都是 dtype:几何。我想用这个我可以很容易地进行比较，但是当我使用我的代码时，我得到的每个点的答案都是错误的。即使有在美国的积分。

代码是:

import geopandas as gp
import pandas as pd
import xlsxwriter
import xlrd
from shapely.geometry import Point, Polygon

df1 = pd.read_excel('PATH')
gdf = gp.GeoDataFrame(df1, geometry= gp.points_from_xy(df1.longitude, df1.latitude))

US = gp.read_file('PATH')

print(gdf['geometry'].contains(US['geometry']))

有谁知道我做错了什么？

Answer 1

contains在 GeoPandas 中，目前以一对一的方式工作，而不是一对多。为此，请使用 sjoin .

points_within = gp.sjoin(gdf, US, op='within')

这将仅返回 US 内的那些点.或者，您可以过滤包含点的多边形。

polygons_contains = gp.sjoin(US, gdf, op='contains')