为何CosmosDB的ST_WITHIN、ST_INTERSECTS查询与经纬度范围查询结果不一致

导致查询结果不一致的核心原因 #

1. 坐标计算模型本质差异
   显式经纬度范围查询是将经纬度作为平面笛卡尔坐标系下的数值直接做大小比对，完全忽略地球曲率影响。而CosmosDB的geography类型默认基于WGS84椭球模型做空间计算，你构造的矩形多边形在球面上的实际边界不是平面坐标系下平行于经纬线的直线，部分在平面数值上刚好落在边界范围内的点，在椭球模型计算中会被判定在多边形外部。
2. 空间索引精度的近似计算误差
   CosmosDB的地理空间索引采用网格分层的近似索引策略，默认精度等级为11。精度等级越低，网格划分颗粒度越粗，边界附近的微小坐标偏差会被索引的近似判定逻辑误判为不在范围内，你遇到的极小数值偏差导致点位被排除就是该机制导致的。
3. 浮点精度与边界判定规则差异
   显式范围查询使用>=/<=的硬规则匹配数值，只要数值满足条件就会返回。而空间函数计算过程中会存在浮点运算舍入误差，哪怕点位经纬度理论上与边界完全重合，在内部转换计算时也可能因为舍入偏差被判定在边界外。另外你存储的geopoint字段如果精度高于你构造Polygon时传入的坐标精度，也会出现匹配偏差。
4. 额外排查项：检查geopoint的存储格式是否符合GeoJSON规范，要求坐标顺序为[经度, 纬度]，如果存储时顺序写反也会出现边界匹配异常。

修复方案 #

• 调整容器的空间索引精度到最高等级12，降低近似判定的误差
• 构造查询用的Polygon时，将边界向外缓冲1e-7度左右（约1厘米级的误差，不会影响业务范围），抵消浮点运算的偏差
• 对边界附近的点位可增加显式经纬度校验做二次兜底匹配
• 如果你的业务场景不需要考虑地球曲率，可将空间类型从geography改为geometry平面类型，计算逻辑就和显式经纬度比对完全一致

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