src/gdcmJpeg12.cxx

   1 // gdcmJpeg12.cxx
   2 //-----------------------------------------------------------------------------
   3 #include <stdio.h>
   4 #include "gdcmFile.h"
   5
   6 #define BITS_IN_JSAMPLE 12
   7
   8 #ifdef GDCM_DEBUG
   9 #define GDCM_jpr_DEBUG 0
  10 #endif   //GDCM_DEBUG
  11
  12 // BITS_IN_JSAMPLE is a compile time defined options.
  13 // We need both 8 an 12;
  14 // To avoid renaming *all* the Jpeg functions,
  15 // we hard code the 'brain damaged liker' option.
  16 // For all the functions, we shall have the 8 and 12 version
  17 // (8 with the 'long' name, 12 with the 'short' name)
  18
  19 #define jpeg_read_header        jReadHeader
  20 #define my_error_exit           myErrorExit
  21 #define jpeg_destroy_decompress jDestDecompress
  22 #define jpeg_stdio_src          jStdSrc
  23 #define jpeg_read_header        jReadHeader
  24 #define jpeg_read_scanlines     jReadScanlines
  25 #define jpeg_finish_decompress  jFinDecompress
  26 //#define jpeg_create_decompress  jCreaDecompress //FIXME
  27
  28 // -----------------
  29 #define jpeg_std_error           jStdError
  30 #define jpeg_CreateCompress      jCreaCompress
  31 #define jpeg_CreateDecompress    jCreaDecompress
  32 #define jpeg_destroy_compress    jDestCompress
  33 #define jpeg_destroy_decompress  jDestDecompress
  34 #define jpeg_stdio_dest          jStdDest
  35 #define jpeg_stdio_src           jStdSrc
  36 #define jpeg_set_defaults        jSetDefaults
  37 #define jpeg_set_colorspace      jSetColorspace
  38 #define jpeg_default_colorspace  jDefColorspace
  39 #define jpeg_set_quality         jSetQuality
  40 #define jpeg_set_linear_quality  jSetLQuality
  41 #define jpeg_add_quant_table     jAddQuantTable
  42 #define jpeg_quality_scaling     jQualityScaling
  43 #define jpeg_simple_progression  jSimProgress
  44 #define jpeg_suppress_tables     jSuppressTables
  45 #define jpeg_alloc_quant_table   jAlcQTable
  46 #define jpeg_alloc_huff_table    jAlcHTable
  47 #define jpeg_start_compress      jStrtCompress
  48 #define jpeg_write_scanlines     jWrtScanlines
  49 #define jpeg_finish_compress     jFinCompress
  50 #define jpeg_write_raw_data      jWrtRawData
  51 #define jpeg_write_marker        jWrtMarker
  52 #define jpeg_write_m_header      jWrtMHeader
  53 #define jpeg_write_m_byte        jWrtMByte
  54 #define jpeg_write_tables        jWrtTables
  55 #define jpeg_read_header         jReadHeader
  56 #define jpeg_start_decompress    jStrtDecompress
  57 #define jpeg_read_scanlines      jReadScanlines
  58 #define jpeg_finish_decompress   jFinDecompress
  59 #define jpeg_read_raw_data       jReadRawData
  60 #define jpeg_has_multiple_scans  jHasMultScn
  61 #define jpeg_start_output        jStrtOutput
  62 #define jpeg_finish_output       jFinOutput
  63 #define jpeg_input_complete      jInComplete
  64 #define jpeg_new_colormap        jNewCMap
  65 #define jpeg_consume_input          jConsumeInput
  66 #define jpeg_calc_output_dimensions jCalcDimensions
  67 #define jpeg_save_markers           jSaveMarkers
  68 #define jpeg_set_marker_processor   jSetMarker
  69 #define jpeg_read_coefficients      jReadCoefs
  70 #define jpeg_write_coefficients     jWrtCoefs
  71 #define jpeg_copy_critical_parameters jCopyCrit
  72 #define jpeg_abort_compress      jAbrtCompress
  73 #define jpeg_abort_decompress    jAbrtDecompress
  74 #define jpeg_abort               jAbort
  75 #define jpeg_destroy             jDestroy
  76 #define jpeg_resync_to_restart   jResyncRestart
  77
  78 /*
  79  * <setjmp.h> is used for the optional error recovery mechanism shown in
  80  * the second part of the example.
  81  */
  82
  83 /*
  84  * Include file for users of JPEG library.
  85  * You will need to have included system headers that define at least
  86  * the typedefs FILE and size_t before you can include jpeglib.h.
  87  * (stdio.h is sufficient on ANSI-conforming systems.)
  88  * You may also wish to include "jerror.h".
  89  */
  90
  91 extern "C" {
  92 #include "jpeg/libijg12/jpeglib12.h"
  93 #include <setjmp.h>
  94 }
  95
  96 /******************** JPEG DECOMPRESSION SAMPLE INTERFACE *******************/
  97
  98 /* This half of the example shows how to read data from the JPEG decompressor.
  99  * It's a bit more refined than the above, in that we show:
 100  *   (a) how to modify the JPEG library's standard error-reporting behavior;
 101  *   (b) how to allocate workspace using the library's memory manager.
 102  *
 103  * Just to make this example a little different from the first one, we'll
 104  * assume that we do not intend to put the whole image into an in-memory
 105  * buffer, but to send it line-by-line someplace else.  We need a one-
 106  * scanline-high JSAMPLE array as a work buffer, and we will let the JPEG
 107  * memory manager allocate it for us.  This approach is actually quite useful
 108  * because we don't need to remember to deallocate the buffer separately: it
 109  * will go away automatically when the JPEG object is cleaned up.
 110  */
 111
 112 /*
 113  * ERROR HANDLING:
 114  *
 115  * The JPEG library's standard error handler (jerror.c) is divided into
 116  * several "methods" which you can override individually.  This lets you
 117  * adjust the behavior without duplicating a lot of code, which you might
 118  * have to update with each future release.
 119  *
 120  * Our example here shows how to override the "error_exit" method so that
 121  * control is returned to the library's caller when a fatal error occurs,
 122  * rather than calling exit() as the standard error_exit method does.
 123  *
 124  * We use C's setjmp/longjmp facility to return control.  This means that the
 125  * routine which calls the JPEG library must first execute a setjmp() call to
 126  * establish the return point.  We want the replacement error_exit to do a
 127  * longjmp().  But we need to make the setjmp buffer accessible to the
 128  * error_exit routine.  To do this, we make a private extension of the
 129  * standard JPEG error handler object.  (If we were using C++, we'd say we
 130  * were making a subclass of the regular error handler.)
 131  *
 132  * Here's the extended error handler struct:
 133  */
 134
 135 //-----------------------------------------------------------------------------
 136 struct my_error_mgr {
 137    struct jpeg_error_mgr pub;   /* "public" fields */
 138    jmp_buf setjmp_buffer;       /* for return to caller */
 139 };
 140
 141 //-----------------------------------------------------------------------------
 142 typedef struct my_error_mgr * my_error_ptr;
 143
 144 /*
 145  * Here's the routine that will replace the standard error_exit method:
 146  */
 147 METHODDEF(void) my_error_exit (j_common_ptr cinfo) {
 148    /* cinfo->err really points to a my_error_mgr struct, so coerce pointer */
 149    my_error_ptr myerr = (my_error_ptr) cinfo->err;
 150
 151    /* Always display the message. */
 152    /* We could postpone this until after returning, if we chose. */
 153    (*cinfo->err->output_message) (cinfo);
 154
 155    /* Return control to the setjmp point */
 156    longjmp(myerr->setjmp_buffer, 1);
 157 }
 158
 159
 160 //-----------------------------------------------------------------------------
 161 /*
 162  * Sample routine for JPEG decompression.  We assume that the source file name
 163  * is passed in.  We want to return 1 on success, 0 on error.
 164  */
 165
 166  /**
 167  * \ingroup gdcmFile
 168  * \brief   routine for JPEG decompression
 169  * @param fp pointer to an already open file descriptor
 170  *                      12 significant bits per pixel
 171  * @param image_buffer to receive uncompressed pixels
 172  * @return 1 on success, 0 on error
 173  */
 174
 175 bool gdcmFile::gdcm_read_JPEG_file12 (FILE *fp,void * image_buffer) {
 176    char *pimage;
 177
 178    /* This struct contains the JPEG decompression parameters and pointers to
 179     * working space (which is allocated as needed by the JPEG library).
 180     */
 181    struct jpeg_decompress_struct cinfo;
 182
 183    /* -------------- inside, we found :
 184     * JDIMENSION image_width;      // input image width
 185     * JDIMENSION image_height;     // input image height
 186     * int input_components;        // nb of color components in input image
 187     * J_COLOR_SPACE in_color_space;// colorspace of input image
 188     * double input_gamma;          // image gamma of input image
 189     * -------------- */
 190
 191    /* We use our private extension JPEG error handler.
 192     * Note that this struct must live as long as the main JPEG parameter
 193     * struct, to avoid dangling-pointer problems.
 194     */
 195    struct my_error_mgr jerr;
 196    /* More stuff */
 197
 198    JSAMPARRAY buffer;/* Output row buffer */
 199
 200    // rappel :
 201    // ------
 202    // typedef unsigned char JSAMPLE;
 203    // typedef JSAMPLE FAR *JSAMPROW;/* ptr to one image row of pixel samples. */
 204    // typedef JSAMPROW *JSAMPARRAY;/* ptr to some rows (a 2-D sample array) */
 205    // typedef JSAMPARRAY *JSAMPIMAGE;/* a 3-D sample array: top index is color */
 206
 207    int row_stride;/* physical row width in output buffer */
 208 #ifdef GDCM_JPG_DEBUG
 209    printf("entree dans gdcmFile::gdcm_read_JPEG_file12, depuis gdcmJpeg\n");
 210 #endif //GDCM_JPG_DEBUG
 211
 212    /* In this example we want to open the input file before doing anything else,
 213     * so that the setjmp() error recovery below can assume the file is open.
 214     * VERY IMPORTANT: use "b" option to fopen() if you are on a machine that
 215     * requires it in order to read binary files.
 216     */
 217
 218    /* Step 1: allocate and initialize JPEG decompression object */
 219 #ifdef GDCM_JPG_DEBUG
 220    printf("Entree Step 1\n");
 221 #endif //GDCM_JPG_DEBUG
 222
 223    /* We set up the normal JPEG error routines, then override error_exit. */
 224
 225    cinfo.err = jpeg_std_error(&jerr.pub);
 226    jerr.pub.error_exit = my_error_exit;
 227
 228    /* Establish the setjmp return context for my_error_exit to use. */
 229    if (setjmp(jerr.setjmp_buffer)) {
 230       /* If we get here, the JPEG code has signaled an error.
 231        * We need to clean up the JPEG object, close the input file, and return.
 232        */
 233       jpeg_destroy_decompress(&cinfo);
 234       return(false);
 235    }
 236
 237    /* Now we can initialize the JPEG decompression object. */
 238    jpeg_create_decompress(&cinfo);
 239
 240    /* Step 2: specify data source (eg, a file) */
 241 #ifdef GDCM_JPG_DEBUG
 242    printf("Entree Step 2\n");
 243 #endif //GDCM_JPG_DEBUG
 244    jpeg_stdio_src(&cinfo, fp);
 245
 246    /* Step 3: read file parameters with jpeg_read_header() */
 247 #ifdef GDCM_JPG_DEBUG
 248    printf("Entree Step 3\n");
 249 #endif //GDCM_JPG_DEBUG
 250    (void) jpeg_read_header(&cinfo, TRUE);
 251
 252    /* We can ignore the return value from jpeg_read_header since
 253     *   (a) suspension is not possible with the stdio data source, and
 254     *   (b) we passed TRUE to reject a tables-only JPEG file as an error.
 255     * See libjpeg.doc for more info.
 256     */
 257
 258 #ifdef GDCM_JPG_DEBUG
 259       printf("--------------Header contents :----------------\n");
 260       printf("image_width %d image_height %d\n",
 261               cinfo.image_width , cinfo.image_height);
 262       printf("bits of precision in image data  %d \n",
 263               cinfo.output_components);
 264       printf("nb of color components returned  %d \n",
 265               cinfo.data_precision);
 266 #endif //GDCM_JPG_DEBUG
 267
 268
 269    /*
 270     * JDIMENSION image_width;       // input image width
 271     * JDIMENSION image_height;      // input image height
 272     * int output_components;        // # of color components returned
 273     * J_COLOR_SPACE in_color_space; // colorspace of input image
 274     * double input_gamma;           // image gamma of input image
 275     * int data_precision;           // bits of precision in image data
 276     */
 277
 278    /* Step 4: set parameters for decompression */
 279 #ifdef GDCM_JPG_DEBUG
 280    printf("Entree Step 4\n");
 281 #endif //GDCM_JPG_DEBUG
 282
 283    /* In this example, we don't need to change any of the defaults set by
 284     * jpeg_read_header(), so we do nothing here.
 285     */
 286
 287    /* Step 5: Start decompressor */
 288 #ifdef GDCM_JPG_DEBUG
 289    printf("Entree Step 5\n");
 290 #endif //GDCM_JPG_DEBUG
 291
 292    (void) jpeg_start_decompress(&cinfo);
 293    /* We can ignore the return value since suspension is not possible
 294     * with the stdio data source.
 295     */
 296
 297    /* We may need to do some setup of our own at this point before reading
 298     * the data.  After jpeg_start_decompress() we have the correct scaled
 299     * output image dimensions available, as well as the output colormap
 300     * if we asked for color quantization.
 301     * In this example, we need to make an output work buffer of the right size.
 302     */
 303
 304    /* JSAMPLEs per row in output buffer */
 305    row_stride = cinfo.output_width * cinfo.output_components;
 306
 307 #ifdef GDCM_JPG_DEBUG
 308       printf ("cinfo.output_width %d cinfo.output_components %d  row_stride %d\n",
 309               cinfo.output_width, cinfo.output_components,row_stride);
 310 #endif //GDCM_JPG_DEBUG
 311
 312    /* Make a one-row-high sample array that will go away when done with image */
 313    buffer = (*cinfo.mem->alloc_sarray)
 314     ((j_common_ptr) &cinfo, JPOOL_IMAGE, row_stride, 1);
 315
 316    /* Step 6: while (scan lines remain to be read) */
 317 #ifdef GDCM_JPG_DEBUG
 318    printf("Entree Step 6\n");
 319 #endif //GDCM_JPG_DEBUG
 320
 321    /*           jpeg_read_scanlines(...); */
 322
 323    /* Here we use the library's state variable cinfo.output_scanline as the
 324     * loop counter, so that we don't have to keep track ourselves.
 325     */
 326
 327 #ifdef GDCM_JPG_DEBUG
 328    printf ("cinfo.output_height %d  cinfo.output_width %d\n",
 329       cinfo.output_height,cinfo.output_width);
 330 #endif //GDCM_JPG_DEBUG
 331
 332    pimage=(char *)image_buffer;
 333
 334    while (cinfo.output_scanline < cinfo.output_height) {
 335       /* jpeg_read_scanlines expects an array of pointers to scanlines.
 336        * Here the array is only one element long, but you could ask for
 337        * more than one scanline at a time if that's more convenient.
 338        */
 339
 340       (void) jpeg_read_scanlines(&cinfo, buffer, 1);
 341
 342       if ( BITS_IN_JSAMPLE == 8) {
 343          memcpy( pimage, buffer[0],row_stride);
 344          pimage+=row_stride;
 345       } else {
 346          memcpy( pimage, buffer[0],row_stride*2 ); // FIXME : *2  car 16 bits?!?
 347          pimage+=row_stride*2;                     // FIXME : *2 car 16 bits?!?
 348       }
 349    }
 350
 351   /* Step 7: Finish decompression */
 352 #ifdef GDCM_JPG_DEBUG
 353    printf("Entree Step 7\n");
 354 #endif //GDCM_JPG_DEBUG
 355    (void) jpeg_finish_decompress(&cinfo);
 356    /* We can ignore the return value since suspension is not possible
 357     * with the stdio data source.
 358     */
 359
 360    /* Step 8: Release JPEG decompression object */
 361 #ifdef GDCM_JPG_DEBUG
 362    printf("Entree Step 8\n");
 363 #endif //GDCM_JPG_DEBUG
 364
 365    /* This is an important step since it will release a good deal of memory. */
 366    jpeg_destroy_decompress(&cinfo);
 367
 368    /* After finish_decompress, we can close the input file.
 369     * Here we postpone it until after no more JPEG errors are possible,
 370     * so as to simplify the setjmp error logic above.  (Actually, I don't
 371     * think that jpeg_destroy can do an error exit, but why assume anything...)
 372     */
 373
 374    /* At this point you may want to check to see whether any corrupt-data
 375     * warnings occurred (test whether jerr.pub.num_warnings is nonzero).
 376     */
 377
 378    /* And we're done! */
 379
 380    return(true);
 381 }
 382
 383 /*
 384  * SOME FINE POINTS:
 385  *
 386  * In the above code, we ignored the return value of jpeg_read_scanlines,
 387  * which is the number of scanlines actually read.  We could get away with
 388  * this because we asked for only one line at a time and we weren't using
 389  * a suspending data source.  See libjpeg.doc for more info.
 390  *
 391  * We cheated a bit by calling alloc_sarray() after jpeg_start_decompress();
 392  * we should have done it beforehand to ensure that the space would be
 393  * counted against the JPEG max_memory setting.  In some systems the above
 394  * code would risk an out-of-memory error.  However, in general we don't
 395  * know the output image dimensions before jpeg_start_decompress(), unless we
 396  * call jpeg_calc_output_dimensions().  See libjpeg.doc for more about this.
 397  *
 398  * Scanlines are returned in the same order as they appear in the JPEG file,
 399  * which is standardly top-to-bottom.  If you must emit data bottom-to-top,
 400  * you can use one of the virtual arrays provided by the JPEG memory manager
 401  * to invert the data.  See wrbmp.c for an example.
 402  *
 403  * As with compression, some operating modes may require temporary files.
 404  * On some systems you may need to set up a signal handler to ensure that
 405  * temporary files are deleted if the program is interrupted.  See libjpeg.doc.
 406  */
 407
 408 //-----------------------------------------------------------------------------