Matlab中色彩适配观察与Bradford变换图像应用技术咨询

一、灰度图像转彩色（观察色彩适配效果） #

如果只是想给灰度图添加伪彩色来观察色彩适配，MATLAB有几种简单且灵活的方式：

• 方法1：使用内置伪彩色映射
  把灰度图转为索引图像，再搭配不同的colormap（比如jet、viridis或者自定义色彩映射），示例代码：

  % 读取灰度图像（如果原图像是彩色，先转灰度）
  I = imread('snimka.jpeg');
  if size(I,3) == 3
      I_gray = rgb2gray(I);
  else
      I_gray = I;
  end
  figure(1), imshow(I_gray); title('原灰度图像');
  
  % 转伪彩色
  [X, map] = gray2ind(I_gray, 256);
  I_color = ind2rgb(X, viridis); % 用viridis色彩映射，也可以换成jet、hot等
  figure(2), imshow(I_color); title('伪彩色适配效果');
  

• 方法2：基于颜色空间的自定义色彩适配
  如果需要更精准的色彩适配（比如模拟特定光照下的色彩），可以把灰度图转到Lab颜色空间，手动调整a/b通道，再转回RGB：

  cform_srgb2lab = makecform('srgb2lab');
  cform_lab2srgb = makecform('lab2srgb');
  
  % 把灰度图扩展为3通道（Lab的L通道就是灰度）
  I_gray_3ch = repmat(I_gray, [1 1 3]);
  I_lab = applycform(I_gray_3ch, cform_srgb2lab);
  
  % 调整a/b通道来添加色彩（这里示例偏红绿色调，你可以自定义数值）
  I_lab(:,:,2) = I_lab(:,:,2) + 30; % a通道：红-绿轴
  I_lab(:,:,3) = I_lab(:,:,3) - 20; % b通道：黄-蓝轴
  
  I_color_adapted = applycform(I_lab, cform_lab2srgb);
  figure(3), imshow(I_color_adapted); title('自定义色彩适配效果');
  

二、Bradford变换的具体实现 #

你提到的Bradford矩阵Mbfd = [.8950 .2664 -.1614; -.7502 1.7135 .0367; .0389 -.0685 1.0296]是用于色彩适应的核心矩阵，操作需要在XYZ颜色空间下进行，具体步骤如下：

步骤拆解： #

1. 将RGB图像转换为XYZ颜色空间（因为Bradford变换是基于XYZ的）
2. 把图像数据转为double类型（避免整数运算的精度损失）
3. 对每个像素的XYZ值与Bradford矩阵进行矩阵乘法运算
4. 将变换后的XYZ数据转回RGB空间
5. 转换回图像的原始数据类型（比如uint8）

完整代码示例： #

% 1. 读取图像并转换为XYZ空间
I = imread('snimka.jpeg');
figure(1), imshow(I); title('原图像');

% 创建SRGB转XYZ的颜色转换结构体
cform_srgb2xyz = makecform('srgb2xyz');
I_xyz = applycform(I, cform_srgb2xyz);

% 2. 转换为double类型以便计算
I_xyz_double = im2double(I_xyz);

% 3. 定义Bradford矩阵并执行变换
Mbfd = [0.8950 0.2664 -0.1614; 
        -0.7502 1.7135 0.0367; 
        0.0389 -0.0685 1.0296];

% 注意：每个像素的XYZ是列向量，所以要做矩阵乘法：Mbfd * XYZ_col
% 先把图像维度调整为N×3（N是像素总数），再相乘，最后恢复维度
[h, w, ~] = size(I_xyz_double);
xyz_pixels = reshape(I_xyz_double, h*w, 3)'; % 转置为3×N的矩阵
transformed_xyz = Mbfd * xyz_pixels;
transformed_xyz = reshape(transformed_xyz', h, w, 3); % 转回h×w×3

% 4. 将变换后的XYZ转回SRGB空间
cform_xyz2srgb = makecform('xyz2srgb');
I_transformed = applycform(transformed_xyz, cform_xyz2srgb);

% 5. 处理超出范围的像素值（确保在[0,1]区间）并转回uint8
I_transformed = max(min(I_transformed, 1), 0);
I_transformed = im2uint8(I_transformed);

figure(2), imshow(I_transformed); title('Bradford变换后的图像');

关键注意点： #

• 矩阵乘法顺序：因为每个像素的XYZ是[X; Y; Z]列向量，所以要用Bradford矩阵左乘这个列向量，也就是Mbfd * [X; Y; Z]，代码中通过转置把像素排成3×N的矩阵，这样可以批量运算，效率更高。
• 数据类型处理：MATLAB图像默认是uint8，直接做矩阵乘法会导致精度丢失，所以必须先转成double类型，运算完成后再转回uint8，同时要把超出[0,1]范围的数值截断，避免图像出现异常色块。
• 颜色转换结构体：makecform函数生成的转换结构体是MATLAB中处理颜色空间转换的标准方式，确保转换符合色彩标准（比如SRGB的色域）。

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