src/gdcmopenjpeg/codec/j2k_to_image.c

   1 /*
   2  * Copyright (c) 2001-2003, David Janssens
   3  * Copyright (c) 2002-2003, Yannick Verschueren
   4  * Copyright (c) 2003-2005, Francois Devaux and Antonin Descampe
   5  * Copyright (c) 2005, Hervé Drolon, FreeImage Team
   6  * Copyright (c) 2002-2005, Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium
   7  * All rights reserved.
   8  *
   9  * Redistribution and use in source and binary forms, with or without
  10  * modification, are permitted provided that the following conditions
  11  * are met:
  12  * 1. Redistributions of source code must retain the above copyright
  13  *    notice, this list of conditions and the following disclaimer.
  14  * 2. Redistributions in binary form must reproduce the above copyright
  15  *    notice, this list of conditions and the following disclaimer in the
  16  *    documentation and/or other materials provided with the distribution.
  17  *
  18  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
  19  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  20  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  21  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  22  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  23  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  24  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  25  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  26  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  27  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  28  * POSSIBILITY OF SUCH DAMAGE.
  29  */
  30 #include <stdio.h>
  31 #include <string.h>
  32 #include <stdlib.h>
  33
  34 #include "openjpeg.h"
  35 #include "compat/getopt.h"
  36 #include "convert.h"
  37
  38 #ifndef WIN32
  39 #define stricmp strcasecmp
  40 #define strnicmp strncasecmp
  41 #endif
  42
  43 /* ----------------------------------------------------------------------- */
  44
  45 #define J2K_CFMT 0
  46 #define JP2_CFMT 1
  47 #define JPT_CFMT 2
  48 #define MJ2_CFMT 3
  49 #define PXM_DFMT 0
  50 #define PGX_DFMT 1
  51 #define BMP_DFMT 2
  52 #define YUV_DFMT 3
  53
  54 /* ----------------------------------------------------------------------- */
  55
  56 void decode_help_display() {
  57   fprintf(stdout,"HELP\n----\n\n");
  58   fprintf(stdout,"- the -h option displays this help information on screen\n\n");
  59
  60   fprintf(stdout,"List of parameters for the JPEG 2000 encoder:\n");
  61   fprintf(stdout,"\n");
  62   fprintf(stdout,"  -i <compressed file>\n");
  63   fprintf(stdout,"    REQUIRED\n");
  64   fprintf(stdout,"    Currently accepts J2K-files, JP2-files and JPT-files. The file type\n");
  65   fprintf(stdout,"    is identified based on its suffix.\n");
  66   fprintf(stdout,"  -o <decompressed file>\n");
  67   fprintf(stdout,"    REQUIRED\n");
  68   fprintf(stdout,"    Currently accepts PGM-files, PPM-files, PNM-files, PGX-files and\n");
  69   fprintf(stdout,"    BMP-files. Binary data is written to the file (not ascii). If a PGX\n");
  70   fprintf(stdout,"    filename is given, there will be as many output files as there are\n");
  71   fprintf(stdout,"    components: an indice starting from 0 will then be appended to the\n");
  72   fprintf(stdout,"    output filename, just before the \"pgx\" extension. If a PGM filename\n");
  73   fprintf(stdout,"    is given and there are more than one component, only the first component\n");
  74   fprintf(stdout,"    will be written to the file.\n");
  75   fprintf(stdout,"  -r <reduce factor>\n");
  76   fprintf(stdout,"    Set the number of highest resolution levels to be discarded. The\n");
  77   fprintf(stdout,"    image resolution is effectively divided by 2 to the power of the\n");
  78   fprintf(stdout,"    number of discarded levels. The reduce factor is limited by the\n");
  79   fprintf(stdout,"    smallest total number of decomposition levels among tiles.\n");
  80   fprintf(stdout,"  -l <number of quality layers to decode>\n");
  81   fprintf(stdout,"    Set the maximum number of quality layers to decode. If there are\n");
  82   fprintf(stdout,"    less quality layers than the specified number, all the quality layers\n");
  83   fprintf(stdout,"    are decoded.\n");
  84   fprintf(stdout,"\n");
  85 }
  86
  87 /* -------------------------------------------------------------------------- */
  88
  89 int get_file_format(char *filename) {
  90   unsigned int i;
  91   static const char *extension[] = {"pgx", "pnm", "pgm", "ppm", "bmp", "j2k", "jp2", "jpt" };
  92   static const int format[] = { PGX_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, BMP_DFMT, J2K_CFMT, JP2_CFMT, JPT_CFMT };
  93   char * ext = strrchr(filename, '.') + 1;
  94   if(ext) {
  95     for(i = 0; i < sizeof(format)/sizeof(*format); i++) {
  96       if(strnicmp(ext, extension[i], 3) == 0) {
  97         return format[i];
  98       }
  99     }
 100   }
 101
 102   return -1;
 103 }
 104
 105 /* -------------------------------------------------------------------------- */
 106
 107 int parse_cmdline_decoder(int argc, char **argv, opj_dparameters_t *parameters) {
 108   /* parse the command line */
 109
 110   while (1) {
 111     int c = getopt(argc, argv, "i:o:r:q:f:t:n:c:b:x:p:s:d:hP:S:E:M:R:T:C:I");
 112     if (c == -1)
 113       break;
 114     switch (c) {
 115       case 'i':      /* input file */
 116       {
 117         char *infile = optarg;
 118         parameters->decod_format = get_file_format(infile);
 119         switch(parameters->decod_format) {
 120           case J2K_CFMT:
 121           case JP2_CFMT:
 122           case JPT_CFMT:
 123             break;
 124           default:
 125             fprintf(stderr,
 126               "!! Unrecognized format for infile : %s [accept only *.j2k, *.jp2, *.jpc or *.jpt] !!\n\n",
 127               infile);
 128             return 1;
 129         }
 130         strncpy(parameters->infile, infile, MAX_PATH);
 131       }
 132       break;
 133
 134         /* ----------------------------------------------------- */
 135
 136       case 'o':      /* output file */
 137       {
 138         char *outfile = optarg;
 139         parameters->cod_format = get_file_format(outfile);
 140         switch(parameters->cod_format) {
 141           case PGX_DFMT:
 142           case PXM_DFMT:
 143           case BMP_DFMT:
 144             break;
 145           default:
 146             fprintf(stderr, "Unknown output format image %s [only *.pnm, *.pgm, *.ppm, *.pgx or *.bmp]!! \n", outfile);
 147             return 1;
 148         }
 149         strncpy(parameters->outfile, outfile, MAX_PATH);
 150       }
 151       break;
 152
 153         /* ----------------------------------------------------- */
 154
 155
 156       case 'r':    /* reduce option */
 157       {
 158         sscanf(optarg, "%d", &parameters->cp_reduce);
 159       }
 160       break;
 161
 162         /* ----------------------------------------------------- */
 163
 164
 165       case 'l':    /* layering option */
 166       {
 167         sscanf(optarg, "%d", &parameters->cp_layer);
 168       }
 169       break;
 170
 171         /* ----------------------------------------------------- */
 172
 173       case 'h':       /* display an help description */
 174         decode_help_display();
 175         return 1;
 176
 177         /* ----------------------------------------------------- */
 178
 179       default:
 180         fprintf(stderr,"WARNING -> this option is not valid \"-%c %s\"\n",c, optarg);
 181         break;
 182     }
 183   }
 184
 185   /* check for possible errors */
 186
 187   if((parameters->infile[0] == 0) || (parameters->outfile[0] == 0)) {
 188     fprintf(stderr,"ERROR -> At least one required argument is missing\nCheck j2k_to_image -h for usage information\n");
 189     return 1;
 190   }
 191
 192   return 0;
 193 }
 194
 195 /* -------------------------------------------------------------------------- */
 196
 197 /**
 198 sample error callback expecting a FILE* client object
 199 */
 200 void error_callback(const char *msg, void *client_data) {
 201   FILE *stream = (FILE*)client_data;
 202   fprintf(stream, "[ERROR] %s", msg);
 203 }
 204 /**
 205 sample warning callback expecting a FILE* client object
 206 */
 207 void warning_callback(const char *msg, void *client_data) {
 208   FILE *stream = (FILE*)client_data;
 209   fprintf(stream, "[WARNING] %s", msg);
 210 }
 211 /**
 212 sample debug callback expecting no client object
 213 */
 214 void info_callback(const char *msg, void *client_data) {
 215   (void)client_data;
 216   fprintf(stdout, "[INFO] %s", msg);
 217 }
 218
 219 /* -------------------------------------------------------------------------- */
 220
 221 int main(int argc, char **argv) {
 222   opj_dparameters_t parameters;  /* decompression parameters */
 223   opj_event_mgr_t event_mgr;    /* event manager */
 224   opj_image_t *image = NULL;
 225   FILE *fsrc = NULL;
 226   unsigned char *src = NULL;
 227   int file_length;
 228
 229   opj_dinfo_t* dinfo = NULL;  /* handle to a decompressor */
 230   opj_cio_t *cio = NULL;
 231
 232   /* configure the event callbacks (not required) */
 233   memset(&event_mgr, 0, sizeof(opj_event_mgr_t));
 234   event_mgr.error_handler = error_callback;
 235   event_mgr.warning_handler = warning_callback;
 236   event_mgr.info_handler = info_callback;
 237
 238   /* set decoding parameters to default values */
 239   opj_set_default_decoder_parameters(&parameters);
 240
 241   /* parse input and get user decoding parameters */
 242   if(parse_cmdline_decoder(argc, argv, &parameters) == 1) {
 243     return 0;
 244   }
 245
 246   /* read the input file and put it in memory */
 247   /* ---------------------------------------- */
 248   fsrc = fopen(parameters.infile, "rb");
 249   if (!fsrc) {
 250     fprintf(stderr, "ERROR -> failed to open %s for reading\n", parameters.infile);
 251     return 1;
 252   }
 253   fseek(fsrc, 0, SEEK_END);
 254   file_length = ftell(fsrc);
 255   fseek(fsrc, 0, SEEK_SET);
 256   src = (unsigned char *) malloc(file_length);
 257   fread(src, 1, file_length, fsrc);
 258   fclose(fsrc);
 259
 260   /* decode the code-stream */
 261   /* ---------------------- */
 262
 263     switch(parameters.decod_format) {
 264     case J2K_CFMT:
 265     {
 266       /* JPEG-2000 codestream */
 267
 268       /* get a decoder handle */
 269       dinfo = opj_create_decompress(CODEC_J2K);
 270
 271       /* catch events using our callbacks and give a local context */
 272       opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);
 273
 274       /* setup the decoder decoding parameters using user parameters */
 275       opj_setup_decoder(dinfo, &parameters);
 276
 277       /* open a byte stream */
 278       cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length);
 279
 280       /* decode the stream and fill the image structure */
 281       image = opj_decode(dinfo, cio);
 282       if(!image) {
 283         fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n");
 284         opj_destroy_decompress(dinfo);
 285         opj_cio_close(cio);
 286         return 1;
 287       }
 288
 289       /* close the byte stream */
 290       opj_cio_close(cio);
 291     }
 292     break;
 293
 294     case JP2_CFMT:
 295     {
 296       /* JPEG 2000 compressed image data */
 297
 298       /* get a decoder handle */
 299       dinfo = opj_create_decompress(CODEC_JP2);
 300
 301       /* catch events using our callbacks and give a local context */
 302       opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);
 303
 304       /* setup the decoder decoding parameters using the current image and using user parameters */
 305       opj_setup_decoder(dinfo, &parameters);
 306
 307       /* open a byte stream */
 308       cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length);
 309
 310       /* decode the stream and fill the image structure */
 311       image = opj_decode(dinfo, cio);
 312       if(!image) {
 313         fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n");
 314         opj_destroy_decompress(dinfo);
 315         opj_cio_close(cio);
 316         return 1;
 317       }
 318
 319       /* close the byte stream */
 320       opj_cio_close(cio);
 321
 322     }
 323     break;
 324
 325     case JPT_CFMT:
 326     {
 327       /* JPEG 2000, JPIP */
 328
 329       /* get a decoder handle */
 330       dinfo = opj_create_decompress(CODEC_JPT);
 331
 332       /* catch events using our callbacks and give a local context */
 333       opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);
 334
 335       /* setup the decoder decoding parameters using user parameters */
 336       opj_setup_decoder(dinfo, &parameters);
 337
 338       /* open a byte stream */
 339       cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length);
 340
 341       /* decode the stream and fill the image structure */
 342       image = opj_decode(dinfo, cio);
 343       if(!image) {
 344         fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n");
 345         opj_destroy_decompress(dinfo);
 346         opj_cio_close(cio);
 347         return 1;
 348       }
 349
 350       /* close the byte stream */
 351       opj_cio_close(cio);
 352     }
 353     break;
 354
 355     default:
 356       fprintf(stderr, "ERROR -> j2k_to_image : Unknown input image format\n");
 357       return 1;
 358   }
 359
 360   /* free the memory containing the code-stream */
 361   free(src);
 362   src = NULL;
 363
 364   /* create output image */
 365   /* ------------------- */
 366
 367   switch (parameters.cod_format) {
 368     case PXM_DFMT:      /* PNM PGM PPM */
 369       imagetopnm(image, parameters.outfile);
 370       break;
 371
 372     case PGX_DFMT:      /* PGX */
 373       imagetopgx(image, parameters.outfile);
 374       break;
 375
 376     case BMP_DFMT:      /* BMP */
 377       imagetobmp(image, parameters.outfile);
 378       break;
 379   }
 380
 381   /* free remaining structures */
 382   if(dinfo) {
 383     opj_destroy_decompress(dinfo);
 384   }
 385
 386   /* free image data structure */
 387   opj_image_destroy(image);
 388
 389   return 0;
 390 }
 391