Home | History | Annotate | Download | only in libjpeg
      1 /*
      2  * jdcolor.c
      3  *
      4  * Copyright (C) 1991-1997, Thomas G. Lane.
      5  * Modified 2011-2012 by Guido Vollbeding.
      6  * This file is part of the Independent JPEG Group's software.
      7  * For conditions of distribution and use, see the accompanying README file.
      8  *
      9  * This file contains output colorspace conversion routines.
     10  */
     11 
     12 #define JPEG_INTERNALS
     13 #include "jinclude.h"
     14 #include "jpeglib.h"
     15 
     16 
     17 /* Private subobject */
     18 
     19 typedef struct {
     20   struct jpeg_color_deconverter pub; /* public fields */
     21 
     22   /* Private state for YCC->RGB conversion */
     23   int * Cr_r_tab;		/* => table for Cr to R conversion */
     24   int * Cb_b_tab;		/* => table for Cb to B conversion */
     25   INT32 * Cr_g_tab;		/* => table for Cr to G conversion */
     26   INT32 * Cb_g_tab;		/* => table for Cb to G conversion */
     27 
     28   /* Private state for RGB->Y conversion */
     29   INT32 * rgb_y_tab;		/* => table for RGB to Y conversion */
     30 } my_color_deconverter;
     31 
     32 typedef my_color_deconverter * my_cconvert_ptr;
     33 
     34 
     35 /**************** YCbCr -> RGB conversion: most common case **************/
     36 /****************   RGB -> Y   conversion: less common case **************/
     37 
     38 /*
     39  * YCbCr is defined per CCIR 601-1, except that Cb and Cr are
     40  * normalized to the range 0..MAXJSAMPLE rather than -0.5 .. 0.5.
     41  * The conversion equations to be implemented are therefore
     42  *
     43  *	R = Y                + 1.40200 * Cr
     44  *	G = Y - 0.34414 * Cb - 0.71414 * Cr
     45  *	B = Y + 1.77200 * Cb
     46  *
     47  *	Y = 0.29900 * R + 0.58700 * G + 0.11400 * B
     48  *
     49  * where Cb and Cr represent the incoming values less CENTERJSAMPLE.
     50  * (These numbers are derived from TIFF 6.0 section 21, dated 3-June-92.)
     51  *
     52  * To avoid floating-point arithmetic, we represent the fractional constants
     53  * as integers scaled up by 2^16 (about 4 digits precision); we have to divide
     54  * the products by 2^16, with appropriate rounding, to get the correct answer.
     55  * Notice that Y, being an integral input, does not contribute any fraction
     56  * so it need not participate in the rounding.
     57  *
     58  * For even more speed, we avoid doing any multiplications in the inner loop
     59  * by precalculating the constants times Cb and Cr for all possible values.
     60  * For 8-bit JSAMPLEs this is very reasonable (only 256 entries per table);
     61  * for 12-bit samples it is still acceptable.  It's not very reasonable for
     62  * 16-bit samples, but if you want lossless storage you shouldn't be changing
     63  * colorspace anyway.
     64  * The Cr=>R and Cb=>B values can be rounded to integers in advance; the
     65  * values for the G calculation are left scaled up, since we must add them
     66  * together before rounding.
     67  */
     68 
     69 #define SCALEBITS	16	/* speediest right-shift on some machines */
     70 #define ONE_HALF	((INT32) 1 << (SCALEBITS-1))
     71 #define FIX(x)		((INT32) ((x) * (1L<<SCALEBITS) + 0.5))
     72 
     73 /* We allocate one big table for RGB->Y conversion and divide it up into
     74  * three parts, instead of doing three alloc_small requests.  This lets us
     75  * use a single table base address, which can be held in a register in the
     76  * inner loops on many machines (more than can hold all three addresses,
     77  * anyway).
     78  */
     79 
     80 #define R_Y_OFF		0			/* offset to R => Y section */
     81 #define G_Y_OFF		(1*(MAXJSAMPLE+1))	/* offset to G => Y section */
     82 #define B_Y_OFF		(2*(MAXJSAMPLE+1))	/* etc. */
     83 #define TABLE_SIZE	(3*(MAXJSAMPLE+1))
     84 
     85 
     86 /*
     87  * Initialize tables for YCC->RGB colorspace conversion.
     88  */
     89 
     90 LOCAL(void)
     91 build_ycc_rgb_table (j_decompress_ptr cinfo)
     92 {
     93   my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
     94   int i;
     95   INT32 x;
     96   SHIFT_TEMPS
     97 
     98   cconvert->Cr_r_tab = (int *)
     99     (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
    100                                 (MAXJSAMPLE+1) * SIZEOF(int));
    101   cconvert->Cb_b_tab = (int *)
    102     (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
    103                                 (MAXJSAMPLE+1) * SIZEOF(int));
    104   cconvert->Cr_g_tab = (INT32 *)
    105     (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
    106                                 (MAXJSAMPLE+1) * SIZEOF(INT32));
    107   cconvert->Cb_g_tab = (INT32 *)
    108     (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
    109                                 (MAXJSAMPLE+1) * SIZEOF(INT32));
    110 
    111   for (i = 0, x = -CENTERJSAMPLE; i <= MAXJSAMPLE; i++, x++) {
    112     /* i is the actual input pixel value, in the range 0..MAXJSAMPLE */
    113     /* The Cb or Cr value we are thinking of is x = i - CENTERJSAMPLE */
    114     /* Cr=>R value is nearest int to 1.40200 * x */
    115     cconvert->Cr_r_tab[i] = (int)
    116                     RIGHT_SHIFT(FIX(1.40200) * x + ONE_HALF, SCALEBITS);
    117     /* Cb=>B value is nearest int to 1.77200 * x */
    118     cconvert->Cb_b_tab[i] = (int)
    119                     RIGHT_SHIFT(FIX(1.77200) * x + ONE_HALF, SCALEBITS);
    120     /* Cr=>G value is scaled-up -0.71414 * x */
    121     cconvert->Cr_g_tab[i] = (- FIX(0.71414)) * x;
    122     /* Cb=>G value is scaled-up -0.34414 * x */
    123     /* We also add in ONE_HALF so that need not do it in inner loop */
    124     cconvert->Cb_g_tab[i] = (- FIX(0.34414)) * x + ONE_HALF;
    125   }
    126 }
    127 
    128 
    129 /*
    130  * Convert some rows of samples to the output colorspace.
    131  *
    132  * Note that we change from noninterleaved, one-plane-per-component format
    133  * to interleaved-pixel format.  The output buffer is therefore three times
    134  * as wide as the input buffer.
    135  * A starting row offset is provided only for the input buffer.  The caller
    136  * can easily adjust the passed output_buf value to accommodate any row
    137  * offset required on that side.
    138  */
    139 
    140 METHODDEF(void)
    141 ycc_rgb_convert (j_decompress_ptr cinfo,
    142                  JSAMPIMAGE input_buf, JDIMENSION input_row,
    143                  JSAMPARRAY output_buf, int num_rows)
    144 {
    145   my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
    146   register int y, cb, cr;
    147   register JSAMPROW outptr;
    148   register JSAMPROW inptr0, inptr1, inptr2;
    149   register JDIMENSION col;
    150   JDIMENSION num_cols = cinfo->output_width;
    151   /* copy these pointers into registers if possible */
    152   register JSAMPLE * range_limit = cinfo->sample_range_limit;
    153   register int * Crrtab = cconvert->Cr_r_tab;
    154   register int * Cbbtab = cconvert->Cb_b_tab;
    155   register INT32 * Crgtab = cconvert->Cr_g_tab;
    156   register INT32 * Cbgtab = cconvert->Cb_g_tab;
    157   SHIFT_TEMPS
    158 
    159   while (--num_rows >= 0) {
    160     inptr0 = input_buf[0][input_row];
    161     inptr1 = input_buf[1][input_row];
    162     inptr2 = input_buf[2][input_row];
    163     input_row++;
    164     outptr = *output_buf++;
    165     for (col = 0; col < num_cols; col++) {
    166       y  = GETJSAMPLE(inptr0[col]);
    167       cb = GETJSAMPLE(inptr1[col]);
    168       cr = GETJSAMPLE(inptr2[col]);
    169       /* Range-limiting is essential due to noise introduced by DCT losses. */
    170       outptr[RGB_RED] =   range_limit[y + Crrtab[cr]];
    171       outptr[RGB_GREEN] = range_limit[y +
    172                               ((int) RIGHT_SHIFT(Cbgtab[cb] + Crgtab[cr],
    173                                                  SCALEBITS))];
    174       outptr[RGB_BLUE] =  range_limit[y + Cbbtab[cb]];
    175       outptr += RGB_PIXELSIZE;
    176     }
    177   }
    178 }
    179 
    180 
    181 /**************** Cases other than YCbCr -> RGB **************/
    182 
    183 
    184 /*
    185  * Initialize for RGB->grayscale colorspace conversion.
    186  */
    187 
    188 LOCAL(void)
    189 build_rgb_y_table (j_decompress_ptr cinfo)
    190 {
    191   my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
    192   INT32 * rgb_y_tab;
    193   INT32 i;
    194 
    195   /* Allocate and fill in the conversion tables. */
    196   cconvert->rgb_y_tab = rgb_y_tab = (INT32 *)
    197     (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
    198                                 (TABLE_SIZE * SIZEOF(INT32)));
    199 
    200   for (i = 0; i <= MAXJSAMPLE; i++) {
    201     rgb_y_tab[i+R_Y_OFF] = FIX(0.29900) * i;
    202     rgb_y_tab[i+G_Y_OFF] = FIX(0.58700) * i;
    203     rgb_y_tab[i+B_Y_OFF] = FIX(0.11400) * i + ONE_HALF;
    204   }
    205 }
    206 
    207 
    208 /*
    209  * Convert RGB to grayscale.
    210  */
    211 
    212 METHODDEF(void)
    213 rgb_gray_convert (j_decompress_ptr cinfo,
    214                   JSAMPIMAGE input_buf, JDIMENSION input_row,
    215                   JSAMPARRAY output_buf, int num_rows)
    216 {
    217   my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
    218   register INT32 * ctab = cconvert->rgb_y_tab;
    219   register int r, g, b;
    220   register JSAMPROW outptr;
    221   register JSAMPROW inptr0, inptr1, inptr2;
    222   register JDIMENSION col;
    223   JDIMENSION num_cols = cinfo->output_width;
    224 
    225   while (--num_rows >= 0) {
    226     inptr0 = input_buf[0][input_row];
    227     inptr1 = input_buf[1][input_row];
    228     inptr2 = input_buf[2][input_row];
    229     input_row++;
    230     outptr = *output_buf++;
    231     for (col = 0; col < num_cols; col++) {
    232       r = GETJSAMPLE(inptr0[col]);
    233       g = GETJSAMPLE(inptr1[col]);
    234       b = GETJSAMPLE(inptr2[col]);
    235       /* Y */
    236       outptr[col] = (JSAMPLE)
    237                 ((ctab[r+R_Y_OFF] + ctab[g+G_Y_OFF] + ctab[b+B_Y_OFF])
    238                  >> SCALEBITS);
    239     }
    240   }
    241 }
    242 
    243 
    244 /*
    245  * [R-G,G,B-G] to [R,G,B] conversion with modulo calculation
    246  * (inverse color transform).
    247  */
    248 
    249 METHODDEF(void)
    250 rgb1_rgb_convert (j_decompress_ptr cinfo,
    251                   JSAMPIMAGE input_buf, JDIMENSION input_row,
    252                   JSAMPARRAY output_buf, int num_rows)
    253 {
    254   register int r, g, b;
    255   register JSAMPROW outptr;
    256   register JSAMPROW inptr0, inptr1, inptr2;
    257   register JDIMENSION col;
    258   JDIMENSION num_cols = cinfo->output_width;
    259 
    260   while (--num_rows >= 0) {
    261     inptr0 = input_buf[0][input_row];
    262     inptr1 = input_buf[1][input_row];
    263     inptr2 = input_buf[2][input_row];
    264     input_row++;
    265     outptr = *output_buf++;
    266     for (col = 0; col < num_cols; col++) {
    267       r = GETJSAMPLE(inptr0[col]);
    268       g = GETJSAMPLE(inptr1[col]);
    269       b = GETJSAMPLE(inptr2[col]);
    270       /* Assume that MAXJSAMPLE+1 is a power of 2, so that the MOD
    271        * (modulo) operator is equivalent to the bitmask operator AND.
    272        */
    273       outptr[RGB_RED]   = (JSAMPLE) ((r + g - CENTERJSAMPLE) & MAXJSAMPLE);
    274       outptr[RGB_GREEN] = (JSAMPLE) g;
    275       outptr[RGB_BLUE]  = (JSAMPLE) ((b + g - CENTERJSAMPLE) & MAXJSAMPLE);
    276       outptr += RGB_PIXELSIZE;
    277     }
    278   }
    279 }
    280 
    281 
    282 /*
    283  * [R-G,G,B-G] to grayscale conversion with modulo calculation
    284  * (inverse color transform).
    285  */
    286 
    287 METHODDEF(void)
    288 rgb1_gray_convert (j_decompress_ptr cinfo,
    289                    JSAMPIMAGE input_buf, JDIMENSION input_row,
    290                    JSAMPARRAY output_buf, int num_rows)
    291 {
    292   my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
    293   register INT32 * ctab = cconvert->rgb_y_tab;
    294   register int r, g, b;
    295   register JSAMPROW outptr;
    296   register JSAMPROW inptr0, inptr1, inptr2;
    297   register JDIMENSION col;
    298   JDIMENSION num_cols = cinfo->output_width;
    299 
    300   while (--num_rows >= 0) {
    301     inptr0 = input_buf[0][input_row];
    302     inptr1 = input_buf[1][input_row];
    303     inptr2 = input_buf[2][input_row];
    304     input_row++;
    305     outptr = *output_buf++;
    306     for (col = 0; col < num_cols; col++) {
    307       r = GETJSAMPLE(inptr0[col]);
    308       g = GETJSAMPLE(inptr1[col]);
    309       b = GETJSAMPLE(inptr2[col]);
    310       /* Assume that MAXJSAMPLE+1 is a power of 2, so that the MOD
    311        * (modulo) operator is equivalent to the bitmask operator AND.
    312        */
    313       r = (r + g - CENTERJSAMPLE) & MAXJSAMPLE;
    314       b = (b + g - CENTERJSAMPLE) & MAXJSAMPLE;
    315       /* Y */
    316       outptr[col] = (JSAMPLE)
    317                 ((ctab[r+R_Y_OFF] + ctab[g+G_Y_OFF] + ctab[b+B_Y_OFF])
    318                  >> SCALEBITS);
    319     }
    320   }
    321 }
    322 
    323 
    324 /*
    325  * No colorspace change, but conversion from separate-planes
    326  * to interleaved representation.
    327  */
    328 
    329 METHODDEF(void)
    330 rgb_convert (j_decompress_ptr cinfo,
    331              JSAMPIMAGE input_buf, JDIMENSION input_row,
    332              JSAMPARRAY output_buf, int num_rows)
    333 {
    334   register JSAMPROW outptr;
    335   register JSAMPROW inptr0, inptr1, inptr2;
    336   register JDIMENSION col;
    337   JDIMENSION num_cols = cinfo->output_width;
    338 
    339   while (--num_rows >= 0) {
    340     inptr0 = input_buf[0][input_row];
    341     inptr1 = input_buf[1][input_row];
    342     inptr2 = input_buf[2][input_row];
    343     input_row++;
    344     outptr = *output_buf++;
    345     for (col = 0; col < num_cols; col++) {
    346       /* We can dispense with GETJSAMPLE() here */
    347       outptr[RGB_RED]   = inptr0[col];
    348       outptr[RGB_GREEN] = inptr1[col];
    349       outptr[RGB_BLUE]  = inptr2[col];
    350       outptr += RGB_PIXELSIZE;
    351     }
    352   }
    353 }
    354 
    355 
    356 /*
    357  * Color conversion for no colorspace change: just copy the data,
    358  * converting from separate-planes to interleaved representation.
    359  */
    360 
    361 METHODDEF(void)
    362 null_convert (j_decompress_ptr cinfo,
    363               JSAMPIMAGE input_buf, JDIMENSION input_row,
    364               JSAMPARRAY output_buf, int num_rows)
    365 {
    366   int ci;
    367   register int nc = cinfo->num_components;
    368   register JSAMPROW outptr;
    369   register JSAMPROW inptr;
    370   register JDIMENSION col;
    371   JDIMENSION num_cols = cinfo->output_width;
    372 
    373   while (--num_rows >= 0) {
    374     for (ci = 0; ci < nc; ci++) {
    375       inptr = input_buf[ci][input_row];
    376       outptr = output_buf[0] + ci;
    377       for (col = 0; col < num_cols; col++) {
    378         *outptr = *inptr++;	/* needn't bother with GETJSAMPLE() here */
    379         outptr += nc;
    380       }
    381     }
    382     input_row++;
    383     output_buf++;
    384   }
    385 }
    386 
    387 
    388 /*
    389  * Color conversion for grayscale: just copy the data.
    390  * This also works for YCbCr -> grayscale conversion, in which
    391  * we just copy the Y (luminance) component and ignore chrominance.
    392  */
    393 
    394 METHODDEF(void)
    395 grayscale_convert (j_decompress_ptr cinfo,
    396                    JSAMPIMAGE input_buf, JDIMENSION input_row,
    397                    JSAMPARRAY output_buf, int num_rows)
    398 {
    399   jcopy_sample_rows(input_buf[0], (int) input_row, output_buf, 0,
    400                     num_rows, cinfo->output_width);
    401 }
    402 
    403 
    404 /*
    405  * Convert grayscale to RGB: just duplicate the graylevel three times.
    406  * This is provided to support applications that don't want to cope
    407  * with grayscale as a separate case.
    408  */
    409 
    410 METHODDEF(void)
    411 gray_rgb_convert (j_decompress_ptr cinfo,
    412                   JSAMPIMAGE input_buf, JDIMENSION input_row,
    413                   JSAMPARRAY output_buf, int num_rows)
    414 {
    415   register JSAMPROW outptr;
    416   register JSAMPROW inptr;
    417   register JDIMENSION col;
    418   JDIMENSION num_cols = cinfo->output_width;
    419 
    420   while (--num_rows >= 0) {
    421     inptr = input_buf[0][input_row++];
    422     outptr = *output_buf++;
    423     for (col = 0; col < num_cols; col++) {
    424       /* We can dispense with GETJSAMPLE() here */
    425       outptr[RGB_RED] = outptr[RGB_GREEN] = outptr[RGB_BLUE] = inptr[col];
    426       outptr += RGB_PIXELSIZE;
    427     }
    428   }
    429 }
    430 
    431 
    432 /*
    433  * Adobe-style YCCK->CMYK conversion.
    434  * We convert YCbCr to R=1-C, G=1-M, and B=1-Y using the same
    435  * conversion as above, while passing K (black) unchanged.
    436  * We assume build_ycc_rgb_table has been called.
    437  */
    438 
    439 METHODDEF(void)
    440 ycck_cmyk_convert (j_decompress_ptr cinfo,
    441                    JSAMPIMAGE input_buf, JDIMENSION input_row,
    442                    JSAMPARRAY output_buf, int num_rows)
    443 {
    444   my_cconvert_ptr cconvert = (my_cconvert_ptr) cinfo->cconvert;
    445   register int y, cb, cr;
    446   register JSAMPROW outptr;
    447   register JSAMPROW inptr0, inptr1, inptr2, inptr3;
    448   register JDIMENSION col;
    449   JDIMENSION num_cols = cinfo->output_width;
    450   /* copy these pointers into registers if possible */
    451   register JSAMPLE * range_limit = cinfo->sample_range_limit;
    452   register int * Crrtab = cconvert->Cr_r_tab;
    453   register int * Cbbtab = cconvert->Cb_b_tab;
    454   register INT32 * Crgtab = cconvert->Cr_g_tab;
    455   register INT32 * Cbgtab = cconvert->Cb_g_tab;
    456   SHIFT_TEMPS
    457 
    458   while (--num_rows >= 0) {
    459     inptr0 = input_buf[0][input_row];
    460     inptr1 = input_buf[1][input_row];
    461     inptr2 = input_buf[2][input_row];
    462     inptr3 = input_buf[3][input_row];
    463     input_row++;
    464     outptr = *output_buf++;
    465     for (col = 0; col < num_cols; col++) {
    466       y  = GETJSAMPLE(inptr0[col]);
    467       cb = GETJSAMPLE(inptr1[col]);
    468       cr = GETJSAMPLE(inptr2[col]);
    469       /* Range-limiting is essential due to noise introduced by DCT losses. */
    470       outptr[0] = range_limit[MAXJSAMPLE - (y + Crrtab[cr])];	/* red */
    471       outptr[1] = range_limit[MAXJSAMPLE - (y +			/* green */
    472                               ((int) RIGHT_SHIFT(Cbgtab[cb] + Crgtab[cr],
    473                                                  SCALEBITS)))];
    474       outptr[2] = range_limit[MAXJSAMPLE - (y + Cbbtab[cb])];	/* blue */
    475       /* K passes through unchanged */
    476       outptr[3] = inptr3[col];	/* don't need GETJSAMPLE here */
    477       outptr += 4;
    478     }
    479   }
    480 }
    481 
    482 
    483 /*
    484  * Empty method for start_pass.
    485  */
    486 
    487 METHODDEF(void)
    488 start_pass_dcolor (j_decompress_ptr cinfo)
    489 {
    490   /* no work needed */
    491 }
    492 
    493 
    494 /*
    495  * Module initialization routine for output colorspace conversion.
    496  */
    497 
    498 GLOBAL(void)
    499 jinit_color_deconverter (j_decompress_ptr cinfo)
    500 {
    501   my_cconvert_ptr cconvert;
    502   int ci;
    503 
    504   cconvert = (my_cconvert_ptr)
    505     (*cinfo->mem->alloc_small) ((j_common_ptr) cinfo, JPOOL_IMAGE,
    506                                 SIZEOF(my_color_deconverter));
    507   cinfo->cconvert = &cconvert->pub;
    508   cconvert->pub.start_pass = start_pass_dcolor;
    509 
    510   /* Make sure num_components agrees with jpeg_color_space */
    511   switch (cinfo->jpeg_color_space) {
    512   case JCS_GRAYSCALE:
    513     if (cinfo->num_components != 1)
    514       ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
    515     break;
    516 
    517   case JCS_RGB:
    518   case JCS_YCbCr:
    519     if (cinfo->num_components != 3)
    520       ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
    521     break;
    522 
    523   case JCS_CMYK:
    524   case JCS_YCCK:
    525     if (cinfo->num_components != 4)
    526       ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
    527     break;
    528 
    529   default:			/* JCS_UNKNOWN can be anything */
    530     if (cinfo->num_components < 1)
    531       ERREXIT(cinfo, JERR_BAD_J_COLORSPACE);
    532     break;
    533   }
    534 
    535   /* Support color transform only for RGB colorspace */
    536   if (cinfo->color_transform && cinfo->jpeg_color_space != JCS_RGB)
    537     ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
    538 
    539   /* Set out_color_components and conversion method based on requested space.
    540    * Also clear the component_needed flags for any unused components,
    541    * so that earlier pipeline stages can avoid useless computation.
    542    */
    543 
    544   switch (cinfo->out_color_space) {
    545   case JCS_GRAYSCALE:
    546     cinfo->out_color_components = 1;
    547     if (cinfo->jpeg_color_space == JCS_GRAYSCALE ||
    548         cinfo->jpeg_color_space == JCS_YCbCr) {
    549       cconvert->pub.color_convert = grayscale_convert;
    550       /* For color->grayscale conversion, only the Y (0) component is needed */
    551       for (ci = 1; ci < cinfo->num_components; ci++)
    552         cinfo->comp_info[ci].component_needed = FALSE;
    553     } else if (cinfo->jpeg_color_space == JCS_RGB) {
    554       switch (cinfo->color_transform) {
    555       case JCT_NONE:
    556         cconvert->pub.color_convert = rgb_gray_convert;
    557         break;
    558       case JCT_SUBTRACT_GREEN:
    559         cconvert->pub.color_convert = rgb1_gray_convert;
    560         break;
    561       default:
    562         ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
    563         break;
    564       }
    565       build_rgb_y_table(cinfo);
    566     } else
    567       ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
    568     break;
    569 
    570   case JCS_RGB:
    571     cinfo->out_color_components = RGB_PIXELSIZE;
    572     if (cinfo->jpeg_color_space == JCS_YCbCr) {
    573       cconvert->pub.color_convert = ycc_rgb_convert;
    574       build_ycc_rgb_table(cinfo);
    575     } else if (cinfo->jpeg_color_space == JCS_GRAYSCALE) {
    576       cconvert->pub.color_convert = gray_rgb_convert;
    577     } else if (cinfo->jpeg_color_space == JCS_RGB) {
    578       switch (cinfo->color_transform) {
    579       case JCT_NONE:
    580         cconvert->pub.color_convert = rgb_convert;
    581         break;
    582       case JCT_SUBTRACT_GREEN:
    583         cconvert->pub.color_convert = rgb1_rgb_convert;
    584         break;
    585       default:
    586         ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
    587         break;
    588       }
    589     } else
    590       ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
    591     break;
    592 
    593   case JCS_CMYK:
    594     cinfo->out_color_components = 4;
    595     if (cinfo->jpeg_color_space == JCS_YCCK) {
    596       cconvert->pub.color_convert = ycck_cmyk_convert;
    597       build_ycc_rgb_table(cinfo);
    598     } else if (cinfo->jpeg_color_space == JCS_CMYK) {
    599       cconvert->pub.color_convert = null_convert;
    600     } else
    601       ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
    602     break;
    603 
    604   default:
    605     /* Permit null conversion to same output space */
    606     if (cinfo->out_color_space == cinfo->jpeg_color_space) {
    607       cinfo->out_color_components = cinfo->num_components;
    608       cconvert->pub.color_convert = null_convert;
    609     } else			/* unsupported non-null conversion */
    610       ERREXIT(cinfo, JERR_CONVERSION_NOTIMPL);
    611     break;
    612   }
    613 
    614   if (cinfo->quantize_colors)
    615     cinfo->output_components = 1; /* single colormapped output component */
    616   else
    617     cinfo->output_components = cinfo->out_color_components;
    618 }
    619