SQL Server STWithin/STIntersects判断不准确及墨卡托投影适配问题

问题根源 #

SQL Server的geography类型严格遵循OGC标准，多边形的顶点环绕方向直接决定了它代表的区域：

• 针对WGS84坐标系（SRID=4326），顺时针顶点顺序的多边形代表的是直观可见的内部区域；如果是逆时针顺序，SQL会把它解析成「整个地球减去该多边形」的剩余区域。
• 你创建的大范围多边形顶点顺序是西南→东南→东北→西北→西南，在球面坐标系下属于逆时针方向，导致实际匹配的是你预期范围之外的区域，所以点不在查询结果里。

解决方案 #

1. 一键反转多边形方向 #

用ReorientObject()方法直接修正环的方向，让SQL正确识别内部区域：

declare @swLng float, @swLat float, @neLng float, @neLat float
declare @g geography, @poly varchar(max)

select @swLng=-147.45514585847684, @swLat=28.841195680455485, @neLng=-59.86599816500535, @neLat=56.14838208179307
set @poly = 'POLYGON((' + cast(@swLng as varchar) + ' ' +  cast(@swLat as varchar) + ',' + cast(@neLng as varchar) + ' ' + cast(@swLat as varchar) + ',' + cast(@neLng as varchar) + ' ' + cast(@neLat as varchar) + ',' + cast(@swLng as varchar) + ' ' + cast(@neLat as varchar) + ',' +  cast(@swLng as varchar) + ' ' +  cast(@swLat as varchar) + '))'
set @g = geography::STGeomFromText(@poly, 4326);
-- 反转环方向，修正区域判定
set @g = @g.ReorientObject();

-- 现在查询就能返回正确结果
select * from #gtest where GeoPoint.STWithin(@g) = 1
select * from #gtest where GeoPoint.STIntersects(@g) = 1

2. 手动调整顶点顺序 #

直接把顶点改成顺时针顺序，比如调整为西南→西北→东北→东南→西南：

set @poly = 'POLYGON((' + cast(@swLng as varchar) + ' ' +  cast(@swLat as varchar) + ',' + cast(@swLng as varchar) + ' ' + cast(@neLat as varchar) + ',' + cast(@neLng as varchar) + ' ' + cast(@neLat as varchar) + ',' + cast(@neLng as varchar) + ' ' + cast(@swLat as varchar) + ',' +  cast(@swLng as varchar) + ' ' +  cast(@swLat as varchar) + '))'
set @g = geography::STGeomFromText(@poly, 4326);

3. Mapbox墨卡托投影适配 #

Mapbox用的是Web墨卡托投影（SRID=3857），如果需要在SQL中处理墨卡托坐标：

• 将WGS84的地理对象转换为墨卡托投影的几何对象（geometry类型）后再判断：

-- 把大范围多边形转成墨卡托投影
declare @g_merc geometry = @g.ReorientObject().STTransform(3857);
-- 点也同步转换后执行判断
select * from #gtest where geography::Point(Lat, Lng, 4326).STTransform(3857).STWithin(@g_merc) = 1

注意：墨卡托投影在高纬度或跨半球大范围场景下会有明显变形，这类场景优先用geography类型并保证环方向正确。

验证方法 #

用STArea()可以快速验证多边形区域是否正确：

-- 未修正的多边形面积会接近地球总表面积（约5.1亿平方公里）
select @g.STArea()
-- 修正后的多边形面积应该是你预期的区域大小
select @g.ReorientObject().STArea()

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