src/gdcmJpeg12.cxx

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