RGB转XYZ再转Lab时B(Z)通道计算结果异常求助

• 检查RGB预处理的伽马校正逻辑
  确认BreakdownRGB中0-255转0-1后的伽马校正是否完全符合标准：

  • 当归一化后的值 ≤ 0.04045时，使用 val / 12.92
  • 当值 > 0.04045时，使用 Math.Pow((val + 0.055) / 1.055, 2.4)
    测试RGB(57,37,18)时，B值归一化后约为0.0706，属于后者分支，核对代码是否正确触发该逻辑——伽马校正错误会直接导致后续Z分量计算偏差。

• 核对D50标准的转换矩阵
  2°视场D50的sRGB转XYZ矩阵必须准确，尤其是Z分量对应的第三行：

  // 正确的D50转换矩阵
  float[,] matrix = new float[,] {
      { 0.4360747f, 0.3850649f, 0.1430804f },
      { 0.2225045f, 0.7168786f, 0.0606169f },
      { 0.0139322f, 0.0971045f, 0.7141733f }
  };
  

  避免混淆D65矩阵（Z分量系数差异极大），同时确认矩阵乘法的行列顺序正确，Z分量计算为 Z = correctedR * matrix[2,0] + correctedG * matrix[2,1] + correctedB * matrix[2,2]。

• 确认Lab转换的D50参考白值
  D50的参考白XYZ值为：X=96.422，Y=100.0，Z=82.521。检查XYZtoLAB中是否用此值做归一化，而非D65的参考白（X=95.047，Z=108.883）。错误的参考白会直接导致b分量偏离（如第一个测试案例中b值偏低，大概率是Z归一化时用错了参考白）。

• 排查数值精度与函数实现细节

  • 尝试将所有计算临时切换为double类型，避免float精度不足导致低亮度B值的计算误差被放大；
  • 核对Lab转换中的非线性变换函数f(t)：

    double f(double t) {
        return t > 0.008856 ? Math.Pow(t, 1.0/3) : (t * 3 * Math.Pow(29.0/6, 2)) + 4.0/29;
    }
    

    确保阈值（0.008856）和幂次（1/3）正确，不要误写为3次方或错误阈值。

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