Proj4 recjection using R

Question

Я пытаюсь перепрограммировать координаты из WGS84 в зону 53 МГА, проекцию UTM на основе данных GDA94. Я получаю бесконечность в качестве своего результата, что определенно неверно. Я использую пакет 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 ::: project по умолчанию равна градусам. Результаты совпадают с ptransform, если вы конвертируете в радианы.

Answer 2

Где находится пакет proj4?

Попробуйте rgdal, если вы можете установить его:

< ДФ - data.frame ("X" = с (131,1, 131.102, 131.1106, 133,34), "Y" = с (-13,23, -13.243, -13,22, -22,66))

библиотека (rgdal)

## проект ожидает, что матрица, предполагает источник является longlat/WGS84

проект (as.matrix (DF), «+ Рго = UTM + зона = 53 + юг + ellps = GRS80 + towgs84 = 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