5

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

我基本上有两个问题：1）有人看到我的代码有任何问题吗？和 2) 是否有人熟悉将 UTM 转换为 lat/long 并仅在 Rgooglemaps 包中使用 UTM？

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

数据：

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

结果

> 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”包来创建一些热图和内核密度图。

4

2 回答 2

6

我正在使用以下代码将 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
}
于 2016-04-10T07:55:25.690 回答
1

我认为您只需要交换 X 和 Y--Northing 为 Y，Easting 为 X。原始代码在此更改中对我有用。

于 2020-05-12T12:13:26.993 回答