jdcolor.c 28 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896
  1. /*
  2. * jdcolor.c
  3. *
  4. * This file was part of the Independent JPEG Group's software:
  5. * Copyright (C) 1991-1997, Thomas G. Lane.
  6. * Modified 2011 by Guido Vollbeding.
  7. * libjpeg-turbo Modifications:
  8. * Copyright 2009 Pierre Ossman <[email protected]> for Cendio AB
  9. * Copyright (C) 2009, 2011-2012, 2014-2015, D. R. Commander.
  10. * Copyright (C) 2013, Linaro Limited.
  11. * For conditions of distribution and use, see the accompanying README file.
  12. *
  13. * This file contains output colorspace conversion routines.
  14. */
  15. #define JPEG_INTERNALS
  16. #include "jinclude.h"
  17. #include "jpeglib.h"
  18. #include "jsimd.h"
  19. #include "jconfigint.h"
  20. /* Private subobject */
  21. typedef struct {
  22. struct jpeg_color_deconverter pub; /* public fields */
  23. /* Private state for YCC->RGB conversion */
  24. int * Cr_r_tab; /* => table for Cr to R conversion */
  25. int * Cb_b_tab; /* => table for Cb to B conversion */
  26. INT32 * Cr_g_tab; /* => table for Cr to G conversion */
  27. INT32 * Cb_g_tab; /* => table for Cb to G conversion */
  28. /* Private state for RGB->Y conversion */
  29. INT32 * rgb_y_tab; /* => table for RGB to Y conversion */
  30. } my_color_deconverter;
  31. typedef my_color_deconverter * my_cconvert_ptr;
  32. /**************** YCbCr -> RGB conversion: most common case **************/
  33. /**************** RGB -> Y conversion: less common case **************/
  34. /*
  35. * YCbCr is defined per CCIR 601-1, except that Cb and Cr are
  36. * normalized to the range 0..MAXJSAMPLE rather than -0.5 .. 0.5.
  37. * The conversion equations to be implemented are therefore
  38. *
  39. * R = Y + 1.40200 * Cr
  40. * G = Y - 0.34414 * Cb - 0.71414 * Cr
  41. * B = Y + 1.77200 * Cb
  42. *
  43. * Y = 0.29900 * R + 0.58700 * G + 0.11400 * B
  44. *
  45. * where Cb and Cr represent the incoming values less CENTERJSAMPLE.
  46. * (These numbers are derived from TIFF 6.0 section 21, dated 3-June-92.)
  47. *
  48. * To avoid floating-point arithmetic, we represent the fractional constants
  49. * as integers scaled up by 2^16 (about 4 digits precision); we have to divide
  50. * the products by 2^16, with appropriate rounding, to get the correct answer.
  51. * Notice that Y, being an integral input, does not contribute any fraction
  52. * so it need not participate in the rounding.
  53. *
  54. * For even more speed, we avoid doing any multiplications in the inner loop
  55. * by precalculating the constants times Cb and Cr for all possible values.
  56. * For 8-bit JSAMPLEs this is very reasonable (only 256 entries per table);
  57. * for 12-bit samples it is still acceptable. It's not very reasonable for
  58. * 16-bit samples, but if you want lossless storage you shouldn't be changing
  59. * colorspace anyway.
  60. * The Cr=>R and Cb=>B values can be rounded to integers in advance; the
  61. * values for the G calculation are left scaled up, since we must add them
  62. * together before rounding.
  63. */
  64. #define SCALEBITS 16 /* speediest right-shift on some machines */
  65. #define ONE_HALF ((INT32) 1 << (SCALEBITS-1))
  66. #define FIX(x) ((INT32) ((x) * (1L<<SCALEBITS) + 0.5))
  67. /* We allocate one big table for RGB->Y conversion and divide it up into
  68. * three parts, instead of doing three alloc_small requests. This lets us
  69. * use a single table base address, which can be held in a register in the
  70. * inner loops on many machines (more than can hold all three addresses,
  71. * anyway).
  72. */
  73. #define R_Y_OFF 0 /* offset to R => Y section */
  74. #define G_Y_OFF (1*(MAXJSAMPLE+1)) /* offset to G => Y section */
  75. #define B_Y_OFF (2*(MAXJSAMPLE+1)) /* etc. */
  76. #define TABLE_SIZE (3*(MAXJSAMPLE+1))
  77. /* Include inline routines for colorspace extensions */
  78. #include "jdcolext.c"
  79. #undef RGB_RED
  80. #undef RGB_GREEN
  81. #undef RGB_BLUE
  82. #undef RGB_PIXELSIZE
  83. #define RGB_RED EXT_RGB_RED
  84. #define RGB_GREEN EXT_RGB_GREEN
  85. #define RGB_BLUE EXT_RGB_BLUE
  86. #define RGB_PIXELSIZE EXT_RGB_PIXELSIZE
  87. #define ycc_rgb_convert_internal ycc_extrgb_convert_internal
  88. #define gray_rgb_convert_internal gray_extrgb_convert_internal
  89. #define rgb_rgb_convert_internal rgb_extrgb_convert_internal
  90. #include "jdcolext.c"
  91. #undef RGB_RED
  92. #undef RGB_GREEN
  93. #undef RGB_BLUE
  94. #undef RGB_PIXELSIZE
  95. #undef ycc_rgb_convert_internal
  96. #undef gray_rgb_convert_internal
  97. #undef rgb_rgb_convert_internal
  98. #define RGB_RED EXT_RGBX_RED
  99. #define RGB_GREEN EXT_RGBX_GREEN
  100. #define RGB_BLUE EXT_RGBX_BLUE
  101. #define RGB_ALPHA 3
  102. #define RGB_PIXELSIZE EXT_RGBX_PIXELSIZE
  103. #define ycc_rgb_convert_internal ycc_extrgbx_convert_internal
  104. #define gray_rgb_convert_internal gray_extrgbx_convert_internal
  105. #define rgb_rgb_convert_internal rgb_extrgbx_convert_internal
  106. #include "jdcolext.c"
  107. #undef RGB_RED
  108. #undef RGB_GREEN
  109. #undef RGB_BLUE
  110. #undef RGB_ALPHA
  111. #undef RGB_PIXELSIZE
  112. #undef ycc_rgb_convert_internal
  113. #undef gray_rgb_convert_internal
  114. #undef rgb_rgb_convert_internal
  115. #define RGB_RED EXT_BGR_RED
  116. #define RGB_GREEN EXT_BGR_GREEN
  117. #define RGB_BLUE EXT_BGR_BLUE
  118. #define RGB_PIXELSIZE EXT_BGR_PIXELSIZE
  119. #define ycc_rgb_convert_internal ycc_extbgr_convert_internal
  120. #define gray_rgb_convert_internal gray_extbgr_convert_internal
  121. #define rgb_rgb_convert_internal rgb_extbgr_convert_internal
  122. #include "jdcolext.c"
  123. #undef RGB_RED
  124. #undef RGB_GREEN
  125. #undef RGB_BLUE
  126. #undef RGB_PIXELSIZE
  127. #undef ycc_rgb_convert_internal
  128. #undef gray_rgb_convert_internal
  129. #undef rgb_rgb_convert_internal
  130. #define RGB_RED EXT_BGRX_RED
  131. #define RGB_GREEN EXT_BGRX_GREEN
  132. #define RGB_BLUE EXT_BGRX_BLUE
  133. #define RGB_ALPHA 3
  134. #define RGB_PIXELSIZE EXT_BGRX_PIXELSIZE
  135. #define ycc_rgb_convert_internal ycc_extbgrx_convert_internal
  136. #define gray_rgb_convert_internal gray_extbgrx_convert_internal
  137. #define rgb_rgb_convert_internal rgb_extbgrx_convert_internal
  138. #include "jdcolext.c"
  139. #undef RGB_RED
  140. #undef RGB_GREEN
  141. #undef RGB_BLUE
  142. #undef RGB_ALPHA
  143. #undef RGB_PIXELSIZE
  144. #undef ycc_rgb_convert_internal
  145. #undef gray_rgb_convert_internal
  146. #undef rgb_rgb_convert_internal
  147. #define RGB_RED EXT_XBGR_RED
  148. #define RGB_GREEN EXT_XBGR_GREEN
  149. #define RGB_BLUE EXT_XBGR_BLUE
  150. #define RGB_ALPHA 0
  151. #define RGB_PIXELSIZE EXT_XBGR_PIXELSIZE
  152. #define ycc_rgb_convert_internal ycc_extxbgr_convert_internal
  153. #define gray_rgb_convert_internal gray_extxbgr_convert_internal
  154. #define rgb_rgb_convert_internal rgb_extxbgr_convert_internal
  155. #include "jdcolext.c"
  156. #undef RGB_RED
  157. #undef RGB_GREEN
  158. #undef RGB_BLUE
  159. #undef RGB_ALPHA
  160. #undef RGB_PIXELSIZE
  161. #undef ycc_rgb_convert_internal
  162. #undef gray_rgb_convert_internal
  163. #undef rgb_rgb_convert_internal
  164. #define RGB_RED EXT_XRGB_RED
  165. #define RGB_GREEN EXT_XRGB_GREEN
  166. #define RGB_BLUE EXT_XRGB_BLUE
  167. #define RGB_ALPHA 0
  168. #define RGB_PIXELSIZE EXT_XRGB_PIXELSIZE
  169. #define ycc_rgb_convert_internal ycc_extxrgb_convert_internal
  170. #define gray_rgb_convert_internal gray_extxrgb_convert_internal
  171. #define rgb_rgb_convert_internal rgb_extxrgb_convert_internal
  172. #include "jdcolext.c"
  173. #undef RGB_RED
  174. #undef RGB_GREEN
  175. #undef RGB_BLUE
  176. #undef RGB_ALPHA
  177. #undef RGB_PIXELSIZE
  178. #undef ycc_rgb_convert_internal
  179. #undef gray_rgb_convert_internal
  180. #undef rgb_rgb_convert_internal
  181. /*
  182. * Initialize tables for YCC->RGB colorspace conversion.
  183. */
  184. LOCAL(void)
  185. build_ycc_rgb_table (j_decompress_ptr cinfo)
  186. {
  187. my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
  188. int i;
  189. INT32 x;
  190. SHIFT_TEMPS
  191. cconvert->Cr_r_tab = (int *)
  192. (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
  193. (MAXJSAMPLE+1) * sizeof(int));
  194. cconvert->Cb_b_tab = (int *)
  195. (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
  196. (MAXJSAMPLE+1) * sizeof(int));
  197. cconvert->Cr_g_tab = (INT32 *)
  198. (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
  199. (MAXJSAMPLE+1) * sizeof(INT32));
  200. cconvert->Cb_g_tab = (INT32 *)
  201. (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
  202. (MAXJSAMPLE+1) * sizeof(INT32));
  203. for (i = 0, x = -CENTERJSAMPLE; i <= MAXJSAMPLE; i++, x++) {
  204. /* i is the actual input pixel value, in the range 0..MAXJSAMPLE */
  205. /* The Cb or Cr value we are thinking of is x = i - CENTERJSAMPLE */
  206. /* Cr=>R value is nearest int to 1.40200 * x */
  207. cconvert->Cr_r_tab[i] = (int)
  208. RIGHT_SHIFT(FIX(1.40200) * x + ONE_HALF, SCALEBITS);
  209. /* Cb=>B value is nearest int to 1.77200 * x */
  210. cconvert->Cb_b_tab[i] = (int)
  211. RIGHT_SHIFT(FIX(1.77200) * x + ONE_HALF, SCALEBITS);
  212. /* Cr=>G value is scaled-up -0.71414 * x */
  213. cconvert->Cr_g_tab[i] = (- FIX(0.71414)) * x;
  214. /* Cb=>G value is scaled-up -0.34414 * x */
  215. /* We also add in ONE_HALF so that need not do it in inner loop */
  216. cconvert->Cb_g_tab[i] = (- FIX(0.34414)) * x + ONE_HALF;
  217. }
  218. }
  219. /*
  220. * Convert some rows of samples to the output colorspace.
  221. */
  222. METHODDEF(void)
  223. ycc_rgb_convert (j_decompress_ptr cinfo,
  224. JSAMPIMAGE input_buf, JDIMENSION input_row,
  225. JSAMPARRAY output_buf, int num_rows)
  226. {
  227. switch (cinfo->out_color_space) {
  228. case JCS_EXT_RGB:
  229. ycc_extrgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  230. num_rows);
  231. break;
  232. case JCS_EXT_RGBX:
  233. case JCS_EXT_RGBA:
  234. ycc_extrgbx_convert_internal(cinfo, input_buf, input_row, output_buf,
  235. num_rows);
  236. break;
  237. case JCS_EXT_BGR:
  238. ycc_extbgr_convert_internal(cinfo, input_buf, input_row, output_buf,
  239. num_rows);
  240. break;
  241. case JCS_EXT_BGRX:
  242. case JCS_EXT_BGRA:
  243. ycc_extbgrx_convert_internal(cinfo, input_buf, input_row, output_buf,
  244. num_rows);
  245. break;
  246. case JCS_EXT_XBGR:
  247. case JCS_EXT_ABGR:
  248. ycc_extxbgr_convert_internal(cinfo, input_buf, input_row, output_buf,
  249. num_rows);
  250. break;
  251. case JCS_EXT_XRGB:
  252. case JCS_EXT_ARGB:
  253. ycc_extxrgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  254. num_rows);
  255. break;
  256. default:
  257. ycc_rgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  258. num_rows);
  259. break;
  260. }
  261. }
  262. /**************** Cases other than YCbCr -> RGB **************/
  263. /*
  264. * Initialize for RGB->grayscale colorspace conversion.
  265. */
  266. LOCAL(void)
  267. build_rgb_y_table (j_decompress_ptr cinfo)
  268. {
  269. my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
  270. INT32 * rgb_y_tab;
  271. INT32 i;
  272. /* Allocate and fill in the conversion tables. */
  273. cconvert->rgb_y_tab = rgb_y_tab = (INT32 *)
  274. (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
  275. (TABLE_SIZE * sizeof(INT32)));
  276. for (i = 0; i <= MAXJSAMPLE; i++) {
  277. rgb_y_tab[i+R_Y_OFF] = FIX(0.29900) * i;
  278. rgb_y_tab[i+G_Y_OFF] = FIX(0.58700) * i;
  279. rgb_y_tab[i+B_Y_OFF] = FIX(0.11400) * i + ONE_HALF;
  280. }
  281. }
  282. /*
  283. * Convert RGB to grayscale.
  284. */
  285. METHODDEF(void)
  286. rgb_gray_convert (j_decompress_ptr cinfo,
  287. JSAMPIMAGE input_buf, JDIMENSION input_row,
  288. JSAMPARRAY output_buf, int num_rows)
  289. {
  290. my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
  291. register int r, g, b;
  292. register INT32 * ctab = cconvert->rgb_y_tab;
  293. register JSAMPROW outptr;
  294. register JSAMPROW inptr0, inptr1, inptr2;
  295. register JDIMENSION col;
  296. JDIMENSION num_cols = cinfo->output_width;
  297. while (--num_rows >= 0) {
  298. inptr0 = input_buf[0][input_row];
  299. inptr1 = input_buf[1][input_row];
  300. inptr2 = input_buf[2][input_row];
  301. input_row++;
  302. outptr = *output_buf++;
  303. for (col = 0; col < num_cols; col++) {
  304. r = GETJSAMPLE(inptr0[col]);
  305. g = GETJSAMPLE(inptr1[col]);
  306. b = GETJSAMPLE(inptr2[col]);
  307. /* Y */
  308. outptr[col] = (JSAMPLE)
  309. ((ctab[r+R_Y_OFF] + ctab[g+G_Y_OFF] + ctab[b+B_Y_OFF])
  310. >> SCALEBITS);
  311. }
  312. }
  313. }
  314. /*
  315. * Color conversion for no colorspace change: just copy the data,
  316. * converting from separate-planes to interleaved representation.
  317. */
  318. METHODDEF(void)
  319. null_convert (j_decompress_ptr cinfo,
  320. JSAMPIMAGE input_buf, JDIMENSION input_row,
  321. JSAMPARRAY output_buf, int num_rows)
  322. {
  323. register JSAMPROW inptr, inptr0, inptr1, inptr2, inptr3, outptr;
  324. register JDIMENSION col;
  325. register int num_components = cinfo->num_components;
  326. JDIMENSION num_cols = cinfo->output_width;
  327. int ci;
  328. if (num_components == 3) {
  329. while (--num_rows >= 0) {
  330. inptr0 = input_buf[0][input_row];
  331. inptr1 = input_buf[1][input_row];
  332. inptr2 = input_buf[2][input_row];
  333. input_row++;
  334. outptr = *output_buf++;
  335. for (col = 0; col < num_cols; col++) {
  336. *outptr++ = inptr0[col];
  337. *outptr++ = inptr1[col];
  338. *outptr++ = inptr2[col];
  339. }
  340. }
  341. } else if (num_components == 4) {
  342. while (--num_rows >= 0) {
  343. inptr0 = input_buf[0][input_row];
  344. inptr1 = input_buf[1][input_row];
  345. inptr2 = input_buf[2][input_row];
  346. inptr3 = input_buf[3][input_row];
  347. input_row++;
  348. outptr = *output_buf++;
  349. for (col = 0; col < num_cols; col++) {
  350. *outptr++ = inptr0[col];
  351. *outptr++ = inptr1[col];
  352. *outptr++ = inptr2[col];
  353. *outptr++ = inptr3[col];
  354. }
  355. }
  356. } else {
  357. while (--num_rows >= 0) {
  358. for (ci = 0; ci < num_components; ci++) {
  359. inptr = input_buf[ci][input_row];
  360. outptr = *output_buf;
  361. for (col = 0; col < num_cols; col++) {
  362. outptr[ci] = inptr[col];
  363. outptr += num_components;
  364. }
  365. }
  366. output_buf++;
  367. input_row++;
  368. }
  369. }
  370. }
  371. /*
  372. * Color conversion for grayscale: just copy the data.
  373. * This also works for YCbCr -> grayscale conversion, in which
  374. * we just copy the Y (luminance) component and ignore chrominance.
  375. */
  376. METHODDEF(void)
  377. grayscale_convert (j_decompress_ptr cinfo,
  378. JSAMPIMAGE input_buf, JDIMENSION input_row,
  379. JSAMPARRAY output_buf, int num_rows)
  380. {
  381. jcopy_sample_rows(input_buf[0], (int) input_row, output_buf, 0,
  382. num_rows, cinfo->output_width);
  383. }
  384. /*
  385. * Convert grayscale to RGB
  386. */
  387. METHODDEF(void)
  388. gray_rgb_convert (j_decompress_ptr cinfo,
  389. JSAMPIMAGE input_buf, JDIMENSION input_row,
  390. JSAMPARRAY output_buf, int num_rows)
  391. {
  392. switch (cinfo->out_color_space) {
  393. case JCS_EXT_RGB:
  394. gray_extrgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  395. num_rows);
  396. break;
  397. case JCS_EXT_RGBX:
  398. case JCS_EXT_RGBA:
  399. gray_extrgbx_convert_internal(cinfo, input_buf, input_row, output_buf,
  400. num_rows);
  401. break;
  402. case JCS_EXT_BGR:
  403. gray_extbgr_convert_internal(cinfo, input_buf, input_row, output_buf,
  404. num_rows);
  405. break;
  406. case JCS_EXT_BGRX:
  407. case JCS_EXT_BGRA:
  408. gray_extbgrx_convert_internal(cinfo, input_buf, input_row, output_buf,
  409. num_rows);
  410. break;
  411. case JCS_EXT_XBGR:
  412. case JCS_EXT_ABGR:
  413. gray_extxbgr_convert_internal(cinfo, input_buf, input_row, output_buf,
  414. num_rows);
  415. break;
  416. case JCS_EXT_XRGB:
  417. case JCS_EXT_ARGB:
  418. gray_extxrgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  419. num_rows);
  420. break;
  421. default:
  422. gray_rgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  423. num_rows);
  424. break;
  425. }
  426. }
  427. /*
  428. * Convert plain RGB to extended RGB
  429. */
  430. METHODDEF(void)
  431. rgb_rgb_convert (j_decompress_ptr cinfo,
  432. JSAMPIMAGE input_buf, JDIMENSION input_row,
  433. JSAMPARRAY output_buf, int num_rows)
  434. {
  435. switch (cinfo->out_color_space) {
  436. case JCS_EXT_RGB:
  437. rgb_extrgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  438. num_rows);
  439. break;
  440. case JCS_EXT_RGBX:
  441. case JCS_EXT_RGBA:
  442. rgb_extrgbx_convert_internal(cinfo, input_buf, input_row, output_buf,
  443. num_rows);
  444. break;
  445. case JCS_EXT_BGR:
  446. rgb_extbgr_convert_internal(cinfo, input_buf, input_row, output_buf,
  447. num_rows);
  448. break;
  449. case JCS_EXT_BGRX:
  450. case JCS_EXT_BGRA:
  451. rgb_extbgrx_convert_internal(cinfo, input_buf, input_row, output_buf,
  452. num_rows);
  453. break;
  454. case JCS_EXT_XBGR:
  455. case JCS_EXT_ABGR:
  456. rgb_extxbgr_convert_internal(cinfo, input_buf, input_row, output_buf,
  457. num_rows);
  458. break;
  459. case JCS_EXT_XRGB:
  460. case JCS_EXT_ARGB:
  461. rgb_extxrgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  462. num_rows);
  463. break;
  464. default:
  465. rgb_rgb_convert_internal(cinfo, input_buf, input_row, output_buf,
  466. num_rows);
  467. break;
  468. }
  469. }
  470. /*
  471. * Adobe-style YCCK->CMYK conversion.
  472. * We convert YCbCr to R=1-C, G=1-M, and B=1-Y using the same
  473. * conversion as above, while passing K (black) unchanged.
  474. * We assume build_ycc_rgb_table has been called.
  475. */
  476. METHODDEF(void)
  477. ycck_cmyk_convert (j_decompress_ptr cinfo,
  478. JSAMPIMAGE input_buf, JDIMENSION input_row,
  479. JSAMPARRAY output_buf, int num_rows)
  480. {
  481. my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
  482. register int y, cb, cr;
  483. register JSAMPROW outptr;
  484. register JSAMPROW inptr0, inptr1, inptr2, inptr3;
  485. register JDIMENSION col;
  486. JDIMENSION num_cols = cinfo->output_width;
  487. /* copy these pointers into registers if possible */
  488. register JSAMPLE * range_limit = cinfo->sample_range_limit;
  489. register int * Crrtab = cconvert->Cr_r_tab;
  490. register int * Cbbtab = cconvert->Cb_b_tab;
  491. register INT32 * Crgtab = cconvert->Cr_g_tab;
  492. register INT32 * Cbgtab = cconvert->Cb_g_tab;
  493. SHIFT_TEMPS
  494. while (--num_rows >= 0) {
  495. inptr0 = input_buf[0][input_row];
  496. inptr1 = input_buf[1][input_row];
  497. inptr2 = input_buf[2][input_row];
  498. inptr3 = input_buf[3][input_row];
  499. input_row++;
  500. outptr = *output_buf++;
  501. for (col = 0; col < num_cols; col++) {
  502. y = GETJSAMPLE(inptr0[col]);
  503. cb = GETJSAMPLE(inptr1[col]);
  504. cr = GETJSAMPLE(inptr2[col]);
  505. /* Range-limiting is essential due to noise introduced by DCT losses. */
  506. outptr[0] = range_limit[MAXJSAMPLE - (y + Crrtab[cr])]; /* red */
  507. outptr[1] = range_limit[MAXJSAMPLE - (y + /* green */
  508. ((int) RIGHT_SHIFT(Cbgtab[cb] + Crgtab[cr],
  509. SCALEBITS)))];
  510. outptr[2] = range_limit[MAXJSAMPLE - (y + Cbbtab[cb])]; /* blue */
  511. /* K passes through unchanged */
  512. outptr[3] = inptr3[col]; /* don't need GETJSAMPLE here */
  513. outptr += 4;
  514. }
  515. }
  516. }
  517. /*
  518. * RGB565 conversion
  519. */
  520. #define PACK_SHORT_565_LE(r, g, b) ((((r) << 8) & 0xF800) | \
  521. (((g) << 3) & 0x7E0) | ((b) >> 3))
  522. #define PACK_SHORT_565_BE(r, g, b) (((r) & 0xF8) | ((g) >> 5) | \
  523. (((g) << 11) & 0xE000) | \
  524. (((b) << 5) & 0x1F00))
  525. #define PACK_TWO_PIXELS_LE(l, r) ((r << 16) | l)
  526. #define PACK_TWO_PIXELS_BE(l, r) ((l << 16) | r)
  527. #define PACK_NEED_ALIGNMENT(ptr) (((size_t)(ptr)) & 3)
  528. #define WRITE_TWO_ALIGNED_PIXELS(addr, pixels) ((*(int *)(addr)) = pixels)
  529. #define DITHER_565_R(r, dither) ((r) + ((dither) & 0xFF))
  530. #define DITHER_565_G(g, dither) ((g) + (((dither) & 0xFF) >> 1))
  531. #define DITHER_565_B(b, dither) ((b) + ((dither) & 0xFF))
  532. /* Declarations for ordered dithering
  533. *
  534. * We use a 4x4 ordered dither array packed into 32 bits. This array is
  535. * sufficent for dithering RGB888 to RGB565.
  536. */
  537. #define DITHER_MASK 0x3
  538. #define DITHER_ROTATE(x) (((x) << 24) | (((x) >> 8) & 0x00FFFFFF))
  539. static const INT32 dither_matrix[4] = {
  540. 0x0008020A,
  541. 0x0C040E06,
  542. 0x030B0109,
  543. 0x0F070D05
  544. };
  545. static INLINE boolean is_big_endian(void)
  546. {
  547. int test_value = 1;
  548. if(*(char *)&test_value != 1)
  549. return TRUE;
  550. return FALSE;
  551. }
  552. /* Include inline routines for RGB565 conversion */
  553. #define PACK_SHORT_565 PACK_SHORT_565_LE
  554. #define PACK_TWO_PIXELS PACK_TWO_PIXELS_LE
  555. #define ycc_rgb565_convert_internal ycc_rgb565_convert_le
  556. #define ycc_rgb565D_convert_internal ycc_rgb565D_convert_le
  557. #define rgb_rgb565_convert_internal rgb_rgb565_convert_le
  558. #define rgb_rgb565D_convert_internal rgb_rgb565D_convert_le
  559. #define gray_rgb565_convert_internal gray_rgb565_convert_le
  560. #define gray_rgb565D_convert_internal gray_rgb565D_convert_le
  561. #include "jdcol565.c"
  562. #undef PACK_SHORT_565
  563. #undef PACK_TWO_PIXELS
  564. #undef ycc_rgb565_convert_internal
  565. #undef ycc_rgb565D_convert_internal
  566. #undef rgb_rgb565_convert_internal
  567. #undef rgb_rgb565D_convert_internal
  568. #undef gray_rgb565_convert_internal
  569. #undef gray_rgb565D_convert_internal
  570. #define PACK_SHORT_565 PACK_SHORT_565_BE
  571. #define PACK_TWO_PIXELS PACK_TWO_PIXELS_BE
  572. #define ycc_rgb565_convert_internal ycc_rgb565_convert_be
  573. #define ycc_rgb565D_convert_internal ycc_rgb565D_convert_be
  574. #define rgb_rgb565_convert_internal rgb_rgb565_convert_be
  575. #define rgb_rgb565D_convert_internal rgb_rgb565D_convert_be
  576. #define gray_rgb565_convert_internal gray_rgb565_convert_be
  577. #define gray_rgb565D_convert_internal gray_rgb565D_convert_be
  578. #include "jdcol565.c"
  579. #undef PACK_SHORT_565
  580. #undef PACK_TWO_PIXELS
  581. #undef ycc_rgb565_convert_internal
  582. #undef ycc_rgb565D_convert_internal
  583. #undef rgb_rgb565_convert_internal
  584. #undef rgb_rgb565D_convert_internal
  585. #undef gray_rgb565_convert_internal
  586. #undef gray_rgb565D_convert_internal
  587. METHODDEF(void)
  588. ycc_rgb565_convert (j_decompress_ptr cinfo,
  589. JSAMPIMAGE input_buf, JDIMENSION input_row,
  590. JSAMPARRAY output_buf, int num_rows)
  591. {
  592. if (is_big_endian())
  593. ycc_rgb565_convert_be(cinfo, input_buf, input_row, output_buf, num_rows);
  594. else
  595. ycc_rgb565_convert_le(cinfo, input_buf, input_row, output_buf, num_rows);
  596. }
  597. METHODDEF(void)
  598. ycc_rgb565D_convert (j_decompress_ptr cinfo,
  599. JSAMPIMAGE input_buf, JDIMENSION input_row,
  600. JSAMPARRAY output_buf, int num_rows)
  601. {
  602. if (is_big_endian())
  603. ycc_rgb565D_convert_be(cinfo, input_buf, input_row, output_buf, num_rows);
  604. else
  605. ycc_rgb565D_convert_le(cinfo, input_buf, input_row, output_buf, num_rows);
  606. }
  607. METHODDEF(void)
  608. rgb_rgb565_convert (j_decompress_ptr cinfo,
  609. JSAMPIMAGE input_buf, JDIMENSION input_row,
  610. JSAMPARRAY output_buf, int num_rows)
  611. {
  612. if (is_big_endian())
  613. rgb_rgb565_convert_be(cinfo, input_buf, input_row, output_buf, num_rows);
  614. else
  615. rgb_rgb565_convert_le(cinfo, input_buf, input_row, output_buf, num_rows);
  616. }
  617. METHODDEF(void)
  618. rgb_rgb565D_convert (j_decompress_ptr cinfo,
  619. JSAMPIMAGE input_buf, JDIMENSION input_row,
  620. JSAMPARRAY output_buf, int num_rows)
  621. {
  622. if (is_big_endian())
  623. rgb_rgb565D_convert_be(cinfo, input_buf, input_row, output_buf, num_rows);
  624. else
  625. rgb_rgb565D_convert_le(cinfo, input_buf, input_row, output_buf, num_rows);
  626. }
  627. METHODDEF(void)
  628. gray_rgb565_convert (j_decompress_ptr cinfo,
  629. JSAMPIMAGE input_buf, JDIMENSION input_row,
  630. JSAMPARRAY output_buf, int num_rows)
  631. {
  632. if (is_big_endian())
  633. gray_rgb565_convert_be(cinfo, input_buf, input_row, output_buf, num_rows);
  634. else
  635. gray_rgb565_convert_le(cinfo, input_buf, input_row, output_buf, num_rows);
  636. }
  637. METHODDEF(void)
  638. gray_rgb565D_convert (j_decompress_ptr cinfo,
  639. JSAMPIMAGE input_buf, JDIMENSION input_row,
  640. JSAMPARRAY output_buf, int num_rows)
  641. {
  642. if (is_big_endian())
  643. gray_rgb565D_convert_be(cinfo, input_buf, input_row, output_buf, num_rows);
  644. else
  645. gray_rgb565D_convert_le(cinfo, input_buf, input_row, output_buf, num_rows);
  646. }
  647. /*
  648. * Empty method for start_pass.
  649. */
  650. METHODDEF(void)
  651. start_pass_dcolor (j_decompress_ptr cinfo)
  652. {
  653. /* no work needed */
  654. }
  655. /*
  656. * Module initialization routine for output colorspace conversion.
  657. */
  658. GLOBAL(void)
  659. jinit_color_deconverter (j_decompress_ptr cinfo)
  660. {
  661. my_cconvert_ptr cconvert;
  662. int ci;
  663. cconvert = (my_cconvert_ptr)
  664. (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
  665. sizeof(my_color_deconverter));
  666. cinfo->cconvert = (struct jpeg_color_deconverter *) cconvert;
  667. cconvert->pub.start_pass = start_pass_dcolor;
  668. /* Make sure num_components agrees with jpeg_color_space */
  669. switch (cinfo->jpeg_color_space) {
  670. case JCS_GRAYSCALE:
  671. if (cinfo->num_components != 1)
  672. ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
  673. break;
  674. case JCS_RGB:
  675. case JCS_YCbCr:
  676. if (cinfo->num_components != 3)
  677. ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
  678. break;
  679. case JCS_CMYK:
  680. case JCS_YCCK:
  681. if (cinfo->num_components != 4)
  682. ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
  683. break;
  684. default: /* JCS_UNKNOWN can be anything */
  685. if (cinfo->num_components < 1)
  686. ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
  687. break;
  688. }
  689. /* Set out_color_components and conversion method based on requested space.
  690. * Also clear the component_needed flags for any unused components,
  691. * so that earlier pipeline stages can avoid useless computation.
  692. */
  693. switch (cinfo->out_color_space) {
  694. case JCS_GRAYSCALE:
  695. cinfo->out_color_components = 1;
  696. if (cinfo->jpeg_color_space == JCS_GRAYSCALE ||
  697. cinfo->jpeg_color_space == JCS_YCbCr) {
  698. cconvert->pub.color_convert = grayscale_convert;
  699. /* For color->grayscale conversion, only the Y (0) component is needed */
  700. for (ci = 1; ci < cinfo->num_components; ci++)
  701. cinfo->comp_info[ci].component_needed = FALSE;
  702. } else if (cinfo->jpeg_color_space == JCS_RGB) {
  703. cconvert->pub.color_convert = rgb_gray_convert;
  704. build_rgb_y_table(cinfo);
  705. } else
  706. ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
  707. break;
  708. case JCS_RGB:
  709. case JCS_EXT_RGB:
  710. case JCS_EXT_RGBX:
  711. case JCS_EXT_BGR:
  712. case JCS_EXT_BGRX:
  713. case JCS_EXT_XBGR:
  714. case JCS_EXT_XRGB:
  715. case JCS_EXT_RGBA:
  716. case JCS_EXT_BGRA:
  717. case JCS_EXT_ABGR:
  718. case JCS_EXT_ARGB:
  719. cinfo->out_color_components = rgb_pixelsize[cinfo->out_color_space];
  720. if (cinfo->jpeg_color_space == JCS_YCbCr) {
  721. if (jsimd_can_ycc_rgb())
  722. cconvert->pub.color_convert = jsimd_ycc_rgb_convert;
  723. else {
  724. cconvert->pub.color_convert = ycc_rgb_convert;
  725. build_ycc_rgb_table(cinfo);
  726. }
  727. } else if (cinfo->jpeg_color_space == JCS_GRAYSCALE) {
  728. cconvert->pub.color_convert = gray_rgb_convert;
  729. } else if (cinfo->jpeg_color_space == JCS_RGB) {
  730. if (rgb_red[cinfo->out_color_space] == 0 &&
  731. rgb_green[cinfo->out_color_space] == 1 &&
  732. rgb_blue[cinfo->out_color_space] == 2 &&
  733. rgb_pixelsize[cinfo->out_color_space] == 3)
  734. cconvert->pub.color_convert = null_convert;
  735. else
  736. cconvert->pub.color_convert = rgb_rgb_convert;
  737. } else
  738. ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
  739. break;
  740. case JCS_RGB565:
  741. cinfo->out_color_components = 3;
  742. if (cinfo->dither_mode == JDITHER_NONE) {
  743. if (cinfo->jpeg_color_space == JCS_YCbCr) {
  744. if (jsimd_can_ycc_rgb565())
  745. cconvert->pub.color_convert = jsimd_ycc_rgb565_convert;
  746. else {
  747. cconvert->pub.color_convert = ycc_rgb565_convert;
  748. build_ycc_rgb_table(cinfo);
  749. }
  750. } else if (cinfo->jpeg_color_space == JCS_GRAYSCALE) {
  751. cconvert->pub.color_convert = gray_rgb565_convert;
  752. } else if (cinfo->jpeg_color_space == JCS_RGB) {
  753. cconvert->pub.color_convert = rgb_rgb565_convert;
  754. } else
  755. ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
  756. } else {
  757. /* only ordered dithering is supported */
  758. if (cinfo->jpeg_color_space == JCS_YCbCr) {
  759. cconvert->pub.color_convert = ycc_rgb565D_convert;
  760. build_ycc_rgb_table(cinfo);
  761. } else if (cinfo->jpeg_color_space == JCS_GRAYSCALE) {
  762. cconvert->pub.color_convert = gray_rgb565D_convert;
  763. } else if (cinfo->jpeg_color_space == JCS_RGB) {
  764. cconvert->pub.color_convert = rgb_rgb565D_convert;
  765. } else
  766. ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
  767. }
  768. break;
  769. case JCS_CMYK:
  770. cinfo->out_color_components = 4;
  771. if (cinfo->jpeg_color_space == JCS_YCCK) {
  772. cconvert->pub.color_convert = ycck_cmyk_convert;
  773. build_ycc_rgb_table(cinfo);
  774. } else if (cinfo->jpeg_color_space == JCS_CMYK) {
  775. cconvert->pub.color_convert = null_convert;
  776. } else
  777. ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
  778. break;
  779. default:
  780. /* Permit null conversion to same output space */
  781. if (cinfo->out_color_space == cinfo->jpeg_color_space) {
  782. cinfo->out_color_components = cinfo->num_components;
  783. cconvert->pub.color_convert = null_convert;
  784. } else /* unsupported non-null conversion */
  785. ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
  786. break;
  787. }
  788. if (cinfo->quantize_colors)
  789. cinfo->output_components = 1; /* single colormapped output component */
  790. else
  791. cinfo->output_components = cinfo->out_color_components;
  792. }