Polygons.contains()不工作(多边形有面积且is_valid = True)

Question

我正在尝试从给定的经度和纬度的.shp文件中获取状态：

我已经加载了.shp (地图)和.csv (纬度经度值)的文件

我设法返回了一张里面有所有点的地图：

    geometry = [Point (xy) for xy in zip(Atrain["lon"],Atrain["lat"])]
    geo_df = gpd.GeoDataFrame(geometry = geometry)

    fig, ax = plt.subplots(figsize = ((15,15)))
    shp_br.plot(ax =ax)
    geo_df.plot(ax = ax,color = "red")

结果是：

https://i.stack.imgur.com/LFGwG.png

(对不起，没有足够的声誉来发布图片，网站的新用户)

但是，我的循环导致了一个应该存在的值的false：

    pt = Point(-10, -50)
    s = 0
    for i in shp_br["geometry"]:
        if i.contains(pt):
            print(shp_br["NM_ESTADO"][s])
        s +=1

有可能是什么原因造成的线索吗？谢谢!

Answer 1

我想明白了，我的观点不应该出现在那里，因为我把它写成了一种格式(Lat，Lon)，而它应该是(Lon，Lat)