Revision: | 1.1 |
Committed: | Sat Feb 22 02:07:21 2003 UTC (22 years, 2 months ago) by greg |
Branch: | MAIN |
CVS Tags: | rad5R2, rad4R2P2, rad5R0, rad5R1, rad3R7P2, rad3R7P1, rad4R2, rad4R1, rad4R0, rad3R5, rad3R6, rad3R6P1, rad3R8, rad3R9, rad4R2P1 |
Log Message: | Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release |
# | Content |
---|---|
1 | { |
2 | Peerless' brightness mapping function |
3 | |
4 | Mapping of Luminance to Brightness for CRT display. |
5 | Hand this file to pcomb(1) with the -f option. |
6 | The picture file should have been run previously through |
7 | the automatic exposure procedure of pfilt(1), and |
8 | pcomb should also be given -o option. Like so: |
9 | |
10 | pfilt input.pic | pcomb -f peerless.cal -o - > output.pic |
11 | |
12 | } |
13 | La_w : 0.5/le(1) * WE; |
14 | Lmax_v : 120; |
15 | La_v : Lmax_v/2; |
16 | b(La) : .338*La^.034; |
17 | k(La) : -1.5*log10(La) + 6.1; |
18 | B_w(L) : k(La_w) * L^b(La_w); |
19 | L_w(B) : (B/k(La_w))^(1/b(La_w)); |
20 | B_v(L) : k(La_v) * L^b(La_v); |
21 | L_v(B) : (B/k(La_v))^(1/b(La_v)); |
22 | Bmax_w : B_w(3*La_w); |
23 | Bmax_v : B_v(Lmax_v); |
24 | |
25 | lum_w = li(1)*WE; |
26 | |
27 | lum_v = L_v(Bmax_v/Bmax_w*B_w(lum_w)); |
28 | |
29 | mult = WE/Lmax_v*lum_v/lum_w; |
30 | |
31 | ro = mult*ri(1); |
32 | go = mult*gi(1); |
33 | bo = mult*bi(1); |