第三节：色彩管理矩阵——矩阵乘法原理与3x3色彩转换矩阵
好，咱们今天聊点硬核的。色彩管理矩阵，说白了就是一组数字，用来把颜色从一个空间搬到另一个空间。我刚开始接触这块时，也觉得不就是乘来乘去嘛，后来才发现——嗯，这里面的门道还真不少。
3.1 矩阵乘法：别怕，就是点积
先说说矩阵乘法。你想想看，一个3x3矩阵乘以一个3x1向量，结果还是一个3x1向量。怎么算的？
拿第一行举例：
[R'] [a11 a12 a13] [R]
[G'] = [a21 a22 a23] * [G]
[B'] [a31 a32 a33] [B]
结果就是：
R' = a11*R + a12*G + a13*B
G' = a21*R + a22*G + a23*B
B' = a31*R + a32*G + a33*B
说白了，每个新通道的值，都是原三个通道的加权和。我在项目中遇到过有人把矩阵写反了，结果颜色全乱套——所以顺序很重要，矩阵在左，向量在右。
小技巧： 我习惯把矩阵乘法想象成「揉面」——三个原色按不同比例揉在一起，得到新颜色。比例就是矩阵系数。
3.2 3x3色彩转换矩阵：从RGB到XYZ
RGB到XYZ的转换，是色彩管理里最基础的一步。为什么？因为RGB是设备相关的，而XYZ是设备无关的。你想想看，同一个RGB值，在不同显示器上颜色不一样，但XYZ值对应的颜色是唯一的。
标准转换矩阵长这样（以D65白点为例）：
[X] [0.4124564 0.3575761 0.1804375] [R_linear]
[Y] = [0.2126729 0.7151522 0.0721750] * [G_linear]
[Z] [0.0193339 0.1191920 0.9503041] [B_linear]
注意！这里的RGB必须是线性值。我见过太多人直接把sRGB的8位值往里塞，结果颜色偏得离谱。记住：先做Gamma校正，再做矩阵乘法。
避坑指南： 我曾经在项目中直接用非线性RGB做矩阵转换，结果Y值算出来全是错的。后来排查了两天才发现——矩阵乘法要求输入是线性光强，不是编码值。
3.3 白点适应：让颜色在不同光源下保持一致
白点适应，说白了就是「换灯泡」。同一个物体，在日光下和钨丝灯下看起来颜色不一样，但人眼会自动校正。色彩管理里也要做这件事。
最常用的方法是Bradford变换。它的核心思路是：先把XYZ转到一个「锥体响应」空间，调整白点，再转回来。
步骤是这样的：
- 把源白点XYZ_src和目标白点XYZ_dst都转到LMS空间
- 计算缩放因子：L_scale = L_dst / L_src，M和S同理
- 构造对角矩阵，乘到LMS值上
- 转回XYZ
举个例子，从D65到D50的白点适应：
# 源白点 D65: X=0.95047, Y=1.0, Z=1.08883
# 目标白点 D50: X=0.96422, Y=1.0, Z=0.82521
# Bradford变换矩阵
M_Bradford = [
[ 0.8951, 0.2664, -0.1614],
[-0.7502, 1.7135, 0.0367],
[ 0.0389, -0.0685, 1.0296]
]
# 先转LMS
LMS_src = M_Bradford * XYZ_src
LMS_dst = M_Bradford * XYZ_dst
# 计算缩放
scale = [LMS_dst[0]/LMS_src[0],
LMS_dst[1]/LMS_src[1],
LMS_dst[2]/LMS_src[2]]
# 构造对角矩阵并应用
M_adapt = M_Bradford_inv * diag(scale) * M_Bradford
核心要点： 白点适应不是可选的。如果你把D65的XYZ数据直接当成D50用，颜色会整体偏蓝。我建议在每次色彩空间转换时，都明确标注白点信息。
3.4 知识体系结构图
下面这张图，是我自己梳理的色彩管理矩阵核心逻辑。你看一遍，应该就能把今天讲的内容串起来了。
3.5 实战中的坑与经验
最后分享几个我踩过的坑：
- 矩阵精度问题： 我建议用double精度，float在多次变换后误差会累积。有一次我用float做三次矩阵连乘，结果颜色偏差肉眼可见。
- 白点标注： 每个XYZ数据都必须带白点信息。我见过一个项目，数据文件里没写白点，结果不同模块对同一数据理解不同，颜色对不上。
- 逆矩阵验证： 写完正向变换，一定要用逆矩阵验证一下。如果正反变换后误差超过0.001，说明矩阵有问题。
我的习惯： 每次做色彩矩阵变换，我都会在代码里加一个「白点检查」——输入一个纯白RGB，看输出的XYZ是否等于白点值。如果不等，说明矩阵或Gamma有问题。
好了，这一节的内容就到这。矩阵乘法、3x3转换矩阵、白点适应，这三块是色彩管理的基石。你把这些搞透了，后面再讲更复杂的色彩空间变换，就轻松多了。
公众号:蓝海资料掘金营，微信deep3321