Question

我在SWEREF 99 TM中有一些协调，无法使用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

我设法使用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')))

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

1 个答案: