NodeJS下GeoJSON点位所属多边形的高效查询实现方案

直接全量遍历多边形做点在多边形（PIP, Point in Polygon）判断的方式，在行政区要素数量多、查询并发高的时候确实性能很差，完全不需要自己从零造轮子，现有成熟的开源能力完全可以满足你匹配GPS点到对应行政区多边形的需求，而且适配你“无重叠多边形”的场景逻辑会更简单。

常用实现依赖选型 #

• turf.js 地理计算生态
  这是目前最适配标准GeoJSON格式（也就是你提到的FeatureCollection结构）的地理计算工具集，本身就内置了点在多边形判断的基础能力，还支持配合空间索引做查询加速：
  • 基础的点判断能力来自@turf/boolean-point-in-polygon方法，单独使用就是纯遍历逻辑，但如果搭配同生态的R树空间索引模块@turf/rbush，可以预先给每个多边形算好最小外接包围盒（MBR），查询时先通过包围盒快速排除90%以上完全不可能覆盖当前点的多边形，只对极少数候选对象做精确的多边形边相交校验，性能比纯遍历高一个数量级以上。
• 轻量独立空间索引实现
  如果不想引入全量turf的依赖包，可以单独使用通用空间索引库rbush，自己为每个多边形计算最小外接矩形存入索引，查询时先按矩形范围过滤候选多边形，再搭配极简的射线法实现点在多边形判断逻辑，包体积更小，适合对依赖大小敏感的场景。
• 网格索引方案
  如果数据量达到几十万级别的行政区规模、或者需要支撑高并发GPS点查询，可以自己实现固定经纬度网格的索引：按固定步长把全球/辖区范围切成网格，每个网格绑定覆盖该范围的多边形ID，查询时直接通过经纬度计算当前点所属的网格，直接获取对应多边形做匹配，查询时间复杂度可以到O(1)，比R树索引的查询响应更快。

参考封装实现逻辑 #

因为你的场景不存在重叠多边形，逻辑不需要处理多结果冲突，只要返回精确命中的多边形即可。

import rbush from '@turf/rbush';
import { booleanPointInPolygon } from '@turf/boolean-point-in-polygon';
import { bbox } from '@turf/bbox';

/**
 * 初始化多边形查询能力（仅在服务启动时执行一次，不要每次查询都调用）
 * @param {FeatureCollection} featureCollection 行政区GeoJSON要素集合
 * @returns {Function} 可复用的坐标查询方法
 */
function initPolygonQuery(featureCollection) {
  // 预构建R树空间索引，所有计算只在初始化阶段跑一次
  const spatialIndex = rbush();
  const indexRecords = featureCollection.features.map(feature => {
    const [minLng, minLat, maxLng, maxLat] = bbox(feature);
    return {
      minX: minLng,
      minY: minLat,
      maxX: maxLng,
      maxY: maxLat,
      sourceFeature: feature
    };
  });
  spatialIndex.load(indexRecords);

  // 暴露对外的查询方法，和你预期的`getPolygonFromPoint(FeatureCollection, x, y)`调用逻辑对齐
  return function getPolygonFromPoint(x, y) {
    // 第一步：通过空间索引快速筛选可能覆盖当前点的候选多边形
    const candidates = spatialIndex.search({
      minX: x,
      minY: y,
      maxX: x,
      maxY: y
    });
    // 第二步：对候选多边形做精确的点归属判断
    const queryPoint = { type: 'Point', coordinates: [x, y] };
    for (const item of candidates) {
      if (booleanPointInPolygon(queryPoint, item.sourceFeature)) {
        return item.sourceFeature;
      }
    }
    // 没有匹配到则返回空
    return null;
  }
}

// 实际使用示例
// const yourAdminDivisionData = require('./your_admin_division.geojson');
// const getMatchedPolygon = initPolygonQuery(yourAdminDivisionData);
// const matched = getMatchedPolygon(116.404, 39.915);

实践注意事项 #

• 空间索引只需要在服务启动、加载完GeoJSON数据时构建一次，不要每次坐标查询都重新计算索引，这是最容易踩的性能坑。
• 注意坐标系统一：GPS原始输出是WGS84坐标系，确保你的行政区GeoJSON数据也使用同一坐标系，不要混用GCJ02之类的偏移坐标系，否则会出现匹配错位的问题。
• 你当前场景不存在重叠多边形，不需要设计复杂的多边形优先级判断逻辑，拿到精确命中的第一个结果直接返回即可。

内容的提问来源于stack exchange，提问作者João Pimentel Ferreira