Преобразование координат с помощью R - SWEREF в WGS84

Question

У меня есть некоторые координаты в SWEREF 99 TM, которые я не могу преобразовать в WGS84 с помощью пакета proj4 в R.

temp <-  data.frame (x = c(598223, 598812, 598824, 598232, 597614, 597629), 
        y = c(7095460, 7095426, 7094827, 7094227, 7094821, 7095433))


ptransform(temp, 
    src.proj = '+proj=utm +zone=33 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs ',
    dst.proj = '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs', silent=TRUE)

Обратите внимание, что выходные данные не эквивалентны координатам в Швеции - широта эквивалентнавозле экватора.

Answer 1

Мне удалось преобразовать координаты, используя SpatialPoints и spTransform из rgdal:

temp <- data.frame (x = c(598223, 598812, 598824, 598232, 597614, 597629), y = c(7095460, 7095426, 7094827, 7094227, 7094821, 7095433))

temp <- SpatialPoints (temp, proj4string = CRS ('+proj=utm +zone=33 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs'))

temp <- as.data.frame (spTransform(temp, CRS('+proj=longlat +datum=WGS84 +no_defs')))