TorchMetrics中ERGAS指标计算方法的修正说明

背景介绍

ERGAS(Erreur Relative Globale Adimensionnelle de Synthèse)是一种广泛应用于遥感图像融合质量评估的指标，特别是在全色锐化(Pansharpening)领域。该指标由Wald等人在1997年提出，用于评估多光谱图像与全色图像融合后的质量。

问题发现

在TorchMetrics项目的实现中，ERGAS指标的计算公式存在一个关键性的错误。原始实现中使用了乘法运算(* ratio)，而实际上应该使用除法运算(/ ratio)。这一错误源于对原始论文中参数含义的误解。

技术分析

ERGAS指标的标准计算公式为：

ERGAS = 100 * (h/l) * sqrt(1/N * sum(RMSE^2/mean^2))

其中：

• h/l表示高分辨率图像与低分辨率图像之间的像素尺寸比
• N是波段数量
• RMSE是均方根误差
• mean是参考图像的平均值

在实际应用中，h/l通常是一个小于1的分数值。例如，对于常见的卫星图像(如Pleiades-1A、WorldView-2、GeoEye-1等)，这个比值通常为1/4。因此，在代码实现时，应该使用整数比例参数ratio的倒数，即1/ratio。

影响评估

这一错误会导致计算出的ERGAS值比实际值大ratio^2倍。例如，当ratio=4时，计算结果会比正确值大16倍。这会严重影响评估结果的准确性，可能导致对图像融合算法性能的错误判断。

修正方案

正确的实现应该将乘法运算改为除法运算：

# 错误实现
ergas = 100 * ratio * torch.sqrt(torch.mean(rmse_per_band / mean_per_band))

# 正确实现
ergas = 100 / ratio * torch.sqrt(torch.mean(rmse_per_band / mean_per_band))

结论

ERGAS指标的正确计算对于遥感图像处理领域至关重要。TorchMetrics项目团队已经确认了这一错误，并将在后续版本中进行修正。用户在使用该指标时应注意这一计算差异，特别是在进行算法性能比较时，应确保使用正确的计算方法以获得准确的评估结果。