理解Shapefile的非经纬度坐标点及空间连接失效原因

嘿，这个情况太常见啦！你看到的那些超大数字坐标，说明你的邮编Shapefile用的是投影坐标系（不是咱们常用的经纬度地理坐标系）——大概率是UTM分区或者墨西哥本地的专用投影，比如墨西哥城常用的EPSG:32614（UTM 14N）或者EPSG:6372（墨西哥官方投影）。别慌，几步就能转成经纬度，顺利做空间连接：

第一步：确认当前数据的坐标系 #

先搞清楚你的数据到底用的是什么投影，不同工具的查看方法：

• QGIS：右键图层 → 属性 → 源标签，直接看CRS（坐标参考系）信息
• Python/Geopandas：用一行代码就能查

import geopandas as gpd

# 加载你的邮编Shapefile
zip_data = gpd.read_file("你的邮编文件路径.shp")
print("当前坐标系:", zip_data.crs)

第二步：转换成经纬度坐标系（WGS84，EPSG:4326） #

咱们目标是转成全球通用的WGS84经纬度（就是GPS用的那种），这样就能和你关注的区域坐标匹配了：

用QGIS可视化操作： #

1. 右键邮编图层 → 导出 → 保存要素为
2. 保存格式选你需要的（比如还是Shapefile），然后在「CRS」下拉框里找到并选择 EPSG:4326 - WGS 84
3. 设置好保存路径，点击确定就行

用Python批量处理： #

# 转换成WGS84经纬度
zip_data_lonlat = zip_data.to_crs(epsg=4326)
# 保存转换后的文件
zip_data_lonlat.to_file("墨西哥城邮编_经纬度版.shp")

第三步：验证转换结果 #

转换完成后，你再看多边形坐标，就会变成类似 POLYGON ((-99.134 19.428, ...)) 这样的经纬度格式了，这时候就可以正常和你关注的区域做空间连接啦！

小提示：如果不知道原数据的EPSG代码怎么办？ #

要是查不到原坐标系，试试墨西哥城常用的几个投影：

• EPSG:32614（UTM 14N）：最常用的，大概率能匹配
• EPSG:6372：墨西哥官方的Lambert投影
  可以先拿其中一个试试转换，看结果的经纬度是否符合墨西哥城的范围（经度大概-99到-100，纬度大概19到20），不对的话再换另一个就行。

内容的提问来源于stack exchange，提问作者J. Yu