30米DEM数据处理实战：SRTM选型、填洼与流域面积计算

简介广东省三十米高清DEM高程数据是一份覆盖广东全域的数字高程模型数据包面向地理信息系统、测绘与城乡规划领域的开发者和研究人员适用于地形分析、坡度坡向计算、汇水区划分及地质灾害风险评估等场景。包内共含二十六个文件以shp和tif为核心载体shp存储高程点、等高线等矢量要素tif提供连续高程栅格影像便于制作三维地形和渲染地表形态同时配套dbf属性表、ovr金字塔、xml元数据、adf辅助文件可在ArcGIS、QGIS中直接读取、显示与缓存加速整体约五十四兆结构清晰。目前已有超过一千六百人浏览学习。数据能细致反映广东省真实地表起伏除绘制高程图与等高线外还可用于制作三维地形模型、提取山谷山脊线、模拟径流路径并结合坡度坡向分析支撑土地规划、道路选线与生态环境评估是区域地理信息项目的基础资料。1. 30米DEM数据源选型为什么我建议优先考虑SRTM而不是ASTER GDEM拿到一个“广东省30米高清DEM高程数据”的需求第一步不是急着下载而是想清楚你要的到底是哪个30米数据源。这个听起来像废话但实际上很多人一开始没区分清楚下载回来的数据用得别扭后面返工成本很高。目前在公开渠道能稳定拿到、且覆盖广东省全境的30米分辨率DEM主流就是两个SRTMShuttle Radar Topography Mission航天飞机雷达地形测绘任务2000年采集目前公开版本有SRTM390米和经过后处理的SRTM 1弧秒约30米版本。国内很多平台标注的“30米DEM”其实就是SRTM1。ASTER GDEMAdvanced Spaceborne Thermal Emission and Reflection RadiometerTerra卫星影像立体像对生成的全球DEM标称分辨率也是30米目前是V3版本。这两个数据源我都在广东省范围内实测对比过。整体感受是SRTM的地形连续性和水文一致性更好ASTER在局部细节上有时候更“碎”。原因在于两者生成原理不同——SRTM是雷达干涉测量不受云层影响对地表连续地形还原好ASTER是光学立体像对受云、阴影、植被覆盖影响大城区和山区容易出现异常塌陷或凸起。对比项SRTM 30米1弧秒ASTER GDEM V3 30米采集方式雷达干涉光学立体像对覆盖范围全球N60°~S56°全球N83°~S83°广东省可用性全覆盖全覆盖垂直精度约5~10米CE90约7~14米常见异常水域高程偏低不稳定云影区空洞、局部尖峰推荐场景水文分析、流域提取、地形建模宏观地形浏览、视觉制图基于这个对比我的建议是做水文分析、流域面积计算、等高线生成这类对地形连续性和精度要求高的任务直接用SRTM 30米如果是做宏观可视化展示ASTER GDEM V3也能用。下文所有实操都以SRTM 30米为例来讲解这也是目前国内“30米高清DEM”语境下最主流的选择。提示国内许多GIS数据平台下载的“30米DEM”需要仔细看元数据部分平台标注的30米是经过插值重采样得到的原始分辨率可能只有90米。务必下载前确认数据来源和原始分辨率信息。2. 广东省30米DEM数据获取实操2.1 下载前的三个准备工作下载广东省30米DEM之前先把三件事准备好否则下载到一半容易卡壳第一个是明确边界范围。“广东省”三个字听上去简单但你要的是全省范围拼接数据还是某个具体地市这决定了你要下载多少个分幅文件。广东省经纬度范围大约在北纬20°13′~25°31′、东经109°39′~117°19′之间涉及SRTM标准分幅大概是十几到二十个分片。如果只做广州、深圳这种单个城市下载对应的分片即可不用全量拉全省。第二个是确认坐标系统。下载回来的SRTM原始数据是WGS84坐标系EPSG:4326平面单位是度。做面积、长度量算时必须投影到合适的平面坐标系广东地区常用的是CGCS2000 3度分带高斯-克吕格投影中央经线一般为111°E第37带或114°E第38带也有用UTM 49N/50N的情况。这个后面处理时会重点说。第三个是准备好数据存放目录。建议按“原始数据→预处理→成果数据”三层目录组织文件命名带来源和分辨率后缀比如GD_SRTM30_orig.tif、GD_SRTM30_fill.tif。这个习惯能帮你省掉后期大量重复劳动。2.2 三种主流下载渠道对比目前国内下载SRTM 30米数据我实测下来有两条免费线路和一条付费备选渠道一OpenTopographyopentopography.org。这是我首推的渠道支持直接框选范围下载可以选SRTM GL130米产品输出GeoTIFF格式还自带一些简单的可视化预览功能。操作门槛最低但国内访问速度不稳定偶尔需要重试。渠道二USGS EarthExplorerearthexplorer.usgs.gov。这是数据源头SRTM数据的原始发布平台搜索功能强大可以按条带号检索也支持矢量边界导入。缺点是界面相对复杂新手第一次操作容易找不到下载按钮。下载需要注册账号免费。渠道三国内地理空间数据云gscloud.cn。对国内用户最友好速度快注册即可下载但需要注意部分数据集可能不是原始SRTM数据而是经过处理的版本下载前要看清楚数据集说明。另外目前平台对批量下载有一定限制。渠道访问速度操作难度是否免费推荐指数OpenTopography一般低是★★★★★USGS EarthExplorer一般中是★★★★地理空间数据云快低是★★★★2.3 OpenTopography下载SRTM 30米数据详细步骤以OpenTopography为例我实际操作的完整流程如下访问opentopography.org首页点击“Find Data”进入数据检索界面。在数据源列表中找到“SRTM GL1”即1弧秒、约30米分辨率有的版本标为“SRTM 1 Arc-Second Global”。使用地图框选工具用矩形选择广东省范围或者点击“Define Polygon”手动导入广东省的边界矢量文件如shp或GeoJSON格式。用矢量边界的好处是下载范围更精确不会多下载不需要的周边区域。设置输出参数。输出格式选GeoTIFF坐标系默认WGS84经纬度像元大小保持原始30米约0.00027778度即可不要勾选重采样功能。填入接收邮箱提交任务。系统会排队处理快则几分钟慢则半小时完成后会给邮箱发下载链接。下载后解压得到整个广东范围的完整GeoTIFF文件命名类似SRTM_GL1_广东范围.tif。注意OpenTopography的任务排队时间不稳定如果急需数据建议一次性把范围框准确不要提交多个重叠任务浪费等待时间。3. 数据预处理填洼是30米DEM最关键的“坑”3.1 为什么要填洼——水文分析的前置条件从平台下载下来的SRTM 30米DEM表面看是一张挺漂亮的高程灰度图但如果你直接拿它去做河流提取、流域面积计算结果大概率惨不忍睹。原因就两个字洼地sink。洼地是DEM中高程值低于周围所有像元的局部低点。它的产生机制有几种SRTM雷达信号在陡峭峡谷、密集林区会产生噪声地面真实情况如喀斯特地貌的天坑、数据处理过程中的插值误差。广东省粤北、粤西地区山地丘陵多植被覆盖率高DEM里的洼地数量相当可观。河流是沿着高差向下流的如果把DEM里的洼地当成自然洼地水流流到洼地就停住了无法继续往下游汇流这样计算出来的河网会断断续续流域面积严重偏差。所以水文分析前必须填洼。3.2 ArcGIS中的填洼操作与参数ArcGIS里操作很简单但有几个参数细节值得注意使用Spatial Analyst工具箱中的**“填洼Fill”工具**输入30米DEM输出参数Z限制Z limit一般填默认值空即无限制填洼。但我个人建议根据研究区地形特点设置一个合理的Z limit值。比如广东省珠三角区域地势平坦填洼量不大Z limit可以设到5~10米而粤北山区地形起伏大Z limit可以设到20~30米防止过度填洼把真实的地形凹陷也抹平了。填洼完成后做一步检查用**“减去Minus”工具**用填洼后的DEM减去原始DEM得到洼地深度栅格。统计这层栅格的最大值和平均值如果最大值超过50米就要警惕是不是填过头了回头检查原始数据是否在特定区域存在明显异常值。3.3 投影转换30米数据必须处理的坐标问题填洼之前的另一个前置步骤是投影转换。OpenTopography下载的数据是WGS84经纬度坐标系而填洼、流向、流量计算这些水文分析在ArcGIS的默认环境里对单位不敏感但当你要计算面积、长度、坡度时就麻烦了——经纬度下的栅格像元不是等面积的直接算出来的面积、坡度是错误的。我推荐的流程是先将原始DEM投影为WGS84 / UTM 49N如果数据范围主要在广东中西部或WGS84 / UTM 50N如果主要是东部或者直接用CGCS2000 / 3-degree Gauss-Kruger zone 38中央经线114°E。重采样方法选双线性插值Bilinear或三次卷积插值Cubic不要用最邻近法因为会导致等高线呈锯齿状。像元大小设置成30米X和Y都填30确保投影前后分辨率一致。投影完成后再填洼这样洼地深度、流向分析的结果才是在等面积平面上计算出来的。3.4 镶嵌与裁剪多分片数据的无缝处理如果你的下载方式是按SRTM分幅下载的那就涉及**镶嵌Mosaic和裁剪Clip**环节。ArcGIS的镶嵌工具Mosaic To New Raster是首选关键选项是波段数填1像素类型选16位有符号整型16 Bit Signed高程有正有负有符号类型才能正确表示负海拔比如珠江口附近低海拔区域。Mosaic运算符选FIRST或LAST都可以但建议勾选“接边羽化”选项来平滑分幅之间的接边痕迹避免后续做坡度、山体阴影分析时出现条带效应。裁剪时建议用广东省行政区划矢量边界裁剪并勾选“使用输入要素裁剪几何”选项让边缘贴合边界线必要时在外扩一个像元避免边缘黑边。4. 30米DEM的实际应用流域面积计算与Python辅助操作4.1 ArcGIS水文分析全流程从填洼到流域面积计算拿到处理干净的30米DEM之后来演示一个高频需求——计算河流流域面积。这对应了热词里那个“arcgis用dem计算河流流域面积”也是很多水土保持、水利规划、环境影响评价项目里最常见的任务。标准的水文分析流程一共五步填洼Fill上面说过了这是第一步。流向计算Flow Direction用D8算法计算每个像元流向哪个相邻像元8个方向。这个步骤的关键是输入必须是填洼后的DEM否则流向会断。流量累积Flow Accumulation统计每个像元有多少上游像元汇入。这一步输出的值就是该点的累积汇水面积像元数乘以像元面积30×30900平方米就是实际汇水面积。河网提取Stream Definition设定累积流量阈值超过阈值的像元即为河流。比如广东省30米DEM累积流量阈值一般设置在1000~5000之间。阈值越低河网越密集。没有绝对标准可以根据想提取的河流等级来调节。我一般会试几个阈值对应生成的河网长度密度对比选一个跟地形图上的水系套合度最高的。分水岭提取Watershed选定某出水口或用倾泻点 Pour Point用分水岭工具框出对应流域范围该工具的输入是流向栅格和倾泻点位置。输出的矢量面就是流域边界面积直接在属性表里就能查到。这五步走完流域面积从DEM上就计算出来了。注意面积结果依赖于投影坐标系地理坐标系下算出来的面积不可信。4.2 Python处理DEM常用技巧裁剪、重采样、DXF叠加热词搜索里出现了“dem dxf 叠加 python”这在实操中是个高频场景——比如你要把设计红线DXF叠到DEM生成的等高线或者地形断面上辅助判断设计对地形的影响。Python处理DEM我推荐GDAL库这是最底层的栅格读写库其他库大多是它的封装。简单演示几个常用操作用Python裁剪DEMfrom osgeo import gdal # 打开DEM ds gdal.Open(GD_SRTM30_fill.tif) # 按shp裁剪 DEM ds_result gdal.Warp(GD_SRTM30_clip.tif, ds, cutlineDSNameguangdong.shp, cropToCutlineTrue, dstSRSEPSG:32649) ds_result None用Python重采样DEMfrom osgeo import gdal # 将30米重采样为90米 ds gdal.Open(GD_SRTM30_clip.tif) ds_result gdal.Warp(GD_SRTM90_resample.tif, ds, xRes90, yRes90, resampleAlgbilinear) ds_result NoneDEM叠加DXF的原理是DXF作为矢量线要素叠加到DEM上通常是查看等高线与设计线的空间关系。实操上推荐用QGIS或ArcGIS先把DXF导入再叠加DEM生成的山体阴影Hillshade和等高线。如果真的需要程序化处理可以用geopandas读取DXF的线要素然后用rasterstats库提取每个线节点的高程值。比如可以用以下代码批量提取折点处的高程import geopandas as gpd from rasterstats import point_query # 读取DXF中的线要素 lines gpd.read_file(design_line.dxf) # 提取线的所有折点 points [] for geom in lines.geometry: for x, y in geom.coords: points.append({geometry: Point(x, y)}) points_gdf gpd.GeoDataFrame(points, crsEPSG:32649) # 提取对应DEM高程值 elevations point_query(points_gdf, GD_SRTM30_clip.tif) points_gdf[elevation] elevations这样每条设计线的每个折点就有高程属性了可以进一步算坡度、填挖方量或者生成纵断面图。这是我在实际项目中反复用到的脚本稳定可靠。4.3 DEM生成等高线的经验参数用30米DEM生成等高线等高距设置是需要动脑子的。原则是等高距不能小于DEM垂直精度的3倍。SRTM 30米数据垂直精度约5~10米所以等高距最好不低于10米20米是广东省山区制图安全选择。如果你非要5米等高距一是等高线会很碎很毛糙二是精度上根本支撑不了属于伪精度。ArcGIS里用“等值线Contour”工具输入填洼后的DEM等高距填20输出的是矢量线。勾选“基于Z值创建等值线”适合生成标注不勾选则只是几何线需要后期加属性标注。生成的等高线在高程突变处可能出现闭合短小圈这是DEM栅格化产生的正常现象可以在ArcGIS里用“消除Eliminate”工具处理。5. 常见问题与避坑技巧实录5.1 下载的数据有明显条带/空洞怎么办SRTM数据在广东省偶尔会出现局部数据空洞尤其是陡峭山地和浓密林区表现为DEM上一块明显的“洞”或异常偏低区域。处理这种问题我实测最有效的方式是下载ASTER GDEM V3的同范围数据做补洞参考。在ArcGIS里用“栅格计算器Raster Calculator”结合条件函数CON做替换或者用GDAL的fill_nodata算法。补洞后的区域边界要做一次3×3窗口的低通滤波均值滤波过渡平滑避免补丁边界生硬。5.2 填洼后河流流向异常、河道弯曲不对这种情况第一怀疑对象是填洼前的流向计算错误。检查一下流程顺序必须是先投影、再填洼、再流向计算。如果先填洼后投影或者先流向后填洼都会导致流向信息错位。第二怀疑对象是DEM噪声导致的平行河道问题parallel flow。解决办法是填洼时设置合理的Z limit然后用ArcGIS的“水文分析——平滑河网”工具对提取的河网做一次平滑处理。记住一点河网提取是分析结果不需要追求特别多的细节关键是主导水系要对。5.3 DEM与实测高程校验偏差过大广东省RTK实测点与SRTM 30米对比平原地区一般偏差在2~5米山区偏差可能到10米以上。如果校验发现系统性偏高或偏低比如所有点位都偏高3米可以用高程纠正方法做整体平移纠正后DEM 原始DEM - 系统偏差值在ArcGIS栅格计算器里直接写GD_SRTM30_fill.tif - 3就完成了一次系统性高程校正。这对于需要和实测数据联合使用的项目来说是非常实用的一招。5.4 批量下载时如何避免遗漏分幅如果按分幅下载我的习惯是先做一张shp网格把整个研究区按SRTM条带号划分出来然后与广东省矢量边界做一遍叠加查询确定需要哪些编号的图幅逐一记录到Excel表格中下载时对照表格核对。等全部下载完成再在ArcGIS里用“Mosaic To New Raster”镶嵌检查四角是否还有空洞有的话回头看是否漏了一幅。6. 30米DEM能否满足项目需求聊聊精度边界最后说说这个30米DEM的适用边界。做省级宏观分析、流域规划、区域地形分类30米完全够用。但具体到单条排水沟、小地块土方量精细化计算、农田平整设计30米就明显不够了这时候要考虑12.5米ALOS PALSAR数据或者1米/2米的高精度航测数据。另外注意“高清”这个词在DEM领域30米不算高清准确说是中等分辨率。所谓“30米高清”更多是相对以前90米SRTM和1公里GDEM数据而言的。如果预算允许、研究区范围不大可以优先选择12.5米ALOS PALSAR DEM处理流程和30米DEM完全一致但地形细节和坡度精度会好不少这在热词里也出现了“alos 12.5米dem数据”说明这确实是一个热度很高的替代选择。ALOS数据在OpenTopography和ASF Data Search上都可以下载。我个人的结论是30米SRTM是精度和覆盖度之间最均衡的免费选择适合做广东省全境这类大范围的DEM底图到了局部精细研究阶段再叠加更高分辨率数据做局部加密也不迟。数据这东西够用就好不必盲目追求“更高”关键是清楚每一级分辨率能支撑什么样的分析结论。这套操作下来从数据获取到预处理再到实际应用基本能覆盖绝大多数用DEM做地形分析的需求。如果一个项目组从一开始就用统一的数据源和处理流程后续所有分析结果的对比性和可复现性都会有保障这也算是我做了这么多DEM数据处理项目后最想提醒大家的一点。本文还有配套的精品资源点击获取