Page 61 - Computational Colour Science Using MATLAB

48 COMPUTING CIE TRISTIMULUS VALUES

disp(’use d65___64 for D65 and 1964 observer’); return;
end

% ignore the extreme wavelengths
cie = cie(8:37,:);

% interpolate to double the sampling rate
% the following three lines may be repeated
% for finer resolution
cie1(:,1) = pinterp(cie(:,1))’;
cie1(:,2) = pinterp(cie(:,2))’;
cie1(:,3) = pinterp(cie(:,3))’;

xy(:,1) = cie1(:,1)./sum(cie1’)’;
xy(:,2) = cie1(:,2)./sum(cie1’)’;