{ RCSid $Id: cielab.cal,v 1.2 2005/05/25 04:44:25 greg Exp $ }
{
	Compute Lab color space and color differences from XYZ primaries.

	5/25/02	G.Ward
}
Xw = 1; Yw = 1; Zw = 1;			{ default white }
sq(x) : x*x;
Ls(Yi) : if(Yi/Yw - .01, 116*(Yi/Yw)^(1/3) - 16, 903.3*Yi/Yw);
as(Xi,Yi,Zi) : 500*((Xi/Xw)^(1/3) - (Yi/Yw)^(1/3));
bs(Xi,Yi,Zi) : 200*((Yi/Yw)^(1/3) - (Zi/Zw)^(1/3));

					{ CIE 1976 Lab delta E formula }
dE76(X1,Y1,Z1,X2,Y2,Z2) : sqrt( sq(Ls(Y1)-Ls(Y2))
				+ sq(as(X1,Y1,Z1)-as(X2,Y2,Z2))
				+ sq(bs(X1,Y1,Z1)-bs(X2,Y2,Z2)) );

					{ CIE 1994 Lab delta E formula }
kL = 1; kC = 1; kH = 1;			{ viewing factors }
Cs(Xi,Yi,Zi) : sqrt( sq(as(Xi,Yi,Zi)) + sq(bs(Xi,Yi,Zi)) );
C1s(X1,Y1,Z1,X2,Y2,Z2) : sqrt(Cs(X1,Y1,Z1)*Cs(X2,Y2,Z2));
dHs(X1,Y1,Z1,X2,Y2,Z2) : sqrt( sq(as(X1,Y1,Z1)-as(X2,Y2,Z2))   
                               + sq(bs(X1,Y1,Z1)-bs(X2,Y2,Z2))
			       - sq(Cs(X1,Y1,Z1)-Cs(X2,Y2,Z2)) );
dE94(X1,Y1,Z1,X2,Y2,Z2) : sqrt( sq((Ls(X1,Y1,Z1)-Ls(X2,Y2,Z2))/kL)
			+ sq( (Cs(X1,Y1,Z1)-Cs(X2,Y2,Z2)) /
				(kC*(1+.048*C1s(X1,Y1,Z1,X2,Y2,Z2))) )
			+ sq( dHs(X1,Y1,Z1,X2,Y2,Z2) /
				(kH*(1+.014*C1s(X1,Y1,Z1,X2,Y2,Z2))) ) )