我遇到一个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 ?
谢谢!