如何用gdal_warp合并多目录TIFF mosaic？结果仅含最后一项

这个问题我之前碰到过不少，核心原因是你的脚本大概率是在循环里每次调用gdalwarp时，都直接把输出写入mergedM.tif——每一次新的处理都会覆盖掉之前的文件，最后自然就只剩最后一个mosaic的内容了。GDAL Warp默认不会自动累加合并，得用正确的姿势传递输入文件才行，给你两种靠谱的解决方案：

方案一：一次性传递所有输入文件（推荐，效率更高） #

这种方法适合文件数量不算特别多、内存足够的场景，直接把所有要合并的mosaic路径一次性传给gdalwarp，让它一次性完成合并，避免覆盖问题：

basepath=/home/rose/Desktop/b # 按需修改顶层目录
cd $basepath

# 用find收集所有mosaic文件（假设是.tif格式，可根据实际调整后缀）
# sort是为了保证合并顺序稳定，按需可去掉
mosaic_list=$(find $PWD -name "*.tif" | sort)

# 调用gdalwarp一次性合并所有文件，启用多线程加速（可选但推荐）
gdalwarp -multi -wo NUM_THREADS=ALL_CPUS $mosaic_list mergedM.tif

小提示： #

• 用find代替ls **更可靠，尤其是路径里有空格、特殊字符的情况，而且不需要额外开启bash的递归 glob 选项。
• 如果你的mosaic_list已经是提前定义好的列表变量，直接把它放在gdalwarp后面作为输入参数就行，不用循环。
• 可以添加-t_srs（指定目标投影）、-tr（指定目标分辨率）等参数，确保所有输入文件的空间参考一致，避免自动重采样带来的意外。

方案二：增量合并（适合大文件/内存不足场景） #

如果要合并的文件数量极大，一次性加载内存不够，可以用增量合并的方式，每次把一个新的mosaic合并到现有输出文件里：

basepath=/home/rose/Desktop/b
cd $basepath

# 先找到第一个mosaic作为初始合并文件
first_mosaic=$(find $PWD -name "*.tif" | sort | head -n1)
cp $first_mosaic mergedM.tif

# 循环处理剩余的mosaic文件
find $PWD -name "*.tif" | sort | tail -n +2 | while read mosaic_file; do
    # 用临时文件存储中间合并结果，避免直接覆盖出错
    temp_output="temp_merged.tif"
    gdalwarp mergedM.tif $mosaic_file $temp_output
    # 替换原输出文件
    mv $temp_output mergedM.tif
done

注意： #

• 每次合并都会生成临时文件，记得如果脚本中断，要手动清理临时文件。
• 这种方式的效率比一次性合并低，因为每次都要读写文件，适合内存紧张的情况。

内容的提问来源于stack exchange，提问作者Rose