基于DEM的ArcGIS水文分析-河网和流域的提取

基于DEM的ArcGIS水文分析—河网和流域的提取在地理信息科学与水文研究领域，数字高程模型（DEM）因其能够精确反映地形起伏特征，成为水文分析的重要基础数据。利用DEM进行河网提取与流域划分，对于水资源管理、洪水预报、生态保护等工作具有至关重要的意义。ArcGIS作为一款功能强大的地理信息系统平台，提供了一系列用于水文分析的工具，能够帮助我们从DEM数据中有效地提取河网结构并划定流域边界。本文将详细阐述基于DEM在ArcGIS中进行河网和流域提取的完整流程与关键技术点。一、数据准备与预处理任何分析工作的成功都离不开高质量的基础数据。DEM数据的质量直接影响后续水文分析结果的准确性。因此，在进行水文分析之前，对DEM数据的准备与预处理是首要步骤。首先，我们需要获取研究区域的DEM数据。这些数据可以来源于国家地理信息中心、相关科研机构或商业数据提供商。DEM数据的分辨率选择需根据研究区域的大小和分析精度要求来决定。一般而言，小区域精细分析宜采用较高分辨率DEM，而大范围宏观分析则可选用较低分辨率数据。获取DEM数据后，首要任务是对其进行质量检查。这包括检查数据是否存在空值、异常值，以及是否有明显的噪声或数据偏移。在ArcGIS中，可以通过“数据管理工具”下的“镶嵌至新栅格”或“复制栅格”等工具确保数据的完整性。同时，需确认DEM数据的坐标系统和投影信息，确保其与研究区域的其他数据保持一致，避免因投影不统一导致的分析误差。二、DEM预处理——核心步骤原始DEM数据往往存在一些不符合水文分析假设的地形特征，例如洼地和尖峰。这些特征会影响水流方向的正确计算，因此必须进行预处理。（一）填洼处理自然界中，水流总是从高处流向低处。但原始DEM中可能存在一些封闭的洼地，这些洼地会导致水流在此汇集而无法流出，形成“死水”区域，这与实际水文过程不符。填洼（Fill）操作的目的就是消除这些洼地，使DEM表面形成一个连续的、无洼地的“平滑”表面，确保水流能够顺畅流动。在ArcGIS的“空间分析工具”->“水文分析”工具箱中，提供了“填洼”工具（Fill）。该工具通过抬高洼地内栅格的高程值，使其与洼地出口的高程一致，从而消除洼地。在使用时，需要注意“Z限制”参数的设置。“Z限制”定义了允许填充的最大深度，如果设置为0，则所有洼地都将被填充；若设置为一个具体值，则只填充深度小于该值的洼地。在实际应用中，若对研究区域的地形特征不十分了解，通常建议先采用默认的0值进行完全填充，以确保后续分析的顺利进行。（二）削峰处理（可选）与洼地相对的是“尖峰”或“伪山脊”，即一些孤立的、高程异常高的栅格单元，它们可能导致水流方向计算出现错误。削峰（Smooth）处理可以平滑这些异常值，使地形更符合自然形态。但削峰操作并非总是必需的，需根据DEM数据的质量和具体研究需求来决定是否执行。ArcGIS中并无专门的“削峰”工具，但可以通过“邻域分析”中的“平均值”或“中值”等工具进行平滑处理，以达到类似效果。三、水流方向计算在经过预处理的DEM基础上，下一步是计算水流方向（FlowDirection）。水流方向是指每个栅格单元的水流将流向其相邻的哪个栅格单元。ArcGIS中采用D8算法（Deterministic8-directionmethod）来计算水流方向。该算法认为每个栅格单元的水流方向只可能流向其8个相邻栅格单元（上、下、左、右及四个对角线方向）中坡度最陡的那个单元。使用“水文分析”工具箱中的“流向”工具（FlowDirection），可以生成一个水流方向栅格。该栅格的每个像元值代表了该像元水流流出的方向，通常用1、2、4、8、16、32、64、128这八个数值分别代表东、东南、南、西南、西、西北、北、东北八个方向。水流方向栅格是后续所有水文分析的基础，其计算的准确性直接关系到整个分析结果的可靠性。四、汇流累积量计算汇流累积量（FlowAccumulation）是指流经每个栅格单元的上游汇水区域的栅格数量（或面积）。它反映了该栅格单元汇聚的水流总量。汇流累积量的大小与河网的形成密切相关，汇流累积量越大的区域，越有可能发育成河流。在ArcGIS中，利用“水文分析”工具箱中的“流量”工具（FlowAccumulation），以水流方向栅格作为输入，可以计算得到汇流累积量栅格。该栅格的每个像元值表示流经该像元的上游栅格数量。理论上，汇流累积量为0的区域表示该区域没有上游来水，而汇流累积量值较大的区域则对应着河流的位置。五、河网提取有了汇流累积量栅格，我们就可以通过设定一个合适的阈值来提取河网。这个阈值代表了能够形成河流的最小汇流累积量。阈值的选择是河网提取的关键，直接影响提取出的河网密度和形态。（一）确定阈值阈值的确定并没有一个统一的标准，需要结合研究区域的地形特征、DEM分辨率以及实际的水文情况进行综合判断。一种常用的方法是通过观察汇流累积量的频率分布，选择一个能够区分河流与非河流区域的拐点值作为阈值。在ArcGIS中，可以通过“统计分析”工具查看汇流累积量的最大值、最小值、平均值等统计信息，并结合经验进行初步设定。此外，也可以参考研究区域已有的河流数据，通过对比调整阈值，使提取的河网与实际河网尽可能吻合。（二）提取河网栅格使用“空间分析工具”->“地图代数”->“栅格计算器”（RasterCalculator），输入表达式“Con(Flow_Accumulation>=阈值,1)”，即可将汇流累积量大于等于阈值的栅格赋值为1（表示河流），其余为NoData，从而得到河网栅格。（三）河网矢量化为了便于后续的分析和应用（如计算河长、河网密度，或与其他矢量数据叠加分析），通常需要将提取出的河网栅格转换为矢量线要素。ArcGIS的“水文分析”工具箱中的“栅格河网转矢量”工具（StreamtoFeature）可以实现这一转换。该工具能够将河网栅格转换为具有拓扑关系的矢量线要素，并可以选择是否简化线要素。六、流域提取流域，也称为集水区，是指由分水岭所包围的、汇集水流到某一出口断面的区域。基于已提取的河网，我们可以进一步划定各个子流域或整个研究区域的流域边界。（一）确定流域出口点流域的提取通常需要指定一个或多个出口点（Outlets）。出口点可以是河流的某一特定位置（如水文站所在位置）、河口或研究区域的边界点。在ArcGIS中，可以通过“创建随机点”工具在河网上生成点，或手动数字化出口点。（二）分水岭计算利用“水文分析”工具箱中的“分水岭”工具（Watershed），以水流方向栅格和出口点数据作为输入，即可计算得到每个出口点所控制的流域范围。该工具会生成一个栅格数据，每个栅格的值代表其所属的流域编号。（三）流域矢量化与后处理与河网矢量化类似，得到的流域栅格也可以通过“栅格转面”工具（RastertoPolygon）转换为矢量面要素。转换后，还可以对矢量流域边界进行平滑、融合等后处理操作，以提高其美观性和实用性。七、关键技术点与经验分享1.DEM分辨率的影响：高分辨率DEM能提供更精细的地形细节，提取的河网也更接近实际，但数据量较大，计算耗时。低分辨率DEM则适用于宏观分析。2.阈值的敏感性：汇流累积量阈值的选择对河网提取结果影响显著。阈值过小，会提取出过多的细小支流，甚至噪声；阈值过大，则可能遗漏一些重要的小河流。建议通过多次试验，并结合实地考察或高分辨率遥感影像进行验证。3.洼地填充的合理性：过度填充可能会改变原有的地形特征，特别是在一些喀斯特地貌区，天然洼地可能是重要的水文单元。因此，在填洼时，若条件允许，可结合区域地质水文知识，合理设置“Z限制”。4.结果的验证：提取的河网和流域边界需要与实际情况进行对比验证。可以利用已有的地形图、水系图或高分辨率卫星影像进行目视检查，或通过计算一些水文参数（如河网密度、流域面积）与已知数据进行比较。八、结论基于DEM的ArcGIS水文分析为河网提取和流域划分提供了一套高