1
2/* pngrtran.c - transforms the data in a row for PNG readers
3 *
4 * Last changed in libpng 1.6.22 [(PENDING RELEASE)]
5 * Copyright (c) 1998-2002,2004,2006-2016 Glenn Randers-Pehrson
6 * (Version 0.96 Copyright (c) 1996, 1997 Andreas Dilger)
7 * (Version 0.88 Copyright (c) 1995, 1996 Guy Eric Schalnat, Group 42, Inc.)
8 *
9 * This code is released under the libpng license.
10 * For conditions of distribution and use, see the disclaimer
11 * and license in png.h
12 *
13 * This file contains functions optionally called by an application
14 * in order to tell libpng how to handle data when reading a PNG.
15 * Transformations that are used in both reading and writing are
16 * in pngtrans.c.
17 */
18
19#include "pngpriv.h"
20
21#ifdef PNG_READ_SUPPORTED
22
23/* Set the action on getting a CRC error for an ancillary or critical chunk. */
24void PNGAPI
25png_set_crc_action(png_structrp png_ptr, int crit_action, int ancil_action)
26{
27 png_debug(1, "in png_set_crc_action");
28
29 if (png_ptr == NULL)
30 return;
31
32 /* Tell libpng how we react to CRC errors in critical chunks */
33 switch (crit_action)
34 {
35 case PNG_CRC_NO_CHANGE: /* Leave setting as is */
36 break;
37
38 case PNG_CRC_WARN_USE: /* Warn/use data */
39 png_ptr->flags &= ~PNG_FLAG_CRC_CRITICAL_MASK;
40 png_ptr->flags |= PNG_FLAG_CRC_CRITICAL_USE;
41 break;
42
43 case PNG_CRC_QUIET_USE: /* Quiet/use data */
44 png_ptr->flags &= ~PNG_FLAG_CRC_CRITICAL_MASK;
45 png_ptr->flags |= PNG_FLAG_CRC_CRITICAL_USE |
46 PNG_FLAG_CRC_CRITICAL_IGNORE;
47 break;
48
49 case PNG_CRC_WARN_DISCARD: /* Not a valid action for critical data */
50 png_warning(png_ptr,
51 "Can't discard critical data on CRC error");
52 case PNG_CRC_ERROR_QUIT: /* Error/quit */
53
54 case PNG_CRC_DEFAULT:
55 default:
56 png_ptr->flags &= ~PNG_FLAG_CRC_CRITICAL_MASK;
57 break;
58 }
59
60 /* Tell libpng how we react to CRC errors in ancillary chunks */
61 switch (ancil_action)
62 {
63 case PNG_CRC_NO_CHANGE: /* Leave setting as is */
64 break;
65
66 case PNG_CRC_WARN_USE: /* Warn/use data */
67 png_ptr->flags &= ~PNG_FLAG_CRC_ANCILLARY_MASK;
68 png_ptr->flags |= PNG_FLAG_CRC_ANCILLARY_USE;
69 break;
70
71 case PNG_CRC_QUIET_USE: /* Quiet/use data */
72 png_ptr->flags &= ~PNG_FLAG_CRC_ANCILLARY_MASK;
73 png_ptr->flags |= PNG_FLAG_CRC_ANCILLARY_USE |
74 PNG_FLAG_CRC_ANCILLARY_NOWARN;
75 break;
76
77 case PNG_CRC_ERROR_QUIT: /* Error/quit */
78 png_ptr->flags &= ~PNG_FLAG_CRC_ANCILLARY_MASK;
79 png_ptr->flags |= PNG_FLAG_CRC_ANCILLARY_NOWARN;
80 break;
81
82 case PNG_CRC_WARN_DISCARD: /* Warn/discard data */
83
84 case PNG_CRC_DEFAULT:
85 default:
86 png_ptr->flags &= ~PNG_FLAG_CRC_ANCILLARY_MASK;
87 break;
88 }
89}
90
91#ifdef PNG_READ_TRANSFORMS_SUPPORTED
92/* Is it OK to set a transformation now? Only if png_start_read_image or
93 * png_read_update_info have not been called. It is not necessary for the IHDR
94 * to have been read in all cases; the need_IHDR parameter allows for this
95 * check too.
96 */
97static int
98png_rtran_ok(png_structrp png_ptr, int need_IHDR)
99{
100 if (png_ptr != NULL)
101 {
102 if ((png_ptr->flags & PNG_FLAG_ROW_INIT) != 0)
103 png_app_error(png_ptr,
104 "invalid after png_start_read_image or png_read_update_info");
105
106 else if (need_IHDR && (png_ptr->mode & PNG_HAVE_IHDR) == 0)
107 png_app_error(png_ptr, "invalid before the PNG header has been read");
108
109 else
110 {
111 /* Turn on failure to initialize correctly for all transforms. */
112 png_ptr->flags |= PNG_FLAG_DETECT_UNINITIALIZED;
113
114 return 1; /* Ok */
115 }
116 }
117
118 return 0; /* no png_error possible! */
119}
120#endif
121
122#ifdef PNG_READ_BACKGROUND_SUPPORTED
123/* Handle alpha and tRNS via a background color */
124void PNGFAPI
125png_set_background_fixed(png_structrp png_ptr,
126 png_const_color_16p background_color, int background_gamma_code,
127 int need_expand, png_fixed_point background_gamma)
128{
129 png_debug(1, "in png_set_background_fixed");
130
131 if (png_rtran_ok(png_ptr, 0) == 0 || background_color == NULL)
132 return;
133
134 if (background_gamma_code == PNG_BACKGROUND_GAMMA_UNKNOWN)
135 {
136 png_warning(png_ptr, "Application must supply a known background gamma");
137 return;
138 }
139
140 png_ptr->transformations |= PNG_COMPOSE | PNG_STRIP_ALPHA;
141 png_ptr->transformations &= ~PNG_ENCODE_ALPHA;
142 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
143
144 png_ptr->background = *background_color;
145 png_ptr->background_gamma = background_gamma;
146 png_ptr->background_gamma_type = (png_byte)(background_gamma_code);
147 if (need_expand != 0)
148 png_ptr->transformations |= PNG_BACKGROUND_EXPAND;
149 else
150 png_ptr->transformations &= ~PNG_BACKGROUND_EXPAND;
151}
152
153# ifdef PNG_FLOATING_POINT_SUPPORTED
154void PNGAPI
155png_set_background(png_structrp png_ptr,
156 png_const_color_16p background_color, int background_gamma_code,
157 int need_expand, double background_gamma)
158{
159 png_set_background_fixed(png_ptr, background_color, background_gamma_code,
160 need_expand, png_fixed(png_ptr, background_gamma, "png_set_background"));
161}
162# endif /* FLOATING_POINT */
163#endif /* READ_BACKGROUND */
164
165/* Scale 16-bit depth files to 8-bit depth. If both of these are set then the
166 * one that pngrtran does first (scale) happens. This is necessary to allow the
167 * TRANSFORM and API behavior to be somewhat consistent, and it's simpler.
168 */
169#ifdef PNG_READ_SCALE_16_TO_8_SUPPORTED
170void PNGAPI
171png_set_scale_16(png_structrp png_ptr)
172{
173 png_debug(1, "in png_set_scale_16");
174
175 if (png_rtran_ok(png_ptr, 0) == 0)
176 return;
177
178 png_ptr->transformations |= PNG_SCALE_16_TO_8;
179}
180#endif
181
182#ifdef PNG_READ_STRIP_16_TO_8_SUPPORTED
183/* Chop 16-bit depth files to 8-bit depth */
184void PNGAPI
185png_set_strip_16(png_structrp png_ptr)
186{
187 png_debug(1, "in png_set_strip_16");
188
189 if (png_rtran_ok(png_ptr, 0) == 0)
190 return;
191
192 png_ptr->transformations |= PNG_16_TO_8;
193}
194#endif
195
196#ifdef PNG_READ_STRIP_ALPHA_SUPPORTED
197void PNGAPI
198png_set_strip_alpha(png_structrp png_ptr)
199{
200 png_debug(1, "in png_set_strip_alpha");
201
202 if (png_rtran_ok(png_ptr, 0) == 0)
203 return;
204
205 png_ptr->transformations |= PNG_STRIP_ALPHA;
206}
207#endif
208
209#if defined(PNG_READ_ALPHA_MODE_SUPPORTED) || defined(PNG_READ_GAMMA_SUPPORTED)
210static png_fixed_point
211translate_gamma_flags(png_structrp png_ptr, png_fixed_point output_gamma,
212 int is_screen)
213{
214 /* Check for flag values. The main reason for having the old Mac value as a
215 * flag is that it is pretty near impossible to work out what the correct
216 * value is from Apple documentation - a working Mac system is needed to
217 * discover the value!
218 */
219 if (output_gamma == PNG_DEFAULT_sRGB ||
220 output_gamma == PNG_FP_1 / PNG_DEFAULT_sRGB)
221 {
222 /* If there is no sRGB support this just sets the gamma to the standard
223 * sRGB value. (This is a side effect of using this function!)
224 */
225# ifdef PNG_READ_sRGB_SUPPORTED
226 png_ptr->flags |= PNG_FLAG_ASSUME_sRGB;
227# else
228 PNG_UNUSED(png_ptr)
229# endif
230 if (is_screen != 0)
231 output_gamma = PNG_GAMMA_sRGB;
232 else
233 output_gamma = PNG_GAMMA_sRGB_INVERSE;
234 }
235
236 else if (output_gamma == PNG_GAMMA_MAC_18 ||
237 output_gamma == PNG_FP_1 / PNG_GAMMA_MAC_18)
238 {
239 if (is_screen != 0)
240 output_gamma = PNG_GAMMA_MAC_OLD;
241 else
242 output_gamma = PNG_GAMMA_MAC_INVERSE;
243 }
244
245 return output_gamma;
246}
247
248# ifdef PNG_FLOATING_POINT_SUPPORTED
249static png_fixed_point
250convert_gamma_value(png_structrp png_ptr, double output_gamma)
251{
252 /* The following silently ignores cases where fixed point (times 100,000)
253 * gamma values are passed to the floating point API. This is safe and it
254 * means the fixed point constants work just fine with the floating point
255 * API. The alternative would just lead to undetected errors and spurious
256 * bug reports. Negative values fail inside the _fixed API unless they
257 * correspond to the flag values.
258 */
259 if (output_gamma > 0 && output_gamma < 128)
260 output_gamma *= PNG_FP_1;
261
262 /* This preserves -1 and -2 exactly: */
263 output_gamma = floor(output_gamma + .5);
264
265 if (output_gamma > PNG_FP_MAX || output_gamma < PNG_FP_MIN)
266 png_fixed_error(png_ptr, "gamma value");
267
268 return (png_fixed_point)output_gamma;
269}
270# endif
271#endif /* READ_ALPHA_MODE || READ_GAMMA */
272
273#ifdef PNG_READ_ALPHA_MODE_SUPPORTED
274void PNGFAPI
275png_set_alpha_mode_fixed(png_structrp png_ptr, int mode,
276 png_fixed_point output_gamma)
277{
278 int compose = 0;
279 png_fixed_point file_gamma;
280
281 png_debug(1, "in png_set_alpha_mode");
282
283 if (png_rtran_ok(png_ptr, 0) == 0)
284 return;
285
286 output_gamma = translate_gamma_flags(png_ptr, output_gamma, 1/*screen*/);
287
288 /* Validate the value to ensure it is in a reasonable range. The value
289 * is expected to be 1 or greater, but this range test allows for some
290 * viewing correction values. The intent is to weed out users of this API
291 * who use the inverse of the gamma value accidentally! Since some of these
292 * values are reasonable this may have to be changed:
293 *
294 * 1.6.x: changed from 0.07..3 to 0.01..100 (to accomodate the optimal 16-bit
295 * gamma of 36, and its reciprocal.)
296 */
297 if (output_gamma < 1000 || output_gamma > 10000000)
298 png_error(png_ptr, "output gamma out of expected range");
299
300 /* The default file gamma is the inverse of the output gamma; the output
301 * gamma may be changed below so get the file value first:
302 */
303 file_gamma = png_reciprocal(output_gamma);
304
305 /* There are really 8 possibilities here, composed of any combination
306 * of:
307 *
308 * premultiply the color channels
309 * do not encode non-opaque pixels
310 * encode the alpha as well as the color channels
311 *
312 * The differences disappear if the input/output ('screen') gamma is 1.0,
313 * because then the encoding is a no-op and there is only the choice of
314 * premultiplying the color channels or not.
315 *
316 * png_set_alpha_mode and png_set_background interact because both use
317 * png_compose to do the work. Calling both is only useful when
318 * png_set_alpha_mode is used to set the default mode - PNG_ALPHA_PNG - along
319 * with a default gamma value. Otherwise PNG_COMPOSE must not be set.
320 */
321 switch (mode)
322 {
323 case PNG_ALPHA_PNG: /* default: png standard */
324 /* No compose, but it may be set by png_set_background! */
325 png_ptr->transformations &= ~PNG_ENCODE_ALPHA;
326 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
327 break;
328
329 case PNG_ALPHA_ASSOCIATED: /* color channels premultiplied */
330 compose = 1;
331 png_ptr->transformations &= ~PNG_ENCODE_ALPHA;
332 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
333 /* The output is linear: */
334 output_gamma = PNG_FP_1;
335 break;
336
337 case PNG_ALPHA_OPTIMIZED: /* associated, non-opaque pixels linear */
338 compose = 1;
339 png_ptr->transformations &= ~PNG_ENCODE_ALPHA;
340 png_ptr->flags |= PNG_FLAG_OPTIMIZE_ALPHA;
341 /* output_gamma records the encoding of opaque pixels! */
342 break;
343
344 case PNG_ALPHA_BROKEN: /* associated, non-linear, alpha encoded */
345 compose = 1;
346 png_ptr->transformations |= PNG_ENCODE_ALPHA;
347 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
348 break;
349
350 default:
351 png_error(png_ptr, "invalid alpha mode");
352 }
353
354 /* Only set the default gamma if the file gamma has not been set (this has
355 * the side effect that the gamma in a second call to png_set_alpha_mode will
356 * be ignored.)
357 */
358 if (png_ptr->colorspace.gamma == 0)
359 {
360 png_ptr->colorspace.gamma = file_gamma;
361 png_ptr->colorspace.flags |= PNG_COLORSPACE_HAVE_GAMMA;
362 }
363
364 /* But always set the output gamma: */
365 png_ptr->screen_gamma = output_gamma;
366
367 /* Finally, if pre-multiplying, set the background fields to achieve the
368 * desired result.
369 */
370 if (compose != 0)
371 {
372 /* And obtain alpha pre-multiplication by composing on black: */
373 memset(&png_ptr->background, 0, (sizeof png_ptr->background));
374 png_ptr->background_gamma = png_ptr->colorspace.gamma; /* just in case */
375 png_ptr->background_gamma_type = PNG_BACKGROUND_GAMMA_FILE;
376 png_ptr->transformations &= ~PNG_BACKGROUND_EXPAND;
377
378 if ((png_ptr->transformations & PNG_COMPOSE) != 0)
379 png_error(png_ptr,
380 "conflicting calls to set alpha mode and background");
381
382 png_ptr->transformations |= PNG_COMPOSE;
383 }
384}
385
386# ifdef PNG_FLOATING_POINT_SUPPORTED
387void PNGAPI
388png_set_alpha_mode(png_structrp png_ptr, int mode, double output_gamma)
389{
390 png_set_alpha_mode_fixed(png_ptr, mode, convert_gamma_value(png_ptr,
391 output_gamma));
392}
393# endif
394#endif
395
396#ifdef PNG_READ_QUANTIZE_SUPPORTED
397/* Dither file to 8-bit. Supply a palette, the current number
398 * of elements in the palette, the maximum number of elements
399 * allowed, and a histogram if possible. If the current number
400 * of colors is greater than the maximum number, the palette will be
401 * modified to fit in the maximum number. "full_quantize" indicates
402 * whether we need a quantizing cube set up for RGB images, or if we
403 * simply are reducing the number of colors in a paletted image.
404 */
405
406typedef struct png_dsort_struct
407{
408 struct png_dsort_struct * next;
409 png_byte left;
410 png_byte right;
411} png_dsort;
412typedef png_dsort * png_dsortp;
413typedef png_dsort * * png_dsortpp;
414
415void PNGAPI
416png_set_quantize(png_structrp png_ptr, png_colorp palette,
417 int num_palette, int maximum_colors, png_const_uint_16p histogram,
418 int full_quantize)
419{
420 png_debug(1, "in png_set_quantize");
421
422 if (png_rtran_ok(png_ptr, 0) == 0)
423 return;
424
425 png_ptr->transformations |= PNG_QUANTIZE;
426
427 if (full_quantize == 0)
428 {
429 int i;
430
431 png_ptr->quantize_index = (png_bytep)png_malloc(png_ptr,
432 (png_uint_32)(num_palette * (sizeof (png_byte))));
433 for (i = 0; i < num_palette; i++)
434 png_ptr->quantize_index[i] = (png_byte)i;
435 }
436
437 if (num_palette > maximum_colors)
438 {
439 if (histogram != NULL)
440 {
441 /* This is easy enough, just throw out the least used colors.
442 * Perhaps not the best solution, but good enough.
443 */
444
445 int i;
446
447 /* Initialize an array to sort colors */
448 png_ptr->quantize_sort = (png_bytep)png_malloc(png_ptr,
449 (png_uint_32)(num_palette * (sizeof (png_byte))));
450
451 /* Initialize the quantize_sort array */
452 for (i = 0; i < num_palette; i++)
453 png_ptr->quantize_sort[i] = (png_byte)i;
454
455 /* Find the least used palette entries by starting a
456 * bubble sort, and running it until we have sorted
457 * out enough colors. Note that we don't care about
458 * sorting all the colors, just finding which are
459 * least used.
460 */
461
462 for (i = num_palette - 1; i >= maximum_colors; i--)
463 {
464 int done; /* To stop early if the list is pre-sorted */
465 int j;
466
467 done = 1;
468 for (j = 0; j < i; j++)
469 {
470 if (histogram[png_ptr->quantize_sort[j]]
471 < histogram[png_ptr->quantize_sort[j + 1]])
472 {
473 png_byte t;
474
475 t = png_ptr->quantize_sort[j];
476 png_ptr->quantize_sort[j] = png_ptr->quantize_sort[j + 1];
477 png_ptr->quantize_sort[j + 1] = t;
478 done = 0;
479 }
480 }
481
482 if (done != 0)
483 break;
484 }
485
486 /* Swap the palette around, and set up a table, if necessary */
487 if (full_quantize != 0)
488 {
489 int j = num_palette;
490
491 /* Put all the useful colors within the max, but don't
492 * move the others.
493 */
494 for (i = 0; i < maximum_colors; i++)
495 {
496 if ((int)png_ptr->quantize_sort[i] >= maximum_colors)
497 {
498 do
499 j--;
500 while ((int)png_ptr->quantize_sort[j] >= maximum_colors);
501
502 palette[i] = palette[j];
503 }
504 }
505 }
506 else
507 {
508 int j = num_palette;
509
510 /* Move all the used colors inside the max limit, and
511 * develop a translation table.
512 */
513 for (i = 0; i < maximum_colors; i++)
514 {
515 /* Only move the colors we need to */
516 if ((int)png_ptr->quantize_sort[i] >= maximum_colors)
517 {
518 png_color tmp_color;
519
520 do
521 j--;
522 while ((int)png_ptr->quantize_sort[j] >= maximum_colors);
523
524 tmp_color = palette[j];
525 palette[j] = palette[i];
526 palette[i] = tmp_color;
527 /* Indicate where the color went */
528 png_ptr->quantize_index[j] = (png_byte)i;
529 png_ptr->quantize_index[i] = (png_byte)j;
530 }
531 }
532
533 /* Find closest color for those colors we are not using */
534 for (i = 0; i < num_palette; i++)
535 {
536 if ((int)png_ptr->quantize_index[i] >= maximum_colors)
537 {
538 int min_d, k, min_k, d_index;
539
540 /* Find the closest color to one we threw out */
541 d_index = png_ptr->quantize_index[i];
542 min_d = PNG_COLOR_DIST(palette[d_index], palette[0]);
543 for (k = 1, min_k = 0; k < maximum_colors; k++)
544 {
545 int d;
546
547 d = PNG_COLOR_DIST(palette[d_index], palette[k]);
548
549 if (d < min_d)
550 {
551 min_d = d;
552 min_k = k;
553 }
554 }
555 /* Point to closest color */
556 png_ptr->quantize_index[i] = (png_byte)min_k;
557 }
558 }
559 }
560 png_free(png_ptr, png_ptr->quantize_sort);
561 png_ptr->quantize_sort = NULL;
562 }
563 else
564 {
565 /* This is much harder to do simply (and quickly). Perhaps
566 * we need to go through a median cut routine, but those
567 * don't always behave themselves with only a few colors
568 * as input. So we will just find the closest two colors,
569 * and throw out one of them (chosen somewhat randomly).
570 * [We don't understand this at all, so if someone wants to
571 * work on improving it, be our guest - AED, GRP]
572 */
573 int i;
574 int max_d;
575 int num_new_palette;
576 png_dsortp t;
577 png_dsortpp hash;
578
579 t = NULL;
580
581 /* Initialize palette index arrays */
582 png_ptr->index_to_palette = (png_bytep)png_malloc(png_ptr,
583 (png_uint_32)(num_palette * (sizeof (png_byte))));
584 png_ptr->palette_to_index = (png_bytep)png_malloc(png_ptr,
585 (png_uint_32)(num_palette * (sizeof (png_byte))));
586
587 /* Initialize the sort array */
588 for (i = 0; i < num_palette; i++)
589 {
590 png_ptr->index_to_palette[i] = (png_byte)i;
591 png_ptr->palette_to_index[i] = (png_byte)i;
592 }
593
594 hash = (png_dsortpp)png_calloc(png_ptr, (png_uint_32)(769 *
595 (sizeof (png_dsortp))));
596
597 num_new_palette = num_palette;
598
599 /* Initial wild guess at how far apart the farthest pixel
600 * pair we will be eliminating will be. Larger
601 * numbers mean more areas will be allocated, Smaller
602 * numbers run the risk of not saving enough data, and
603 * having to do this all over again.
604 *
605 * I have not done extensive checking on this number.
606 */
607 max_d = 96;
608
609 while (num_new_palette > maximum_colors)
610 {
611 for (i = 0; i < num_new_palette - 1; i++)
612 {
613 int j;
614
615 for (j = i + 1; j < num_new_palette; j++)
616 {
617 int d;
618
619 d = PNG_COLOR_DIST(palette[i], palette[j]);
620
621 if (d <= max_d)
622 {
623
624 t = (png_dsortp)png_malloc_warn(png_ptr,
625 (png_uint_32)(sizeof (png_dsort)));
626
627 if (t == NULL)
628 break;
629
630 t->next = hash[d];
631 t->left = (png_byte)i;
632 t->right = (png_byte)j;
633 hash[d] = t;
634 }
635 }
636 if (t == NULL)
637 break;
638 }
639
640 if (t != NULL)
641 for (i = 0; i <= max_d; i++)
642 {
643 if (hash[i] != NULL)
644 {
645 png_dsortp p;
646
647 for (p = hash[i]; p; p = p->next)
648 {
649 if ((int)png_ptr->index_to_palette[p->left]
650 < num_new_palette &&
651 (int)png_ptr->index_to_palette[p->right]
652 < num_new_palette)
653 {
654 int j, next_j;
655
656 if (num_new_palette & 0x01)
657 {
658 j = p->left;
659 next_j = p->right;
660 }
661 else
662 {
663 j = p->right;
664 next_j = p->left;
665 }
666
667 num_new_palette--;
668 palette[png_ptr->index_to_palette[j]]
669 = palette[num_new_palette];
670 if (full_quantize == 0)
671 {
672 int k;
673
674 for (k = 0; k < num_palette; k++)
675 {
676 if (png_ptr->quantize_index[k] ==
677 png_ptr->index_to_palette[j])
678 png_ptr->quantize_index[k] =
679 png_ptr->index_to_palette[next_j];
680
681 if ((int)png_ptr->quantize_index[k] ==
682 num_new_palette)
683 png_ptr->quantize_index[k] =
684 png_ptr->index_to_palette[j];
685 }
686 }
687
688 png_ptr->index_to_palette[png_ptr->palette_to_index
689 [num_new_palette]] = png_ptr->index_to_palette[j];
690
691 png_ptr->palette_to_index[png_ptr->index_to_palette[j]]
692 = png_ptr->palette_to_index[num_new_palette];
693
694 png_ptr->index_to_palette[j] =
695 (png_byte)num_new_palette;
696
697 png_ptr->palette_to_index[num_new_palette] =
698 (png_byte)j;
699 }
700 if (num_new_palette <= maximum_colors)
701 break;
702 }
703 if (num_new_palette <= maximum_colors)
704 break;
705 }
706 }
707
708 for (i = 0; i < 769; i++)
709 {
710 if (hash[i] != NULL)
711 {
712 png_dsortp p = hash[i];
713 while (p)
714 {
715 t = p->next;
716 png_free(png_ptr, p);
717 p = t;
718 }
719 }
720 hash[i] = 0;
721 }
722 max_d += 96;
723 }
724 png_free(png_ptr, hash);
725 png_free(png_ptr, png_ptr->palette_to_index);
726 png_free(png_ptr, png_ptr->index_to_palette);
727 png_ptr->palette_to_index = NULL;
728 png_ptr->index_to_palette = NULL;
729 }
730 num_palette = maximum_colors;
731 }
732 if (png_ptr->palette == NULL)
733 {
734 png_ptr->palette = palette;
735 }
736 png_ptr->num_palette = (png_uint_16)num_palette;
737
738 if (full_quantize != 0)
739 {
740 int i;
741 png_bytep distance;
742 int total_bits = PNG_QUANTIZE_RED_BITS + PNG_QUANTIZE_GREEN_BITS +
743 PNG_QUANTIZE_BLUE_BITS;
744 int num_red = (1 << PNG_QUANTIZE_RED_BITS);
745 int num_green = (1 << PNG_QUANTIZE_GREEN_BITS);
746 int num_blue = (1 << PNG_QUANTIZE_BLUE_BITS);
747 png_size_t num_entries = ((png_size_t)1 << total_bits);
748
749 png_ptr->palette_lookup = (png_bytep)png_calloc(png_ptr,
750 (png_uint_32)(num_entries * (sizeof (png_byte))));
751
752 distance = (png_bytep)png_malloc(png_ptr, (png_uint_32)(num_entries *
753 (sizeof (png_byte))));
754
755 memset(distance, 0xff, num_entries * (sizeof (png_byte)));
756
757 for (i = 0; i < num_palette; i++)
758 {
759 int ir, ig, ib;
760 int r = (palette[i].red >> (8 - PNG_QUANTIZE_RED_BITS));
761 int g = (palette[i].green >> (8 - PNG_QUANTIZE_GREEN_BITS));
762 int b = (palette[i].blue >> (8 - PNG_QUANTIZE_BLUE_BITS));
763
764 for (ir = 0; ir < num_red; ir++)
765 {
766 /* int dr = abs(ir - r); */
767 int dr = ((ir > r) ? ir - r : r - ir);
768 int index_r = (ir << (PNG_QUANTIZE_BLUE_BITS +
769 PNG_QUANTIZE_GREEN_BITS));
770
771 for (ig = 0; ig < num_green; ig++)
772 {
773 /* int dg = abs(ig - g); */
774 int dg = ((ig > g) ? ig - g : g - ig);
775 int dt = dr + dg;
776 int dm = ((dr > dg) ? dr : dg);
777 int index_g = index_r | (ig << PNG_QUANTIZE_BLUE_BITS);
778
779 for (ib = 0; ib < num_blue; ib++)
780 {
781 int d_index = index_g | ib;
782 /* int db = abs(ib - b); */
783 int db = ((ib > b) ? ib - b : b - ib);
784 int dmax = ((dm > db) ? dm : db);
785 int d = dmax + dt + db;
786
787 if (d < (int)distance[d_index])
788 {
789 distance[d_index] = (png_byte)d;
790 png_ptr->palette_lookup[d_index] = (png_byte)i;
791 }
792 }
793 }
794 }
795 }
796
797 png_free(png_ptr, distance);
798 }
799}
800#endif /* READ_QUANTIZE */
801
802#ifdef PNG_READ_GAMMA_SUPPORTED
803void PNGFAPI
804png_set_gamma_fixed(png_structrp png_ptr, png_fixed_point scrn_gamma,
805 png_fixed_point file_gamma)
806{
807 png_debug(1, "in png_set_gamma_fixed");
808
809 if (png_rtran_ok(png_ptr, 0) == 0)
810 return;
811
812 /* New in libpng-1.5.4 - reserve particular negative values as flags. */
813 scrn_gamma = translate_gamma_flags(png_ptr, scrn_gamma, 1/*screen*/);
814 file_gamma = translate_gamma_flags(png_ptr, file_gamma, 0/*file*/);
815
816 /* Checking the gamma values for being >0 was added in 1.5.4 along with the
817 * premultiplied alpha support; this actually hides an undocumented feature
818 * of the previous implementation which allowed gamma processing to be
819 * disabled in background handling. There is no evidence (so far) that this
820 * was being used; however, png_set_background itself accepted and must still
821 * accept '0' for the gamma value it takes, because it isn't always used.
822 *
823 * Since this is an API change (albeit a very minor one that removes an
824 * undocumented API feature) the following checks were only enabled in
825 * libpng-1.6.0.
826 */
827 if (file_gamma <= 0)
828 png_error(png_ptr, "invalid file gamma in png_set_gamma");
829
830 if (scrn_gamma <= 0)
831 png_error(png_ptr, "invalid screen gamma in png_set_gamma");
832
833 /* Set the gamma values unconditionally - this overrides the value in the PNG
834 * file if a gAMA chunk was present. png_set_alpha_mode provides a
835 * different, easier, way to default the file gamma.
836 */
837 png_ptr->colorspace.gamma = file_gamma;
838 png_ptr->colorspace.flags |= PNG_COLORSPACE_HAVE_GAMMA;
839 png_ptr->screen_gamma = scrn_gamma;
840}
841
842# ifdef PNG_FLOATING_POINT_SUPPORTED
843void PNGAPI
844png_set_gamma(png_structrp png_ptr, double scrn_gamma, double file_gamma)
845{
846 png_set_gamma_fixed(png_ptr, convert_gamma_value(png_ptr, scrn_gamma),
847 convert_gamma_value(png_ptr, file_gamma));
848}
849# endif /* FLOATING_POINT */
850#endif /* READ_GAMMA */
851
852#ifdef PNG_READ_EXPAND_SUPPORTED
853/* Expand paletted images to RGB, expand grayscale images of
854 * less than 8-bit depth to 8-bit depth, and expand tRNS chunks
855 * to alpha channels.
856 */
857void PNGAPI
858png_set_expand(png_structrp png_ptr)
859{
860 png_debug(1, "in png_set_expand");
861
862 if (png_rtran_ok(png_ptr, 0) == 0)
863 return;
864
865 png_ptr->transformations |= (PNG_EXPAND | PNG_EXPAND_tRNS);
866}
867
868/* GRR 19990627: the following three functions currently are identical
869 * to png_set_expand(). However, it is entirely reasonable that someone
870 * might wish to expand an indexed image to RGB but *not* expand a single,
871 * fully transparent palette entry to a full alpha channel--perhaps instead
872 * convert tRNS to the grayscale/RGB format (16-bit RGB value), or replace
873 * the transparent color with a particular RGB value, or drop tRNS entirely.
874 * IOW, a future version of the library may make the transformations flag
875 * a bit more fine-grained, with separate bits for each of these three
876 * functions.
877 *
878 * More to the point, these functions make it obvious what libpng will be
879 * doing, whereas "expand" can (and does) mean any number of things.
880 *
881 * GRP 20060307: In libpng-1.2.9, png_set_gray_1_2_4_to_8() was modified
882 * to expand only the sample depth but not to expand the tRNS to alpha
883 * and its name was changed to png_set_expand_gray_1_2_4_to_8().
884 */
885
886/* Expand paletted images to RGB. */
887void PNGAPI
888png_set_palette_to_rgb(png_structrp png_ptr)
889{
890 png_debug(1, "in png_set_palette_to_rgb");
891
892 if (png_rtran_ok(png_ptr, 0) == 0)
893 return;
894
895 png_ptr->transformations |= (PNG_EXPAND | PNG_EXPAND_tRNS);
896}
897
898/* Expand grayscale images of less than 8-bit depth to 8 bits. */
899void PNGAPI
900png_set_expand_gray_1_2_4_to_8(png_structrp png_ptr)
901{
902 png_debug(1, "in png_set_expand_gray_1_2_4_to_8");
903
904 if (png_rtran_ok(png_ptr, 0) == 0)
905 return;
906
907 png_ptr->transformations |= PNG_EXPAND;
908}
909
910/* Expand tRNS chunks to alpha channels. */
911void PNGAPI
912png_set_tRNS_to_alpha(png_structrp png_ptr)
913{
914 png_debug(1, "in png_set_tRNS_to_alpha");
915
916 if (png_rtran_ok(png_ptr, 0) == 0)
917 return;
918
919 png_ptr->transformations |= (PNG_EXPAND | PNG_EXPAND_tRNS);
920}
921#endif /* READ_EXPAND */
922
923#ifdef PNG_READ_EXPAND_16_SUPPORTED
924/* Expand to 16-bit channels, expand the tRNS chunk too (because otherwise
925 * it may not work correctly.)
926 */
927void PNGAPI
928png_set_expand_16(png_structrp png_ptr)
929{
930 png_debug(1, "in png_set_expand_16");
931
932 if (png_rtran_ok(png_ptr, 0) == 0)
933 return;
934
935 png_ptr->transformations |= (PNG_EXPAND_16 | PNG_EXPAND | PNG_EXPAND_tRNS);
936}
937#endif
938
939#ifdef PNG_READ_GRAY_TO_RGB_SUPPORTED
940void PNGAPI
941png_set_gray_to_rgb(png_structrp png_ptr)
942{
943 png_debug(1, "in png_set_gray_to_rgb");
944
945 if (png_rtran_ok(png_ptr, 0) == 0)
946 return;
947
948 /* Because rgb must be 8 bits or more: */
949 png_set_expand_gray_1_2_4_to_8(png_ptr);
950 png_ptr->transformations |= PNG_GRAY_TO_RGB;
951}
952#endif
953
954#ifdef PNG_READ_RGB_TO_GRAY_SUPPORTED
955void PNGFAPI
956png_set_rgb_to_gray_fixed(png_structrp png_ptr, int error_action,
957 png_fixed_point red, png_fixed_point green)
958{
959 png_debug(1, "in png_set_rgb_to_gray");
960
961 /* Need the IHDR here because of the check on color_type below. */
962 /* TODO: fix this */
963 if (png_rtran_ok(png_ptr, 1) == 0)
964 return;
965
966 switch (error_action)
967 {
968 case PNG_ERROR_ACTION_NONE:
969 png_ptr->transformations |= PNG_RGB_TO_GRAY;
970 break;
971
972 case PNG_ERROR_ACTION_WARN:
973 png_ptr->transformations |= PNG_RGB_TO_GRAY_WARN;
974 break;
975
976 case PNG_ERROR_ACTION_ERROR:
977 png_ptr->transformations |= PNG_RGB_TO_GRAY_ERR;
978 break;
979
980 default:
981 png_error(png_ptr, "invalid error action to rgb_to_gray");
982 }
983
984 if (png_ptr->color_type == PNG_COLOR_TYPE_PALETTE)
985#ifdef PNG_READ_EXPAND_SUPPORTED
986 png_ptr->transformations |= PNG_EXPAND;
987#else
988 {
989 /* Make this an error in 1.6 because otherwise the application may assume
990 * that it just worked and get a memory overwrite.
991 */
992 png_error(png_ptr,
993 "Cannot do RGB_TO_GRAY without EXPAND_SUPPORTED");
994
995 /* png_ptr->transformations &= ~PNG_RGB_TO_GRAY; */
996 }
997#endif
998 {
999 if (red >= 0 && green >= 0 && red + green <= PNG_FP_1)
1000 {
1001 png_uint_16 red_int, green_int;
1002
1003 /* NOTE: this calculation does not round, but this behavior is retained
1004 * for consistency; the inaccuracy is very small. The code here always
1005 * overwrites the coefficients, regardless of whether they have been
1006 * defaulted or set already.
1007 */
1008 red_int = (png_uint_16)(((png_uint_32)red*32768)/100000);
1009 green_int = (png_uint_16)(((png_uint_32)green*32768)/100000);
1010
1011 png_ptr->rgb_to_gray_red_coeff = red_int;
1012 png_ptr->rgb_to_gray_green_coeff = green_int;
1013 png_ptr->rgb_to_gray_coefficients_set = 1;
1014 }
1015
1016 else
1017 {
1018 if (red >= 0 && green >= 0)
1019 png_app_warning(png_ptr,
1020 "ignoring out of range rgb_to_gray coefficients");
1021
1022 /* Use the defaults, from the cHRM chunk if set, else the historical
1023 * values which are close to the sRGB/HDTV/ITU-Rec 709 values. See
1024 * png_do_rgb_to_gray for more discussion of the values. In this case
1025 * the coefficients are not marked as 'set' and are not overwritten if
1026 * something has already provided a default.
1027 */
1028 if (png_ptr->rgb_to_gray_red_coeff == 0 &&
1029 png_ptr->rgb_to_gray_green_coeff == 0)
1030 {
1031 png_ptr->rgb_to_gray_red_coeff = 6968;
1032 png_ptr->rgb_to_gray_green_coeff = 23434;
1033 /* png_ptr->rgb_to_gray_blue_coeff = 2366; */
1034 }
1035 }
1036 }
1037}
1038
1039#ifdef PNG_FLOATING_POINT_SUPPORTED
1040/* Convert a RGB image to a grayscale of the same width. This allows us,
1041 * for example, to convert a 24 bpp RGB image into an 8 bpp grayscale image.
1042 */
1043
1044void PNGAPI
1045png_set_rgb_to_gray(png_structrp png_ptr, int error_action, double red,
1046 double green)
1047{
1048 png_set_rgb_to_gray_fixed(png_ptr, error_action,
1049 png_fixed(png_ptr, red, "rgb to gray red coefficient"),
1050 png_fixed(png_ptr, green, "rgb to gray green coefficient"));
1051}
1052#endif /* FLOATING POINT */
1053
1054#endif /* RGB_TO_GRAY */
1055
1056#if defined(PNG_READ_USER_TRANSFORM_SUPPORTED) || \
1057 defined(PNG_WRITE_USER_TRANSFORM_SUPPORTED)
1058void PNGAPI
1059png_set_read_user_transform_fn(png_structrp png_ptr, png_user_transform_ptr
1060 read_user_transform_fn)
1061{
1062 png_debug(1, "in png_set_read_user_transform_fn");
1063
1064#ifdef PNG_READ_USER_TRANSFORM_SUPPORTED
1065 png_ptr->transformations |= PNG_USER_TRANSFORM;
1066 png_ptr->read_user_transform_fn = read_user_transform_fn;
1067#endif
1068}
1069#endif
1070
1071#ifdef PNG_READ_TRANSFORMS_SUPPORTED
1072#ifdef PNG_READ_GAMMA_SUPPORTED
1073/* In the case of gamma transformations only do transformations on images where
1074 * the [file] gamma and screen_gamma are not close reciprocals, otherwise it
1075 * slows things down slightly, and also needlessly introduces small errors.
1076 */
1077static int /* PRIVATE */
1078png_gamma_threshold(png_fixed_point screen_gamma, png_fixed_point file_gamma)
1079{
1080 /* PNG_GAMMA_THRESHOLD is the threshold for performing gamma
1081 * correction as a difference of the overall transform from 1.0
1082 *
1083 * We want to compare the threshold with s*f - 1, if we get
1084 * overflow here it is because of wacky gamma values so we
1085 * turn on processing anyway.
1086 */
1087 png_fixed_point gtest;
1088 return !png_muldiv(&gtest, screen_gamma, file_gamma, PNG_FP_1) ||
1089 png_gamma_significant(gtest);
1090}
1091#endif
1092
1093/* Initialize everything needed for the read. This includes modifying
1094 * the palette.
1095 */
1096
1097/* For the moment 'png_init_palette_transformations' and
1098 * 'png_init_rgb_transformations' only do some flag canceling optimizations.
1099 * The intent is that these two routines should have palette or rgb operations
1100 * extracted from 'png_init_read_transformations'.
1101 */
1102static void /* PRIVATE */
1103png_init_palette_transformations(png_structrp png_ptr)
1104{
1105 /* Called to handle the (input) palette case. In png_do_read_transformations
1106 * the first step is to expand the palette if requested, so this code must
1107 * take care to only make changes that are invariant with respect to the
1108 * palette expansion, or only do them if there is no expansion.
1109 *
1110 * STRIP_ALPHA has already been handled in the caller (by setting num_trans
1111 * to 0.)
1112 */
1113 int input_has_alpha = 0;
1114 int input_has_transparency = 0;
1115
1116 if (png_ptr->num_trans > 0)
1117 {
1118 int i;
1119
1120 /* Ignore if all the entries are opaque (unlikely!) */
1121 for (i=0; i<png_ptr->num_trans; ++i)
1122 {
1123 if (png_ptr->trans_alpha[i] == 255)
1124 continue;
1125 else if (png_ptr->trans_alpha[i] == 0)
1126 input_has_transparency = 1;
1127 else
1128 {
1129 input_has_transparency = 1;
1130 input_has_alpha = 1;
1131 break;
1132 }
1133 }
1134 }
1135
1136 /* If no alpha we can optimize. */
1137 if (input_has_alpha == 0)
1138 {
1139 /* Any alpha means background and associative alpha processing is
1140 * required, however if the alpha is 0 or 1 throughout OPTIMIZE_ALPHA
1141 * and ENCODE_ALPHA are irrelevant.
1142 */
1143 png_ptr->transformations &= ~PNG_ENCODE_ALPHA;
1144 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
1145
1146 if (input_has_transparency == 0)
1147 png_ptr->transformations &= ~(PNG_COMPOSE | PNG_BACKGROUND_EXPAND);
1148 }
1149
1150#if defined(PNG_READ_EXPAND_SUPPORTED) && defined(PNG_READ_BACKGROUND_SUPPORTED)
1151 /* png_set_background handling - deals with the complexity of whether the
1152 * background color is in the file format or the screen format in the case
1153 * where an 'expand' will happen.
1154 */
1155
1156 /* The following code cannot be entered in the alpha pre-multiplication case
1157 * because PNG_BACKGROUND_EXPAND is cancelled below.
1158 */
1159 if ((png_ptr->transformations & PNG_BACKGROUND_EXPAND) != 0 &&
1160 (png_ptr->transformations & PNG_EXPAND) != 0)
1161 {
1162 {
1163 png_ptr->background.red =
1164 png_ptr->palette[png_ptr->background.index].red;
1165 png_ptr->background.green =
1166 png_ptr->palette[png_ptr->background.index].green;
1167 png_ptr->background.blue =
1168 png_ptr->palette[png_ptr->background.index].blue;
1169
1170#ifdef PNG_READ_INVERT_ALPHA_SUPPORTED
1171 if ((png_ptr->transformations & PNG_INVERT_ALPHA) != 0)
1172 {
1173 if ((png_ptr->transformations & PNG_EXPAND_tRNS) == 0)
1174 {
1175 /* Invert the alpha channel (in tRNS) unless the pixels are
1176 * going to be expanded, in which case leave it for later
1177 */
1178 int i, istop = png_ptr->num_trans;
1179
1180 for (i=0; i<istop; i++)
1181 png_ptr->trans_alpha[i] = (png_byte)(255 -
1182 png_ptr->trans_alpha[i]);
1183 }
1184 }
1185#endif /* READ_INVERT_ALPHA */
1186 }
1187 } /* background expand and (therefore) no alpha association. */
1188#endif /* READ_EXPAND && READ_BACKGROUND */
1189}
1190
1191static void /* PRIVATE */
1192png_init_rgb_transformations(png_structrp png_ptr)
1193{
1194 /* Added to libpng-1.5.4: check the color type to determine whether there
1195 * is any alpha or transparency in the image and simply cancel the
1196 * background and alpha mode stuff if there isn't.
1197 */
1198 int input_has_alpha = (png_ptr->color_type & PNG_COLOR_MASK_ALPHA) != 0;
1199 int input_has_transparency = png_ptr->num_trans > 0;
1200
1201 /* If no alpha we can optimize. */
1202 if (input_has_alpha == 0)
1203 {
1204 /* Any alpha means background and associative alpha processing is
1205 * required, however if the alpha is 0 or 1 throughout OPTIMIZE_ALPHA
1206 * and ENCODE_ALPHA are irrelevant.
1207 */
1208# ifdef PNG_READ_ALPHA_MODE_SUPPORTED
1209 png_ptr->transformations &= ~PNG_ENCODE_ALPHA;
1210 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
1211# endif
1212
1213 if (input_has_transparency == 0)
1214 png_ptr->transformations &= ~(PNG_COMPOSE | PNG_BACKGROUND_EXPAND);
1215 }
1216
1217#if defined(PNG_READ_EXPAND_SUPPORTED) && defined(PNG_READ_BACKGROUND_SUPPORTED)
1218 /* png_set_background handling - deals with the complexity of whether the
1219 * background color is in the file format or the screen format in the case
1220 * where an 'expand' will happen.
1221 */
1222
1223 /* The following code cannot be entered in the alpha pre-multiplication case
1224 * because PNG_BACKGROUND_EXPAND is cancelled below.
1225 */
1226 if ((png_ptr->transformations & PNG_BACKGROUND_EXPAND) != 0 &&
1227 (png_ptr->transformations & PNG_EXPAND) != 0 &&
1228 (png_ptr->color_type & PNG_COLOR_MASK_COLOR) == 0)
1229 /* i.e., GRAY or GRAY_ALPHA */
1230 {
1231 {
1232 /* Expand background and tRNS chunks */
1233 int gray = png_ptr->background.gray;
1234 int trans_gray = png_ptr->trans_color.gray;
1235
1236 switch (png_ptr->bit_depth)
1237 {
1238 case 1:
1239 gray *= 0xff;
1240 trans_gray *= 0xff;
1241 break;
1242
1243 case 2:
1244 gray *= 0x55;
1245 trans_gray *= 0x55;
1246 break;
1247
1248 case 4:
1249 gray *= 0x11;
1250 trans_gray *= 0x11;
1251 break;
1252
1253 default:
1254
1255 case 8:
1256 /* FALL THROUGH (Already 8 bits) */
1257
1258 case 16:
1259 /* Already a full 16 bits */
1260 break;
1261 }
1262
1263 png_ptr->background.red = png_ptr->background.green =
1264 png_ptr->background.blue = (png_uint_16)gray;
1265
1266 if ((png_ptr->transformations & PNG_EXPAND_tRNS) == 0)
1267 {
1268 png_ptr->trans_color.red = png_ptr->trans_color.green =
1269 png_ptr->trans_color.blue = (png_uint_16)trans_gray;
1270 }
1271 }
1272 } /* background expand and (therefore) no alpha association. */
1273#endif /* READ_EXPAND && READ_BACKGROUND */
1274}
1275
1276void /* PRIVATE */
1277png_init_read_transformations(png_structrp png_ptr)
1278{
1279 png_debug(1, "in png_init_read_transformations");
1280
1281 /* This internal function is called from png_read_start_row in pngrutil.c
1282 * and it is called before the 'rowbytes' calculation is done, so the code
1283 * in here can change or update the transformations flags.
1284 *
1285 * First do updates that do not depend on the details of the PNG image data
1286 * being processed.
1287 */
1288
1289#ifdef PNG_READ_GAMMA_SUPPORTED
1290 /* Prior to 1.5.4 these tests were performed from png_set_gamma, 1.5.4 adds
1291 * png_set_alpha_mode and this is another source for a default file gamma so
1292 * the test needs to be performed later - here. In addition prior to 1.5.4
1293 * the tests were repeated for the PALETTE color type here - this is no
1294 * longer necessary (and doesn't seem to have been necessary before.)
1295 */
1296 {
1297 /* The following temporary indicates if overall gamma correction is
1298 * required.
1299 */
1300 int gamma_correction = 0;
1301
1302 if (png_ptr->colorspace.gamma != 0) /* has been set */
1303 {
1304 if (png_ptr->screen_gamma != 0) /* screen set too */
1305 gamma_correction = png_gamma_threshold(png_ptr->colorspace.gamma,
1306 png_ptr->screen_gamma);
1307
1308 else
1309 /* Assume the output matches the input; a long time default behavior
1310 * of libpng, although the standard has nothing to say about this.
1311 */
1312 png_ptr->screen_gamma = png_reciprocal(png_ptr->colorspace.gamma);
1313 }
1314
1315 else if (png_ptr->screen_gamma != 0)
1316 /* The converse - assume the file matches the screen, note that this
1317 * perhaps undesireable default can (from 1.5.4) be changed by calling
1318 * png_set_alpha_mode (even if the alpha handling mode isn't required
1319 * or isn't changed from the default.)
1320 */
1321 png_ptr->colorspace.gamma = png_reciprocal(png_ptr->screen_gamma);
1322
1323 else /* neither are set */
1324 /* Just in case the following prevents any processing - file and screen
1325 * are both assumed to be linear and there is no way to introduce a
1326 * third gamma value other than png_set_background with 'UNIQUE', and,
1327 * prior to 1.5.4
1328 */
1329 png_ptr->screen_gamma = png_ptr->colorspace.gamma = PNG_FP_1;
1330
1331 /* We have a gamma value now. */
1332 png_ptr->colorspace.flags |= PNG_COLORSPACE_HAVE_GAMMA;
1333
1334 /* Now turn the gamma transformation on or off as appropriate. Notice
1335 * that PNG_GAMMA just refers to the file->screen correction. Alpha
1336 * composition may independently cause gamma correction because it needs
1337 * linear data (e.g. if the file has a gAMA chunk but the screen gamma
1338 * hasn't been specified.) In any case this flag may get turned off in
1339 * the code immediately below if the transform can be handled outside the
1340 * row loop.
1341 */
1342 if (gamma_correction != 0)
1343 png_ptr->transformations |= PNG_GAMMA;
1344
1345 else
1346 png_ptr->transformations &= ~PNG_GAMMA;
1347 }
1348#endif
1349
1350 /* Certain transformations have the effect of preventing other
1351 * transformations that happen afterward in png_do_read_transformations;
1352 * resolve the interdependencies here. From the code of
1353 * png_do_read_transformations the order is:
1354 *
1355 * 1) PNG_EXPAND (including PNG_EXPAND_tRNS)
1356 * 2) PNG_STRIP_ALPHA (if no compose)
1357 * 3) PNG_RGB_TO_GRAY
1358 * 4) PNG_GRAY_TO_RGB iff !PNG_BACKGROUND_IS_GRAY
1359 * 5) PNG_COMPOSE
1360 * 6) PNG_GAMMA
1361 * 7) PNG_STRIP_ALPHA (if compose)
1362 * 8) PNG_ENCODE_ALPHA
1363 * 9) PNG_SCALE_16_TO_8
1364 * 10) PNG_16_TO_8
1365 * 11) PNG_QUANTIZE (converts to palette)
1366 * 12) PNG_EXPAND_16
1367 * 13) PNG_GRAY_TO_RGB iff PNG_BACKGROUND_IS_GRAY
1368 * 14) PNG_INVERT_MONO
1369 * 15) PNG_INVERT_ALPHA
1370 * 16) PNG_SHIFT
1371 * 17) PNG_PACK
1372 * 18) PNG_BGR
1373 * 19) PNG_PACKSWAP
1374 * 20) PNG_FILLER (includes PNG_ADD_ALPHA)
1375 * 21) PNG_SWAP_ALPHA
1376 * 22) PNG_SWAP_BYTES
1377 * 23) PNG_USER_TRANSFORM [must be last]
1378 */
1379#ifdef PNG_READ_STRIP_ALPHA_SUPPORTED
1380 if ((png_ptr->transformations & PNG_STRIP_ALPHA) != 0 &&
1381 (png_ptr->transformations & PNG_COMPOSE) == 0)
1382 {
1383 /* Stripping the alpha channel happens immediately after the 'expand'
1384 * transformations, before all other transformation, so it cancels out
1385 * the alpha handling. It has the side effect negating the effect of
1386 * PNG_EXPAND_tRNS too:
1387 */
1388 png_ptr->transformations &= ~(PNG_BACKGROUND_EXPAND | PNG_ENCODE_ALPHA |
1389 PNG_EXPAND_tRNS);
1390 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
1391
1392 /* Kill the tRNS chunk itself too. Prior to 1.5.4 this did not happen
1393 * so transparency information would remain just so long as it wasn't
1394 * expanded. This produces unexpected API changes if the set of things
1395 * that do PNG_EXPAND_tRNS changes (perfectly possible given the
1396 * documentation - which says ask for what you want, accept what you
1397 * get.) This makes the behavior consistent from 1.5.4:
1398 */
1399 png_ptr->num_trans = 0;
1400 }
1401#endif /* STRIP_ALPHA supported, no COMPOSE */
1402
1403#ifdef PNG_READ_ALPHA_MODE_SUPPORTED
1404 /* If the screen gamma is about 1.0 then the OPTIMIZE_ALPHA and ENCODE_ALPHA
1405 * settings will have no effect.
1406 */
1407 if (png_gamma_significant(png_ptr->screen_gamma) == 0)
1408 {
1409 png_ptr->transformations &= ~PNG_ENCODE_ALPHA;
1410 png_ptr->flags &= ~PNG_FLAG_OPTIMIZE_ALPHA;
1411 }
1412#endif
1413
1414#ifdef PNG_READ_RGB_TO_GRAY_SUPPORTED
1415 /* Make sure the coefficients for the rgb to gray conversion are set
1416 * appropriately.
1417 */
1418 if ((png_ptr->transformations & PNG_RGB_TO_GRAY) != 0)
1419 png_colorspace_set_rgb_coefficients(png_ptr);
1420#endif
1421
1422#ifdef PNG_READ_GRAY_TO_RGB_SUPPORTED
1423#if defined(PNG_READ_EXPAND_SUPPORTED) && defined(PNG_READ_BACKGROUND_SUPPORTED)
1424 /* Detect gray background and attempt to enable optimization for
1425 * gray --> RGB case.
1426 *
1427 * Note: if PNG_BACKGROUND_EXPAND is set and color_type is either RGB or
1428 * RGB_ALPHA (in which case need_expand is superfluous anyway), the
1429 * background color might actually be gray yet not be flagged as such.
1430 * This is not a problem for the current code, which uses
1431 * PNG_BACKGROUND_IS_GRAY only to decide when to do the
1432 * png_do_gray_to_rgb() transformation.
1433 *
1434 * TODO: this code needs to be revised to avoid the complexity and
1435 * interdependencies. The color type of the background should be recorded in
1436 * png_set_background, along with the bit depth, then the code has a record
1437 * of exactly what color space the background is currently in.
1438 */
1439 if ((png_ptr->transformations & PNG_BACKGROUND_EXPAND) != 0)
1440 {
1441 /* PNG_BACKGROUND_EXPAND: the background is in the file color space, so if
1442 * the file was grayscale the background value is gray.
1443 */
1444 if ((png_ptr->color_type & PNG_COLOR_MASK_COLOR) == 0)
1445 png_ptr->mode |= PNG_BACKGROUND_IS_GRAY;
1446 }
1447
1448 else if ((png_ptr->transformations & PNG_COMPOSE) != 0)
1449 {
1450 /* PNG_COMPOSE: png_set_background was called with need_expand false,
1451 * so the color is in the color space of the output or png_set_alpha_mode
1452 * was called and the color is black. Ignore RGB_TO_GRAY because that
1453 * happens before GRAY_TO_RGB.
1454 */
1455 if ((png_ptr->transformations & PNG_GRAY_TO_RGB) != 0)
1456 {
1457 if (png_ptr->background.red == png_ptr->background.green &&
1458 png_ptr->background.red == png_ptr->background.blue)
1459 {
1460 png_ptr->mode |= PNG_BACKGROUND_IS_GRAY;
1461 png_ptr->background.gray = png_ptr->background.red;
1462 }
1463 }
1464 }
1465#endif /* READ_EXPAND && READ_BACKGROUND */
1466#endif /* READ_GRAY_TO_RGB */
1467
1468 /* For indexed PNG data (PNG_COLOR_TYPE_PALETTE) many of the transformations
1469 * can be performed directly on the palette, and some (such as rgb to gray)
1470 * can be optimized inside the palette. This is particularly true of the
1471 * composite (background and alpha) stuff, which can be pretty much all done
1472 * in the palette even if the result is expanded to RGB or gray afterward.
1473 *
1474 * NOTE: this is Not Yet Implemented, the code behaves as in 1.5.1 and
1475 * earlier and the palette stuff is actually handled on the first row. This
1476 * leads to the reported bug that the palette returned by png_get_PLTE is not
1477 * updated.
1478 */
1479 if (png_ptr->color_type == PNG_COLOR_TYPE_PALETTE)
1480 png_init_palette_transformations(png_ptr);
1481
1482 else
1483 png_init_rgb_transformations(png_ptr);
1484
1485#if defined(PNG_READ_BACKGROUND_SUPPORTED) && \
1486 defined(PNG_READ_EXPAND_16_SUPPORTED)
1487 if ((png_ptr->transformations & PNG_EXPAND_16) != 0 &&
1488 (png_ptr->transformations & PNG_COMPOSE) != 0 &&
1489 (png_ptr->transformations & PNG_BACKGROUND_EXPAND) == 0 &&
1490 png_ptr->bit_depth != 16)
1491 {
1492 /* TODO: fix this. Because the expand_16 operation is after the compose
1493 * handling the background color must be 8, not 16, bits deep, but the
1494 * application will supply a 16-bit value so reduce it here.
1495 *
1496 * The PNG_BACKGROUND_EXPAND code above does not expand to 16 bits at
1497 * present, so that case is ok (until do_expand_16 is moved.)
1498 *
1499 * NOTE: this discards the low 16 bits of the user supplied background
1500 * color, but until expand_16 works properly there is no choice!
1501 */
1502# define CHOP(x) (x)=((png_uint_16)PNG_DIV257(x))
1503 CHOP(png_ptr->background.red);
1504 CHOP(png_ptr->background.green);
1505 CHOP(png_ptr->background.blue);
1506 CHOP(png_ptr->background.gray);
1507# undef CHOP
1508 }
1509#endif /* READ_BACKGROUND && READ_EXPAND_16 */
1510
1511#if defined(PNG_READ_BACKGROUND_SUPPORTED) && \
1512 (defined(PNG_READ_SCALE_16_TO_8_SUPPORTED) || \
1513 defined(PNG_READ_STRIP_16_TO_8_SUPPORTED))
1514 if ((png_ptr->transformations & (PNG_16_TO_8|PNG_SCALE_16_TO_8)) != 0 &&
1515 (png_ptr->transformations & PNG_COMPOSE) != 0 &&
1516 (png_ptr->transformations & PNG_BACKGROUND_EXPAND) == 0 &&
1517 png_ptr->bit_depth == 16)
1518 {
1519 /* On the other hand, if a 16-bit file is to be reduced to 8-bits per
1520 * component this will also happen after PNG_COMPOSE and so the background
1521 * color must be pre-expanded here.
1522 *
1523 * TODO: fix this too.
1524 */
1525 png_ptr->background.red = (png_uint_16)(png_ptr->background.red * 257);
1526 png_ptr->background.green =
1527 (png_uint_16)(png_ptr->background.green * 257);
1528 png_ptr->background.blue = (png_uint_16)(png_ptr->background.blue * 257);
1529 png_ptr->background.gray = (png_uint_16)(png_ptr->background.gray * 257);
1530 }
1531#endif
1532
1533 /* NOTE: below 'PNG_READ_ALPHA_MODE_SUPPORTED' is presumed to also enable the
1534 * background support (see the comments in scripts/pnglibconf.dfa), this
1535 * allows pre-multiplication of the alpha channel to be implemented as
1536 * compositing on black. This is probably sub-optimal and has been done in
1537 * 1.5.4 betas simply to enable external critique and testing (i.e. to
1538 * implement the new API quickly, without lots of internal changes.)
1539 */
1540
1541#ifdef PNG_READ_GAMMA_SUPPORTED
1542# ifdef PNG_READ_BACKGROUND_SUPPORTED
1543 /* Includes ALPHA_MODE */
1544 png_ptr->background_1 = png_ptr->background;
1545# endif
1546
1547 /* This needs to change - in the palette image case a whole set of tables are
1548 * built when it would be quicker to just calculate the correct value for
1549 * each palette entry directly. Also, the test is too tricky - why check
1550 * PNG_RGB_TO_GRAY if PNG_GAMMA is not set? The answer seems to be that
1551 * PNG_GAMMA is cancelled even if the gamma is known? The test excludes the
1552 * PNG_COMPOSE case, so apparently if there is no *overall* gamma correction
1553 * the gamma tables will not be built even if composition is required on a
1554 * gamma encoded value.
1555 *
1556 * In 1.5.4 this is addressed below by an additional check on the individual
1557 * file gamma - if it is not 1.0 both RGB_TO_GRAY and COMPOSE need the
1558 * tables.
1559 */
1560 if ((png_ptr->transformations & PNG_GAMMA) != 0 ||
1561 ((png_ptr->transformations & PNG_RGB_TO_GRAY) != 0 &&
1562 (png_gamma_significant(png_ptr->colorspace.gamma) != 0 ||
1563 png_gamma_significant(png_ptr->screen_gamma) != 0)) ||
1564 ((png_ptr->transformations & PNG_COMPOSE) != 0 &&
1565 (png_gamma_significant(png_ptr->colorspace.gamma) != 0 ||
1566 png_gamma_significant(png_ptr->screen_gamma) != 0
1567# ifdef PNG_READ_BACKGROUND_SUPPORTED
1568 || (png_ptr->background_gamma_type == PNG_BACKGROUND_GAMMA_UNIQUE &&
1569 png_gamma_significant(png_ptr->background_gamma) != 0)
1570# endif
1571 )) || ((png_ptr->transformations & PNG_ENCODE_ALPHA) != 0 &&
1572 png_gamma_significant(png_ptr->screen_gamma) != 0))
1573 {
1574 png_build_gamma_table(png_ptr, png_ptr->bit_depth);
1575
1576#ifdef PNG_READ_BACKGROUND_SUPPORTED
1577 if ((png_ptr->transformations & PNG_COMPOSE) != 0)
1578 {
1579 /* Issue a warning about this combination: because RGB_TO_GRAY is
1580 * optimized to do the gamma transform if present yet do_background has
1581 * to do the same thing if both options are set a
1582 * double-gamma-correction happens. This is true in all versions of
1583 * libpng to date.
1584 */
1585 if ((png_ptr->transformations & PNG_RGB_TO_GRAY) != 0)
1586 png_warning(png_ptr,
1587 "libpng does not support gamma+background+rgb_to_gray");
1588
1589 if ((png_ptr->color_type == PNG_COLOR_TYPE_PALETTE) != 0)
1590 {
1591 /* We don't get to here unless there is a tRNS chunk with non-opaque
1592 * entries - see the checking code at the start of this function.
1593 */
1594 png_color back, back_1;
1595 png_colorp palette = png_ptr->palette;
1596 int num_palette = png_ptr->num_palette;
1597 int i;
1598 if (png_ptr->background_gamma_type == PNG_BACKGROUND_GAMMA_FILE)
1599 {
1600
1601 back.red = png_ptr->gamma_table[png_ptr->background.red];
1602 back.green = png_ptr->gamma_table[png_ptr->background.green];
1603 back.blue = png_ptr->gamma_table[png_ptr->background.blue];
1604
1605 back_1.red = png_ptr->gamma_to_1[png_ptr->background.red];
1606 back_1.green = png_ptr->gamma_to_1[png_ptr->background.green];
1607 back_1.blue = png_ptr->gamma_to_1[png_ptr->background.blue];
1608 }
1609 else
1610 {
1611 png_fixed_point g, gs;
1612
1613 switch (png_ptr->background_gamma_type)
1614 {
1615 case PNG_BACKGROUND_GAMMA_SCREEN:
1616 g = (png_ptr->screen_gamma);
1617 gs = PNG_FP_1;
1618 break;
1619
1620 case PNG_BACKGROUND_GAMMA_FILE:
1621 g = png_reciprocal(png_ptr->colorspace.gamma);
1622 gs = png_reciprocal2(png_ptr->colorspace.gamma,
1623 png_ptr->screen_gamma);
1624 break;
1625
1626 case PNG_BACKGROUND_GAMMA_UNIQUE:
1627 g = png_reciprocal(png_ptr->background_gamma);
1628 gs = png_reciprocal2(png_ptr->background_gamma,
1629 png_ptr->screen_gamma);
1630 break;
1631 default:
1632 g = PNG_FP_1; /* back_1 */
1633 gs = PNG_FP_1; /* back */
1634 break;
1635 }
1636
1637 if (png_gamma_significant(gs) != 0)
1638 {
1639 back.red = png_gamma_8bit_correct(png_ptr->background.red,
1640 gs);
1641 back.green = png_gamma_8bit_correct(png_ptr->background.green,
1642 gs);
1643 back.blue = png_gamma_8bit_correct(png_ptr->background.blue,
1644 gs);
1645 }
1646
1647 else
1648 {
1649 back.red = (png_byte)png_ptr->background.red;
1650 back.green = (png_byte)png_ptr->background.green;
1651 back.blue = (png_byte)png_ptr->background.blue;
1652 }
1653
1654 if (png_gamma_significant(g) != 0)
1655 {
1656 back_1.red = png_gamma_8bit_correct(png_ptr->background.red,
1657 g);
1658 back_1.green = png_gamma_8bit_correct(
1659 png_ptr->background.green, g);
1660 back_1.blue = png_gamma_8bit_correct(png_ptr->background.blue,
1661 g);
1662 }
1663
1664 else
1665 {
1666 back_1.red = (png_byte)png_ptr->background.red;
1667 back_1.green = (png_byte)png_ptr->background.green;
1668 back_1.blue = (png_byte)png_ptr->background.blue;
1669 }
1670 }
1671
1672 for (i = 0; i < num_palette; i++)
1673 {
1674 if (i < (int)png_ptr->num_trans &&
1675 png_ptr->trans_alpha[i] != 0xff)
1676 {
1677 if (png_ptr->trans_alpha[i] == 0)
1678 {
1679 palette[i] = back;
1680 }
1681 else /* if (png_ptr->trans_alpha[i] != 0xff) */
1682 {
1683 png_byte v, w;
1684
1685 v = png_ptr->gamma_to_1[palette[i].red];
1686 png_composite(w, v, png_ptr->trans_alpha[i], back_1.red);
1687 palette[i].red = png_ptr->gamma_from_1[w];
1688
1689 v = png_ptr->gamma_to_1[palette[i].green];
1690 png_composite(w, v, png_ptr->trans_alpha[i], back_1.green);
1691 palette[i].green = png_ptr->gamma_from_1[w];
1692
1693 v = png_ptr->gamma_to_1[palette[i].blue];
1694 png_composite(w, v, png_ptr->trans_alpha[i], back_1.blue);
1695 palette[i].blue = png_ptr->gamma_from_1[w];
1696 }
1697 }
1698 else
1699 {
1700 palette[i].red = png_ptr->gamma_table[palette[i].red];
1701 palette[i].green = png_ptr->gamma_table[palette[i].green];
1702 palette[i].blue = png_ptr->gamma_table[palette[i].blue];
1703 }
1704 }
1705
1706 /* Prevent the transformations being done again.
1707 *
1708 * NOTE: this is highly dubious; it removes the transformations in
1709 * place. This seems inconsistent with the general treatment of the
1710 * transformations elsewhere.
1711 */
1712 png_ptr->transformations &= ~(PNG_COMPOSE | PNG_GAMMA);
1713 } /* color_type == PNG_COLOR_TYPE_PALETTE */
1714
1715 /* if (png_ptr->background_gamma_type!=PNG_BACKGROUND_GAMMA_UNKNOWN) */
1716 else /* color_type != PNG_COLOR_TYPE_PALETTE */
1717 {
1718 int gs_sig, g_sig;
1719 png_fixed_point g = PNG_FP_1; /* Correction to linear */
1720 png_fixed_point gs = PNG_FP_1; /* Correction to screen */
1721
1722 switch (png_ptr->background_gamma_type)
1723 {
1724 case PNG_BACKGROUND_GAMMA_SCREEN:
1725 g = png_ptr->screen_gamma;
1726 /* gs = PNG_FP_1; */
1727 break;
1728
1729 case PNG_BACKGROUND_GAMMA_FILE:
1730 g = png_reciprocal(png_ptr->colorspace.gamma);
1731 gs = png_reciprocal2(png_ptr->colorspace.gamma,
1732 png_ptr->screen_gamma);
1733 break;
1734
1735 case PNG_BACKGROUND_GAMMA_UNIQUE:
1736 g = png_reciprocal(png_ptr->background_gamma);
1737 gs = png_reciprocal2(png_ptr->background_gamma,
1738 png_ptr->screen_gamma);
1739 break;
1740
1741 default:
1742 png_error(png_ptr, "invalid background gamma type");
1743 }
1744
1745 g_sig = png_gamma_significant(g);
1746 gs_sig = png_gamma_significant(gs);
1747
1748 if (g_sig != 0)
1749 png_ptr->background_1.gray = png_gamma_correct(png_ptr,
1750 png_ptr->background.gray, g);
1751
1752 if (gs_sig != 0)
1753 png_ptr->background.gray = png_gamma_correct(png_ptr,
1754 png_ptr->background.gray, gs);
1755
1756 if ((png_ptr->background.red != png_ptr->background.green) ||
1757 (png_ptr->background.red != png_ptr->background.blue) ||
1758 (png_ptr->background.red != png_ptr->background.gray))
1759 {
1760 /* RGB or RGBA with color background */
1761 if (g_sig != 0)
1762 {
1763 png_ptr->background_1.red = png_gamma_correct(png_ptr,
1764 png_ptr->background.red, g);
1765
1766 png_ptr->background_1.green = png_gamma_correct(png_ptr,
1767 png_ptr->background.green, g);
1768
1769 png_ptr->background_1.blue = png_gamma_correct(png_ptr,
1770 png_ptr->background.blue, g);
1771 }
1772
1773 if (gs_sig != 0)
1774 {
1775 png_ptr->background.red = png_gamma_correct(png_ptr,
1776 png_ptr->background.red, gs);
1777
1778 png_ptr->background.green = png_gamma_correct(png_ptr,
1779 png_ptr->background.green, gs);
1780
1781 png_ptr->background.blue = png_gamma_correct(png_ptr,
1782 png_ptr->background.blue, gs);
1783 }
1784 }
1785
1786 else
1787 {
1788 /* GRAY, GRAY ALPHA, RGB, or RGBA with gray background */
1789 png_ptr->background_1.red = png_ptr->background_1.green
1790 = png_ptr->background_1.blue = png_ptr->background_1.gray;
1791
1792 png_ptr->background.red = png_ptr->background.green
1793 = png_ptr->background.blue = png_ptr->background.gray;
1794 }
1795
1796 /* The background is now in screen gamma: */
1797 png_ptr->background_gamma_type = PNG_BACKGROUND_GAMMA_SCREEN;
1798 } /* color_type != PNG_COLOR_TYPE_PALETTE */
1799 }/* png_ptr->transformations & PNG_BACKGROUND */
1800
1801 else
1802 /* Transformation does not include PNG_BACKGROUND */
1803#endif /* READ_BACKGROUND */
1804 if (png_ptr->color_type == PNG_COLOR_TYPE_PALETTE
1805#ifdef PNG_READ_RGB_TO_GRAY_SUPPORTED
1806 /* RGB_TO_GRAY needs to have non-gamma-corrected values! */
1807 && ((png_ptr->transformations & PNG_EXPAND) == 0 ||
1808 (png_ptr->transformations & PNG_RGB_TO_GRAY) == 0)
1809#endif
1810 )
1811 {
1812 png_colorp palette = png_ptr->palette;
1813 int num_palette = png_ptr->num_palette;
1814 int i;
1815
1816 /* NOTE: there are other transformations that should probably be in
1817 * here too.
1818 */
1819 for (i = 0; i < num_palette; i++)
1820 {
1821 palette[i].red = png_ptr->gamma_table[palette[i].red];
1822 palette[i].green = png_ptr->gamma_table[palette[i].green];
1823 palette[i].blue = png_ptr->gamma_table[palette[i].blue];
1824 }
1825
1826 /* Done the gamma correction. */
1827 png_ptr->transformations &= ~PNG_GAMMA;
1828 } /* color_type == PALETTE && !PNG_BACKGROUND transformation */
1829 }
1830#ifdef PNG_READ_BACKGROUND_SUPPORTED
1831 else
1832#endif
1833#endif /* READ_GAMMA */
1834
1835#ifdef PNG_READ_BACKGROUND_SUPPORTED
1836 /* No GAMMA transformation (see the hanging else 4 lines above) */
1837 if ((png_ptr->transformations & PNG_COMPOSE) != 0 &&
1838 (png_ptr->color_type == PNG_COLOR_TYPE_PALETTE))
1839 {
1840 int i;
1841 int istop = (int)png_ptr->num_trans;
1842 png_color back;
1843 png_colorp palette = png_ptr->palette;
1844
1845 back.red = (png_byte)png_ptr->background.red;
1846 back.green = (png_byte)png_ptr->background.green;
1847 back.blue = (png_byte)png_ptr->background.blue;
1848
1849 for (i = 0; i < istop; i++)
1850 {
1851 if (png_ptr->trans_alpha[i] == 0)
1852 {
1853 palette[i] = back;
1854 }
1855
1856 else if (png_ptr->trans_alpha[i] != 0xff)
1857 {
1858 /* The png_composite() macro is defined in png.h */
1859 png_composite(palette[i].red, palette[i].red,
1860 png_ptr->trans_alpha[i], back.red);
1861
1862 png_composite(palette[i].green, palette[i].green,
1863 png_ptr->trans_alpha[i], back.green);
1864
1865 png_composite(palette[i].blue, palette[i].blue,
1866 png_ptr->trans_alpha[i], back.blue);
1867 }
1868 }
1869
1870 png_ptr->transformations &= ~PNG_COMPOSE;
1871 }
1872#endif /* READ_BACKGROUND */
1873
1874#ifdef PNG_READ_SHIFT_SUPPORTED
1875 if ((png_ptr->transformations & PNG_SHIFT) != 0 &&
1876 (png_ptr->transformations & PNG_EXPAND) == 0 &&
1877 (png_ptr->color_type == PNG_COLOR_TYPE_PALETTE))
1878 {
1879 int i;
1880 int istop = png_ptr->num_palette;
1881 int shift = 8 - png_ptr->sig_bit.red;
1882
1883 png_ptr->transformations &= ~PNG_SHIFT;
1884
1885 /* significant bits can be in the range 1 to 7 for a meaninful result, if
1886 * the number of significant bits is 0 then no shift is done (this is an
1887 * error condition which is silently ignored.)
1888 */
1889 if (shift > 0 && shift < 8)
1890 for (i=0; i<istop; ++i)
1891 {
1892 int component = png_ptr->palette[i].red;
1893
1894 component >>= shift;
1895 png_ptr->palette[i].red = (png_byte)component;
1896 }
1897
1898 shift = 8 - png_ptr->sig_bit.green;
1899 if (shift > 0 && shift < 8)
1900 for (i=0; i<istop; ++i)
1901 {
1902 int component = png_ptr->palette[i].green;
1903
1904 component >>= shift;
1905 png_ptr->palette[i].green = (png_byte)component;
1906 }
1907
1908 shift = 8 - png_ptr->sig_bit.blue;
1909 if (shift > 0 && shift < 8)
1910 for (i=0; i<istop; ++i)
1911 {
1912 int component = png_ptr->palette[i].blue;
1913
1914 component >>= shift;
1915 png_ptr->palette[i].blue = (png_byte)component;
1916 }
1917 }
1918#endif /* READ_SHIFT */
1919}
1920
1921/* Modify the info structure to reflect the transformations. The
1922 * info should be updated so a PNG file could be written with it,
1923 * assuming the transformations result in valid PNG data.
1924 */
1925void /* PRIVATE */
1926png_read_transform_info(png_structrp png_ptr, png_inforp info_ptr)
1927{
1928 png_debug(1, "in png_read_transform_info");
1929
1930#ifdef PNG_READ_EXPAND_SUPPORTED
1931 if ((png_ptr->transformations & PNG_EXPAND) != 0)
1932 {
1933 if (info_ptr->color_type == PNG_COLOR_TYPE_PALETTE)
1934 {
1935 /* This check must match what actually happens in
1936 * png_do_expand_palette; if it ever checks the tRNS chunk to see if
1937 * it is all opaque we must do the same (at present it does not.)
1938 */
1939 if (png_ptr->num_trans > 0)
1940 info_ptr->color_type = PNG_COLOR_TYPE_RGB_ALPHA;
1941
1942 else
1943 info_ptr->color_type = PNG_COLOR_TYPE_RGB;
1944
1945 info_ptr->bit_depth = 8;
1946 info_ptr->num_trans = 0;
1947
1948 if (png_ptr->palette == NULL)
1949 png_error (png_ptr, "Palette is NULL in indexed image");
1950 }
1951 else
1952 {
1953 if (png_ptr->num_trans != 0)
1954 {
1955 if ((png_ptr->transformations & PNG_EXPAND_tRNS) != 0)
1956 info_ptr->color_type |= PNG_COLOR_MASK_ALPHA;
1957 }
1958 if (info_ptr->bit_depth < 8)
1959 info_ptr->bit_depth = 8;
1960
1961 info_ptr->num_trans = 0;
1962 }
1963 }
1964#endif
1965
1966#if defined(PNG_READ_BACKGROUND_SUPPORTED) ||\
1967 defined(PNG_READ_ALPHA_MODE_SUPPORTED)
1968 /* The following is almost certainly wrong unless the background value is in
1969 * the screen space!
1970 */
1971 if ((png_ptr->transformations & PNG_COMPOSE) != 0)
1972 info_ptr->background = png_ptr->background;
1973#endif
1974
1975#ifdef PNG_READ_GAMMA_SUPPORTED
1976 /* The following used to be conditional on PNG_GAMMA (prior to 1.5.4),
1977 * however it seems that the code in png_init_read_transformations, which has
1978 * been called before this from png_read_update_info->png_read_start_row
1979 * sometimes does the gamma transform and cancels the flag.
1980 *
1981 * TODO: this looks wrong; the info_ptr should end up with a gamma equal to
1982 * the screen_gamma value. The following probably results in weirdness if
1983 * the info_ptr is used by the app after the rows have been read.
1984 */
1985 info_ptr->colorspace.gamma = png_ptr->colorspace.gamma;
1986#endif
1987
1988 if (info_ptr->bit_depth == 16)
1989 {
1990# ifdef PNG_READ_16BIT_SUPPORTED
1991# ifdef PNG_READ_SCALE_16_TO_8_SUPPORTED
1992 if ((png_ptr->transformations & PNG_SCALE_16_TO_8) != 0)
1993 info_ptr->bit_depth = 8;
1994# endif
1995
1996# ifdef PNG_READ_STRIP_16_TO_8_SUPPORTED
1997 if ((png_ptr->transformations & PNG_16_TO_8) != 0)
1998 info_ptr->bit_depth = 8;
1999# endif
2000
2001# else
2002 /* No 16-bit support: force chopping 16-bit input down to 8, in this case
2003 * the app program can chose if both APIs are available by setting the
2004 * correct scaling to use.
2005 */
2006# ifdef PNG_READ_STRIP_16_TO_8_SUPPORTED
2007 /* For compatibility with previous versions use the strip method by
2008 * default. This code works because if PNG_SCALE_16_TO_8 is already
2009 * set the code below will do that in preference to the chop.
2010 */
2011 png_ptr->transformations |= PNG_16_TO_8;
2012 info_ptr->bit_depth = 8;
2013# else
2014
2015# ifdef PNG_READ_SCALE_16_TO_8_SUPPORTED
2016 png_ptr->transformations |= PNG_SCALE_16_TO_8;
2017 info_ptr->bit_depth = 8;
2018# else
2019
2020 CONFIGURATION ERROR: you must enable at least one 16 to 8 method
2021# endif
2022# endif
2023#endif /* !READ_16BIT */
2024 }
2025
2026#ifdef PNG_READ_GRAY_TO_RGB_SUPPORTED
2027 if ((png_ptr->transformations & PNG_GRAY_TO_RGB) != 0)
2028 info_ptr->color_type = (png_byte)(info_ptr->color_type |
2029 PNG_COLOR_MASK_COLOR);
2030#endif
2031
2032#ifdef PNG_READ_RGB_TO_GRAY_SUPPORTED
2033 if ((png_ptr->transformations & PNG_RGB_TO_GRAY) != 0)
2034 info_ptr->color_type = (png_byte)(info_ptr->color_type &
2035 ~PNG_COLOR_MASK_COLOR);
2036#endif
2037
2038#ifdef PNG_READ_QUANTIZE_SUPPORTED
2039 if ((png_ptr->transformations & PNG_QUANTIZE) != 0)
2040 {
2041 if (((info_ptr->color_type == PNG_COLOR_TYPE_RGB) ||
2042 (info_ptr->color_type == PNG_COLOR_TYPE_RGB_ALPHA)) &&
2043 png_ptr->palette_lookup != 0 && info_ptr->bit_depth == 8)
2044 {
2045 info_ptr->color_type = PNG_COLOR_TYPE_PALETTE;
2046 }
2047 }
2048#endif
2049
2050#ifdef PNG_READ_EXPAND_16_SUPPORTED
2051 if ((png_ptr->transformations & PNG_EXPAND_16) != 0 &&
2052 info_ptr->bit_depth == 8 &&
2053 info_ptr->color_type != PNG_COLOR_TYPE_PALETTE)
2054 {
2055 info_ptr->bit_depth = 16;
2056 }
2057#endif
2058
2059#ifdef PNG_READ_PACK_SUPPORTED
2060 if ((png_ptr->transformations & PNG_PACK) != 0 &&
2061 (info_ptr->bit_depth < 8))
2062 info_ptr->bit_depth = 8;
2063#endif
2064
2065 if (info_ptr->color_type == PNG_COLOR_TYPE_PALETTE)
2066 info_ptr->channels = 1;
2067
2068 else if ((info_ptr->color_type & PNG_COLOR_MASK_COLOR) != 0)
2069 info_ptr->channels = 3;
2070
2071 else
2072 info_ptr->channels = 1;
2073
2074#ifdef PNG_READ_STRIP_ALPHA_SUPPORTED
2075 if ((png_ptr->transformations & PNG_STRIP_ALPHA) != 0)
2076 {
2077 info_ptr->color_type = (png_byte)(info_ptr->color_type &
2078 ~PNG_COLOR_MASK_ALPHA);
2079 info_ptr->num_trans = 0;
2080 }
2081#endif
2082
2083 if ((info_ptr->color_type & PNG_COLOR_MASK_ALPHA) != 0)
2084 info_ptr->channels++;
2085
2086#ifdef PNG_READ_FILLER_SUPPORTED
2087 /* STRIP_ALPHA and FILLER allowed: MASK_ALPHA bit stripped above */
2088 if ((png_ptr->transformations & PNG_FILLER) != 0 &&
2089 (info_ptr->color_type == PNG_COLOR_TYPE_RGB ||
2090 info_ptr->color_type == PNG_COLOR_TYPE_GRAY))
2091 {
2092 info_ptr->channels++;
2093 /* If adding a true alpha channel not just filler */
2094 if ((png_ptr->transformations & PNG_ADD_ALPHA) != 0)
2095 info_ptr->color_type |= PNG_COLOR_MASK_ALPHA;
2096 }
2097#endif
2098
2099#if defined(PNG_USER_TRANSFORM_PTR_SUPPORTED) && \
2100defined(PNG_READ_USER_TRANSFORM_SUPPORTED)
2101 if ((png_ptr->transformations & PNG_USER_TRANSFORM) != 0)
2102 {
2103 if (png_ptr->user_transform_depth != 0)
2104 info_ptr->bit_depth = png_ptr->user_transform_depth;
2105
2106 if (png_ptr->user_transform_channels != 0)
2107 info_ptr->channels = png_ptr->user_transform_channels;
2108 }
2109#endif
2110
2111 info_ptr->pixel_depth = (png_byte)(info_ptr->channels *
2112 info_ptr->bit_depth);
2113
2114 info_ptr->rowbytes = PNG_ROWBYTES(info_ptr->pixel_depth, info_ptr->width);
2115
2116 /* Adding in 1.5.4: cache the above value in png_struct so that we can later
2117 * check in png_rowbytes that the user buffer won't get overwritten. Note
2118 * that the field is not always set - if png_read_update_info isn't called
2119 * the application has to either not do any transforms or get the calculation
2120 * right itself.
2121 */
2122 png_ptr->info_rowbytes = info_ptr->rowbytes;
2123
2124#ifndef PNG_READ_EXPAND_SUPPORTED
2125 if (png_ptr != NULL)
2126 return;
2127#endif
2128}
2129
2130#ifdef PNG_READ_PACK_SUPPORTED
2131/* Unpack pixels of 1, 2, or 4 bits per pixel into 1 byte per pixel,
2132 * without changing the actual values. Thus, if you had a row with
2133 * a bit depth of 1, you would end up with bytes that only contained
2134 * the numbers 0 or 1. If you would rather they contain 0 and 255, use
2135 * png_do_shift() after this.
2136 */
2137static void
2138png_do_unpack(png_row_infop row_info, png_bytep row)
2139{
2140 png_debug(1, "in png_do_unpack");
2141
2142 if (row_info->bit_depth < 8)
2143 {
2144 png_uint_32 i;
2145 png_uint_32 row_width=row_info->width;
2146
2147 switch (row_info->bit_depth)
2148 {
2149 case 1:
2150 {
2151 png_bytep sp = row + (png_size_t)((row_width - 1) >> 3);
2152 png_bytep dp = row + (png_size_t)row_width - 1;
2153 png_uint_32 shift = 7 - (int)((row_width + 7) & 0x07);
2154 for (i = 0; i < row_width; i++)
2155 {
2156 *dp = (png_byte)((*sp >> shift) & 0x01);
2157
2158 if (shift == 7)
2159 {
2160 shift = 0;
2161 sp--;
2162 }
2163
2164 else
2165 shift++;
2166
2167 dp--;
2168 }
2169 break;
2170 }
2171
2172 case 2:
2173 {
2174
2175 png_bytep sp = row + (png_size_t)((row_width - 1) >> 2);
2176 png_bytep dp = row + (png_size_t)row_width - 1;
2177 png_uint_32 shift = (int)((3 - ((row_width + 3) & 0x03)) << 1);
2178 for (i = 0; i < row_width; i++)
2179 {
2180 *dp = (png_byte)((*sp >> shift) & 0x03);
2181
2182 if (shift == 6)
2183 {
2184 shift = 0;
2185 sp--;
2186 }
2187
2188 else
2189 shift += 2;
2190
2191 dp--;
2192 }
2193 break;
2194 }
2195
2196 case 4:
2197 {
2198 png_bytep sp = row + (png_size_t)((row_width - 1) >> 1);
2199 png_bytep dp = row + (png_size_t)row_width - 1;
2200 png_uint_32 shift = (int)((1 - ((row_width + 1) & 0x01)) << 2);
2201 for (i = 0; i < row_width; i++)
2202 {
2203 *dp = (png_byte)((*sp >> shift) & 0x0f);
2204
2205 if (shift == 4)
2206 {
2207 shift = 0;
2208 sp--;
2209 }
2210
2211 else
2212 shift = 4;
2213
2214 dp--;
2215 }
2216 break;
2217 }
2218
2219 default:
2220 break;
2221 }
2222 row_info->bit_depth = 8;
2223 row_info->pixel_depth = (png_byte)(8 * row_info->channels);
2224 row_info->rowbytes = row_width * row_info->channels;
2225 }
2226}
2227#endif
2228
2229#ifdef PNG_READ_SHIFT_SUPPORTED
2230/* Reverse the effects of png_do_shift. This routine merely shifts the
2231 * pixels back to their significant bits values. Thus, if you have
2232 * a row of bit depth 8, but only 5 are significant, this will shift
2233 * the values back to 0 through 31.
2234 */
2235static void
2236png_do_unshift(png_row_infop row_info, png_bytep row,
2237 png_const_color_8p sig_bits)
2238{
2239 int color_type;
2240
2241 png_debug(1, "in png_do_unshift");
2242
2243 /* The palette case has already been handled in the _init routine. */
2244 color_type = row_info->color_type;
2245
2246 if (color_type != PNG_COLOR_TYPE_PALETTE)
2247 {
2248 int shift[4];
2249 int channels = 0;
2250 int bit_depth = row_info->bit_depth;
2251
2252 if ((color_type & PNG_COLOR_MASK_COLOR) != 0)
2253 {
2254 shift[channels++] = bit_depth - sig_bits->red;
2255 shift[channels++] = bit_depth - sig_bits->green;
2256 shift[channels++] = bit_depth - sig_bits->blue;
2257 }
2258
2259 else
2260 {
2261 shift[channels++] = bit_depth - sig_bits->gray;
2262 }
2263
2264 if ((color_type & PNG_COLOR_MASK_ALPHA) != 0)
2265 {
2266 shift[channels++] = bit_depth - sig_bits->alpha;
2267 }
2268
2269 {
2270 int c, have_shift;
2271
2272 for (c = have_shift = 0; c < channels; ++c)
2273 {
2274 /* A shift of more than the bit depth is an error condition but it
2275 * gets ignored here.
2276 */
2277 if (shift[c] <= 0 || shift[c] >= bit_depth)
2278 shift[c] = 0;
2279
2280 else
2281 have_shift = 1;
2282 }
2283
2284 if (have_shift == 0)
2285 return;
2286 }
2287
2288 switch (bit_depth)
2289 {
2290 default:
2291 /* Must be 1bpp gray: should not be here! */
2292 /* NOTREACHED */
2293 break;
2294
2295 case 2:
2296 /* Must be 2bpp gray */
2297 /* assert(channels == 1 && shift[0] == 1) */
2298 {
2299 png_bytep bp = row;
2300 png_bytep bp_end = bp + row_info->rowbytes;
2301
2302 while (bp < bp_end)
2303 {
2304 int b = (*bp >> 1) & 0x55;
2305 *bp++ = (png_byte)b;
2306 }
2307 break;
2308 }
2309
2310 case 4:
2311 /* Must be 4bpp gray */
2312 /* assert(channels == 1) */
2313 {
2314 png_bytep bp = row;
2315 png_bytep bp_end = bp + row_info->rowbytes;
2316 int gray_shift = shift[0];
2317 int mask = 0xf >> gray_shift;
2318
2319 mask |= mask << 4;
2320
2321 while (bp < bp_end)
2322 {
2323 int b = (*bp >> gray_shift) & mask;
2324 *bp++ = (png_byte)b;
2325 }
2326 break;
2327 }
2328
2329 case 8:
2330 /* Single byte components, G, GA, RGB, RGBA */
2331 {
2332 png_bytep bp = row;
2333 png_bytep bp_end = bp + row_info->rowbytes;
2334 int channel = 0;
2335
2336 while (bp < bp_end)
2337 {
2338 int b = *bp >> shift[channel];
2339 if (++channel >= channels)
2340 channel = 0;
2341 *bp++ = (png_byte)b;
2342 }
2343 break;
2344 }
2345
2346#ifdef PNG_READ_16BIT_SUPPORTED
2347 case 16:
2348 /* Double byte components, G, GA, RGB, RGBA */
2349 {
2350 png_bytep bp = row;
2351 png_bytep bp_end = bp + row_info->rowbytes;
2352 int channel = 0;
2353
2354 while (bp < bp_end)
2355 {
2356 int value = (bp[0] << 8) + bp[1];
2357
2358 value >>= shift[channel];
2359 if (++channel >= channels)
2360 channel = 0;
2361 *bp++ = (png_byte)(value >> 8);
2362 *bp++ = (png_byte)value;
2363 }
2364 break;
2365 }
2366#endif
2367 }
2368 }
2369}
2370#endif
2371
2372#ifdef PNG_READ_SCALE_16_TO_8_SUPPORTED
2373/* Scale rows of bit depth 16 down to 8 accurately */
2374static void
2375png_do_scale_16_to_8(png_row_infop row_info, png_bytep row)
2376{
2377 png_debug(1, "in png_do_scale_16_to_8");
2378
2379 if (row_info->bit_depth == 16)
2380 {
2381 png_bytep sp = row; /* source */
2382 png_bytep dp = row; /* destination */
2383 png_bytep ep = sp + row_info->rowbytes; /* end+1 */
2384
2385 while (sp < ep)
2386 {
2387 /* The input is an array of 16-bit components, these must be scaled to
2388 * 8 bits each. For a 16-bit value V the required value (from the PNG
2389 * specification) is:
2390 *
2391 * (V * 255) / 65535
2392 *
2393 * This reduces to round(V / 257), or floor((V + 128.5)/257)
2394 *
2395 * Represent V as the two byte value vhi.vlo. Make a guess that the
2396 * result is the top byte of V, vhi, then the correction to this value
2397 * is:
2398 *
2399 * error = floor(((V-vhi.vhi) + 128.5) / 257)
2400 * = floor(((vlo-vhi) + 128.5) / 257)
2401 *
2402 * This can be approximated using integer arithmetic (and a signed
2403 * shift):
2404 *
2405 * error = (vlo-vhi+128) >> 8;
2406 *
2407 * The approximate differs from the exact answer only when (vlo-vhi) is
2408 * 128; it then gives a correction of +1 when the exact correction is
2409 * 0. This gives 128 errors. The exact answer (correct for all 16-bit
2410 * input values) is:
2411 *
2412 * error = (vlo-vhi+128)*65535 >> 24;
2413 *
2414 * An alternative arithmetic calculation which also gives no errors is:
2415 *
2416 * (V * 255 + 32895) >> 16
2417 */
2418
2419 png_int_32 tmp = *sp++; /* must be signed! */
2420 tmp += (((int)*sp++ - tmp + 128) * 65535) >> 24;
2421 *dp++ = (png_byte)tmp;
2422 }
2423
2424 row_info->bit_depth = 8;
2425 row_info->pixel_depth = (png_byte)(8 * row_info->channels);
2426 row_info->rowbytes = row_info->width * row_info->channels;
2427 }
2428}
2429#endif
2430
2431#ifdef PNG_READ_STRIP_16_TO_8_SUPPORTED
2432static void
2433/* Simply discard the low byte. This was the default behavior prior
2434 * to libpng-1.5.4.
2435 */
2436png_do_chop(png_row_infop row_info, png_bytep row)
2437{
2438 png_debug(1, "in png_do_chop");
2439
2440 if (row_info->bit_depth == 16)
2441 {
2442 png_bytep sp = row; /* source */
2443 png_bytep dp = row; /* destination */
2444 png_bytep ep = sp + row_info->rowbytes; /* end+1 */
2445
2446 while (sp < ep)
2447 {
2448 *dp++ = *sp;
2449 sp += 2; /* skip low byte */
2450 }
2451
2452 row_info->bit_depth = 8;
2453 row_info->pixel_depth = (png_byte)(8 * row_info->channels);
2454 row_info->rowbytes = row_info->width * row_info->channels;
2455 }
2456}
2457#endif
2458
2459#ifdef PNG_READ_SWAP_ALPHA_SUPPORTED
2460static void
2461png_do_read_swap_alpha(png_row_infop row_info, png_bytep row)
2462{
2463 png_debug(1, "in png_do_read_swap_alpha");
2464
2465 {
2466 png_uint_32 row_width = row_info->width;
2467 if (row_info->color_type == PNG_COLOR_TYPE_RGB_ALPHA)
2468 {
2469 /* This converts from RGBA to ARGB */
2470 if (row_info->bit_depth == 8)
2471 {
2472 png_bytep sp = row + row_info->rowbytes;
2473 png_bytep dp = sp;
2474 png_byte save;
2475 png_uint_32 i;
2476
2477 for (i = 0; i < row_width; i++)
2478 {
2479 save = *(--sp);
2480 *(--dp) = *(--sp);
2481 *(--dp) = *(--sp);
2482 *(--dp) = *(--sp);
2483 *(--dp) = save;
2484 }
2485 }
2486
2487#ifdef PNG_READ_16BIT_SUPPORTED
2488 /* This converts from RRGGBBAA to AARRGGBB */
2489 else
2490 {
2491 png_bytep sp = row + row_info->rowbytes;
2492 png_bytep dp = sp;
2493 png_byte save[2];
2494 png_uint_32 i;
2495
2496 for (i = 0; i < row_width; i++)
2497 {
2498 save[0] = *(--sp);
2499 save[1] = *(--sp);
2500 *(--dp) = *(--sp);
2501 *(--dp) = *(--sp);
2502 *(--dp) = *(--sp);
2503 *(--dp) = *(--sp);
2504 *(--dp) = *(--sp);
2505 *(--dp) = *(--sp);
2506 *(--dp) = save[0];
2507 *(--dp) = save[1];
2508 }
2509 }
2510#endif
2511 }
2512
2513 else if (row_info->color_type == PNG_COLOR_TYPE_GRAY_ALPHA)
2514 {
2515 /* This converts from GA to AG */
2516 if (row_info->bit_depth == 8)
2517 {
2518 png_bytep sp = row + row_info->rowbytes;
2519 png_bytep dp = sp;
2520 png_byte save;
2521 png_uint_32 i;
2522
2523 for (i = 0; i < row_width; i++)
2524 {
2525 save = *(--sp);
2526 *(--dp) = *(--sp);
2527 *(--dp) = save;
2528 }
2529 }
2530
2531#ifdef PNG_READ_16BIT_SUPPORTED
2532 /* This converts from GGAA to AAGG */
2533 else
2534 {
2535 png_bytep sp = row + row_info->rowbytes;
2536 png_bytep dp = sp;
2537 png_byte save[2];
2538 png_uint_32 i;
2539
2540 for (i = 0; i < row_width; i++)
2541 {
2542 save[0] = *(--sp);
2543 save[1] = *(--sp);
2544 *(--dp) = *(--sp);
2545 *(--dp) = *(--sp);
2546 *(--dp) = save[0];
2547 *(--dp) = save[1];
2548 }
2549 }
2550#endif
2551 }
2552 }
2553}
2554#endif
2555
2556#ifdef PNG_READ_INVERT_ALPHA_SUPPORTED
2557static void
2558png_do_read_invert_alpha(png_row_infop row_info, png_bytep row)
2559{
2560 png_uint_32 row_width;
2561 png_debug(1, "in png_do_read_invert_alpha");
2562
2563 row_width = row_info->width;
2564 if (row_info->color_type == PNG_COLOR_TYPE_RGB_ALPHA)
2565 {
2566 if (row_info->bit_depth == 8)
2567 {
2568 /* This inverts the alpha channel in RGBA */
2569 png_bytep sp = row + row_info->rowbytes;
2570 png_bytep dp = sp;
2571 png_uint_32 i;
2572
2573 for (i = 0; i < row_width; i++)
2574 {
2575 *(--dp) = (png_byte)(255 - *(--sp));
2576
2577/* This does nothing:
2578 *(--dp) = *(--sp);
2579 *(--dp) = *(--sp);
2580 *(--dp) = *(--sp);
2581 We can replace it with:
2582*/
2583 sp-=3;
2584 dp=sp;
2585 }
2586 }
2587
2588#ifdef PNG_READ_16BIT_SUPPORTED
2589 /* This inverts the alpha channel in RRGGBBAA */
2590 else
2591 {
2592 png_bytep sp = row + row_info->rowbytes;
2593 png_bytep dp = sp;
2594 png_uint_32 i;
2595
2596 for (i = 0; i < row_width; i++)
2597 {
2598 *(--dp) = (png_byte)(255 - *(--sp));
2599 *(--dp) = (png_byte)(255 - *(--sp));
2600
2601/* This does nothing:
2602 *(--dp) = *(--sp);
2603 *(--dp) = *(--sp);
2604 *(--dp) = *(--sp);
2605 *(--dp) = *(--sp);
2606 *(--dp) = *(--sp);
2607 *(--dp) = *(--sp);
2608 We can replace it with:
2609*/
2610 sp-=6;
2611 dp=sp;
2612 }
2613 }
2614#endif
2615 }
2616 else if (row_info->color_type == PNG_COLOR_TYPE_GRAY_ALPHA)
2617 {
2618 if (row_info->bit_depth == 8)
2619 {
2620 /* This inverts the alpha channel in GA */
2621 png_bytep sp = row + row_info->rowbytes;
2622 png_bytep dp = sp;
2623 png_uint_32 i;
2624
2625 for (i = 0; i < row_width; i++)
2626 {
2627 *(--dp) = (png_byte)(255 - *(--sp));
2628 *(--dp) = *(--sp);
2629 }
2630 }
2631
2632#ifdef PNG_READ_16BIT_SUPPORTED
2633 else
2634 {
2635 /* This inverts the alpha channel in GGAA */
2636 png_bytep sp = row + row_info->rowbytes;
2637 png_bytep dp = sp;
2638 png_uint_32 i;
2639
2640 for (i = 0; i < row_width; i++)
2641 {
2642 *(--dp) = (png_byte)(255 - *(--sp));
2643 *(--dp) = (png_byte)(255 - *(--sp));
2644/*
2645 *(--dp) = *(--sp);
2646 *(--dp) = *(--sp);
2647*/
2648 sp-=2;
2649 dp=sp;
2650 }
2651 }
2652#endif
2653 }
2654}
2655#endif
2656
2657#ifdef PNG_READ_FILLER_SUPPORTED
2658/* Add filler channel if we have RGB color */
2659static void
2660png_do_read_filler(png_row_infop row_info, png_bytep row,
2661 png_uint_32 filler, png_uint_32 flags)
2662{
2663 png_uint_32 i;
2664 png_uint_32 row_width = row_info->width;
2665
2666#ifdef PNG_READ_16BIT_SUPPORTED
2667 png_byte hi_filler = (png_byte)(filler>>8);
2668#endif
2669 png_byte lo_filler = (png_byte)filler;
2670
2671 png_debug(1, "in png_do_read_filler");
2672
2673 if (
2674 row_info->color_type == PNG_COLOR_TYPE_GRAY)
2675 {
2676 if (row_info->bit_depth == 8)
2677 {
2678 if ((flags & PNG_FLAG_FILLER_AFTER) != 0)
2679 {
2680 /* This changes the data from G to GX */
2681 png_bytep sp = row + (png_size_t)row_width;
2682 png_bytep dp = sp + (png_size_t)row_width;
2683 for (i = 1; i < row_width; i++)
2684 {
2685 *(--dp) = lo_filler;
2686 *(--dp) = *(--sp);
2687 }
2688 *(--dp) = lo_filler;
2689 row_info->channels = 2;
2690 row_info->pixel_depth = 16;
2691 row_info->rowbytes = row_width * 2;
2692 }
2693
2694 else
2695 {
2696 /* This changes the data from G to XG */
2697 png_bytep sp = row + (png_size_t)row_width;
2698 png_bytep dp = sp + (png_size_t)row_width;
2699 for (i = 0; i < row_width; i++)
2700 {
2701 *(--dp) = *(--sp);
2702 *(--dp) = lo_filler;
2703 }
2704 row_info->channels = 2;
2705 row_info->pixel_depth = 16;
2706 row_info->rowbytes = row_width * 2;
2707 }
2708 }
2709
2710#ifdef PNG_READ_16BIT_SUPPORTED
2711 else if (row_info->bit_depth == 16)
2712 {
2713 if ((flags & PNG_FLAG_FILLER_AFTER) != 0)
2714 {
2715 /* This changes the data from GG to GGXX */
2716 png_bytep sp = row + (png_size_t)row_width * 2;
2717 png_bytep dp = sp + (png_size_t)row_width * 2;
2718 for (i = 1; i < row_width; i++)
2719 {
2720 *(--dp) = lo_filler;
2721 *(--dp) = hi_filler;
2722 *(--dp) = *(--sp);
2723 *(--dp) = *(--sp);
2724 }
2725 *(--dp) = lo_filler;
2726 *(--dp) = hi_filler;
2727 row_info->channels = 2;
2728 row_info->pixel_depth = 32;
2729 row_info->rowbytes = row_width * 4;
2730 }
2731
2732 else
2733 {
2734 /* This changes the data from GG to XXGG */
2735 png_bytep sp = row + (png_size_t)row_width * 2;
2736 png_bytep dp = sp + (png_size_t)row_width * 2;
2737 for (i = 0; i < row_width; i++)
2738 {
2739 *(--dp) = *(--sp);
2740 *(--dp) = *(--sp);
2741 *(--dp) = lo_filler;
2742 *(--dp) = hi_filler;
2743 }
2744 row_info->channels = 2;
2745 row_info->pixel_depth = 32;
2746 row_info->rowbytes = row_width * 4;
2747 }
2748 }
2749#endif
2750 } /* COLOR_TYPE == GRAY */
2751 else if (row_info->color_type == PNG_COLOR_TYPE_RGB)
2752 {
2753 if (row_info->bit_depth == 8)
2754 {
2755 if ((flags & PNG_FLAG_FILLER_AFTER) != 0)
2756 {
2757 /* This changes the data from RGB to RGBX */
2758 png_bytep sp = row + (png_size_t)row_width * 3;
2759 png_bytep dp = sp + (png_size_t)row_width;
2760 for (i = 1; i < row_width; i++)
2761 {
2762 *(--dp) = lo_filler;
2763 *(--dp) = *(--sp);
2764 *(--dp) = *(--sp);
2765 *(--dp) = *(--sp);
2766 }
2767 *(--dp) = lo_filler;
2768 row_info->channels = 4;
2769 row_info->pixel_depth = 32;
2770 row_info->rowbytes = row_width * 4;
2771 }
2772
2773 else
2774 {
2775 /* This changes the data from RGB to XRGB */
2776 png_bytep sp = row + (png_size_t)row_width * 3;
2777 png_bytep dp = sp + (png_size_t)row_width;
2778 for (i = 0; i < row_width; i++)
2779 {
2780 *(--dp) = *(--sp);
2781 *(--dp) = *(--sp);
2782 *(--dp) = *(--sp);
2783 *(--dp) = lo_filler;
2784 }
2785 row_info->channels = 4;
2786 row_info->pixel_depth = 32;
2787 row_info->rowbytes = row_width * 4;
2788 }
2789 }
2790
2791#ifdef PNG_READ_16BIT_SUPPORTED
2792 else if (row_info->bit_depth == 16)
2793 {
2794 if ((flags & PNG_FLAG_FILLER_AFTER) != 0)
2795 {
2796 /* This changes the data from RRGGBB to RRGGBBXX */
2797 png_bytep sp = row + (png_size_t)row_width * 6;
2798 png_bytep dp = sp + (png_size_t)row_width * 2;
2799 for (i = 1; i < row_width; i++)
2800 {
2801 *(--dp) = lo_filler;
2802 *(--dp) = hi_filler;
2803 *(--dp) = *(--sp);
2804 *(--dp) = *(--sp);
2805 *(--dp) = *(--sp);
2806 *(--dp) = *(--sp);
2807 *(--dp) = *(--sp);
2808 *(--dp) = *(--sp);
2809 }
2810 *(--dp) = lo_filler;
2811 *(--dp) = hi_filler;
2812 row_info->channels = 4;
2813 row_info->pixel_depth = 64;
2814 row_info->rowbytes = row_width * 8;
2815 }
2816
2817 else
2818 {
2819 /* This changes the data from RRGGBB to XXRRGGBB */
2820 png_bytep sp = row + (png_size_t)row_width * 6;
2821 png_bytep dp = sp + (png_size_t)row_width * 2;
2822 for (i = 0; i < row_width; i++)
2823 {
2824 *(--dp) = *(--sp);
2825 *(--dp) = *(--sp);
2826 *(--dp) = *(--sp);
2827 *(--dp) = *(--sp);
2828 *(--dp) = *(--sp);
2829 *(--dp) = *(--sp);
2830 *(--dp) = lo_filler;
2831 *(--dp) = hi_filler;
2832 }
2833
2834 row_info->channels = 4;
2835 row_info->pixel_depth = 64;
2836 row_info->rowbytes = row_width * 8;
2837 }
2838 }
2839#endif
2840 } /* COLOR_TYPE == RGB */
2841}
2842#endif
2843
2844#ifdef PNG_READ_GRAY_TO_RGB_SUPPORTED
2845/* Expand grayscale files to RGB, with or without alpha */
2846static void
2847png_do_gray_to_rgb(png_row_infop row_info, png_bytep row)
2848{
2849 png_uint_32 i;
2850 png_uint_32 row_width = row_info->width;
2851
2852 png_debug(1, "in png_do_gray_to_rgb");
2853
2854 if (row_info->bit_depth >= 8 &&
2855 (row_info->color_type & PNG_COLOR_MASK_COLOR) == 0)
2856 {
2857 if (row_info->color_type == PNG_COLOR_TYPE_GRAY)
2858 {
2859 if (row_info->bit_depth == 8)
2860 {
2861 /* This changes G to RGB */
2862 png_bytep sp = row + (png_size_t)row_width - 1;
2863 png_bytep dp = sp + (png_size_t)row_width * 2;
2864 for (i = 0; i < row_width; i++)
2865 {
2866 *(dp--) = *sp;
2867 *(dp--) = *sp;
2868 *(dp--) = *(sp--);
2869 }
2870 }
2871
2872 else
2873 {
2874 /* This changes GG to RRGGBB */
2875 png_bytep sp = row + (png_size_t)row_width * 2 - 1;
2876 png_bytep dp = sp + (png_size_t)row_width * 4;
2877 for (i = 0; i < row_width; i++)
2878 {
2879 *(dp--) = *sp;
2880 *(dp--) = *(sp - 1);
2881 *(dp--) = *sp;
2882 *(dp--) = *(sp - 1);
2883 *(dp--) = *(sp--);
2884 *(dp--) = *(sp--);
2885 }
2886 }
2887 }
2888
2889 else if (row_info->color_type == PNG_COLOR_TYPE_GRAY_ALPHA)
2890 {
2891 if (row_info->bit_depth == 8)
2892 {
2893 /* This changes GA to RGBA */
2894 png_bytep sp = row + (png_size_t)row_width * 2 - 1;
2895 png_bytep dp = sp + (png_size_t)row_width * 2;
2896 for (i = 0; i < row_width; i++)
2897 {
2898 *(dp--) = *(sp--);
2899 *(dp--) = *sp;
2900 *(dp--) = *sp;
2901 *(dp--) = *(sp--);
2902 }
2903 }
2904
2905 else
2906 {
2907 /* This changes GGAA to RRGGBBAA */
2908 png_bytep sp = row + (png_size_t)row_width * 4 - 1;
2909 png_bytep dp = sp + (png_size_t)row_width * 4;
2910 for (i = 0; i < row_width; i++)
2911 {
2912 *(dp--) = *(sp--);
2913 *(dp--) = *(sp--);
2914 *(dp--) = *sp;
2915 *(dp--) = *(sp - 1);
2916 *(dp--) = *sp;
2917 *(dp--) = *(sp - 1);
2918 *(dp--) = *(sp--);
2919 *(dp--) = *(sp--);
2920 }
2921 }
2922 }
2923 row_info->channels = (png_byte)(row_info->channels + 2);
2924 row_info->color_type |= PNG_COLOR_MASK_COLOR;
2925 row_info->pixel_depth = (png_byte)(row_info->channels *
2926 row_info->bit_depth);
2927 row_info->rowbytes = PNG_ROWBYTES(row_info->pixel_depth, row_width);
2928 }
2929}
2930#endif
2931
2932#ifdef PNG_READ_RGB_TO_GRAY_SUPPORTED
2933/* Reduce RGB files to grayscale, with or without alpha
2934 * using the equation given in Poynton's ColorFAQ of 1998-01-04 at
2935 * <http://www.inforamp.net/~poynton/> (THIS LINK IS DEAD June 2008 but
2936 * versions dated 1998 through November 2002 have been archived at
2937 * http://web.archive.org/web/20000816232553/http://www.inforamp.net/
2938 * ~poynton/notes/colour_and_gamma/ColorFAQ.txt )
2939 * Charles Poynton poynton at poynton.com
2940 *
2941 * Y = 0.212671 * R + 0.715160 * G + 0.072169 * B
2942 *
2943 * which can be expressed with integers as
2944 *
2945 * Y = (6969 * R + 23434 * G + 2365 * B)/32768
2946 *
2947 * Poynton's current link (as of January 2003 through July 2011):
2948 * <http://www.poynton.com/notes/colour_and_gamma/>
2949 * has changed the numbers slightly:
2950 *
2951 * Y = 0.2126*R + 0.7152*G + 0.0722*B
2952 *
2953 * which can be expressed with integers as
2954 *
2955 * Y = (6966 * R + 23436 * G + 2366 * B)/32768
2956 *
2957 * Historically, however, libpng uses numbers derived from the ITU-R Rec 709
2958 * end point chromaticities and the D65 white point. Depending on the
2959 * precision used for the D65 white point this produces a variety of different
2960 * numbers, however if the four decimal place value used in ITU-R Rec 709 is
2961 * used (0.3127,0.3290) the Y calculation would be:
2962 *
2963 * Y = (6968 * R + 23435 * G + 2366 * B)/32768
2964 *
2965 * While this is correct the rounding results in an overflow for white, because
2966 * the sum of the rounded coefficients is 32769, not 32768. Consequently
2967 * libpng uses, instead, the closest non-overflowing approximation:
2968 *
2969 * Y = (6968 * R + 23434 * G + 2366 * B)/32768
2970 *
2971 * Starting with libpng-1.5.5, if the image being converted has a cHRM chunk
2972 * (including an sRGB chunk) then the chromaticities are used to calculate the
2973 * coefficients. See the chunk handling in pngrutil.c for more information.
2974 *
2975 * In all cases the calculation is to be done in a linear colorspace. If no
2976 * gamma information is available to correct the encoding of the original RGB
2977 * values this results in an implicit assumption that the original PNG RGB
2978 * values were linear.
2979 *
2980 * Other integer coefficents can be used via png_set_rgb_to_gray(). Because
2981 * the API takes just red and green coefficients the blue coefficient is
2982 * calculated to make the sum 32768. This will result in different rounding
2983 * to that used above.
2984 */
2985static int
2986png_do_rgb_to_gray(png_structrp png_ptr, png_row_infop row_info, png_bytep row)
2987
2988{
2989 int rgb_error = 0;
2990
2991 png_debug(1, "in png_do_rgb_to_gray");
2992
2993 if ((row_info->color_type & PNG_COLOR_MASK_PALETTE) == 0 &&
2994 (row_info->color_type & PNG_COLOR_MASK_COLOR) != 0)
2995 {
2996 PNG_CONST png_uint_32 rc = png_ptr->rgb_to_gray_red_coeff;
2997 PNG_CONST png_uint_32 gc = png_ptr->rgb_to_gray_green_coeff;
2998 PNG_CONST png_uint_32 bc = 32768 - rc - gc;
2999 PNG_CONST png_uint_32 row_width = row_info->width;
3000 PNG_CONST int have_alpha =
3001 (row_info->color_type & PNG_COLOR_MASK_ALPHA) != 0;
3002
3003 if (row_info->bit_depth == 8)
3004 {
3005#ifdef PNG_READ_GAMMA_SUPPORTED
3006 /* Notice that gamma to/from 1 are not necessarily inverses (if
3007 * there is an overall gamma correction). Prior to 1.5.5 this code
3008 * checked the linearized values for equality; this doesn't match
3009 * the documentation, the original values must be checked.
3010 */
3011 if (png_ptr->gamma_from_1 != NULL && png_ptr->gamma_to_1 != NULL)
3012 {
3013 png_bytep sp = row;
3014 png_bytep dp = row;
3015 png_uint_32 i;
3016
3017 for (i = 0; i < row_width; i++)
3018 {
3019 png_byte red = *(sp++);
3020 png_byte green = *(sp++);
3021 png_byte blue = *(sp++);
3022
3023 if (red != green || red != blue)
3024 {
3025 red = png_ptr->gamma_to_1[red];
3026 green = png_ptr->gamma_to_1[green];
3027 blue = png_ptr->gamma_to_1[blue];
3028
3029 rgb_error |= 1;
3030 *(dp++) = png_ptr->gamma_from_1[
3031 (rc*red + gc*green + bc*blue + 16384)>>15];
3032 }
3033
3034 else
3035 {
3036 /* If there is no overall correction the table will not be
3037 * set.
3038 */
3039 if (png_ptr->gamma_table != NULL)
3040 red = png_ptr->gamma_table[red];
3041
3042 *(dp++) = red;
3043 }
3044
3045 if (have_alpha != 0)
3046 *(dp++) = *(sp++);
3047 }
3048 }
3049 else
3050#endif
3051 {
3052 png_bytep sp = row;
3053 png_bytep dp = row;
3054 png_uint_32 i;
3055
3056 for (i = 0; i < row_width; i++)
3057 {
3058 png_byte red = *(sp++);
3059 png_byte green = *(sp++);
3060 png_byte blue = *(sp++);
3061
3062 if (red != green || red != blue)
3063 {
3064 rgb_error |= 1;
3065 /* NOTE: this is the historical approach which simply
3066 * truncates the results.
3067 */
3068 *(dp++) = (png_byte)((rc*red + gc*green + bc*blue)>>15);
3069 }
3070
3071 else
3072 *(dp++) = red;
3073
3074 if (have_alpha != 0)
3075 *(dp++) = *(sp++);
3076 }
3077 }
3078 }
3079
3080 else /* RGB bit_depth == 16 */
3081 {
3082#ifdef PNG_READ_GAMMA_SUPPORTED
3083 if (png_ptr->gamma_16_to_1 != NULL && png_ptr->gamma_16_from_1 != NULL)
3084 {
3085 png_bytep sp = row;
3086 png_bytep dp = row;
3087 png_uint_32 i;
3088
3089 for (i = 0; i < row_width; i++)
3090 {
3091 png_uint_16 red, green, blue, w;
3092 png_byte hi,lo;
3093
3094 hi=*(sp)++; lo=*(sp)++; red = (png_uint_16)((hi << 8) | (lo));
3095 hi=*(sp)++; lo=*(sp)++; green = (png_uint_16)((hi << 8) | (lo));
3096 hi=*(sp)++; lo=*(sp)++; blue = (png_uint_16)((hi << 8) | (lo));
3097
3098 if (red == green && red == blue)
3099 {
3100 if (png_ptr->gamma_16_table != NULL)
3101 w = png_ptr->gamma_16_table[(red & 0xff)
3102 >> png_ptr->gamma_shift][red >> 8];
3103
3104 else
3105 w = red;
3106 }
3107
3108 else
3109 {
3110 png_uint_16 red_1 = png_ptr->gamma_16_to_1[(red & 0xff)
3111 >> png_ptr->gamma_shift][red>>8];
3112 png_uint_16 green_1 =
3113 png_ptr->gamma_16_to_1[(green & 0xff) >>
3114 png_ptr->gamma_shift][green>>8];
3115 png_uint_16 blue_1 = png_ptr->gamma_16_to_1[(blue & 0xff)
3116 >> png_ptr->gamma_shift][blue>>8];
3117 png_uint_16 gray16 = (png_uint_16)((rc*red_1 + gc*green_1
3118 + bc*blue_1 + 16384)>>15);
3119 w = png_ptr->gamma_16_from_1[(gray16 & 0xff) >>
3120 png_ptr->gamma_shift][gray16 >> 8];
3121 rgb_error |= 1;
3122 }
3123
3124 *(dp++) = (png_byte)((w>>8) & 0xff);
3125 *(dp++) = (png_byte)(w & 0xff);
3126
3127 if (have_alpha != 0)
3128 {
3129 *(dp++) = *(sp++);
3130 *(dp++) = *(sp++);
3131 }
3132 }
3133 }
3134 else
3135#endif
3136 {
3137 png_bytep sp = row;
3138 png_bytep dp = row;
3139 png_uint_32 i;
3140
3141 for (i = 0; i < row_width; i++)
3142 {
3143 png_uint_16 red, green, blue, gray16;
3144 png_byte hi,lo;
3145
3146 hi=*(sp)++; lo=*(sp)++; red = (png_uint_16)((hi << 8) | (lo));
3147 hi=*(sp)++; lo=*(sp)++; green = (png_uint_16)((hi << 8) | (lo));
3148 hi=*(sp)++; lo=*(sp)++; blue = (png_uint_16)((hi << 8) | (lo));
3149
3150 if (red != green || red != blue)
3151 rgb_error |= 1;
3152
3153 /* From 1.5.5 in the 16-bit case do the accurate conversion even
3154 * in the 'fast' case - this is because this is where the code
3155 * ends up when handling linear 16-bit data.
3156 */
3157 gray16 = (png_uint_16)((rc*red + gc*green + bc*blue + 16384) >>
3158 15);
3159 *(dp++) = (png_byte)((gray16 >> 8) & 0xff);
3160 *(dp++) = (png_byte)(gray16 & 0xff);
3161
3162 if (have_alpha != 0)
3163 {
3164 *(dp++) = *(sp++);
3165 *(dp++) = *(sp++);
3166 }
3167 }
3168 }
3169 }
3170
3171 row_info->channels = (png_byte)(row_info->channels - 2);
3172 row_info->color_type = (png_byte)(row_info->color_type &
3173 ~PNG_COLOR_MASK_COLOR);
3174 row_info->pixel_depth = (png_byte)(row_info->channels *
3175 row_info->bit_depth);
3176 row_info->rowbytes = PNG_ROWBYTES(row_info->pixel_depth, row_width);
3177 }
3178 return rgb_error;
3179}
3180#endif
3181
3182#if defined(PNG_READ_BACKGROUND_SUPPORTED) ||\
3183 defined(PNG_READ_ALPHA_MODE_SUPPORTED)
3184/* Replace any alpha or transparency with the supplied background color.
3185 * "background" is already in the screen gamma, while "background_1" is
3186 * at a gamma of 1.0. Paletted files have already been taken care of.
3187 */
3188static void
3189png_do_compose(png_row_infop row_info, png_bytep row, png_structrp png_ptr)
3190{
3191#ifdef PNG_READ_GAMMA_SUPPORTED
3192 png_const_bytep gamma_table = png_ptr->gamma_table;
3193 png_const_bytep gamma_from_1 = png_ptr->gamma_from_1;
3194 png_const_bytep gamma_to_1 = png_ptr->gamma_to_1;
3195 png_const_uint_16pp gamma_16 = png_ptr->gamma_16_table;
3196 png_const_uint_16pp gamma_16_from_1 = png_ptr->gamma_16_from_1;
3197 png_const_uint_16pp gamma_16_to_1 = png_ptr->gamma_16_to_1;
3198 int gamma_shift = png_ptr->gamma_shift;
3199 int optimize = (png_ptr->flags & PNG_FLAG_OPTIMIZE_ALPHA) != 0;
3200#endif
3201
3202 png_bytep sp;
3203 png_uint_32 i;
3204 png_uint_32 row_width = row_info->width;
3205 int shift;
3206
3207 png_debug(1, "in png_do_compose");
3208
3209 {
3210 switch (row_info->color_type)
3211 {
3212 case PNG_COLOR_TYPE_GRAY:
3213 {
3214 switch (row_info->bit_depth)
3215 {
3216 case 1:
3217 {
3218 sp = row;
3219 shift = 7;
3220 for (i = 0; i < row_width; i++)
3221 {
3222 if ((png_uint_16)((*sp >> shift) & 0x01)
3223 == png_ptr->trans_color.gray)
3224 {
3225 unsigned int tmp = *sp & (0x7f7f >> (7 - shift));
3226 tmp |= png_ptr->background.gray << shift;
3227 *sp = (png_byte)(tmp & 0xff);
3228 }
3229
3230 if (shift == 0)
3231 {
3232 shift = 7;
3233 sp++;
3234 }
3235
3236 else
3237 shift--;
3238 }
3239 break;
3240 }
3241
3242 case 2:
3243 {
3244#ifdef PNG_READ_GAMMA_SUPPORTED
3245 if (gamma_table != NULL)
3246 {
3247 sp = row;
3248 shift = 6;
3249 for (i = 0; i < row_width; i++)
3250 {
3251 if ((png_uint_16)((*sp >> shift) & 0x03)
3252 == png_ptr->trans_color.gray)
3253 {
3254 unsigned int tmp = *sp & (0x3f3f >> (6 - shift));
3255 tmp |= png_ptr->background.gray << shift;
3256 *sp = (png_byte)(tmp & 0xff);
3257 }
3258
3259 else
3260 {
3261 unsigned int p = (*sp >> shift) & 0x03;
3262 unsigned int g = (gamma_table [p | (p << 2) |
3263 (p << 4) | (p << 6)] >> 6) & 0x03;
3264 unsigned int tmp = *sp & (0x3f3f >> (6 - shift));
3265 tmp |= g << shift;
3266 *sp = (png_byte)(tmp & 0xff);
3267 }
3268
3269 if (shift == 0)
3270 {
3271 shift = 6;
3272 sp++;
3273 }
3274
3275 else
3276 shift -= 2;
3277 }
3278 }
3279
3280 else
3281#endif
3282 {
3283 sp = row;
3284 shift = 6;
3285 for (i = 0; i < row_width; i++)
3286 {
3287 if ((png_uint_16)((*sp >> shift) & 0x03)
3288 == png_ptr->trans_color.gray)
3289 {
3290 unsigned int tmp = *sp & (0x3f3f >> (6 - shift));
3291 tmp |= png_ptr->background.gray << shift;
3292 *sp = (png_byte)(tmp & 0xff);
3293 }
3294
3295 if (shift == 0)
3296 {
3297 shift = 6;
3298 sp++;
3299 }
3300
3301 else
3302 shift -= 2;
3303 }
3304 }
3305 break;
3306 }
3307
3308 case 4:
3309 {
3310#ifdef PNG_READ_GAMMA_SUPPORTED
3311 if (gamma_table != NULL)
3312 {
3313 sp = row;
3314 shift = 4;
3315 for (i = 0; i < row_width; i++)
3316 {
3317 if ((png_uint_16)((*sp >> shift) & 0x0f)
3318 == png_ptr->trans_color.gray)
3319 {
3320 unsigned int tmp = *sp & (0x0f0f >> (4 - shift));
3321 tmp |= png_ptr->background.gray << shift;
3322 *sp = (png_byte)(tmp & 0xff);
3323 }
3324
3325 else
3326 {
3327 unsigned int p = (*sp >> shift) & 0x0f;
3328 unsigned int g = (gamma_table[p | (p << 4)] >> 4) &
3329 0x0f;
3330 unsigned int tmp = *sp & (0x0f0f >> (4 - shift));
3331 tmp |= g << shift;
3332 *sp = (png_byte)(tmp & 0xff);
3333 }
3334
3335 if (shift == 0)
3336 {
3337 shift = 4;
3338 sp++;
3339 }
3340
3341 else
3342 shift -= 4;
3343 }
3344 }
3345
3346 else
3347#endif
3348 {
3349 sp = row;
3350 shift = 4;
3351 for (i = 0; i < row_width; i++)
3352 {
3353 if ((png_uint_16)((*sp >> shift) & 0x0f)
3354 == png_ptr->trans_color.gray)
3355 {
3356 unsigned int tmp = *sp & (0x0f0f >> (4 - shift));
3357 tmp |= png_ptr->background.gray << shift;
3358 *sp = (png_byte)(tmp & 0xff);
3359 }
3360
3361 if (shift == 0)
3362 {
3363 shift = 4;
3364 sp++;
3365 }
3366
3367 else
3368 shift -= 4;
3369 }
3370 }
3371 break;
3372 }
3373
3374 case 8:
3375 {
3376#ifdef PNG_READ_GAMMA_SUPPORTED
3377 if (gamma_table != NULL)
3378 {
3379 sp = row;
3380 for (i = 0; i < row_width; i++, sp++)
3381 {
3382 if (*sp == png_ptr->trans_color.gray)
3383 *sp = (png_byte)png_ptr->background.gray;
3384
3385 else
3386 *sp = gamma_table[*sp];
3387 }
3388 }
3389 else
3390#endif
3391 {
3392 sp = row;
3393 for (i = 0; i < row_width; i++, sp++)
3394 {
3395 if (*sp == png_ptr->trans_color.gray)
3396 *sp = (png_byte)png_ptr->background.gray;
3397 }
3398 }
3399 break;
3400 }
3401
3402 case 16:
3403 {
3404#ifdef PNG_READ_GAMMA_SUPPORTED
3405 if (gamma_16 != NULL)
3406 {
3407 sp = row;
3408 for (i = 0; i < row_width; i++, sp += 2)
3409 {
3410 png_uint_16 v;
3411
3412 v = (png_uint_16)(((*sp) << 8) + *(sp + 1));
3413
3414 if (v == png_ptr->trans_color.gray)
3415 {
3416 /* Background is already in screen gamma */
3417 *sp = (png_byte)((png_ptr->background.gray >> 8)
3418 & 0xff);
3419 *(sp + 1) = (png_byte)(png_ptr->background.gray
3420 & 0xff);
3421 }
3422
3423 else
3424 {
3425 v = gamma_16[*(sp + 1) >> gamma_shift][*sp];
3426 *sp = (png_byte)((v >> 8) & 0xff);
3427 *(sp + 1) = (png_byte)(v & 0xff);
3428 }
3429 }
3430 }
3431 else
3432#endif
3433 {
3434 sp = row;
3435 for (i = 0; i < row_width; i++, sp += 2)
3436 {
3437 png_uint_16 v;
3438
3439 v = (png_uint_16)(((*sp) << 8) + *(sp + 1));
3440
3441 if (v == png_ptr->trans_color.gray)
3442 {
3443 *sp = (png_byte)((png_ptr->background.gray >> 8)
3444 & 0xff);
3445 *(sp + 1) = (png_byte)(png_ptr->background.gray
3446 & 0xff);
3447 }
3448 }
3449 }
3450 break;
3451 }
3452
3453 default:
3454 break;
3455 }
3456 break;
3457 }
3458
3459 case PNG_COLOR_TYPE_RGB:
3460 {
3461 if (row_info->bit_depth == 8)
3462 {
3463#ifdef PNG_READ_GAMMA_SUPPORTED
3464 if (gamma_table != NULL)
3465 {
3466 sp = row;
3467 for (i = 0; i < row_width; i++, sp += 3)
3468 {
3469 if (*sp == png_ptr->trans_color.red &&
3470 *(sp + 1) == png_ptr->trans_color.green &&
3471 *(sp + 2) == png_ptr->trans_color.blue)
3472 {
3473 *sp = (png_byte)png_ptr->background.red;
3474 *(sp + 1) = (png_byte)png_ptr->background.green;
3475 *(sp + 2) = (png_byte)png_ptr->background.blue;
3476 }
3477
3478 else
3479 {
3480 *sp = gamma_table[*sp];
3481 *(sp + 1) = gamma_table[*(sp + 1)];
3482 *(sp + 2) = gamma_table[*(sp + 2)];
3483 }
3484 }
3485 }
3486 else
3487#endif
3488 {
3489 sp = row;
3490 for (i = 0; i < row_width; i++, sp += 3)
3491 {
3492 if (*sp == png_ptr->trans_color.red &&
3493 *(sp + 1) == png_ptr->trans_color.green &&
3494 *(sp + 2) == png_ptr->trans_color.blue)
3495 {
3496 *sp = (png_byte)png_ptr->background.red;
3497 *(sp + 1) = (png_byte)png_ptr->background.green;
3498 *(sp + 2) = (png_byte)png_ptr->background.blue;
3499 }
3500 }
3501 }
3502 }
3503 else /* if (row_info->bit_depth == 16) */
3504 {
3505#ifdef PNG_READ_GAMMA_SUPPORTED
3506 if (gamma_16 != NULL)
3507 {
3508 sp = row;
3509 for (i = 0; i < row_width; i++, sp += 6)
3510 {
3511 png_uint_16 r = (png_uint_16)(((*sp) << 8) + *(sp + 1));
3512
3513 png_uint_16 g = (png_uint_16)(((*(sp + 2)) << 8)
3514 + *(sp + 3));
3515
3516 png_uint_16 b = (png_uint_16)(((*(sp + 4)) << 8)
3517 + *(sp + 5));
3518
3519 if (r == png_ptr->trans_color.red &&
3520 g == png_ptr->trans_color.green &&
3521 b == png_ptr->trans_color.blue)
3522 {
3523 /* Background is already in screen gamma */
3524 *sp = (png_byte)((png_ptr->background.red >> 8) & 0xff);
3525 *(sp + 1) = (png_byte)(png_ptr->background.red & 0xff);
3526 *(sp + 2) = (png_byte)((png_ptr->background.green >> 8)
3527 & 0xff);
3528 *(sp + 3) = (png_byte)(png_ptr->background.green
3529 & 0xff);
3530 *(sp + 4) = (png_byte)((png_ptr->background.blue >> 8)
3531 & 0xff);
3532 *(sp + 5) = (png_byte)(png_ptr->background.blue & 0xff);
3533 }
3534
3535 else
3536 {
3537 png_uint_16 v = gamma_16[*(sp + 1) >> gamma_shift][*sp];
3538 *sp = (png_byte)((v >> 8) & 0xff);
3539 *(sp + 1) = (png_byte)(v & 0xff);
3540
3541 v = gamma_16[*(sp + 3) >> gamma_shift][*(sp + 2)];
3542 *(sp + 2) = (png_byte)((v >> 8) & 0xff);
3543 *(sp + 3) = (png_byte)(v & 0xff);
3544
3545 v = gamma_16[*(sp + 5) >> gamma_shift][*(sp + 4)];
3546 *(sp + 4) = (png_byte)((v >> 8) & 0xff);
3547 *(sp + 5) = (png_byte)(v & 0xff);
3548 }
3549 }
3550 }
3551
3552 else
3553#endif
3554 {
3555 sp = row;
3556 for (i = 0; i < row_width; i++, sp += 6)
3557 {
3558 png_uint_16 r = (png_uint_16)(((*sp) << 8) + *(sp + 1));
3559
3560 png_uint_16 g = (png_uint_16)(((*(sp + 2)) << 8)
3561 + *(sp + 3));
3562
3563 png_uint_16 b = (png_uint_16)(((*(sp + 4)) << 8)
3564 + *(sp + 5));
3565
3566 if (r == png_ptr->trans_color.red &&
3567 g == png_ptr->trans_color.green &&
3568 b == png_ptr->trans_color.blue)
3569 {
3570 *sp = (png_byte)((png_ptr->background.red >> 8) & 0xff);
3571 *(sp + 1) = (png_byte)(png_ptr->background.red & 0xff);
3572 *(sp + 2) = (png_byte)((png_ptr->background.green >> 8)
3573 & 0xff);
3574 *(sp + 3) = (png_byte)(png_ptr->background.green
3575 & 0xff);
3576 *(sp + 4) = (png_byte)((png_ptr->background.blue >> 8)
3577 & 0xff);
3578 *(sp + 5) = (png_byte)(png_ptr->background.blue & 0xff);
3579 }
3580 }
3581 }
3582 }
3583 break;
3584 }
3585
3586 case PNG_COLOR_TYPE_GRAY_ALPHA:
3587 {
3588 if (row_info->bit_depth == 8)
3589 {
3590#ifdef PNG_READ_GAMMA_SUPPORTED
3591 if (gamma_to_1 != NULL && gamma_from_1 != NULL &&
3592 gamma_table != NULL)
3593 {
3594 sp = row;
3595 for (i = 0; i < row_width; i++, sp += 2)
3596 {
3597 png_uint_16 a = *(sp + 1);
3598
3599 if (a == 0xff)
3600 *sp = gamma_table[*sp];
3601
3602 else if (a == 0)
3603 {
3604 /* Background is already in screen gamma */
3605 *sp = (png_byte)png_ptr->background.gray;
3606 }
3607
3608 else
3609 {
3610 png_byte v, w;
3611
3612 v = gamma_to_1[*sp];
3613 png_composite(w, v, a, png_ptr->background_1.gray);
3614 if (optimize == 0)
3615 w = gamma_from_1[w];
3616 *sp = w;
3617 }
3618 }
3619 }
3620 else
3621#endif
3622 {
3623 sp = row;
3624 for (i = 0; i < row_width; i++, sp += 2)
3625 {
3626 png_byte a = *(sp + 1);
3627
3628 if (a == 0)
3629 *sp = (png_byte)png_ptr->background.gray;
3630
3631 else if (a < 0xff)
3632 png_composite(*sp, *sp, a, png_ptr->background.gray);
3633 }
3634 }
3635 }
3636 else /* if (png_ptr->bit_depth == 16) */
3637 {
3638#ifdef PNG_READ_GAMMA_SUPPORTED
3639 if (gamma_16 != NULL && gamma_16_from_1 != NULL &&
3640 gamma_16_to_1 != NULL)
3641 {
3642 sp = row;
3643 for (i = 0; i < row_width; i++, sp += 4)
3644 {
3645 png_uint_16 a = (png_uint_16)(((*(sp + 2)) << 8)
3646 + *(sp + 3));
3647
3648 if (a == (png_uint_16)0xffff)
3649 {
3650 png_uint_16 v;
3651
3652 v = gamma_16[*(sp + 1) >> gamma_shift][*sp];
3653 *sp = (png_byte)((v >> 8) & 0xff);
3654 *(sp + 1) = (png_byte)(v & 0xff);
3655 }
3656
3657 else if (a == 0)
3658 {
3659 /* Background is already in screen gamma */
3660 *sp = (png_byte)((png_ptr->background.gray >> 8)
3661 & 0xff);
3662 *(sp + 1) = (png_byte)(png_ptr->background.gray & 0xff);
3663 }
3664
3665 else
3666 {
3667 png_uint_16 g, v, w;
3668
3669 g = gamma_16_to_1[*(sp + 1) >> gamma_shift][*sp];
3670 png_composite_16(v, g, a, png_ptr->background_1.gray);
3671 if (optimize != 0)
3672 w = v;
3673 else
3674 w = gamma_16_from_1[(v & 0xff) >>
3675 gamma_shift][v >> 8];
3676 *sp = (png_byte)((w >> 8) & 0xff);
3677 *(sp + 1) = (png_byte)(w & 0xff);
3678 }
3679 }
3680 }
3681 else
3682#endif
3683 {
3684 sp = row;
3685 for (i = 0; i < row_width; i++, sp += 4)
3686 {
3687 png_uint_16 a = (png_uint_16)(((*(sp + 2)) << 8)
3688 + *(sp + 3));
3689
3690 if (a == 0)
3691 {
3692 *sp = (png_byte)((png_ptr->background.gray >> 8)
3693 & 0xff);
3694 *(sp + 1) = (png_byte)(png_ptr->background.gray & 0xff);
3695 }
3696
3697 else if (a < 0xffff)
3698 {
3699 png_uint_16 g, v;
3700
3701 g = (png_uint_16)(((*sp) << 8) + *(sp + 1));
3702 png_composite_16(v, g, a, png_ptr->background.gray);
3703 *sp = (png_byte)((v >> 8) & 0xff);
3704 *(sp + 1) = (png_byte)(v & 0xff);
3705 }
3706 }
3707 }
3708 }
3709 break;
3710 }
3711
3712 case PNG_COLOR_TYPE_RGB_ALPHA:
3713 {
3714 if (row_info->bit_depth == 8)
3715 {
3716#ifdef PNG_READ_GAMMA_SUPPORTED
3717 if (gamma_to_1 != NULL && gamma_from_1 != NULL &&
3718 gamma_table != NULL)
3719 {
3720 sp = row;
3721 for (i = 0; i < row_width; i++, sp += 4)
3722 {
3723 png_byte a = *(sp + 3);
3724
3725 if (a == 0xff)
3726 {
3727 *sp = gamma_table[*sp];
3728 *(sp + 1) = gamma_table[*(sp + 1)];
3729 *(sp + 2) = gamma_table[*(sp + 2)];
3730 }
3731
3732 else if (a == 0)
3733 {
3734 /* Background is already in screen gamma */
3735 *sp = (png_byte)png_ptr->background.red;
3736 *(sp + 1) = (png_byte)png_ptr->background.green;
3737 *(sp + 2) = (png_byte)png_ptr->background.blue;
3738 }
3739
3740 else
3741 {
3742 png_byte v, w;
3743
3744 v = gamma_to_1[*sp];
3745 png_composite(w, v, a, png_ptr->background_1.red);
3746 if (optimize == 0) w = gamma_from_1[w];
3747 *sp = w;
3748
3749 v = gamma_to_1[*(sp + 1)];
3750 png_composite(w, v, a, png_ptr->background_1.green);
3751 if (optimize == 0) w = gamma_from_1[w];
3752 *(sp + 1) = w;
3753
3754 v = gamma_to_1[*(sp + 2)];
3755 png_composite(w, v, a, png_ptr->background_1.blue);
3756 if (optimize == 0) w = gamma_from_1[w];
3757 *(sp + 2) = w;
3758 }
3759 }
3760 }
3761 else
3762#endif
3763 {
3764 sp = row;
3765 for (i = 0; i < row_width; i++, sp += 4)
3766 {
3767 png_byte a = *(sp + 3);
3768
3769 if (a == 0)
3770 {
3771 *sp = (png_byte)png_ptr->background.red;
3772 *(sp + 1) = (png_byte)png_ptr->background.green;
3773 *(sp + 2) = (png_byte)png_ptr->background.blue;
3774 }
3775
3776 else if (a < 0xff)
3777 {
3778 png_composite(*sp, *sp, a, png_ptr->background.red);
3779
3780 png_composite(*(sp + 1), *(sp + 1), a,
3781 png_ptr->background.green);
3782
3783 png_composite(*(sp + 2), *(sp + 2), a,
3784 png_ptr->background.blue);
3785 }
3786 }
3787 }
3788 }
3789 else /* if (row_info->bit_depth == 16) */
3790 {
3791#ifdef PNG_READ_GAMMA_SUPPORTED
3792 if (gamma_16 != NULL && gamma_16_from_1 != NULL &&
3793 gamma_16_to_1 != NULL)
3794 {
3795 sp = row;
3796 for (i = 0; i < row_width; i++, sp += 8)
3797 {
3798 png_uint_16 a = (png_uint_16)(((png_uint_16)(*(sp + 6))
3799 << 8) + (png_uint_16)(*(sp + 7)));
3800
3801 if (a == (png_uint_16)0xffff)
3802 {
3803 png_uint_16 v;
3804
3805 v = gamma_16[*(sp + 1) >> gamma_shift][*sp];
3806 *sp = (png_byte)((v >> 8) & 0xff);
3807 *(sp + 1) = (png_byte)(v & 0xff);
3808
3809 v = gamma_16[*(sp + 3) >> gamma_shift][*(sp + 2)];
3810 *(sp + 2) = (png_byte)((v >> 8) & 0xff);
3811 *(sp + 3) = (png_byte)(v & 0xff);
3812
3813 v = gamma_16[*(sp + 5) >> gamma_shift][*(sp + 4)];
3814 *(sp + 4) = (png_byte)((v >> 8) & 0xff);
3815 *(sp + 5) = (png_byte)(v & 0xff);
3816 }
3817
3818 else if (a == 0)
3819 {
3820 /* Background is already in screen gamma */
3821 *sp = (png_byte)((png_ptr->background.red >> 8) & 0xff);
3822 *(sp + 1) = (png_byte)(png_ptr->background.red & 0xff);
3823 *(sp + 2) = (png_byte)((png_ptr->background.green >> 8)
3824 & 0xff);
3825 *(sp + 3) = (png_byte)(png_ptr->background.green
3826 & 0xff);
3827 *(sp + 4) = (png_byte)((png_ptr->background.blue >> 8)
3828 & 0xff);
3829 *(sp + 5) = (png_byte)(png_ptr->background.blue & 0xff);
3830 }
3831
3832 else
3833 {
3834 png_uint_16 v, w;
3835
3836 v = gamma_16_to_1[*(sp + 1) >> gamma_shift][*sp];
3837 png_composite_16(w, v, a, png_ptr->background_1.red);
3838 if (optimize == 0)
3839 w = gamma_16_from_1[((w & 0xff) >> gamma_shift)][w >>
3840 8];
3841 *sp = (png_byte)((w >> 8) & 0xff);
3842 *(sp + 1) = (png_byte)(w & 0xff);
3843
3844 v = gamma_16_to_1[*(sp + 3) >> gamma_shift][*(sp + 2)];
3845 png_composite_16(w, v, a, png_ptr->background_1.green);
3846 if (optimize == 0)
3847 w = gamma_16_from_1[((w & 0xff) >> gamma_shift)][w >>
3848 8];
3849
3850 *(sp + 2) = (png_byte)((w >> 8) & 0xff);
3851 *(sp + 3) = (png_byte)(w & 0xff);
3852
3853 v = gamma_16_to_1[*(sp + 5) >> gamma_shift][*(sp + 4)];
3854 png_composite_16(w, v, a, png_ptr->background_1.blue);
3855 if (optimize == 0)
3856 w = gamma_16_from_1[((w & 0xff) >> gamma_shift)][w >>
3857 8];
3858
3859 *(sp + 4) = (png_byte)((w >> 8) & 0xff);
3860 *(sp + 5) = (png_byte)(w & 0xff);
3861 }
3862 }
3863 }
3864
3865 else
3866#endif
3867 {
3868 sp = row;
3869 for (i = 0; i < row_width; i++, sp += 8)
3870 {
3871 png_uint_16 a = (png_uint_16)(((png_uint_16)(*(sp + 6))
3872 << 8) + (png_uint_16)(*(sp + 7)));
3873
3874 if (a == 0)
3875 {
3876 *sp = (png_byte)((png_ptr->background.red >> 8) & 0xff);
3877 *(sp + 1) = (png_byte)(png_ptr->background.red & 0xff);
3878 *(sp + 2) = (png_byte)((png_ptr->background.green >> 8)
3879 & 0xff);
3880 *(sp + 3) = (png_byte)(png_ptr->background.green
3881 & 0xff);
3882 *(sp + 4) = (png_byte)((png_ptr->background.blue >> 8)
3883 & 0xff);
3884 *(sp + 5) = (png_byte)(png_ptr->background.blue & 0xff);
3885 }
3886
3887 else if (a < 0xffff)
3888 {
3889 png_uint_16 v;
3890
3891 png_uint_16 r = (png_uint_16)(((*sp) << 8) + *(sp + 1));
3892 png_uint_16 g = (png_uint_16)(((*(sp + 2)) << 8)
3893 + *(sp + 3));
3894 png_uint_16 b = (png_uint_16)(((*(sp + 4)) << 8)
3895 + *(sp + 5));
3896
3897 png_composite_16(v, r, a, png_ptr->background.red);
3898 *sp = (png_byte)((v >> 8) & 0xff);
3899 *(sp + 1) = (png_byte)(v & 0xff);
3900
3901 png_composite_16(v, g, a, png_ptr->background.green);
3902 *(sp + 2) = (png_byte)((v >> 8) & 0xff);
3903 *(sp + 3) = (png_byte)(v & 0xff);
3904
3905 png_composite_16(v, b, a, png_ptr->background.blue);
3906 *(sp + 4) = (png_byte)((v >> 8) & 0xff);
3907 *(sp + 5) = (png_byte)(v & 0xff);
3908 }
3909 }
3910 }
3911 }
3912 break;
3913 }
3914
3915 default:
3916 break;
3917 }
3918 }
3919}
3920#endif /* READ_BACKGROUND || READ_ALPHA_MODE */
3921
3922#ifdef PNG_READ_GAMMA_SUPPORTED
3923/* Gamma correct the image, avoiding the alpha channel. Make sure
3924 * you do this after you deal with the transparency issue on grayscale
3925 * or RGB images. If your bit depth is 8, use gamma_table, if it
3926 * is 16, use gamma_16_table and gamma_shift. Build these with
3927 * build_gamma_table().
3928 */
3929static void
3930png_do_gamma(png_row_infop row_info, png_bytep row, png_structrp png_ptr)
3931{
3932 png_const_bytep gamma_table = png_ptr->gamma_table;
3933 png_const_uint_16pp gamma_16_table = png_ptr->gamma_16_table;
3934 int gamma_shift = png_ptr->gamma_shift;
3935
3936 png_bytep sp;
3937 png_uint_32 i;
3938 png_uint_32 row_width=row_info->width;
3939
3940 png_debug(1, "in png_do_gamma");
3941
3942 if (((row_info->bit_depth <= 8 && gamma_table != NULL) ||
3943 (row_info->bit_depth == 16 && gamma_16_table != NULL)))
3944 {
3945 switch (row_info->color_type)
3946 {
3947 case PNG_COLOR_TYPE_RGB:
3948 {
3949 if (row_info->bit_depth == 8)
3950 {
3951 sp = row;
3952 for (i = 0; i < row_width; i++)
3953 {
3954 *sp = gamma_table[*sp];
3955 sp++;
3956 *sp = gamma_table[*sp];
3957 sp++;
3958 *sp = gamma_table[*sp];
3959 sp++;
3960 }
3961 }
3962
3963 else /* if (row_info->bit_depth == 16) */
3964 {
3965 sp = row;
3966 for (i = 0; i < row_width; i++)
3967 {
3968 png_uint_16 v;
3969
3970 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
3971 *sp = (png_byte)((v >> 8) & 0xff);
3972 *(sp + 1) = (png_byte)(v & 0xff);
3973 sp += 2;
3974
3975 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
3976 *sp = (png_byte)((v >> 8) & 0xff);
3977 *(sp + 1) = (png_byte)(v & 0xff);
3978 sp += 2;
3979
3980 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
3981 *sp = (png_byte)((v >> 8) & 0xff);
3982 *(sp + 1) = (png_byte)(v & 0xff);
3983 sp += 2;
3984 }
3985 }
3986 break;
3987 }
3988
3989 case PNG_COLOR_TYPE_RGB_ALPHA:
3990 {
3991 if (row_info->bit_depth == 8)
3992 {
3993 sp = row;
3994 for (i = 0; i < row_width; i++)
3995 {
3996 *sp = gamma_table[*sp];
3997 sp++;
3998
3999 *sp = gamma_table[*sp];
4000 sp++;
4001
4002 *sp = gamma_table[*sp];
4003 sp++;
4004
4005 sp++;
4006 }
4007 }
4008
4009 else /* if (row_info->bit_depth == 16) */
4010 {
4011 sp = row;
4012 for (i = 0; i < row_width; i++)
4013 {
4014 png_uint_16 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
4015 *sp = (png_byte)((v >> 8) & 0xff);
4016 *(sp + 1) = (png_byte)(v & 0xff);
4017 sp += 2;
4018
4019 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
4020 *sp = (png_byte)((v >> 8) & 0xff);
4021 *(sp + 1) = (png_byte)(v & 0xff);
4022 sp += 2;
4023
4024 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
4025 *sp = (png_byte)((v >> 8) & 0xff);
4026 *(sp + 1) = (png_byte)(v & 0xff);
4027 sp += 4;
4028 }
4029 }
4030 break;
4031 }
4032
4033 case PNG_COLOR_TYPE_GRAY_ALPHA:
4034 {
4035 if (row_info->bit_depth == 8)
4036 {
4037 sp = row;
4038 for (i = 0; i < row_width; i++)
4039 {
4040 *sp = gamma_table[*sp];
4041 sp += 2;
4042 }
4043 }
4044
4045 else /* if (row_info->bit_depth == 16) */
4046 {
4047 sp = row;
4048 for (i = 0; i < row_width; i++)
4049 {
4050 png_uint_16 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
4051 *sp = (png_byte)((v >> 8) & 0xff);
4052 *(sp + 1) = (png_byte)(v & 0xff);
4053 sp += 4;
4054 }
4055 }
4056 break;
4057 }
4058
4059 case PNG_COLOR_TYPE_GRAY:
4060 {
4061 if (row_info->bit_depth == 2)
4062 {
4063 sp = row;
4064 for (i = 0; i < row_width; i += 4)
4065 {
4066 int a = *sp & 0xc0;
4067 int b = *sp & 0x30;
4068 int c = *sp & 0x0c;
4069 int d = *sp & 0x03;
4070
4071 *sp = (png_byte)(
4072 ((((int)gamma_table[a|(a>>2)|(a>>4)|(a>>6)]) ) & 0xc0)|
4073 ((((int)gamma_table[(b<<2)|b|(b>>2)|(b>>4)])>>2) & 0x30)|
4074 ((((int)gamma_table[(c<<4)|(c<<2)|c|(c>>2)])>>4) & 0x0c)|
4075 ((((int)gamma_table[(d<<6)|(d<<4)|(d<<2)|d])>>6) ));
4076 sp++;
4077 }
4078 }
4079
4080 if (row_info->bit_depth == 4)
4081 {
4082 sp = row;
4083 for (i = 0; i < row_width; i += 2)
4084 {
4085 int msb = *sp & 0xf0;
4086 int lsb = *sp & 0x0f;
4087
4088 *sp = (png_byte)((((int)gamma_table[msb | (msb >> 4)]) & 0xf0)
4089 | (((int)gamma_table[(lsb << 4) | lsb]) >> 4));
4090 sp++;
4091 }
4092 }
4093
4094 else if (row_info->bit_depth == 8)
4095 {
4096 sp = row;
4097 for (i = 0; i < row_width; i++)
4098 {
4099 *sp = gamma_table[*sp];
4100 sp++;
4101 }
4102 }
4103
4104 else if (row_info->bit_depth == 16)
4105 {
4106 sp = row;
4107 for (i = 0; i < row_width; i++)
4108 {
4109 png_uint_16 v = gamma_16_table[*(sp + 1) >> gamma_shift][*sp];
4110 *sp = (png_byte)((v >> 8) & 0xff);
4111 *(sp + 1) = (png_byte)(v & 0xff);
4112 sp += 2;
4113 }
4114 }
4115 break;
4116 }
4117
4118 default:
4119 break;
4120 }
4121 }
4122}
4123#endif
4124
4125#ifdef PNG_READ_ALPHA_MODE_SUPPORTED
4126/* Encode the alpha channel to the output gamma (the input channel is always
4127 * linear.) Called only with color types that have an alpha channel. Needs the
4128 * from_1 tables.
4129 */
4130static void
4131png_do_encode_alpha(png_row_infop row_info, png_bytep row, png_structrp png_ptr)
4132{
4133 png_uint_32 row_width = row_info->width;
4134
4135 png_debug(1, "in png_do_encode_alpha");
4136
4137 if ((row_info->color_type & PNG_COLOR_MASK_ALPHA) != 0)
4138 {
4139 if (row_info->bit_depth == 8)
4140 {
4141 PNG_CONST png_bytep table = png_ptr->gamma_from_1;
4142
4143 if (table != NULL)
4144 {
4145 PNG_CONST int step =
4146 (row_info->color_type & PNG_COLOR_MASK_COLOR) ? 4 : 2;
4147
4148 /* The alpha channel is the last component: */
4149 row += step - 1;
4150
4151 for (; row_width > 0; --row_width, row += step)
4152 *row = table[*row];
4153
4154 return;
4155 }
4156 }
4157
4158 else if (row_info->bit_depth == 16)
4159 {
4160 PNG_CONST png_uint_16pp table = png_ptr->gamma_16_from_1;
4161 PNG_CONST int gamma_shift = png_ptr->gamma_shift;
4162
4163 if (table != NULL)
4164 {
4165 PNG_CONST int step =
4166 (row_info->color_type & PNG_COLOR_MASK_COLOR) ? 8 : 4;
4167
4168 /* The alpha channel is the last component: */
4169 row += step - 2;
4170
4171 for (; row_width > 0; --row_width, row += step)
4172 {
4173 png_uint_16 v;
4174
4175 v = table[*(row + 1) >> gamma_shift][*row];
4176 *row = (png_byte)((v >> 8) & 0xff);
4177 *(row + 1) = (png_byte)(v & 0xff);
4178 }
4179
4180 return;
4181 }
4182 }
4183 }
4184
4185 /* Only get to here if called with a weird row_info; no harm has been done,
4186 * so just issue a warning.
4187 */
4188 png_warning(png_ptr, "png_do_encode_alpha: unexpected call");
4189}
4190#endif
4191
4192#ifdef PNG_READ_EXPAND_SUPPORTED
4193/* Expands a palette row to an RGB or RGBA row depending
4194 * upon whether you supply trans and num_trans.
4195 */
4196static void
4197png_do_expand_palette(png_row_infop row_info, png_bytep row,
4198 png_const_colorp palette, png_const_bytep trans_alpha, int num_trans)
4199{
4200 int shift, value;
4201 png_bytep sp, dp;
4202 png_uint_32 i;
4203 png_uint_32 row_width=row_info->width;
4204
4205 png_debug(1, "in png_do_expand_palette");
4206
4207 if (row_info->color_type == PNG_COLOR_TYPE_PALETTE)
4208 {
4209 if (row_info->bit_depth < 8)
4210 {
4211 switch (row_info->bit_depth)
4212 {
4213 case 1:
4214 {
4215 sp = row + (png_size_t)((row_width - 1) >> 3);
4216 dp = row + (png_size_t)row_width - 1;
4217 shift = 7 - (int)((row_width + 7) & 0x07);
4218 for (i = 0; i < row_width; i++)
4219 {
4220 if ((*sp >> shift) & 0x01)
4221 *dp = 1;
4222
4223 else
4224 *dp = 0;
4225
4226 if (shift == 7)
4227 {
4228 shift = 0;
4229 sp--;
4230 }
4231
4232 else
4233 shift++;
4234
4235 dp--;
4236 }
4237 break;
4238 }
4239
4240 case 2:
4241 {
4242 sp = row + (png_size_t)((row_width - 1) >> 2);
4243 dp = row + (png_size_t)row_width - 1;
4244 shift = (int)((3 - ((row_width + 3) & 0x03)) << 1);
4245 for (i = 0; i < row_width; i++)
4246 {
4247 value = (*sp >> shift) & 0x03;
4248 *dp = (png_byte)value;
4249 if (shift == 6)
4250 {
4251 shift = 0;
4252 sp--;
4253 }
4254
4255 else
4256 shift += 2;
4257
4258 dp--;
4259 }
4260 break;
4261 }
4262
4263 case 4:
4264 {
4265 sp = row + (png_size_t)((row_width - 1) >> 1);
4266 dp = row + (png_size_t)row_width - 1;
4267 shift = (int)((row_width & 0x01) << 2);
4268 for (i = 0; i < row_width; i++)
4269 {
4270 value = (*sp >> shift) & 0x0f;
4271 *dp = (png_byte)value;
4272 if (shift == 4)
4273 {
4274 shift = 0;
4275 sp--;
4276 }
4277
4278 else
4279 shift += 4;
4280
4281 dp--;
4282 }
4283 break;
4284 }
4285
4286 default:
4287 break;
4288 }
4289 row_info->bit_depth = 8;
4290 row_info->pixel_depth = 8;
4291 row_info->rowbytes = row_width;
4292 }
4293
4294 if (row_info->bit_depth == 8)
4295 {
4296 {
4297 if (num_trans > 0)
4298 {
4299 sp = row + (png_size_t)row_width - 1;
4300 dp = row + (png_size_t)(row_width << 2) - 1;
4301
4302 for (i = 0; i < row_width; i++)
4303 {
4304 if ((int)(*sp) >= num_trans)
4305 *dp-- = 0xff;
4306
4307 else
4308 *dp-- = trans_alpha[*sp];
4309
4310 *dp-- = palette[*sp].blue;
4311 *dp-- = palette[*sp].green;
4312 *dp-- = palette[*sp].red;
4313 sp--;
4314 }
4315 row_info->bit_depth = 8;
4316 row_info->pixel_depth = 32;
4317 row_info->rowbytes = row_width * 4;
4318 row_info->color_type = 6;
4319 row_info->channels = 4;
4320 }
4321
4322 else
4323 {
4324 sp = row + (png_size_t)row_width - 1;
4325 dp = row + (png_size_t)(row_width * 3) - 1;
4326
4327 for (i = 0; i < row_width; i++)
4328 {
4329 *dp-- = palette[*sp].blue;
4330 *dp-- = palette[*sp].green;
4331 *dp-- = palette[*sp].red;
4332 sp--;
4333 }
4334
4335 row_info->bit_depth = 8;
4336 row_info->pixel_depth = 24;
4337 row_info->rowbytes = row_width * 3;
4338 row_info->color_type = 2;
4339 row_info->channels = 3;
4340 }
4341 }
4342 }
4343 }
4344}
4345
4346/* If the bit depth < 8, it is expanded to 8. Also, if the already
4347 * expanded transparency value is supplied, an alpha channel is built.
4348 */
4349static void
4350png_do_expand(png_row_infop row_info, png_bytep row,
4351 png_const_color_16p trans_color)
4352{
4353 int shift, value;
4354 png_bytep sp, dp;
4355 png_uint_32 i;
4356 png_uint_32 row_width=row_info->width;
4357
4358 png_debug(1, "in png_do_expand");
4359
4360 {
4361 if (row_info->color_type == PNG_COLOR_TYPE_GRAY)
4362 {
4363 unsigned int gray = trans_color != NULL ? trans_color->gray : 0;
4364
4365 if (row_info->bit_depth < 8)
4366 {
4367 switch (row_info->bit_depth)
4368 {
4369 case 1:
4370 {
4371 gray = (gray & 0x01) * 0xff;
4372 sp = row + (png_size_t)((row_width - 1) >> 3);
4373 dp = row + (png_size_t)row_width - 1;
4374 shift = 7 - (int)((row_width + 7) & 0x07);
4375 for (i = 0; i < row_width; i++)
4376 {
4377 if ((*sp >> shift) & 0x01)
4378 *dp = 0xff;
4379
4380 else
4381 *dp = 0;
4382
4383 if (shift == 7)
4384 {
4385 shift = 0;
4386 sp--;
4387 }
4388
4389 else
4390 shift++;
4391
4392 dp--;
4393 }
4394 break;
4395 }
4396
4397 case 2:
4398 {
4399 gray = (gray & 0x03) * 0x55;
4400 sp = row + (png_size_t)((row_width - 1) >> 2);
4401 dp = row + (png_size_t)row_width - 1;
4402 shift = (int)((3 - ((row_width + 3) & 0x03)) << 1);
4403 for (i = 0; i < row_width; i++)
4404 {
4405 value = (*sp >> shift) & 0x03;
4406 *dp = (png_byte)(value | (value << 2) | (value << 4) |
4407 (value << 6));
4408 if (shift == 6)
4409 {
4410 shift = 0;
4411 sp--;
4412 }
4413
4414 else
4415 shift += 2;
4416
4417 dp--;
4418 }
4419 break;
4420 }
4421
4422 case 4:
4423 {
4424 gray = (gray & 0x0f) * 0x11;
4425 sp = row + (png_size_t)((row_width - 1) >> 1);
4426 dp = row + (png_size_t)row_width - 1;
4427 shift = (int)((1 - ((row_width + 1) & 0x01)) << 2);
4428 for (i = 0; i < row_width; i++)
4429 {
4430 value = (*sp >> shift) & 0x0f;
4431 *dp = (png_byte)(value | (value << 4));
4432 if (shift == 4)
4433 {
4434 shift = 0;
4435 sp--;
4436 }
4437
4438 else
4439 shift = 4;
4440
4441 dp--;
4442 }
4443 break;
4444 }
4445
4446 default:
4447 break;
4448 }
4449
4450 row_info->bit_depth = 8;
4451 row_info->pixel_depth = 8;
4452 row_info->rowbytes = row_width;
4453 }
4454
4455 if (trans_color != NULL)
4456 {
4457 if (row_info->bit_depth == 8)
4458 {
4459 gray = gray & 0xff;
4460 sp = row + (png_size_t)row_width - 1;
4461 dp = row + (png_size_t)(row_width << 1) - 1;
4462
4463 for (i = 0; i < row_width; i++)
4464 {
4465 if ((*sp & 0xffU) == gray)
4466 *dp-- = 0;
4467
4468 else
4469 *dp-- = 0xff;
4470
4471 *dp-- = *sp--;
4472 }
4473 }
4474
4475 else if (row_info->bit_depth == 16)
4476 {
4477 unsigned int gray_high = (gray >> 8) & 0xff;
4478 unsigned int gray_low = gray & 0xff;
4479 sp = row + row_info->rowbytes - 1;
4480 dp = row + (row_info->rowbytes << 1) - 1;
4481 for (i = 0; i < row_width; i++)
4482 {
4483 if ((*(sp - 1) & 0xffU) == gray_high &&
4484 (*(sp) & 0xffU) == gray_low)
4485 {
4486 *dp-- = 0;
4487 *dp-- = 0;
4488 }
4489
4490 else
4491 {
4492 *dp-- = 0xff;
4493 *dp-- = 0xff;
4494 }
4495
4496 *dp-- = *sp--;
4497 *dp-- = *sp--;
4498 }
4499 }
4500
4501 row_info->color_type = PNG_COLOR_TYPE_GRAY_ALPHA;
4502 row_info->channels = 2;
4503 row_info->pixel_depth = (png_byte)(row_info->bit_depth << 1);
4504 row_info->rowbytes = PNG_ROWBYTES(row_info->pixel_depth,
4505 row_width);
4506 }
4507 }
4508 else if (row_info->color_type == PNG_COLOR_TYPE_RGB &&
4509 trans_color != NULL)
4510 {
4511 if (row_info->bit_depth == 8)
4512 {
4513 png_byte red = (png_byte)(trans_color->red & 0xff);
4514 png_byte green = (png_byte)(trans_color->green & 0xff);
4515 png_byte blue = (png_byte)(trans_color->blue & 0xff);
4516 sp = row + (png_size_t)row_info->rowbytes - 1;
4517 dp = row + (png_size_t)(row_width << 2) - 1;
4518 for (i = 0; i < row_width; i++)
4519 {
4520 if (*(sp - 2) == red && *(sp - 1) == green && *(sp) == blue)
4521 *dp-- = 0;
4522
4523 else
4524 *dp-- = 0xff;
4525
4526 *dp-- = *sp--;
4527 *dp-- = *sp--;
4528 *dp-- = *sp--;
4529 }
4530 }
4531 else if (row_info->bit_depth == 16)
4532 {
4533 png_byte red_high = (png_byte)((trans_color->red >> 8) & 0xff);
4534 png_byte green_high = (png_byte)((trans_color->green >> 8) & 0xff);
4535 png_byte blue_high = (png_byte)((trans_color->blue >> 8) & 0xff);
4536 png_byte red_low = (png_byte)(trans_color->red & 0xff);
4537 png_byte green_low = (png_byte)(trans_color->green & 0xff);
4538 png_byte blue_low = (png_byte)(trans_color->blue & 0xff);
4539 sp = row + row_info->rowbytes - 1;
4540 dp = row + (png_size_t)(row_width << 3) - 1;
4541 for (i = 0; i < row_width; i++)
4542 {
4543 if (*(sp - 5) == red_high &&
4544 *(sp - 4) == red_low &&
4545 *(sp - 3) == green_high &&
4546 *(sp - 2) == green_low &&
4547 *(sp - 1) == blue_high &&
4548 *(sp ) == blue_low)
4549 {
4550 *dp-- = 0;
4551 *dp-- = 0;
4552 }
4553
4554 else
4555 {
4556 *dp-- = 0xff;
4557 *dp-- = 0xff;
4558 }
4559
4560 *dp-- = *sp--;
4561 *dp-- = *sp--;
4562 *dp-- = *sp--;
4563 *dp-- = *sp--;
4564 *dp-- = *sp--;
4565 *dp-- = *sp--;
4566 }
4567 }
4568 row_info->color_type = PNG_COLOR_TYPE_RGB_ALPHA;
4569 row_info->channels = 4;
4570 row_info->pixel_depth = (png_byte)(row_info->bit_depth << 2);
4571 row_info->rowbytes = PNG_ROWBYTES(row_info->pixel_depth, row_width);
4572 }
4573 }
4574}
4575#endif
4576
4577#ifdef PNG_READ_EXPAND_16_SUPPORTED
4578/* If the bit depth is 8 and the color type is not a palette type expand the
4579 * whole row to 16 bits. Has no effect otherwise.
4580 */
4581static void
4582png_do_expand_16(png_row_infop row_info, png_bytep row)
4583{
4584 if (row_info->bit_depth == 8 &&
4585 row_info->color_type != PNG_COLOR_TYPE_PALETTE)
4586 {
4587 /* The row have a sequence of bytes containing [0..255] and we need
4588 * to turn it into another row containing [0..65535], to do this we
4589 * calculate:
4590 *
4591 * (input / 255) * 65535
4592 *
4593 * Which happens to be exactly input * 257 and this can be achieved
4594 * simply by byte replication in place (copying backwards).
4595 */
4596 png_byte *sp = row + row_info->rowbytes; /* source, last byte + 1 */
4597 png_byte *dp = sp + row_info->rowbytes; /* destination, end + 1 */
4598 while (dp > sp)
4599 dp[-2] = dp[-1] = *--sp, dp -= 2;
4600
4601 row_info->rowbytes *= 2;
4602 row_info->bit_depth = 16;
4603 row_info->pixel_depth = (png_byte)(row_info->channels * 16);
4604 }
4605}
4606#endif
4607
4608#ifdef PNG_READ_QUANTIZE_SUPPORTED
4609static void
4610png_do_quantize(png_row_infop row_info, png_bytep row,
4611 png_const_bytep palette_lookup, png_const_bytep quantize_lookup)
4612{
4613 png_bytep sp, dp;
4614 png_uint_32 i;
4615 png_uint_32 row_width=row_info->width;
4616
4617 png_debug(1, "in png_do_quantize");
4618
4619 if (row_info->bit_depth == 8)
4620 {
4621 if (row_info->color_type == PNG_COLOR_TYPE_RGB && palette_lookup)
4622 {
4623 int r, g, b, p;
4624 sp = row;
4625 dp = row;
4626 for (i = 0; i < row_width; i++)
4627 {
4628 r = *sp++;
4629 g = *sp++;
4630 b = *sp++;
4631
4632 /* This looks real messy, but the compiler will reduce
4633 * it down to a reasonable formula. For example, with
4634 * 5 bits per color, we get:
4635 * p = (((r >> 3) & 0x1f) << 10) |
4636 * (((g >> 3) & 0x1f) << 5) |
4637 * ((b >> 3) & 0x1f);
4638 */
4639 p = (((r >> (8 - PNG_QUANTIZE_RED_BITS)) &
4640 ((1 << PNG_QUANTIZE_RED_BITS) - 1)) <<
4641 (PNG_QUANTIZE_GREEN_BITS + PNG_QUANTIZE_BLUE_BITS)) |
4642 (((g >> (8 - PNG_QUANTIZE_GREEN_BITS)) &
4643 ((1 << PNG_QUANTIZE_GREEN_BITS) - 1)) <<
4644 (PNG_QUANTIZE_BLUE_BITS)) |
4645 ((b >> (8 - PNG_QUANTIZE_BLUE_BITS)) &
4646 ((1 << PNG_QUANTIZE_BLUE_BITS) - 1));
4647
4648 *dp++ = palette_lookup[p];
4649 }
4650
4651 row_info->color_type = PNG_COLOR_TYPE_PALETTE;
4652 row_info->channels = 1;
4653 row_info->pixel_depth = row_info->bit_depth;
4654 row_info->rowbytes = PNG_ROWBYTES(row_info->pixel_depth, row_width);
4655 }
4656
4657 else if (row_info->color_type == PNG_COLOR_TYPE_RGB_ALPHA &&
4658 palette_lookup != NULL)
4659 {
4660 int r, g, b, p;
4661 sp = row;
4662 dp = row;
4663 for (i = 0; i < row_width; i++)
4664 {
4665 r = *sp++;
4666 g = *sp++;
4667 b = *sp++;
4668 sp++;
4669
4670 p = (((r >> (8 - PNG_QUANTIZE_RED_BITS)) &
4671 ((1 << PNG_QUANTIZE_RED_BITS) - 1)) <<
4672 (PNG_QUANTIZE_GREEN_BITS + PNG_QUANTIZE_BLUE_BITS)) |
4673 (((g >> (8 - PNG_QUANTIZE_GREEN_BITS)) &
4674 ((1 << PNG_QUANTIZE_GREEN_BITS) - 1)) <<
4675 (PNG_QUANTIZE_BLUE_BITS)) |
4676 ((b >> (8 - PNG_QUANTIZE_BLUE_BITS)) &
4677 ((1 << PNG_QUANTIZE_BLUE_BITS) - 1));
4678
4679 *dp++ = palette_lookup[p];
4680 }
4681
4682 row_info->color_type = PNG_COLOR_TYPE_PALETTE;
4683 row_info->channels = 1;
4684 row_info->pixel_depth = row_info->bit_depth;
4685 row_info->rowbytes = PNG_ROWBYTES(row_info->pixel_depth, row_width);
4686 }
4687
4688 else if (row_info->color_type == PNG_COLOR_TYPE_PALETTE &&
4689 quantize_lookup)
4690 {
4691 sp = row;
4692
4693 for (i = 0; i < row_width; i++, sp++)
4694 {
4695 *sp = quantize_lookup[*sp];
4696 }
4697 }
4698 }
4699}
4700#endif /* READ_QUANTIZE */
4701
4702/* Transform the row. The order of transformations is significant,
4703 * and is very touchy. If you add a transformation, take care to
4704 * decide how it fits in with the other transformations here.
4705 */
4706void /* PRIVATE */
4707png_do_read_transformations(png_structrp png_ptr, png_row_infop row_info)
4708{
4709 png_debug(1, "in png_do_read_transformations");
4710
4711 if (png_ptr->row_buf == NULL)
4712 {
4713 /* Prior to 1.5.4 this output row/pass where the NULL pointer is, but this
4714 * error is incredibly rare and incredibly easy to debug without this
4715 * information.
4716 */
4717 png_error(png_ptr, "NULL row buffer");
4718 }
4719
4720 /* The following is debugging; prior to 1.5.4 the code was never compiled in;
4721 * in 1.5.4 PNG_FLAG_DETECT_UNINITIALIZED was added and the macro
4722 * PNG_WARN_UNINITIALIZED_ROW removed. In 1.6 the new flag is set only for
4723 * all transformations, however in practice the ROW_INIT always gets done on
4724 * demand, if necessary.
4725 */
4726 if ((png_ptr->flags & PNG_FLAG_DETECT_UNINITIALIZED) != 0 &&
4727 (png_ptr->flags & PNG_FLAG_ROW_INIT) == 0)
4728 {
4729 /* Application has failed to call either png_read_start_image() or
4730 * png_read_update_info() after setting transforms that expand pixels.
4731 * This check added to libpng-1.2.19 (but not enabled until 1.5.4).
4732 */
4733 png_error(png_ptr, "Uninitialized row");
4734 }
4735
4736#ifdef PNG_READ_EXPAND_SUPPORTED
4737 if ((png_ptr->transformations & PNG_EXPAND) != 0)
4738 {
4739 if (row_info->color_type == PNG_COLOR_TYPE_PALETTE)
4740 {
4741 png_do_expand_palette(row_info, png_ptr->row_buf + 1,
4742 png_ptr->palette, png_ptr->trans_alpha, png_ptr->num_trans);
4743 }
4744
4745 else
4746 {
4747 if (png_ptr->num_trans != 0 &&
4748 (png_ptr->transformations & PNG_EXPAND_tRNS) != 0)
4749 png_do_expand(row_info, png_ptr->row_buf + 1,
4750 &(png_ptr->trans_color));
4751
4752 else
4753 png_do_expand(row_info, png_ptr->row_buf + 1,
4754 NULL);
4755 }
4756 }
4757#endif
4758
4759#ifdef PNG_READ_STRIP_ALPHA_SUPPORTED
4760 if ((png_ptr->transformations & PNG_STRIP_ALPHA) != 0 &&
4761 (png_ptr->transformations & PNG_COMPOSE) == 0 &&
4762 (row_info->color_type == PNG_COLOR_TYPE_RGB_ALPHA ||
4763 row_info->color_type == PNG_COLOR_TYPE_GRAY_ALPHA))
4764 png_do_strip_channel(row_info, png_ptr->row_buf + 1,
4765 0 /* at_start == false, because SWAP_ALPHA happens later */);
4766#endif
4767
4768#ifdef PNG_READ_RGB_TO_GRAY_SUPPORTED
4769 if ((png_ptr->transformations & PNG_RGB_TO_GRAY) != 0)
4770 {
4771 int rgb_error =
4772 png_do_rgb_to_gray(png_ptr, row_info,
4773 png_ptr->row_buf + 1);
4774
4775 if (rgb_error != 0)
4776 {
4777 png_ptr->rgb_to_gray_status=1;
4778 if ((png_ptr->transformations & PNG_RGB_TO_GRAY) ==
4779 PNG_RGB_TO_GRAY_WARN)
4780 png_warning(png_ptr, "png_do_rgb_to_gray found nongray pixel");
4781
4782 if ((png_ptr->transformations & PNG_RGB_TO_GRAY) ==
4783 PNG_RGB_TO_GRAY_ERR)
4784 png_error(png_ptr, "png_do_rgb_to_gray found nongray pixel");
4785 }
4786 }
4787#endif
4788
4789/* From Andreas Dilger e-mail to png-implement, 26 March 1998:
4790 *
4791 * In most cases, the "simple transparency" should be done prior to doing
4792 * gray-to-RGB, or you will have to test 3x as many bytes to check if a
4793 * pixel is transparent. You would also need to make sure that the
4794 * transparency information is upgraded to RGB.
4795 *
4796 * To summarize, the current flow is:
4797 * - Gray + simple transparency -> compare 1 or 2 gray bytes and composite
4798 * with background "in place" if transparent,
4799 * convert to RGB if necessary
4800 * - Gray + alpha -> composite with gray background and remove alpha bytes,
4801 * convert to RGB if necessary
4802 *
4803 * To support RGB backgrounds for gray images we need:
4804 * - Gray + simple transparency -> convert to RGB + simple transparency,
4805 * compare 3 or 6 bytes and composite with
4806 * background "in place" if transparent
4807 * (3x compare/pixel compared to doing
4808 * composite with gray bkgrnd)
4809 * - Gray + alpha -> convert to RGB + alpha, composite with background and
4810 * remove alpha bytes (3x float
4811 * operations/pixel compared with composite
4812 * on gray background)
4813 *
4814 * Greg's change will do this. The reason it wasn't done before is for
4815 * performance, as this increases the per-pixel operations. If we would check
4816 * in advance if the background was gray or RGB, and position the gray-to-RGB
4817 * transform appropriately, then it would save a lot of work/time.
4818 */
4819
4820#ifdef PNG_READ_GRAY_TO_RGB_SUPPORTED
4821 /* If gray -> RGB, do so now only if background is non-gray; else do later
4822 * for performance reasons
4823 */
4824 if ((png_ptr->transformations & PNG_GRAY_TO_RGB) != 0 &&
4825 (png_ptr->mode & PNG_BACKGROUND_IS_GRAY) == 0)
4826 png_do_gray_to_rgb(row_info, png_ptr->row_buf + 1);
4827#endif
4828
4829#if defined(PNG_READ_BACKGROUND_SUPPORTED) ||\
4830 defined(PNG_READ_ALPHA_MODE_SUPPORTED)
4831 if ((png_ptr->transformations & PNG_COMPOSE) != 0)
4832 png_do_compose(row_info, png_ptr->row_buf + 1, png_ptr);
4833#endif
4834
4835#ifdef PNG_READ_GAMMA_SUPPORTED
4836 if ((png_ptr->transformations & PNG_GAMMA) != 0 &&
4837#ifdef PNG_READ_RGB_TO_GRAY_SUPPORTED
4838 /* Because RGB_TO_GRAY does the gamma transform. */
4839 (png_ptr->transformations & PNG_RGB_TO_GRAY) == 0 &&
4840#endif
4841#if defined(PNG_READ_BACKGROUND_SUPPORTED) ||\
4842 defined(PNG_READ_ALPHA_MODE_SUPPORTED)
4843 /* Because PNG_COMPOSE does the gamma transform if there is something to
4844 * do (if there is an alpha channel or transparency.)
4845 */
4846 !((png_ptr->transformations & PNG_COMPOSE) != 0 &&
4847 ((png_ptr->num_trans != 0) ||
4848 (png_ptr->color_type & PNG_COLOR_MASK_ALPHA) != 0)) &&
4849#endif
4850 /* Because png_init_read_transformations transforms the palette, unless
4851 * RGB_TO_GRAY will do the transform.
4852 */
4853 (png_ptr->color_type != PNG_COLOR_TYPE_PALETTE))
4854 png_do_gamma(row_info, png_ptr->row_buf + 1, png_ptr);
4855#endif
4856
4857#ifdef PNG_READ_STRIP_ALPHA_SUPPORTED
4858 if ((png_ptr->transformations & PNG_STRIP_ALPHA) != 0 &&
4859 (png_ptr->transformations & PNG_COMPOSE) != 0 &&
4860 (row_info->color_type == PNG_COLOR_TYPE_RGB_ALPHA ||
4861 row_info->color_type == PNG_COLOR_TYPE_GRAY_ALPHA))
4862 png_do_strip_channel(row_info, png_ptr->row_buf + 1,
4863 0 /* at_start == false, because SWAP_ALPHA happens later */);
4864#endif
4865
4866#ifdef PNG_READ_ALPHA_MODE_SUPPORTED
4867 if ((png_ptr->transformations & PNG_ENCODE_ALPHA) != 0 &&
4868 (row_info->color_type & PNG_COLOR_MASK_ALPHA) != 0)
4869 png_do_encode_alpha(row_info, png_ptr->row_buf + 1, png_ptr);
4870#endif
4871
4872#ifdef PNG_READ_SCALE_16_TO_8_SUPPORTED
4873 if ((png_ptr->transformations & PNG_SCALE_16_TO_8) != 0)
4874 png_do_scale_16_to_8(row_info, png_ptr->row_buf + 1);
4875#endif
4876
4877#ifdef PNG_READ_STRIP_16_TO_8_SUPPORTED
4878 /* There is no harm in doing both of these because only one has any effect,
4879 * by putting the 'scale' option first if the app asks for scale (either by
4880 * calling the API or in a TRANSFORM flag) this is what happens.
4881 */
4882 if ((png_ptr->transformations & PNG_16_TO_8) != 0)
4883 png_do_chop(row_info, png_ptr->row_buf + 1);
4884#endif
4885
4886#ifdef PNG_READ_QUANTIZE_SUPPORTED
4887 if ((png_ptr->transformations & PNG_QUANTIZE) != 0)
4888 {
4889 png_do_quantize(row_info, png_ptr->row_buf + 1,
4890 png_ptr->palette_lookup, png_ptr->quantize_index);
4891
4892 if (row_info->rowbytes == 0)
489