ST_Intersects对不相交Geography返回True的原因排查

地理数据示例 #

select 'SRID=4326;MULTIPOLYGON (((105.67151103718783 26.882512088629863, 115.99865947468783 31.74650295647645, 116.37219463093783 43.48451247459262, 147.02954326375033 62.362161196046245, -177.2690102030465 66.81825017074695, -168.77863849894493 66.26641372873476, -169.78045429484337 60.98223428128946, 165.38220439656283 44.47269646435982, 172.99572978718783 23.382219202927942, 168.86487041218783 -11.63112025581975, 140.82776103718783 -11.372743681342062, 130.19299541218783 -9.644484126039586, 124.04065166218783 -12.662180464457489, 103.91369853718783 -12.23305845779793, 90.37854228718783 6.9914492544889555, 88.79651103718783 23.945718564044213, 105.67151103718783 26.882512088629863)))'::geography
union all
select 'SRID=4326;POINT (0 80)'::geography

数据可视化情况 #

该数据在DBeaver中显示时，由于反子午线（180°经线）的存在，多边形出现显示变形，但它实际代表的是东南亚区域。从视觉上看，点(0 80)与多边形在几何和地理上均不相交。

异常查询结果 #

执行以下ST_Intersects查询时，却返回True：

select st_intersects(
'SRID=4326;MULTIPOLYGON (((105.67151103718783 26.882512088629863, 115.99865947468783 31.74650295647645, 116.37219463093783 43.48451247459262, 147.02954326375033 62.362161196046245, -177.2690102030465 66.81825017074695, -168.77863849894493 66.26641372873476, -169.78045429484337 60.98223428128946, 165.38220439656283 44.47269646435982, 172.99572978718783 23.382219202927942, 168.86487041218783 -11.63112025581975, 140.82776103718783 -11.372743681342062, 130.19299541218783 -9.644484126039586, 124.04065166218783 -12.662180464457489, 103.91369853718783 -12.23305845779793, 90.37854228718783 6.9914492544889555, 88.79651103718783 23.945718564044213, 105.67151103718783 26.882512088629863)))'::geography,
'SRID=4326;POINT (0 80)'::geography
)

原因分析 #

这是因为PostGIS处理geography类型多边形时，默认遵循"多边形内部是地球表面面积较小的区域"规则。你的多边形跨越了反子午线，当前坐标序列实际定义的是包含北极点附近的"大区域多边形"，而你预期的东南亚区域反而成了这个多边形的外部。

当用ST_Intersects判断点(0 80)与该多边形的关系时，PostGIS认为这个点位于多边形内部，因此返回True。

要解决这个问题，需要调整多边形的顶点顺序：

• 对于geography类型，多边形顶点的顺时针/逆时针顺序决定了内部区域范围，尤其是跨越180°经线或接近极点时。
• 可以使用ST_Reverse函数反转顶点顺序，或者重新整理坐标序列，确保定义的是你需要的较小区域（东南亚）。

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