更改shapefile投影

时间:2019-04-04 10:47:15

标签: r

我遇到一个shapefile问题,需要为其更改投影,但是它不能完全按照我的意愿工作。

这是我要转换的shapefile:

library(rgdal)
tmp <- tempdir()
url <- "https://www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_admin_0_countries.zip"
file <- basename(url)
download.file(url, file)
unzip(file, exdir = tmp )
world <- readOGR(dsn = tmp, layer = "ne_10m_admin_0_countries")

实际上,我只想要欧洲:

Europe <- c("Northern Europe", "Western Europe", "Southern Europe", "Eastern Europe")
europe <- world[world$SUBREGION %in% Europe,]

我想将投影更改为LAEA:

europe_laea <- spTransform(europe, CRS("+proj=laea")) 

这可行,但是坐标是UTM,尽管我要DMS。即使我将CRS("+proj=laea")更改为CRS("+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"),坐标仍然在UTM中。

所以我的问题是:是否可以通过DMS坐标将项目更改为LAEA but

谢谢!

0 个答案:

没有答案