6 |
|
} |
7 |
|
Xw = 1; Yw = 1; Zw = 1; { default white } |
8 |
|
sq(x) : x*x; |
9 |
< |
Ls(Yi) : if(Yi/Yw - .01, 116*(Yi/Yw)^(1/3) - 16, 903.3*Yi/Yw); |
10 |
< |
as(Xi,Yi,Zi) : 500*((Xi/Xw)^(1/3) - (Yi/Yw)^(1/3)); |
11 |
< |
bs(Xi,Yi,Zi) : 200*((Yi/Yw)^(1/3) - (Zi/Zw)^(1/3)); |
9 |
> |
DEL : 6/29; |
10 |
> |
bfunc(x) : if(x - DEL*DEL*DEL, x^(1/3), x/(3*DEL*DEL) + 2/3*DEL); |
11 |
> |
Ls(Yi) : 116*bfunc(Yi/Yw) - 16; |
12 |
> |
as(Xi,Yi,Zi) : 500*(bfunc(Xi/Xw) - bfunc(Yi/Yw)); |
13 |
> |
bs(Xi,Yi,Zi) : 200*(bfunc(Yi/Yw) - bfunc(Zi/Zw)); |
14 |
|
|
15 |
|
{ CIE 1976 Lab delta E formula } |
16 |
|
dE76(X1,Y1,Z1,X2,Y2,Z2) : sqrt( sq(Ls(Y1)-Ls(Y2)) |