如何在 r 中打开 USGS DEM 下载？答案

【问题标题】：How does one open a USGS DEM download in r?如何在 r 中打开 USGS DEM 下载？
【发布时间】：2019-06-18 11:54:54
【问题描述】：

我从 USGS https://dds.cr.usgs.gov/srtm/version2_1/SRTM30/ 下载了 30 弧秒 SRTM（航天飞机雷达地形任务）.DEM 切片，并希望在 r 中打开它们，最好是作为栅格。

gdalinfo([file path and name].dem)

返回消息“'[文件路径和名称].dem'不被识别为受支持的文件格式。”我对 gdal 不太了解，但错误消息暗示安装了一个相当最新的版本：“运行命令'”/Library/Frameworks/GDAL.framework/Versions/2.1/Programs/gdalinfo”...”当然这意味着gdal_translate 也失败了。

用普通格式打开 .dem 也不起作用

raster([file path and name])

我应该怎么做？我无法想象这是一个独特的问题，但我发现的问题仅基于 DEM 已经采用 .tif 或其他标准栅格格式，或者涉及使用 ArcGIS 或 r 以外的其他软件。

我没有发布大的 .dem 文件，但如果有好的方法我可以提供一个可重现的示例，请告诉我。

【问题讨论】：

标签： r r-raster rgdal

【解决方案1】：

.DEM 文件仅包含高度数据（作为简单的整数矩阵）。您可以使用 readBin("filename.DEM", "integer", size = 2, signed = TRUE, n = NROWS * NCOLS, endian = "big") 将矩阵导入 R。但是，如果要将文件作为空间栅格导入，还需要从同一 Web 文件夹位置下载 .hdr 文件，因为这些文件包含空间参考。一旦将两个文件都放在同一个目录中，您就可以简单地使用

x = raster("filename.DEM")

【讨论】：