PageRenderTime 62ms CodeModel.GetById 22ms RepoModel.GetById 0ms app.codeStats 0ms

/src/FreeImage/Source/LibTIFF/tif_getimage.c

https://bitbucket.org/cabalistic/ogredeps/
C | 2554 lines | 2069 code | 183 blank | 302 comment | 389 complexity | 01a416068600961142938173e5984966 MD5 | raw file
Possible License(s): LGPL-3.0, BSD-3-Clause, CPL-1.0, Unlicense, GPL-2.0, GPL-3.0, LGPL-2.0, MPL-2.0-no-copyleft-exception, BSD-2-Clause, LGPL-2.1 
    

  1. /* $Id: tif_getimage.c,v 1.39 2011/04/10 17:14:09 drolon Exp $ */
    2. 

  2. 

  3. /*
    4. 

  4.  * Copyright (c) 1991-1997 Sam Leffler
    5. 

  5.  * Copyright (c) 1991-1997 Silicon Graphics, Inc.
    6. 

  6.  *
    7. 

  7.  * Permission to use, copy, modify, distribute, and sell this software and 
    8. 

  8.  * its documentation for any purpose is hereby granted without fee, provided
    9. 

  9.  * that (i) the above copyright notices and this permission notice appear in
    10. 

  10.  * all copies of the software and related documentation, and (ii) the names of
    11. 

  11.  * Sam Leffler and Silicon Graphics may not be used in any advertising or
    12. 

  12.  * publicity relating to the software without the specific, prior written
    13. 

  13.  * permission of Sam Leffler and Silicon Graphics.
    14. 

  14.  * 
    15. 

  15.  * THE SOFTWARE IS PROVIDED "AS-IS" AND WITHOUT WARRANTY OF ANY KIND, 
    16. 

  16.  * EXPRESS, IMPLIED OR OTHERWISE, INCLUDING WITHOUT LIMITATION, ANY 
    17. 

  17.  * WARRANTY OF MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE.  
    18. 

  18.  * 
    19. 

  19.  * IN NO EVENT SHALL SAM LEFFLER OR SILICON GRAPHICS BE LIABLE FOR
    20. 

  20.  * ANY SPECIAL, INCIDENTAL, INDIRECT OR CONSEQUENTIAL DAMAGES OF ANY KIND,
    21. 

  21.  * OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS,
    22. 

  22.  * WHETHER OR NOT ADVISED OF THE POSSIBILITY OF DAMAGE, AND ON ANY THEORY OF 
    23. 

  23.  * LIABILITY, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE 
    24. 

  24.  * OF THIS SOFTWARE.
    25. 

  25.  */
    26. 

  26. 

  27. /*
    28. 

  28.  * TIFF Library
    29. 

  29.  *
    30. 

  30.  * Read and return a packed RGBA image.
    31. 

  31.  */
    32. 

  32. #include "tiffiop.h"
    33. 

  33. #include <stdio.h>
    34. 

  34. 

  35. static int gtTileContig(TIFFRGBAImage*, uint32*, uint32, uint32);
    36. 

  36. static int gtTileSeparate(TIFFRGBAImage*, uint32*, uint32, uint32);
    37. 

  37. static int gtStripContig(TIFFRGBAImage*, uint32*, uint32, uint32);
    38. 

  38. static int gtStripSeparate(TIFFRGBAImage*, uint32*, uint32, uint32);
    39. 

  39. static int PickContigCase(TIFFRGBAImage*);
    40. 

  40. static int PickSeparateCase(TIFFRGBAImage*);
    41. 

  41. static const char photoTag[] = "PhotometricInterpretation";
    42. 

  42. 

  43. /* 
    44. 

  44.  * Helper constants used in Orientation tag handling
    45. 

  45.  */
    46. 

  46. #define FLIP_VERTICALLY 0x01
    47. 

  47. #define FLIP_HORIZONTALLY 0x02
    48. 

  48. 

  49. /*
    50. 

  50.  * Color conversion constants. We will define display types here.
    51. 

  51.  */
    52. 

  52. 

  53. TIFFDisplay display_sRGB = {
    54. 

  54. 	{			/* XYZ -> luminance matrix */
    55. 

  55. 		{  3.2410F, -1.5374F, -0.4986F },
    56. 

  56. 		{  -0.9692F, 1.8760F, 0.0416F },
    57. 

  57. 		{  0.0556F, -0.2040F, 1.0570F }
    58. 

  58. 	},	
    59. 

  59. 	100.0F, 100.0F, 100.0F,	/* Light o/p for reference white */
    60. 

  60. 	255, 255, 255,		/* Pixel values for ref. white */
    61. 

  61. 	1.0F, 1.0F, 1.0F,	/* Residual light o/p for black pixel */
    62. 

  62. 	2.4F, 2.4F, 2.4F,	/* Gamma values for the three guns */
    63. 

  63. };
    64. 

  64. 

  65. /*
    66. 

  66.  * Check the image to see if TIFFReadRGBAImage can deal with it.
    67. 

  67.  * 1/0 is returned according to whether or not the image can
    68. 

  68.  * be handled.  If 0 is returned, emsg contains the reason
    69. 

  69.  * why it is being rejected.
    70. 

  70.  */
    71. 

  71. int
    72. 

  72. TIFFRGBAImageOK(TIFF* tif, char emsg[1024])
    73. 

  73. {
    74. 

  74. 	TIFFDirectory* td = &tif->tif_dir;
    75. 

  75. 	uint16 photometric;
    76. 

  76. 	int colorchannels;
    77. 

  77. 

  78. 	if (!tif->tif_decodestatus) {
    79. 

  79. 		sprintf(emsg, "Sorry, requested compression method is not configured");
    80. 

  80. 		return (0);
    81. 

  81. 	}
    82. 

  82. 	switch (td->td_bitspersample) {
    83. 

  83. 		case 1:
    84. 

  84. 		case 2:
    85. 

  85. 		case 4:
    86. 

  86. 		case 8:
    87. 

  87. 		case 16:
    88. 

  88. 			break;
    89. 

  89. 		default:
    90. 

  90. 			sprintf(emsg, "Sorry, can not handle images with %d-bit samples",
    91. 

  91. 			    td->td_bitspersample);
    92. 

  92. 			return (0);
    93. 

  93. 	}
    94. 

  94. 	colorchannels = td->td_samplesperpixel - td->td_extrasamples;
    95. 

  95. 	if (!TIFFGetField(tif, TIFFTAG_PHOTOMETRIC, &photometric)) {
    96. 

  96. 		switch (colorchannels) {
    97. 

  97. 			case 1:
    98. 

  98. 				photometric = PHOTOMETRIC_MINISBLACK;
    99. 

  99. 				break;
    100. 

  100. 			case 3:
    101. 

  101. 				photometric = PHOTOMETRIC_RGB;
    102. 

  102. 				break;
    103. 

  103. 			default:
    104. 

  104. 				sprintf(emsg, "Missing needed %s tag", photoTag);
    105. 

  105. 				return (0);
    106. 

  106. 		}
    107. 

  107. 	}
    108. 

  108. 	switch (photometric) {
    109. 

  109. 		case PHOTOMETRIC_MINISWHITE:
    110. 

  110. 		case PHOTOMETRIC_MINISBLACK:
    111. 

  111. 		case PHOTOMETRIC_PALETTE:
    112. 

  112. 			if (td->td_planarconfig == PLANARCONFIG_CONTIG
    113. 

  113. 			    && td->td_samplesperpixel != 1
    114. 

  114. 			    && td->td_bitspersample < 8 ) {
    115. 

  115. 				sprintf(emsg,
    116. 

  116. 				    "Sorry, can not handle contiguous data with %s=%d, "
    117. 

  117. 				    "and %s=%d and Bits/Sample=%d",
    118. 

  118. 				    photoTag, photometric,
    119. 

  119. 				    "Samples/pixel", td->td_samplesperpixel,
    120. 

  120. 				    td->td_bitspersample);
    121. 

  121. 				return (0);
    122. 

  122. 			}
    123. 

  123. 			/*
    124. 

  124. 			 * We should likely validate that any extra samples are either
    125. 

  125. 			 * to be ignored, or are alpha, and if alpha we should try to use
    126. 

  126. 			 * them.  But for now we won't bother with this.
    127. 

  127. 			*/
    128. 

  128. 			break;
    129. 

  129. 		case PHOTOMETRIC_YCBCR:
    130. 

  130. 			/*
    131. 

  131. 			 * TODO: if at all meaningful and useful, make more complete
    132. 

  132. 			 * support check here, or better still, refactor to let supporting
    133. 

  133. 			 * code decide whether there is support and what meaningfull
    134. 

  134. 			 * error to return
    135. 

  135. 			 */
    136. 

  136. 			break;
    137. 

  137. 		case PHOTOMETRIC_RGB:
    138. 

  138. 			if (colorchannels < 3) {
    139. 

  139. 				sprintf(emsg, "Sorry, can not handle RGB image with %s=%d",
    140. 

  140. 				    "Color channels", colorchannels);
    141. 

  141. 				return (0);
    142. 

  142. 			}
    143. 

  143. 			break;
    144. 

  144. 		case PHOTOMETRIC_SEPARATED:
    145. 

  145. 			{
    146. 

  146. 				uint16 inkset;
    147. 

  147. 				TIFFGetFieldDefaulted(tif, TIFFTAG_INKSET, &inkset);
    148. 

  148. 				if (inkset != INKSET_CMYK) {
    149. 

  149. 					sprintf(emsg,
    150. 

  150. 					    "Sorry, can not handle separated image with %s=%d",
    151. 

  151. 					    "InkSet", inkset);
    152. 

  152. 					return 0;
    153. 

  153. 				}
    154. 

  154. 				if (td->td_samplesperpixel < 4) {
    155. 

  155. 					sprintf(emsg,
    156. 

  156. 					    "Sorry, can not handle separated image with %s=%d",
    157. 

  157. 					    "Samples/pixel", td->td_samplesperpixel);
    158. 

  158. 					return 0;
    159. 

  159. 				}
    160. 

  160. 				break;
    161. 

  161. 			}
    162. 

  162. 		case PHOTOMETRIC_LOGL:
    163. 

  163. 			if (td->td_compression != COMPRESSION_SGILOG) {
    164. 

  164. 				sprintf(emsg, "Sorry, LogL data must have %s=%d",
    165. 

  165. 				    "Compression", COMPRESSION_SGILOG);
    166. 

  166. 				return (0);
    167. 

  167. 			}
    168. 

  168. 			break;
    169. 

  169. 		case PHOTOMETRIC_LOGLUV:
    170. 

  170. 			if (td->td_compression != COMPRESSION_SGILOG &&
    171. 

  171. 			    td->td_compression != COMPRESSION_SGILOG24) {
    172. 

  172. 				sprintf(emsg, "Sorry, LogLuv data must have %s=%d or %d",
    173. 

  173. 				    "Compression", COMPRESSION_SGILOG, COMPRESSION_SGILOG24);
    174. 

  174. 				return (0);
    175. 

  175. 			}
    176. 

  176. 			if (td->td_planarconfig != PLANARCONFIG_CONTIG) {
    177. 

  177. 				sprintf(emsg, "Sorry, can not handle LogLuv images with %s=%d",
    178. 

  178. 				    "Planarconfiguration", td->td_planarconfig);
    179. 

  179. 				return (0);
    180. 

  180. 			}
    181. 

  181. 			break;
    182. 

  182. 		case PHOTOMETRIC_CIELAB:
    183. 

  183. 			break;
    184. 

  184. 		default:
    185. 

  185. 			sprintf(emsg, "Sorry, can not handle image with %s=%d",
    186. 

  186. 			    photoTag, photometric);
    187. 

  187. 			return (0);
    188. 

  188. 	}
    189. 

  189. 	return (1);
    190. 

  190. }
    191. 

  191. 

  192. void
    193. 

  193. TIFFRGBAImageEnd(TIFFRGBAImage* img)
    194. 

  194. {
    195. 

  195. 	if (img->Map)
    196. 

  196. 		_TIFFfree(img->Map), img->Map = NULL;
    197. 

  197. 	if (img->BWmap)
    198. 

  198. 		_TIFFfree(img->BWmap), img->BWmap = NULL;
    199. 

  199. 	if (img->PALmap)
    200. 

  200. 		_TIFFfree(img->PALmap), img->PALmap = NULL;
    201. 

  201. 	if (img->ycbcr)
    202. 

  202. 		_TIFFfree(img->ycbcr), img->ycbcr = NULL;
    203. 

  203. 	if (img->cielab)
    204. 

  204. 		_TIFFfree(img->cielab), img->cielab = NULL;
    205. 

  205. 	if( img->redcmap ) {
    206. 

  206. 		_TIFFfree( img->redcmap );
    207. 

  207. 		_TIFFfree( img->greencmap );
    208. 

  208. 		_TIFFfree( img->bluecmap );
    209. 

  209. 	}
    210. 

  210. }
    211. 

  211. 

  212. static int
    213. 

  213. isCCITTCompression(TIFF* tif)
    214. 

  214. {
    215. 

  215.     uint16 compress;
    216. 

  216.     TIFFGetField(tif, TIFFTAG_COMPRESSION, &compress);
    217. 

  217.     return (compress == COMPRESSION_CCITTFAX3 ||
    218. 

  218. 	    compress == COMPRESSION_CCITTFAX4 ||
    219. 

  219. 	    compress == COMPRESSION_CCITTRLE ||
    220. 

  220. 	    compress == COMPRESSION_CCITTRLEW);
    221. 

  221. }
    222. 

  222. 

  223. int
    224. 

  224. TIFFRGBAImageBegin(TIFFRGBAImage* img, TIFF* tif, int stop, char emsg[1024])
    225. 

  225. {
    226. 

  226. 	uint16* sampleinfo;
    227. 

  227. 	uint16 extrasamples;
    228. 

  228. 	uint16 planarconfig;
    229. 

  229. 	uint16 compress;
    230. 

  230. 	int colorchannels;
    231. 

  231. 	uint16 *red_orig, *green_orig, *blue_orig;
    232. 

  232. 	int n_color;
    233. 

  233. 

  234. 	/* Initialize to normal values */
    235. 

  235. 	img->row_offset = 0;
    236. 

  236. 	img->col_offset = 0;
    237. 

  237. 	img->redcmap = NULL;
    238. 

  238. 	img->greencmap = NULL;
    239. 

  239. 	img->bluecmap = NULL;
    240. 

  240. 	img->req_orientation = ORIENTATION_BOTLEFT;     /* It is the default */
    241. 

  241. 

  242. 	img->tif = tif;
    243. 

  243. 	img->stoponerr = stop;
    244. 

  244. 	TIFFGetFieldDefaulted(tif, TIFFTAG_BITSPERSAMPLE, &img->bitspersample);
    245. 

  245. 	switch (img->bitspersample) {
    246. 

  246. 		case 1:
    247. 

  247. 		case 2:
    248. 

  248. 		case 4:
    249. 

  249. 		case 8:
    250. 

  250. 		case 16:
    251. 

  251. 			break;
    252. 

  252. 		default:
    253. 

  253. 			sprintf(emsg, "Sorry, can not handle images with %d-bit samples",
    254. 

  254. 			    img->bitspersample);
    255. 

  255. 			return (0);
    256. 

  256. 	}
    257. 

  257. 	img->alpha = 0;
    258. 

  258. 	TIFFGetFieldDefaulted(tif, TIFFTAG_SAMPLESPERPIXEL, &img->samplesperpixel);
    259. 

  259. 	TIFFGetFieldDefaulted(tif, TIFFTAG_EXTRASAMPLES,
    260. 

  260. 	    &extrasamples, &sampleinfo);
    261. 

  261. 	if (extrasamples >= 1)
    262. 

  262. 	{
    263. 

  263. 		switch (sampleinfo[0]) {
    264. 

  264. 			case EXTRASAMPLE_UNSPECIFIED:          /* Workaround for some images without */
    265. 

  265. 				if (img->samplesperpixel > 3)  /* correct info about alpha channel */
    266. 

  266. 					img->alpha = EXTRASAMPLE_ASSOCALPHA;
    267. 

  267. 				break;
    268. 

  268. 			case EXTRASAMPLE_ASSOCALPHA:           /* data is pre-multiplied */
    269. 

  269. 			case EXTRASAMPLE_UNASSALPHA:           /* data is not pre-multiplied */
    270. 

  270. 				img->alpha = sampleinfo[0];
    271. 

  271. 				break;
    272. 

  272. 		}
    273. 

  273. 	}
    274. 

  274. 

  275. #ifdef DEFAULT_EXTRASAMPLE_AS_ALPHA
    276. 

  276. 	if( !TIFFGetField(tif, TIFFTAG_PHOTOMETRIC, &img->photometric))
    277. 

  277. 		img->photometric = PHOTOMETRIC_MINISWHITE;
    278. 

  278. 

  279. 	if( extrasamples == 0
    280. 

  280. 	    && img->samplesperpixel == 4
    281. 

  281. 	    && img->photometric == PHOTOMETRIC_RGB )
    282. 

  282. 	{
    283. 

  283. 		img->alpha = EXTRASAMPLE_ASSOCALPHA;
    284. 

  284. 		extrasamples = 1;
    285. 

  285. 	}
    286. 

  286. #endif
    287. 

  287. 

  288. 	colorchannels = img->samplesperpixel - extrasamples;
    289. 

  289. 	TIFFGetFieldDefaulted(tif, TIFFTAG_COMPRESSION, &compress);
    290. 

  290. 	TIFFGetFieldDefaulted(tif, TIFFTAG_PLANARCONFIG, &planarconfig);
    291. 

  291. 	if (!TIFFGetField(tif, TIFFTAG_PHOTOMETRIC, &img->photometric)) {
    292. 

  292. 		switch (colorchannels) {
    293. 

  293. 			case 1:
    294. 

  294. 				if (isCCITTCompression(tif))
    295. 

  295. 					img->photometric = PHOTOMETRIC_MINISWHITE;
    296. 

  296. 				else
    297. 

  297. 					img->photometric = PHOTOMETRIC_MINISBLACK;
    298. 

  298. 				break;
    299. 

  299. 			case 3:
    300. 

  300. 				img->photometric = PHOTOMETRIC_RGB;
    301. 

  301. 				break;
    302. 

  302. 			default:
    303. 

  303. 				sprintf(emsg, "Missing needed %s tag", photoTag);
    304. 

  304. 				return (0);
    305. 

  305. 		}
    306. 

  306. 	}
    307. 

  307. 	switch (img->photometric) {
    308. 

  308. 		case PHOTOMETRIC_PALETTE:
    309. 

  309. 			if (!TIFFGetField(tif, TIFFTAG_COLORMAP,
    310. 

  310. 			    &red_orig, &green_orig, &blue_orig)) {
    311. 

  311. 				sprintf(emsg, "Missing required \"Colormap\" tag");
    312. 

  312. 				return (0);
    313. 

  313. 			}
    314. 

  314. 

  315. 			/* copy the colormaps so we can modify them */
    316. 

  316. 			n_color = (1L << img->bitspersample);
    317. 

  317. 			img->redcmap = (uint16 *) _TIFFmalloc(sizeof(uint16)*n_color);
    318. 

  318. 			img->greencmap = (uint16 *) _TIFFmalloc(sizeof(uint16)*n_color);
    319. 

  319. 			img->bluecmap = (uint16 *) _TIFFmalloc(sizeof(uint16)*n_color);
    320. 

  320. 			if( !img->redcmap || !img->greencmap || !img->bluecmap ) {
    321. 

  321. 				sprintf(emsg, "Out of memory for colormap copy");
    322. 

  322. 				return (0);
    323. 

  323. 			}
    324. 

  324. 

  325. 			_TIFFmemcpy( img->redcmap, red_orig, n_color * 2 );
    326. 

  326. 			_TIFFmemcpy( img->greencmap, green_orig, n_color * 2 );
    327. 

  327. 			_TIFFmemcpy( img->bluecmap, blue_orig, n_color * 2 );
    328. 

  328. 

  329. 			/* fall thru... */
    330. 

  330. 		case PHOTOMETRIC_MINISWHITE:
    331. 

  331. 		case PHOTOMETRIC_MINISBLACK:
    332. 

  332. 			if (planarconfig == PLANARCONFIG_CONTIG
    333. 

  333. 			    && img->samplesperpixel != 1
    334. 

  334. 			    && img->bitspersample < 8 ) {
    335. 

  335. 				sprintf(emsg,
    336. 

  336. 				    "Sorry, can not handle contiguous data with %s=%d, "
    337. 

  337. 				    "and %s=%d and Bits/Sample=%d",
    338. 

  338. 				    photoTag, img->photometric,
    339. 

  339. 				    "Samples/pixel", img->samplesperpixel,
    340. 

  340. 				    img->bitspersample);
    341. 

  341. 				return (0);
    342. 

  342. 			}
    343. 

  343. 			break;
    344. 

  344. 		case PHOTOMETRIC_YCBCR:
    345. 

  345. 			/* It would probably be nice to have a reality check here. */
    346. 

  346. 			if (planarconfig == PLANARCONFIG_CONTIG)
    347. 

  347. 				/* can rely on libjpeg to convert to RGB */
    348. 

  348. 				/* XXX should restore current state on exit */
    349. 

  349. 				switch (compress) {
    350. 

  350. 					case COMPRESSION_JPEG:
    351. 

  351. 						/*
    352. 

  352. 						 * TODO: when complete tests verify complete desubsampling
    353. 

  353. 						 * and YCbCr handling, remove use of TIFFTAG_JPEGCOLORMODE in
    354. 

  354. 						 * favor of tif_getimage.c native handling
    355. 

  355. 						 */
    356. 

  356. 						TIFFSetField(tif, TIFFTAG_JPEGCOLORMODE, JPEGCOLORMODE_RGB);
    357. 

  357. 						img->photometric = PHOTOMETRIC_RGB;
    358. 

  358. 						break;
    359. 

  359. 					default:
    360. 

  360. 						/* do nothing */;
    361. 

  361. 						break;
    362. 

  362. 				}
    363. 

  363. 			/*
    364. 

  364. 			 * TODO: if at all meaningful and useful, make more complete
    365. 

  365. 			 * support check here, or better still, refactor to let supporting
    366. 

  366. 			 * code decide whether there is support and what meaningfull
    367. 

  367. 			 * error to return
    368. 

  368. 			 */
    369. 

  369. 			break;
    370. 

  370. 		case PHOTOMETRIC_RGB:
    371. 

  371. 			if (colorchannels < 3) {
    372. 

  372. 				sprintf(emsg, "Sorry, can not handle RGB image with %s=%d",
    373. 

  373. 				    "Color channels", colorchannels);
    374. 

  374. 				return (0);
    375. 

  375. 			}
    376. 

  376. 			break;
    377. 

  377. 		case PHOTOMETRIC_SEPARATED:
    378. 

  378. 			{
    379. 

  379. 				uint16 inkset;
    380. 

  380. 				TIFFGetFieldDefaulted(tif, TIFFTAG_INKSET, &inkset);
    381. 

  381. 				if (inkset != INKSET_CMYK) {
    382. 

  382. 					sprintf(emsg, "Sorry, can not handle separated image with %s=%d",
    383. 

  383. 					    "InkSet", inkset);
    384. 

  384. 					return (0);
    385. 

  385. 				}
    386. 

  386. 				if (img->samplesperpixel < 4) {
    387. 

  387. 					sprintf(emsg, "Sorry, can not handle separated image with %s=%d",
    388. 

  388. 					    "Samples/pixel", img->samplesperpixel);
    389. 

  389. 					return (0);
    390. 

  390. 				}
    391. 

  391. 			}
    392. 

  392. 			break;
    393. 

  393. 		case PHOTOMETRIC_LOGL:
    394. 

  394. 			if (compress != COMPRESSION_SGILOG) {
    395. 

  395. 				sprintf(emsg, "Sorry, LogL data must have %s=%d",
    396. 

  396. 				    "Compression", COMPRESSION_SGILOG);
    397. 

  397. 				return (0);
    398. 

  398. 			}
    399. 

  399. 			TIFFSetField(tif, TIFFTAG_SGILOGDATAFMT, SGILOGDATAFMT_8BIT);
    400. 

  400. 			img->photometric = PHOTOMETRIC_MINISBLACK;	/* little white lie */
    401. 

  401. 			img->bitspersample = 8;
    402. 

  402. 			break;
    403. 

  403. 		case PHOTOMETRIC_LOGLUV:
    404. 

  404. 			if (compress != COMPRESSION_SGILOG && compress != COMPRESSION_SGILOG24) {
    405. 

  405. 				sprintf(emsg, "Sorry, LogLuv data must have %s=%d or %d",
    406. 

  406. 				    "Compression", COMPRESSION_SGILOG, COMPRESSION_SGILOG24);
    407. 

  407. 				return (0);
    408. 

  408. 			}
    409. 

  409. 			if (planarconfig != PLANARCONFIG_CONTIG) {
    410. 

  410. 				sprintf(emsg, "Sorry, can not handle LogLuv images with %s=%d",
    411. 

  411. 				    "Planarconfiguration", planarconfig);
    412. 

  412. 				return (0);
    413. 

  413. 			}
    414. 

  414. 			TIFFSetField(tif, TIFFTAG_SGILOGDATAFMT, SGILOGDATAFMT_8BIT);
    415. 

  415. 			img->photometric = PHOTOMETRIC_RGB;		/* little white lie */
    416. 

  416. 			img->bitspersample = 8;
    417. 

  417. 			break;
    418. 

  418. 		case PHOTOMETRIC_CIELAB:
    419. 

  419. 			break;
    420. 

  420. 		default:
    421. 

  421. 			sprintf(emsg, "Sorry, can not handle image with %s=%d",
    422. 

  422. 			    photoTag, img->photometric);
    423. 

  423. 			return (0);
    424. 

  424. 	}
    425. 

  425. 	img->Map = NULL;
    426. 

  426. 	img->BWmap = NULL;
    427. 

  427. 	img->PALmap = NULL;
    428. 

  428. 	img->ycbcr = NULL;
    429. 

  429. 	img->cielab = NULL;
    430. 

  430. 	TIFFGetField(tif, TIFFTAG_IMAGEWIDTH, &img->width);
    431. 

  431. 	TIFFGetField(tif, TIFFTAG_IMAGELENGTH, &img->height);
    432. 

  432. 	TIFFGetFieldDefaulted(tif, TIFFTAG_ORIENTATION, &img->orientation);
    433. 

  433. 	img->isContig =
    434. 

  434. 	    !(planarconfig == PLANARCONFIG_SEPARATE && colorchannels > 1);
    435. 

  435. 	if (img->isContig) {
    436. 

  436. 		if (!PickContigCase(img)) {
    437. 

  437. 			sprintf(emsg, "Sorry, can not handle image");
    438. 

  438. 			return 0;
    439. 

  439. 		}
    440. 

  440. 	} else {
    441. 

  441. 		if (!PickSeparateCase(img)) {
    442. 

  442. 			sprintf(emsg, "Sorry, can not handle image");
    443. 

  443. 			return 0;
    444. 

  444. 		}
    445. 

  445. 	}
    446. 

  446. 	return 1;
    447. 

  447. }
    448. 

  448. 

  449. int
    450. 

  450. TIFFRGBAImageGet(TIFFRGBAImage* img, uint32* raster, uint32 w, uint32 h)
    451. 

  451. {
    452. 

  452.     if (img->get == NULL) {
    453. 

  453. 		TIFFErrorExt(img->tif->tif_clientdata, TIFFFileName(img->tif), "No \"get\" routine setup");
    454. 

  454. 		return (0);
    455. 

  455. 	}
    456. 

  456. 	if (img->put.any == NULL) {
    457. 

  457. 		TIFFErrorExt(img->tif->tif_clientdata, TIFFFileName(img->tif),
    458. 

  458. 		"No \"put\" routine setupl; probably can not handle image format");
    459. 

  459. 		return (0);
    460. 

  460.     }
    461. 

  461.     return (*img->get)(img, raster, w, h);
    462. 

  462. }
    463. 

  463. 

  464. /*
    465. 

  465.  * Read the specified image into an ABGR-format rastertaking in account
    466. 

  466.  * specified orientation.
    467. 

  467.  */
    468. 

  468. int
    469. 

  469. TIFFReadRGBAImageOriented(TIFF* tif,
    470. 

  470. 			  uint32 rwidth, uint32 rheight, uint32* raster,
    471. 

  471. 			  int orientation, int stop)
    472. 

  472. {
    473. 

  473.     char emsg[1024] = "";
    474. 

  474.     TIFFRGBAImage img;
    475. 

  475.     int ok;
    476. 

  476. 

  477. 	if (TIFFRGBAImageOK(tif, emsg) && TIFFRGBAImageBegin(&img, tif, stop, emsg)) {
    478. 

  478. 		img.req_orientation = orientation;
    479. 

  479. 		/* XXX verify rwidth and rheight against width and height */
    480. 

  480. 		ok = TIFFRGBAImageGet(&img, raster+(rheight-img.height)*rwidth,
    481. 

  481. 			rwidth, img.height);
    482. 

  482. 		TIFFRGBAImageEnd(&img);
    483. 

  483. 	} else {
    484. 

  484. 		TIFFErrorExt(tif->tif_clientdata, TIFFFileName(tif), "%s", emsg);
    485. 

  485. 		ok = 0;
    486. 

  486.     }
    487. 

  487.     return (ok);
    488. 

  488. }
    489. 

  489. 

  490. /*
    491. 

  491.  * Read the specified image into an ABGR-format raster. Use bottom left
    492. 

  492.  * origin for raster by default.
    493. 

  493.  */
    494. 

  494. int
    495. 

  495. TIFFReadRGBAImage(TIFF* tif,
    496. 

  496. 		  uint32 rwidth, uint32 rheight, uint32* raster, int stop)
    497. 

  497. {
    498. 

  498. 	return TIFFReadRGBAImageOriented(tif, rwidth, rheight, raster,
    499. 

  499. 					 ORIENTATION_BOTLEFT, stop);
    500. 

  500. }
    501. 

  501. 

  502. static int 
    503. 

  503. setorientation(TIFFRGBAImage* img)
    504. 

  504. {
    505. 

  505. 	switch (img->orientation) {
    506. 

  506. 		case ORIENTATION_TOPLEFT:
    507. 

  507. 		case ORIENTATION_LEFTTOP:
    508. 

  508. 			if (img->req_orientation == ORIENTATION_TOPRIGHT ||
    509. 

  509. 			    img->req_orientation == ORIENTATION_RIGHTTOP)
    510. 

  510. 				return FLIP_HORIZONTALLY;
    511. 

  511. 			else if (img->req_orientation == ORIENTATION_BOTRIGHT ||
    512. 

  512. 			    img->req_orientation == ORIENTATION_RIGHTBOT)
    513. 

  513. 				return FLIP_HORIZONTALLY | FLIP_VERTICALLY;
    514. 

  514. 			else if (img->req_orientation == ORIENTATION_BOTLEFT ||
    515. 

  515. 			    img->req_orientation == ORIENTATION_LEFTBOT)
    516. 

  516. 				return FLIP_VERTICALLY;
    517. 

  517. 			else
    518. 

  518. 				return 0;
    519. 

  519. 		case ORIENTATION_TOPRIGHT:
    520. 

  520. 		case ORIENTATION_RIGHTTOP:
    521. 

  521. 			if (img->req_orientation == ORIENTATION_TOPLEFT ||
    522. 

  522. 			    img->req_orientation == ORIENTATION_LEFTTOP)
    523. 

  523. 				return FLIP_HORIZONTALLY;
    524. 

  524. 			else if (img->req_orientation == ORIENTATION_BOTRIGHT ||
    525. 

  525. 			    img->req_orientation == ORIENTATION_RIGHTBOT)
    526. 

  526. 				return FLIP_VERTICALLY;
    527. 

  527. 			else if (img->req_orientation == ORIENTATION_BOTLEFT ||
    528. 

  528. 			    img->req_orientation == ORIENTATION_LEFTBOT)
    529. 

  529. 				return FLIP_HORIZONTALLY | FLIP_VERTICALLY;
    530. 

  530. 			else
    531. 

  531. 				return 0;
    532. 

  532. 		case ORIENTATION_BOTRIGHT:
    533. 

  533. 		case ORIENTATION_RIGHTBOT:
    534. 

  534. 			if (img->req_orientation == ORIENTATION_TOPLEFT ||
    535. 

  535. 			    img->req_orientation == ORIENTATION_LEFTTOP)
    536. 

  536. 				return FLIP_HORIZONTALLY | FLIP_VERTICALLY;
    537. 

  537. 			else if (img->req_orientation == ORIENTATION_TOPRIGHT ||
    538. 

  538. 			    img->req_orientation == ORIENTATION_RIGHTTOP)
    539. 

  539. 				return FLIP_VERTICALLY;
    540. 

  540. 			else if (img->req_orientation == ORIENTATION_BOTLEFT ||
    541. 

  541. 			    img->req_orientation == ORIENTATION_LEFTBOT)
    542. 

  542. 				return FLIP_HORIZONTALLY;
    543. 

  543. 			else
    544. 

  544. 				return 0;
    545. 

  545. 		case ORIENTATION_BOTLEFT:
    546. 

  546. 		case ORIENTATION_LEFTBOT:
    547. 

  547. 			if (img->req_orientation == ORIENTATION_TOPLEFT ||
    548. 

  548. 			    img->req_orientation == ORIENTATION_LEFTTOP)
    549. 

  549. 				return FLIP_VERTICALLY;
    550. 

  550. 			else if (img->req_orientation == ORIENTATION_TOPRIGHT ||
    551. 

  551. 			    img->req_orientation == ORIENTATION_RIGHTTOP)
    552. 

  552. 				return FLIP_HORIZONTALLY | FLIP_VERTICALLY;
    553. 

  553. 			else if (img->req_orientation == ORIENTATION_BOTRIGHT ||
    554. 

  554. 			    img->req_orientation == ORIENTATION_RIGHTBOT)
    555. 

  555. 				return FLIP_HORIZONTALLY;
    556. 

  556. 			else
    557. 

  557. 				return 0;
    558. 

  558. 		default:	/* NOTREACHED */
    559. 

  559. 			return 0;
    560. 

  560. 	}
    561. 

  561. }
    562. 

  562. 

  563. /*
    564. 

  564.  * Get an tile-organized image that has
    565. 

  565.  *	PlanarConfiguration contiguous if SamplesPerPixel > 1
    566. 

  566.  * or
    567. 

  567.  *	SamplesPerPixel == 1
    568. 

  568.  */	
    569. 

  569. static int
    570. 

  570. gtTileContig(TIFFRGBAImage* img, uint32* raster, uint32 w, uint32 h)
    571. 

  571. {
    572. 

  572.     TIFF* tif = img->tif;
    573. 

  573.     tileContigRoutine put = img->put.contig;
    574. 

  574.     uint32 col, row, y, rowstoread;
    575. 

  575.     uint32 pos;
    576. 

  576.     uint32 tw, th;
    577. 

  577.     unsigned char* buf;
    578. 

  578.     int32 fromskew, toskew;
    579. 

  579.     uint32 nrow;
    580. 

  580.     int ret = 1, flip;
    581. 

  581. 

  582.     buf = (unsigned char*) _TIFFmalloc(TIFFTileSize(tif));
    583. 

  583.     if (buf == 0) {
    584. 

  584. 		TIFFErrorExt(tif->tif_clientdata, TIFFFileName(tif), "No space for tile buffer");
    585. 

  585. 		return (0);
    586. 

  586.     }
    587. 

  587.     _TIFFmemset(buf, 0, TIFFTileSize(tif));
    588. 

  588.     TIFFGetField(tif, TIFFTAG_TILEWIDTH, &tw);
    589. 

  589.     TIFFGetField(tif, TIFFTAG_TILELENGTH, &th);
    590. 

  590. 

  591.     flip = setorientation(img);
    592. 

  592.     if (flip & FLIP_VERTICALLY) {
    593. 

  593. 	    y = h - 1;
    594. 

  594. 	    toskew = -(int32)(tw + w);
    595. 

  595.     }
    596. 

  596.     else {
    597. 

  597. 	    y = 0;
    598. 

  598. 	    toskew = -(int32)(tw - w);
    599. 

  599.     }
    600. 

  600.      
    601. 

  601.     for (row = 0; row < h; row += nrow)
    602. 

  602.     {
    603. 

  603.         rowstoread = th - (row + img->row_offset) % th;
    604. 

  604.     	nrow = (row + rowstoread > h ? h - row : rowstoread);
    605. 

  605. 	for (col = 0; col < w; col += tw) 
    606. 

  606.         {
    607. 

  607.             if (TIFFReadTile(tif, buf, col+img->col_offset,
    608. 

  608.                              row+img->row_offset, 0, 0) < 0 && img->stoponerr)
    609. 

  609.             {
    610. 

  610.                 ret = 0;
    611. 

  611.                 break;
    612. 

  612.             }
    613. 

  613. 	    
    614. 

  614.             pos = ((row+img->row_offset) % th) * TIFFTileRowSize(tif);
    615. 

  615. 

  616.     	    if (col + tw > w) 
    617. 

  617.             {
    618. 

  618.                 /*
    619. 

  619.                  * Tile is clipped horizontally.  Calculate
    620. 

  620.                  * visible portion and skewing factors.
    621. 

  621.                  */
    622. 

  622.                 uint32 npix = w - col;
    623. 

  623.                 fromskew = tw - npix;
    624. 

  624.                 (*put)(img, raster+y*w+col, col, y,
    625. 

  625.                        npix, nrow, fromskew, toskew + fromskew, buf + pos);
    626. 

  626.             }
    627. 

  627.             else 
    628. 

  628.             {
    629. 

  629.                 (*put)(img, raster+y*w+col, col, y, tw, nrow, 0, toskew, buf + pos);
    630. 

  630.             }
    631. 

  631.         }
    632. 

  632. 

  633.         y += (flip & FLIP_VERTICALLY ? -(int32) nrow : (int32) nrow);
    634. 

  634.     }
    635. 

  635.     _TIFFfree(buf);
    636. 

  636. 

  637.     if (flip & FLIP_HORIZONTALLY) {
    638. 

  638. 	    uint32 line;
    639. 

  639. 

  640. 	    for (line = 0; line < h; line++) {
    641. 

  641. 		    uint32 *left = raster + (line * w);
    642. 

  642. 		    uint32 *right = left + w - 1;
    643. 

  643. 		    
    644. 

  644. 		    while ( left < right ) {
    645. 

  645. 			    uint32 temp = *left;
    646. 

  646. 			    *left = *right;
    647. 

  647. 			    *right = temp;
    648. 

  648. 			    left++, right--;
    649. 

  649. 		    }
    650. 

  650. 	    }
    651. 

  651.     }
    652. 

  652. 

  653.     return (ret);
    654. 

  654. }
    655. 

  655. 

  656. /*
    657. 

  657.  * Get an tile-organized image that has
    658. 

  658.  *	 SamplesPerPixel > 1
    659. 

  659.  *	 PlanarConfiguration separated
    660. 

  660.  * We assume that all such images are RGB.
    661. 

  661.  */	
    662. 

  662. static int
    663. 

  663. gtTileSeparate(TIFFRGBAImage* img, uint32* raster, uint32 w, uint32 h)
    664. 

  664. {
    665. 

  665. 	TIFF* tif = img->tif;
    666. 

  666. 	tileSeparateRoutine put = img->put.separate;
    667. 

  667. 	uint32 col, row, y, rowstoread;
    668. 

  668. 	uint32 pos;
    669. 

  669. 	uint32 tw, th;
    670. 

  670. 	unsigned char* buf;
    671. 

  671. 	unsigned char* p0;
    672. 

  672. 	unsigned char* p1;
    673. 

  673. 	unsigned char* p2;
    674. 

  674. 	unsigned char* pa;
    675. 

  675. 	tsize_t tilesize;
    676. 

  676. 	int32 fromskew, toskew;
    677. 

  677. 	int alpha = img->alpha;
    678. 

  678. 	uint32 nrow;
    679. 

  679. 	int ret = 1, flip;
    680. 

  680. 

  681. 	tilesize = TIFFTileSize(tif);
    682. 

  682. 	buf = (unsigned char*) _TIFFmalloc((alpha?4:3)*tilesize);
    683. 

  683. 	if (buf == 0) {
    684. 

  684. 		TIFFErrorExt(tif->tif_clientdata, TIFFFileName(tif), "No space for tile buffer");
    685. 

  685. 		return (0);
    686. 

  686. 	}
    687. 

  687. 	_TIFFmemset(buf, 0, (alpha?4:3)*tilesize);
    688. 

  688. 	p0 = buf;
    689. 

  689. 	p1 = p0 + tilesize;
    690. 

  690. 	p2 = p1 + tilesize;
    691. 

  691. 	pa = (alpha?(p2+tilesize):NULL);
    692. 

  692. 	TIFFGetField(tif, TIFFTAG_TILEWIDTH, &tw);
    693. 

  693. 	TIFFGetField(tif, TIFFTAG_TILELENGTH, &th);
    694. 

  694. 

  695. 	flip = setorientation(img);
    696. 

  696. 	if (flip & FLIP_VERTICALLY) {
    697. 

  697. 		y = h - 1;
    698. 

  698. 		toskew = -(int32)(tw + w);
    699. 

  699. 	}
    700. 

  700. 	else {
    701. 

  701. 		y = 0;
    702. 

  702. 		toskew = -(int32)(tw - w);
    703. 

  703. 	}
    704. 

  704. 

  705. 	for (row = 0; row < h; row += nrow)
    706. 

  706. 	{
    707. 

  707. 		rowstoread = th - (row + img->row_offset) % th;
    708. 

  708. 		nrow = (row + rowstoread > h ? h - row : rowstoread);
    709. 

  709. 		for (col = 0; col < w; col += tw)
    710. 

  710. 		{
    711. 

  711. 			if (TIFFReadTile(tif, p0, col+img->col_offset,
    712. 

  712. 			    row+img->row_offset,0,0) < 0 && img->stoponerr)
    713. 

  713. 			{
    714. 

  714. 				ret = 0;
    715. 

  715. 				break;
    716. 

  716. 			}
    717. 

  717. 			if (TIFFReadTile(tif, p1, col+img->col_offset,
    718. 

  718. 			    row+img->row_offset,0,1) < 0 && img->stoponerr)
    719. 

  719. 			{
    720. 

  720. 				ret = 0;
    721. 

  721. 				break;
    722. 

  722. 			}
    723. 

  723. 			if (TIFFReadTile(tif, p2, col+img->col_offset,
    724. 

  724. 			    row+img->row_offset,0,2) < 0 && img->stoponerr)
    725. 

  725. 			{
    726. 

  726. 				ret = 0;
    727. 

  727. 				break;
    728. 

  728. 			}
    729. 

  729. 			if (alpha)
    730. 

  730. 			{
    731. 

  731. 				if (TIFFReadTile(tif,pa,col+img->col_offset,
    732. 

  732. 				    row+img->row_offset,0,3) < 0 && img->stoponerr)
    733. 

  733. 				{
    734. 

  734. 					ret = 0;
    735. 

  735. 					break;
    736. 

  736. 				}
    737. 

  737. 			}
    738. 

  738. 

  739. 			pos = ((row+img->row_offset) % th) * TIFFTileRowSize(tif);
    740. 

  740. 

  741. 			if (col + tw > w)
    742. 

  742. 			{
    743. 

  743. 				/*
    744. 

  744. 				 * Tile is clipped horizontally.  Calculate
    745. 

  745. 				 * visible portion and skewing factors.
    746. 

  746. 				 */
    747. 

  747. 				uint32 npix = w - col;
    748. 

  748. 				fromskew = tw - npix;
    749. 

  749. 				(*put)(img, raster+y*w+col, col, y,
    750. 

  750. 				    npix, nrow, fromskew, toskew + fromskew,
    751. 

  751. 				    p0 + pos, p1 + pos, p2 + pos, (alpha?(pa+pos):NULL));
    752. 

  752. 			} else {
    753. 

  753. 				(*put)(img, raster+y*w+col, col, y,
    754. 

  754. 				    tw, nrow, 0, toskew, p0 + pos, p1 + pos, p2 + pos, (alpha?(pa+pos):NULL));
    755. 

  755. 			}
    756. 

  756. 		}
    757. 

  757. 

  758. 		y += (flip & FLIP_VERTICALLY ?-(int32) nrow : (int32) nrow);
    759. 

  759. 	}
    760. 

  760. 

  761. 	if (flip & FLIP_HORIZONTALLY) {
    762. 

  762. 		uint32 line;
    763. 

  763. 

  764. 		for (line = 0; line < h; line++) {
    765. 

  765. 			uint32 *left = raster + (line * w);
    766. 

  766. 			uint32 *right = left + w - 1;
    767. 

  767. 

  768. 			while ( left < right ) {
    769. 

  769. 				uint32 temp = *left;
    770. 

  770. 				*left = *right;
    771. 

  771. 				*right = temp;
    772. 

  772. 				left++, right--;
    773. 

  773. 			}
    774. 

  774. 		}
    775. 

  775. 	}
    776. 

  776. 

  777. 	_TIFFfree(buf);
    778. 

  778. 	return (ret);
    779. 

  779. }
    780. 

  780. 

  781. /*
    782. 

  782.  * Get a strip-organized image that has
    783. 

  783.  *	PlanarConfiguration contiguous if SamplesPerPixel > 1
    784. 

  784.  * or
    785. 

  785.  *	SamplesPerPixel == 1
    786. 

  786.  */	
    787. 

  787. static int
    788. 

  788. gtStripContig(TIFFRGBAImage* img, uint32* raster, uint32 w, uint32 h)
    789. 

  789. {
    790. 

  790. 	TIFF* tif = img->tif;
    791. 

  791. 	tileContigRoutine put = img->put.contig;
    792. 

  792. 	uint32 row, y, nrow, nrowsub, rowstoread;
    793. 

  793. 	uint32 pos;
    794. 

  794. 	unsigned char* buf;
    795. 

  795. 	uint32 rowsperstrip;
    796. 

  796. 	uint16 subsamplinghor,subsamplingver;
    797. 

  797. 	uint32 imagewidth = img->width;
    798. 

  798. 	tsize_t scanline;
    799. 

  799. 	int32 fromskew, toskew;
    800. 

  800. 	int ret = 1, flip;
    801. 

  801. 

  802. 	buf = (unsigned char*) _TIFFmalloc(TIFFStripSize(tif));
    803. 

  803. 	if (buf == 0) {
    804. 

  804. 		TIFFErrorExt(tif->tif_clientdata, TIFFFileName(tif), "No space for strip buffer");
    805. 

  805. 		return (0);
    806. 

  806. 	}
    807. 

  807. 	_TIFFmemset(buf, 0, TIFFStripSize(tif));
    808. 

  808. 

  809. 	flip = setorientation(img);
    810. 

  810. 	if (flip & FLIP_VERTICALLY) {
    811. 

  811. 		y = h - 1;
    812. 

  812. 		toskew = -(int32)(w + w);
    813. 

  813. 	} else {
    814. 

  814. 		y = 0;
    815. 

  815. 		toskew = -(int32)(w - w);
    816. 

  816. 	}
    817. 

  817. 

  818. 	TIFFGetFieldDefaulted(tif, TIFFTAG_ROWSPERSTRIP, &rowsperstrip);
    819. 

  819. 	TIFFGetFieldDefaulted(tif, TIFFTAG_YCBCRSUBSAMPLING, &subsamplinghor, &subsamplingver);
    820. 

  820. 	scanline = TIFFNewScanlineSize(tif);
    821. 

  821. 	fromskew = (w < imagewidth ? imagewidth - w : 0);
    822. 

  822. 	for (row = 0; row < h; row += nrow)
    823. 

  823. 	{
    824. 

  824. 		rowstoread = rowsperstrip - (row + img->row_offset) % rowsperstrip;
    825. 

  825. 		nrow = (row + rowstoread > h ? h - row : rowstoread);
    826. 

  826. 		nrowsub = nrow;
    827. 

  827. 		if ((nrowsub%subsamplingver)!=0)
    828. 

  828. 			nrowsub+=subsamplingver-nrowsub%subsamplingver;
    829. 

  829. 		if (TIFFReadEncodedStrip(tif,
    830. 

  830. 		    TIFFComputeStrip(tif,row+img->row_offset, 0),
    831. 

  831. 		    buf,
    832. 

  832. 		    ((row + img->row_offset)%rowsperstrip + nrowsub) * scanline) < 0
    833. 

  833. 		    && img->stoponerr)
    834. 

  834. 		{
    835. 

  835. 			ret = 0;
    836. 

  836. 			break;
    837. 

  837. 		}
    838. 

  838. 

  839. 		pos = ((row + img->row_offset) % rowsperstrip) * scanline;
    840. 

  840. 		(*put)(img, raster+y*w, 0, y, w, nrow, fromskew, toskew, buf + pos);
    841. 

  841. 		y += (flip & FLIP_VERTICALLY ? -(int32) nrow : (int32) nrow);
    842. 

  842. 	}
    843. 

  843. 

  844. 	if (flip & FLIP_HORIZONTALLY) {
    845. 

  845. 		uint32 line;
    846. 

  846. 

  847. 		for (line = 0; line < h; line++) {
    848. 

  848. 			uint32 *left = raster + (line * w);
    849. 

  849. 			uint32 *right = left + w - 1;
    850. 

  850. 

  851. 			while ( left < right ) {
    852. 

  852. 				uint32 temp = *left;
    853. 

  853. 				*left = *right;
    854. 

  854. 				*right = temp;
    855. 

  855. 				left++, right--;
    856. 

  856. 			}
    857. 

  857. 		}
    858. 

  858. 	}
    859. 

  859. 

  860. 	_TIFFfree(buf);
    861. 

  861. 	return (ret);
    862. 

  862. }
    863. 

  863. 

  864. /*
    865. 

  865.  * Get a strip-organized image with
    866. 

  866.  *	 SamplesPerPixel > 1
    867. 

  867.  *	 PlanarConfiguration separated
    868. 

  868.  * We assume that all such images are RGB.
    869. 

  869.  */
    870. 

  870. static int
    871. 

  871. gtStripSeparate(TIFFRGBAImage* img, uint32* raster, uint32 w, uint32 h)
    872. 

  872. {
    873. 

  873. 	TIFF* tif = img->tif;
    874. 

  874. 	tileSeparateRoutine put = img->put.separate;
    875. 

  875. 	unsigned char *buf;
    876. 

  876. 	unsigned char *p0, *p1, *p2, *pa;
    877. 

  877. 	uint32 row, y, nrow, rowstoread;
    878. 

  878. 	uint32 pos;
    879. 

  879. 	tsize_t scanline;
    880. 

  880. 	uint32 rowsperstrip, offset_row;
    881. 

  881. 	uint32 imagewidth = img->width;
    882. 

  882. 	tsize_t stripsize;
    883. 

  883. 	int32 fromskew, toskew;
    884. 

  884. 	int alpha = img->alpha;
    885. 

  885. 	int ret = 1, flip;
    886. 

  886. 

  887. 	stripsize = TIFFStripSize(tif);
    888. 

  888. 	p0 = buf = (unsigned char *)_TIFFmalloc((alpha?4:3)*stripsize);
    889. 

  889. 	if (buf == 0) {
    890. 

  890. 		TIFFErrorExt(tif->tif_clientdata, TIFFFileName(tif), "No space for tile buffer");
    891. 

  891. 		return (0);
    892. 

  892. 	}
    893. 

  893. 	_TIFFmemset(buf, 0, (alpha?4:3)*stripsize);
    894. 

  894. 	p1 = p0 + stripsize;
    895. 

  895. 	p2 = p1 + stripsize;
    896. 

  896. 	pa = (alpha?(p2+stripsize):NULL);
    897. 

  897. 

  898. 	flip = setorientation(img);
    899. 

  899. 	if (flip & FLIP_VERTICALLY) {
    900. 

  900. 		y = h - 1;
    901. 

  901. 		toskew = -(int32)(w + w);
    902. 

  902. 	}
    903. 

  903. 	else {
    904. 

  904. 		y = 0;
    905. 

  905. 		toskew = -(int32)(w - w);
    906. 

  906. 	}
    907. 

  907. 

  908. 	TIFFGetFieldDefaulted(tif, TIFFTAG_ROWSPERSTRIP, &rowsperstrip);
    909. 

  909. 	scanline = TIFFScanlineSize(tif);
    910. 

  910. 	fromskew = (w < imagewidth ? imagewidth - w : 0);
    911. 

  911. 	for (row = 0; row < h; row += nrow)
    912. 

  912. 	{
    913. 

  913. 		rowstoread = rowsperstrip - (row + img->row_offset) % rowsperstrip;
    914. 

  914. 		nrow = (row + rowstoread > h ? h - row : rowstoread);
    915. 

  915. 		offset_row = row + img->row_offset;
    916. 

  916. 		if (TIFFReadEncodedStrip(tif, TIFFComputeStrip(tif, offset_row, 0),
    917. 

  917. 		    p0, ((row + img->row_offset)%rowsperstrip + nrow) * scanline) < 0
    918. 

  918. 		    && img->stoponerr)
    919. 

  919. 		{
    920. 

  920. 			ret = 0;
    921. 

  921. 			break;
    922. 

  922. 		}
    923. 

  923. 		if (TIFFReadEncodedStrip(tif, TIFFComputeStrip(tif, offset_row, 1),
    924. 

  924. 		    p1, ((row + img->row_offset)%rowsperstrip + nrow) * scanline) < 0
    925. 

  925. 		    && img->stoponerr)
    926. 

  926. 		{
    927. 

  927. 			ret = 0;
    928. 

  928. 			break;
    929. 

  929. 		}
    930. 

  930. 		if (TIFFReadEncodedStrip(tif, TIFFComputeStrip(tif, offset_row, 2),
    931. 

  931. 		    p2, ((row + img->row_offset)%rowsperstrip + nrow) * scanline) < 0
    932. 

  932. 		    && img->stoponerr)
    933. 

  933. 		{
    934. 

  934. 			ret = 0;
    935. 

  935. 			break;
    936. 

  936. 		}
    937. 

  937. 		if (alpha)
    938. 

  938. 		{
    939. 

  939. 			if (TIFFReadEncodedStrip(tif, TIFFComputeStrip(tif, offset_row, 3),
    940. 

  940. 			    pa, ((row + img->row_offset)%rowsperstrip + nrow) * scanline) < 0
    941. 

  941. 			    && img->stoponerr)
    942. 

  942. 			{
    943. 

  943. 				ret = 0;
    944. 

  944. 				break;
    945. 

  945. 			}
    946. 

  946. 		}
    947. 

  947. 

  948. 		pos = ((row + img->row_offset) % rowsperstrip) * scanline;
    949. 

  949. 		(*put)(img, raster+y*w, 0, y, w, nrow, fromskew, toskew, p0 + pos, p1 + pos,
    950. 

  950. 		    p2 + pos, (alpha?(pa+pos):NULL));
    951. 

  951. 		y += (flip & FLIP_VERTICALLY ? -(int32) nrow : (int32) nrow);
    952. 

  952. 	}
    953. 

  953. 

  954. 	if (flip & FLIP_HORIZONTALLY) {
    955. 

  955. 		uint32 line;
    956. 

  956. 

  957. 		for (line = 0; line < h; line++) {
    958. 

  958. 			uint32 *left = raster + (line * w);
    959. 

  959. 			uint32 *right = left + w - 1;
    960. 

  960. 

  961. 			while ( left < right ) {
    962. 

  962. 				uint32 temp = *left;
    963. 

  963. 				*left = *right;
    964. 

  964. 				*right = temp;
    965. 

  965. 				left++, right--;
    966. 

  966. 			}
    967. 

  967. 		}
    968. 

  968. 	}
    969. 

  969. 

  970. 	_TIFFfree(buf);
    971. 

  971. 	return (ret);
    972. 

  972. }
    973. 

  973. 

  974. /*
    975. 

  975.  * The following routines move decoded data returned
    976. 

  976.  * from the TIFF library into rasters filled with packed
    977. 

  977.  * ABGR pixels (i.e. suitable for passing to lrecwrite.)
    978. 

  978.  *
    979. 

  979.  * The routines have been created according to the most
    980. 

  980.  * important cases and optimized.  PickContigCase and
    981. 

  981.  * PickSeparateCase analyze the parameters and select
    982. 

  982.  * the appropriate "get" and "put" routine to use.
    983. 

  983.  */
    984. 

  984. #define	REPEAT8(op)	REPEAT4(op); REPEAT4(op)
    985. 

  985. #define	REPEAT4(op)	REPEAT2(op); REPEAT2(op)
    986. 

  986. #define	REPEAT2(op)	op; op
    987. 

  987. #define	CASE8(x,op)			\
    988. 

  988.     switch (x) {			\
    989. 

  989.     case 7: op; case 6: op; case 5: op;	\
    990. 

  990.     case 4: op; case 3: op; case 2: op;	\
    991. 

  991.     case 1: op;				\
    992. 

  992.     }
    993. 

  993. #define	CASE4(x,op)	switch (x) { case 3: op; case 2: op; case 1: op; }
    994. 

  994. #define	NOP
    995. 

  995. 

  996. #define	UNROLL8(w, op1, op2) {		\
    997. 

  997.     uint32 _x;				\
    998. 

  998.     for (_x = w; _x >= 8; _x -= 8) {	\
    999. 

  999. 	op1;				\
    1000. 

  1000. 	REPEAT8(op2);			\
    1001. 

  1001.     }					\
    1002. 

  1002.     if (_x > 0) {			\
    1003. 

  1003. 	op1;				\
    1004. 

  1004. 	CASE8(_x,op2);			\
    1005. 

  1005.     }					\
    1006. 

  1006. }
    1007. 

  1007. #define	UNROLL4(w, op1, op2) {		\
    1008. 

  1008.     uint32 _x;				\
    1009. 

  1009.     for (_x = w; _x >= 4; _x -= 4) {	\
    1010. 

  1010. 	op1;				\
    1011. 

  1011. 	REPEAT4(op2);			\
    1012. 

  1012.     }					\
    1013. 

  1013.     if (_x > 0) {			\
    1014. 

  1014. 	op1;				\
    1015. 

  1015. 	CASE4(_x,op2);			\
    1016. 

  1016.     }					\
    1017. 

  1017. }
    1018. 

  1018. #define	UNROLL2(w, op1, op2) {		\
    1019. 

  1019.     uint32 _x;				\
    1020. 

  1020.     for (_x = w; _x >= 2; _x -= 2) {	\
    1021. 

  1021. 	op1;				\
    1022. 

  1022. 	REPEAT2(op2);			\
    1023. 

  1023.     }					\
    1024. 

  1024.     if (_x) {				\
    1025. 

  1025. 	op1;				\
    1026. 

  1026. 	op2;				\
    1027. 

  1027.     }					\
    1028. 

  1028. }
    1029. 

  1029.     
    1030. 

  1030. #define	SKEW(r,g,b,skew)	{ r += skew; g += skew; b += skew; }
    1031. 

  1031. #define	SKEW4(r,g,b,a,skew)	{ r += skew; g += skew; b += skew; a+= skew; }
    1032. 

  1032. 

  1033. #define A1 (((uint32)0xffL)<<24)
    1034. 

  1034. #define	PACK(r,g,b)	\
    1035. 

  1035. 	((uint32)(r)|((uint32)(g)<<8)|((uint32)(b)<<16)|A1)
    1036. 

  1036. #define	PACK4(r,g,b,a)	\
    1037. 

  1037. 	((uint32)(r)|((uint32)(g)<<8)|((uint32)(b)<<16)|((uint32)(a)<<24))
    1038. 

  1038. #define W2B(v) (((v)>>8)&0xff)
    1039. 

  1039. #define	PACKW(r,g,b)	\
    1040. 

  1040. 	((uint32)W2B(r)|((uint32)W2B(g)<<8)|((uint32)W2B(b)<<16)|A1)
    1041. 

  1041. #define	PACKW4(r,g,b,a)	\
    1042. 

  1042. 	((uint32)W2B(r)|((uint32)W2B(g)<<8)|((uint32)W2B(b)<<16)|((uint32)W2B(a)<<24))
    1043. 

  1043. 

  1044. #define	DECLAREContigPutFunc(name) \
    1045. 

  1045. static void name(\
    1046. 

  1046.     TIFFRGBAImage* img, \
    1047. 

  1047.     uint32* cp, \
    1048. 

  1048.     uint32 x, uint32 y, \
    1049. 

  1049.     uint32 w, uint32 h, \
    1050. 

  1050.     int32 fromskew, int32 toskew, \
    1051. 

  1051.     unsigned char* pp \
    1052. 

  1052. )
    1053. 

  1053. 

  1054. /*
    1055. 

  1055.  * 8-bit palette => colormap/RGB
    1056. 

  1056.  */
    1057. 

  1057. DECLAREContigPutFunc(put8bitcmaptile)
    1058. 

  1058. {
    1059. 

  1059.     uint32** PALmap = img->PALmap;
    1060. 

  1060.     int samplesperpixel = img->samplesperpixel;
    1061. 

  1061. 

  1062.     (void) y;
    1063. 

  1063.     while (h-- > 0) {
    1064. 

  1064. 	for (x = w; x-- > 0;)
    1065. 

  1065.         {
    1066. 

  1066. 	    *cp++ = PALmap[*pp][0];
    1067. 

  1067.             pp += samplesperpixel;
    1068. 

  1068.         }
    1069. 

  1069. 	cp += toskew;
    1070. 

  1070. 	pp += fromskew;
    1071. 

  1071.     }
    1072. 

  1072. }
    1073. 

  1073. 

  1074. /*
    1075. 

  1075.  * 4-bit palette => colormap/RGB
    1076. 

  1076.  */
    1077. 

  1077. DECLAREContigPutFunc(put4bitcmaptile)
    1078. 

  1078. {
    1079. 

  1079.     uint32** PALmap = img->PALmap;
    1080. 

  1080. 

  1081.     (void) x; (void) y;
    1082. 

  1082.     fromskew /= 2;
    1083. 

  1083.     while (h-- > 0) {
    1084. 

  1084. 	uint32* bw;
    1085. 

  1085. 	UNROLL2(w, bw = PALmap[*pp++], *cp++ = *bw++);
    1086. 

  1086. 	cp += toskew;
    1087. 

  1087. 	pp += fromskew;
    1088. 

  1088.     }
    1089. 

  1089. }
    1090. 

  1090. 

  1091. /*
    1092. 

  1092.  * 2-bit palette => colormap/RGB
    1093. 

  1093.  */
    1094. 

  1094. DECLAREContigPutFunc(put2bitcmaptile)
    1095. 

  1095. {
    1096. 

  1096.     uint32** PALmap = img->PALmap;
    1097. 

  1097. 

  1098.     (void) x; (void) y;
    1099. 

  1099.     fromskew /= 4;
    1100. 

  1100.     while (h-- > 0) {
    1101. 

  1101. 	uint32* bw;
    1102. 

  1102. 	UNROLL4(w, bw = PALmap[*pp++], *cp++ = *bw++);
    1103. 

  1103. 	cp += toskew;
    1104. 

  1104. 	pp += fromskew;
    1105. 

  1105.     }
    1106. 

  1106. }
    1107. 

  1107. 

  1108. /*
    1109. 

  1109.  * 1-bit palette => colormap/RGB
    1110. 

  1110.  */
    1111. 

  1111. DECLAREContigPutFunc(put1bitcmaptile)
    1112. 

  1112. {
    1113. 

  1113.     uint32** PALmap = img->PALmap;
    1114. 

  1114. 

  1115.     (void) x; (void) y;
    1116. 

  1116.     fromskew /= 8;
    1117. 

  1117.     while (h-- > 0) {
    1118. 

  1118. 	uint32* bw;
    1119. 

  1119. 	UNROLL8(w, bw = PALmap[*pp++], *cp++ = *bw++);
    1120. 

  1120. 	cp += toskew;
    1121. 

  1121. 	pp += fromskew;
    1122. 

  1122.     }
    1123. 

  1123. }
    1124. 

  1124. 

  1125. /*
    1126. 

  1126.  * 8-bit greyscale => colormap/RGB
    1127. 

  1127.  */
    1128. 

  1128. DECLAREContigPutFunc(putgreytile)
    1129. 

  1129. {
    1130. 

  1130.     int samplesperpixel = img->samplesperpixel;
    1131. 

  1131.     uint32** BWmap = img->BWmap;
    1132. 

  1132. 

  1133.     (void) y;
    1134. 

  1134.     while (h-- > 0) {
    1135. 

  1135. 	for (x = w; x-- > 0;)
    1136. 

  1136.         {
    1137. 

  1137. 	    *cp++ = BWmap[*pp][0];
    1138. 

  1138.             pp += samplesperpixel;
    1139. 

  1139.         }
    1140. 

  1140. 	cp += toskew;
    1141. 

  1141. 	pp += fromskew;
    1142. 

  1142.     }
    1143. 

  1143. }
    1144. 

  1144. 

  1145. /*
    1146. 

  1146.  * 16-bit greyscale => colormap/RGB
    1147. 

  1147.  */
    1148. 

  1148. DECLAREContigPutFunc(put16bitbwtile)
    1149. 

  1149. {
    1150. 

  1150.     int samplesperpixel = img->samplesperpixel;
    1151. 

  1151.     uint32** BWmap = img->BWmap;
    1152. 

  1152. 

  1153.     (void) y;
    1154. 

  1154.     while (h-- > 0) {
    1155. 

  1155.         uint16 *wp = (uint16 *) pp;
    1156. 

  1156. 

  1157. 	for (x = w; x-- > 0;)
    1158. 

  1158.         {
    1159. 

  1159.             /* use high order byte of 16bit value */
    1160. 

  1160. 

  1161. 	    *cp++ = BWmap[*wp >> 8][0];
    1162. 

  1162.             pp += 2 * samplesperpixel;
    1163. 

  1163.             wp += samplesperpixel;
    1164. 

  1164.         }
    1165. 

  1165. 	cp += toskew;
    1166. 

  1166. 	pp += fromskew;
    1167. 

  1167.     }
    1168. 

  1168. }
    1169. 

  1169. 

  1170. /*
    1171. 

  1171.  * 1-bit bilevel => colormap/RGB
    1172. 

  1172.  */
    1173. 

  1173. DECLAREContigPutFunc(put1bitbwtile)
    1174. 

  1174. {
    1175. 

  1175.     uint32** BWmap = img->BWmap;
    1176. 

  1176. 

  1177.     (void) x; (void) y;
    1178. 

  1178.     fromskew /= 8;
    1179. 

  1179.     while (h-- > 0) {
    1180. 

  1180. 	uint32* bw;
    1181. 

  1181. 	UNROLL8(w, bw = BWmap[*pp++], *cp++ = *bw++);
    1182. 

  1182. 	cp += toskew;
    1183. 

  1183. 	pp += fromskew;
    1184. 

  1184.     }
    1185. 

  1185. }
    1186. 

  1186. 

  1187. /*
    1188. 

  1188.  * 2-bit greyscale => colormap/RGB
    1189. 

  1189.  */
    1190. 

  1190. DECLAREContigPutFunc(put2bitbwtile)
    1191. 

  1191. {
    1192. 

  1192.     uint32** BWmap = img->BWmap;
    1193. 

  1193. 

  1194.     (void) x; (void) y;
    1195. 

  1195.     fromskew /= 4;
    1196. 

  1196.     while (h-- > 0) {
    1197. 

  1197. 	uint32* bw;
    1198. 

  1198. 	UNROLL4(w, bw = BWmap[*pp++], *cp++ = *bw++);
    1199. 

  1199. 	cp += toskew;
    1200. 

  1200. 	pp += fromskew;
    1201. 

  1201.     }
    1202. 

  1202. }
    1203. 

  1203. 

  1204. /*
    1205. 

  1205.  * 4-bit greyscale => colormap/RGB
    1206. 

  1206.  */
    1207. 

  1207. DECLAREContigPutFunc(put4bitbwtile)
    1208. 

  1208. {
    1209. 

  1209.     uint32** BWmap = img->BWmap;
    1210. 

  1210. 

  1211.     (void) x; (void) y;
    1212. 

  1212.     fromskew /= 2;
    1213. 

  1213.     while (h-- > 0) {
    1214. 

  1214. 	uint32* bw;
    1215. 

  1215. 	UNROLL2(w, bw = BWmap[*pp++], *cp++ = *bw++);
    1216. 

  1216. 	cp += toskew;
    1217. 

  1217. 	pp += fromskew;
    1218. 

  1218.     }
    1219. 

  1219. }
    1220. 

  1220. 

  1221. /*
    1222. 

  1222.  * 8-bit packed samples, no Map => RGB
    1223. 

  1223.  */
    1224. 

  1224. DECLAREContigPutFunc(putRGBcontig8bittile)
    1225. 

  1225. {
    1226. 

  1226.     int samplesperpixel = img->samplesperpixel;
    1227. 

  1227. 

  1228.     (void) x; (void) y;
    1229. 

  1229.     fromskew *= samplesperpixel;
    1230. 

  1230.     while (h-- > 0) {
    1231. 

  1231. 	UNROLL8(w, NOP,
    1232. 

  1232. 	    *cp++ = PACK(pp[0], pp[1], pp[2]);
    1233. 

  1233. 	    pp += samplesperpixel);
    1234. 

  1234. 	cp += toskew;
    1235. 

  1235. 	pp += fromskew;
    1236. 

  1236.     }
    1237. 

  1237. }
    1238. 

  1238. 

  1239. /*
    1240. 

  1240.  * 8-bit packed samples => RGBA w/ associated alpha
    1241. 

  1241.  * (known to have Map == NULL)
    1242. 

  1242.  */
    1243. 

  1243. DECLAREContigPutFunc(putRGBAAcontig8bittile)
    1244. 

  1244. {
    1245. 

  1245.     int samplesperpixel = img->samplesperpixel;
    1246. 

  1246. 

  1247.     (void) x; (void) y;
    1248. 

  1248.     fromskew *= samplesperpixel;
    1249. 

  1249.     while (h-- > 0) {
    1250. 

  1250. 	UNROLL8(w, NOP,
    1251. 

  1251. 	    *cp++ = PACK4(pp[0], pp[1], pp[2], pp[3]);
    1252. 

  1252. 	    pp += samplesperpixel);
    1253. 

  1253. 	cp += toskew;
    1254. 

  1254. 	pp += fromskew;
    1255. 

  1255.     }
    1256. 

  1256. }
    1257. 

  1257. 

  1258. /*
    1259. 

  1259.  * 8-bit packed samples => RGBA w/ unassociated alpha
    1260. 

  1260.  * (known to have Map == NULL)
    1261. 

  1261.  */
    1262. 

  1262. DECLAREContigPutFunc(putRGBUAcontig8bittile)
    1263. 

  1263. {
    1264. 

  1264. 	int samplesperpixel = img->samplesperpixel;
    1265. 

  1265. 	(void) y;
    1266. 

  1266. 	fromskew *= samplesperpixel;
    1267. 

  1267. 	while (h-- > 0) {
    1268. 

  1268. 		uint32 r, g, b, a;
    1269. 

  1269. 		for (x = w; x-- > 0;) {
    1270. 

  1270. 			a = pp[3];
    1271. 

  1271.                         r = (a*pp[0] + 127) / 255;
    1272. 

  1272.                         g = (a*pp[1] + 127) / 255;
    1273. 

  1273.                         b = (a*pp[2] + 127) / 255;
    1274. 

  1274. 			*cp++ = PACK4(r,g,b,a);
    1275. 

  1275. 			pp += samplesperpixel;
    1276. 

  1276. 		}
    1277. 

  1277. 		cp += toskew;
    1278. 

  1278. 		pp += fromskew;
    1279. 

  1279. 	}
    1280. 

  1280. }
    1281. 

  1281. 

  1282. /*
    1283. 

  1283.  * 16-bit packed samples => RGB
    1284. 

  1284.  */
    1285. 

  1285. DECLAREContigPutFunc(putRGBcontig16bittile)
    1286. 

  1286. {
    1287. 

  1287. 	int samplesperpixel = img->samplesperpixel;
    1288. 

  1288. 	uint16 *wp = (uint16 *)pp;
    1289. 

  1289. 	(void) y;
    1290. 

  1290. 	fromskew *= samplesperpixel;
    1291. 

  1291. 	while (h-- > 0) {
    1292. 

  1292. 		for (x = w; x-- > 0;) {
    1293. 

  1293.                     *cp++ = PACKW(wp[0],wp[1],wp[2]);
    1294. 

  1294.                     wp += samplesperpixel;
    1295. 

  1295. 		}
    1296. 

  1296. 		cp += toskew;
    1297. 

  1297. 		wp += fromskew;
    1298. 

  1298. 	}
    1299. 

  1299. }
    1300. 

  1300. 

  1301. /*
    1302. 

  1302.  * 16-bit packed samples => RGBA w/ associated alpha
    1303. 

  1303.  * (known to have Map == NULL)
    1304. 

  1304.  */
    1305. 

  1305. DECLAREContigPutFunc(putRGBAAcontig16bittile)
    1306. 

  1306. {
    1307. 

  1307. 	int samplesperpixel = img->samplesperpixel;
    1308. 

  1308. 	uint16 *wp = (uint16 *)pp;
    1309. 

  1309. 	(void) y;
    1310. 

  1310. 	fromskew *= samplesperpixel;
    1311. 

  1311. 	while (h-- > 0) {
    1312. 

  1312. 		for (x = w; x-- > 0;) {
    1313. 

  1313.                     *cp++ = PACKW4(wp[0],wp[1],wp[2],wp[3]);
    1314. 

  1314.                     wp += samplesperpixel;
    1315. 

  1315. 		}
    1316. 

  1316. 		cp += toskew;
    1317. 

  1317. 		wp += fromskew;
    1318. 

  1318. 	}
    1319. 

  1319. }
    1320. 

  1320. 

  1321. /*
    1322. 

  1322.  * 16-bit packed samples => RGBA w/ unassociated alpha
    1323. 

  1323.  * (known to have Map == NULL)
    1324. 

  1324.  */
    1325. 

  1325. DECLAREContigPutFunc(putRGBUAcontig16bittile)
    1326. 

  1326. {
    1327. 

  1327. 	int samplesperpixel = img->samplesperpixel;
    1328. 

  1328. 	uint16 *wp = (uint16 *)pp;
    1329. 

  1329. 	(void) y;
    1330. 

  1330. 	fromskew *= samplesperpixel;
    1331. 

  1331. 	while (h-- > 0) {
    1332. 

  1332. 		uint32 r,g,b,a;
    1333. 

  1333. 		for (x = w; x-- > 0;) {
    1334. 

  1334.                     a = W2B(wp[3]);
    1335. 

  1335.                     r = (a*W2B(wp[0]) + 127) / 255;
    1336. 

  1336.                     g = (a*W2B(wp[1]) + 127) / 255;
    1337. 

  1337.                     b = (a*W2B(wp[2]) + 127) / 255;
    1338. 

  1338.                     *cp++ = PACK4(r,g,b,a);
    1339. 

  1339.                     wp += samplesperpixel;
    1340. 

  1340. 		}
    1341. 

  1341. 		cp += toskew;
    1342. 

  1342. 		wp += fromskew;
    1343. 

  1343. 	}
    1344. 

  1344. }
    1345. 

  1345. 

  1346. /*
    1347. 

  1347.  * 8-bit packed CMYK samples w/o Map => RGB
    1348. 

  1348.  *
    1349. 

  1349.  * NB: The conversion of CMYK->RGB is *very* crude.
    1350. 

  1350.  */
    1351. 

  1351. DECLAREContigPutFunc(putRGBcontig8bitCMYKtile)
    1352. 

  1352. {
    1353. 

  1353.     int samplesperpixel = img->samplesperpixel;
    1354. 

  1354.     uint16 r, g, b, k;
    1355. 

  1355. 

  1356.     (void) x; (void) y;
    1357. 

  1357.     fromskew *= samplesperpixel;
    1358. 

  1358.     while (h-- > 0) {
    1359. 

  1359. 	UNROLL8(w, NOP,
    1360. 

  1360. 	    k = 255 - pp[3];
    1361. 

  1361. 	    r = (k*(255-pp[0]))/255;
    1362. 

  1362. 	    g = (k*(255-pp[1]))/255;
    1363. 

  1363. 	    b = (k*(255-pp[2]))/255;
    1364. 

  1364. 	    *cp++ = PACK(r, g, b);
    1365. 

  1365. 	    pp += samplesperpixel);
    1366. 

  1366. 	cp += toskew;
    1367. 

  1367. 	pp += fromskew;
    1368. 

  1368.     }
    1369. 

  1369. }
    1370. 

  1370. 

  1371. /*
    1372. 

  1372.  * 8-bit packed CMYK samples w/Map => RGB
    1373. 

  1373.  *
    1374. 

  1374.  * NB: The conversion of CMYK->RGB is *very* crude.
    1375. 

  1375.  */
    1376. 

  1376. DECLAREContigPutFunc(putRGBcontig8bitCMYKMaptile)
    1377. 

  1377. {
    1378. 

  1378.     int samplesperpixel = img->samplesperpixel;
    1379. 

  1379.     TIFFRGBValue* Map = img->Map;
    1380. 

  1380.     uint16 r, g, b, k;
    1381. 

  1381. 

  1382.     (void) y;
    1383. 

  1383.     fromskew *= samplesperpixel;
    1384. 

  1384.     while (h-- > 0) {
    1385. 

  1385. 	for (x = w; x-- > 0;) {
    1386. 

  1386. 	    k = 255 - pp[3];
    1387. 

  1387. 	    r = (k*(255-pp[0]))/255;
    1388. 

  1388. 	    g = (k*(255-pp[1]))/255;
    1389. 

  1389. 	    b = (k*(255-pp[2]))/255;
    1390. 

  1390. 	    *cp++ = PACK(Map[r], Map[g], Map[b]);
    1391. 

  1391. 	    pp += samplesperpixel;
    1392. 

  1392. 	}
    1393. 

  1393. 	pp += fromskew;
    1394. 

  1394. 	cp += toskew;
    1395. 

  1395.     }
    1396. 

  1396. }
    1397. 

  1397. 

  1398. #define	DECLARESepPutFunc(name) \
    1399. 

  1399. static void name(\
    1400. 

  1400.     TIFFRGBAImage* img,\
    1401. 

  1401.     uint32* cp,\
    1402. 

  1402.     uint32 x, uint32 y, \
    1403. 

  1403.     uint32 w, uint32 h,\
    1404. 

  1404.     int32 fromskew, int32 toskew,\
    1405. 

  1405.     unsigned char* r, unsigned char* g, unsigned char* b, unsigned char* a\
    1406. 

  1406. )
    1407. 

  1407. 

  1408. /*
    1409. 

  1409.  * 8-bit unpacked samples => RGB
    1410. 

  1410.  */
    1411. 

  1411. DECLARESepPutFunc(putRGBseparate8bittile)
    1412. 

  1412. {
    1413. 

  1413.     (void) img; (void) x; (void) y; (void) a;
    1414. 

  1414.     while (h-- > 0) {
    1415. 

  1415. 	UNROLL8(w, NOP, *cp++ = PACK(*r++, *g++, *b++));
    1416. 

  1416. 	SKEW(r, g, b, fromskew);
    1417. 

  1417. 	cp += toskew;
    1418. 

  1418.     }
    1419. 

  1419. }
    1420. 

  1420. 

  1421. /*
    1422. 

  1422.  * 8-bit unpacked samples => RGBA w/ associated alpha
    1423. 

  1423.  */
    1424. 

  1424. DECLARESepPutFunc(putRGBAAseparate8bittile)
    1425. 

  1425. {
    1426. 

  1426. 	(void) img; (void) x; (void) y;
    1427. 

  1427. 	while (h-- > 0) {
    1428. 

  1428. 		UNROLL8(w, NOP, *cp++ = PACK4(*r++, *g++, *b++, *a++));
    1429. 

  1429. 		SKEW4(r, g, b, a, fromskew);
    1430. 

  1430. 		cp += toskew;
    1431. 

  1431. 	}
    1432. 

  1432. }
    1433. 

  1433. 

  1434. /*
    1435. 

  1435.  * 8-bit unpacked samples => RGBA w/ unassociated alpha
    1436. 

  1436.  */
    1437. 

  1437. DECLARESepPutFunc(putRGBUAseparate8bittile)
    1438. 

  1438. {
    1439. 

  1439. 	(void) img; (void) y;
    1440. 

  1440. 	while (h-- > 0) {
    1441. 

  1441. 		uint32 rv, gv, bv, av;
    1442. 

  1442. 		for (x = w; x-- > 0;) {
    1443. 

  1443. 			av = *a++;
    1444. 

  1444.                         rv = (av* *r++ + 127) / 255;
    1445. 

  1445.                         gv = (av* *g++ + 127) / 255;
    1446. 

  1446.                         bv = (av* *b++ + 127) / 255;
    1447. 

  1447. 			*cp++ = PACK4(rv,gv,bv,av);
    1448. 

  1448. 		}
    1449. 

  1449. 		SKEW4(r, g, b, a, fromskew);
    1450. 

  1450. 		cp += toskew;
    1451. 

  1451. 	}
    1452. 

  1452. }
    1453. 

  1453. 

  1454. /*
    1455. 

  1455.  * 16-bit unpacked samples => RGB
    1456. 

  1456.  */
    1457. 

  1457. DECLARESepPutFunc(putRGBseparate16bittile)
    1458. 

  1458. {
    1459. 

  1459. 	uint16 *wr = (uint16*) r;
    1460. 

  1460. 	uint16 *wg = (uint16*) g;
    1461. 

  1461. 	uint16 *wb = (uint16*) b;
    1462. 

  1462. 	(void) img; (void) y; (void) a;
    1463. 

  1463. 	while (h-- > 0) {
    1464. 

  1464. 		for (x = 0; x < w; x++)
    1465. 

  1465.                     *cp++ = PACKW(*wr++,*wg++,*wb++);
    1466. 

  1466. 		SKEW(wr, wg, wb, fromskew);
    1467. 

  1467. 		cp += toskew;
    1468. 

  1468. 	}
    1469. 

  1469. }
    1470. 

  1470. 

  1471. /*
    1472. 

  1472.  * 16-bit unpacked samples => RGBA w/ associated alpha
    1473. 

  1473.  */
    1474. 

  1474. DECLARESepPutFunc(putRGBAAseparate16bittile)
    1475. 

  1475. {
    1476. 

  1476. 	uint16 *wr = (uint16*) r;
    1477. 

  1477. 	uint16 *wg = (uint16*) g;
    1478. 

  1478. 	uint16 *wb = (uint16*) b;
    1479. 

  1479. 	uint16 *wa = (uint16*) a;
    1480. 

  1480. 	(void) img; (void) y;
    1481. 

  1481. 	while (h-- > 0) {
    1482. 

  1482. 		for (x = 0; x < w; x++)
    1483. 

  1483.                     *cp++ = PACKW4(*wr++,*wg++,*wb++,*wa++);
    1484. 

  1484. 		SKEW4(wr, wg, wb, wa, fromskew);
    1485. 

  1485. 		cp += toskew;
    1486. 

  1486. 	}
    1487. 

  1487. }
    1488. 

  1488. 

  1489. /*
    1490. 

  1490.  * 16-bit unpacked samples => RGBA w/ unassociated alpha
    1491. 

  1491.  */
    1492. 

  1492. DECLARESepPutFunc(putRGBUAseparate16bittile)
    1493. 

  1493. {
    1494. 

  1494. 	uint16 *wr = (uint16*) r;
    1495. 

  1495. 	uint16 *wg = (uint16*) g;
    1496. 

  1496. 	uint16 *wb = (uint16*) b;
    1497. 

  1497. 	uint16 *wa = (uint16*) a;
    1498. 

  1498. 	(void) img; (void) y;
    1499. 

  1499. 	while (h-- > 0) {
    1500. 

  1500. 		uint32 r,g,b,a;
    1501. 

  1501. 		for (x = w; x-- > 0;) {
    1502. 

  1502.                     a = W2B(*wa++);
    1503. 

  1503.                     r = (a*W2B(*wr++) + 127) / 255;
    1504. 

  1504.                     g = (a*W2B(*wg++) + 127) / 255;
    1505. 

  1505.                     b = (a*W2B(*wb++) + 127) / 255;
    1506. 

  1506.                     *cp++ = PACK4(r,g,b,a);
    1507. 

  1507. 		}
    1508. 

  1508. 		SKEW4(wr, wg, wb, wa, fromskew);
    1509. 

  1509. 		cp += toskew;
    1510. 

  1510. 	}
    1511. 

  1511. }
    1512. 

  1512. 

  1513. /*
    1514. 

  1514.  * 8-bit packed CIE L*a*b 1976 samples => RGB
    1515. 

  1515.  */
    1516. 

  1516. DECLAREContigPutFunc(putcontig8bitCIELab)
    1517. 

  1517. {
    1518. 

  1518. 	float X, Y, Z;
    1519. 

  1519. 	uint32 r, g, b;
    1520. 

  1520. 	(void) y;
    1521. 

  1521. 	fromskew *= 3;
    1522. 

  1522. 	while (h-- > 0) {
    1523. 

  1523. 		for (x = w; x-- > 0;) {
    1524. 

  1524. 			TIFFCIELabToXYZ(img->cielab,
    1525. 

  1525. 					(unsigned char)pp[0],
    1526. 

  1526. 					(signed char)pp[1],
    1527. 

  1527. 					(signed char)pp[2],
    1528. 

  1528. 					&X, &Y, &Z);
    1529. 

  1529. 			TIFFXYZToRGB(img->cielab, X, Y, Z, &r, &g, &b);
    1530. 

  1530. 			*cp++ = PACK(r, g, b);
    1531. 

  1531. 			pp += 3;
    1532. 

  1532. 		}
    1533. 

  1533. 		cp += toskew;
    1534. 

  1534. 		pp += fromskew;
    1535. 

  1535. 	}
    1536. 

  1536. }
    1537. 

  1537. 

  1538. /*
    1539. 

  1539.  * YCbCr -> RGB conversion and packing routines.
    1540. 

  1540.  */
    1541. 

  1541. 

  1542. #define	YCbCrtoRGB(dst, Y) {						\
    1543. 

  1543. 	uint32 r, g, b;							\
    1544. 

  1544. 	TIFFYCbCrtoRGB(img->ycbcr, (Y), Cb, Cr, &r, &g, &b);		\
    1545. 

  1545. 	dst = PACK(r, g, b);						\
    1546. 

  1546. }
    1547. 

  1547. 

  1548. /*
    1549. 

  1549.  * 8-bit packed YCbCr samples => RGB 
    1550. 

  1550.  * This function is generic for different sampling sizes, 
    1551. 

  1551.  * and can handle blocks sizes that aren't multiples of the
    1552. 

  1552.  * sampling size.  However, it is substantially less optimized
    1553. 

  1553.  * than the specific sampling cases.  It is used as a fallback
    1554. 

  1554.  * for difficult blocks.
    1555. 

  1555.  */
    1556. 

  1556. #ifdef notdef
    1557. 

  1557. static void putcontig8bitYCbCrGenericTile( 
    1558. 

  1558.     TIFFRGBAImage* img, 
    1559. 

  1559.     uint32* cp, 
    1560. 

  1560.     uint32 x, uint32 y, 
    1561. 

  1561.     uint32 w, uint32 h, 
    1562. 

  1562.     int32 fromskew, int32 toskew, 
    1563. 

  1563.     unsigned char* pp,
    1564. 

  1564.     int h_group, 
    1565. 

  1565.     int v_group )
    1566. 

  1566. 

  1567. {
    1568. 

  1568.     uint32* cp1 = cp+w+toskew;
    1569. 

  1569.     uint32* cp2 = cp1+w+toskew;
    1570. 

  1570.     uint32* cp3 = cp2+w+toskew;
    1571. 

  1571.     int32 incr = 3*w+4*toskew;
    1572. 

  1572.     int32   Cb, Cr;
    1573. 

  1573.     int     group_size = v_group * h_group + 2;
    1574. 

  1574. 

  1575.     (void) y;
    1576. 

  1576.     fromskew = (fromskew * group_size) / h_group;
    1577. 

  1577. 

  1578.     for( yy = 0; yy < h; yy++ )
    1579. 

  1579.     {
    1580. 

  1580.         unsigned char *pp_line;
    1581. 

  1581.         int     y_line_group = yy / v_group;
    1582. 

  1582.         int     y_remainder = yy - y_line_group * v_group;
    1583. 

  1583. 

  1584.         pp_line = pp + v_line_group * 
    1585. 

  1585. 

  1586.         
    1587. 

  1587.         for( xx = 0; xx < w; xx++ )
    1588. 

  1588.         {
    1589. 

  1589.             Cb = pp
    1590. 

  1590.         }
    1591. 

  1591.     }
    1592. 

  1592.     for (; h >= 4; h -= 4) {
    1593. 

  1593. 	x = w>>2;
    1594. 

  1594. 	do {
    1595. 

  1595. 	    Cb = pp[16];
    1596. 

  1596. 	    Cr = pp[17];
    1597. 

  1597. 

  1598. 	    YCbCrtoRGB(cp [0], pp[ 0]);
    1599. 

  1599. 	    YCbCrtoRGB(cp [1], pp[ 1]);
    1600. 

  1600. 	    YCbCrtoRGB(cp [2], pp[ 2]);
    1601. 

  1601. 	    YCbCrtoRGB(cp [3], pp[ 3]);
    1602. 

  1602. 	    YCbCrtoRGB(cp1[0], pp[ 4]);
    1603. 

  1603. 	    YCbCrtoRGB(cp1[1], pp[ 5]);
    1604. 

  1604. 	    YCbCrtoRGB(cp1[2], pp[ 6]);
    1605. 

  1605. 	    YCbCrtoRGB(cp1[3], pp[ 7]);
    1606. 

  1606. 	    YCbCrtoRGB(cp2[0], pp[ 8]);
    1607. 

  1607. 	    YCbCrtoRGB(cp2[1], pp[ 9]);
    1608. 

  1608. 	    YCbCrtoRGB(cp2[2], pp[10]);
    1609. 

  1609. 	    YCbCrtoRGB(cp2[3], pp[11]);
    1610. 

  1610. 	    YCbCrtoRGB(cp3[0], pp[12]);
    1611. 

  1611. 	    YCbCrtoRGB(cp3[1], pp[13]);
    1612. 

  1612. 	    YCbCrtoRGB(cp3[2], pp[14]);
    1613. 

  1613. 	    YCbCrtoRGB(cp3[3], pp[15]);
    1614. 

  1614. 

  1615. 	    cp += 4, cp1 += 4, cp2 += 4, cp3 += 4;
    1616. 

  1616. 	    pp += 18;
    1617. 

  1617. 	} while (--x);
    1618. 

  1618. 	cp += incr, cp1 += incr, cp2 += incr, cp3 += incr;
    1619. 

  1619. 	pp += fromskew;
    1620. 

  1620.     }
    1621. 

  1621. }
    1622. 

  1622. #endif
    1623. 

  1623. 

  1624. /*
    1625. 

  1625.  * 8-bit packed YCbCr samples w/ 4,4 subsampling => RGB
    1626. 

  1626.  */
    1627. 

  1627. DECLAREContigPutFunc(putcontig8bitYCbCr44tile)
    1628. 

  1628. {
    1629. 

  1629.     uint32* cp1 = cp+w+toskew;
    1630. 

  1630.     uint32* cp2 = cp1+w+toskew;
    1631. 

  1631.     uint32* cp3 = cp2+w+toskew;
    1632. 

  1632.     int32 incr = 3*w+4*toskew;
    1633. 

  1633. 

  1634.     (void) y;
    1635. 

  1635.     /* adjust fromskew */
    1636. 

  1636.     fromskew = (fromskew * 18) / 4;
    1637. 

  1637.     if ((h & 3) == 0 && (w & 3) == 0) {				        
    1638. 

  1638.         for (; h >= 4; h -= 4) {
    1639. 

  1639.             x = w>>2;
    1640. 

  1640.             do {
    1641. 

  1641.                 int32 Cb = pp[16];
    1642. 

  1642.                 int32 Cr = pp[17];
    1643. 

  1643. 

  1644.                 YCbCrtoRGB(cp [0], pp[ 0]);
    1645. 

  1645.                 YCbCrtoRGB(cp [1], pp[ 1]);
    1646. 

  1646.                 YCbCrtoRGB(cp [2], pp[ 2]);
    1647. 

  1647.                 YCbCrtoRGB(cp [3], pp[ 3]);
    1648. 

  1648.                 YCbCrtoRGB(cp1[0], pp[ 4]);
    1649. 

  1649.                 YCbCrtoRGB(cp1[1], pp[ 5]);
    1650. 

  1650.                 YCbCrtoRGB(cp1[2], pp[ 6]);
    1651. 

  1651.                 YCbCrtoRGB(cp1[3], pp[ 7]);
    1652. 

  1652.                 YCbCrtoRGB(cp2[0], pp[ 8]);
    1653. 

  1653.                 YCbCrtoRGB(cp2[1], pp[ 9]);
    1654. 

  1654.                 YCbCrtoRGB(cp2[2], pp[10]);
    1655. 

  1655.                 YCbCrtoRGB(cp2[3], pp[11]);
    1656. 

  1656.                 YCbCrtoRGB(cp3[0], pp[12]);
    1657. 

  1657.                 YCbCrtoRGB(cp3[1], pp[13]);
    1658. 

  1658.                 YCbCrtoRGB(cp3[2], pp[14]);
    1659. 

  1659.                 YCbCrtoRGB(cp3[3], pp[15]);
    1660. 

  1660. 

  1661.                 cp += 4, cp1 += 4, cp2 += 4, cp3 += 4;
    1662. 

  1662.                 pp += 18;
    1663. 

  1663.             } while (--x);
    1664. 

  1664.             cp += incr, cp1 += incr, cp2 += incr, cp3 += incr;
    1665. 

  1665.             pp += fromskew;
    1666. 

  1666.         }
    1667. 

  1667.     } else {
    1668. 

  1668.         while (h > 0) {
    1669. 

  1669.             for (x = w; x > 0;) {
    1670. 

  1670.                 int32 Cb = pp[16];
    1671. 

  1671.                 int32 Cr = pp[17];
    1672. 

  1672.                 switch (x) {
    1673. 

  1673.                 default:
    1674. 

  1674.                     switch (h) {
    1675. 

  1675.                     default: YCbCrtoRGB(cp3[3], pp[15]); /* FALLTHROUGH */
    1676. 

  1676.                     case 3:  YCbCrtoRGB(cp2[3], pp[11]); /* FALLTHROUGH */
    1677. 

  1677.                     case 2:  YCbCrtoRGB(cp1[3], pp[ 7]); /* FALLTHROUGH */
    1678. 

  1678.                     case 1:  YCbCrtoRGB(cp [3], pp[ 3]); /* FALLTHROUGH */
    1679. 

  1679.                     }                                    /* FALLTHROUGH */
    1680. 

  1680.                 case 3:
    1681. 

  1681.                     switch (h) {
    1682. 

  1682.                     default: YCbCrtoRGB(cp3[2], pp[14]); /* FALLTHROUGH */
    1683. 

  1683.                     case 3:  YCbCrtoRGB(cp2[2], pp[10]); /* FALLTHROUGH */
    1684. 

  1684.                     case 2:  YCbCrtoRGB(cp1[2], pp[ 6]); /* FALLTHROUGH */
    1685. 

  1685.                     case 1:  YCbCrtoRGB(cp [2], pp[ 2]); /* FALLTHROUGH */
    1686. 

  1686.                     }                                    /* FALLTHROUGH */
    1687. 

  1687.                 case 2:
    1688. 

  1688.                     switch (h) {
    1689. 

  1689.                     default: YCbCrtoRGB(cp3[1], pp[13]); /* FALLTHROUGH */
    1690. 

  1690.                     case 3:  YCbCrtoRGB(cp2[1], pp[ 9]); /* FALLTHROUGH */
    1691. 

  1691.                     case 2:  YCbCrtoRGB(cp1[1], pp[ 5]); /* FALLTHROUGH */
    1692. 

  1692.                     case 1:  YCbCrtoRGB(cp [1], pp[ 1]); /* FALLTHROUGH */
    1693. 

  1693.                     }                                    /* FALLTHROUGH */
    1694. 

  1694.                 case 1:
    1695. 

  1695.                     switch (h) {
    1696. 

  1696.                     default: YCbCrtoRGB(cp3[0], pp[12]); /* FALLTHROUGH */
    1697. 

  1697.                     case 3:  YCbCrtoRGB(cp2[0], pp[ 8]); /* FALLTHROUGH */
    1698. 

  1698.                     case 2:  YCbCrtoRGB(cp1[0], pp[ 4]); /* FALLTHROUGH */
    1699. 

  1699.                     case 1:  YCbCrtoRGB(cp [0], pp[ 0]); /* FALLTHROUGH */
    1700. 

  1700.                     }                                    /* FALLTHROUGH */
    1701. 

  1701.                 }
    1702. 

  1702.                 if (x < 4) {
    1703. 

  1703.                     cp += x; cp1 += x; cp2 += x; cp3 += x;
    1704. 

  1704.                     x = 0;
    1705. 

  1705.                 }
    1706. 

  1706.                 else {
    1707. 

  1707.                     cp += 4; cp1 += 4; cp2 += 4; cp3 += 4;
    1708. 

  1708.                     x -= 4;
    1709. 

  1709.                 }
    1710. 

  1710.                 pp += 18;
    1711. 

  1711.             }
    1712. 

  1712.             if (h <= 4)
    1713. 

  1713.                 break;
    1714. 

  1714.             h -= 4;
    1715. 

  1715.             cp += incr, cp1 += incr, cp2 += incr, cp3 += incr;
    1716. 

  1716.             pp += fromskew;
    1717. 

  1717.         }
    1718. 

  1718.     }
    1719. 

  1719. }
    1720. 

  1720. 

  1721. /*
    1722. 

  1722.  * 8-bit packed YCbCr samples w/ 4,2 subsampling => RGB
    1723. 

  1723.  */
    1724. 

  1724. DECLAREContigPutFunc(putcontig8bitYCbCr42tile)
    1725. 

  1725. {
    1726. 

  1726.     uint32* cp1 = cp+w+toskew;
    1727. 

  1727.     int32 incr = 2*toskew+w;
    1728. 

  1728. 

  1729.     (void) y;
    1730. 

  1730.     fromskew = (fromskew * 10) / 4;
    1731. 

  1731.     if ((h & 3) == 0 && (w & 1) == 0) {
    1732. 

  1732.         for (; h >= 2; h -= 2) {
    1733. 

  1733.             x = w>>2;
    1734. 

  1734.             do {
    1735. 

  1735.                 int32 Cb = pp[8];
    1736. 

  1736.                 int32 Cr = pp[9];
    1737. 

  1737.                 
    1738. 

  1738.                 YCbCrtoRGB(cp [0], pp[0]);
    1739. 

  1739.                 YCbCrtoRGB(cp [1], pp[1]);
    1740. 

  1740.                 YCbCrtoRGB(cp [2], pp[2]);
    1741. 

  1741.                 YCbCrtoRGB(cp [3], pp[3]);
    1742. 

  1742.                 YCbCrtoRGB(cp1[0], pp[4]);
    1743. 

  1743.                 YCbCrtoRGB(cp1[1], pp[5]);
    1744. 

  1744.                 YCbCrtoRGB(cp1[2], pp[6]);
    1745. 

  1745.                 YCbCrtoRGB(cp1[3], pp[7]);
    1746. 

  1746.                 
    1747. 

  1747.                 cp += 4, cp1 += 4;
    1748. 

  1748.                 pp += 10;
    1749. 

  1749.             } while (--x);
    1750. 

  1750.             cp += incr, cp1 += incr;
    1751. 

  1751.             pp += fromskew;
    1752. 

  1752.         }
    1753. 

  1753.     } else {
    1754. 

  1754.         while (h > 0) {
    1755. 

  1755.             for (x = w; x > 0;) {
    1756. 

  1756.                 int32 Cb = pp[8];
    1757. 

  1757.                 int32 Cr = pp[9];
    1758. 

  1758.                 switch (x) {
    1759. 

  1759.                 default:
    1760. 

  1760.                     switch (h) {
    1761. 

  1761.                     default: YCbCrtoRGB(cp1[3], pp[ 7]); /* FALLTHROUGH */
    1762. 

  1762.                     case 1:  YCbCrtoRGB(cp [3], pp[ 3]); /* FALLTHROUGH */
    1763. 

  1763.                     }                                    /* FALLTHROUGH */
    1764. 

  1764.                 case 3:
    1765. 

  1765.                     switch (h) {
    1766. 

  1766.                     default: YCbCrtoRGB(cp1[2], pp[ 6]); /* FALLTHROUGH */
    1767. 

  1767.                     case 1:  YCbCrtoRGB(cp [2], pp[ 2]); /* FALLTHROUGH */
    1768. 

  1768.                     }                                    /* FALLTHROUGH */
    1769. 

  1769.                 case 2:
    1770. 

  1770.                     switch (h) {
    1771. 

  1771.                     default: YCbCrtoRGB(cp1[1], pp[ 5]); /* FALLTHROUGH */
    1772. 

  1772.                     case 1:  YCbCrtoRGB(cp [1], pp[ 1]); /* FALLTHROUGH */
    1773. 

  1773.                     }                                    /* FALLTHROUGH */
    1774. 

  1774.                 case 1:
    1775. 

  1775.                     switch (h) {
    1776. 

  1776.                     default: YCbCrtoRGB(cp1[0], pp[ 4]); /* FALLTHROUGH */
    1777. 

  1777.                     case 1:  YCbCrtoRGB(cp [0], pp[ 0]); /* FALLTHROUGH */
    1778. 

  1778.                     }                                    /* FALLTHROUGH */
    1779. 

  1779.                 }
    1780. 

  1780.                 if (x < 4) {
    1781. 

  1781.                     cp += x; cp1 += x;
    1782. 

  1782.                     x = 0;
    1783. 

  1783.                 }
    1784. 

  1784.                 else {
    1785. 

  1785.                     cp += 4; cp1 += 4;
    1786. 

  1786.                     x -= 4;
    1787. 

  1787.                 }
    1788. 

  1788.                 pp += 10;
    1789. 

  1789.             }
    1790. 

  1790.             if (h <= 2)
    1791. 

  1791.                 break;
    1792. 

  1792.             h -= 2;
    1793. 

  1793.             cp += incr, cp1 += incr;
    1794. 

  1794.             pp += fromskew;
    1795. 

  1795.         }
    1796. 

  1796.     }
    1797. 

  1797. }
    1798. 

  1798. 

  1799. /*
    1800. 

  1800.  * 8-bit packed YCbCr samples w/ 4,1 subsampling => RGB
    1801. 

  1801.  */
    1802. 

  1802. DECLAREContigPutFunc(putcontig8bitYCbCr41tile)
    1803. 

  1803. {
    1804. 

  1804.     (void) y;
    1805. 

  1805.     /* XXX adjust fromskew */
    1806. 

  1806.     do {
    1807. 

  1807. 	x = w>>2;
    1808. 

  1808. 	do {
    1809. 

  1809. 	    int32 Cb = pp[4];
    1810. 

  1810. 	    int32 Cr = pp[5];
    1811. 

  1811. 

  1812. 	    YCbCrtoRGB(cp [0], pp[0]);
    1813. 

  1813. 	    YCbCrtoRGB(cp [1], pp[1]);
    1814. 

  1814. 	    YCbCrtoRGB(cp [2], pp[2]);
    1815. 

  1815. 	    YCbCrtoRGB(cp [3], pp[3]);
    1816. 

  1816. 

  1817. 	    cp += 4;
    1818. 

  1818. 	    pp += 6;
    1819. 

  1819. 	} while (--x);
    1820. 

  1820. 

  1821.         if( (w&3) != 0 )
    1822. 

  1822.         {
    1823. 

  1823. 	    int32 Cb = pp[4];
    1824. 

  1824. 	    int32 Cr = pp[5];
    1825. 

  1825. 

  1826.             switch( (w&3) ) {
    1827. 

  1827.               case 3: YCbCrtoRGB(cp [2], pp[2]);
    1828. 

  1828.               case 2: YCbCrtoRGB(cp [1], pp[1]);
    1829. 

  1829.               case 1: YCbCrtoRGB(cp [0], pp[0]);
    1830. 

  1830.               case 0: break;
    1831. 

  1831.             }
    1832. 

  1832. 

  1833.             cp += (w&3);
    1834. 

  1834.             pp += 6;
    1835. 

  1835.         }
    1836. 

  1836. 

  1837. 	cp += toskew;
    1838. 

  1838. 	pp += fromskew;
    1839. 

  1839.     } while (--h);
    1840. 

  1840. 

  1841. }
    1842. 

  1842. 

  1843. /*
    1844. 

  1844.  * 8-bit packed YCbCr samples w/ 2,2 subsampling => RGB
    1845. 

  1845.  */
    1846. 

  1846. DECLAREContigPutFunc(putcontig8bitYCbCr22tile)
    1847. 

  1847. {
    1848. 

  1848. 	uint32* cp2;
    1849. 

  1849. 	int32 incr = 2*toskew+w;
    1850. 

  1850. 	(void) y;
    1851. 

  1851. 	fromskew = (fromskew / 2) * 6;
    1852. 

  1852. 	cp2 = cp+w+toskew;
    1853. 

  1853. 	while (h>=2) {
    1854. 

  1854. 		x = w;
    1855. 

  1855. 		while (x>=2) {
    1856. 

  1856. 			uint32 Cb = pp[4];
    1857. 

  1857. 			uint32 Cr = pp[5];
    1858. 

  1858. 			YCbCrtoRGB(cp[0], pp[0]);
    1859. 

  1859. 			YCbCrtoRGB(cp[1], pp[1]);
    1860. 

  1860. 			YCbCrtoRGB(cp2[0], pp[2]);
    1861. 

  1861. 			YCbCrtoRGB(cp2[1], pp[3]);
    1862. 

  1862. 			cp += 2;
    1863. 

  1863. 			cp2 += 2;
    1864. 

  1864. 			pp += 6;
    1865. 

  1865. 			x -= 2;
    1866. 

  1866. 		}
    1867. 

  1867. 		if (x==1) {
    1868. 

  1868. 			uint32 Cb = pp[4];
    1869. 

  1869. 			uint32 Cr = pp[5];
    1870. 

  1870. 			YCbCrtoRGB(cp[0], pp[0]);
    1871. 

  1871. 			YCbCrtoRGB(cp2[0], pp[2]);
    1872. 

  1872. 			cp ++ ;
    1873. 

  1873. 			cp2 ++ ;
    1874. 

  1874. 			pp += 6;
    1875. 

  1875. 		}
    1876. 

  1876. 		cp += incr;
    1877. 

  1877. 		cp2 += incr;
    1878. 

  1878. 		pp += fromskew;
    1879. 

  1879. 		h-=2;
    1880. 

  1880. 	}
    1881. 

  1881. 	if (h==1) {
    1882. 

  1882. 		x = w;
    1883. 

  1883. 		while (x>=2) {
    1884. 

  1884. 			uint32 Cb = pp[4];
    1885. 

  1885. 			uint32 Cr = pp[5];
    1886. 

  1886. 			YCbCrtoRGB(cp[0], pp[0]);
    1887. 

  1887. 			YCbCrtoRGB(cp[1], pp[1]);
    1888. 

  1888. 			cp += 2;
    1889. 

  1889. 			cp2 += 2;
    1890. 

  1890. 			pp += 6;
    1891. 

  1891. 			x -= 2;
    1892. 

  1892. 		}
    1893. 

  1893. 		if (x==1) {
    1894. 

  1894. 			uint32 Cb = pp[4];
    1895. 

  1895. 			uint32 Cr = pp[5];
    1896. 

  1896. 			YCbCrtoRGB(cp[0], pp[0]);
    1897. 

  1897. 		}
    1898. 

  1898. 	}
    1899. 

  1899. }
    1900. 

  1900. 

  1901. /*
    1902. 

  1902.  * 8-bit packed YCbCr samples w/ 2,1 subsampling => RGB
    1903. 

  1903.  */
    1904. 

  1904. DECLAREContigPutFunc(putcontig8bitYCbCr21tile)
    1905. 

  1905. {
    1906. 

  1906. 	(void) y;
    1907. 

  1907. 	fromskew = (fromskew * 4) / 2;
    1908. 

  1908. 	do {
    1909. 

  1909. 		x = w>>1;
    1910. 

  1910. 		do {
    1911. 

  1911. 			int32 Cb = pp[2];
    1912. 

  1912. 			int32 Cr = pp[3];
    1913. 

  1913. 

  1914. 			YCbCrtoRGB(cp[0], pp[0]);
    1915. 

  1915. 			YCbCrtoRGB(cp[1], pp[1]);
    1916. 

  1916. 

  1917. 			cp += 2;
    1918. 

  1918. 			pp += 4;
    1919. 

  1919. 		} while (--x);
    1920. 

  1920. 

  1921. 		if( (w&1) != 0 )
    1922. 

  1922. 		{
    1923. 

  1923. 			int32 Cb = pp[2];
    1924. 

  1924. 			int32 Cr = pp[3];
    1925. 

  1925. 

  1926. 			YCbCrtoRGB(cp[0], pp[0]);
    1927. 

  1927. 

  1928. 			cp += 1;
    1929. 

  1929. 			pp += 4;
    1930. 

  1930. 		}
    1931. 

  1931. 

  1932. 		cp += toskew;
    1933. 

  1933. 		pp += fromskew;
    1934. 

  1934. 	} while (--h);
    1935. 

  1935. }
    1936. 

  1936. 

  1937. /*
    1938. 

  1938.  * 8-bit packed YCbCr samples w/ 1,2 subsampling => RGB
    1939. 

  1939.  */
    1940. 

  1940. DECLAREContigPutFunc(putcontig8bitYCbCr12tile)
    1941. 

  1941. {
    1942. 

  1942. 	uint32* cp2;
    1943. 

  1943. 	int32 incr = 2*toskew+w;
    1944. 

  1944. 	(void) y;
    1945. 

  1945. 	fromskew = (fromskew / 2) * 4;
    1946. 

  1946. 	cp2 = cp+w+toskew;
    1947. 

  1947. 	while (h>=2) {
    1948. 

  1948. 		x = w;
    1949. 

  1949. 		do {
    1950. 

  1950. 			uint32 Cb = pp[2];
    1951. 

  1951. 			uint32 Cr = pp[3];
    1952. 

  1952. 			YCbCrtoRGB(cp[0], pp[0]);
    1953. 

  1953. 			YCbCrtoRGB(cp2[0], pp[1]);
    1954. 

  1954. 			cp ++;
    1955. 

  1955. 			cp2 ++;
    1956. 

  1956. 			pp += 4;
    1957. 

  1957. 		} while (--x);
    1958. 

  1958. 		cp += incr;
    1959. 

  1959. 		cp2 += incr;
    1960. 

  1960. 		pp += fromskew;
    1961. 

  1961. 		h-=2;
    1962. 

  1962. 	}
    1963. 

  1963. 	if (h==1) {
    1964. 

  1964. 		x = w;
    1965. 

  1965. 		do {
    1966. 

  1966. 			uint32 Cb = pp[2];
    1967. 

  1967. 			uint32 Cr = pp[3];
    1968. 

  1968. 			YCbCrtoRGB(cp[0], pp[0]);
    1969. 

  1969. 			cp ++;
    1970. 

  1970. 			pp += 4;
    1971. 

  1971. 		} while (--x);
    1972. 

  1972. 	}
    1973. 

  1973. }
    1974. 

  1974. 

  1975. /*
    1976. 

  1976.  * 8-bit packed YCbCr samples w/ no subsampling => RGB
    1977. 

  1977.  */
    1978. 

  1978. DECLAREContigPutFunc(putcontig8bitYCbCr11tile)
    1979. 

  1979. {
    1980. 

  1980. 	(void) y;
    1981. 

  1981. 	fromskew *= 3;
    1982. 

  1982. 	do {
    1983. 

  1983. 		x = w; /* was x = w>>1; patched 2000/09/25 warmerda@home.com */
    1984. 

  1984. 		do {
    1985. 

  1985. 			int32 Cb = pp[1];
    1986. 

  1986. 			int32 Cr = pp[2];
    1987. 

  1987. 

  1988. 			YCbCrtoRGB(*cp++, pp[0]);
    1989. 

  1989. 

  1990. 			pp += 3;
    1991. 

  1991. 		} while (--x);
    1992. 

  1992. 		cp += toskew;
    1993. 

  1993. 		pp += fromskew;
    1994. 

  1994. 	} while (--h);
    1995. 

  1995. }
    1996. 

  1996. 

  1997. /*
    1998. 

  1998.  * 8-bit packed YCbCr samples w/ no subsampling => RGB
    1999. 

  1999.  */
    2000. 

  2000. DECLARESepPutFunc(putseparate8bitYCbCr11tile)
    2001. 

  2001. {
    2002. 

  2002. 	(void) y;
    2003. 

  2003. 	(void) a;
    2004. 

  2004. 	/* TODO: naming of input vars is still off, change obfuscating declaration inside define, or resolve obfuscation */
    2005. 

  2005. 	while (h-- > 0) {
    2006. 

  2006. 		x = w;
    2007. 

  2007. 		do {
    2008. 

  2008. 			uint32 dr, dg, db;
    2009. 

  2009. 			TIFFYCbCrtoRGB(img->ycbcr,*r++,*g++,*b++,&dr,&dg,&db);
    2010. 

  2010. 			*cp++ = PACK(dr,dg,db);
    2011. 

  2011. 		} while (--x);
    2012. 

  2012. 		SKEW(r, g, b, fromskew);
    2013. 

  2013. 		cp += toskew;
    2014. 

  2014. 	}
    2015. 

  2015. }
    2016. 

  2016. #undef YCbCrtoRGB
    2017. 

  2017. 

  2018. static int
    2019. 

  2019. initYCbCrConversion(TIFFRGBAImage* img)
    2020. 

  2020. {
    2021. 

  2021. 	static char module[] = "initYCbCrConversion";
    2022. 

  2022. 

  2023. 	float *luma, *refBlackWhite;
    2024. 

  2024. 

  2025. 	if (img->ycbcr == NULL) {
    2026. 

  2026. 		img->ycbcr = (TIFFYCbCrToRGB*) _TIFFmalloc(
    2027. 

  2027. 		    TIFFroundup(sizeof (TIFFYCbCrToRGB), sizeof (long))
    2028. 

  2028. 		    + 4*256*sizeof (TIFFRGBValue)
    2029. 

  2029. 		    + 2*256*sizeof (int)
    2030. 

  2030. 		    + 3*256*sizeof (int32)
    2031. 

  2031. 		    );
    2032. 

  2032. 		if (img->ycbcr == NULL) {
    2033. 

  2033. 			TIFFErrorExt(img->tif->tif_clientdata, module,
    2034. 

  2034. 			    "No space for YCbCr->RGB conversion state");
    2035. 

  2035. 			return (0);
    2036. 

  2036. 		}
    2037. 

  2037. 	}
    2038. 

  2038. 

  2039. 	TIFFGetFieldDefaulted(img->tif, TIFFTAG_YCBCRCOEFFICIENTS, &luma);
    2040. 

  2040. 	TIFFGetFieldDefaulted(img->tif, TIFFTAG_REFERENCEBLACKWHITE,
    2041. 

  2041. 	    &refBlackWhite);
    2042. 

  2042. 	if (TIFFYCbCrToRGBInit(img->ycbcr, luma, refBlackWhite) < 0)
    2043. 

  2043. 		return(0);
    2044. 

  2044. 	return (1);
    2045. 

  2045. }
    2046. 

  2046. 

  2047. static tileContigRoutine
    2048. 

  2048. initCIELabConversion(TIFFRGBAImage* img)
    2049. 

  2049. {
    2050. 

  2050. 	static char module[] = "initCIELabConversion";
    2051. 

  2051. 

  2052. 	float   *whitePoint;
    2053. 

  2053. 	float   refWhite[3];
    2054. 

  2054. 

  2055. 	if (!img->cielab) {
    2056. 

  2056. 		img->cielab = (TIFFCIELabToRGB *)
    2057. 

  2057. 			_TIFFmalloc(sizeof(TIFFCIELabToRGB));
    2058. 

  2058. 		if (!img->cielab) {
    2059. 

  2059. 			TIFFErrorExt(img->tif->tif_clientdata, module,
    2060. 

  2060. 			    "No space for CIE L*a*b*->RGB conversion state.");
    2061. 

  2061. 			return NULL;
    2062. 

  2062. 		}
    2063. 

  2063. 	}
    2064. 

  2064. 

  2065. 	TIFFGetFieldDefaulted(img->tif, TIFFTAG_WHITEPOINT, &whitePoint);
    2066. 

  2066. 	refWhite[1] = 100.0F;
    2067. 

  2067. 	refWhite[0] = whitePoint[0] / whitePoint[1] * refWhite[1];
    2068. 

  2068. 	refWhite[2] = (1.0F - whitePoint[0] - whitePoint[1])
    2069. 

  2069. 		      / whitePoint[1] * refWhite[1];
    2070. 

  2070. 	if (TIFFCIELabToRGBInit(img->cielab, &display_sRGB, refWhite) < 0) {
    2071. 

  2071. 		TIFFErrorExt(img->tif->tif_clientdata, module,
    2072. 

  2072. 		    "Failed to initialize CIE L*a*b*->RGB conversion state.");
    2073. 

  2073. 		_TIFFfree(img->cielab);
    2074. 

  2074. 		return NULL;
    2075. 

  2075. 	}
    2076. 

  2076. 

  2077. 	return putcontig8bitCIELab;
    2078. 

  2078. }
    2079. 

  2079. 

  2080. /*
    2081. 

  2081.  * Greyscale images with less than 8 bits/sample are handled
    2082. 

  2082.  * with a table to avoid lots of shifts and masks.  The table
    2083. 

  2083.  * is setup so that put*bwtile (below) can retrieve 8/bitspersample
    2084. 

  2084.  * pixel values simply by indexing into the table with one
    2085. 

  2085.  * number.
    2086. 

  2086.  */
    2087. 

  2087. static int
    2088. 

  2088. makebwmap(TIFFRGBAImage* img)
    2089. 

  2089. {
    2090. 

  2090.     TIFFRGBValue* Map = img->Map;
    2091. 

  2091.     int bitspersample = img->bitspersample;
    2092. 

  2092.     int nsamples = 8 / bitspersample;
    2093. 

  2093.     int i;
    2094. 

  2094.     uint32* p;
    2095. 

  2095. 

  2096.     if( nsamples == 0 )
    2097. 

  2097.         nsamples = 1;
    2098. 

  2098. 

  2099.     img->BWmap = (uint32**) _TIFFmalloc(
    2100. 

  2100. 	256*sizeof (uint32 *)+(256*nsamples*sizeof(uint32)));
    2101. 

  2101.     if (img->BWmap == NULL) {
    2102. 

  2102. 		TIFFErrorExt(img->tif->tif_clientdata, TIFFFileName(img->tif), "No space for B&W mapping table");
    2103. 

  2103. 		return (0);
    2104. 

  2104.     }
    2105. 

  2105.     p = (uint32*)(img->BWmap + 256);
    2106. 

  2106.     for (i = 0; i < 256; i++) {
    2107. 

  2107. 	TIFFRGBValue c;
    2108. 

  2108. 	img->BWmap[i] = p;
    2109. 

  2109. 	switch (bitspersample) {
    2110. 

  2110. #define	GREY(x)	c = Map[x]; *p++ = PACK(c,c,c);
    2111. 

  2111. 	case 1:
    2112. 

  2112. 	    GREY(i>>7);
    2113. 

  2113. 	    GREY((i>>6)&1);
    2114. 

  2114. 	    GREY((i>>5)&1);
    2115. 

  2115. 	    GREY((i>>4)&1);
    2116. 

  2116. 	    GREY((i>>3)&1);
    2117. 

  2117. 	    GREY((i>>2)&1);
    2118. 

  2118. 	    GREY((i>>1)&1);
    2119. 

  2119. 	    GREY(i&1);
    2120. 

  2120. 	    break;
    2121. 

  2121. 	case 2:
    2122. 

  2122. 	    GREY(i>>6);
    2123. 

  2123. 	    GREY((i>>4)&3);
    2124. 

  2124. 	    GREY((i>>2)&3);
    2125. 

  2125. 	    GREY(i&3);
    2126. 

  2126. 	    break;
    2127. 

  2127. 	case 4:
    2128. 

  2128. 	    GREY(i>>4);
    2129. 

  2129. 	    GREY(i&0xf);
    2130. 

  2130. 	    break;
    2131. 

  2131. 	case 8:
    2132. 

  2132.         case 16:
    2133. 

  2133. 	    GREY(i);
    2134. 

  2134. 	    break;
    2135. 

  2135. 	}
    2136. 

  2136. #undef	GREY
    2137. 

  2137.     }
    2138. 

  2138.     return (1);
    2139. 

  2139. }
    2140. 

  2140. 

  2141. /*
    2142. 

  2142.  * Construct a mapping table to convert from the range
    2143. 

  2143.  * of the data samples to [0,255] --for display.  This
    2144. 

  2144.  * process also handles inverting B&W images when needed.
    2145. 

  2145.  */ 
    2146. 

  2146. static int
    2147. 

  2147. setupMap(TIFFRGBAImage* img)
    2148. 

  2148. {
    2149. 

  2149.     int32 x, range;
    2150. 

  2150. 

  2151.     range = (int32)((1L<<img->bitspersample)-1);
    2152. 

  2152.     
    2153. 

  2153.     /* treat 16 bit the same as eight bit */
    2154. 

  2154.     if( img->bitspersample == 16 )
    2155. 

  2155.         range = (int32) 255;
    2156. 

  2156. 

  2157.     img->Map = (TIFFRGBValue*) _TIFFmalloc((range+1) * sizeof (TIFFRGBValue));
    2158. 

  2158.     if (img->Map == NULL) {
    2159. 

  2159. 		TIFFErrorExt(img->tif->tif_clientdata, TIFFFileName(img->tif),
    2160. 

  2160. 			"No space for photometric conversion table");
    2161. 

  2161. 		return (0);
    2162. 

  2162.     }
    2163. 

  2163.     if (img->photometric == PHOTOMETRIC_MINISWHITE) {
    2164. 

  2164. 	for (x = 0; x <= range; x++)
    2165. 

  2165. 	    img->Map[x] = (TIFFRGBValue) (((range - x) * 255) / range);
    2166. 

  2166.     } else {
    2167. 

  2167. 	for (x = 0; x <= range; x++)
    2168. 

  2168. 	    img->Map[x] = (TIFFRGBValue) ((x * 255) / range);
    2169. 

  2169.     }
    2170. 

  2170.     if (img->bitspersample <= 16 &&
    2171. 

  2171. 	(img->photometric == PHOTOMETRIC_MINISBLACK ||
    2172. 

  2172. 	 img->photometric == PHOTOMETRIC_MINISWHITE)) {
    2173. 

  2173. 	/*
    2174. 

  2174. 	 * Use photometric mapping table to construct
    2175. 

  2175. 	 * unpacking tables for samples <= 8 bits.
    2176. 

  2176. 	 */
    2177. 

  2177. 	if (!makebwmap(img))
    2178. 

  2178. 	    return (0);
    2179. 

  2179. 	/* no longer need Map, free it */
    2180. 

  2180. 	_TIFFfree(img->Map), img->Map = NULL;
    2181. 

  2181.     }
    2182. 

  2182.     return (1);
    2183. 

  2183. }
    2184. 

  2184. 

  2185. static int
    2186. 

  2186. checkcmap(TIFFRGBAImage* img)
    2187. 

  2187. {
    2188. 

  2188.     uint16* r = img->redcmap;
    2189. 

  2189.     uint16* g = img->greencmap;
    2190. 

  2190.     uint16* b = img->bluecmap;
    2191. 

  2191.     long n = 1L<<img->bitspersample;
    2192. 

  2192. 

  2193.     while (n-- > 0)
    2194. 

  2194. 	if (*r++ >= 256 || *g++ >= 256 || *b++ >= 256)
    2195. 

  2195. 	    return (16);
    2196. 

  2196.     return (8);
    2197. 

  2197. }
    2198. 

  2198. 

  2199. static void
    2200. 

  2200. cvtcmap(TIFFRGBAImage* img)
    2201. 

  2201. {
    2202. 

  2202.     uint16* r = img->redcmap;
    2203. 

  2203.     uint16* g = img->greencmap;
    2204. 

  2204.     uint16* b = img->bluecmap;
    2205. 

  2205.     long i;
    2206. 

  2206. 

  2207.     for (i = (1L<<img->bitspersample)-1; i >= 0; i--) {
    2208. 

  2208. #define	CVT(x)		((uint16)((x)>>8))
    2209. 

  2209. 	r[i] = CVT(r[i]);
    2210. 

  2210. 	g[i] = CVT(g[i]);
    2211. 

  2211. 	b[i] = CVT(b[i]);
    2212. 

  2212. #undef	CVT
    2213. 

  2213.     }
    2214. 

  2214. }
    2215. 

  2215. 

  2216. /*
    2217. 

  2217.  * Palette images with <= 8 bits/sample are handled
    2218. 

  2218.  * with a table to avoid lots of shifts and masks.  The table
    2219. 

  2219.  * is setup so that put*cmaptile (below) can retrieve 8/bitspersample
    2220. 

  2220.  * pixel values simply by indexing into the table with one
    2221. 

  2221.  * number.
    2222. 

  2222.  */
    2223. 

  2223. static int
    2224. 

  2224. makecmap(TIFFRGBAImage* img)
    2225. 

  2225. {
    2226. 

  2226.     int bitspersample = img->bitspersample;
    2227. 

  2227.     int nsamples = 8 / bitspersample;
    2228. 

  2228.     uint16* r = img->redcmap;
    2229. 

  2229.     uint16* g = img->greencmap;
    2230. 

  2230.     uint16* b = img->bluecmap;
    2231. 

  2231.     uint32 *p;
    2232. 

  2232.     int i;
    2233. 

  2233. 

  2234.     img->PALmap = (uint32**) _TIFFmalloc(
    2235. 

  2235. 	256*sizeof (uint32 *)+(256*nsamples*sizeof(uint32)));
    2236. 

  2236.     if (img->PALmap == NULL) {
    2237. 

  2237. 		TIFFErrorExt(img->tif->tif_clientdata, TIFFFileName(img->tif), "No space for Palette mapping table");
    2238. 

  2238. 		return (0);
    2239. 

  2239. 	}
    2240. 

  2240.     p = (uint32*)(img->PALmap + 256);
    2241. 

  2241.     for (i = 0; i < 256; i++) {
    2242. 

  2242. 	TIFFRGBValue c;
    2243. 

  2243. 	img->PALmap[i] = p;
    2244. 

  2244. #define	CMAP(x)	c = (TIFFRGBValue) x; *p++ = PACK(r[c]&0xff, g[c]&0xff, b[c]&0xff);
    2245. 

  2245. 	switch (bitspersample) {
    2246. 

  2246. 	case 1:
    2247. 

  2247. 	    CMAP(i>>7);
    2248. 

  2248. 	    CMAP((i>>6)&1);
    2249. 

  2249. 	    CMAP((i>>5)&1);
    2250. 

  2250. 	    CMAP((i>>4)&1);
    2251. 

  2251. 	    CMAP((i>>3)&1);
    2252. 

  2252. 	    CMAP((i>>2)&1);
    2253. 

  2253. 	    CMAP((i>>1)&1);
    2254. 

  2254. 	    CMAP(i&1);
    2255. 

  2255. 	    break;
    2256. 

  2256. 	case 2:
    2257. 

  2257. 	    CMAP(i>>6);
    2258. 

  2258. 	    CMAP((i>>4)&3);
    2259. 

  2259. 	    CMAP((i>>2)&3);
    2260. 

  2260. 	    CMAP(i&3);
    2261. 

  2261. 	    break;
    2262. 

  2262. 	case 4:
    2263. 

  2263. 	    CMAP(i>>4);
    2264. 

  2264. 	    CMAP(i&0xf);
    2265. 

  2265. 	    break;
    2266. 

  2266. 	case 8:
    2267. 

  2267. 	    CMAP(i);
    2268. 

  2268. 	    break;
    2269. 

  2269. 	}
    2270. 

  2270. #undef CMAP
    2271. 

  2271.     }
    2272. 

  2272.     return (1);
    2273. 

  2273. }
    2274. 

  2274. 

  2275. /* 
    2276. 

  2276.  * Construct any mapping table used
    2277. 

  2277.  * by the associated put routine.
    2278. 

  2278.  */
    2279. 

  2279. static int
    2280. 

  2280. buildMap(TIFFRGBAImage* img)
    2281. 

  2281. {
    2282. 

  2282.     switch (img->photometric) {
    2283. 

  2283.     case PHOTOMETRIC_RGB:
    2284. 

  2284.     case PHOTOMETRIC_YCBCR:
    2285. 

  2285.     case PHOTOMETRIC_SEPARATED:
    2286. 

  2286. 	if (img->bitspersample == 8)
    2287. 

  2287. 	    break;
    2288. 

  2288. 	/* fall thru... */
    2289. 

  2289.     case PHOTOMETRIC_MINISBLACK:
    2290. 

  2290.     case PHOTOMETRIC_MINISWHITE:
    2291. 

  2291. 	if (!setupMap(img))
    2292. 

  2292. 	    return (0);
    2293. 

  2293. 	break;
    2294. 

  2294.     case PHOTOMETRIC_PALETTE:
    2295. 

  2295. 	/*
    2296. 

  2296. 	 * Convert 16-bit colormap to 8-bit (unless it looks
    2297. 

  2297. 	 * like an old-style 8-bit colormap).
    2298. 

  2298. 	 */
    2299. 

  2299. 	if (checkcmap(img) == 16)
    2300. 

  2300. 	    cvtcmap(img);
    2301. 

  2301. 	else
    2302. 

  2302. 	    TIFFWarningExt(img->tif->tif_clientdata, TIFFFileName(img->tif), "Assuming 8-bit colormap");
    2303. 

  2303. 	/*
    2304. 

  2304. 	 * Use mapping table and colormap to construct
    2305. 

  2305. 	 * unpacking tables for samples < 8 bits.
    2306. 

  2306. 	 */
    2307. 

  2307. 	if (img->bitspersample <= 8 && !makecmap(img))
    2308. 

  2308. 	    return (0);
    2309. 

  2309. 	break;
    2310. 

  2310.     }
    2311. 

  2311.     return (1);
    2312. 

  2312. }
    2313. 

  2313. 

  2314. /*
    2315. 

  2315.  * Select the appropriate conversion routine for packed data.
    2316. 

  2316.  */
    2317. 

  2317. static int
    2318. 

  2318. PickContigCase(TIFFRGBAImage* img)
    2319. 

  2319. {
    2320. 

  2320. 	img->get = TIFFIsTiled(img->tif) ? gtTileContig : gtStripContig;
    2321. 

  2321. 	img->put.contig = NULL;
    2322. 

  2322. 	switch (img->photometric) {
    2323. 

  2323. 		case PHOTOMETRIC_RGB:
    2324. 

  2324. 			switch (img->bitspersample) {
    2325. 

  2325. 				case 8:
    2326. 

  2326. 					if (img->alpha == EXTRASAMPLE_ASSOCALPHA)
    2327. 

  2327. 						img->put.contig = putRGBAAcontig8bittile;
    2328. 

  2328. 					else if (img->alpha == EXTRASAMPLE_UNASSALPHA)
    2329. 

  2329. 					{
    2330. 

  2330.                                             img->put.contig = putRGBUAcontig8bittile;
    2331. 

  2331. 					}
    2332. 

  2332. 					else
    2333. 

  2333.                                             img->put.contig = putRGBcontig8bittile;
    2334. 

  2334. 					break;
    2335. 

  2335. 				case 16:
    2336. 

  2336. 					if (img->alpha == EXTRASAMPLE_ASSOCALPHA)
    2337. 

  2337. 					{
    2338. 

  2338.                                             img->put.contig = putRGBAAcontig16bittile;
    2339. 

  2339. 					}
    2340. 

  2340. 					else if (img->alpha == EXTRASAMPLE_UNASSALPHA)
    2341. 

  2341. 					{
    2342. 

  2342.                                             img->put.contig = putRGBUAcontig16bittile;
    2343. 

  2343. 					}
    2344. 

  2344. 					else
    2345. 

  2345. 					{
    2346. 

  2346.                                             img->put.contig = putRGBcontig16bittile;
    2347. 

  2347. 					}
    2348. 

  2348. 					break;
    2349. 

  2349. 			}
    2350. 

  2350. 			break;
    2351. 

  2351. 		case PHOTOMETRIC_SEPARATED:
    2352. 

  2352. 			if (buildMap(img)) {
    2353. 

  2353. 				if (img->bitspersample == 8) {
    2354. 

  2354. 					if (!img->Map)
    2355. 

  2355. 						img->put.contig = putRGBcontig8bitCMYKtile;
    2356. 

  2356. 					else
    2357. 

  2357. 						img->put.contig = putRGBcontig8bitCMYKMaptile;
    2358. 

  2358. 				}
    2359. 

  2359. 			}
    2360. 

  2360. 			break;
    2361. 

  2361. 		case PHOTOMETRIC_PALETTE:
    2362. 

  2362. 			if (buildMap(img)) {
    2363. 

  2363. 				switch (img->bitspersample) {
    2364. 

  2364. 					case 8:
    2365. 

  2365. 						img->put.contig = put8bitcmaptile;
    2366. 

  2366. 						break;
    2367. 

  2367. 					case 4:
    2368. 

  2368. 						img->put.contig = put4bitcmaptile;
    2369. 

  2369. 						break;
    2370. 

  2370. 					case 2:
    2371. 

  2371. 						img->put.contig = put2bitcmaptile;
    2372. 

  2372. 						break;
    2373. 

  2373. 					case 1:
    2374. 

  2374. 						img->put.contig = put1bitcmaptile;
    2375. 

  2375. 						break;
    2376. 

  2376. 				}
    2377. 

  2377. 			}
    2378. 

  2378. 			break;
    2379. 

  2379. 		case PHOTOMETRIC_MINISWHITE:
    2380. 

  2380. 		case PHOTOMETRIC_MINISBLACK:
    2381. 

  2381. 			if (buildMap(img)) {
    2382. 

  2382. 				switch (img->bitspersample) {
    2383. 

  2383. 					case 16:
    2384. 

  2384. 						img->put.contig = put16bitbwtile;
    2385. 

  2385. 						break;
    2386. 

  2386. 					case 8:
    2387. 

  2387. 						img->put.contig = putgreytile;
    2388. 

  2388. 						break;
    2389. 

  2389. 					case 4:
    2390. 

  2390. 						img->put.contig = put4bitbwtile;
    2391. 

  2391. 						break;
    2392. 

  2392. 					case 2:
    2393. 

  2393. 						img->put.contig = put2bitbwtile;
    2394. 

  2394. 						break;
    2395. 

  2395. 					case 1:
    2396. 

  2396. 						img->put.contig = put1bitbwtile;
    2397. 

  2397. 						break;
    2398. 

  2398. 				}
    2399. 

  2399. 			}
    2400. 

  2400. 			break;
    2401. 

  2401. 		case PHOTOMETRIC_YCBCR:
    2402. 

  2402. 			if ((img->bitspersample==8) && (img->samplesperpixel==3))
    2403. 

  2403. 			{
    2404. 

  2404. 				if (initYCbCrConversion(img)!=0)
    2405. 

  2405. 				{
    2406. 

  2406. 					/*
    2407. 

  2407. 					 * The 6.0 spec says that subsampling must be
    2408. 

  2408. 					 * one of 1, 2, or 4, and that vertical subsampling
    2409. 

  2409. 					 * must always be <= horizontal subsampling; so
    2410. 

  2410. 					 * there are only a few possibilities and we just
    2411. 

  2411. 					 * enumerate the cases.
    2412. 

  2412. 					 * Joris: added support for the [1,2] case, nonetheless, to accomodate
    2413. 

  2413. 					 * some OJPEG files
    2414. 

  2414. 					 */
    2415. 

  2415. 					uint16 SubsamplingHor;
    2416. 

  2416. 					uint16 SubsamplingVer;
    2417. 

  2417. 					TIFFGetFieldDefaulted(img->tif, TIFFTAG_YCBCRSUBSAMPLING, &SubsamplingHor, &SubsamplingVer);
    2418. 

  2418. 					switch ((SubsamplingHor<<4)|SubsamplingVer) {
    2419. 

  2419. 						case 0x44:
    2420. 

  2420. 							img->put.contig = putcontig8bitYCbCr44tile;
    2421. 

  2421. 							break;
    2422. 

  2422. 						case 0x42:
    2423. 

  2423. 							img->put.contig = putcontig8bitYCbCr42tile;
    2424. 

  2424. 							break;
    2425. 

  2425. 						case 0x41:
    2426. 

  2426. 							img->put.contig = putcontig8bitYCbCr41tile;
    2427. 

  2427. 							break;
    2428. 

  2428. 						case 0x22:
    2429. 

  2429. 							img->put.contig = putcontig8bitYCbCr22tile;
    2430. 

  2430. 							break;
    2431. 

  2431. 						case 0x21:
    2432. 

  2432. 							img->put.contig = putcontig8bitYCbCr21tile;
    2433. 

  2433. 							break;
    2434. 

  2434. 						case 0x12:
    2435. 

  2435. 							img->put.contig = putcontig8bitYCbCr12tile;
    2436. 

  2436. 							break;
    2437. 

  2437. 						case 0x11:
    2438. 

  2438. 							img->put.contig = putcontig8bitYCbCr11tile;
    2439. 

  2439. 							break;
    2440. 

  2440. 					}
    2441. 

  2441. 				}
    2442. 

  2442. 			}
    2443. 

  2443. 			break;
    2444. 

  2444. 		case PHOTOMETRIC_CIELAB:
    2445. 

  2445. 			if (buildMap(img)) {
    2446. 

  2446. 				if (img->bitspersample == 8)
    2447. 

  2447. 					img->put.contig = initCIELabConversion(img);
    2448. 

  2448. 				break;
    2449. 

  2449. 			}
    2450. 

  2450. 	}
    2451. 

  2451. 	return ((img->get!=NULL) && (img->put.contig!=NULL));
    2452. 

  2452. }
    2453. 

  2453. 

  2454. /*
    2455. 

  2455.  * Select the appropriate conversion routine for unpacked data.
    2456. 

  2456.  *
    2457. 

  2457.  * NB: we assume that unpacked single channel data is directed
    2458. 

  2458.  *	 to the "packed routines.
    2459. 

  2459.  */
    2460. 

  2460. static int
    2461. 

  2461. PickSeparateCase(TIFFRGBAImage* img)
    2462. 

  2462. {
    2463. 

  2463. 	img->get = TIFFIsTiled(img->tif) ? gtTileSeparate : gtStripSeparate;
    2464. 

  2464. 	img->put.separate = NULL;
    2465. 

  2465. 	switch (img->photometric) {
    2466. 

  2466. 		case PHOTOMETRIC_RGB:
    2467. 

  2467. 			switch (img->bitspersample) {
    2468. 

  2468. 				case 8:
    2469. 

  2469. 					if (img->alpha == EXTRASAMPLE_ASSOCALPHA)
    2470. 

  2470. 						img->put.separate = putRGBAAseparate8bittile;
    2471. 

  2471. 					else if (img->alpha == EXTRASAMPLE_UNASSALPHA)
    2472. 

  2472. 					{
    2473. 

  2473.                                             img->put.separate = putRGBUAseparate8bittile;
    2474. 

  2474. 					}
    2475. 

  2475. 					else
    2476. 

  2476. 						img->put.separate = putRGBseparate8bittile;
    2477. 

  2477. 					break;
    2478. 

  2478. 				case 16:
    2479. 

  2479. 					if (img->alpha == EXTRASAMPLE_ASSOCALPHA)
    2480. 

  2480. 					{
    2481. 

  2481.                                             img->put.separate = putRGBAAseparate16bittile;
    2482. 

  2482. 					}
    2483. 

  2483. 					else if (img->alpha == EXTRASAMPLE_UNASSALPHA)
    2484. 

  2484. 					{
    2485. 

  2485.                                             img->put.separate = putRGBUAseparate16bittile;
    2486. 

  2486. 					}
    2487. 

  2487. 					else
    2488. 

  2488. 					{
    2489. 

  2489.                                             img->put.separate = putRGBseparate16bittile;
    2490. 

  2490. 					}
    2491. 

  2491. 					break;
    2492. 

  2492. 			}
    2493. 

  2493. 			break;
    2494. 

  2494. 		case PHOTOMETRIC_YCBCR:
    2495. 

  2495. 			if ((img->bitspersample==8) && (img->samplesperpixel==3))
    2496. 

  2496. 			{
    2497. 

  2497. 				if (initYCbCrConversion(img)!=0)
    2498. 

  2498. 				{
    2499. 

  2499. 					uint16 hs, vs;
    2500. 

  2500. 					TIFFGetFieldDefaulted(img->tif, TIFFTAG_YCBCRSUBSAMPLING, &hs, &vs);
    2501. 

  2501. 					switch ((hs<<4)|vs) {
    2502. 

  2502. 						case 0x11:
    2503. 

  2503. 							img->put.separate = putseparate8bitYCbCr11tile;
    2504. 

  2504. 							break;
    2505. 

  2505. 						/* TODO: add other cases here */
    2506. 

  2506. 					}
    2507. 

  2507. 				}
    2508. 

  2508. 			}
    2509. 

  2509. 			break;
    2510. 

  2510. 	}
    2511. 

  2511. 	return ((img->get!=NULL) && (img->put.separate!=NULL));
    2512. 

  2512. }
    2513. 

  2513. 

  2514. /*
    2515. 

  2515.  * Read a whole strip off data from the file, and convert to RGBA form.
    2516. 

  2516.  * If this is the last strip, then it will only contain the portion of
    2517. 

  2517.  * the strip that is actually within the image space.  The result is
    2518. 

  2518.  * organized in bottom to top form.
    2519. 

  2519.  */
    2520. 

  2520. 

  2521. 

  2522. int
    2523. 

  2523. TIFFReadRGBAStrip(TIFF* tif, uint32 row, uint32 * raster )
    2524. 

  2524. 

  2525. {
    2526. 

  2526.     char 	emsg[1024] = "";
    2527. 

  2527.     TIFFRGBAImage img;
    2528. 

  2528.     int 	ok;
    2529. 

  2529.     uint32	rowsperstrip, rows_to_read;
    2530. 

  2530. 

  2531.     if( TIFFIsTiled( tif ) )
    2532. 

  2532.     {
    2533. 

  2533. 		TIFFErrorExt(tif->tif_clientdata, TIFFFileName(tif),
    2534. 

  2534.                   "Can't use TIFFReadRGBAStrip() with tiled file.");
    2535. 

  2535. 	return (0);
    2536. 

  2536.     }
    2537. 

  2537.     
    2538. 

  2538.     TIFFGetFieldDefaulted(tif, TIFFTAG_ROWSPERSTRIP, &rowsperstrip);
    2539. 

  2539.     if( (row % rowsperstrip) != 0 )
    2540. 

  2540.     {
    2541. 

  2541. 		TIFFErrorExt(tif->tif_clientdata, TIFFFileName(tif),
    2542. 

  2542. 				"Row passed to TIFFReadRGBAStrip() must be first in a strip.");
    2543. 

  2543. 		return (0);
    2544. 

  2544.     }
    2545. 

  2545. 

  2546.     if (TIFFRGBAImageOK(tif, emsg) && TIFFRGBAImageBegin(&img, tif, 0, emsg)) {
    2547. 

  2547. 

  2548.         img.row_offset = row;
    2549. 

  2549.         img.col_offset = 0;
    2550. 

  2550. 

  2551.         if( row + rowsperstrip > img.height )
    2552. 

  2552.             rows_to_read = img.height - row;
    2553. 

  2553.         else
    2554. 

  2554.             rows_to_read 
    2555.