如何用ImageMagick将RGB卫星影像叠加到单波段卫星影像上？

问题背景 #

手头有两张Himawari 8的GeoTIFF影像：

• IDE00435.202307240340.tif：1km分辨率可见光RGB影像
• IDE00409.202307240340.tif：500m分辨率单波段（全色）影像

使用gdal_pansharpen融合后效果过暗，转而尝试ImageMagick，已完成RGB影像的放大与对比度调整：

convert -level 20%,100%,2.0 -resize 200% IDE00435.202307240340.tif rgb-corrected.tif

但尝试用-colorize和-composite融合全色与RGB影像时失败，无法同时保留RGB色彩和全色影像的亮度细节。

解决方案 #

你之前用-colorize的思路不对，-colorize是给灰度图上色的工具，不适合做高分辨率全色与RGB的融合。正确的思路是替换RGB的亮度通道为全色影像，或者用合适的混合模式叠加细节，以下是两种可行方案：

方案一：HSV通道替换（精准保留色彩+全色细节） #

全色影像本质是亮度信息，我们可以把RGB转为HSV色彩空间（H=色相，S=饱和度，V=亮度），将V通道替换为全色影像后再转回RGB：

1. 先把全色影像调整到与RGB一致的亮度范围（避免亮度不匹配）：

# 先通过`identify -verbose rgb-corrected.tif`查看RGB影像位深，8位则用-depth 8，16位则用-depth 16
convert IDE00409.202307240340.tif -level 0%,100%,1.0 -depth 8 pan-normalized.tif

2. 执行通道替换与融合：

convert rgb-corrected.tif pan-normalized.tif \
  -colorspace HSV -channel R -separate +channel \
  \( -clone 0-1 pan-normalized.tif \) \
  -combine -colorspace RGB pansharpen-final.tif

• 解释：-colorspace HSV将RGB转为HSV，-separate拆分出H、S、V三个通道；然后把前两个通道（H、S）和归一化后的全色影像（作为新的V通道）重新组合，最后转回RGB色彩空间。

方案二：HardLight混合模式（快速叠加纹理细节） #

如果不需要精准的亮度替换，用HardLight混合模式可以快速把全色影像的纹理叠加到RGB上，同时保留色彩：

1. 调整全色影像的亮度范围：

convert IDE00409.202307240340.tif -level 0%,100%,1.0 pan-adjusted.tif

2. 执行混合融合：

convert pan-adjusted.tif rgb-corrected.tif -compose HardLight -composite pansharpen-hardlight.tif

• 解释：HardLight模式会根据全色影像的灰度值，对RGB影像做提亮或加深处理，从而保留全色的高分辨率纹理。

关键说明 #

• 确保两张影像地理配准完全对齐：可以用gdalinfo查看两张图的地理范围和投影，避免融合后错位；
• 保持位深一致：如果是16位卫星影像，务必在命令中加入对应-depth参数，避免色彩精度损失；
• 对比度调整：如果融合后亮度仍不理想，可以在全色影像的-level参数中微调，比如-level 5%,95%,1.2来增强对比度。

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