Proj4 Reprojektion R unter Verwendung von
-
09-10-2019 - |
Frage
Ich versuche, Koordinaten von WGS84 zu MGA Zone 53, eine UTM-Projektion auf der Grundlage des GDA94 Datum zu projizieren. Ich bekomme unendlich wie mein Ergebnis, was definitiv falsch ist. Ich bin mit R proj4
Paket wie folgt:
> 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
>
Was läuft falsch hier?
Lösung
Das Problem ist, dass Ptransform erwartet Radiant nicht Grad. Die Funktion proj4 ::: Projektvorgaben in Grad. Die Ergebnisse sind gleich mit Ptransform, wenn Sie in Radiant umwandeln.
Andere Tipps
Wo ist das proj4 Paket erhalten aus?
Versuchen rgdal, wenn Sie es installieren können:
df <- data.frame ( "x" = c (131,1, 131,102, 131,1106, 133.34), "y" = c (-13,23, -13,243, -13,22, -22,66))
Bibliothek (rgdal)
## Projekt eine Matrix erwartet nimmt Quelle ist longlat / WGS84
Projekt (as.matrix (df), "+ proj = utm + zone = 53 + Süden + ellps = GRS80 + towgs84 = 0,0,0,0,0,0,0 + Einheiten = m + no_defs" )
[,1] [,2]
[1] 77.177,18 8534132
[2] 77.416,79 8532695
[3] 78.310,75 8535258
[4] 329.440,68 7493165