基于SRTM-DEM的甘河小流域精细化提取与特征分析研究

基于SRTM-DEM的甘河小流域精细化提取与特征分析研究一、引言1.1研究背景与意义甘河小流域在区域生态和经济发展中占据着举足轻重的地位。它不仅是区域水资源的重要涵养地，对维持当地的生态平衡、保障生态系统的稳定运行起着关键作用，还是当地居民生活用水以及农业灌溉用水的重要来源，直接关系到当地居民的生产生活和区域农业的可持续发展。然而，要实现对甘河小流域的科学管理和有效保护，深入了解其水文地理信息是必不可少的前提条件。SRTM-DEM（ShuttleRadarTopographyMission-DigitalElevationModel，航天飞机雷达地形测绘使命数字高程模型）数据在流域提取研究中发挥着关键作用。作为一种高精度的地形数据，SRTM-DEM能够精确地反映地表的高程信息，为流域提取提供了坚实的数据基础。利用SRTM-DEM数据进行流域提取，具有传统方法无可比拟的优势。它能够快速、准确地获取流域的地形特征、河网分布以及流域边界等关键信息，极大地提高了研究效率，并且能够避免传统实地测量和手工绘制方法可能出现的人为误差，提高了研究结果的准确性和可靠性。此外，SRTM-DEM数据覆盖范围广泛，能够为不同区域的流域提取研究提供数据支持，具有很强的通用性和实用性。本研究基于SRTM-DEM数据对甘河小流域进行提取研究，旨在深入挖掘甘河小流域的地形和水文特征，为该流域的水文模拟、水资源管理、生态保护以及防洪减灾等工作提供科学依据和数据支持，具有重要的理论和实践意义。通过准确提取甘河小流域的河网和流域边界，能够更精准地掌握流域内的水流路径和汇水区域，为水资源的合理开发利用和优化配置提供决策支持；通过分析流域的地形特征，如坡度、坡向等，能够更好地了解流域内的水土流失情况和生态环境状况，为生态保护和修复提供科学指导；通过对流域水文过程的研究，能够提高对洪水等自然灾害的预测和预警能力，为防洪减灾工作提供有力保障。1.2国内外研究现状在国外，利用SRTM-DEM数据进行流域提取的研究开展得较早，并且取得了一系列丰硕的成果。许多学者运用先进的算法和模型，结合SRTM-DEM数据，对不同类型的流域进行了深入研究。例如，一些学者通过改进地表径流漫流模型，提高了河网提取的精度和效率；还有一些学者利用SRTM-DEM数据与其他地理信息数据相结合，对流域的生态环境、水资源变化等进行了综合分析。然而，这些研究在面对复杂地形和特殊地理条件的流域时，仍然存在一定的局限性，如在山区等地形起伏较大的区域，河网提取的准确性有待进一步提高；在数据处理和分析过程中，对于一些细节特征的把握还不够精准。在国内，随着地理信息技术的快速发展，基于SRTM-DEM数据的流域提取研究也得到了广泛关注和深入开展。众多学者针对不同流域的特点，开展了大量的实证研究，在数据处理、算法优化以及应用拓展等方面取得了显著进展。一些研究通过对比不同分辨率的SRTM-DEM数据对流域提取结果的影响，为数据选择提供了参考依据；还有一些研究将流域提取结果应用于水资源规划、水利工程建设等实际领域，取得了良好的效果。但目前国内研究在不同数据源的融合利用、多尺度分析以及模型的适应性等方面仍存在不足，需要进一步加强研究。综合国内外研究现状，虽然利用SRTM-DEM数据进行流域提取的研究已经取得了一定的成果，但仍存在一些问题和挑战有待解决。例如，如何进一步提高河网提取的精度和可靠性，尤其是在复杂地形条件下；如何更好地融合多源数据，提高流域特征提取的全面性和准确性；如何加强流域提取结果在实际应用中的有效性和针对性等。本研究将在前人研究的基础上，针对甘河小流域的特点，探索更加有效的流域提取方法和技术，以期为该领域的研究提供新的思路和方法，填补相关研究的空白或不足。1.3研究目标与内容本研究旨在基于SRTM-DEM数据，运用先进的地理信息技术和方法，准确、高效地提取甘河小流域的地形和水文特征，为流域的科学管理和可持续发展提供全面、可靠的基础数据和科学依据。具体研究内容包括：数据处理：全面收集与甘河小流域相关的SRTM-DEM数据，并对其进行严格的预处理和格式转换。通过细致的检查和处理，确保数据的质量和完整性，消除数据中的噪声和错误，使其能够满足后续深入分析和处理的要求。地形分析：借助专业的GIS软件，深入分析SRTM-DEM数据，精确提取甘河小流域的关键地形特征，如高程、坡度、坡向等。通过对这些地形特征的深入研究，全面了解流域的地形地貌格局，为后续的流域边界提取和河网提取提供坚实的地形基础。河网提取：依据地形分析的结果，运用科学合理的方法，从SRTM-DEM数据中准确提取甘河小流域的河网信息。详细确定河流的具体流向、精确测量河流的长度，并对河网的分级和结构进行深入分析，全面掌握流域内的水系分布和水流路径。流域边界提取：利用先进的技术手段和算法，基于SRTM-DEM数据精确提取甘河小流域的边界。通过对流域边界的准确界定，明确流域的范围和边界条件，为后续的水文分析和模拟提供准确的边界数据。验证与分析：将提取得到的河网和流域边界与实际情况进行细致的对比和验证，运用多种评估指标和方法，对提取结果的精度进行全面、客观的评价。深入分析提取过程中可能存在的误差来源和影响因素，提出针对性的改进措施和优化方案，不断提高提取结果的准确性和可靠性。1.4研究方法与技术路线本研究主要采用以下方法进行数据获取与处理：通过中国科学院计算机网络信息中心国际科学数据镜像网站等权威渠道，获取覆盖甘河小流域的SRTM-DEM数据。利用ArcGIS软件强大的数据处理功能，对获取的原始数据进行全面的预处理，包括格式转换，将数据转换为适合ArcGIS软件处理的格式；坐标系统统一，确保数据在同一坐标系统下，避免因坐标差异导致的误差；以及数据裁剪，去除与研究区域无关的数据，减少数据量，提高处理效率。运用ArcGIS软件中的水文分析模块，对预处理后的数据进行深入分析。首先，通过填洼处理，消除DEM数据中的洼地，确保水流方向的准确性；然后，基于D8算法，精确计算流向，确定每个栅格单元的水流方向；接着，根据流向数据，准确计算汇流累积量，为河网提取提供关键数据；最后，通过设定合理的汇流累积阈值，成功提取河网，并利用水流追踪算法，准确提取流域边界。本研究的技术路线如下：首先，进行数据获取，收集甘河小流域的SRTM-DEM数据。接着，在ArcGIS软件中对数据进行预处理，包括格式转换、坐标系统统一和裁剪等操作。然后，利用水文分析模块，依次进行填洼、流向计算、汇流累积量计算、河网提取和流域边界提取等步骤。在河网提取过程中，通过不断调整汇流累积阈值，分析不同阈值下的河网提取结果，选择最佳的阈值。在流域边界提取后，将提取结果与实际情况进行对比验证，对精度进行评估。最后，根据验证和评估结果，对提取结果进行优化和完善，得出最终的甘河小流域提取结果。二、研究区概况与数据源2.1甘河小流域概况甘河小流域位于内蒙古自治区东北部，地处大兴安岭山脉东侧，其经纬度范围大致为东经[具体起始经度]-[具体结束经度]，北纬[具体起始纬度]-[具体结束纬度]。该流域呈西北-东南走向，发源于海拔1197米的沃其山麓，最终在黑蒙交界处附近汇入嫩江，是嫩江的重要一级支流。从地形地貌来看，甘河小流域内地形以丘陵和平原为主，地势西北高、东南低。流域内山峦起伏，河谷纵横，山地海拔多在500-1000米之间，相对高差较大，坡度一般在15°-35°之间，局部地区可达40°以上。在河流的长期侵蚀和堆积作用下，形成了较为宽阔的河谷平原，谷底宽为3-4公里，呈梯形或浅槽形河谷。甘河谷坡为20°-40°，河道沿右岸而行，有较多的石砬子，水流较急，河道曲折，呈S形或Ω形接连不断，旧河道和牛轭湖时时可见。从两岸注入甘河的小溪较多，水泡子星罗棋布，砂滩、砂咀出现甚多，一般在0.4-1.0平方公里左右，最大的在老道口附近，长达2.5公里，宽约500米，为砂卵石。河流凹岸因水流冲刷坍塌形成土岸，凸岸则为砂、卵石滩淤积。甘河小流域属于温带大陆性季风气候，夏季受副热带海洋气团影响，冬季受极地大陆气团控制。其主要气候特点为春季干旱少雨，多大风天气，气温回升较快；夏季炎热、降雨集中，6-8月的降水量约占全年降水量的70%-80%，年降水量为500毫米左右；秋季凉爽，降温快，常有早霜；冬季漫长、寒冷，封冻期长达6个月，最大冰厚可达2.0米以上，11月上、中旬封河，4月中、下旬开河。加格达奇站多年平均气温为-0.8℃，最高气温出现在6月份，月平均气温15.0℃，最低气温出现在1月份，月平均气温-23.3℃，最高气温37.3℃（1969年6月21日），最低气温-45.4℃（1980年1月13日），平均无霜期106天。在水文条件方面，甘河全长490公里，河宽50-200米，河深为1-3米，流域面积1.97万平方公里，多年平均流量130立方米/秒。据水文站记载，多年水位平均在98.24米，最高洪水位为102.01米（1988年8月9日测），最低水位为97.56米（1988年4月4日测）；多年平均流量为66.7立方米/秒，最大流量为3200立方米/秒（1988年8月9日测），最小流量为0.43立方米/秒（1985年2月12日测）。水中含盐度极值平均为209毫克当量/升，最大时为260毫克当量/升（1983年7月20日测）；河水总硬度极值从1975年至1984年平均为0.7毫克当量/升，最大为1.03毫克当量/升（1976年2月16日测）；酸碱度（PH）极值多年平均为7.0，最大为8.1（1976年7月11日测）。其主要支流有克一河、阿里河、吉文河、奎勒河等，这些支流的水文特征也各有差异，共同构成了甘河小流域复杂的水系网络。克一河发源于大兴安岭北段东麓鄂伦春自治旗境内伊图里河林业局库布春林场北沟，河长79.2公里，河深0.65米，河宽15.5米，流域面积1358平方公里，流速0.85米/秒，流量9.89立方米/秒，年平均径流量3.12亿立方米，河道平均比降1.9‰，年结冰期约5-6个月；阿里河发源于鄂伦春自治旗东北部伊勒呼里山西南麓，河长106公里（一说124公里），流域面积2150平方公里（一说2183平方公里），多年平均流量16.1立方米/秒，自然落差415米，水能理论蕴藏量1.75万千瓦；奎勒河发源于伦春自治旗中部大兴安岭东麓，河长189公里，流域面积4723平方公里，多年平均流量28.1立方米/秒，自然落差415米。2.2SRTM-DEM数据介绍本研究使用的SRTM-DEM数据是由美国国家航空航天局（NASA）和国家影像与制图局（NIMA）联合实施的航天飞机雷达地形测绘使命获取的数据产品。该数据在2000年2月由“奋进”号航天飞机搭载SRTM雷达系统进行了全球范围的地形测量，其数据覆盖范围广泛，能够为全球大部分地区提供高精度的地形信息。在获取方式上，通过中国科学院计算机网络信息中心国际科学数据镜像网站，利用该网站提供的下载工具和数据搜索功能，输入甘河小流域的地理位置范围等参数，成功下载了覆盖甘河小流域的SRTM-DEM数据。这种获取方式方便快捷，且数据来源权威可靠，能够满足本研究对数据准确性和完整性的要求。SRTM-DEM数据具有较高的精度，其垂直精度在平原地区可达±16米，在山区的垂直精度约为±20米。水平精度方面，数据的定位精度能够达到一定的标准，确保了地形信息在平面位置上的准确性。该数据分辨率为30米，即每个栅格单元代表地面上30米×30米的区域，这种分辨率能够较为细致地反映地表的地形起伏变化，对于提取甘河小流域的地形和水文特征具有重要意义。在地形起伏较大的山区，30米的分辨率可以清晰地分辨出山谷、山脊等地形要素；在河谷平原地区，也能够准确地捕捉到地形的微地貌变化，为河网提取和流域边界界定提供了精确的数据基础。SRTM-DEM数据在本研究中具有良好的适用性。首先，其高分辨率和精度能够准确地反映甘河小流域复杂的地形地貌特征，为地形分析提供了可靠的数据支持。通过对该数据的处理和分析，可以精确地提取流域的高程、坡度、坡向等地形参数，这些参数对于理解流域的地形格局和水流运动规律至关重要。其次，该数据的广泛覆盖范围确保了甘河小流域能够完整地被数据覆盖，不存在数据缺失或空白区域，保证了研究的全面性和完整性。再者，SRTM-DEM数据是目前应用较为广泛的地形数据之一，相关的数据处理和分析方法已经相对成熟，在ArcGIS等地理信息系统软件中能够方便地进行处理和分析，这为研究的顺利开展提供了便利条件。三、数据预处理3.1数据格式转换本研究获取的原始SRTM-DEM数据为GeoTIFF格式，为了便于在ArcGIS软件中进行后续的处理和分析，需要将其转换为ArcGIS软件能够直接读取和处理的格式，如EsriGrid格式。在格式转换过程中，选用ArcGIS软件中的“转换工具”模块下的“从栅格转出”工具集中的“栅格转ASCII”工具，将GeoTIFF格式的SRTM-DEM数据转换为ASCII格式。然后，再利用“数据管理工具”模块下的“栅格”工具集中的“创建栅格数据集”工具，将ASCII格式的数据转换为EsriGrid格式。在进行格式转换时，为了确保数据的准确性和完整性，采取了一系列质量控制措施。首先，在转换前仔细检查原始数据的元数据信息，包括坐标系、分辨率、数据范围等，确保这些信息在转换过程中不会丢失或发生错误。其次，在转换过程中，密切关注软件的提示信息和日志文件，及时发现并处理可能出现的错误和警告。最后，在转换完成后，对转换后的数据进行可视化检查，对比转换前后的数据图像，确保数据的空间位置、地形特征等没有发生明显变化。同时，还对转换后的数据进行了统计分析，检查数据的最小值、最大值、均值等统计指标是否合理，以进一步验证数据的质量。通过这些质量控制措施，有效地保证了数据格式转换的准确性和可靠性，为后续的数据分析和处理奠定了良好的基础。3.2数据拼接与裁剪由于甘河小流域的范围较大，获取的SRTM-DEM数据可能由多景数据组成。为了获得完整的覆盖甘河小流域的地形数据，需要对这些多景数据进行拼接处理。在ArcGIS软件中，使用“数据管理工具”模块下的“栅格”工具集中的“镶嵌至新栅格”工具进行数据拼接。在使用该工具时，将所有需要拼接的SRTM-DEM数据添加到“输入栅格”列表中，设置好输出路径和文件名，选择合适的镶嵌运算符，如“最大值”“最小值”“平均值”等。本研究根据数据的特点和研究需求，选择“平均值”作为镶嵌运算符，以确保拼接后的地形数据能够更加准确地反映实际地形情况。同时，还设置了输出数据的像元大小、数据类型等参数，使其与原始数据保持一致，避免因参数设置不当而导致的数据失真或精度损失。完成数据拼接后，得到的是覆盖范围较大的地形数据，其中包含了许多与甘河小流域研究区域无关的数据。为了减少数据量，提高后续处理的效率和精度，需要根据甘河小流域的研究区域范围对拼接后的数据进行裁剪。在ArcGIS软件中，利用“数据管理工具”模块下的“栅格”工具集中的“裁剪”工具进行数据裁剪。首先，准备好甘河小流域的研究区域边界矢量数据，将其加载到ArcGIS软件中。然后，在“裁剪”工具中，将拼接后的SRTM-DEM数据作为输入栅格，将研究区域边界矢量数据作为裁剪要素，设置好输出路径和文件名，点击“确定”按钮即可完成数据裁剪操作。在裁剪过程中，同样需要关注输出数据的质量，检查裁剪后的边界是否准确，是否存在数据丢失或异常值等问题。通过对裁剪后的数据进行可视化检查和统计分析，确保裁剪后的SRTM-DEM数据能够准确地覆盖甘河小流域的研究区域，为后续的地形分析和流域提取提供准确的数据支持。3.3洼地填充与平滑处理在SRTM-DEM数据中，由于数据采集误差、地形测量方法的局限性以及内插算法的影响，常常会存在一些洼地（被较高高程区域围绕的低地）。这些洼地的存在会对水流模拟产生显著的影响，导致流域水流不畅，不能形成完整的流域网络。在计算水流方向时，洼地中的栅格会使得水流方向出现错误，从而影响汇流累积量的计算和河网提取的准确性；在进行水文分析时，洼地会导致水流在这些区域聚集，无法按照实际的水流路径流动，使得模拟的水文过程与实际情况不符。因此，在利用模拟法进行流域地形分析时，需要首先对DEM数据中的洼地进行处理。填充洼地的方法主要有两种：一种是基于局部算法的填充方法，另一种是基于全局算法的填充方法。本研究采用基于局部算法的填充方法，即利用ArcGIS软件中的“水文分析”工具集中的“填洼”工具进行洼地填充。该方法的原理是通过计算每个洼地的最低出水口高程，将洼地内的栅格高程值提高到最低出水口高程，从而消除洼地。在使用“填洼”工具时，设置合适的填充阈值非常关键。如果填充阈值设置过小，可能无法完全填充所有的洼地；如果填充阈值设置过大，可能会过度填充，导致一些真实的地形特征被破坏。本研究通过多次试验，根据甘河小流域的地形特点和数据精度，确定了合适的填充阈值，使得洼地填充既能有效地消除洼地对水流模拟的影响，又能最大程度地保留真实的地形特征。除了洼地填充，还对SRTM-DEM数据进行了平滑处理。平滑处理的目的是消除数据中的噪声和微小的地形起伏，使地形表面更加光滑，从而提高水流模拟的准确性。在ArcGIS软件中，使用“空间分析”工具集中的“邻域分析”工具进行平滑处理。具体操作是选择“均值”作为邻域统计类型，设置合适的邻域大小，如3×3、5×5等。本研究经过对比不同邻域大小的平滑效果，选择了5×5的邻域大小，在这个邻域大小下，能够在有效地消除噪声的同时，较好地保留地形的主要特征。平滑处理的原理是通过计算邻域内栅格的平均值，用该平均值代替中心栅格的值，从而使地形表面更加平滑。经过洼地填充和平滑处理后，得到了更加准确和光滑的SRTM-DEM数据，为后续的河网提取和流域边界提取提供了高质量的数据基础。四、基于SRTM-DEM的地形特征提取4.1高程信息提取与分析利用ArcGIS软件强大的空间分析功能，从经过预处理的SRTM-DEM数据中成功提取出甘河小流域的高程数据。通过对提取的高程数据进行深入分析，绘制出甘河小流域的高程分布图，清晰地展示了该流域的高程分布特征。甘河小流域的高程呈现出明显的西北高、东南低的分布格局。流域西北部的山地地区，高程多在500-1000米之间，其中最高峰海拔达到了[具体山峰名称及海拔高度]，形成了高耸的地形屏障，阻挡了来自北方的冷空气，对流域内的气候和降水分布产生了重要影响。这些山地地区地势陡峭，地形起伏较大，是河流的发源地和水源涵养区，丰富的降水和高山积雪融水为河流提供了充足的水源。而流域东南部则以平原和丘陵为主，高程普遍在200-500米之间，地势较为平坦开阔，有利于农业生产和人类活动的开展。这里土地肥沃，灌溉水源充足，是当地重要的农业产区，主要种植小麦、玉米、大豆等农作物。为了更直观地展示甘河小流域的高程分布情况，制作了高程频率直方图。从直方图中可以看出，高程在300-400米之间的区域面积占比最大，约为[X]%，这表明该高程范围在流域内分布最为广泛，是流域的主要地形区域。在这个高程范围内，地形相对平缓，土壤肥沃，适宜发展农业和建设居民点。而高程大于800米的区域面积占比较小，仅为[X]%，这些高海拔地区主要分布在流域的西北部山区，地形复杂，交通不便，生态环境脆弱，人类活动相对较少，但对于维护区域生态平衡具有重要意义。通过对高程数据的统计分析，得出甘河小流域的平均高程为[具体平均高程数值]米，标准差为[具体标准差数值]。平均高程反映了流域整体的地形高度水平，而标准差则体现了高程数据的离散程度，即地形的起伏变化情况。较大的标准差说明流域内地形起伏较大，存在明显的高差变化，这与实际的地形地貌特征相符。在山区，由于山峰和山谷的交错分布，高程变化剧烈，导致标准差增大；而在平原地区，地形相对平坦，高程变化较小，标准差也相应较小。4.2坡度与坡向计算坡度和坡向是重要的地形参数，对于理解地表的形态、水流运动以及生态系统的分布具有重要意义。在ArcGIS软件中，利用“空间分析”工具箱中的“坡度”和“坡向”工具，对甘河小流域的SRTM-DEM数据进行计算，得到了该流域的坡度和坡向数据。坡度的计算基于地表的高程变化，通过对相邻栅格单元之间的高程差进行分析，确定每个栅格单元的倾斜程度。计算公式为：[坡度计算公式]，其中，[公式中各参数的含义解释]。计算结果以度为单位表示，范围从0°（表示水平地面）到90°（表示垂直地面）。坡向的计算则是确定每个栅格单元的下坡方向，以度为单位表示，范围从0°（表示正北方向）到360°，按照顺时针方向旋转。例如，90°表示正东方向，180°表示正南方向，270°表示正西方向。经过计算，甘河小流域的坡度分布呈现出明显的区域差异。在山地地区，坡度普遍较大，多在15°-35°之间，局部地区可达40°以上。这些陡峭的山坡地形起伏大，水流速度快，容易造成水土流失和山体滑坡等地质灾害。在强降雨天气下，大量的雨水迅速汇聚成地表径流，对山坡上的土壤和植被产生强烈的冲刷作用，导致土壤肥力下降，植被破坏，进而影响生态环境的稳定性。而在平原和丘陵地区，坡度相对较小，一般在5°-15°之间，地形较为平缓，有利于农业生产和基础设施建设。在这些地区，土地平整，灌溉条件良好，适合大规模种植农作物和建设农田水利设施，同时也便于交通线路的铺设和居民点的布局。坡向方面，甘河小流域的坡向分布较为均匀。各个方向的坡向都有一定的面积占比，但总体上南坡和东南坡的面积相对较大。南坡和东南坡由于朝向太阳，光照条件较好，热量充足，有利于植被的生长和光合作用。在这些坡向上，植被种类丰富，生长茂盛，形成了较为稳定的生态系统。同时，南坡和东南坡的气温相对较高，土壤温度也较高，有利于农作物的生长和发育，是农业生产的理想区域。相比之下，北坡和西北坡由于光照时间较短，气温较低，植被生长相对缓慢，生态系统相对脆弱。在冬季，北坡和西北坡容易受到冷空气的侵袭，气温骤降，对植被和农作物造成冻害，影响农业生产和生态环境的稳定。坡度和坡向对甘河小流域的水文和生态过程有着重要的影响。在水文方面，坡度直接影响水流的速度和方向。较大的坡度会使水流速度加快，增加地表径流的流量和流速，容易导致洪水的发生；而较小的坡度则使水流速度减缓，有利于水分的下渗和土壤的涵养水源。坡向也会影响降水的分布和蒸发量。南坡和东南坡由于光照充足，气温较高，蒸发量相对较大；而北坡和西北坡则相反，蒸发量较小。这些因素都会影响流域内的水资源分布和水循环过程。在生态方面，坡度和坡向影响植被的分布和生长。不同的坡度和坡向提供了不同的光照、温度和水分条件，适合不同类型的植被生长。例如，喜光植物多分布在南坡和东南坡，而耐阴植物则多分布在北坡和西北坡。此外，坡度还会影响土壤的侵蚀和肥力，进而影响植被的生长和生态系统的稳定性。4.3地形起伏度与粗糙度分析地形起伏度和粗糙度是衡量地形复杂程度的重要指标，对于研究流域的地貌特征和人类活动具有重要意义。地形起伏度是指一定区域内最高点与最低点的高差，它反映了地形的垂直变化程度；地形粗糙度则是通过对坡度和坡向的综合分析，衡量地形表面的不规则程度。在计算地形起伏度时，使用ArcGIS软件中的“焦点统计”工具，设置合适的分析窗口大小，如3×3、5×5等。本研究选择了5×5的分析窗口，在这个窗口大小下，能够较好地反映地形的局部起伏特征。通过计算每个分析窗口内的高程最大值和最小值，然后用最大值减去最小值，得到该窗口内的地形起伏度。将所有分析窗口的地形起伏度结果进行汇总，得到甘河小流域的地形起伏度分布图。对于地形粗糙度的计算，首先利用前面计算得到的坡度和坡向数据，通过特定的算法将其转换为地形粗糙度指标。具体的计算方法为：[地形粗糙度计算公式]，其中，[公式中各参数的含义解释]。计算结果越大，表示地形表面越粗糙，不规则程度越高；反之，则表示地形表面越平坦，不规则程度越低。经过计算，得到甘河小流域的地形粗糙度分布图。从地形起伏度分布图可以看出，甘河小流域的地形起伏度呈现出明显的区域差异。在山地地区，地形起伏度较大，多在200-500米之间，部分区域甚至超过500米。这些地区山峰林立，山谷幽深，地形复杂多变，是地质构造运动和长期的外力侵蚀作用的结果。强烈的地壳运动使地层发生褶皱和断裂，形成了高山和深谷；而长期的流水侵蚀、风力侵蚀等外力作用则进一步塑造了复杂的地形地貌，使得地形起伏度增大。在这些地区，地形起伏度大导致交通不便，工程建设难度大，人类活动相对较少，但也为旅游业的发展提供了独特的自然景观资源，如奇峰异石、峡谷飞瀑等，吸引了众多游客前来观赏。而在平原和丘陵地区，地形起伏度相对较小，一般在50-200米之间，地形较为平缓，地势相对开阔。这些地区是人类活动的主要区域，交通便利，农业发达，人口密集。由于地形起伏度小，有利于大规模的农业生产和基础设施建设，如修建道路、桥梁、农田水利设施等，促进了区域经济的发展。地形粗糙度分布图显示，甘河小流域的地形粗糙度也存在明显的区域差异。在山地和河流附近，地形粗糙度较高。山地地区由于地形起伏大，坡度和坡向变化复杂，导致地形表面不规则程度高；河流附近则由于水流的侵蚀和堆积作用，形成了复杂的河漫滩、阶地等地形，使得地形粗糙度增大。在这些地区，地形粗糙度高增加了土地利用的难度，限制了一些大规模的开发活动。例如，在山地地区，由于地形崎岖，土地破碎，不利于大规模的机械化农业生产，只能发展一些小规模的特色农业或林业；在河流附近，由于地形复杂，容易受到洪水的威胁，在进行建设活动时需要考虑防洪和地质稳定性等因素。而在平原地区，地形粗糙度相对较低，地形表面较为平坦，土地利用效率高，适合进行大规模的农业生产和城市建设。在平原地区，地势平坦开阔，有利于大规模的农田开垦和灌溉设施建设，提高农业生产效率；同时，也便于城市的规划和布局，建设各类基础设施和建筑物，促进城市化进程的发展。地形起伏度和粗糙度对甘河小流域的地貌和人类活动产生了重要影响。在地貌方面，地形起伏度和粗糙度是塑造流域地貌形态的重要因素。较大的地形起伏度和粗糙度使得地表形态复杂多样，形成了高山、峡谷、丘陵、平原等不同的地貌类型，这些地貌类型相互交错，构成了甘河小流域独特的自然景观。在人类活动方面，地形起伏度和粗糙度影响着土地利用方式、交通建设和工程建设等。在地形起伏度和粗糙度较大的地区，土地利用受到限制，交通建设和工程建设难度大，成本高；而在地形起伏度和粗糙度较小的地区，土地利用相对容易，交通建设和工程建设条件较好，有利于经济的发展和人类活动的开展。因此，在进行流域规划和开发时，需要充分考虑地形起伏度和粗糙度的影响，合理布局各类人类活动，实现区域的可持续发展。五、甘河小流域河网提取5.1流向确定在河网提取过程中，确定水流方向是至关重要的一步。本研究采用D8算法来确定水流方向，该算法基于最陡坡度法，假设单个栅格中的水流只能流入与之相邻的8个栅格中。其原理是在3×3的DEM栅格上，计算中心栅格与各相邻栅格间的距离权落差，即栅格中心点落差除以栅格中心点之间的距离，取距离权落差最大的栅格为中心栅格的流出栅格，从而确定水流方向。在ArcGIS中，通过“空间分析工具”模块下的“水文分析”工具集中的“流向”工具来实现D8算法计算水流方向。具体操作步骤如下：首先，打开“流向”工具对话框，将经过预处理和地形分析后的SRTM-DEM数据作为输入表面栅格；然后，设置好输出表面栅格的保存路径和文件名；在参数设置中，选择默认的D8流向类型，其他参数保持默认设置；最后，点击“确定”按钮，即可得到甘河小流域的水流方向数据。在这个过程中，需要注意输入数据的准确性和完整性，以及参数设置的合理性，以确保计算结果的可靠性。水流方向数据以栅格形式呈现，每个栅格的值代表了该栅格单元的水流方向，取值范围为1-255，其中1表示东，2表示东南，4表示南，8表示西南，16表示西，32表示西北，64表示北，128表示东北，这些数值对应的方向编码为后续的汇流累积量计算和河网提取提供了关键依据。5.2汇流累积量计算汇流累积量的计算是基于前面确定的水流方向数据，其计算原理是模拟地表水流的汇流过程，统计每个栅格的上游集水区贡献的累积水量，通常以像元数量来表示。在ArcGIS中，利用“空间分析工具”模块下的“水文分析”工具集中的“流量”工具进行汇流累积量计算。具体操作时，在“流量”工具对话框中，将前面计算得到的水流方向数据作为输入流向栅格；设置好输出蓄积栅格数据的保存路径和文件名；其他参数采用默认设置，点击“确定”按钮即可完成计算。汇流累积量的大小代表了上游有多少栅格的水流方向最终流过该栅格，它反映了流域内水流的汇聚程度。在甘河小流域中，汇流累积量较大的区域通常是河流的主干道和较大的支流所在位置，这些区域接收了来自上游众多栅格的水流，表明水流在此处大量汇聚；而汇流累积量较小的区域则多为流域的上游源头和地势较高的地区，水流相对分散。通过对汇流累积量数据的分析，可以清晰地了解流域内水流的分布和汇聚特征，为后续的河网提取提供重要的数据支持，准确地识别出河流的潜在位置和走向。计算得到的汇流累积量数据结果如图[具体图号]所示，从图中可以直观地看出汇流累积量在流域内的分布情况，颜色越深表示汇流累积量越大，颜色越浅表示汇流累积量越小，这种可视化的表达方式有助于更直观地理解流域内水流的汇聚模式。5.3河网生成与阈值确定河网的生成是基于汇流累积量数据，通过设定合适的汇流累积阈值，将汇流累积量大于阈值的栅格视为河流，从而提取出河网。在ArcGIS中，使用“地图代数”工具集中的“栅格计算器”工具进行河网生成。在“栅格计算器”中，输入表达式“Con(汇流累积量栅格>阈值，1,NoData)”，其中“汇流累积量栅格”为前面计算得到的汇流累积量数据，“阈值”为设定的汇流累积阈值，“1”表示满足条件（汇流累积量大于阈值）的栅格赋值为1，代表河流，“NoData”表示不满足条件的栅格赋值为空值，即非河流区域。通过调整阈值的大小，可以得到不同密度的河网。为了确定最佳阈值，本研究进行了一系列实验，分别设置汇流累积阈值为100、200、300、400、500等，并对不同阈值下提取的河网进行分析。当阈值设置为100时，提取出的河网非常密集，包含了许多细小的支流和可能的伪河道，这些细小的支流和伪河道可能是由于数据噪声或地形微小起伏导致的，它们在实际的水系中可能并不存在，或者对流域的整体水文特征影响较小。随着阈值逐渐增大到200时，河网密度有所降低，一些细小的支流和伪河道被去除，河网结构更加清晰，与实际水系的主要河道更为接近，能够较好地反映流域的主要水系格局。当阈值增大到300时，河网进一步简化，一些较小的支流也被舍去，虽然河网更加简洁，但可能会丢失一些对流域水文分析有一定价值的信息。当阈值继续增大到400和500时，河网变得过于稀疏，许多重要的支流都被忽略，无法完整地反映流域的水系特征。通过对不同阈值下提取的河网与实际水系的对比分析，综合考虑河网的完整性和准确性，确定在甘河小流域中，汇流累积阈值为200时提取的河网效果最佳，能够在保留主要河道和重要支流的同时，有效地去除噪声和伪河道，准确地反映甘河小流域的实际水系情况。5.4河网精度验证为了评估提取河网的准确性和精度，将提取得到的河网与实际水系图进行对比验证。实际水系图来源于当地的水利部门或相关地理信息数据库，具有较高的准确性和权威性。首先，将提取的河网数据和实际水系图在ArcGIS软件中进行叠加显示，通过直观的视觉对比，初步观察提取河网与实际水系在空间位置、走向和河流分布上的一致性。从叠加图中可以看出，大部分主要河流的位置和走向与实际水系较为吻合，但在一些细节部分，如部分支流的长度和位置，仍存在一定的差异。为了更准确地评估河网提取的精度，采用定量分析的方法，计算提取河网与实际水系的吻合度指标，如河流长度吻合率、河流位置偏差等。河流长度吻合率的计算方法是：将提取河网中每条河流的长度与实际水系中对应河流的长度进行对比，计算两者长度的比值，然后对所有河流的长度比值求平均值，得到河流长度吻合率。河流位置偏差的计算则是通过计算提取河网中河流中心线与实际水系中对应河流中心线之间的平均距离来衡量。经过计算，提取河网与实际水系的河流长度吻合率达到了[具体吻合率数值]，表明在河流长度方面，提取结果与实际情况较为接近；河流位置偏差的平均值为[具体偏差数值]米，说明在河流位置上，提取河网与实际水系存在一定的偏差，但总体偏差在可接受范围内。通过与实际水系图的对比验证，虽然提取的河网在一些细节上存在一定的误差，但总体上能够较为准确地反映甘河小流域的水系特征，满足了本研究对河网提取精度的要求。对于存在的误差，可能是由于SRTM-DEM数据的精度限制、地形复杂导致的水流模拟误差以及阈值选择的不确定性等因素引起的。在后续的研究中，可以进一步优化数据处理方法和算法，提高河网提取的精度，以更好地服务于甘河小流域的水文分析和管理决策。六、甘河小流域边界提取6.1流域边界提取方法本研究基于水流方向和汇流累积量确定流域边界。在完成流向确定和汇流累积量计算后，利用ArcGIS软件中的“水文分析”工具集中的“分水岭”工具进行流域边界提取。其原理是根据水流方向数据，从每个栅格单元出发，按照水流方向追踪，直到遇到分水岭为止，这些分水岭的连线就构成了流域边界。具体操作步骤为：打开“分水岭”工具对话框，将流向数据作为输入流向栅格，设置好输出栅格数据的保存路径和文件名，其他参数采用默认设置，点击“确定”按钮即可得到甘河小流域的边界数据。通过这种方法提取的流域边界，能够准确地反映流域内水流的汇集范围，将所有流向同一出口的区域划分在同一个流域内。该方法基于地形数据和水流模拟，充分考虑了地形对水流的影响，具有较高的科学性和准确性。在甘河小流域中，这种方法能够清晰地界定流域的范围，为后续的水文分析和管理提供了准确的边界数据。例如，在一些山区，通过该方法能够准确地识别出山脊等地形特征作为流域边界，避免了因人为判断不准确而导致的边界误差。同时，这种基于数据和算法的提取方法，具有较高的自动化程度，能够快速地完成流域边界的提取工作，提高了研究效率。6.2子流域划分在确定甘河小流域的边界后，根据河网和地形特征进行子流域划分。利用ArcGIS软件中的“流域分析”工具集，其中的“河网流域划分”功能可以根据水流方向和汇流区域自动识别子流域。在进行子流域划分时，以提取的河网作为基础，将河网中的每条河流作为子流域的划分依据，从河流的源头开始，沿着水流方向，将流向该河流的区域划分为一个子流域。同时，考虑地形特征，如分水岭、山谷等，确保子流域的划分符合地形地貌特征和水文学规律。例如，在地势较高的山区，以分水岭为界划分不同的子流域；在山谷地区，根据水流的汇聚情况确定子流域的范围。子流域划分在流域研究中具有重要意义。首先，不同的子流域可能具有不同的地形、土壤、植被等特征，这些特征会影响到流域内的水文过程、水土流失情况以及生态系统的分布。通过划分不同的子流域，可以更细致地研究这些因素对流域生态环境和水资源的影响，为制定针对性的保护和管理措施提供依据。其次，在水资源管理方面，不同子流域的水资源量、用水需求和开发利用程度可能存在差异。通过子流域划分，可以对每个子流域的水资源进行单独评估和管理，实现水资源的合理分配和高效利用。再者，在防洪减灾方面，了解不同子流域的洪水发生规律和汇流特性，有助于制定更加科学合理的防洪预案，提高防洪减灾能力。例如，对于洪水风险较高的子流域，可以加强防洪工程建设和洪水监测预警；对于水资源短缺的子流域，可以采取节水措施和水资源调配方案。6.3流域边界验证与修正为了确保提取的甘河小流域边界的准确性，利用实地调查和其他数据对其进行验证和修正。组织专业人员对甘河小流域进行实地考察，记录流域边界处的地形特征、河流走向以及土地利用情况等信息。将实地调查得到的信息与提取的流域边界数据进行对比分析，检查边界是否与实际地形相符。同时，收集当地的地形图、水系图以及高分辨率遥感影像等其他数据，将这些数据与提取的流域边界进行叠加分析，从多个角度验证边界的准确性。在验证过程中，发现一些区域的提取边界与实际情况存在偏差。例如，在部分山区，由于地形复杂，数据分辨率有限，导致提取的边界与实际的分水岭位置存在一定的偏移；在一些河流交汇处，由于水流情况复杂，提取的边界未能准确反映实际的汇水区域。针对这些偏差，根据实地调查和其他数据提供的信息，对提取的流域边界进行修正。在ArcGIS软件中，使用“编辑”工具对边界进行手动调整，使其更加符合实际地形和水流情况。通过多次验证和修正，最终得到了更加准确的甘河小流域边界数据，提高了研究结果的可靠性，为后续的流域分析和管理工作提供了坚实的数据基础。七、结果分析与讨论7.1小流域地形与河网特征分析通过对甘河小流域地形与河网特征的提取和分析，发现二者之间存在着紧密的相互关系。从地形对河网的影响来看，地势高低和坡度大小直接决定了水流的方向和速度。在甘河小流域，地势西北高东南低，这使得水流总体上自西北向东南流动，从而塑造了河网的基本走向。在坡度较大的山地地区，水流速度快，侵蚀作用强烈，容易形成深切的河谷和狭窄的河道，这些区域的河网往往较为密集，且河流的弯曲度较小，以适应快速的水流。而在坡度较小的平原和丘陵地区，水流速度减缓，河流的侧向侵蚀作用增强，容易形成宽阔的河谷和弯曲的河道，河网的密度相对较低，但河流的连通性较好。例如，在甘河的上游山区，由于地势陡峭，坡度较大，河流的流速较快，形成了众多短小而湍急的支流，这些支流相互交织，构成了密集的河网；而在下游的平原地区，河流流速减缓，河道逐渐变宽，河网的密度相对降低，但河流之间的连通性更好，形成了较为稳定的水系格局。河网对地形也有着重要的塑造作用。河流的侵蚀、搬运和堆积作用不断改变着地形地貌。在河流的上游，强烈的下切侵蚀作用使得河谷不断加深，形成了峡谷等地形；在中游，河流的侧向侵蚀作用使得河谷逐渐拓宽，形成了河漫滩和阶地等地形；在下游，河流携带的大量泥沙在河口地区堆积，形成了三角洲等地形。甘河小流域的河谷平原就是在河流长期的侵蚀和堆积作用下形成的，这些平原地区地势平坦，土壤肥沃，是人类活动的主要区域。此外，河网的分布还影响着地形的稳定性。河流的存在使得地下水水位升高，增加了土壤的含水量，从而影响了土壤的力学性质和抗侵蚀能力。在一些河流附近，由于地下水位较高，土壤容易发生软化和滑动，导致山体滑坡等地质灾害的发生。而在河网稀疏的地区，由于缺乏水流的作用，地形的变化相对较为缓慢，地貌形态相对较为稳定。7.2与其他研究结果对比将本研究基于SRTM-DEM数据提取甘河小流域的结果与同类研究进行对比分析，发现存在一定的差异。在河网提取方面，与一些利用高分辨率航空影像或LiDAR数据提取河网的研究相比，本研究提取的河网在细节上存在一定差距。高分辨率航空影像和LiDAR数据能够提供更详细的地表信息，对于一些细小的支流和河网的细微特征能够更准确地识别和提取。而SRTM-DEM数据由于分辨率的限制，可能会遗漏一些细小的支流，对于河网的一些细微变化也难以准确捕捉。例如，在某些高分辨率影像研究中，能够清晰地识别出一些长度较短、宽度较窄的小溪流，而在本研究中，这些小溪流可能由于汇流累积量未达到设定阈值而未被提取出来。然而，本研究基于SRTM-DEM数据提取河网的优势在于能够快速、全面地获取大范围的河网信息，对于宏观的水系格局把握较为准确，并且成本相对较低，适用于对流域河网进行初步的分析和研究。在流域边界提取方面，与一些采用实地调查结合传统测量方法确定流域边界的研究相比，本研究利用SRTM-DEM数据提取的流域边界在准确性上存在一定误差。实地调查和传统测量方法能够直接获取地面的实际情况，对于一些地形复杂、难以通过数据准确判断的区域，能够更准确地确定流域边界。而基于SRTM-DEM数据提取流域边界时，可能会受到数据精度、地形复杂性以及算法局限性等因素的影响，导致边界的确定存在一定偏差。在山区，由于地形起伏较大，数据分辨率有限，可能会使提取的流域边界与实际分水岭位置存在一定偏移；在河流交汇处，水流情况复杂，也可能导致提取的边界未能准确反映实际的汇水区域。然而，本研究方法具有高效、快速的特点，能够在短时间内获取整个流域的边界信息，为流域的宏观分析和管理提供了便利。这些差异的产生主要是由于数据来源、研究方法和研究区域的不同所导致的。不同的数据具有不同的精度和分辨率，不同的研究方法也有其各自的优缺点和适用范围，而研究区域的地形地貌、气候条件等因素也会对提取结果产生影响。7.3不确定性分析数据精度是影响提取结果不确定性的重要因素之一。SRTM-DEM数据虽然具有较高的精度，但仍然存在一定的误差。其垂直精度在平原地区可达±16米，在山区约为±20米，水平精度也存在一定的误差范围。这些误差会在地形分析、河网提取和流域边界提取等过程中逐渐积累和传播，导致最终提取结果的不确定性。在地形分析中，高程数据的误差会影响坡度、坡向等地形参数的计算精度，进而影响河网提取和流域边界提取的准确性。如果高程数据存在误差，可能会使计算得到的坡度和坡向与实际情况不