使用GeoJSON与Shapely判断点是否在多边形内：排查检测失败问题

嘿，我看了你的代码和问题，发现有两个关键问题导致检测一直返回False，下面给你拆解一下解决方案：

1. 坐标顺序搞反了！ #

这是地理空间处理里最容易踩的坑：Shapely的Point构造函数参数顺序是**(经度, 纬度)**，但你传入的(40.712776, -74.005974)是纬度在前、经度在后的顺序。而GeoJSON本身的坐标格式就是[经度, 纬度]，所以你的点其实不在美国的多边形范围内，自然返回False。

赶紧把点的创建代码改成这样：

point = Point(-74.005974, 40.712776)

2. 处理美国GeoJSON里的复杂几何类型 #

美国的GeoJSON包含阿拉斯加、夏威夷这些飞地，对应的几何类型是MultiPolygon（偶尔还会有GeometryCollection）。虽然Shapely的contains方法对MultiPolygon能直接生效，但为了避免遗漏子多边形，最好增加对子几何的遍历逻辑：

修正后的完整代码：

import json
from shapely.geometry import shape, Point

with open('USA.geo.json', 'r') as f:
    js = json.load(f)

# 修正坐标顺序：经度在前，纬度在后
point = Point(-74.005974, 40.712776)

for feature in js['features']:
    geom = shape(feature['geometry'])
    # 检查是否是包含多个子几何的类型（比如MultiPolygon）
    if hasattr(geom, 'geoms'):
        for sub_geom in geom.geoms:
            if sub_geom.contains(point):
                print('Found containing polygon:', feature)
                break
    else:
        if geom.contains(point):
            print('Found containing polygon:', feature)

额外验证技巧 #

如果还是不确定问题在哪，可以先打印出多边形的边界范围，看看你的点是不是在这个区间里：

print(polygon.bounds)  # 输出格式：(最小经度, 最小纬度, 最大经度, 最大纬度)

如果你的点的经纬度不在这个区间里，那肯定不会匹配成功，这时候优先检查坐标顺序哦。

内容的提问来源于stack exchange，提问作者Josh