如何利用R/ArcGIS将Excel经纬度匹配至新加坡子区Shapefile？

嘿，这个需求我经常碰到，给你两种靠谱的解决方案——一种适合批量自动化处理，另一种适合不想写代码的朋友：

这是效率最高的方式，能快速完成批量点的区域匹配：

• 首先安装必要的依赖库：

  pip install geopandas pandas openpyxl
  

• 然后运行以下代码（记得替换成你的文件路径和列名）：

  import geopandas as gpd
  import pandas as pd
  
  # 1. 加载新加坡子区域的Shapefile（确保.shp/.shx/.dbf等配套文件在同一文件夹）
  subregions = gpd.read_file("singapore_subregions.shp")
  # 可以先打印坐标系确认：print(subregions.crs)，新加坡常用EPSG:3414或EPSG:4326
  
  # 2. 加载带经纬度的Excel数据
  points_df = pd.read_excel("location_coords.xlsx")
  # 将经纬度转换为地理点对象，注意替换成你Excel里的经度/纬度列名
  points_gdf = gpd.GeoDataFrame(
      points_df,
      geometry=gpd.points_from_xy(points_df["经度"], points_df["纬度"]),
      crs="EPSG:4326"  # 大多数GPS/定位数据用WGS84坐标系
  )
  
  # 3. 统一坐标系：如果Shapefile和点数据的坐标系不一样，必须转换对齐
  if subregions.crs != points_gdf.crs:
      points_gdf = points_gdf.to_crs(subregions.crs)
  
  # 4. 执行空间连接：判断每个点属于哪个子区域
  # predicate="within" 表示找包含该点的子区域
  result_gdf = gpd.sjoin(points_gdf, subregions, how="left", predicate="within")
  
  # 5. 导出结果到Excel（去掉几何列，保留原始数据+匹配的子区域信息）
  result_gdf.drop(columns="geometry").to_excel("matched_locations.xlsx", index=False)
  

• 小提示：

  • 如果你的Excel里经纬度列名不是"经度"/"纬度"，一定要对应修改代码里的字段名
  • 如果有些点不在任何子区域内（比如海上），结果里的子区域字段会显示NaN，可以后续手动处理
  • 如果子区域有重叠，how="left"会返回第一个匹配的区域，需要的话可以调整参数获取所有重叠区域

如果你不太熟悉Python，用QGIS这个开源GIS工具也能轻松完成：

• 步骤1：打开QGIS，导入新加坡子区域Shapefile：
  点击菜单栏「图层」→「添加矢量图层」，选择你的.shp文件导入

• 步骤2：导入Excel的经纬度数据：
  点击菜单栏「图层」→「添加图层」→「添加分隔文本图层」，选择你的Excel文件；在弹出的窗口里，指定「X字段」为经度列，「Y字段」为纬度列，坐标系选择EPSG:4326，点击确定导入

• 步骤3：执行空间匹配：
  右键点击导入的点图层→「属性」→「连接/关联」→「添加连接」；在连接设置里：

  • 连接图层选择你的子区域Shapefile图层
  • 连接类型选择「取第一个匹配的要素」
  • 匹配条件选择「点在面内」
    设置完成后点击确定，点图层就会带上对应的子区域信息

• 步骤4：导出结果：
  右键点击点图层→「导出」→「保存要素为」，选择保存格式为「Excel」，设置保存路径即可

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