R:如何将EPSG 25832坐标转换为EPSG 4326坐标

Question

我有一个csv文件，其中三列都是value x y格式。我知道x和y指的是EPSG 25832。我需要将这些坐标转换为EPSG 4326，因为当我在shapefile上调用summary()来绘制数据时，我得到了这一行：

proj4string :
[+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0]

WGS84对应于EPSG 4326，对吗？

通过阅读其他答案，我了解到我需要使用rgdal包中的spTransform函数。然而，我没有找到任何关于函数用法的详细解释。请帮帮我！

文档对我来说很神秘(我对R和空间数据不熟悉)，所以它帮不了我。

编辑：添加dput(head(data))的输出

structure(list(Value = c(10L, 9L, 17L, 13L, 10L, 6L), X = c(687199.0608, 
687199.0608, 687199.0608, 687199.0608, 687199.0608, 687199.0608
), Y = c(4928179.721, 4928179.721, 4928179.721, 4928179.721, 
4928179.721, 4928179.721)), .Names = c("Value", "X", "Y"), row.names = c(NA, 
-6L), class = c("tbl_df", "tbl", "data.frame"))

Answer 1

这应该就行了！

问题是..。你需要用CRS定义两个坐标系统..。这些都是由他们的EPSG代码定义的。然后，您只需使用spTransform :)来转换它们，但这当然只有在为具有坐标和proj4string的数据分配正确的坐标系时才有可能。

library(rgdal)
library(data.table)
d <- structure(list(Value = c(10L, 9L, 17L, 13L, 10L, 6L), X = c(687199.0608, 
687199.0608, 687199.0608, 687199.0608, 687199.0608, 687199.0608
), Y = c(4928179.721, 4928179.721, 4928179.721, 4928179.721, 
4928179.721, 4928179.721)), .Names = c("Value", "X", "Y"), row.names = c(NA, 
-6L), class = c("tbl_df", "tbl", "data.frame"))

d = as.data.table(d)
d = d[,.(X,Y)]
coordinates(d) <- c("X","Y")
proj4string(d) <- CRS("+init=epsg:25832") 
CRS.new <- CRS("+init=epsg:4326") # WGS 84
dnew <- spTransform(d, CRS.new)