使用r重新投影ProJ4

Question

我正在尝试将基于GDA94基于基于GDA94基于的UTM投影的MGA区域53重新投影。我的结果是无限的，这绝对是不正确的。我正在使用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

> 

这里出了什么问题？

Answer 1

问题是Ptransform期望弧度，而不是学位。函数proj4 :::项目默认为学位。如果您转换为弧度，则Ptransform的结果是相同的。

Answer 2

ProJ4软件包从哪里获得？

如果可以安装rgdal，请尝试：

df < - data.frame（“ x” = C（131.1，131.102，131.1106，133.34），“ y” = C（-13.23，-13.243，-13.22，-22.66））

图书馆（rgdal）

##项目期望矩阵，假设源为longlat/wgs84

项目（as.matrix（df），“ +proj = utm +Zone = 53 +South +Ellps = GRS80 +TOWGS84 = 0,0,0,0,0,0,0,0,0,0 +单位= m +no_defs”）

     [,1]    [,2]

[1,] 77177.18 8534132

[2,] 77416.79 8532695

[3,] 78310.75 8535258

[4,] 329440.68 7493165