如何使用 R 将 WGS84 转换为纬度/经度

Question

我有一个包含 WGS84 坐标的大型数据集，我想将其转换为 Long/Lat 坐标。 WGS84 坐标目前分为自己的列（即“北”“东”），另外还有“区域”和相应的区域。总共有几个区域。

关于将这些列转换为 Lat/Long 并将它们保存到新列中的任何建议？

数据集示例：

Longitude Latitude zone  
233243    6571770   33  
262706    6653520   33  
195348    6573696   33   
256880    6645020   33  
260610    6654042   32  
13799     6505840   33  

Answer 1

您可能需要阅读sp 包、坐标和spTransform 函数。除其他事项外。

对于d 作为数据框：

> d
  Longitude Latitude zone
1    233243  6571770   33
2    262706  6653520   33
3    195348  6573696   33
4    256880  6645020   33
5    260610  6654042   32
6     13799  6505840   33

计划是：通过zone id 拆分成一个列表，然后为每个列表元素创建一个空间点数据框，使用坐标参考系统“+init=epsg:326”+“zone”（例如“ +init=epsg:32633" for zone 33)，然后转换为 epsg:4326 lat-long，然后将地块合并到单个空间点数据框：

需要这些：

> library(sp) ; library(raster)

按区域 ID 创建列表：

> byzone = split(d,d$zone)

在所有部分上，设置坐标、CRS和变换：

> zdll = lapply(byzone, 
    function(zd){
      coordinates(zd)=~Longitude+Latitude
      proj4string(zd)=paste0("+init=epsg:326",zd$zone[1])
      spTransform(zd, CRS("+init=epsg:4326"))
     })

现在加入他们：

> dll = do.call(rbind.SpatialPointsDataFrame, zdll)

如果您现在绘制dll，您将看到经纬度上的点。如果您想要坐标并与原始坐标进行比较，请执行以下操作：

> cbind(d, coordinates(dll))
  Longitude Latitude zone Longitude Latitude
5    233243  6571770   33  4.712098 59.95395
1    262706  6653520   33 10.327202 59.20086
2    195348  6573696   33 10.750120 59.95049
3    256880  6645020   33  9.663825 59.19259
4    260610  6654042   32 10.656096 59.87099
6     13799  6505840   33  6.663177 58.42136

他们在预期的地方吗？