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