arcgis获取面DEM

### 如何在 ArcGIS 中获取和处理面 DEM 数据

#### 加载 USGS 格式的 DEM 文件并转换为栅格
对于 USGS 格式的 DEM 文件，在 ArcGIS 的操作相对简单。可以利用 `ArcToolBox` 下的 “转换工具”—“转到栅格”—“DEM 转到栅格” 工具来完成这一过程[^1]。

import arcpy
from arcpy import env

env.workspace = "C:/data"
arcpy.DEMToRaster_conversion("input_dem", "output_raster", "FLOAT")

#### 合并多个 DEM 文件
当需要合并若干个 DEM 文件时，可以通过 `ArcToolbox` 提供的功能实现。具体路径为 `Spatial Analysis Tools`—`Extraction`—`Extract by Mask`。需要注意的是，执行此功能前要确认已安装 Spatial Analyst 扩展模块，并且确保用于裁剪的矢量图层为多边形类型以及其空间参考与目标 DEM 一致；如果不一致，则需提前使用 PROJECT 进行转换[^2]。

inMaskData = "mask.shp"
inputRaster = "elevation"
outExtractByMask = ExtractByMask(inputRaster, inMaskData)
outExtractByMask.save("C:/sapyexamples/output/outextmask")

#### 计算坡度
为了基于 DEM 数据计算坡度，首先要保证输入的数据具有合适的坐标系统。如果原始数据仅具备地理坐标系而非投影坐标系，那么应当先将其转换至适合的投影坐标系下再进行后续运算。这一步骤可通过菜单栏中的 【数据管理】—【投影与变换】–【投影栅格】 来达成。接着便能运用 Slope 函数来进行实际的坡度分析工作了[^4]。

inRaster = "C:/data/inraster.tif"
outMeasurement = "DEGREE"
zFactor = 1
outSlope = Slope(inRaster, outMeasurement, zFactor)
outSlope.save("C:/data/outslope01")