使用GDAL合并GeoTIFF遇黑白问题，求正确方法或Python插件

一、正确使用GDAL合并GeoTIFF的方法 #

你遇到的黑白问题，大概率是合并时默认只保留了第一个波段，或是丢失了调色板（Colormap）信息。以下是两种可靠的解决方式：

1. 命令行工具 gdal_merge.py #

GDAL自带的gdal_merge.py是专门处理TIFF合并的工具，通过添加参数可完整保留色彩：

• 若输入是单波段带调色板的TIFF：

  gdal_merge.py -o merged.tif -pct input1.tif input2.tif ...
  

  -pct参数会自动合并并保留所有输入文件的调色板，避免黑白化。
• 若输入是**多波段（如RGB）**的TIFF：

  gdal_merge.py -o merged.tif input1.tif input2.tif ...
  

  默认会保留所有波段，无需额外参数。如果出现只保留单波段的情况，检查输入文件是否确实为多波段，或添加-separate强制按波段顺序合并。

2. Python调用GDAL库 #

通过gdal.Warp实现合并，需确保保留元数据和色彩信息：

from osgeo import gdal

input_files = ["input1.tif", "input2.tif", "input3.tif"]
output_file = "merged.tif"

# 配置合并选项
warp_options = gdal.WarpOptions(
    format="GTiff",
    creationOptions=["COMPRESS=LZW"],
    copyMetadata=True,
    copyColorTable=True,  # 显式保留调色板
    resampleAlg="near"
)

# 执行合并
merged_ds = gdal.Warp(output_file, input_files, options=warp_options)
merged_ds = None  # 关闭数据集，确保写入完成

注意：若输入文件的CRS不一致，需在warp_options中添加dstSRS参数统一坐标系，否则合并结果会错位。

二、推荐的Python合并GeoTIFF插件 #

1. Rasterio #

Rasterio是基于GDAL的轻量化Python库，API更简洁，对色彩信息的处理更友好：

import rasterio
from rasterio.merge import merge

input_files = ["input1.tif", "input2.tif", "input3.tif"]
output_file = "merged.tif"

# 读取所有源文件
src_list = [rasterio.open(fp) for fp in input_files]

# 执行合并
mosaic_data, mosaic_transform = merge(src_list)

# 更新输出元数据
out_meta = src_list[0].meta.copy()
out_meta.update({
    "height": mosaic_data.shape[1],
    "width": mosaic_data.shape[2],
    "transform": mosaic_transform,
    "compress": "LZW"
})

# 写入输出文件
with rasterio.open(output_file, "w", **out_meta) as dest:
    dest.write(mosaic_data)

# 关闭源文件
for src in src_list:
    src.close()

Rasterio会自动继承输入文件的波段数、调色板及色彩空间，无需额外配置即可避免黑白问题。

2. GDAL Python绑定（原生） #

即前面提到的gdal.Warp方式，适合需要深度定制GDAL参数的场景，比如处理复杂投影或重采样需求。

常见问题排查 #

• 合并后仍为黑白：检查输入文件是否为单波段且无调色板，或合并时误删了波段；可通过gdalinfo input.tif查看文件的波段数和色彩信息。
• 合并结果错位：确认所有输入文件的CRS一致，不一致时需先通过gdalwarp统一坐标系再合并。

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