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