我有一个包含 9,000 多个 UTM 坐标的 .csv 文件，我想将其转换为十进制，但遇到了一些麻烦。我搜索了此处和其他地方发布的几篇帖子，但似乎找不到将我的 UTM 集转换为可用且准确的纬度/经度的解决方案。

我基本上有两个问题：1）有人看到我的代码有任何问题吗？ 2）是否有人熟悉放弃将 UTM 转换为纬度/经度并仅在 Rgooglemaps 包中使用 UTM 的情况？

以下是我的代码和数据的一些示例：

Data:

>head(utm)
-Northing Easting
1  4236576  615805
2  4236576  615805
3  4236576  615805
4  4236576  615805
5  4236576  615805
6  4236576  615805

到目前为止的代码：

utm <- read.csv(file="utm.csv", header=TRUE, sep=",")
library(rgdal)
utm <- utm[complete.cases(utm),]
utm1 <- data.frame(x=utm$Northing,y=utm$Easting) 
coordinates(utm1) <- ~x+y 
class(utm1)
proj4string(utm1) <- CRS("+proj=utm +zone=10 +datum=WGS84 +units=m +ellps=WGS84") 
utm2 <- spTransform(utm1,CRS("+proj=longlat +datum=WGS84"))

Results

> head(utm2)
SpatialPoints:
             x        y
[1,] -91.08516 4.727323
[2,] -91.08516 4.727323
[3,] -91.08516 4.727323
[4,] -91.08516 4.727323
[5,] -91.08516 4.727323
[6,] -91.08516 4.727323
Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84 +ellps=WGS84
+towgs84=0,0,0 

所以，我得到了一些输出，但没有得到合理的输出。我在这里缺少什么吗？另外，就其价值而言，我计划使用“Rgooglemaps”包来创建一些热图和核密度图。

我使用以下代码将 UTM 转换为纬度/经度。它适用于伦敦地区

wgs84 = "+init=epsg:4326"
bng = '+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 
+ellps=airy +datum=OSGB36 +units=m +no_defs'

ConvertCoordinates <- function(easting,northing) {
out = cbind(easting,northing)
mask = !is.na(easting)
sp <-  sp::spTransform(sp::SpatialPoints(list(easting[mask],northing[mask]),proj4string=sp::CRS(bng)),sp::CRS(wgs84))
out[mask,]=sp@coords
out
}