我正在尝试将坐标从 WGS84 重新投影到 MGA Zone 53,这是一个基于 GDA94 数据的 UTM 投影。我的结果是无穷大,这绝对是不正确的。我正在使用 R 的 proj4
像这样打包:
> library(proj4)
> df <- data.frame("x" = c(131.1, 131.102, 131.1106, 133.34), "y" = c(-13.23, -13.243, -13.22, -22.66))
> df
x y
1 131.1000 -13.230
2 131.1020 -13.243
3 131.1106 -13.220
4 133.3400 -22.660
> ptransform(data = df, src.proj = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", dst.proj = "+proj=utm +zone=53 +south +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")
$x
[1] Inf Inf Inf Inf
$y
[1] Inf Inf Inf Inf
$z
[1] 0 0 0 0
>
这里出了什么问题?
最佳答案
问题是 ptransform 需要弧度,而不是度数。函数 proj4:::project 默认为度。如果转换为弧度,结果与 ptransform 相同。
关于r - 使用 R 的 Proj4 重投影,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/3984385/