UTM转经纬度异常：R语言terra包转换结果与预期不符求助

你的转换结果错误的核心原因是错误指定了原始坐标的坐标系。你提供的X=791423.36、Y=6341151.30并非UTM 31N（EPSG:32631）的坐标，而是法国官方坐标系**Lambert-93（EPSG:2154）**的坐标，这也是法国Salindres地区常用的坐标系。

当你用UTM 31N去解析Lambert-93坐标时，必然得到错误的经纬度结果（对应北欧区域的57.12°N、7.81°E），而非预期的Salindres地区坐标。

解决方案 #

将原始坐标的CRS指定为EPSG:2154，再转换到WGS84经纬度（EPSG:4326）即可得到正确结果：

library(terra)
# 用Lambert-93（EPSG:2154）定义原始坐标的CRS
v <- vect(as.matrix(cbind(791423.36, 6341151.30)), crs = "EPSG:2154")
# 转换到WGS84经纬度
y <- project(v, "EPSG:4326")
# 查看转换后的坐标值
geom(y)

执行上述代码后，输出的经纬度将为：

x        y
1 4.1397 44.3608

与你预期的4.14（经度）和44.36（纬度）完全吻合。

补充说明 #

• Lambert-93是法国本土的官方投影坐标系，多数法国境内的地理数据会采用该坐标系，而非UTM。
• 若你需要快速判断坐标所属坐标系，可参考范围：Lambert-93的X范围约在600000到1200000之间，Y范围约在6000000到7200000之间，与你提供的坐标完全匹配。

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