xyz2rgb

function Irgb = dos_xyz2rgb(Ixyz, primaries, white, tf_params)

Ixyz = im2double(Ixyz);

X = Ixyz(:,:,1);
Y = Ixyz(:,:,2);
Z = Ixyz(:,:,3);

R = primaries.R;
G = primaries.G;
B = primaries.B;

Irgb = zeros(size(Ixyz));
Irgb(:,:,1) = R(1)*X + R(2)*Y + R(3)*Z;
Irgb(:,:,2) = G(1)*X + G(2)*Y + G(3)*Z;
Irgb(:,:,3) = B(1)*X + B(2)*Y + B(3)*Z;

end