1 | /*****************************************************************************/ |
2 | // Copyright 2006-2007 Adobe Systems Incorporated |
3 | // All Rights Reserved. |
4 | // |
5 | // NOTICE: Adobe permits you to use, modify, and distribute this file in |
6 | // accordance with the terms of the Adobe license agreement accompanying it. |
7 | /*****************************************************************************/ |
8 | |
9 | /* $Id: //mondo/dng_sdk_1_4/dng_sdk/source/dng_render.cpp#1 $ */ |
10 | /* $DateTime: 2012/05/30 13:28:51 $ */ |
11 | /* $Change: 832332 $ */ |
12 | /* $Author: tknoll $ */ |
13 | |
14 | /*****************************************************************************/ |
15 | |
16 | #include "dng_render.h" |
17 | |
18 | #include "dng_1d_table.h" |
19 | #include "dng_bottlenecks.h" |
20 | #include "dng_camera_profile.h" |
21 | #include "dng_color_space.h" |
22 | #include "dng_color_spec.h" |
23 | #include "dng_filter_task.h" |
24 | #include "dng_host.h" |
25 | #include "dng_image.h" |
26 | #include "dng_negative.h" |
27 | #include "dng_resample.h" |
28 | #include "dng_safe_arithmetic.h" |
29 | #include "dng_utils.h" |
30 | |
31 | /*****************************************************************************/ |
32 | |
33 | dng_function_exposure_ramp::dng_function_exposure_ramp (real64 white, |
34 | real64 black, |
35 | real64 minBlack) |
36 | |
37 | : fSlope ((white == black) ? 0.0f : 1.0 / (white - black)) |
38 | , fBlack (black) |
39 | |
40 | , fRadius (0.0) |
41 | , fQScale (0.0) |
42 | |
43 | { |
44 | if (fSlope == 0.0) |
45 | { |
46 | ThrowBadFormat (); |
47 | } |
48 | |
49 | const real64 kMaxCurveX = 0.5; // Fraction of minBlack. |
50 | |
51 | const real64 kMaxCurveY = 1.0 / 16.0; // Fraction of white. |
52 | |
53 | fRadius = Min_real64 (kMaxCurveX * minBlack, |
54 | kMaxCurveY / fSlope); |
55 | |
56 | if (fRadius > 0.0) |
57 | fQScale= fSlope / (4.0 * fRadius); |
58 | else |
59 | fQScale = 0.0; |
60 | |
61 | } |
62 | |
63 | /*****************************************************************************/ |
64 | |
65 | real64 dng_function_exposure_ramp::Evaluate (real64 x) const |
66 | { |
67 | |
68 | if (x <= fBlack - fRadius) |
69 | return 0.0; |
70 | |
71 | if (x >= fBlack + fRadius) |
72 | return Min_real64 ((x - fBlack) * fSlope, 1.0); |
73 | |
74 | real64 y = x - (fBlack - fRadius); |
75 | |
76 | return fQScale * y * y; |
77 | |
78 | } |
79 | |
80 | /*****************************************************************************/ |
81 | |
82 | dng_function_exposure_tone::dng_function_exposure_tone (real64 exposure) |
83 | |
84 | : fIsNOP (exposure >= 0.0) |
85 | |
86 | , fSlope (0.0) |
87 | |
88 | , a (0.0) |
89 | , b (0.0) |
90 | , c (0.0) |
91 | |
92 | { |
93 | |
94 | if (!fIsNOP) |
95 | { |
96 | |
97 | // Find slope to use for the all except the highest two f-stops. |
98 | |
99 | fSlope = pow (2.0, exposure); |
100 | |
101 | // Find quadradic parameters that match this darking at the crossover |
102 | // point, yet still map pure white to pure white. |
103 | |
104 | a = 16.0 / 9.0 * (1.0 - fSlope); |
105 | |
106 | b = fSlope - 0.5 * a; |
107 | |
108 | c = 1.0 - a - b; |
109 | |
110 | } |
111 | |
112 | } |
113 | |
114 | /*****************************************************************************/ |
115 | |
116 | real64 dng_function_exposure_tone::Evaluate (real64 x) const |
117 | { |
118 | |
119 | if (!fIsNOP) |
120 | { |
121 | |
122 | if (x <= 0.25) |
123 | x = x * fSlope; |
124 | |
125 | else |
126 | x = (a * x + b) * x + c; |
127 | |
128 | } |
129 | |
130 | return x; |
131 | |
132 | } |
133 | |
134 | /*****************************************************************************/ |
135 | |
136 | real64 dng_tone_curve_acr3_default::Evaluate (real64 x) const |
137 | { |
138 | |
139 | static const real32 kTable [] = |
140 | { |
141 | 0.00000f, 0.00078f, 0.00160f, 0.00242f, |
142 | 0.00314f, 0.00385f, 0.00460f, 0.00539f, |
143 | 0.00623f, 0.00712f, 0.00806f, 0.00906f, |
144 | 0.01012f, 0.01122f, 0.01238f, 0.01359f, |
145 | 0.01485f, 0.01616f, 0.01751f, 0.01890f, |
146 | 0.02033f, 0.02180f, 0.02331f, 0.02485f, |
147 | 0.02643f, 0.02804f, 0.02967f, 0.03134f, |
148 | 0.03303f, 0.03475f, 0.03648f, 0.03824f, |
149 | 0.04002f, 0.04181f, 0.04362f, 0.04545f, |
150 | 0.04730f, 0.04916f, 0.05103f, 0.05292f, |
151 | 0.05483f, 0.05675f, 0.05868f, 0.06063f, |
152 | 0.06259f, 0.06457f, 0.06655f, 0.06856f, |
153 | 0.07057f, 0.07259f, 0.07463f, 0.07668f, |
154 | 0.07874f, 0.08081f, 0.08290f, 0.08499f, |
155 | 0.08710f, 0.08921f, 0.09134f, 0.09348f, |
156 | 0.09563f, 0.09779f, 0.09996f, 0.10214f, |
157 | 0.10433f, 0.10652f, 0.10873f, 0.11095f, |
158 | 0.11318f, 0.11541f, 0.11766f, 0.11991f, |
159 | 0.12218f, 0.12445f, 0.12673f, 0.12902f, |
160 | 0.13132f, 0.13363f, 0.13595f, 0.13827f, |
161 | 0.14061f, 0.14295f, 0.14530f, 0.14765f, |
162 | 0.15002f, 0.15239f, 0.15477f, 0.15716f, |
163 | 0.15956f, 0.16197f, 0.16438f, 0.16680f, |
164 | 0.16923f, 0.17166f, 0.17410f, 0.17655f, |
165 | 0.17901f, 0.18148f, 0.18395f, 0.18643f, |
166 | 0.18891f, 0.19141f, 0.19391f, 0.19641f, |
167 | 0.19893f, 0.20145f, 0.20398f, 0.20651f, |
168 | 0.20905f, 0.21160f, 0.21416f, 0.21672f, |
169 | 0.21929f, 0.22185f, 0.22440f, 0.22696f, |
170 | 0.22950f, 0.23204f, 0.23458f, 0.23711f, |
171 | 0.23963f, 0.24215f, 0.24466f, 0.24717f, |
172 | 0.24967f, 0.25216f, 0.25465f, 0.25713f, |
173 | 0.25961f, 0.26208f, 0.26454f, 0.26700f, |
174 | 0.26945f, 0.27189f, 0.27433f, 0.27676f, |
175 | 0.27918f, 0.28160f, 0.28401f, 0.28641f, |
176 | 0.28881f, 0.29120f, 0.29358f, 0.29596f, |
177 | 0.29833f, 0.30069f, 0.30305f, 0.30540f, |
178 | 0.30774f, 0.31008f, 0.31241f, 0.31473f, |
179 | 0.31704f, 0.31935f, 0.32165f, 0.32395f, |
180 | 0.32623f, 0.32851f, 0.33079f, 0.33305f, |
181 | 0.33531f, 0.33756f, 0.33981f, 0.34205f, |
182 | 0.34428f, 0.34650f, 0.34872f, 0.35093f, |
183 | 0.35313f, 0.35532f, 0.35751f, 0.35969f, |
184 | 0.36187f, 0.36404f, 0.36620f, 0.36835f, |
185 | 0.37050f, 0.37264f, 0.37477f, 0.37689f, |
186 | 0.37901f, 0.38112f, 0.38323f, 0.38533f, |
187 | 0.38742f, 0.38950f, 0.39158f, 0.39365f, |
188 | 0.39571f, 0.39777f, 0.39982f, 0.40186f, |
189 | 0.40389f, 0.40592f, 0.40794f, 0.40996f, |
190 | 0.41197f, 0.41397f, 0.41596f, 0.41795f, |
191 | 0.41993f, 0.42191f, 0.42388f, 0.42584f, |
192 | 0.42779f, 0.42974f, 0.43168f, 0.43362f, |
193 | 0.43554f, 0.43747f, 0.43938f, 0.44129f, |
194 | 0.44319f, 0.44509f, 0.44698f, 0.44886f, |
195 | 0.45073f, 0.45260f, 0.45447f, 0.45632f, |
196 | 0.45817f, 0.46002f, 0.46186f, 0.46369f, |
197 | 0.46551f, 0.46733f, 0.46914f, 0.47095f, |
198 | 0.47275f, 0.47454f, 0.47633f, 0.47811f, |
199 | 0.47989f, 0.48166f, 0.48342f, 0.48518f, |
200 | 0.48693f, 0.48867f, 0.49041f, 0.49214f, |
201 | 0.49387f, 0.49559f, 0.49730f, 0.49901f, |
202 | 0.50072f, 0.50241f, 0.50410f, 0.50579f, |
203 | 0.50747f, 0.50914f, 0.51081f, 0.51247f, |
204 | 0.51413f, 0.51578f, 0.51742f, 0.51906f, |
205 | 0.52069f, 0.52232f, 0.52394f, 0.52556f, |
206 | 0.52717f, 0.52878f, 0.53038f, 0.53197f, |
207 | 0.53356f, 0.53514f, 0.53672f, 0.53829f, |
208 | 0.53986f, 0.54142f, 0.54297f, 0.54452f, |
209 | 0.54607f, 0.54761f, 0.54914f, 0.55067f, |
210 | 0.55220f, 0.55371f, 0.55523f, 0.55673f, |
211 | 0.55824f, 0.55973f, 0.56123f, 0.56271f, |
212 | 0.56420f, 0.56567f, 0.56715f, 0.56861f, |
213 | 0.57007f, 0.57153f, 0.57298f, 0.57443f, |
214 | 0.57587f, 0.57731f, 0.57874f, 0.58017f, |
215 | 0.58159f, 0.58301f, 0.58443f, 0.58583f, |
216 | 0.58724f, 0.58864f, 0.59003f, 0.59142f, |
217 | 0.59281f, 0.59419f, 0.59556f, 0.59694f, |
218 | 0.59830f, 0.59966f, 0.60102f, 0.60238f, |
219 | 0.60373f, 0.60507f, 0.60641f, 0.60775f, |
220 | 0.60908f, 0.61040f, 0.61173f, 0.61305f, |
221 | 0.61436f, 0.61567f, 0.61698f, 0.61828f, |
222 | 0.61957f, 0.62087f, 0.62216f, 0.62344f, |
223 | 0.62472f, 0.62600f, 0.62727f, 0.62854f, |
224 | 0.62980f, 0.63106f, 0.63232f, 0.63357f, |
225 | 0.63482f, 0.63606f, 0.63730f, 0.63854f, |
226 | 0.63977f, 0.64100f, 0.64222f, 0.64344f, |
227 | 0.64466f, 0.64587f, 0.64708f, 0.64829f, |
228 | 0.64949f, 0.65069f, 0.65188f, 0.65307f, |
229 | 0.65426f, 0.65544f, 0.65662f, 0.65779f, |
230 | 0.65897f, 0.66013f, 0.66130f, 0.66246f, |
231 | 0.66362f, 0.66477f, 0.66592f, 0.66707f, |
232 | 0.66821f, 0.66935f, 0.67048f, 0.67162f, |
233 | 0.67275f, 0.67387f, 0.67499f, 0.67611f, |
234 | 0.67723f, 0.67834f, 0.67945f, 0.68055f, |
235 | 0.68165f, 0.68275f, 0.68385f, 0.68494f, |
236 | 0.68603f, 0.68711f, 0.68819f, 0.68927f, |
237 | 0.69035f, 0.69142f, 0.69249f, 0.69355f, |
238 | 0.69461f, 0.69567f, 0.69673f, 0.69778f, |
239 | 0.69883f, 0.69988f, 0.70092f, 0.70196f, |
240 | 0.70300f, 0.70403f, 0.70506f, 0.70609f, |
241 | 0.70711f, 0.70813f, 0.70915f, 0.71017f, |
242 | 0.71118f, 0.71219f, 0.71319f, 0.71420f, |
243 | 0.71520f, 0.71620f, 0.71719f, 0.71818f, |
244 | 0.71917f, 0.72016f, 0.72114f, 0.72212f, |
245 | 0.72309f, 0.72407f, 0.72504f, 0.72601f, |
246 | 0.72697f, 0.72794f, 0.72890f, 0.72985f, |
247 | 0.73081f, 0.73176f, 0.73271f, 0.73365f, |
248 | 0.73460f, 0.73554f, 0.73647f, 0.73741f, |
249 | 0.73834f, 0.73927f, 0.74020f, 0.74112f, |
250 | 0.74204f, 0.74296f, 0.74388f, 0.74479f, |
251 | 0.74570f, 0.74661f, 0.74751f, 0.74842f, |
252 | 0.74932f, 0.75021f, 0.75111f, 0.75200f, |
253 | 0.75289f, 0.75378f, 0.75466f, 0.75555f, |
254 | 0.75643f, 0.75730f, 0.75818f, 0.75905f, |
255 | 0.75992f, 0.76079f, 0.76165f, 0.76251f, |
256 | 0.76337f, 0.76423f, 0.76508f, 0.76594f, |
257 | 0.76679f, 0.76763f, 0.76848f, 0.76932f, |
258 | 0.77016f, 0.77100f, 0.77183f, 0.77267f, |
259 | 0.77350f, 0.77432f, 0.77515f, 0.77597f, |
260 | 0.77680f, 0.77761f, 0.77843f, 0.77924f, |
261 | 0.78006f, 0.78087f, 0.78167f, 0.78248f, |
262 | 0.78328f, 0.78408f, 0.78488f, 0.78568f, |
263 | 0.78647f, 0.78726f, 0.78805f, 0.78884f, |
264 | 0.78962f, 0.79040f, 0.79118f, 0.79196f, |
265 | 0.79274f, 0.79351f, 0.79428f, 0.79505f, |
266 | 0.79582f, 0.79658f, 0.79735f, 0.79811f, |
267 | 0.79887f, 0.79962f, 0.80038f, 0.80113f, |
268 | 0.80188f, 0.80263f, 0.80337f, 0.80412f, |
269 | 0.80486f, 0.80560f, 0.80634f, 0.80707f, |
270 | 0.80780f, 0.80854f, 0.80926f, 0.80999f, |
271 | 0.81072f, 0.81144f, 0.81216f, 0.81288f, |
272 | 0.81360f, 0.81431f, 0.81503f, 0.81574f, |
273 | 0.81645f, 0.81715f, 0.81786f, 0.81856f, |
274 | 0.81926f, 0.81996f, 0.82066f, 0.82135f, |
275 | 0.82205f, 0.82274f, 0.82343f, 0.82412f, |
276 | 0.82480f, 0.82549f, 0.82617f, 0.82685f, |
277 | 0.82753f, 0.82820f, 0.82888f, 0.82955f, |
278 | 0.83022f, 0.83089f, 0.83155f, 0.83222f, |
279 | 0.83288f, 0.83354f, 0.83420f, 0.83486f, |
280 | 0.83552f, 0.83617f, 0.83682f, 0.83747f, |
281 | 0.83812f, 0.83877f, 0.83941f, 0.84005f, |
282 | 0.84069f, 0.84133f, 0.84197f, 0.84261f, |
283 | 0.84324f, 0.84387f, 0.84450f, 0.84513f, |
284 | 0.84576f, 0.84639f, 0.84701f, 0.84763f, |
285 | 0.84825f, 0.84887f, 0.84949f, 0.85010f, |
286 | 0.85071f, 0.85132f, 0.85193f, 0.85254f, |
287 | 0.85315f, 0.85375f, 0.85436f, 0.85496f, |
288 | 0.85556f, 0.85615f, 0.85675f, 0.85735f, |
289 | 0.85794f, 0.85853f, 0.85912f, 0.85971f, |
290 | 0.86029f, 0.86088f, 0.86146f, 0.86204f, |
291 | 0.86262f, 0.86320f, 0.86378f, 0.86435f, |
292 | 0.86493f, 0.86550f, 0.86607f, 0.86664f, |
293 | 0.86720f, 0.86777f, 0.86833f, 0.86889f, |
294 | 0.86945f, 0.87001f, 0.87057f, 0.87113f, |
295 | 0.87168f, 0.87223f, 0.87278f, 0.87333f, |
296 | 0.87388f, 0.87443f, 0.87497f, 0.87552f, |
297 | 0.87606f, 0.87660f, 0.87714f, 0.87768f, |
298 | 0.87821f, 0.87875f, 0.87928f, 0.87981f, |
299 | 0.88034f, 0.88087f, 0.88140f, 0.88192f, |
300 | 0.88244f, 0.88297f, 0.88349f, 0.88401f, |
301 | 0.88453f, 0.88504f, 0.88556f, 0.88607f, |
302 | 0.88658f, 0.88709f, 0.88760f, 0.88811f, |
303 | 0.88862f, 0.88912f, 0.88963f, 0.89013f, |
304 | 0.89063f, 0.89113f, 0.89163f, 0.89212f, |
305 | 0.89262f, 0.89311f, 0.89360f, 0.89409f, |
306 | 0.89458f, 0.89507f, 0.89556f, 0.89604f, |
307 | 0.89653f, 0.89701f, 0.89749f, 0.89797f, |
308 | 0.89845f, 0.89892f, 0.89940f, 0.89987f, |
309 | 0.90035f, 0.90082f, 0.90129f, 0.90176f, |
310 | 0.90222f, 0.90269f, 0.90316f, 0.90362f, |
311 | 0.90408f, 0.90454f, 0.90500f, 0.90546f, |
312 | 0.90592f, 0.90637f, 0.90683f, 0.90728f, |
313 | 0.90773f, 0.90818f, 0.90863f, 0.90908f, |
314 | 0.90952f, 0.90997f, 0.91041f, 0.91085f, |
315 | 0.91130f, 0.91173f, 0.91217f, 0.91261f, |
316 | 0.91305f, 0.91348f, 0.91392f, 0.91435f, |
317 | 0.91478f, 0.91521f, 0.91564f, 0.91606f, |
318 | 0.91649f, 0.91691f, 0.91734f, 0.91776f, |
319 | 0.91818f, 0.91860f, 0.91902f, 0.91944f, |
320 | 0.91985f, 0.92027f, 0.92068f, 0.92109f, |
321 | 0.92150f, 0.92191f, 0.92232f, 0.92273f, |
322 | 0.92314f, 0.92354f, 0.92395f, 0.92435f, |
323 | 0.92475f, 0.92515f, 0.92555f, 0.92595f, |
324 | 0.92634f, 0.92674f, 0.92713f, 0.92753f, |
325 | 0.92792f, 0.92831f, 0.92870f, 0.92909f, |
326 | 0.92947f, 0.92986f, 0.93025f, 0.93063f, |
327 | 0.93101f, 0.93139f, 0.93177f, 0.93215f, |
328 | 0.93253f, 0.93291f, 0.93328f, 0.93366f, |
329 | 0.93403f, 0.93440f, 0.93478f, 0.93515f, |
330 | 0.93551f, 0.93588f, 0.93625f, 0.93661f, |
331 | 0.93698f, 0.93734f, 0.93770f, 0.93807f, |
332 | 0.93843f, 0.93878f, 0.93914f, 0.93950f, |
333 | 0.93986f, 0.94021f, 0.94056f, 0.94092f, |
334 | 0.94127f, 0.94162f, 0.94197f, 0.94231f, |
335 | 0.94266f, 0.94301f, 0.94335f, 0.94369f, |
336 | 0.94404f, 0.94438f, 0.94472f, 0.94506f, |
337 | 0.94540f, 0.94573f, 0.94607f, 0.94641f, |
338 | 0.94674f, 0.94707f, 0.94740f, 0.94774f, |
339 | 0.94807f, 0.94839f, 0.94872f, 0.94905f, |
340 | 0.94937f, 0.94970f, 0.95002f, 0.95035f, |
341 | 0.95067f, 0.95099f, 0.95131f, 0.95163f, |
342 | 0.95194f, 0.95226f, 0.95257f, 0.95289f, |
343 | 0.95320f, 0.95351f, 0.95383f, 0.95414f, |
344 | 0.95445f, 0.95475f, 0.95506f, 0.95537f, |
345 | 0.95567f, 0.95598f, 0.95628f, 0.95658f, |
346 | 0.95688f, 0.95718f, 0.95748f, 0.95778f, |
347 | 0.95808f, 0.95838f, 0.95867f, 0.95897f, |
348 | 0.95926f, 0.95955f, 0.95984f, 0.96013f, |
349 | 0.96042f, 0.96071f, 0.96100f, 0.96129f, |
350 | 0.96157f, 0.96186f, 0.96214f, 0.96242f, |
351 | 0.96271f, 0.96299f, 0.96327f, 0.96355f, |
352 | 0.96382f, 0.96410f, 0.96438f, 0.96465f, |
353 | 0.96493f, 0.96520f, 0.96547f, 0.96574f, |
354 | 0.96602f, 0.96629f, 0.96655f, 0.96682f, |
355 | 0.96709f, 0.96735f, 0.96762f, 0.96788f, |
356 | 0.96815f, 0.96841f, 0.96867f, 0.96893f, |
357 | 0.96919f, 0.96945f, 0.96971f, 0.96996f, |
358 | 0.97022f, 0.97047f, 0.97073f, 0.97098f, |
359 | 0.97123f, 0.97149f, 0.97174f, 0.97199f, |
360 | 0.97223f, 0.97248f, 0.97273f, 0.97297f, |
361 | 0.97322f, 0.97346f, 0.97371f, 0.97395f, |
362 | 0.97419f, 0.97443f, 0.97467f, 0.97491f, |
363 | 0.97515f, 0.97539f, 0.97562f, 0.97586f, |
364 | 0.97609f, 0.97633f, 0.97656f, 0.97679f, |
365 | 0.97702f, 0.97725f, 0.97748f, 0.97771f, |
366 | 0.97794f, 0.97817f, 0.97839f, 0.97862f, |
367 | 0.97884f, 0.97907f, 0.97929f, 0.97951f, |
368 | 0.97973f, 0.97995f, 0.98017f, 0.98039f, |
369 | 0.98061f, 0.98082f, 0.98104f, 0.98125f, |
370 | 0.98147f, 0.98168f, 0.98189f, 0.98211f, |
371 | 0.98232f, 0.98253f, 0.98274f, 0.98295f, |
372 | 0.98315f, 0.98336f, 0.98357f, 0.98377f, |
373 | 0.98398f, 0.98418f, 0.98438f, 0.98458f, |
374 | 0.98478f, 0.98498f, 0.98518f, 0.98538f, |
375 | 0.98558f, 0.98578f, 0.98597f, 0.98617f, |
376 | 0.98636f, 0.98656f, 0.98675f, 0.98694f, |
377 | 0.98714f, 0.98733f, 0.98752f, 0.98771f, |
378 | 0.98789f, 0.98808f, 0.98827f, 0.98845f, |
379 | 0.98864f, 0.98882f, 0.98901f, 0.98919f, |
380 | 0.98937f, 0.98955f, 0.98973f, 0.98991f, |
381 | 0.99009f, 0.99027f, 0.99045f, 0.99063f, |
382 | 0.99080f, 0.99098f, 0.99115f, 0.99133f, |
383 | 0.99150f, 0.99167f, 0.99184f, 0.99201f, |
384 | 0.99218f, 0.99235f, 0.99252f, 0.99269f, |
385 | 0.99285f, 0.99302f, 0.99319f, 0.99335f, |
386 | 0.99351f, 0.99368f, 0.99384f, 0.99400f, |
387 | 0.99416f, 0.99432f, 0.99448f, 0.99464f, |
388 | 0.99480f, 0.99495f, 0.99511f, 0.99527f, |
389 | 0.99542f, 0.99558f, 0.99573f, 0.99588f, |
390 | 0.99603f, 0.99619f, 0.99634f, 0.99649f, |
391 | 0.99664f, 0.99678f, 0.99693f, 0.99708f, |
392 | 0.99722f, 0.99737f, 0.99751f, 0.99766f, |
393 | 0.99780f, 0.99794f, 0.99809f, 0.99823f, |
394 | 0.99837f, 0.99851f, 0.99865f, 0.99879f, |
395 | 0.99892f, 0.99906f, 0.99920f, 0.99933f, |
396 | 0.99947f, 0.99960f, 0.99974f, 0.99987f, |
397 | 1.00000f |
398 | }; |
399 | |
400 | const uint32 kTableSize = sizeof (kTable ) / |
401 | sizeof (kTable [0]); |
402 | |
403 | real32 y = (real32) x * (real32) (kTableSize - 1); |
404 | |
405 | int32 index = Pin_int32 (0, (int32) y, kTableSize - 2); |
406 | |
407 | real32 fract = y - (real32) index; |
408 | |
409 | return kTable [index ] * (1.0f - fract) + |
410 | kTable [index + 1] * ( fract); |
411 | |
412 | } |
413 | |
414 | /*****************************************************************************/ |
415 | |
416 | real64 dng_tone_curve_acr3_default::EvaluateInverse (real64 x) const |
417 | { |
418 | |
419 | static const real32 kTable [] = |
420 | { |
421 | 0.00000f, 0.00121f, 0.00237f, 0.00362f, |
422 | 0.00496f, 0.00621f, 0.00738f, 0.00848f, |
423 | 0.00951f, 0.01048f, 0.01139f, 0.01227f, |
424 | 0.01312f, 0.01393f, 0.01471f, 0.01547f, |
425 | 0.01620f, 0.01692f, 0.01763f, 0.01831f, |
426 | 0.01899f, 0.01965f, 0.02030f, 0.02094f, |
427 | 0.02157f, 0.02218f, 0.02280f, 0.02340f, |
428 | 0.02399f, 0.02458f, 0.02517f, 0.02574f, |
429 | 0.02631f, 0.02688f, 0.02744f, 0.02800f, |
430 | 0.02855f, 0.02910f, 0.02965f, 0.03019f, |
431 | 0.03072f, 0.03126f, 0.03179f, 0.03232f, |
432 | 0.03285f, 0.03338f, 0.03390f, 0.03442f, |
433 | 0.03493f, 0.03545f, 0.03596f, 0.03647f, |
434 | 0.03698f, 0.03749f, 0.03799f, 0.03849f, |
435 | 0.03899f, 0.03949f, 0.03998f, 0.04048f, |
436 | 0.04097f, 0.04146f, 0.04195f, 0.04244f, |
437 | 0.04292f, 0.04341f, 0.04389f, 0.04437f, |
438 | 0.04485f, 0.04533f, 0.04580f, 0.04628f, |
439 | 0.04675f, 0.04722f, 0.04769f, 0.04816f, |
440 | 0.04863f, 0.04910f, 0.04956f, 0.05003f, |
441 | 0.05049f, 0.05095f, 0.05141f, 0.05187f, |
442 | 0.05233f, 0.05278f, 0.05324f, 0.05370f, |
443 | 0.05415f, 0.05460f, 0.05505f, 0.05551f, |
444 | 0.05595f, 0.05640f, 0.05685f, 0.05729f, |
445 | 0.05774f, 0.05818f, 0.05863f, 0.05907f, |
446 | 0.05951f, 0.05995f, 0.06039f, 0.06083f, |
447 | 0.06126f, 0.06170f, 0.06214f, 0.06257f, |
448 | 0.06301f, 0.06344f, 0.06388f, 0.06431f, |
449 | 0.06474f, 0.06517f, 0.06560f, 0.06602f, |
450 | 0.06645f, 0.06688f, 0.06731f, 0.06773f, |
451 | 0.06815f, 0.06858f, 0.06900f, 0.06943f, |
452 | 0.06985f, 0.07027f, 0.07069f, 0.07111f, |
453 | 0.07152f, 0.07194f, 0.07236f, 0.07278f, |
454 | 0.07319f, 0.07361f, 0.07402f, 0.07444f, |
455 | 0.07485f, 0.07526f, 0.07567f, 0.07608f, |
456 | 0.07650f, 0.07691f, 0.07732f, 0.07772f, |
457 | 0.07813f, 0.07854f, 0.07895f, 0.07935f, |
458 | 0.07976f, 0.08016f, 0.08057f, 0.08098f, |
459 | 0.08138f, 0.08178f, 0.08218f, 0.08259f, |
460 | 0.08299f, 0.08339f, 0.08379f, 0.08419f, |
461 | 0.08459f, 0.08499f, 0.08539f, 0.08578f, |
462 | 0.08618f, 0.08657f, 0.08697f, 0.08737f, |
463 | 0.08776f, 0.08816f, 0.08855f, 0.08894f, |
464 | 0.08934f, 0.08973f, 0.09012f, 0.09051f, |
465 | 0.09091f, 0.09130f, 0.09169f, 0.09208f, |
466 | 0.09247f, 0.09286f, 0.09324f, 0.09363f, |
467 | 0.09402f, 0.09440f, 0.09479f, 0.09518f, |
468 | 0.09556f, 0.09595f, 0.09633f, 0.09672f, |
469 | 0.09710f, 0.09749f, 0.09787f, 0.09825f, |
470 | 0.09863f, 0.09901f, 0.09939f, 0.09978f, |
471 | 0.10016f, 0.10054f, 0.10092f, 0.10130f, |
472 | 0.10167f, 0.10205f, 0.10243f, 0.10281f, |
473 | 0.10319f, 0.10356f, 0.10394f, 0.10432f, |
474 | 0.10469f, 0.10507f, 0.10544f, 0.10582f, |
475 | 0.10619f, 0.10657f, 0.10694f, 0.10731f, |
476 | 0.10768f, 0.10806f, 0.10843f, 0.10880f, |
477 | 0.10917f, 0.10954f, 0.10991f, 0.11029f, |
478 | 0.11066f, 0.11103f, 0.11141f, 0.11178f, |
479 | 0.11215f, 0.11253f, 0.11290f, 0.11328f, |
480 | 0.11365f, 0.11403f, 0.11440f, 0.11478f, |
481 | 0.11516f, 0.11553f, 0.11591f, 0.11629f, |
482 | 0.11666f, 0.11704f, 0.11742f, 0.11780f, |
483 | 0.11818f, 0.11856f, 0.11894f, 0.11932f, |
484 | 0.11970f, 0.12008f, 0.12046f, 0.12084f, |
485 | 0.12122f, 0.12161f, 0.12199f, 0.12237f, |
486 | 0.12276f, 0.12314f, 0.12352f, 0.12391f, |
487 | 0.12429f, 0.12468f, 0.12506f, 0.12545f, |
488 | 0.12583f, 0.12622f, 0.12661f, 0.12700f, |
489 | 0.12738f, 0.12777f, 0.12816f, 0.12855f, |
490 | 0.12894f, 0.12933f, 0.12972f, 0.13011f, |
491 | 0.13050f, 0.13089f, 0.13129f, 0.13168f, |
492 | 0.13207f, 0.13247f, 0.13286f, 0.13325f, |
493 | 0.13365f, 0.13404f, 0.13444f, 0.13483f, |
494 | 0.13523f, 0.13563f, 0.13603f, 0.13642f, |
495 | 0.13682f, 0.13722f, 0.13762f, 0.13802f, |
496 | 0.13842f, 0.13882f, 0.13922f, 0.13962f, |
497 | 0.14003f, 0.14043f, 0.14083f, 0.14124f, |
498 | 0.14164f, 0.14204f, 0.14245f, 0.14285f, |
499 | 0.14326f, 0.14366f, 0.14407f, 0.14448f, |
500 | 0.14489f, 0.14530f, 0.14570f, 0.14611f, |
501 | 0.14652f, 0.14693f, 0.14734f, 0.14776f, |
502 | 0.14817f, 0.14858f, 0.14900f, 0.14941f, |
503 | 0.14982f, 0.15024f, 0.15065f, 0.15107f, |
504 | 0.15148f, 0.15190f, 0.15232f, 0.15274f, |
505 | 0.15316f, 0.15357f, 0.15399f, 0.15441f, |
506 | 0.15483f, 0.15526f, 0.15568f, 0.15610f, |
507 | 0.15652f, 0.15695f, 0.15737f, 0.15779f, |
508 | 0.15822f, 0.15864f, 0.15907f, 0.15950f, |
509 | 0.15992f, 0.16035f, 0.16078f, 0.16121f, |
510 | 0.16164f, 0.16207f, 0.16250f, 0.16293f, |
511 | 0.16337f, 0.16380f, 0.16423f, 0.16467f, |
512 | 0.16511f, 0.16554f, 0.16598f, 0.16641f, |
513 | 0.16685f, 0.16729f, 0.16773f, 0.16816f, |
514 | 0.16860f, 0.16904f, 0.16949f, 0.16993f, |
515 | 0.17037f, 0.17081f, 0.17126f, 0.17170f, |
516 | 0.17215f, 0.17259f, 0.17304f, 0.17349f, |
517 | 0.17393f, 0.17438f, 0.17483f, 0.17528f, |
518 | 0.17573f, 0.17619f, 0.17664f, 0.17709f, |
519 | 0.17754f, 0.17799f, 0.17845f, 0.17890f, |
520 | 0.17936f, 0.17982f, 0.18028f, 0.18073f, |
521 | 0.18119f, 0.18165f, 0.18211f, 0.18257f, |
522 | 0.18303f, 0.18350f, 0.18396f, 0.18442f, |
523 | 0.18489f, 0.18535f, 0.18582f, 0.18629f, |
524 | 0.18676f, 0.18723f, 0.18770f, 0.18817f, |
525 | 0.18864f, 0.18911f, 0.18958f, 0.19005f, |
526 | 0.19053f, 0.19100f, 0.19147f, 0.19195f, |
527 | 0.19243f, 0.19291f, 0.19339f, 0.19387f, |
528 | 0.19435f, 0.19483f, 0.19531f, 0.19579f, |
529 | 0.19627f, 0.19676f, 0.19724f, 0.19773f, |
530 | 0.19821f, 0.19870f, 0.19919f, 0.19968f, |
531 | 0.20017f, 0.20066f, 0.20115f, 0.20164f, |
532 | 0.20214f, 0.20263f, 0.20313f, 0.20362f, |
533 | 0.20412f, 0.20462f, 0.20512f, 0.20561f, |
534 | 0.20611f, 0.20662f, 0.20712f, 0.20762f, |
535 | 0.20812f, 0.20863f, 0.20913f, 0.20964f, |
536 | 0.21015f, 0.21066f, 0.21117f, 0.21168f, |
537 | 0.21219f, 0.21270f, 0.21321f, 0.21373f, |
538 | 0.21424f, 0.21476f, 0.21527f, 0.21579f, |
539 | 0.21631f, 0.21683f, 0.21735f, 0.21787f, |
540 | 0.21839f, 0.21892f, 0.21944f, 0.21997f, |
541 | 0.22049f, 0.22102f, 0.22155f, 0.22208f, |
542 | 0.22261f, 0.22314f, 0.22367f, 0.22420f, |
543 | 0.22474f, 0.22527f, 0.22581f, 0.22634f, |
544 | 0.22688f, 0.22742f, 0.22796f, 0.22850f, |
545 | 0.22905f, 0.22959f, 0.23013f, 0.23068f, |
546 | 0.23123f, 0.23178f, 0.23232f, 0.23287f, |
547 | 0.23343f, 0.23398f, 0.23453f, 0.23508f, |
548 | 0.23564f, 0.23620f, 0.23675f, 0.23731f, |
549 | 0.23787f, 0.23843f, 0.23899f, 0.23956f, |
550 | 0.24012f, 0.24069f, 0.24125f, 0.24182f, |
551 | 0.24239f, 0.24296f, 0.24353f, 0.24410f, |
552 | 0.24468f, 0.24525f, 0.24582f, 0.24640f, |
553 | 0.24698f, 0.24756f, 0.24814f, 0.24872f, |
554 | 0.24931f, 0.24989f, 0.25048f, 0.25106f, |
555 | 0.25165f, 0.25224f, 0.25283f, 0.25342f, |
556 | 0.25401f, 0.25460f, 0.25520f, 0.25579f, |
557 | 0.25639f, 0.25699f, 0.25759f, 0.25820f, |
558 | 0.25880f, 0.25940f, 0.26001f, 0.26062f, |
559 | 0.26122f, 0.26183f, 0.26244f, 0.26306f, |
560 | 0.26367f, 0.26429f, 0.26490f, 0.26552f, |
561 | 0.26614f, 0.26676f, 0.26738f, 0.26800f, |
562 | 0.26863f, 0.26925f, 0.26988f, 0.27051f, |
563 | 0.27114f, 0.27177f, 0.27240f, 0.27303f, |
564 | 0.27367f, 0.27431f, 0.27495f, 0.27558f, |
565 | 0.27623f, 0.27687f, 0.27751f, 0.27816f, |
566 | 0.27881f, 0.27945f, 0.28011f, 0.28076f, |
567 | 0.28141f, 0.28207f, 0.28272f, 0.28338f, |
568 | 0.28404f, 0.28470f, 0.28536f, 0.28602f, |
569 | 0.28669f, 0.28736f, 0.28802f, 0.28869f, |
570 | 0.28937f, 0.29004f, 0.29071f, 0.29139f, |
571 | 0.29207f, 0.29274f, 0.29342f, 0.29410f, |
572 | 0.29479f, 0.29548f, 0.29616f, 0.29685f, |
573 | 0.29754f, 0.29823f, 0.29893f, 0.29962f, |
574 | 0.30032f, 0.30102f, 0.30172f, 0.30242f, |
575 | 0.30312f, 0.30383f, 0.30453f, 0.30524f, |
576 | 0.30595f, 0.30667f, 0.30738f, 0.30809f, |
577 | 0.30881f, 0.30953f, 0.31025f, 0.31097f, |
578 | 0.31170f, 0.31242f, 0.31315f, 0.31388f, |
579 | 0.31461f, 0.31534f, 0.31608f, 0.31682f, |
580 | 0.31755f, 0.31829f, 0.31904f, 0.31978f, |
581 | 0.32053f, 0.32127f, 0.32202f, 0.32277f, |
582 | 0.32353f, 0.32428f, 0.32504f, 0.32580f, |
583 | 0.32656f, 0.32732f, 0.32808f, 0.32885f, |
584 | 0.32962f, 0.33039f, 0.33116f, 0.33193f, |
585 | 0.33271f, 0.33349f, 0.33427f, 0.33505f, |
586 | 0.33583f, 0.33662f, 0.33741f, 0.33820f, |
587 | 0.33899f, 0.33978f, 0.34058f, 0.34138f, |
588 | 0.34218f, 0.34298f, 0.34378f, 0.34459f, |
589 | 0.34540f, 0.34621f, 0.34702f, 0.34783f, |
590 | 0.34865f, 0.34947f, 0.35029f, 0.35111f, |
591 | 0.35194f, 0.35277f, 0.35360f, 0.35443f, |
592 | 0.35526f, 0.35610f, 0.35694f, 0.35778f, |
593 | 0.35862f, 0.35946f, 0.36032f, 0.36117f, |
594 | 0.36202f, 0.36287f, 0.36372f, 0.36458f, |
595 | 0.36545f, 0.36631f, 0.36718f, 0.36805f, |
596 | 0.36891f, 0.36979f, 0.37066f, 0.37154f, |
597 | 0.37242f, 0.37331f, 0.37419f, 0.37507f, |
598 | 0.37596f, 0.37686f, 0.37775f, 0.37865f, |
599 | 0.37955f, 0.38045f, 0.38136f, 0.38227f, |
600 | 0.38317f, 0.38409f, 0.38500f, 0.38592f, |
601 | 0.38684f, 0.38776f, 0.38869f, 0.38961f, |
602 | 0.39055f, 0.39148f, 0.39242f, 0.39335f, |
603 | 0.39430f, 0.39524f, 0.39619f, 0.39714f, |
604 | 0.39809f, 0.39904f, 0.40000f, 0.40097f, |
605 | 0.40193f, 0.40289f, 0.40386f, 0.40483f, |
606 | 0.40581f, 0.40679f, 0.40777f, 0.40875f, |
607 | 0.40974f, 0.41073f, 0.41172f, 0.41272f, |
608 | 0.41372f, 0.41472f, 0.41572f, 0.41673f, |
609 | 0.41774f, 0.41875f, 0.41977f, 0.42079f, |
610 | 0.42181f, 0.42284f, 0.42386f, 0.42490f, |
611 | 0.42594f, 0.42697f, 0.42801f, 0.42906f, |
612 | 0.43011f, 0.43116f, 0.43222f, 0.43327f, |
613 | 0.43434f, 0.43540f, 0.43647f, 0.43754f, |
614 | 0.43862f, 0.43970f, 0.44077f, 0.44186f, |
615 | 0.44295f, 0.44404f, 0.44514f, 0.44624f, |
616 | 0.44734f, 0.44845f, 0.44956f, 0.45068f, |
617 | 0.45179f, 0.45291f, 0.45404f, 0.45516f, |
618 | 0.45630f, 0.45744f, 0.45858f, 0.45972f, |
619 | 0.46086f, 0.46202f, 0.46318f, 0.46433f, |
620 | 0.46550f, 0.46667f, 0.46784f, 0.46901f, |
621 | 0.47019f, 0.47137f, 0.47256f, 0.47375f, |
622 | 0.47495f, 0.47615f, 0.47735f, 0.47856f, |
623 | 0.47977f, 0.48099f, 0.48222f, 0.48344f, |
624 | 0.48467f, 0.48590f, 0.48714f, 0.48838f, |
625 | 0.48963f, 0.49088f, 0.49213f, 0.49340f, |
626 | 0.49466f, 0.49593f, 0.49721f, 0.49849f, |
627 | 0.49977f, 0.50106f, 0.50236f, 0.50366f, |
628 | 0.50496f, 0.50627f, 0.50758f, 0.50890f, |
629 | 0.51023f, 0.51155f, 0.51289f, 0.51422f, |
630 | 0.51556f, 0.51692f, 0.51827f, 0.51964f, |
631 | 0.52100f, 0.52237f, 0.52374f, 0.52512f, |
632 | 0.52651f, 0.52790f, 0.52930f, 0.53070f, |
633 | 0.53212f, 0.53353f, 0.53495f, 0.53638f, |
634 | 0.53781f, 0.53925f, 0.54070f, 0.54214f, |
635 | 0.54360f, 0.54506f, 0.54653f, 0.54800f, |
636 | 0.54949f, 0.55098f, 0.55247f, 0.55396f, |
637 | 0.55548f, 0.55699f, 0.55851f, 0.56003f, |
638 | 0.56156f, 0.56310f, 0.56464f, 0.56621f, |
639 | 0.56777f, 0.56933f, 0.57091f, 0.57248f, |
640 | 0.57407f, 0.57568f, 0.57727f, 0.57888f, |
641 | 0.58050f, 0.58213f, 0.58376f, 0.58541f, |
642 | 0.58705f, 0.58871f, 0.59037f, 0.59204f, |
643 | 0.59373f, 0.59541f, 0.59712f, 0.59882f, |
644 | 0.60052f, 0.60226f, 0.60399f, 0.60572f, |
645 | 0.60748f, 0.60922f, 0.61099f, 0.61276f, |
646 | 0.61455f, 0.61635f, 0.61814f, 0.61996f, |
647 | 0.62178f, 0.62361f, 0.62545f, 0.62730f, |
648 | 0.62917f, 0.63104f, 0.63291f, 0.63480f, |
649 | 0.63671f, 0.63862f, 0.64054f, 0.64249f, |
650 | 0.64443f, 0.64638f, 0.64835f, 0.65033f, |
651 | 0.65232f, 0.65433f, 0.65633f, 0.65836f, |
652 | 0.66041f, 0.66245f, 0.66452f, 0.66660f, |
653 | 0.66868f, 0.67078f, 0.67290f, 0.67503f, |
654 | 0.67717f, 0.67932f, 0.68151f, 0.68368f, |
655 | 0.68587f, 0.68809f, 0.69033f, 0.69257f, |
656 | 0.69482f, 0.69709f, 0.69939f, 0.70169f, |
657 | 0.70402f, 0.70634f, 0.70869f, 0.71107f, |
658 | 0.71346f, 0.71587f, 0.71829f, 0.72073f, |
659 | 0.72320f, 0.72567f, 0.72818f, 0.73069f, |
660 | 0.73323f, 0.73579f, 0.73838f, 0.74098f, |
661 | 0.74360f, 0.74622f, 0.74890f, 0.75159f, |
662 | 0.75429f, 0.75704f, 0.75979f, 0.76257f, |
663 | 0.76537f, 0.76821f, 0.77109f, 0.77396f, |
664 | 0.77688f, 0.77982f, 0.78278f, 0.78579f, |
665 | 0.78883f, 0.79187f, 0.79498f, 0.79809f, |
666 | 0.80127f, 0.80445f, 0.80767f, 0.81095f, |
667 | 0.81424f, 0.81757f, 0.82094f, 0.82438f, |
668 | 0.82782f, 0.83133f, 0.83488f, 0.83847f, |
669 | 0.84210f, 0.84577f, 0.84951f, 0.85328f, |
670 | 0.85713f, 0.86103f, 0.86499f, 0.86900f, |
671 | 0.87306f, 0.87720f, 0.88139f, 0.88566f, |
672 | 0.89000f, 0.89442f, 0.89891f, 0.90350f, |
673 | 0.90818f, 0.91295f, 0.91780f, 0.92272f, |
674 | 0.92780f, 0.93299f, 0.93828f, 0.94369f, |
675 | 0.94926f, 0.95493f, 0.96082f, 0.96684f, |
676 | 0.97305f, 0.97943f, 0.98605f, 0.99291f, |
677 | 1.00000f |
678 | }; |
679 | |
680 | const uint32 kTableSize = sizeof (kTable ) / |
681 | sizeof (kTable [0]); |
682 | |
683 | real32 y = (real32) x * (real32) (kTableSize - 1); |
684 | |
685 | int32 index = Pin_int32 (0, (int32) y, kTableSize - 2); |
686 | |
687 | real32 fract = y - (real32) index; |
688 | |
689 | return kTable [index ] * (1.0f - fract) + |
690 | kTable [index + 1] * ( fract); |
691 | |
692 | } |
693 | |
694 | /*****************************************************************************/ |
695 | |
696 | const dng_1d_function & dng_tone_curve_acr3_default::Get () |
697 | { |
698 | |
699 | static dng_tone_curve_acr3_default static_dng_tone_curve_acr3_default; |
700 | |
701 | return static_dng_tone_curve_acr3_default; |
702 | |
703 | } |
704 | |
705 | /*****************************************************************************/ |
706 | |
707 | class dng_render_task: public dng_filter_task |
708 | { |
709 | |
710 | protected: |
711 | |
712 | const dng_negative &fNegative; |
713 | |
714 | const dng_render &fParams; |
715 | |
716 | dng_point fSrcOffset; |
717 | |
718 | dng_vector fCameraWhite; |
719 | dng_matrix fCameraToRGB; |
720 | |
721 | AutoPtr<dng_hue_sat_map> fHueSatMap; |
722 | |
723 | dng_1d_table fExposureRamp; |
724 | |
725 | AutoPtr<dng_hue_sat_map> fLookTable; |
726 | |
727 | dng_1d_table fToneCurve; |
728 | |
729 | dng_matrix fRGBtoFinal; |
730 | |
731 | dng_1d_table fEncodeGamma; |
732 | |
733 | AutoPtr<dng_1d_table> fHueSatMapEncode; |
734 | AutoPtr<dng_1d_table> fHueSatMapDecode; |
735 | |
736 | AutoPtr<dng_1d_table> fLookTableEncode; |
737 | AutoPtr<dng_1d_table> fLookTableDecode; |
738 | |
739 | AutoPtr<dng_memory_block> fTempBuffer [kMaxMPThreads]; |
740 | |
741 | public: |
742 | |
743 | dng_render_task (const dng_image &srcImage, |
744 | dng_image &dstImage, |
745 | const dng_negative &negative, |
746 | const dng_render ¶ms, |
747 | const dng_point &srcOffset); |
748 | |
749 | virtual dng_rect SrcArea (const dng_rect &dstArea); |
750 | |
751 | virtual void Start (uint32 threadCount, |
752 | const dng_point &tileSize, |
753 | dng_memory_allocator *allocator, |
754 | dng_abort_sniffer *sniffer); |
755 | |
756 | virtual void ProcessArea (uint32 threadIndex, |
757 | dng_pixel_buffer &srcBuffer, |
758 | dng_pixel_buffer &dstBuffer); |
759 | |
760 | }; |
761 | |
762 | /*****************************************************************************/ |
763 | |
764 | dng_render_task::dng_render_task (const dng_image &srcImage, |
765 | dng_image &dstImage, |
766 | const dng_negative &negative, |
767 | const dng_render ¶ms, |
768 | const dng_point &srcOffset) |
769 | |
770 | : dng_filter_task (srcImage, |
771 | dstImage) |
772 | |
773 | , fNegative (negative ) |
774 | , fParams (params ) |
775 | , fSrcOffset (srcOffset) |
776 | |
777 | , fCameraWhite () |
778 | , fCameraToRGB () |
779 | |
780 | , fHueSatMap () |
781 | |
782 | , fExposureRamp () |
783 | |
784 | , fLookTable () |
785 | |
786 | , fToneCurve () |
787 | |
788 | , fRGBtoFinal () |
789 | |
790 | , fEncodeGamma () |
791 | |
792 | , fHueSatMapEncode () |
793 | , fHueSatMapDecode () |
794 | |
795 | , fLookTableEncode () |
796 | , fLookTableDecode () |
797 | |
798 | { |
799 | |
800 | fSrcPixelType = ttFloat; |
801 | fDstPixelType = ttFloat; |
802 | |
803 | } |
804 | |
805 | /*****************************************************************************/ |
806 | |
807 | dng_rect dng_render_task::SrcArea (const dng_rect &dstArea) |
808 | { |
809 | |
810 | return dstArea + fSrcOffset; |
811 | |
812 | } |
813 | |
814 | /*****************************************************************************/ |
815 | |
816 | void dng_render_task::Start (uint32 threadCount, |
817 | const dng_point &tileSize, |
818 | dng_memory_allocator *allocator, |
819 | dng_abort_sniffer *sniffer) |
820 | { |
821 | |
822 | dng_filter_task::Start (threadCount, |
823 | tileSize, |
824 | allocator, |
825 | sniffer); |
826 | |
827 | // Compute camera space to linear ProPhoto RGB parameters. |
828 | |
829 | dng_camera_profile_id profileID; // Default profile ID. |
830 | |
831 | if (!fNegative.IsMonochrome ()) |
832 | { |
833 | |
834 | AutoPtr<dng_color_spec> spec (fNegative.MakeColorSpec (profileID)); |
835 | |
836 | if (fParams.WhiteXY ().IsValid ()) |
837 | { |
838 | |
839 | spec->SetWhiteXY (fParams.WhiteXY ()); |
840 | |
841 | } |
842 | |
843 | else if (fNegative.HasCameraNeutral ()) |
844 | { |
845 | |
846 | spec->SetWhiteXY (spec->NeutralToXY (fNegative.CameraNeutral ())); |
847 | |
848 | } |
849 | |
850 | else if (fNegative.HasCameraWhiteXY ()) |
851 | { |
852 | |
853 | spec->SetWhiteXY (fNegative.CameraWhiteXY ()); |
854 | |
855 | } |
856 | |
857 | else |
858 | { |
859 | |
860 | spec->SetWhiteXY (D55_xy_coord ()); |
861 | |
862 | } |
863 | |
864 | fCameraWhite = spec->CameraWhite (); |
865 | |
866 | fCameraToRGB = dng_space_ProPhoto::Get ().MatrixFromPCS () * |
867 | spec->CameraToPCS (); |
868 | |
869 | // Find Hue/Sat table, if any. |
870 | |
871 | const dng_camera_profile *profile = fNegative.ProfileByID (profileID); |
872 | |
873 | if (profile) |
874 | { |
875 | |
876 | fHueSatMap.Reset (profile->HueSatMapForWhite (spec->WhiteXY ())); |
877 | |
878 | if (profile->HasLookTable ()) |
879 | { |
880 | |
881 | fLookTable.Reset (new dng_hue_sat_map (profile->LookTable ())); |
882 | |
883 | } |
884 | |
885 | if (profile->HueSatMapEncoding () != encoding_Linear) |
886 | { |
887 | |
888 | BuildHueSatMapEncodingTable (*allocator, |
889 | profile->HueSatMapEncoding (), |
890 | fHueSatMapEncode, |
891 | fHueSatMapDecode, |
892 | false); |
893 | |
894 | } |
895 | |
896 | if (profile->LookTableEncoding () != encoding_Linear) |
897 | { |
898 | |
899 | BuildHueSatMapEncodingTable (*allocator, |
900 | profile->LookTableEncoding (), |
901 | fLookTableEncode, |
902 | fLookTableDecode, |
903 | false); |
904 | |
905 | } |
906 | |
907 | } |
908 | |
909 | } |
910 | |
911 | // Compute exposure/shadows ramp. |
912 | |
913 | real64 exposure = fParams.Exposure () + |
914 | fNegative.TotalBaselineExposure (profileID) - |
915 | (log (fNegative.Stage3Gain ()) / log (2.0)); |
916 | |
917 | { |
918 | |
919 | real64 white = 1.0 / pow (2.0, Max_real64 (0.0, exposure)); |
920 | |
921 | real64 black = fParams.Shadows () * |
922 | fNegative.ShadowScale () * |
923 | fNegative.Stage3Gain () * |
924 | 0.001; |
925 | |
926 | black = Min_real64 (black, 0.99 * white); |
927 | |
928 | dng_function_exposure_ramp rampFunction (white, |
929 | black, |
930 | black); |
931 | |
932 | fExposureRamp.Initialize (*allocator, rampFunction); |
933 | |
934 | } |
935 | |
936 | // Compute tone curve. |
937 | |
938 | { |
939 | |
940 | // If there is any negative exposure compenation to perform |
941 | // (beyond what the camera provides for with its baseline exposure), |
942 | // we fake this by darkening the tone curve. |
943 | |
944 | dng_function_exposure_tone exposureTone (exposure); |
945 | |
946 | dng_1d_concatenate totalTone (exposureTone, |
947 | fParams.ToneCurve ()); |
948 | |
949 | fToneCurve.Initialize (*allocator, totalTone); |
950 | |
951 | } |
952 | |
953 | // Compute linear ProPhoto RGB to final space parameters. |
954 | |
955 | { |
956 | |
957 | const dng_color_space &finalSpace = fParams.FinalSpace (); |
958 | |
959 | fRGBtoFinal = finalSpace.MatrixFromPCS () * |
960 | dng_space_ProPhoto::Get ().MatrixToPCS (); |
961 | |
962 | fEncodeGamma.Initialize (*allocator, finalSpace.GammaFunction ()); |
963 | |
964 | } |
965 | |
966 | // Allocate temp buffer to hold one row of RGB data. |
967 | |
968 | uint32 tempBufferSize = 0; |
969 | |
970 | if (!SafeUint32Mult (tileSize.h, (uint32) sizeof (real32), &tempBufferSize) || |
971 | !SafeUint32Mult (tempBufferSize, 3, &tempBufferSize)) |
972 | { |
973 | |
974 | ThrowMemoryFull("Arithmetic overflow computing buffer size." ); |
975 | |
976 | } |
977 | |
978 | for (uint32 threadIndex = 0; threadIndex < threadCount; threadIndex++) |
979 | { |
980 | |
981 | fTempBuffer [threadIndex] . Reset (allocator->Allocate (tempBufferSize)); |
982 | |
983 | } |
984 | |
985 | } |
986 | |
987 | /*****************************************************************************/ |
988 | |
989 | void dng_render_task::ProcessArea (uint32 threadIndex, |
990 | dng_pixel_buffer &srcBuffer, |
991 | dng_pixel_buffer &dstBuffer) |
992 | { |
993 | |
994 | dng_rect srcArea = srcBuffer.fArea; |
995 | dng_rect dstArea = dstBuffer.fArea; |
996 | |
997 | uint32 srcCols = srcArea.W (); |
998 | |
999 | real32 *tPtrR = fTempBuffer [threadIndex]->Buffer_real32 (); |
1000 | |
1001 | real32 *tPtrG = tPtrR + srcCols; |
1002 | real32 *tPtrB = tPtrG + srcCols; |
1003 | |
1004 | for (int32 srcRow = srcArea.t; srcRow < srcArea.b; srcRow++) |
1005 | { |
1006 | |
1007 | // First convert from camera native space to linear PhotoRGB, |
1008 | // applying the white balance and camera profile. |
1009 | |
1010 | { |
1011 | |
1012 | const real32 *sPtrA = (const real32 *) |
1013 | srcBuffer.ConstPixel (srcRow, |
1014 | srcArea.l, |
1015 | 0); |
1016 | |
1017 | if (fSrcPlanes == 1) |
1018 | { |
1019 | |
1020 | // For monochrome cameras, this just requires copying |
1021 | // the data into all three color channels. |
1022 | |
1023 | DoCopyBytes (sPtrA, tPtrR, srcCols * (uint32) sizeof (real32)); |
1024 | DoCopyBytes (sPtrA, tPtrG, srcCols * (uint32) sizeof (real32)); |
1025 | DoCopyBytes (sPtrA, tPtrB, srcCols * (uint32) sizeof (real32)); |
1026 | |
1027 | } |
1028 | |
1029 | else |
1030 | { |
1031 | |
1032 | const real32 *sPtrB = sPtrA + srcBuffer.fPlaneStep; |
1033 | const real32 *sPtrC = sPtrB + srcBuffer.fPlaneStep; |
1034 | |
1035 | if (fSrcPlanes == 3) |
1036 | { |
1037 | |
1038 | DoBaselineABCtoRGB (sPtrA, |
1039 | sPtrB, |
1040 | sPtrC, |
1041 | tPtrR, |
1042 | tPtrG, |
1043 | tPtrB, |
1044 | srcCols, |
1045 | fCameraWhite, |
1046 | fCameraToRGB); |
1047 | |
1048 | } |
1049 | |
1050 | else |
1051 | { |
1052 | |
1053 | const real32 *sPtrD = sPtrC + srcBuffer.fPlaneStep; |
1054 | |
1055 | DoBaselineABCDtoRGB (sPtrA, |
1056 | sPtrB, |
1057 | sPtrC, |
1058 | sPtrD, |
1059 | tPtrR, |
1060 | tPtrG, |
1061 | tPtrB, |
1062 | srcCols, |
1063 | fCameraWhite, |
1064 | fCameraToRGB); |
1065 | |
1066 | } |
1067 | |
1068 | // Apply Hue/Sat map, if any. |
1069 | |
1070 | if (fHueSatMap.Get ()) |
1071 | { |
1072 | |
1073 | DoBaselineHueSatMap (tPtrR, |
1074 | tPtrG, |
1075 | tPtrB, |
1076 | tPtrR, |
1077 | tPtrG, |
1078 | tPtrB, |
1079 | srcCols, |
1080 | *fHueSatMap.Get (), |
1081 | fHueSatMapEncode.Get (), |
1082 | fHueSatMapDecode.Get ()); |
1083 | |
1084 | } |
1085 | |
1086 | } |
1087 | |
1088 | } |
1089 | |
1090 | // Apply exposure curve. |
1091 | |
1092 | DoBaseline1DTable (tPtrR, |
1093 | tPtrR, |
1094 | srcCols, |
1095 | fExposureRamp); |
1096 | |
1097 | DoBaseline1DTable (tPtrG, |
1098 | tPtrG, |
1099 | srcCols, |
1100 | fExposureRamp); |
1101 | |
1102 | DoBaseline1DTable (tPtrB, |
1103 | tPtrB, |
1104 | srcCols, |
1105 | fExposureRamp); |
1106 | |
1107 | // Apply look table, if any. |
1108 | |
1109 | if (fLookTable.Get ()) |
1110 | { |
1111 | |
1112 | DoBaselineHueSatMap (tPtrR, |
1113 | tPtrG, |
1114 | tPtrB, |
1115 | tPtrR, |
1116 | tPtrG, |
1117 | tPtrB, |
1118 | srcCols, |
1119 | *fLookTable.Get (), |
1120 | fLookTableEncode.Get (), |
1121 | fLookTableDecode.Get ()); |
1122 | |
1123 | } |
1124 | |
1125 | // Apply baseline tone curve. |
1126 | |
1127 | DoBaselineRGBTone (tPtrR, |
1128 | tPtrG, |
1129 | tPtrB, |
1130 | tPtrR, |
1131 | tPtrG, |
1132 | tPtrB, |
1133 | srcCols, |
1134 | fToneCurve); |
1135 | |
1136 | // Convert to final color space. |
1137 | |
1138 | int32 dstRow = srcRow + (dstArea.t - srcArea.t); |
1139 | |
1140 | if (fDstPlanes == 1) |
1141 | { |
1142 | |
1143 | real32 *dPtrG = dstBuffer.DirtyPixel_real32 (dstRow, |
1144 | dstArea.l, |
1145 | 0); |
1146 | |
1147 | DoBaselineRGBtoGray (tPtrR, |
1148 | tPtrG, |
1149 | tPtrB, |
1150 | dPtrG, |
1151 | srcCols, |
1152 | fRGBtoFinal); |
1153 | |
1154 | DoBaseline1DTable (dPtrG, |
1155 | dPtrG, |
1156 | srcCols, |
1157 | fEncodeGamma); |
1158 | |
1159 | } |
1160 | |
1161 | else |
1162 | { |
1163 | |
1164 | real32 *dPtrR = dstBuffer.DirtyPixel_real32 (dstRow, |
1165 | dstArea.l, |
1166 | 0); |
1167 | |
1168 | real32 *dPtrG = dPtrR + dstBuffer.fPlaneStep; |
1169 | real32 *dPtrB = dPtrG + dstBuffer.fPlaneStep; |
1170 | |
1171 | DoBaselineRGBtoRGB (tPtrR, |
1172 | tPtrG, |
1173 | tPtrB, |
1174 | dPtrR, |
1175 | dPtrG, |
1176 | dPtrB, |
1177 | srcCols, |
1178 | fRGBtoFinal); |
1179 | |
1180 | DoBaseline1DTable (dPtrR, |
1181 | dPtrR, |
1182 | srcCols, |
1183 | fEncodeGamma); |
1184 | |
1185 | DoBaseline1DTable (dPtrG, |
1186 | dPtrG, |
1187 | srcCols, |
1188 | fEncodeGamma); |
1189 | |
1190 | DoBaseline1DTable (dPtrB, |
1191 | dPtrB, |
1192 | srcCols, |
1193 | fEncodeGamma); |
1194 | |
1195 | } |
1196 | |
1197 | } |
1198 | |
1199 | } |
1200 | |
1201 | /*****************************************************************************/ |
1202 | |
1203 | dng_render::dng_render (dng_host &host, |
1204 | const dng_negative &negative) |
1205 | |
1206 | : fHost (host) |
1207 | , fNegative (negative) |
1208 | |
1209 | , fWhiteXY () |
1210 | |
1211 | , fExposure (0.0) |
1212 | , fShadows (5.0) |
1213 | |
1214 | , fToneCurve (&dng_tone_curve_acr3_default::Get ()) |
1215 | |
1216 | , fFinalSpace (&dng_space_sRGB::Get ()) |
1217 | , fFinalPixelType (ttByte) |
1218 | |
1219 | , fMaximumSize (0) |
1220 | |
1221 | , fProfileToneCurve () |
1222 | |
1223 | { |
1224 | |
1225 | // Switch to NOP default parameters for non-scene referred data. |
1226 | |
1227 | if (fNegative.ColorimetricReference () != crSceneReferred) |
1228 | { |
1229 | |
1230 | fShadows = 0.0; |
1231 | |
1232 | fToneCurve = &dng_1d_identity::Get (); |
1233 | |
1234 | } |
1235 | |
1236 | // Use default tone curve from profile if any. |
1237 | |
1238 | const dng_camera_profile *profile = fNegative.ProfileByID (dng_camera_profile_id ()); |
1239 | |
1240 | if (profile && profile->ToneCurve ().IsValid ()) |
1241 | { |
1242 | |
1243 | fProfileToneCurve.Reset (new dng_spline_solver); |
1244 | |
1245 | profile->ToneCurve ().Solve (*fProfileToneCurve.Get ()); |
1246 | |
1247 | fToneCurve = fProfileToneCurve.Get (); |
1248 | |
1249 | } |
1250 | |
1251 | // Turn off default shadow mapping if requested by profile. |
1252 | |
1253 | if (profile && (profile->DefaultBlackRender () == defaultBlackRender_None)) |
1254 | { |
1255 | |
1256 | fShadows = 0.0; |
1257 | |
1258 | } |
1259 | |
1260 | } |
1261 | |
1262 | /*****************************************************************************/ |
1263 | |
1264 | dng_image * dng_render::Render () |
1265 | { |
1266 | |
1267 | const dng_image *srcImage = fNegative.Stage3Image (); |
1268 | |
1269 | dng_rect srcBounds = fNegative.DefaultCropArea (); |
1270 | |
1271 | dng_point dstSize; |
1272 | |
1273 | dstSize.h = fNegative.DefaultFinalWidth (); |
1274 | dstSize.v = fNegative.DefaultFinalHeight (); |
1275 | |
1276 | if (MaximumSize ()) |
1277 | { |
1278 | |
1279 | if (Max_uint32 (dstSize.h, dstSize.v) > MaximumSize ()) |
1280 | { |
1281 | |
1282 | real64 ratio = fNegative.AspectRatio (); |
1283 | |
1284 | if (ratio >= 1.0) |
1285 | { |
1286 | dstSize.h = MaximumSize (); |
1287 | dstSize.v = Max_uint32 (1, Round_uint32 (dstSize.h / ratio)); |
1288 | } |
1289 | |
1290 | else |
1291 | { |
1292 | dstSize.v = MaximumSize (); |
1293 | dstSize.h = Max_uint32 (1, Round_uint32 (dstSize.v * ratio)); |
1294 | } |
1295 | |
1296 | } |
1297 | |
1298 | } |
1299 | |
1300 | AutoPtr<dng_image> tempImage; |
1301 | |
1302 | if (srcBounds.Size () != dstSize) |
1303 | { |
1304 | |
1305 | tempImage.Reset (fHost.Make_dng_image (dstSize, |
1306 | srcImage->Planes (), |
1307 | srcImage->PixelType ())); |
1308 | |
1309 | ResampleImage (fHost, |
1310 | *srcImage, |
1311 | *tempImage.Get (), |
1312 | srcBounds, |
1313 | tempImage->Bounds (), |
1314 | dng_resample_bicubic::Get ()); |
1315 | |
1316 | srcImage = tempImage.Get (); |
1317 | |
1318 | srcBounds = tempImage->Bounds (); |
1319 | |
1320 | } |
1321 | |
1322 | uint32 dstPlanes = FinalSpace ().IsMonochrome () ? 1 : 3; |
1323 | |
1324 | AutoPtr<dng_image> dstImage (fHost.Make_dng_image (srcBounds.Size (), |
1325 | dstPlanes, |
1326 | FinalPixelType ())); |
1327 | |
1328 | dng_render_task task (*srcImage, |
1329 | *dstImage.Get (), |
1330 | fNegative, |
1331 | *this, |
1332 | srcBounds.TL ()); |
1333 | |
1334 | fHost.PerformAreaTask (task, |
1335 | dstImage->Bounds ()); |
1336 | |
1337 | return dstImage.Release (); |
1338 | |
1339 | } |
1340 | |
1341 | /*****************************************************************************/ |
1342 | |