如何在 r 中打开 USGS DEM 下载?

How does one open a USGS DEM download in r?

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

gdalinfo([file path and name].dem)

returns 消息“'[文件路径和名称].dem' 未被识别为受支持的文件格式。”我对 gdal 了解不多,但错误消息暗示安装了相当新的版本:“运行 command '"/Library/Frameworks/GDAL.framework/Versions/2.1/Programs/gdalinfo” ……”当然这意味着 gdal_translate 也失败了。

用普通

打开 .dem 也没有用
raster([file path and name])

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

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

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

x = raster("filename.DEM")