如何将Paraview中Lab空间的BuRd色卡导入MATLAB？

问题背景 #

需要将Paraview中标记为space="Lab"的BuRd色卡导入MATLAB作为自定义色卡，但直接使用lab2rgb转换会出现负值，直接当作RGB使用效果不符合预期。

关键原因 #

Paraview的Lab色卡中，<Point>的r/g/b值是归一化后的CIELAB分量，对应关系为：

• r → L分量：归一化到[0,1]（对应CIELAB的[0,100]）
• g → a分量：归一化到[0,1]（对应CIELAB的[-127,127]）
• b → b分量：归一化到[0,1]（对应CIELAB的[-127,127]）

MATLAB的lab2rgb函数需要标准CIELAB范围的数值，因此需要先对Paraview的归一化值进行反转换。

实现步骤 #

1. 提取色卡数据：从XML中提取所有<Point>的r/g/b值，整理为n×3的矩阵。
2. 反归一化到标准CIELAB范围：
   • L = r × 100
   • a = (g × 254) - 127
   • b = (b_val × 254) - 127
3. 转换为RGB颜色空间：使用lab2rgb函数将CIELAB值转换为MATLAB支持的RGB值（范围[0,1]）。
4. 生成自定义色卡：将转换后的RGB矩阵作为colormap使用。

MATLAB代码示例 #

% 提取Paraview BuRd色卡的r/g/b数据
paraview_lab_data = [
    0.019608, 0.188235, 0.380392;
    0.088504, 0.321107, 0.564937;
    0.163399, 0.444983, 0.697501;
    0.247059, 0.555709, 0.754095;
    0.420684, 0.676432, 0.818685;
    0.606459, 0.789773, 0.880277;
    0.761476, 0.868512, 0.924567;
    0.878047, 0.925721, 0.951942;
    0.969089, 0.966474, 0.964937;
    0.983852, 0.897578, 0.846828;
    0.982468, 0.800692, 0.706113;
    0.960323, 0.667820, 0.536332;
    0.894579, 0.503806, 0.399769;
    0.817070, 0.332180, 0.281046;
    0.728489, 0.155017, 0.197386;
    0.576932, 0.055363, 0.149250;
    0.403922, 0.000000, 0.121569;
];

% 反归一化到标准CIELAB范围
lab_data = zeros(size(paraview_lab_data));
lab_data(:,1) = paraview_lab_data(:,1) * 100;          % L分量：0-100
lab_data(:,2) = paraview_lab_data(:,2) * 254 - 127;    % a分量：-127到127
lab_data(:,3) = paraview_lab_data(:,3) * 254 - 127;    % b分量：-127到127

% 转换为RGB颜色空间
rgb_colormap = lab2rgb(lab_data);

% 可选：对RGB值进行裁剪，确保在[0,1]范围内（避免因精度问题出现微小负值）
rgb_colormap = clamp(rgb_colormap, 0, 1);

% 设置为当前色卡
colormap(rgb_colormap);

% 辅助函数：裁剪数值到指定范围
function out = clamp(in, min_val, max_val)
    out = max(min(in, max_val), min_val);
end

验证效果 #

运行上述代码后，通过colorbar查看生成的色卡，会与Paraview中的BuRd色卡一致：从深蓝色过渡到白色，再过渡到深红色。

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