| 1 |
{ RCSid $Id$ } |
| 2 |
{ |
| 3 |
Rcalc run to compute CRI and other factors from light source SPD |
| 4 |
convolved with standard curves. |
| 5 |
Used by scripts illumcal, reflcal. |
| 6 |
|
| 7 |
Input is summed wavelength convolution used below. |
| 8 |
} |
| 9 |
Lp = $1; { photopic lumens or luminance } |
| 10 |
Ls = $2; { scotopic lumens or luminance } |
| 11 |
x_k = $3/($3+$4+$5); { source x } |
| 12 |
y_k = $4/($3+$4+$5); { source y } |
| 13 |
u_k = u_XYZ($3,$4,$5); { source u } |
| 14 |
v_k = v_XYZ($3,$4,$5); { source v } |
| 15 |
c_k = c_uv(u_k,v_k); { source c } |
| 16 |
d_k = d_uv(u_k,v_k); { source d } |
| 17 |
x10_k = $6/($6+$7+$8); { source x10 } |
| 18 |
y10_k = $7/($6+$7+$8); { source y10 } |
| 19 |
Yi(i) = select(i, $10, $13, $16, $19, $22, $25, $28, $31) * 100 / $4; |
| 20 |
ui(i) = select(i, |
| 21 |
u_XYZ($9,$10,$11), { u for reference color 1 } |
| 22 |
u_XYZ($12,$13,$14), { u for reference color 2 } |
| 23 |
u_XYZ($15,$16,$17), { u for reference color 3 } |
| 24 |
u_XYZ($18,$19,$20), { u for reference color 4 } |
| 25 |
u_XYZ($21,$22,$23), { u for reference color 5 } |
| 26 |
u_XYZ($24,$25,$26), { u for reference color 6 } |
| 27 |
u_XYZ($27,$28,$29), { u for reference color 7 } |
| 28 |
u_XYZ($30,$31,$32) { u for reference color 8 } |
| 29 |
); |
| 30 |
vi(i) = select(i, |
| 31 |
v_XYZ($9,$10,$11), { v for reference color 1 } |
| 32 |
v_XYZ($12,$13,$14), { v for reference color 2 } |
| 33 |
v_XYZ($15,$16,$17), { v for reference color 3 } |
| 34 |
v_XYZ($18,$19,$20), { v for reference color 4 } |
| 35 |
v_XYZ($21,$22,$23), { v for reference color 5 } |
| 36 |
v_XYZ($24,$25,$26), { v for reference color 6 } |
| 37 |
v_XYZ($27,$28,$29), { v for reference color 7 } |
| 38 |
v_XYZ($30,$31,$32) { v for reference color 8 } |
| 39 |
); |
| 40 |
|
| 41 |
ci(i) = c_uv(ui(i),vi(i)); |
| 42 |
di(i) = d_uv(ui(i),vi(i)); |
| 43 |
|
| 44 |
u_XYZ(X,Y,Z) = 4*X/(X+15*Y+3*Z); |
| 45 |
v_XYZ(X,Y,Z) = 9*Y/(X+15*Y+3*Z); |
| 46 |
c_uv(u,v) = (4-u-10*v)/v; |
| 47 |
d_uv(u,v) = (1.708*v+.404-1.481*u)/v; |
