GDAL项目中的gdal_pansharpen工具栅格对齐问题解析

问题背景

在GDAL 3.7.x及以上版本中，用户报告gdal_pansharpen工具在处理某些高分辨率影像时出现异常。该工具用于将全色影像与多光谱影像融合，但在新版本中会产生大量"Access window out of range"错误，最终生成无效的全零值输出文件。

问题重现与分析

通过创建测试数据，可以稳定重现该问题：

1. 创建多光谱测试影像(6877x6516像素)
2. 创建全色测试影像(27511x26068像素)
3. 执行gdal_pansharpen融合处理

关键发现：

• 全色影像尺寸27511不是多光谱影像尺寸6877的整数倍(27511/4≈6878≠6877)
• 两幅影像的地理范围不完全一致
• GDAL 3.6.x版本能正确处理这种非整数倍关系，但3.7.x及以上版本会出现问题

技术原理

gdal_pansharpen工具的核心是将低分辨率多光谱影像与高分辨率全色影像进行融合，需要精确的栅格对齐和重采样。当输入影像的尺寸关系不是严格的整数倍时：

1. 工具需要计算重采样窗口
2. 在新版本中，窗口计算可能存在边界条件处理缺陷
3. 导致尝试读取超出源影像范围的像素
4. 产生"Access window out of range"错误
5. 最终输出无效数据

解决方案

目前确认的有效解决方案是确保输入影像满足：

1. 全色影像尺寸是多光谱影像的整数倍
2. 两幅影像的地理范围完全对齐

具体操作：

# 使用gdalwarp调整全色影像范围
gdalwarp -te [xmin ymin xmax ymax] pan.tif pan_adjusted.tif

深入建议

1. 输入检查：在使用gdal_pansharpen前，应验证输入影像的尺寸比例和地理范围
2. 版本选择：如需处理非标准尺寸比例的影像，可暂时使用GDAL 3.6.x版本
3. 输出验证：处理完成后应检查输出影像的统计值和直方图
4. 波段设置：确保多光谱影像的第4波段设置为"Undefined"而非"Alpha"

未来展望

该问题反映了GDAL在影像融合算法中对非标准尺寸比例处理的健壮性需要改进。建议开发团队：

1. 增强窗口计算的边界条件处理
2. 添加更严格的输入验证
3. 提供更明确的错误提示
4. 保持与旧版本的兼容性

对于专业用户，理解影像融合过程中的几何对齐要求至关重要，这有助于预防类似问题的发生。