Django GEOS构建大范围多边形的查询异常问题

问题原因 #

核心问题出在多边形顶点的环绕方向上：

• GEOS遵循OGC标准，在EPSG:4326地理坐标系中，多边形的环绕方向决定了它代表的是实际覆盖区域还是补集：逆时针环绕对应正常区域，顺时针环绕会被解析为“挖空”的补集区域。
• 你创建的east多边形顶点顺序为顺时针，单独使用时GEOS会自动修正方向，但组合成MultiPolygon后，自动修正逻辑未触发，导致east被解析为“除东经0°到179°之外的区域”，因此东部的点无法匹配。
• 拆分为四个小多边形时，每个图形的范围足够小，环绕方向的偏差不会触发补集解析，所以查询正常。

解决方案 #

给每个多边形调用normalize()方法，自动调整顶点环绕方向为符合标准的逆时针：

from django.contrib.gis.geos import MultiPolygon, Polygon

# 调整west多边形的环绕方向
west = Polygon(((0, 85), (-179, 85), (-179, -85), (0, -85), (0, 85)), srid=4326).normalize()
# 调整east多边形的环绕方向
east = Polygon(((0, 85), (179, 85), (179, -85), (0, -85), (0, 85)), srid=4326).normalize()
# 组合为MultiPolygon
big = MultiPolygon(east, west)
feature.geography = big
feature.save()

验证 #

调整后，所有查询都应返回匹配结果：

from django.contrib.gis.geos import Point

queryset.filter(feature__geography__intersects=Point(-2.5, 2.5, srid=4326))
queryset.filter(feature__geography__intersects=Point(-2.5, -2.5, srid=4326))
queryset.filter(feature__geography__intersects=Point(2.5, -2.5, srid=4326))
queryset.filter(feature__geography__intersects=Point(2.5, 2.5, srid=4326))

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