Мне известно о загрузке данных матрицы высот с использованием пакетов raster
и elevatr
. Однако эти пакеты обеспечивают разрешение SRTM 90 м.
Я хочу загрузить данные матрицы высот высокого разрешения, в идеале разрешение SRTM 30 м для ~ 120 стран. У меня есть буферная зона в 100 км от береговой линии, и мне нужны только данные для этих районов.
Существуют ли какие-либо пакеты или функции, которые могли бы предоставлять функции, аналогичные пакетам raster
и elevatr
для данных SRTM с разрешением 30 м?
Мой рабочий процесс в R для текущей загрузки 90 м Данные SRTM, как показано ниже:
#setting projection
prj_dd <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
#splitting buffered data into 2 rows i.e. two countries
df_1 <- buffer_100km%>%dplyr::slice(1:2)
##download elevation data
elevation_df_hi <- get_elev_raster(df_1, z = 14, prj = prj_dd, clip ="locations")
# on cluster download ~5h
#merging DEMs still running ~18 hours running to present
#Get country geometry
ABW <- raster::getData("GADM",
country = "ABW",
level=0)
AIA<- raster::getData("GADM",
country = "AIA",
level=0)
#setting as sf objects
ABW<-st_as_sf(ABW)
AIA<-st_as_sf(AIA)
#merging data
df_1_country<-rbind(ABW, AIA)
#mapping country and downloaded elevation data
elev_map_df_1<-tm_shape(df_1_country)+
tm_graticules()+
tm_polygons()+
tm_shape(elevation_df_hi)+
tm_raster()+
tm_facets(by = "NAME_0")