Как открыть загрузку DEM USGS в r?

Question

Я загрузил 30-секундные дуги SRTM (топографическая миссия Shuttle Radar) .DEM тайлы из USGS https://dds.cr.usgs.gov/srtm/version2_1/SRTM30/ и хочу открыть их в r, в идеале как растры.

gdalinfo([file path and name].dem)

возвращает сообщение «[[путь к файлу и имя] .dem» не распознан как поддерживаемый формат файла ». Я не очень разбираюсь в gdal, но из сообщения об ошибке следует, что установлена ​​довольно актуальная версия:« выполняется команда » '"/Library/Frameworks/GDAL.framework/Versions/2.1/Programs/gdalinfovation…”. Конечно, это означает, что gdal_translate тоже не работает.

Также не удалось открыть .dem с помощью простого

raster([file path and name])

Как мне это сделать? Я не могу представить, что это уникальный вопрос, но я нашел только вопросы, основанные на том, что DEM уже находится в .tif или другом стандартном растровом формате, или которые связаны с использованием ArcGIS или другого программного обеспечения, отличного от r.

Я не разместил большой файл .dem, но если есть хороший способ, я могу привести воспроизводимый пример, пожалуйста, дайте мне знать.

Answer 1

Файлы .DEM содержат только данные о высоте (в виде простой целочисленной матрицы).Вы можете импортировать матрицу в R, используя readBin("filename.DEM", "integer", size = 2, signed = TRUE, n = NROWS * NCOLS, endian = "big").Но если вы хотите импортировать файл как пространственный растр, вам также необходимо загрузить файлы .hdr из того же места в веб-папке, поскольку они содержат пространственную привязку.Если у вас есть оба файла в одном каталоге, вы можете просто использовать

x = raster("filename.DEM")