Question

我从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])

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

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

Answer 1

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

x = raster("filename.DEM")

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

1 个答案: