如何查询MongoDB中先后经过两个指定GeoJSON区域的空间轨迹

原有方案问题说明 #

你当前使用的连续$match仅能筛选出同时经过R1、R2两个区域的轨迹，无法判断两个区域的先后经过顺序；如果未给地理字段建立空间索引，两次$match都会触发全表扫描，数据量大时查询效率极低。

前置优化（必须做，否则性能无法达标） #

给存储经纬度的字段建立2dsphere空间索引，大幅降低地理匹配的计算开销：

db.你的集合名.createIndex({"location.lonlat": "2dsphere"})

支持方向判断的高效聚合管道 #

核心逻辑是：先用地理索引快速过滤出同时经过两个区域的轨迹（干掉99%以上不相关数据），再仅对少量符合基础条件的轨迹做顺序校验，避免全量遍历所有轨迹的坐标数组。
聚合管道完整配置如下：

[
  // 第一阶段：用地理索引快速过滤同时经过R1、R2的轨迹，减少后续计算量
  {
    "$match": {
      "$and": [
        {"location.lonlat": {"$geoIntersects": {"$geometry": gjs[0]}}},
        {"location.lonlat": {"$geoIntersects": {"$geometry": gjs[1]}}}
      ]
    }
  },
  // 第二阶段：计算轨迹首次经过R1、R2的位置索引（默认你的轨迹数组按时间正序存储）
  {
    "$addFields": {
      "firstR1Index": {
        "$reduce": {
          "input": "$location.lonlat",
          "initialValue": -1,
          "in": {
            "$cond": [
              {"$and": [{"$eq": ["$$value", -1]}, {"$geoWithin": {"$geometry": gjs[0], "point": "$$this"}}]},
              "$$index",
              "$$value"
            ]
          }
        }
      },
      "firstR2Index": {
        "$reduce": {
          "input": "$location.lonlat",
          "initialValue": -1,
          "in": {
            "$cond": [
              {"$and": [{"$eq": ["$$value", -1]}, {"$geoWithin": {"$geometry": gjs[1], "point": "$$this"}}]},
              "$$index",
              "$$value"
            ]
          }
        }
      }
    }
  },
  // 第三阶段：过滤出先经过R1、后经过R2的轨迹
  {
    "$match": {
      "$expr": {"$lt": ["$firstR1Index", "$firstR2Index"]}
    }
  },
  // 可选阶段：移除临时计算的辅助字段
  {
    "$project": {"firstR1Index": 0, "firstR2Index": 0}
  }
]

额外性能优化建议 #

• 如果单条轨迹的坐标点数量超过1000，可在写入轨迹时提前给每个坐标点加上采集时间戳/顺序序号，校验顺序时直接对比首次命中点的时间戳即可，无需计算数组索引
• 如果查询频率极高，可将常见的区域组合的查询结果做缓存，进一步降低数据库压力

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