У меня проблема с шейп-файлом, для которого мне нужно изменить проекцию, но она работает не совсем так, как мне хотелось бы.
Вот шейп-файл, который я хотел бы преобразовать:
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 ...
Итак, мой вопрос: есть ли способ изменить proj на LAEA , но с координатами DMS?
Спасибо!