r - 如何在 R 中读取 .HGT 文件

Question

我正在尝试将存储在 .HGT 文件中的 NASA 海拔数据读取到 R 中。我已经检查过 rgdal库，但显然它不读取这种类型的文件。
附: Stack Overflow 社区在 Python 中展示了如何读取这种类型的文件和 C++ ，但我一直在寻找纯 R 解决方案。
有关 .HGT 文件中地形数据的更多信息:
2014年，美国宇航局航天飞机雷达地形任务(SRTM)生成的地形数据在全局发布。美国以外地区的 SRTM 数据分辨率为 1 弧秒，或约 30 米(98 英尺)。您可以阅读更多信息here并下载数据 here .
目前，这些数据可用于世界这些地区:

以下是来自 www2.jpl.nasa.gov/srtm/faq.html 的描述:

The SRTM data files have names like "N34W119.hgt". What do the letters and numbers refer to, and what is ".hgt" format?

Each data file covers a one-degree-of-latitude by one-degree-of-longitude block of Earth's surface. The first seven characters indicate the southwest corner of the block, with N, S, E, and W referring to north, south, east, and west. Thus, the "N34W119.hgt" file covers latitudes 34 to 35 North and longitudes 118-119 West (this file includes downtown Los Angeles, California). The filename extension ".hgt" simply stands for the word "height", meaning elevation. It is NOT a format type. These files are in "raw" format (no headers and not compressed), 16-bit signed integers, elevation measured in meters above sea level, in a "geographic" (latitude and longitude array) projection, with data voids indicated by -32768. International 3-arc-second files have 1201 columns and 1201 rows of data, with a total filesize of 2,884,802 bytes ( = 1201 x 1201 x 2). United States 1-arc-second files have 3601 columns and 3601 rows of data, with a total filesize of 25,934,402 bytes ( = 3601 x 3601 x 2). For more information read the text file "SRTM_Topo.txt" at http://edcftp.cr.usgs.gov/pub/data/srtm/Readme.html

Answer 1

答案很简单，使用 raster包(感谢@Pascal 和@hrbrmstr 的评论)。 rgdal包也必须安装。

# Load libraries
  library(raster)
  library(rgdal)

# read file
  elevation <- raster("S23W044.hgt") 
# view image
  image(elevation)