MongoDB地理查询问题：带圆形挖空的正方形区域随机点查询报错

场景与需求 #

我的MongoDB数据库中lootPoint集合存储的Point类型地理数据格式如下：

{
    "_id" : "1856132932",
    "point" : {
        "type" : "Point",
        "coordinates" : [ 21.6550408, 50.1034841 ]
    }
}
{
    "_id" : "2184534902",
    "point" : {
        "type" : "Point",
        "coordinates" : [ 21.6560194, 50.1037961 ]
    }
}

需求为：查询位于指定正方形区域内，但排除中间20米圆形挖空范围的随机点。

原查询问题 #

最初编写的聚合查询运行时报错“无法解析$not字段”：

db.lootPoint.aggregate([
  {
    $match: {
      "point.coordinates": {
        $geoWithin: {
          $geometry: {
            type: "Polygon",
            coordinates: [
              [
                  [21.654939727403415,50.10333291665514],
                  [21.656736672596585,50.10333291665514],
                  [21.656736672596585,50.10453288334485],
                  [21.654939727403415,50.10453288334485],
                  [21.654939727403415,50.10333291665514]
              ]
            ]
          }
        },
        $not: {
          $geoWithin: {
            $centerSphere: [[21.6558382, 50.1039329], 20 / 3963.2]
          }
        }
      }
    }
  },
  { $sample: { size: 1 } }
])

改为$and组合形式后返回0条结果：

db.lootPoint.aggregate([
  {
    $match: {
      $and:[
        { "point.coordinates": {...} },
        { "point.coordinates": { $not:[ ... ] } },
      ]
    }
  },
  ...
])

错误原因分析 #

1. 第一个查询错误：$not不能直接与$geoWithin并列写在同一字段的查询条件中，MongoDB语法要求同一字段的多条件逻辑组合需通过$and或正确嵌套实现，不能将$not作为同层级条件。
2. 第二个查询错误：$not的用法错误，它后面应跟查询条件文档而非数组；同时原查询针对point.coordinates而非完整的GeoJSON point字段查询，可能导致地理空间逻辑识别异常。

正确查询语句 #

方案一：$and组合双条件 #

db.lootPoint.aggregate([
  {
    $match: {
      $and: [
        {
          "point": {
            $geoWithin: {
              $geometry: {
                type: "Polygon",
                coordinates: [
                  [
                    [21.654939727403415,50.10333291665514],
                    [21.656736672596585,50.10333291665514],
                    [21.656736672596585,50.10453288334485],
                    [21.654939727403415,50.10453288334485],
                    [21.654939727403415,50.10333291665514]
                  ]
                ]
              }
            }
          }
        },
        {
          "point": {
            $not: {
              $geoWithin: {
                $centerSphere: [[21.6558382, 50.1039329], 20 / 6378.1]
              }
            }
          }
        }
      ]
    }
  },
  { $sample: { size: 1 } }
])

关键修正：

• 针对完整的GeoJSON point字段查询，而非point.coordinates，确保MongoDB正确识别地理空间类型。
• $not后跟随正确的$geoWithin条件文档。
• 用地球平均半径公里数6378.1计算弧度，与米单位更匹配。

方案二：$nor反向排除 #

db.lootPoint.aggregate([
  {
    $match: {
      "point": {
        $geoWithin: {
          $geometry: {
            type: "Polygon",
            coordinates: [
              [
                [21.654939727403415,50.10333291665514],
                [21.656736672596585,50.10333291665514],
                [21.656736672596585,50.10453288334485],
                [21.654939727403415,50.10453288334485],
                [21.654939727403415,50.10333291665514]
              ]
            ]
          }
        },
        $nor: [
          {
            $geoWithin: {
              $centerSphere: [[21.6558382, 50.1039329], 20 / 6378.1]
            }
          }
        ]
      }
    }
  },
  { $sample: { size: 1 } }
])

替代优化方案 #

若数据量较大，先为point字段创建2dsphere索引提升查询效率：

db.lootPoint.createIndex({ "point": "2dsphere" })

索引创建后，上述查询会利用索引快速过滤文档，再执行$sample随机选取。

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