SRTM 30 м разрешение скачать по стране функции или пакета в R? - PullRequest
0 голосов
/ 07 января 2020

Мне известно о загрузке данных матрицы высот с использованием пакетов 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")
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...