PostGIS ST_Within与MongoDB $geoWithin功能及结果一致性问询

我正在对比PostGIS与MongoDB的空间查询功能，以挑选适配工作需求的产品。我在两个数据库中导入了完全相同的模拟数据，随后用地理围栏查询统计指定多边形范围内的元素数量，但返回结果存在明显差异：MongoDB未查询到范围内的Polygon类型数据，且POINT类型数据的统计数量比PostGIS多。想请教这是查询代码的问题，还是两款数据库的空间函数存在底层差异？

PostGIS 查询语句及结果 #

首先你的PostGIS语句存在笔误：FERE应为FROM，修正后的完整语句如下：

-- # For Postgis
SELECT 
    type
    , count(id)
FROM your_table_name
WHERE ST_Within(geojson, ST_GeomFromText('POLYGON((100.95725876372013 16.84256184570458,101.08808706513712 18.77654822388623,102.55103520605854 16.583603557822975,100.95725876372013 16.84256184570458))'))
GROUP BY type

查询结果：

type        count
"POINT"     38990
"POLYGON"       1

MongoDB 查询语句及结果 #

# for Mongodb
db.geometryRoot.aggregate([
    {
        $match: {
            geojson: {
                $geoWithin: {
                    $geometry: {
                        type: "Polygon",
                        coordinates: [[
                            [100.95725876372013, 16.84256184570458], 
                            [101.08808706513712, 18.77654822388623], 
                            [102.55103520605854, 16.583603557822975], 
                            [100.95725876372013, 16.84256184570458]
                        ]]
                    }
                }
            }
        }
    },
    {
        $group: {
            _id: "$type",
            count: { $count: { } }
        }
    }
])

查询结果：

{
  _id: 'POINT',
  count: 39175
}

差异原因拆解 #

1. 代码与配置问题 #

• PostGIS笔误：原语句的FERE是明显的拼写错误，必须修正为FROM才能正常执行，推测你实际运行时已修正，否则无法得到结果。
• MongoDB索引缺失：如果geojson字段未创建2dsphere索引，MongoDB的$geoWithin无法正确处理非点几何（比如Polygon），直接导致这类数据被过滤。这是MongoDB空间查询的硬性要求，必须先创建索引：

  db.geometryRoot.createIndex({ geojson: "2dsphere" })
  

2. 底层空间逻辑差异 #

• Polygon匹配规则不同：
  • PostGIS的ST_Within(A,B)遵循OGC标准，判断A的所有点是否都在B的内部或边界上，只要满足就会返回。
  • MongoDB的$geoWithin对于非点几何，默认要求整个几何对象完全被包含在查询多边形内，且查询多边形必须符合GeoJSON的逆时针方向要求（顺时针会被视为内部区域反转）。如果你的目标Polygon只是部分重叠，或者查询多边形方向不对，MongoDB都不会返回。
• 点数据计数差异：
  • 两者对边界点的判定逻辑存在精度差异：PostGIS对浮点坐标的精度处理更严格，而MongoDB的$geoWithin可能将刚好落在多边形边界上的点计入结果，或者反之。
  • 坐标顺序验证：MongoDB严格遵循GeoJSON标准（[经度, 纬度]），而PostGIS的WKT格式是(经度 纬度)，如果你的数据存储时坐标顺序不一致（比如PostGIS存的是[纬度, 经度]），会导致点的匹配范围偏差。

验证步骤 #

1. 给MongoDB的geojson字段创建2dsphere索引后重新执行查询，看是否能返回Polygon数据。
2. 反转MongoDB查询多边形的坐标顺序（改成逆时针），验证是否影响Polygon的匹配结果。
3. 抽取几个PostGIS未计数、但MongoDB计数的点，分别用两个数据库的空间函数单独查询，定位边界点的判定差异。
4. 确认两个数据库中存储的GeoJSON坐标顺序完全一致（均为[经度, 纬度]）。

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