如何基于GeoJSON城市选区文件，通过经纬度判断所在选区？（Node.js优先）

核心流程 #

• 解析GeoJSON文件：将费城选区的GeoJSON文件下载到本地后，解析为可操作的对象，提取其中的features数组——每个元素对应一个选区，包含边界几何信息（geometry字段）和选区属性（properties字段，如ID、名称）。
• 调整坐标格式：GeoJSON的坐标遵循**[经度, 纬度]**顺序，所以要把输入的市政厅坐标（纬度39.953306°、经度-75.163577°）转换为[-75.163577, 39.953306]的格式。
• 空间匹配判断：遍历每个选区的几何边界，判断目标坐标是否处于该多边形内部，这是核心环节。
• 返回匹配结果：找到包含目标坐标的选区后，输出其属性信息；无匹配则返回无结果。

Node.js实现方案 #

推荐使用turf.js地理空间库，它封装了成熟的空间运算方法，无需手动实现复杂算法。

1. 安装依赖 #

npm install @turf/turf

2. 代码实现 #

const fs = require('fs');
const turf = require('@turf/turf');

// 读取本地GeoJSON文件（假设文件名为philly-wards.geojson）
const geojsonData = JSON.parse(fs.readFileSync('./philly-wards.geojson', 'utf8'));

// 构造目标点（GeoJSON标准：[经度, 纬度]）
const targetPoint = turf.point([-75.163577, 39.953306]);

// 遍历选区匹配
let matchedWard = null;
for (const feature of geojsonData.features) {
  // 只处理多边形/多类型多边形选区
  if (['Polygon', 'MultiPolygon'].includes(feature.geometry.type)) {
    const isInside = turf.booleanPointInPolygon(targetPoint, feature);
    if (isInside) {
      matchedWard = feature.properties;
      break; // 找到匹配后终止遍历
    }
  }
}

// 输出结果
if (matchedWard) {
  console.log('目标坐标所属选区：', matchedWard);
} else {
  console.log('未找到匹配的选区');
}

优化提示 #

如果GeoJSON文件数据量较大，可使用turf.quadtree创建空间索引，减少不必要的遍历，提升查询效率。

学习资源推荐 #

基础概念 #

• GeoJSON规范：重点掌握Feature、Polygon、MultiPolygon的结构，牢记坐标顺序要求。
• 点-in-多边形算法：学习射线法的核心原理，了解边界点、孔洞等特殊场景的处理逻辑。

工具与库 #

• Turf.js官方文档：系统学习各类空间运算方法，覆盖点、线、多边形的常见操作。
• Node.js地理空间开发专题：关注技术博客中的相关实战内容，了解实际项目中的应用场景。

实战练习 #

• 用不同城市的GeoJSON文件测试代码，熟悉不规则多边形、带孔洞选区的处理。
• 学习空间索引概念，优化大数据量下的查询性能。

内容的提问来源于stack exchange，提问作者God's Drunkest Driver