如何处理Shapefile中超标经纬度坐标？QGIS/GEE如何适配无效坐标？

开发者阿泽

2026-6-22

解决Shapefile超标坐标映射问题

首先得明确：QGIS和GEE能正常处理这些文件，大概率是因为这些Shapefile用的是投影坐标系（比如UTM、高斯-克吕格），而非WGS84这类地理坐标系。你的解析器错误地把投影坐标当成经纬度来校验，才会出现“数值超标”的误判。

正确的处理步骤

第一步：解析PRJ文件判断坐标系
Shapefile的.prj文件里存储了坐标系的具体参数，你需要先解析这个文件：
- 如果是投影坐标系：直接跳过经纬度范围校验，改用对应投影的坐标范围（比如UTM坐标系的东向坐标通常在100000-999999之间，北向坐标在0-10000000之间）。
- 如果确实是地理坐标系但坐标超标：按以下方式映射到有效范围
第二步：地理坐标系下的超标坐标映射
你之前的取模没奏效，是因为用错了计算逻辑：
- 经度（-180~180）：
  先把经度转换为0_{360范围，再转成-180}180：
```
def normalize_lon(lon):
    lon_mod = (lon % 360 + 360) % 360  # 确保结果在0-360之间
    return lon_mod - 360 if lon_mod > 180 else lon_mod
```
  比如1000度：1000%360=280，280>180，所以280-360=-80，对应西经80度。
- 纬度（-90~90）：
  纬度不能直接取模，因为超过90度后是极地的镜像（比如91度相当于南纬89度）：
```
def normalize_lat(lat):
    lat_mod = lat % 180
    if lat_mod > 90:
        return -(180 - lat_mod) if lat >=0 else (180 - lat_mod)
    elif lat_mod < -90:
        return (-180 - lat_mod) if lat <0 else -(-180 - lat_mod)
    else:
        return lat_mod
```
  比如1000度：1000%180=100，100>90，所以返回-(180-100)=-80（南纬80度）；-1000度：-1000%180=80，直接返回80（北纬80度）。
第三步：验证数据
用QGIS打开问题Shapefile，右键图层→属性→信息，查看“坐标系”字段，确认是投影还是地理坐标系，再对应调整你的解析逻辑。