坐标、分辨率、像素类型五花八门几十个tif铺在GIS里连颜色都对不上想导出一张整图根本无从下手。我基本隔三差五就要处理一批这样的高程碎图所以把常用的镶嵌合并方法整理一下Global Mapper、ArcGIS、QGIS三套方案都有覆盖应急处理和批量化生产两种需求后面再聊聊镶嵌完之后高程接边断层线的识别和处理。这内容适合测绘、林业、水利、规划、道路设计以及所有被分幅DEM折磨过的人。1. 先说清楚高程碎图为什么要做镶嵌合并1.1 高程碎图的来源与形态大家手里批量下载的高程数据来源无非这么几个地理空间数据云、USGS EarthExplorer、NASA Earthdata、ASF Data Search或者一些省级地理信息平台。这些平台下出来的DEM大多是按全球标准分幅切的比如SRTM是按经纬度1度一格AW3D30是0.5度一格地理空间数据云提供的ASTER GDEM有时候是1度有时候是你自定义范围内自动切片。下载选项里一旦选了“批量下载”或者“按范围下载”到手的就是一堆几千乘几千像素的tif小文件。这些小文件的格式相当统一但不完全一致。最常见的坑是同一批下载的DEM有的是整型存储有的是浮点型存储有的坐标系是WGS84地理坐标有的平台默认给你加了UTM投影有的中心点经纬度不同像元大小在N37和N45两个纬度带里还不太一样——因为地理坐标系下同是1弧秒在不同纬度代表的实地距离不同。如果直接把这些文件堆进ArcScene或者Global Mapper当普通图层浏览没问题但一做水文分析、坡度计算、等高线提取就全露馅了。1.2 不镶嵌直接用的后果好多第一次接触高程数据的朋友会问几十个小tif我用的时候单独加载不行吗理论上每个tif都能单独出图但实际根本没法用。第一水文分析会有大量“断头”现象。流向分析是以栅格为单位的每个小tif边缘那一圈的像元因为看不到外面的数据软件不知道水该往哪里流结果就是沿着图幅接缝出现一圈人工河道后续算汇流累积量全是假的。第二坡度坡向分析会出现明显的“接缝条带”。即使两个相邻图幅地形是连续的由于数值精度或重采样差异接缝两侧的坡度值会在边缘上下跳动打开坡度图一看整个研究区像被切成了豆腐块。第三提取等高线时相邻图幅的同一根等高线错位断开。这个问题很隐蔽。两边DEM分别生成的等高线在某些陡坡段会错开半个到一个等高距拼起来像地形断层一样。单看每个图幅都正常合起来就不对。第四3D显示完全没法看。Global Mapper里加载几十个小tif做三维地形颜色会自动拉伸每个文件一个色阶整个就是一个大花脸。尤其是跨海、跨山区的项目研究区越大花得越厉害。所以镶嵌合并不是可选项而是必须做的预处理。1.3 镶嵌合并在做什么——栅格数据集的融合逻辑镶嵌合并的学术叫法是Mosaic本质上就是把多个空间相邻、部分区域还可能重叠的栅格数据集在统一的坐标系、统一的像元大小、统一的像素类型之下计算出一个覆盖完整目标区域的新栅格。逻辑上说镶嵌分两种情况一种是多个文件之间完全不重叠只是边界相接那合并就是个“拼接”动作把边缝对齐即可另一种是下载时平台为了覆盖整个矩形范围相邻文件之间会有几行几列的重叠或者在测区边缘不同来源的数据重叠了这时候就要确定重叠区每个像元最终取哪个值。后续第4部分我会专门展开讲重叠区的取值算法这里先记住一个概念镶嵌的难点不在“拼”而在“融”。2. 工具选型哪种方式最适合你的场景2.1 三大主流工具对比处理高程碎图的工具用得最多的其实是这三类Global Mapper、ArcGIS工具箱、QGIS加GDAL命令行。我放在一起对比一下方便你按自己的条件选。工具学习曲线费用擅长场景核心功能关键词Global Mapper平缓半天上手商业授权有试用版中小批量快速合并、三维预览、应急处理Combine/Convert工具feather融合多种重采样方式ArcGISDesktop/Pro较陡需理解工具箱逻辑商业授权较贵生产级批处理、Mosaic Dataset、建模自动化Mosaic to New RasterMosaic DatasetTopo to RasterQGIS GDAL中短命令行需要基础完全免费开源大批量自动化、脚本批量处理、服务器环境gdal_merge.pygdalbuildvrtgdalwarp如果你只用一两次手头也只有几十个文件那Global Mapper确实是最快的。它不挑数据拖进去就能看合并完还能顺手把格式转成你要的GeoTIFF或IMG。如果以后要经常做或者要接单位的生产流程ArcGIS的Mosaic to New Raster更正式能纳入ModelBuilder。至于免费方案QGIS里跑一条gdal_merge.py命令就行连界面都不用开。2.2 我的选择逻辑与实际案例说一个我自己的习惯。拿到一批高程碎图先看数量。少于30景我直接开Global Mapper一边看三维一边检查哪些图幅有问题合并完顺手导出。量再大一点或者要周期性更新我就写个QGIS的批处理脚本用gdal_merge.py一条命令全跑完。如果单位有ArcGIS许可且数据要入库给多个部门用那就建一个Mosaic Dataset把文件夹里所有DEM直接灌进去ArcGIS服务器端可以做动态镶嵌客户端看到的是一整张连续的高程影像但底层不需要物理合并成一个巨大文件这个在生产环境里很好用。有朋友问是不是只要用Mosaic to New Raster就能万事大吉。我觉得不能迷信。Mosaic to New Raster是物理合并一次跑完输出一个文件后续再做裁剪、投影、填挖都方便而Mosaic Dataset是逻辑镶嵌管理灵活但不适合直接扔给其他软件做二次处理。我们这里讲“合并”主推物理合并最后会提虚拟镶嵌做补充。2.3 准备工作的共同点不管选哪套方案合并前请一定做三件事。第一整理文件清单。把下载的tif全部放在同一个文件夹下命名最好有序比如按纬度排序或者按图幅编号。然后用查看器Global Mapper/QGIS/ArcMap都行一个个确认每个文件能正常打开有些平台下到一半断了文件会损坏这种文件在合并时可能整块报错。第二检查坐标系。右键看每个tif的属性确认坐标系字段是否齐全。如果某个文件连坐标系都是未定义的状态或者两个文件来源不一样一个WGS84一个西安80那必须先把它们统一投影到同一个坐标系再合并。绝对不能在坐标系不一致的情况下镶嵌否则会出现数公里的偏移肉眼看着对齐了机器算起来全错。第三统一分辨率。如果这批碎图中既有30米分辨率的SRTM又有12米分辨率的ALOS就不能直接合并因为像元尺寸不同输出网格对不上。必须先通过重采样把低分辨率数据升到和高分辨率一样的像元尺寸或者反过来全部降到统一的较低分辨率。别担心重采样丢精度后面会讲怎么选重采样方法。3. 实操三套方案完整跑一遍3.1 方案AGlobal Mapper 快速镶嵌适合新手和应急Global Mapper是我电脑里最常用的栅格预处理工具尤其是应急处理时特别省事。它的镶嵌功能藏在左侧“Overlay Control Center”叠加控制中心里。第一步把需要合并的tif全部拖进Global Mapper窗口。软件会自动按空间位置把它摆放好如果有重叠区域会有提示画面显示的是当前叠在最上面那一层。第二步在Overlay Control Center里按住Ctrl选中所有需要合并的图层单击右键选择“Combine/Convert”选项。这个命令是全套功能里最高频的它既能合并栅格也能做格式转换、重投影、裁剪。第三步在弹出的设置窗口里在“Options”标签下设置“Output Format”为GeoTIFF或者其他格式如IMG。切记在“Raster Options”里勾选正确的位深——高程数据一般选32-bit floating point如果文件是整数形式且有小数精度需求也必须选浮点型。若数据本身是SRTM等整米级精度选16-bit signed integer也可但普通情况下我建议无脑32位浮点型。第四步关键操作在“Import/Export Options”里设置NoData值。不同来源的高程数据无效值五花八门有些是0有些是-9999有些是32767。如果不统一设置合并后这些无效值会参与计算把地形挖出千米大坑。Global Mapper在Combine时可以在“Output Projection”下方看到“No Data Value”设置一个统一的阈值比如-9999并勾选“Treat as No Data”软件会在输出时把所有等于该值的像元统一标记为NoData。第五步处理重叠区。在“Options”下的“Resampling”中选择合适的重采样方式在“Feathering”区域勾选“Use Feathering”并设置羽化宽度。羽化是Global Mapper处理重叠区不错的功能它让重叠区两侧的值做加权渐变而不是生硬地在一条边界上切换这在多景DEM拼接时能显著削弱接缝。如果你的数据重叠区本来就很小羽化宽度可以设小一点比如10个像元重叠区多且宽就设成30到50个像元。最后点击“OK”软件会弹出一个进度条等它跑完再在Overlay Control Center里右键输出文件“Zoom to Layer”检查有没有异常空洞然后就可以另存或直接使用。整个过程大概几分钟无需写任何代码新手也能独立完成。3.2 方案BArcGIS Mosaic to New Raster适合批量和生产ArcGIS里常用的物理镶嵌工具是“Mosaic to New Raster”位置在ArcToolbox - Data Management Tools - Raster - Mosaic to New Raster。ArcGIS Pro的路径略有不同但参数一样。打开工具后界面参数如下Input Rasters把要合并的所有tif都加进来也可以拖入一个文件夹后逐个添加。Output Location输出目录注意这里只选文件夹不是完整文件名。Raster Dataset Name输出文件名建议加.tif后缀没有后缀时默认存成Esri GRID格式后续给别人传数据会很不方便。Coordinate System可以留空它会用第一个输入文件的坐标系也可以在这里显式选择一个坐标系强行统一输出投影。如果输入文件坐标系不一致强烈建议在这里选目标坐标系。Pixel Type下拉菜单里有一堆选项建议选32_BIT_FLOAT理由前面说过。如果用16_BIT_SIGNED高程大于32767米的地方基本只有极高山才可能会溢出而且在填挖、计算坡度时会有整数运算精度问题。Number of Bands高程是单波段数据填1即可。Mosaic Method这个很关键它决定重叠区域如何取值。常用的有FIRST、LAST、BLEND、MAX、MIN。FIRST和LAST是按输入列表顺序取上表面值适用于不同时相数据需要优先某一块的场合MAX适合取最高的值比如做建筑高度、树冠表面DSM时要保证选中最高的物体MIN则相反。如果只是普通地形DEM推荐BLEND它会在重叠区做加权平均接缝相对自然。Mosaic Colormap Mode如果是高程数据不用管保持默认如果是彩色影像再考虑MATCH等。设置完成后点击运行。输出单文件tif即是合并完成的DEM。ArcGIS的优势在于批处理和自动化。可以把上述工具拖进ModelBuilder然后使用“迭代栅格”控件遍历整个文件夹甚至用Python脚本循环处理每个月更新的多批次数据。生产单位如果每天都有新到的图幅可以直接把它做成一个Geoprocessing tool长期使用。3.3 方案CQGIS GDAL 命令行免费可自动化如果不想花钱买商业软件或者要在大批量文件上重复执行GDAL命令行是神器。QGIS安装目录下自带OSGeo4W Shell打开后就能直接用gdal_merge.py。最朴素的合并命令如下gdal_merge.py -o merged_dem.tif -ot Float32 -a_nodata -9999 -co COMPRESSLZW -co TILEDYES input1.tif input2.tif input3.tif这条命令把input1到input3合并成一个merged_dem.tif位深设为32位浮点NoData设置为-9999使用LZW无损压缩并开启内部瓦片化。如果文件特别多可以用通配符gdal_merge.py -o merged_dem.tif -ot Float32 -a_nodata -9999 -co COMPRESSLZW ./dem_tiles/*.tifGDAL还有一个更聪明的工具gdalbuildvrt。它不实际拷贝像素数据而是生成一个虚拟栅格VRT文件记录所有tif的空间位置和文件路径。从软件层面看VRT就是一个完整的DEM文件可以随时用于分析每次读取时它动态从源文件中取数据。这个方案几乎瞬间完成且不占用额外磁盘空间适合需要反复调整范围或临时预览的场景。gdalbuildvrt merged.vrt ./dem_tiles/*.tif要注意gdalbuildvrt默认不处理数据重叠问题——多个文件重叠时它会选第一个文件的值。如果重叠区很小且数据质量一致这没问题但重叠区大且值有差异最好先统一规范后再用gdalbuildvrt做虚拟镶嵌或者直接用gdal_merge.py做物理合并。另外如果在合并前需要统一投影先用gdalwarp把每个tif转到目标坐标系再合并gdalwarp -t_srs EPSG:32650 -r bilinear input.tif output_32650.tif这里的EPSG:32650是UTM 50N的EPSG代码按你的测区纬度替换即可。重采样方法选bilinear双线性因为高程是连续表面双线性和三次卷积都可用但三次卷积在陡地形上可能出现超调值超出邻域实际范围。我呢平时最常用的组合就是先gdalinfo批量检查所有文件的坐标系和像素类型发现不一致的先处理然后一条gdal_merge.py搞定。再配合Python脚本几秒钟能处理上百个文件。4. 关键细节重叠区处理、NoData 与精度问题4.1 重叠区到底取哪个值——融合权重原理很多人以为镶嵌就是把两个tif像拼图一样“啪”地合到一起接缝两边一贴就行。真实情况是相邻DEM之间通常有重叠区两个文件在重叠区对同一地物的高程值往往不一致差几米到几十米都有可能。为什么同一个地形两片数据测出来不一样因为它们的原始数据源可能不同、处理算法不同、时相不同。比如一个片区用的是SRTM 1弧秒另一个片区用的是ASTER GDEM两者在山区系统误差就能差15米以上。如果镶嵌时只是硬拼接重叠区就会出现“重影”或“台阶”。这时候就要选重叠区的取值策略。FIRST / LAST很好理解按文件在输入列表里的顺序取最上面那层的值。这种策略适合有明显优先级的数据比如最新的航测数据覆盖了旧数据的一角你就可以把新数据放最前面选择FIRST让新数据覆盖旧数据。MAX策略适合DSM场景比如要保留最高地物重叠区取多个来源的较大值可以有效保留建筑物、树冠顶部。MIN则适合只关心最低地形面的场景比如提取山谷线时。但普通地形DEM这两个策略都会在重叠区留下比较生硬的痕迹。BLEND是最推荐的通用策略。它把重叠区按距离接缝计算权重接缝一侧取一个源文件的权重高越往另一侧权重越低中间渐变过渡。实现方式本质上是线性加权平均新值 权重A × 值A 权重B × 值B其中权重A 权重B 1且权重随重叠区内像元到各自图幅边缘的距离线性变化。这解决了硬拼接的台阶感但也有代价重叠区的细节被平均了如果两块数据系统性偏差太大BLEND会得到一条过渡带高精度应用时仍然不可用。此时应该做的是在镶嵌前先把两个数据源做高程基准校正比如用“shift”工具把其中一景整体抬高或降低而不是靠融合策略硬抹。4.2 NoData 黑边是镶嵌第一大坑高程碎图的边缘经常有一圈无效值在SRTM的海洋区域、ASTER GDEM的云遮挡区域、或者数据裁切的边界上尤其常见。如果没有统一处理镶嵌后输出文件里会参杂多种NoData值加载到GIS里就显示成黑块或者颜色怪异的大洞。处理原则合并之前确定一个NoData值把这个值在输入文件中全部统一。Global Mapper的Combine设置里有“No Data Value”选项ArcGIS的Mosaic to New Raster会自动把输入文件的NoData转为输出NoDataGDAL则用-a_nodata参数强制指定。做完镶嵌后一定要检查输出文件是否还有黑边。QGIS里加载tif右键属性——透明度将需要设为透明的值设为“NoData值”ArcGIS里在符号系统里把NoData显示为无色Global Mapper则直接在图层属性里取消背景显示。如果发现仍有异常黑色区域大概率是原始文件中某处出现了0值而其他文件的NoData是-99990被当成了真实高程参与计算。处理办法是合并前逐个文件用Raster Calculator把0值重分类为-9999再统一NoData后重新合并。4.3 高程精度与位深选择高程数据的位深选择直接关系到输出精度。很多教程里随手选了8_BIT_UNSIGNED或者16_BIT_SIGNED结果高程值出现明显的阶梯化。原理很简单8位无符号整型只能存0到255连我国东部平原的高程都存不下16位有符号整型能存-32768到32767看起来够了但它是整数高程超过32767米会溢出而且任何小数位都被截断。假如原始数据高程是123.45米存成16位整型就变成了123山区斜坡上0.5米的精度差异全没了。所以高程数据一律用32_BIT_FLOAT32位浮点它能表达很大的数值范围和足够多的小数位。虽然输出文件会变大但高程精度是后续分析的生命线不值得为省磁盘空间牺牲它。同时高程数据在压缩时选LZW它属于无损压缩对浮点高程的压缩率虽然不如对影像的JPEG压缩那么夸张但通常能把文件减小到原来的60%到70%而且没有任何信息损失。在GDAL里使用-co COMPRESSLZW在Global Mapper里输出GeoTIFF时选择压缩类型为LZW在ArcGIS环境选项里设置合适的压缩参数即可。5. 进阶识别被断层线隔开的等高线的高程问题5.1 镶嵌后的“隐形杀手”接缝断层线说实话单纯把栅格合并完不等于大功告成。更隐蔽的问题是——镶嵌后的DEM沿某些线会出现同一条等高线被隔开、两侧无法对应的情况这种线我习惯叫它“断层线”。断层线产生的原因通常有三个。第一相邻两景DEM的高程基准不一致或有系统误差导致等高线在两景交界处上下错动十几米看起来像地形被横向切断。第二在陡崖、峡谷等地形突变带上不同源数据的自动匹配算法自身会产生跳跃一条山谷在接缝两侧被“切”开了。第三原始数据在采集时存在不同期次比如一个片区是雪季影像一个是非雪季影像积雪表面差异会导致高程连续性问题。识别这种断层线靠肉眼盯着渲染图看也是一种办法但不系统容易漏。更可靠的是用坡度分析和曲率分析。5.2 三步定位断层线第一步整幅DEM计算Slope坡度栅格计算器或ArcGIS的Surface工具都行。断层线在坡度图上显示为一条非常突兀的线性高值条带因为两侧地形如果错开局部坡度会瞬间变得很大而它并不对应真实地形。如果坡度图上有线性亮线先圈出来。第二步提取等高线并加载到ArcGIS里对比断层线两侧的等值线走势。如果某条等高线在这条线附近出现了明显位移或者一根完整的等高线在这里断开为两根不同高程的线基本可以确认断层线存在。随后的关键问题是断裂线两侧的高程值到底哪个对、错的一方差多少以及要不要修正。第三步针对每一条可疑断层线用“高程剖面”工具在断层线垂直方向上拉一条剖面线看高程值是突变还是渐变。如果是突变说明两侧数据有系统偏差此时可以在断裂线两侧取一系列样本点分别记录两边的高程值计算差值均值然后用该均值对其中一侧做整体校正。如果研究对象是等高线这个操作可以保证断裂线两侧同一根等高线的高程值一致恢复正确的等值线序列。5.3 实操修复方案修复断层线有三种常见做法按复杂度从低到高排列。方法一镶嵌时开启羽化Feathering这是最省事的办法。在Global Mapper里勾选Feathering选项后重叠区的值会被线性混合断层现象在多数情况下会明显减轻。但注意如果两侧系统误差有几十米羽化出来的过渡带是假地形依然不能用于精度较高的分析。方法二用ArcGIS的Topo to Raster工具重新构建高程表面。Topo to Raster支持输入等高线、河流、湖泊以及“断线”breakline数据它能把已知的等高线、水流线、断层约束等综合起来生成新的DEM。做法是先从原DEM提取等高线再把你识别出的断层线作为breakline要素加载进来并设置相应的高程约束然后生成一张更平滑、断线处合理过渡的新DEM。这个方案适合做区域性的地形修复不太适合严格保留原始采样值。方法三针对断层带局部用“低通滤波”或“中值滤波”处理后重新赋高程。ArcGIS的Focal Statistics工具里选Median以3x3或5x5窗口处理断层带能有效消除条带状噪声。但中值滤波同样会模糊真实细节所以只建议在断层带范围内处理用一个掩膜把断层带圈出来只对圈内像元滤波圈外不动。我个人实际使用的经验是断层修复要以“不破坏真实地形”为前提。我曾经为了消除一条接缝把一整个山区的高程大幅平滑结果后期算坡度、流向时整个流域边界都变了。正确的思路是能识别出断层线再针对断层线周围有限带宽做局部处理而且处理完一定要叠加原始等高线检查是否恢复连续。关于“识别被断层线隔开的等高线的高程”这个热搜词我的理解是在实际地形成图中等高线本质是DEM等值面与水平面的交线。如果DEM内部存在画不出来的断裂线提取等高线时必然会把连续的一条线切成两段其中一段在断裂线左侧显示为一个高程值另一段在断裂线右侧显示为另一个值。识别这类被隔开的等高线高程就是在修复DEM之前先通过裂缝两侧的等值线数值差去推断出该处应该统一为多少米。这个思路放到操作层面就是第二节和第三节提到的绘制剖面、对比两侧值、整体校准。要真正解决断层线导致的等高线错位根源仍在于高程数据本身的一致性。所以我一般在镶嵌结束后做一次“自查”生成等高线叠加在3D地形上如果发现沿某条线等高线发生密集抖动或断开就立即回到原始文件查看是不是混入了坐标系或基准不同的图幅。6. 常见问题与排查技巧实录把我在实操中遇到的典型问题整理成一个速查表按“现象—原因—解决方案”来列遇到问题可以直接对着查。现象可能原因解决方案镶嵌后沿接缝出现明显棱线或台阶两景数据高程基准/系统误差不一致或重叠区用了FIRST/LAST改成BLEND先整体校准其中一景的高程偏移量再合并输出文件边缘有大片黑色区域NoData值未统一部分文件用0标识无效值重新设置NoData为统一值把0重分类为NoData后重跑等高线沿一条线断裂错位断层线未处理原始数据不同期或不同源按第5章的断层定位流程识别断层线做局部中值滤波或Topo to Raster重建镶嵌后文件巨大加载非常卡输出位深过高且未压缩或数据量大使用LZW压缩、TILED参数对Win平台可设置Overview金字塔合并后整体位置偏移几百米以上输入数据坐标系不一致或某个tif缺投影文件先统一投影必要时用QGIS检查每个文件的空间范围边界部分区域显示为平直的0值地形而非真实山谷输入文件中0被当真实高程参与BLEND平均先用条件函数把0替换为NoData重新镶嵌整幅DEM渲染出现密集条纹像拉链重采样方法设置不当或输入文件分辨率不一致统一分辨率后用bilinear重采样如果分辨率差距过大先升尺度再合并使用Global Mapper输出后坐标体系识别不了输出格式选错了或者GeoTIFF里没有写入坐标系信息输出时指定Output Projection并在Export Options里勾选Embed Projection Info几个我踩过并印象深刻的坑具体说下。第一次处理山东省一个县的多个图幅。下载时看着都是SRTM坐标系也都是WGS84于是直接扔进Global Mapper合并。结果出图后发现中间有一条斜向的台阶来回排查了很久最后发现其中两景的高程值整体偏低约8米因为它下载的是另一个处理版本的SRTMSRTMGL3标称高程参考面与主流的SRTMGL1不同。解决办法是把这两景整体加8米后重新合并。这类系统性偏移用BLEND也救不回来只能先做数据级校准再镶嵌。另一次给某单位跑一个流域范围的DEM合并。我用Mosaic to New RasterMosaic Method选了MAX因为想着MAX可以让山谷里可能存在的低值伪影被覆盖掉。结果合并后山区河谷的高程整体偏高了几米——因为MAX把不同图幅在河谷处较高的一侧全取了导致下游坡降失真。后来改成BLEND问题立刻改善。从此我对MAX策略非常谨慎只在DSM类似的表面数据的场景中推荐。还有一次合并完检查成果时发现一个圆形的白色空洞放大后恰好是ASTER GDEM常见的云遮挡空洞。处理思路是先用邻域插值或者是Nibble工具把空洞补上再加入镶嵌流程。如果先镶嵌再补镶嵌时NoData传播到整圈接缝会影响周围像元。顺序很重要先修补单个文件再做镶嵌。最后关于流程顺序我强烈建议不要在原文件上直接操作。下载的高程碎图永远保留一份原始文件合并输出放在另一个目录。原因很简单你要反复试验Mosaic Method、NoData设置、重采样方式每试一次都可能产出不同的结果原始文件一旦被处理工具覆盖修改再想回到起点就得重新下载。7. 最后再分享一点小经验上面这些坑我基本都完整走了一遍才摸清套路。现在拿到一批高程碎图我的常规流程是先gdalinfo批量检查坐标系、像素类型、分辨率有不一致的先统一再快速在Global Mapper里浏览一遍找出有异常的图幅然后按数据量选择GDAL命令或ArcGIS工具合并NoData统一设-9999位深用32位浮点压缩开LZW最后生成等高线和坡度图做双重校验看有没有接缝断层。整个过程半小时左右几十个图幅就能变成一张干净连续的DEM。镶嵌合并这件事难度不在工具操作本身而在于你知不知道数据里藏了什么问题。高程数据不像普通遥感影像颜色错了肉眼能看出来高程错了可能直到算完坡度、提完等高线才发现满屏错误到那时返工成本就高了。所以把“合并前检查、合并后校验”变成习惯能省掉后期大量时间。这内容如果对你有用可以顺着这条线继续往下做合并完的DEM接下来可以做流域提取、坡度分级、可视域分析也可以跟遥感影像叠加做三维地形浏览。工具选Global Mapper还是GDAL都无所谓关键是每一步想清楚为什么要这么做。