MongoDB无孔GEOJSON多边形$geoWithin查询触发环包含校验错误

问题触发场景 #

使用PyMongo执行带$geoWithin条件的地理空间检索时，传入预期无孔洞的GeoJSON多边形作为查询范围，在遍历返回的查询游标时抛出运行时错误，错误原文为：
OperationFailure: Secondary loops not contained by first exterior loop - secondary loops must be holes
对应的复现逻辑代码如下：

import pymongo
import json

def loadjson(filename):
    with open(filename) as f:
        mydict = json.load(f)
    return mydict

client = pymongo.MongoClient('mongodb://' \
                                    + 'username' \
                                    + ':' \
                                    + "password" + \
                                    '@localhost:num')

gjOr=loadjson('polygon.json')

database = client.database
collection = database["name"]
filParam = {"field": {"$geoWithin": {"$geometry": gjOr["geometry"]}}}
exParam = {"parameter": 1} 

cursorItineraries = collection.find(filParam, exParam).batch_size(2)
for i in cursorItineraries:
    print(i)

报错本质原因 #

MongoDB内置的GeoJSON几何校验规则为：

• 单个Polygon类型的几何对象中，coordinates数组的第一个子环，被强制识别为该面的唯一外部边界
• 排在第一个环之后的所有子环，都会被校验为「嵌入外部环内部的孔洞」，MongoDB不支持单个Polygon内存在多个互不相连的独立外部环

用户感知到的“多边形无孔洞”和MongoDB报的“次级环必须是孔洞”不冲突，本质是你传入的几何对象结构不符合规范，才会触发这类校验失败，常见触发场景：

• 多块互不相连的独立查询范围，被错误封装为单个Polygon类型，每一块独立面的边界环被当成了Polygon的次级环。这类场景本应使用MultiPolygon类型承载多块独立面，而非单个Polygon
• 坐标数组存在层级嵌套错误：比如单个无孔洞多边形的coordinates多套了一层数组，导致本来应该作为首个环的边界坐标被挤到次级环位置，MongoDB读到的第一个环是无效的空/零散坐标，自然判定后续环不被其包含
• 环的闭合逻辑异常：本该属于同一个环的坐标序列，因为首尾点不匹配被几何引擎拆成两个独立环，进而触发次级环校验失败

修复方案 #

• 几何类型匹配修正：如果查询范围是多块不相连的独立面，将GeoJSON的type字段从Polygon改为MultiPolygon，把每一块独立面的坐标结构单独嵌套在MultiPolygon的coordinates子层级中
• 结构层级校验：如果确实是单个无孔洞的多边形，检查coordinates的嵌套层级，正确的无孔洞多边形结构应为coordinates: [[[lon1, lat1], [lon2, lat2], ..., [lon1, lat1]]]，即数组内仅保留1个闭合的外环坐标序列，不要额外嵌套冗余数组
• 环有效性校验：确认所有环的首尾坐标完全一致，满足GeoJSON对闭合环的要求，无孔洞的Polygon内不要塞入其他独立的坐标环

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