Как читать файлы .HGT в R

Я пытаюсь прочитать данные о высоте из NASA, хранящиеся в файлах .HGT, в R. Я проверил rgdal библиотеку, но, очевидно, она не читает этот тип файла.

пс. Сообщество Stack Overflow показывает, как читать этот тип файла в Python и C ++, но я искал чистый Раствор R.

Дополнительная информация о топографических данных в файлах .HGT:

В 2014 году топографические данные, полученные в ходе миссии НАСА по радиолокационной топографии шаттла (SRTM), были опубликованы во всем мире. Разрешение данных SRTM для регионов за пределами США составляет 1 угловую секунду или около 30 метров (98 футов). Вы можете прочитать дополнительную информацию здесь и загрузить данные здесь.

В настоящий момент доступны данные по этим регионам мира: введите описание изображения здесь

Вот описание с www2.jpl.nasa.gov/srtm/faq.html:

Файлы данных SRTM имеют такие имена, как N34W119.hgt. Что означают буквы и цифры и что такое формат .hgt?

Каждый файл данных покрывает блок земной поверхности размером один градус широты и один градус долготы. Первые семь символов обозначают юго-западный угол блока, а N, S, E и W относятся к северу, югу, востоку и западу. Таким образом, файл N34W119.hgt охватывает 34–35 северной широты и 118–119 западной долготы (этот файл включает центр Лос-Анджелеса, Калифорния). Расширение имени файла .hgt просто означает высоту слова, означающую высоту. Это НЕ тип формата. Эти файлы представлены в необработанном формате (без заголовков и без сжатия), 16-битных целых числах со знаком, высота, измеренная в метрах над уровнем моря, в географической проекции (массив широты и долготы), с пустыми данными, обозначенными -32768. Международные файлы с 3 угловыми секундами содержат 1201 столбец и 1201 строку данных с общим размером файла 2 884 802 байта (= 1201 x 1201 x 2). Файлы США в 1 угловую секунду содержат 3601 столбец и 3601 строку данных с общим размером файла 25 934 402 байта (= 3601 x 3601 x 2). Для получения дополнительной информации прочтите текстовый файл SRTM_Topo.txt по адресу http://edcftp.cr.usgs.gov/pub/data/srtm/Readme.html.


r binary gis elevation topography
person rafa.pereira    schedule 28.08.2015    source источник
comment
Я могу открыть эти файлы с помощью raster package.   -  person    schedule 28.08.2015
comment
то же самое. hgt <- raster("N38W112.hgt") ; image(hgt)   -  person hrbrmstr    schedule 28.08.2015

Ответы (1)


arrow_upward
5
arrow_downward

Ответ довольно прост, используя пакет raster (благодаря комментариям @Pascal и @hrbrmstr). Также должен быть установлен пакет rgdal.

# Load libraries
  library(raster)
  library(rgdal)

# read file
  elevation <- raster("S23W044.hgt") 
# view image
  image(elevation)
person rafa.pereira    schedule 29.08.2015