| Revision: | 1.1 |
| Committed: | Sat Feb 22 02:07:21 2003 UTC (21 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); |