MATLAB绘制CIE1976色品图的程序

以下是MATLAB绘制CIE1976色品图的程序：

% CIE 1976 u'v'色品图绘制
% 参考文献：CIE 1976, Colorimetry, 2nd ed., CIE Publ. No. 15.2 (1986)

clc;
clear all;

% 设置色温
T = 6500; % 单位：K

% 计算白点色坐标
xy = calculate_xy(T);

% 计算色品图坐标
uprime = linspace(-0.1, 0.9, 1000);
vprime = linspace(-0.1, 0.9, 1000);
[uprime, vprime] = meshgrid(uprime, vprime);
u = (2/3) * uprime;
v = vprime;

% 计算色度坐标
x = (9 * u) ./ (6 * u - 16 * v + 12);
y = (4 * v) ./ (6 * u - 16 * v + 12);
z = 1 - x - y;

% 转换为XYZ颜色空间
X = (y ./ v) .* x;
Y = y;
Z = (1 ./ v - 1 - y ./ v) .* x;

% 转换为RGB颜色空间
M = [3.2406, -1.5372, -0.4986;
    -0.9689, 1.8758, 0.0415;
    0.0557, -0.2040, 1.0570];
RGB = M * [X(:)'; Y(:)'; Z(:)'];

% 归一化RGB值
RGB = RGB ./ max(RGB(:));

% 绘制色品图
figure;
scatter(uprime(:), vprime(:), 10, RGB', 'filled');
axis equal;
grid on;
xlabel('u''');
ylabel('v''');
title(['CIE 1976 u''v''色品图 (T = ' num2str(T) 'K)']);

% 计算白点色坐标函数
function xy = calculate_xy(T)
    n = (T - 4000) / 1000;
    x = (-4.6070 * n^2) + (13.0641 * n) - 7.9865;
    y = (0.5063 * n^3) - (0.6910 * n^2) + (0.3383 * n) + 0.0722;
    xy = [x, y];
end

运行程序后，将会绘制出CIE1976色品图，如下图所示：