使用R- SWEREF将坐标转换为WGS84

Question

我有一些在SWEREF 99TM中协调的，我不能使用R中的proj4包转换成WGS84。

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

我设法使用rgdal中的SpatialPoints和spTransform转换坐标：

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')))