[2.0] Backport all changes since r2798 (included) from trunk
[openjpeg.git] / src / bin / jp2 / opj_decompress.c
1 /*
2  * The copyright in this software is being made available under the 2-clauses 
3  * BSD License, included below. This software may be subject to other third 
4  * party and contributor rights, including patent rights, and no such rights
5  * are granted under this license.
6  *
7  * Copyright (c) 2002-2014, Universite catholique de Louvain (UCL), Belgium
8  * Copyright (c) 2002-2014, Professor Benoit Macq
9  * Copyright (c) 2001-2003, David Janssens
10  * Copyright (c) 2002-2003, Yannick Verschueren
11  * Copyright (c) 2003-2007, Francois-Olivier Devaux 
12  * Copyright (c) 2003-2014, Antonin Descampe
13  * Copyright (c) 2005, Herve Drolon, FreeImage Team
14  * Copyright (c) 2006-2007, Parvatha Elangovan
15  * Copyright (c) 2008, 2011-2012, Centre National d'Etudes Spatiales (CNES), FR 
16  * Copyright (c) 2012, CS Systemes d'Information, France
17  * All rights reserved.
18  *
19  * Redistribution and use in source and binary forms, with or without
20  * modification, are permitted provided that the following conditions
21  * are met:
22  * 1. Redistributions of source code must retain the above copyright
23  *    notice, this list of conditions and the following disclaimer.
24  * 2. Redistributions in binary form must reproduce the above copyright
25  *    notice, this list of conditions and the following disclaimer in the
26  *    documentation and/or other materials provided with the distribution.
27  *
28  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
29  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
30  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
31  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
32  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
33  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
34  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
35  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
36  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
37  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
38  * POSSIBILITY OF SUCH DAMAGE.
39  */
40 #include "opj_apps_config.h"
41
42 #include <stdio.h>
43 #include <string.h>
44 #include <stdlib.h>
45 #include <math.h>
46
47 #ifdef _WIN32
48 #include "windirent.h"
49 #else
50 #include <dirent.h>
51 #endif /* _WIN32 */
52
53 #ifdef _WIN32
54 #include <windows.h>
55 #define strcasecmp _stricmp
56 #define strncasecmp _strnicmp
57 #else
58 #include <strings.h>
59 #endif /* _WIN32 */
60
61 #include "openjpeg.h"
62 #include "opj_getopt.h"
63 #include "convert.h"
64 #include "index.h"
65
66 #ifdef OPJ_HAVE_LIBLCMS2
67 #include <lcms2.h>
68 #endif
69 #ifdef OPJ_HAVE_LIBLCMS1
70 #include <lcms.h>
71 #endif
72 #include "color.h"
73
74 #include "format_defs.h"
75
76 typedef struct dircnt{
77         /** Buffer for holding images read from Directory*/
78         char *filename_buf;
79         /** Pointer to the buffer*/
80         char **filename;
81 }dircnt_t;
82
83
84 typedef struct img_folder{
85         /** The directory path of the folder containing input images*/
86         char *imgdirpath;
87         /** Output format*/
88         const char *out_format;
89         /** Enable option*/
90         char set_imgdir;
91         /** Enable Cod Format for output*/
92         char set_out_format;
93
94 }img_fol_t;
95
96 /* -------------------------------------------------------------------------- */
97 /* Declarations                                                               */
98 int get_num_images(char *imgdirpath);
99 int load_images(dircnt_t *dirptr, char *imgdirpath);
100 int get_file_format(const char *filename);
101 char get_next_file(int imageno,dircnt_t *dirptr,img_fol_t *img_fol, opj_dparameters_t *parameters);
102 static int infile_format(const char *fname);
103
104 int parse_cmdline_decoder(int argc, char **argv, opj_dparameters_t *parameters,img_fol_t *img_fol, char *indexfilename);
105 int parse_DA_values( char* inArg, unsigned int *DA_x0, unsigned int *DA_y0, unsigned int *DA_x1, unsigned int *DA_y1);
106
107 /* -------------------------------------------------------------------------- */
108 static void decode_help_display(void) {
109         fprintf(stdout,"HELP for opj_decompress\n----\n\n");
110         fprintf(stdout,"- the -h option displays this help information on screen\n\n");
111
112 /* UniPG>> */
113         fprintf(stdout,"List of parameters for the JPEG 2000 "
114 #ifdef USE_JPWL
115                 "+ JPWL "
116 #endif /* USE_JPWL */
117                 "decoder:\n");
118 /* <<UniPG */
119         fprintf(stdout,"\n");
120         fprintf(stdout,"\n");
121         fprintf(stdout,"  -ImgDir \n");
122         fprintf(stdout,"        Image file Directory path \n");
123         fprintf(stdout,"  -OutFor \n");
124         fprintf(stdout,"    REQUIRED only if -ImgDir is used\n");
125         fprintf(stdout,"          Need to specify only format without filename <BMP>  \n");
126         fprintf(stdout,"    Currently accepts PGM, PPM, PNM, PGX, PNG, BMP, TIF, RAW and TGA formats\n");
127         fprintf(stdout,"  -i <compressed file>\n");
128         fprintf(stdout,"    REQUIRED only if an Input image directory not specified\n");
129         fprintf(stdout,"    Currently accepts J2K-files, JP2-files and JPT-files. The file type\n");
130         fprintf(stdout,"    is identified based on its suffix.\n");
131         fprintf(stdout,"  -o <decompressed file>\n");
132         fprintf(stdout,"    REQUIRED\n");
133         fprintf(stdout,"    Currently accepts PGM, PPM, PNM, PGX, PNG, BMP, TIF, RAW and TGA files\n");
134         fprintf(stdout,"    Binary data is written to the file (not ascii). If a PGX\n");
135         fprintf(stdout,"    filename is given, there will be as many output files as there are\n");
136         fprintf(stdout,"    components: an indice starting from 0 will then be appended to the\n");
137         fprintf(stdout,"    output filename, just before the \"pgx\" extension. If a PGM filename\n");
138         fprintf(stdout,"    is given and there are more than one component, only the first component\n");
139         fprintf(stdout,"    will be written to the file.\n");
140         fprintf(stdout,"  -r <reduce factor>\n");
141         fprintf(stdout,"    Set the number of highest resolution levels to be discarded. The\n");
142         fprintf(stdout,"    image resolution is effectively divided by 2 to the power of the\n");
143         fprintf(stdout,"    number of discarded levels. The reduce factor is limited by the\n");
144         fprintf(stdout,"    smallest total number of decomposition levels among tiles.\n");
145         fprintf(stdout,"  -l <number of quality layers to decode>\n");
146         fprintf(stdout,"    Set the maximum number of quality layers to decode. If there are\n");
147         fprintf(stdout,"    less quality layers than the specified number, all the quality layers\n");
148         fprintf(stdout,"    are decoded.\n");
149         fprintf(stdout,"  -x  \n"); 
150         fprintf(stdout,"    Create an index file *.Idx (-x index_name.Idx) \n");
151         fprintf(stdout,"  -d <x0,y0,x1,y1>\n");
152         fprintf(stdout,"    OPTIONAL\n");
153         fprintf(stdout,"    Decoding area\n");
154         fprintf(stdout,"    By default all the image is decoded.\n");
155         fprintf(stdout,"  -t <tile_number>\n");
156         fprintf(stdout,"    OPTIONAL\n");
157         fprintf(stdout,"    Set the tile number of the decoded tile. Follow the JPEG2000 convention from left-up to bottom-up\n");
158         fprintf(stdout,"    By default all tiles are decoded.\n");
159         fprintf(stdout,"\n");
160 /* UniPG>> */
161 #ifdef USE_JPWL
162         fprintf(stdout,"  -W <options>\n");
163         fprintf(stdout,"    Activates the JPWL correction capability, if the codestream complies.\n");
164         fprintf(stdout,"    Options can be a comma separated list of <param=val> tokens:\n");
165         fprintf(stdout,"    c, c=numcomps\n");
166         fprintf(stdout,"       numcomps is the number of expected components in the codestream\n");
167         fprintf(stdout,"       (search of first EPB rely upon this, default is %d)\n", JPWL_EXPECTED_COMPONENTS);
168 #endif /* USE_JPWL */
169 /* <<UniPG */
170         fprintf(stdout,"\n");
171 }
172
173 /* -------------------------------------------------------------------------- */
174
175 int get_num_images(char *imgdirpath){
176         DIR *dir;
177         struct dirent* content; 
178         int num_images = 0;
179
180         /*Reading the input images from given input directory*/
181
182         dir= opendir(imgdirpath);
183         if(!dir){
184                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
185                 return 0;
186         }
187         
188         while((content=readdir(dir))!=NULL){
189                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
190                         continue;
191                 num_images++;
192         }
193         return num_images;
194 }
195
196 /* -------------------------------------------------------------------------- */
197 int load_images(dircnt_t *dirptr, char *imgdirpath){
198         DIR *dir;
199         struct dirent* content; 
200         int i = 0;
201
202         /*Reading the input images from given input directory*/
203
204         dir= opendir(imgdirpath);
205         if(!dir){
206                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
207                 return 1;
208         }else   {
209                 fprintf(stderr,"Folder opened successfully\n");
210         }
211         
212         while((content=readdir(dir))!=NULL){
213                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
214                         continue;
215
216                 strcpy(dirptr->filename[i],content->d_name);
217                 i++;
218         }
219         return 0;       
220 }
221
222 /* -------------------------------------------------------------------------- */
223 int get_file_format(const char *filename) {
224         unsigned int i;
225         static const char *extension[] = {"pgx", "pnm", "pgm", "ppm", "bmp","tif", "raw", "rawl", "tga", "png", "j2k", "jp2", "jpt", "j2c", "jpc" };
226         static const int format[] = { PGX_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, BMP_DFMT, TIF_DFMT, RAW_DFMT, RAWL_DFMT, TGA_DFMT, PNG_DFMT, J2K_CFMT, JP2_CFMT, JPT_CFMT, J2K_CFMT, J2K_CFMT };
227         char * ext = strrchr(filename, '.');
228         if (ext == NULL)
229                 return -1;
230         ext++;
231         if(*ext) {
232                 for(i = 0; i < sizeof(format)/sizeof(*format); i++) {
233                         if(strcasecmp(ext, extension[i]) == 0) {
234                                 return format[i];
235                         }
236                 }
237         }
238
239         return -1;
240 }
241
242 /* -------------------------------------------------------------------------- */
243 char get_next_file(int imageno,dircnt_t *dirptr,img_fol_t *img_fol, opj_dparameters_t *parameters){
244         char image_filename[OPJ_PATH_LEN], infilename[OPJ_PATH_LEN],outfilename[OPJ_PATH_LEN],temp_ofname[OPJ_PATH_LEN];
245         char *temp_p, temp1[OPJ_PATH_LEN]="";
246
247         strcpy(image_filename,dirptr->filename[imageno]);
248         fprintf(stderr,"File Number %d \"%s\"\n",imageno,image_filename);
249         parameters->decod_format = infile_format(image_filename);
250         if (parameters->decod_format == -1)
251                 return 1;
252         sprintf(infilename,"%s/%s",img_fol->imgdirpath,image_filename);
253         strncpy(parameters->infile, infilename, sizeof(infilename));
254
255         /*Set output file*/
256         strcpy(temp_ofname,strtok(image_filename,"."));
257         while((temp_p = strtok(NULL,".")) != NULL){
258                 strcat(temp_ofname,temp1);
259                 sprintf(temp1,".%s",temp_p);
260         }
261         if(img_fol->set_out_format==1){
262                 sprintf(outfilename,"%s/%s.%s",img_fol->imgdirpath,temp_ofname,img_fol->out_format);
263                 strncpy(parameters->outfile, outfilename, sizeof(outfilename));
264         }
265         return 0;
266 }
267
268 /* -------------------------------------------------------------------------- */
269 #define JP2_RFC3745_MAGIC "\x00\x00\x00\x0c\x6a\x50\x20\x20\x0d\x0a\x87\x0a"
270 #define JP2_MAGIC "\x0d\x0a\x87\x0a"
271 /* position 45: "\xff\x52" */
272 #define J2K_CODESTREAM_MAGIC "\xff\x4f\xff\x51"
273
274 static int infile_format(const char *fname)
275 {
276         FILE *reader;
277         const char *s, *magic_s;
278         int ext_format, magic_format;
279         unsigned char buf[12];
280         OPJ_SIZE_T l_nb_read;
281
282         reader = fopen(fname, "rb");
283
284         if (reader == NULL)
285                 return -2;
286
287         memset(buf, 0, 12);
288         l_nb_read = fread(buf, 1, 12, reader);
289         fclose(reader);
290         if (l_nb_read != 12)
291                 return -1;
292
293
294
295         ext_format = get_file_format(fname);
296
297         if (ext_format == JPT_CFMT)
298                 return JPT_CFMT;
299
300         if (memcmp(buf, JP2_RFC3745_MAGIC, 12) == 0 || memcmp(buf, JP2_MAGIC, 4) == 0) {
301                 magic_format = JP2_CFMT;
302                 magic_s = ".jp2";
303         }
304         else if (memcmp(buf, J2K_CODESTREAM_MAGIC, 4) == 0) {
305                 magic_format = J2K_CFMT;
306                 magic_s = ".j2k or .jpc or .j2c";
307         }
308         else
309                 return -1;
310
311         if (magic_format == ext_format)
312                 return ext_format;
313
314         s = fname + strlen(fname) - 4;
315
316         fputs("\n===========================================\n", stderr);
317         fprintf(stderr, "The extension of this file is incorrect.\n"
318                                         "FOUND %s. SHOULD BE %s\n", s, magic_s);
319         fputs("===========================================\n", stderr);
320
321         return magic_format;
322 }
323
324 /* -------------------------------------------------------------------------- */
325 /**
326  * Parse the command line
327  */
328 /* -------------------------------------------------------------------------- */
329 int parse_cmdline_decoder(int argc, char **argv, opj_dparameters_t *parameters,img_fol_t *img_fol, char *indexfilename) {
330         /* parse the command line */
331         int totlen, c;
332         opj_option_t long_option[]={
333                 {"ImgDir",REQ_ARG, NULL ,'y'},
334                 {"OutFor",REQ_ARG, NULL ,'O'},
335         };
336
337         const char optlist[] = "i:o:r:l:x:d:t:"
338
339 /* UniPG>> */
340 #ifdef USE_JPWL
341                                         "W:"
342 #endif /* USE_JPWL */
343 /* <<UniPG */
344                         "h"             ;
345         totlen=sizeof(long_option);
346         img_fol->set_out_format = 0;
347         do {
348                 c = opj_getopt_long(argc, argv,optlist,long_option,totlen);
349                 if (c == -1)
350                         break;
351                 switch (c) {
352                         case 'i':                       /* input file */
353                         {
354                                 char *infile = opj_optarg;
355                                 parameters->decod_format = infile_format(infile);
356                                 switch(parameters->decod_format) {
357                                         case J2K_CFMT:
358                                                 break;
359                                         case JP2_CFMT:
360                                                 break;
361                                         case JPT_CFMT:
362                                                 break;
363                                         case -2:
364                                                 fprintf(stderr, 
365                                                         "!! infile cannot be read: %s !!\n\n", 
366                                                         infile);
367                                                 return 1;
368                                         default:
369                                                 fprintf(stderr, 
370                                                         "!! Unrecognized format for infile: %s [accept only *.j2k, *.jp2, *.jpc or *.jpt] !!\n\n", 
371                                                         infile);
372                                                 return 1;
373                                 }
374                                 strncpy(parameters->infile, infile, sizeof(parameters->infile)-1);
375                         }
376                         break;
377                                 
378                                 /* ----------------------------------------------------- */
379
380                         case 'o':                       /* output file */
381                         {
382                                 char *outfile = opj_optarg;
383                                 parameters->cod_format = get_file_format(outfile);
384                                 switch(parameters->cod_format) {
385                                         case PGX_DFMT:
386                                                 break;
387                                         case PXM_DFMT:
388                                                 break;
389                                         case BMP_DFMT:
390                                                 break;
391                                         case TIF_DFMT:
392                                                 break;
393                                         case RAW_DFMT:
394                                                 break;
395                                         case RAWL_DFMT:
396                                                 break;
397                                         case TGA_DFMT:
398                                                 break;
399                                         case PNG_DFMT:
400                                                 break;
401                                         default:
402                                                 fprintf(stderr, "Unknown output format image %s [only *.pnm, *.pgm, *.ppm, *.pgx, *.bmp, *.tif, *.raw or *.tga]!! \n", outfile);
403                                                 return 1;
404                                 }
405                                 strncpy(parameters->outfile, outfile, sizeof(parameters->outfile)-1);
406                         }
407                         break;
408                         
409                                 /* ----------------------------------------------------- */
410
411                         case 'O':                       /* output format */
412                         {
413                                 char outformat[50];
414                                 char *of = opj_optarg;
415                                 sprintf(outformat,".%s",of);
416                                 img_fol->set_out_format = 1;
417                                 parameters->cod_format = get_file_format(outformat);
418                                 switch(parameters->cod_format) {
419                                         case PGX_DFMT:
420                                                 img_fol->out_format = "pgx";
421                                                 break;
422                                         case PXM_DFMT:
423                                                 img_fol->out_format = "ppm";
424                                                 break;
425                                         case BMP_DFMT:
426                                                 img_fol->out_format = "bmp";
427                                                 break;
428                                         case TIF_DFMT:
429                                                 img_fol->out_format = "tif";
430                                                 break;
431                                         case RAW_DFMT:
432                                                 img_fol->out_format = "raw";
433                                                 break;
434                                         case RAWL_DFMT:
435                                                 img_fol->out_format = "rawl";
436                                                 break;
437                                         case TGA_DFMT:
438                                                 img_fol->out_format = "raw";
439                                                 break;
440                                         case PNG_DFMT:
441                                                 img_fol->out_format = "png";
442                                                 break;
443                                         default:
444                                                 fprintf(stderr, "Unknown output format image %s [only *.pnm, *.pgm, *.ppm, *.pgx, *.bmp, *.tif, *.raw or *.tga]!! \n", outformat);
445                                                 return 1;
446                                                 break;
447                                 }
448                         }
449                         break;
450
451                                 /* ----------------------------------------------------- */
452
453
454                         case 'r':               /* reduce option */
455                         {
456                                 sscanf(opj_optarg, "%ud", &parameters->cp_reduce);
457                         }
458                         break;
459                         
460                                 /* ----------------------------------------------------- */
461       
462
463                         case 'l':               /* layering option */
464                         {
465                                 sscanf(opj_optarg, "%ud", &parameters->cp_layer);
466                         }
467                         break;
468                         
469                                 /* ----------------------------------------------------- */
470
471                         case 'h':                       /* display an help description */
472                                 decode_help_display();
473                                 return 1;                               
474
475                                 /* ------------------------------------------------------ */
476
477                         case 'y':                       /* Image Directory path */
478                                 {
479                                         img_fol->imgdirpath = (char*)malloc(strlen(opj_optarg) + 1);
480                                         strcpy(img_fol->imgdirpath,opj_optarg);
481                                         img_fol->set_imgdir=1;
482                                 }
483                                 break;
484
485                                 /* ----------------------------------------------------- */
486
487                         case 'd':               /* Input decode ROI */
488                         {
489                                 int size_optarg = (int)strlen(opj_optarg) + 1;
490                                 char *ROI_values = (char*) malloc((size_t)size_optarg);
491                                 ROI_values[0] = '\0';
492                                 strncpy(ROI_values, opj_optarg, strlen(opj_optarg));
493                                 ROI_values[strlen(opj_optarg)] = '\0';
494                                 /*printf("ROI_values = %s [%d / %d]\n", ROI_values, strlen(ROI_values), size_optarg ); */
495                                 parse_DA_values( ROI_values, &parameters->DA_x0, &parameters->DA_y0, &parameters->DA_x1, &parameters->DA_y1);
496
497                                 free(ROI_values);
498                         }
499                         break;
500
501                         /* ----------------------------------------------------- */
502
503                         case 't':               /* Input tile index */
504                         {
505                                 sscanf(opj_optarg, "%ud", &parameters->tile_index);
506                                 parameters->nb_tile_to_decode = 1;
507                         }
508                         break;
509
510                                 /* ----------------------------------------------------- */                                                             
511
512                         case 'x':                       /* Creation of index file */
513                                 {
514                                         char *index = opj_optarg;
515                                         strncpy(indexfilename, index, OPJ_PATH_LEN);
516                                 }
517                                 break;
518                                 /* ----------------------------------------------------- */
519                                 /* UniPG>> */
520 #ifdef USE_JPWL
521                         
522                         case 'W':                       /* activate JPWL correction */
523                         {
524                                 char *token = NULL;
525
526                                 token = strtok(opj_optarg, ",");
527                                 while(token != NULL) {
528
529                                         /* search expected number of components */
530                                         if (*token == 'c') {
531
532                                                 static int compno;
533
534                                                 compno = JPWL_EXPECTED_COMPONENTS; /* predefined no. of components */
535
536                                                 if(sscanf(token, "c=%d", &compno) == 1) {
537                                                         /* Specified */
538                                                         if ((compno < 1) || (compno > 256)) {
539                                                                 fprintf(stderr, "ERROR -> invalid number of components c = %d\n", compno);
540                                                                 return 1;
541                                                         }
542                                                         parameters->jpwl_exp_comps = compno;
543
544                                                 } else if (!strcmp(token, "c")) {
545                                                         /* default */
546                                                         parameters->jpwl_exp_comps = compno; /* auto for default size */
547
548                                                 } else {
549                                                         fprintf(stderr, "ERROR -> invalid components specified = %s\n", token);
550                                                         return 1;
551                                                 };
552                                         }
553
554                                         /* search maximum number of tiles */
555                                         if (*token == 't') {
556
557                                                 static int tileno;
558
559                                                 tileno = JPWL_MAXIMUM_TILES; /* maximum no. of tiles */
560
561                                                 if(sscanf(token, "t=%d", &tileno) == 1) {
562                                                         /* Specified */
563                                                         if ((tileno < 1) || (tileno > JPWL_MAXIMUM_TILES)) {
564                                                                 fprintf(stderr, "ERROR -> invalid number of tiles t = %d\n", tileno);
565                                                                 return 1;
566                                                         }
567                                                         parameters->jpwl_max_tiles = tileno;
568
569                                                 } else if (!strcmp(token, "t")) {
570                                                         /* default */
571                                                         parameters->jpwl_max_tiles = tileno; /* auto for default size */
572
573                                                 } else {
574                                                         fprintf(stderr, "ERROR -> invalid tiles specified = %s\n", token);
575                                                         return 1;
576                                                 };
577                                         }
578
579                                         /* next token or bust */
580                                         token = strtok(NULL, ",");
581                                 };
582                                 parameters->jpwl_correct = OPJ_TRUE;
583                                 fprintf(stdout, "JPWL correction capability activated\n");
584                                 fprintf(stdout, "- expecting %d components\n", parameters->jpwl_exp_comps);
585                         }
586                         break;  
587 #endif /* USE_JPWL */
588 /* <<UniPG */            
589
590                                 /* ----------------------------------------------------- */
591                         
592                         default:
593                                 fprintf(stderr,"WARNING -> this option is not valid \"-%c %s\"\n",c, opj_optarg);
594                                 break;
595                 }
596         }while(c != -1);
597
598         /* check for possible errors */
599         if(img_fol->set_imgdir==1){
600                 if(!(parameters->infile[0]==0)){
601                         fprintf(stderr, "Error: options -ImgDir and -i cannot be used together !!\n");
602                         return 1;
603                 }
604                 if(img_fol->set_out_format == 0){
605                         fprintf(stderr, "Error: When -ImgDir is used, -OutFor <FORMAT> must be used !!\n");
606                         fprintf(stderr, "Only one format allowed! Valid format PGM, PPM, PNM, PGX, BMP, TIF, RAW and TGA!!\n");
607                         return 1;
608                 }
609                 if(!((parameters->outfile[0] == 0))){
610                         fprintf(stderr, "Error: options -ImgDir and -o cannot be used together !!\n");
611                         return 1;
612                 }
613         }else{
614                 if((parameters->infile[0] == 0) || (parameters->outfile[0] == 0)) {
615                         fprintf(stderr, "Example: %s -i image.j2k -o image.pgm\n",argv[0]);
616                         fprintf(stderr, "    Try: %s -h\n",argv[0]);
617                         return 1;
618                 }
619         }
620
621         return 0;
622 }
623
624 /* -------------------------------------------------------------------------- */
625 /**
626  * Parse decoding area input values
627  * separator = ","
628  */
629 /* -------------------------------------------------------------------------- */
630 int parse_DA_values( char* inArg, unsigned int *DA_x0, unsigned int *DA_y0, unsigned int *DA_x1, unsigned int *DA_y1)
631 {
632         int it = 0;
633         int values[4];
634         char delims[] = ",";
635         char *result = NULL;
636         result = strtok( inArg, delims );
637
638         while( (result != NULL) && (it < 4 ) ) {
639                 values[it] = atoi(result);
640                 result = strtok( NULL, delims );
641                 it++;
642         }
643
644         if (it != 4) {
645                 return EXIT_FAILURE;
646         }
647         else{
648                 *DA_x0 = (OPJ_UINT32)values[0]; *DA_y0 = (OPJ_UINT32)values[1];
649                 *DA_x1 = (OPJ_UINT32)values[2]; *DA_y1 = (OPJ_UINT32)values[3];
650                 return EXIT_SUCCESS;
651         }
652 }
653
654 /* -------------------------------------------------------------------------- */
655
656 /**
657 sample error callback expecting a FILE* client object
658 */
659 static void error_callback(const char *msg, void *client_data) {
660         (void)client_data;
661         fprintf(stdout, "[ERROR] %s", msg);
662 }
663 /**
664 sample warning callback expecting a FILE* client object
665 */
666 static void warning_callback(const char *msg, void *client_data) {
667         (void)client_data;
668         fprintf(stdout, "[WARNING] %s", msg);
669 }
670 /**
671 sample debug callback expecting no client object
672 */
673 static void info_callback(const char *msg, void *client_data) {
674         (void)client_data;
675         fprintf(stdout, "[INFO] %s", msg);
676 }
677
678 /* -------------------------------------------------------------------------- */
679 /**
680  * OPJ_DECOMPRESS MAIN
681  */
682 /* -------------------------------------------------------------------------- */
683 int main(int argc, char **argv)
684 {
685         FILE *fsrc = NULL;
686
687         opj_dparameters_t parameters;                   /* decompression parameters */
688         opj_image_t* image = NULL;
689         opj_stream_t *l_stream = NULL;                          /* Stream */
690         opj_codec_t* l_codec = NULL;                            /* Handle to a decompressor */
691         opj_codestream_index_t* cstr_index = NULL;
692
693         char indexfilename[OPJ_PATH_LEN];       /* index file name */
694
695         OPJ_INT32 num_images, imageno;
696         img_fol_t img_fol;
697         dircnt_t *dirptr = NULL;
698   int failed = 0;
699
700         /* set decoding parameters to default values */
701         opj_set_default_decoder_parameters(&parameters);
702
703         /* FIXME Initialize indexfilename and img_fol */
704         *indexfilename = 0;
705
706         /* Initialize img_fol */
707         memset(&img_fol,0,sizeof(img_fol_t));
708
709         /* parse input and get user encoding parameters */
710         if(parse_cmdline_decoder(argc, argv, &parameters,&img_fol, indexfilename) == 1) {
711                 return EXIT_FAILURE;
712         }
713
714         /* Initialize reading of directory */
715         if(img_fol.set_imgdir==1){      
716                 int it_image;
717                 num_images=get_num_images(img_fol.imgdirpath);
718
719                 dirptr=(dircnt_t*)malloc(sizeof(dircnt_t));
720                 if(dirptr){
721                         dirptr->filename_buf = (char*)malloc((size_t)num_images*OPJ_PATH_LEN*sizeof(char));     /* Stores at max 10 image file names*/
722                         dirptr->filename = (char**) malloc((size_t)num_images*sizeof(char*));
723
724                         if(!dirptr->filename_buf){
725                                 return EXIT_FAILURE;
726                         }
727                         for(it_image=0;it_image<num_images;it_image++){
728                                 dirptr->filename[it_image] = dirptr->filename_buf + it_image*OPJ_PATH_LEN;
729                         }
730                 }
731                 if(load_images(dirptr,img_fol.imgdirpath)==1){
732                         return EXIT_FAILURE;
733                 }
734                 if (num_images==0){
735                         fprintf(stdout,"Folder is empty\n");
736                         return EXIT_FAILURE;
737                 }
738         }else{
739                 num_images=1;
740         }
741
742         /*Decoding image one by one*/
743         for(imageno = 0; imageno < num_images ; imageno++)      {
744
745                 fprintf(stderr,"\n");
746
747                 if(img_fol.set_imgdir==1){
748                         if (get_next_file(imageno, dirptr,&img_fol, &parameters)) {
749                                 fprintf(stderr,"skipping file...\n");
750                                 continue;
751                         }
752                 }
753
754                 /* read the input file and put it in memory */
755                 /* ---------------------------------------- */
756                 fsrc = fopen(parameters.infile, "rb");
757                 if (!fsrc) {
758                         fprintf(stderr, "ERROR -> failed to open %s for reading\n", parameters.infile);
759                         return EXIT_FAILURE;
760                 }
761
762                 l_stream = opj_stream_create_default_file_stream(fsrc,1);
763                 if (!l_stream){
764                         fclose(fsrc);
765                         fprintf(stderr, "ERROR -> failed to create the stream from the file\n");
766                         return EXIT_FAILURE;
767                 }
768
769                 /* decode the JPEG2000 stream */
770                 /* ---------------------- */
771
772                 switch(parameters.decod_format) {
773                         case J2K_CFMT:  /* JPEG-2000 codestream */
774                         {
775                                 /* Get a decoder handle */
776                                 l_codec = opj_create_decompress(OPJ_CODEC_J2K);
777                                 break;
778                         }
779                         case JP2_CFMT:  /* JPEG 2000 compressed image data */
780                         {
781                                 /* Get a decoder handle */
782                                 l_codec = opj_create_decompress(OPJ_CODEC_JP2);
783                                 break;
784                         }
785                         case JPT_CFMT:  /* JPEG 2000, JPIP */
786                         {
787                                 /* Get a decoder handle */
788                                 l_codec = opj_create_decompress(OPJ_CODEC_JPT);
789                                 break;
790                         }
791                         default:
792                                 fprintf(stderr, "skipping file..\n");
793                                 opj_stream_destroy(l_stream);
794                                 continue;
795                 }
796
797                 /* catch events using our callbacks and give a local context */         
798                 opj_set_info_handler(l_codec, info_callback,00);
799                 opj_set_warning_handler(l_codec, warning_callback,00);
800                 opj_set_error_handler(l_codec, error_callback,00);
801
802                 /* Setup the decoder decoding parameters using user parameters */
803                 if ( !opj_setup_decoder(l_codec, &parameters) ){
804                         fprintf(stderr, "ERROR -> opj_compress: failed to setup the decoder\n");
805                         opj_stream_destroy(l_stream);
806                         fclose(fsrc);
807                         opj_destroy_codec(l_codec);
808                         return EXIT_FAILURE;
809                 }
810
811
812                 /* Read the main header of the codestream and if necessary the JP2 boxes*/
813                 if(! opj_read_header(l_stream, l_codec, &image)){
814                         fprintf(stderr, "ERROR -> opj_decompress: failed to read the header\n");
815                         opj_stream_destroy(l_stream);
816                         fclose(fsrc);
817                         opj_destroy_codec(l_codec);
818                         opj_image_destroy(image);
819                         return EXIT_FAILURE;
820                 }
821
822                 if (!parameters.nb_tile_to_decode) {
823                         /* Optional if you want decode the entire image */
824                         if (!opj_set_decode_area(l_codec, image, (OPJ_INT32)parameters.DA_x0,
825                                         (OPJ_INT32)parameters.DA_y0, (OPJ_INT32)parameters.DA_x1, (OPJ_INT32)parameters.DA_y1)){
826                                 fprintf(stderr, "ERROR -> opj_decompress: failed to set the decoded area\n");
827                                 opj_stream_destroy(l_stream);
828                                 opj_destroy_codec(l_codec);
829                                 opj_image_destroy(image);
830                                 fclose(fsrc);
831                                 return EXIT_FAILURE;
832                         }
833
834                         /* Get the decoded image */
835                         if (!(opj_decode(l_codec, l_stream, image) && opj_end_decompress(l_codec,       l_stream))) {
836                                 fprintf(stderr,"ERROR -> opj_decompress: failed to decode image!\n");
837                                 opj_destroy_codec(l_codec);
838                                 opj_stream_destroy(l_stream);
839                                 opj_image_destroy(image);
840                                 fclose(fsrc);
841                                 return EXIT_FAILURE;
842                         }
843                 }
844                 else {
845
846                         /* It is just here to illustrate how to use the resolution after set parameters */
847                         /*if (!opj_set_decoded_resolution_factor(l_codec, 5)) {
848                                 fprintf(stderr, "ERROR -> opj_decompress: failed to set the resolution factor tile!\n");
849                                 opj_destroy_codec(l_codec);
850                                 opj_stream_destroy(l_stream);
851                                 opj_image_destroy(image);
852                                 fclose(fsrc);
853                                 return EXIT_FAILURE;
854                         }*/
855
856                         if (!opj_get_decoded_tile(l_codec, l_stream, image, parameters.tile_index)) {
857                                 fprintf(stderr, "ERROR -> opj_decompress: failed to decode tile!\n");
858                                 opj_destroy_codec(l_codec);
859                                 opj_stream_destroy(l_stream);
860                                 opj_image_destroy(image);
861                                 fclose(fsrc);
862                                 return EXIT_FAILURE;
863                         }
864                         fprintf(stdout, "tile %d is decoded!\n\n", parameters.tile_index);
865                 }
866
867                 /* Close the byte stream */
868                 opj_stream_destroy(l_stream);
869                 fclose(fsrc);
870
871                 if(image->color_space == OPJ_CLRSPC_SYCC){
872                         color_sycc_to_rgb(image); /* FIXME */
873                 }
874                 
875                 if( image->color_space != OPJ_CLRSPC_SYCC 
876                         && image->numcomps == 3 && image->comps[0].dx == image->comps[0].dy
877                         && image->comps[1].dx != 1 )
878                         image->color_space = OPJ_CLRSPC_SYCC;
879                 else if (image->numcomps <= 2)
880                         image->color_space = OPJ_CLRSPC_GRAY;
881
882                 if(image->icc_profile_buf) {
883 #if defined(OPJ_HAVE_LIBLCMS1) || defined(OPJ_HAVE_LIBLCMS2)
884                         color_apply_icc_profile(image); /* FIXME */
885 #endif
886                         free(image->icc_profile_buf);
887                         image->icc_profile_buf = NULL; image->icc_profile_len = 0;
888                 }
889
890                 /* create output image */
891                 /* ------------------- */
892                 switch (parameters.cod_format) {
893                 case PXM_DFMT:                  /* PNM PGM PPM */
894                         if (imagetopnm(image, parameters.outfile)) {
895                                 fprintf(stderr,"Outfile %s not generated\n",parameters.outfile);
896         failed = 1;
897                         }
898                         else {
899                                 fprintf(stdout,"Generated Outfile %s\n",parameters.outfile);
900                         }
901                         break;
902
903                 case PGX_DFMT:                  /* PGX */
904                         if(imagetopgx(image, parameters.outfile)){
905                                 fprintf(stderr,"Outfile %s not generated\n",parameters.outfile);
906         failed = 1;
907                         }
908                         else {
909                                 fprintf(stdout,"Generated Outfile %s\n",parameters.outfile);
910                         }
911                         break;
912
913                 case BMP_DFMT:                  /* BMP */
914                         if(imagetobmp(image, parameters.outfile)){
915                                 fprintf(stderr,"Outfile %s not generated\n",parameters.outfile);
916         failed = 1;
917                         }
918                         else {
919                                 fprintf(stdout,"Generated Outfile %s\n",parameters.outfile);
920                         }
921                         break;
922 #ifdef OPJ_HAVE_LIBTIFF
923                 case TIF_DFMT:                  /* TIFF */
924                         if(imagetotif(image, parameters.outfile)){
925                                 fprintf(stderr,"Outfile %s not generated\n",parameters.outfile);
926         failed = 1;
927                         }
928                         else {
929                                 fprintf(stdout,"Generated Outfile %s\n",parameters.outfile);
930                         }
931                         break;
932 #endif /* OPJ_HAVE_LIBTIFF */
933                 case RAW_DFMT:                  /* RAW */
934                         if(imagetoraw(image, parameters.outfile)){
935                                 fprintf(stderr,"Error generating raw file. Outfile %s not generated\n",parameters.outfile);
936         failed = 1;
937                         }
938                         else {
939                                 fprintf(stdout,"Successfully generated Outfile %s\n",parameters.outfile);
940                         }
941                         break;
942
943                 case RAWL_DFMT:                 /* RAWL */
944                         if(imagetorawl(image, parameters.outfile)){
945                                 fprintf(stderr,"Error generating rawl file. Outfile %s not generated\n",parameters.outfile);
946         failed = 1;
947                         }
948                         else {
949                                 fprintf(stdout,"Successfully generated Outfile %s\n",parameters.outfile);
950                         }
951                         break;
952
953                 case TGA_DFMT:                  /* TGA */
954                         if(imagetotga(image, parameters.outfile)){
955                                 fprintf(stderr,"Error generating tga file. Outfile %s not generated\n",parameters.outfile);
956         failed = 1;
957                         }
958                         else {
959                                 fprintf(stdout,"Successfully generated Outfile %s\n",parameters.outfile);
960                         }
961                         break;
962 #ifdef OPJ_HAVE_LIBPNG
963                 case PNG_DFMT:                  /* PNG */
964                         if(imagetopng(image, parameters.outfile)){
965                                 fprintf(stderr,"Error generating png file. Outfile %s not generated\n",parameters.outfile);
966         failed = 1;
967                         }
968                         else {
969                                 fprintf(stdout,"Successfully generated Outfile %s\n",parameters.outfile);
970                         }
971                         break;
972 #endif /* OPJ_HAVE_LIBPNG */
973 /* Can happen if output file is TIFF or PNG
974  * and OPJ_HAVE_LIBTIF or OPJ_HAVE_LIBPNG is undefined
975 */
976                         default:
977                                 fprintf(stderr,"Outfile %s not generated\n",parameters.outfile);
978         failed = 1;
979                 }
980
981                 /* free remaining structures */
982                 if (l_codec) {
983                         opj_destroy_codec(l_codec);
984                 }
985
986
987                 /* free image data structure */
988                 opj_image_destroy(image);
989
990                 /* destroy the codestream index */
991                 opj_destroy_cstr_index(&cstr_index);
992
993         }
994         return failed ? EXIT_FAILURE : EXIT_SUCCESS;
995 }
996 /*end main*/
997
998
999
1000