2. DEM数据获取与预处理实战

上一节我们理解了DEM的栅格本质与高程含义，这一节要把“数据”真正变成“能用的数据”。一个常见误区是：从公开数据源下载到的DEM压缩包，解压后直接拖进GIS就开始算流向。实际上，原始下载产品往往带着投影不统一、分辨率过粗、范围过大、存在洼地等问题，直接使用会让后续的流向与汇流分析产生大量伪路径。本节以“为一个小流域准备可用DEM”为目标，把下载、投影转换、重采样、裁剪、填洼这条链路讲清楚。

一、数据源选择：先明确分辨率与覆盖需求

公开DEM数据源中，USGS EarthExplorer 提供 SRTM、ASTER GDEM、Copernicus DEM 等产品，OpenTopography 则提供经过处理的多种全球与区域DEM，并支持按范围裁剪下载。选择时先回答两个问题：研究区需要多高的水平分辨率，以及垂直精度能否满足坡度、汇流分析。

以一个小流域（例如面积约 50 km² 的山区小流域）为例，若只做宏观地形起伏判读，30 m 分辨率足够；若要识别坡面细沟与小型沟道，则倾向 12 m 或更高。需要注意，Copernicus DEM 的 GLO-30 为 30 m 级，GLO-90 为 90 m 级；SRTM 常见为 1 弧秒（约 30 m）。分辨率数字是格网间距，不是精度承诺，植被覆盖区的高程仍可能偏差数米。

下载时建议直接按流域外接矩形略加缓冲（例如外扩 1～2 km）选取范围，避免后续裁剪把边界处的汇流路径截断。若数据源支持选择投影与格式，优先选地理坐标（经纬度）的 GeoTIFF，便于后续统一处理。

二、投影转换：让“米”成为分析单位

下载的DEM多为地理坐标系（如 WGS 84，EPSG:4326），其像元大小以度为单位。直接在这种数据上计算坡度、流向和汇流面积，会因经度方向的实际地面距离随纬度变化而产生系统性偏差。因此第一步是把DEM转换到投影坐标系，使像元在东西、南北方向都对应近似相等的实际距离。

对小流域，常用 UTM 投影（按所在带选择 EPSG 编码，如中国东部某流域可能落在 EPSG:32650）。转换时使用“投影栅格”工具，重采样方法选双线性插值或三次卷积，不要用最近邻——最近邻会保留原始高程值但产生阶梯状像元，对连续地形不利。转换后检查像元大小是否接近设计值（如 30 m），并确认范围没有明显偏移。

三、重采样与裁剪：统一格网、锁定范围

若手头数据分辨率与目标不一致（例如拿到 90 m 数据但需要 30 m 格网），需要重采样。这里要清醒：重采样不能创造真实细节。把 90 m 升采样到 30 m 只是插值出更密的格网，地形细节仍是 90 m 的水平；反之降采样会平滑掉小沟谷。因此重采样应服务于“统一格网”，而不是“提升精度”。

裁剪则用流域边界（或外接矩形）提取子区。推荐用“按掩膜提取”，掩膜为流域面要素；若暂时没有精确边界，可先用外接矩形裁剪，待边界确定后再精裁。裁剪后务必检查：行列数、像元大小、NoData 值是否合理，边界处是否出现整行整列的 NoData 导致后续填洼异常。

四、填补洼地：让水流“走得出去”

DEM中的洼地分两类：真实洼地（如岩溶漏斗、冰川湖盆）和伪洼地（由数据噪声、插值误差、低分辨率造成）。基于D8等单流向算法时，水流进入洼地后无更低邻域可去，会形成断流或错误汇聚。因此填洼是流向分析前的必要步骤。

常用做法是运行“填洼”工具，它把每个洼地抬升到其出口最低点的高程，形成平坦或微倾斜面。操作时注意两点：一是设置合理的填洼阈值，例如只填平深度不超过 2 m 的洼地，避免把真实湖盆也抹平；二是填洼后对比填洼前后的高程差栅格，检查被修改区域是否集中在数据边缘或噪声区，若大片区域被抬升，说明数据质量或分辨率可能有问题。

填洼不是终点。填洼后应再检查一次流向栅格：是否仍有大量像元指向自身或形成短小闭合环；若有，可能是边界NoData或平坦区处理不当，需要结合“平坦区流向处理”进一步修正。

五、一份可执行的预处理检查清单

为便于复现，下面给出一份针对小流域DEM的检查清单，字段与判断依据均围绕本节任务设计，可直接对照操作记录使用。

完成上述步骤后，得到的DEM才具备进入流向计算的条件：投影一致、格网统一、范围完整、洼地受控。需要强调的是，本节示例中的 30 m、2 m 阈值、UTM 带号均为教学假设，实际项目应根据流域尺度、数据源精度和分析目的调整，并以填洼前后对比与流向检查结果作为判断依据，而不是机械套用参数。