我试图将存储在.HGT文件中的海拔数据读入R。我已经检查了rgdal库,但显然它没有读取这种类型的文件。
ps。Stack溢出社区展示了如何在Python和C++中读取这个文件的typee,但是我正在寻找一个纯R解决方案。
有关.HGT文件中地形数据的更多信息:
2014年,NASA的航天飞机雷达地形任务(SRTM)生成的地形数据在全球发布。美国以外地区的SRTM数据分辨率为1弧秒,约30米(98英尺)。您可以阅读更多信息这里并下载数据这里。
目前,世界上这些区域的数据如下:

以下是www.jpl.nasa.gov/srtm/faq.html中的描述:
SRTM数据文件的名称类似于"N34W119.hgt“。字母和数字所指的是什么?".hgt“格式是什么? 每个数据文件覆盖地球表面的一个经度分块的一个纬度。前7个字符表示地块的西南角,N、S、E和W分别指北、南、东和西。因此,"N34W119.hgt“文件涵盖北纬34至35度和西经118至119度(该文件包括加利福尼亚州洛杉矶市中心)。文件扩展名".hgt“只是表示”高度“一词,意思是高度。它不是格式类型。这些文件采用“原始”格式(没有标头,没有压缩),16位有符号整数,海拔以米为单位,在“地理”(纬度和经度阵列)投影中,数据空洞由-32768表示。国际3弧秒文件有1201列和1201行数据,总文件大小为2,884,802字节(=1201x1201x2).美国1弧秒文件有3601列和3601行数据,总文件大小为25,934,402字节(= 3601 x 3601x2)。有关更多信息,请阅读http://edcftp.cr.usgs.gov/pub/data/srtm/Readme.html上的文本文件“http://edcftp.cr.usgs.gov/pub/data/srtm/Readme.html”
发布于 2015-08-29 18:38:42
使用raster包回答非常简单(感谢@Pascal和@hrbrmstr的注释)。还必须安装rgdal包。
# Load libraries
library(raster)
library(rgdal)
# read file
elevation <- raster("S23W044.hgt")
# view image
image(elevation)https://stackoverflow.com/questions/32268234
复制相似问题