rを使用したproj4再注射
-
09-10-2019 - |
質問
私は、WGS84から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
>
ここで何が問題になっているのですか?
解決
問題は、Ptransformが度ではなくラジアンを期待していることです。関数proj4 :::プロジェクトはデフォルトであることがあります。結果は、ラジアンに変換する場合、ptransformで同じです。
他のヒント
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 +ゾーン= 53 +ellps = grs80 +towgs84 = 0,0,0,0,0,0 +unit = m +no_defs"))
[,1] [,2]
[1,] 77177.18 8534132
[2,] 77416.79 8532695
[3,] 78310.75 8535258
[4,] 329440.68 7493165
所属していません StackOverflow