src/bin/jp2/opj_compress.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, Jerome Fimes, Communications & Systemes <jerome.fimes@c-s.fr>
  10  * Copyright (c) 2011-2012, Centre National d'Etudes Spatiales (CNES), France
  11  * Copyright (c) 2012, CS Systemes d'Information, France
  12  * All rights reserved.
  13  *
  14  * Redistribution and use in source and binary forms, with or without
  15  * modification, are permitted provided that the following conditions
  16  * are met:
  17  * 1. Redistributions of source code must retain the above copyright
  18  *    notice, this list of conditions and the following disclaimer.
  19  * 2. Redistributions in binary form must reproduce the above copyright
  20  *    notice, this list of conditions and the following disclaimer in the
  21  *    documentation and/or other materials provided with the distribution.
  22  *
  23  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
  24  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  25  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  26  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  27  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  28  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  29  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  30  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  31  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  32  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  33  * POSSIBILITY OF SUCH DAMAGE.
  34  */
  35
  36 #include <stdio.h>
  37 #include <string.h>
  38 #include <stdlib.h>
  39 #include <math.h>
  40 #include <assert.h>
  41
  42 #ifdef _WIN32
  43 #include "windirent.h"
  44 #else
  45 #include <dirent.h>
  46 #endif /* _WIN32 */
  47
  48 #ifdef _WIN32
  49 #include <windows.h>
  50 #define strcasecmp _stricmp
  51 #define strncasecmp _strnicmp
  52 #else
  53 #include <strings.h>
  54 #endif /* _WIN32 */
  55
  56 #include "opj_apps_config.h"
  57 #include "openjpeg.h"
  58 #include "opj_getopt.h"
  59 #include "convert.h"
  60 #include "index.h"
  61
  62 #include "format_defs.h"
  63
  64 #define CINEMA_24_CS 1302083    /*Codestream length for 24fps*/
  65 #define CINEMA_48_CS 651041             /*Codestream length for 48fps*/
  66 #define COMP_24_CS 1041666              /*Maximum size per color component for 2K & 4K @ 24fps*/
  67 #define COMP_48_CS 520833               /*Maximum size per color component for 2K @ 48fps*/
  68
  69 typedef struct dircnt{
  70     /** Buffer for holding images read from Directory*/
  71     char *filename_buf;
  72     /** Pointer to the buffer*/
  73     char **filename;
  74 }dircnt_t;
  75
  76 typedef struct img_folder{
  77     /** The directory path of the folder containing input images*/
  78     char *imgdirpath;
  79     /** Output format*/
  80     char *out_format;
  81     /** Enable option*/
  82     char set_imgdir;
  83     /** Enable Cod Format for output*/
  84     char set_out_format;
  85 }img_fol_t;
  86
  87 static void encode_help_display(void) {
  88     fprintf(stdout,"HELP for opj_compress\n----\n\n");
  89     fprintf(stdout,"- the -h option displays this help information on screen\n\n");
  90
  91     /* UniPG>> */
  92     fprintf(stdout,"List of parameters for the JPEG 2000 "
  93         #ifdef USE_JPWL
  94             "+ JPWL "
  95         #endif /* USE_JPWL */
  96             "encoder:\n");
  97     /* <<UniPG */
  98     fprintf(stdout,"\n");
  99     fprintf(stdout,"REMARKS:\n");
 100     fprintf(stdout,"---------\n");
 101     fprintf(stdout,"\n");
 102     fprintf(stdout,"The markers written to the main_header are : SOC SIZ COD QCD COM.\n");
 103     fprintf(stdout,"COD and QCD never appear in the tile_header.\n");
 104     fprintf(stdout,"\n");
 105     fprintf(stdout,"By default:\n");
 106     fprintf(stdout,"------------\n");
 107     fprintf(stdout,"\n");
 108     fprintf(stdout," * Lossless\n");
 109     fprintf(stdout," * 1 tile\n");
 110     fprintf(stdout," * Size of precinct : 2^15 x 2^15 (means 1 precinct)\n");
 111     fprintf(stdout," * Size of code-block : 64 x 64\n");
 112     fprintf(stdout," * Number of resolutions: 6\n");
 113     fprintf(stdout," * No SOP marker in the codestream\n");
 114     fprintf(stdout," * No EPH marker in the codestream\n");
 115     fprintf(stdout," * No sub-sampling in x or y direction\n");
 116     fprintf(stdout," * No mode switch activated\n");
 117     fprintf(stdout," * Progression order: LRCP\n");
 118     fprintf(stdout," * No index file\n");
 119     fprintf(stdout," * No ROI upshifted\n");
 120     fprintf(stdout," * No offset of the origin of the image\n");
 121     fprintf(stdout," * No offset of the origin of the tiles\n");
 122     fprintf(stdout," * Reversible DWT 5-3\n");
 123     /* UniPG>> */
 124 #ifdef USE_JPWL
 125     fprintf(stdout," * No JPWL protection\n");
 126 #endif /* USE_JPWL */
 127     /* <<UniPG */
 128     fprintf(stdout,"\n");
 129     fprintf(stdout,"Parameters:\n");
 130     fprintf(stdout,"------------\n");
 131     fprintf(stdout,"\n");
 132     fprintf(stdout,"Required Parameters (except with -h):\n");
 133     fprintf(stdout,"One of the two options -ImgDir or -i must be used\n");
 134     fprintf(stdout,"\n");
 135     fprintf(stdout,"-ImgDir      : Image file Directory path (example ../Images) \n");
 136     fprintf(stdout,"    When using this option -OutFor must be used\n");
 137     fprintf(stdout,"\n");
 138     fprintf(stdout,"-OutFor \n");
 139     fprintf(stdout,"    REQUIRED only if -ImgDir is used\n");
 140     fprintf(stdout,"      Need to specify only format without filename <BMP>  \n");
 141     fprintf(stdout,"    Currently accepts PBM, PGM, PPM, PNM, PAM, PGX, PNG, BMP, TIF, RAW, RAWL and TGA formats\n");
 142     fprintf(stdout,"\n");
 143     fprintf(stdout,"-i           : source file  (-i source.pnm also *pbm, *.pgm, *.ppm, *.pam, *.pgx, *png, *.bmp, *.tif, *.raw, *.tga) \n");
 144     fprintf(stdout,"    When using this option -o must be used\n");
 145     fprintf(stdout,"\n");
 146     fprintf(stdout,"-o           : destination file (-o dest.j2k or .jp2) \n");
 147     fprintf(stdout,"\n");
 148     fprintf(stdout,"Optional Parameters:\n");
 149     fprintf(stdout,"\n");
 150     fprintf(stdout,"-h           : display the help information \n ");
 151     fprintf(stdout,"\n");
 152     fprintf(stdout,"-cinema2K    : Digital Cinema 2K profile compliant codestream for 2K resolution.(-CINEMA2K 24 or 48) \n");
 153     fprintf(stdout,"      Need to specify the frames per second for a 2K resolution. Only 24 or 48 fps is allowed\n");
 154     fprintf(stdout,"\n");
 155     fprintf(stdout,"-cinema4K    : Digital Cinema 4K profile compliant codestream for 4K resolution \n");
 156     fprintf(stdout,"      Frames per second not required. Default value is 24fps\n");
 157     fprintf(stdout,"\n");
 158     fprintf(stdout,"-r           : different compression ratios for successive layers (-r 20,10,5)\n ");
 159     fprintf(stdout,"             - The rate specified for each quality level is the desired \n");
 160     fprintf(stdout,"               compression factor.\n");
 161     fprintf(stdout,"               Example: -r 20,10,1 means quality 1: compress 20x, \n");
 162     fprintf(stdout,"                 quality 2: compress 10x and quality 3: compress lossless\n");
 163     fprintf(stdout,"\n");
 164     fprintf(stdout,"               (options -r and -q cannot be used together)\n ");
 165     fprintf(stdout,"\n");
 166
 167     fprintf(stdout,"-q           : different psnr for successive layers (-q 30,40,50) \n ");
 168
 169     fprintf(stdout,"               (options -r and -q cannot be used together)\n ");
 170
 171     fprintf(stdout,"\n");
 172     fprintf(stdout,"-n           : number of resolutions (-n 3) \n");
 173     fprintf(stdout,"\n");
 174     fprintf(stdout,"-b           : size of code block (-b 32,32) \n");
 175     fprintf(stdout,"\n");
 176     fprintf(stdout,"-c           : size of precinct (-c 128,128) \n");
 177     fprintf(stdout,"\n");
 178     fprintf(stdout,"-t           : size of tile (-t 512,512) \n");
 179     fprintf(stdout,"\n");
 180     fprintf(stdout,"-p           : progression order (-p LRCP) [LRCP, RLCP, RPCL, PCRL, CPRL] \n");
 181     fprintf(stdout,"\n");
 182     fprintf(stdout,"-s           : subsampling factor (-s 2,2) [-s X,Y] \n");
 183     fprintf(stdout,"         Remark: subsampling bigger than 2 can produce error\n");
 184     fprintf(stdout,"\n");
 185     fprintf(stdout,"-POC         : Progression order change (-POC T1=0,0,1,5,3,CPRL/T1=5,0,1,6,3,CPRL) \n");
 186     fprintf(stdout,"      Example: T1=0,0,1,5,3,CPRL \n");
 187     fprintf(stdout,"                     : Ttilenumber=Resolution num start,Component num start,Layer num end,Resolution num end,Component num end,Progression order\n");
 188     fprintf(stdout,"\n");
 189     fprintf(stdout,"-SOP         : write SOP marker before each packet \n");
 190     fprintf(stdout,"\n");
 191     fprintf(stdout,"-EPH         : write EPH marker after each header packet \n");
 192     fprintf(stdout,"\n");
 193     fprintf(stdout,"-M           : mode switch (-M 3) [1=BYPASS(LAZY) 2=RESET 4=RESTART(TERMALL)\n");
 194     fprintf(stdout,"                 8=VSC 16=ERTERM(SEGTERM) 32=SEGMARK(SEGSYM)] \n");
 195     fprintf(stdout,"                 Indicate multiple modes by adding their values. \n");
 196     fprintf(stdout,"                 ex: RESTART(4) + RESET(2) + SEGMARK(32) = -M 38\n");
 197     fprintf(stdout,"\n");
 198     fprintf(stdout,"-TP          : divide packets of every tile into tile-parts (-TP R) [R, L, C]\n");
 199     fprintf(stdout,"\n");
 200     fprintf(stdout,"-x           : create an index file *.Idx (-x index_name.Idx) \n");
 201     fprintf(stdout,"\n");
 202     fprintf(stdout,"-ROI         : c=%%d,U=%%d : quantization indices upshifted \n");
 203     fprintf(stdout,"               for component c=%%d [%%d = 0,1,2]\n");
 204     fprintf(stdout,"               with a value of U=%%d [0 <= %%d <= 37] (i.e. -ROI c=0,U=25) \n");
 205     fprintf(stdout,"\n");
 206     fprintf(stdout,"-d           : offset of the origin of the image (-d 150,300) \n");
 207     fprintf(stdout,"\n");
 208     fprintf(stdout,"-T           : offset of the origin of the tiles (-T 100,75) \n");
 209     fprintf(stdout,"\n");
 210     fprintf(stdout,"-I           : use the irreversible DWT 9-7 (-I) \n");
 211     fprintf(stdout,"\n");
 212     fprintf(stdout,"-F           : characteristics of the raw input image\n");
 213     fprintf(stdout,"               -F rawWidth,rawHeight,rawComp,rawBitDepth,s/u (Signed/Unsigned)\n");
 214     fprintf(stdout,"               Example: -i lena.raw -o lena.j2k -F 512,512,3,8,u\n");
 215     fprintf(stdout,"\n");
 216     fprintf(stdout,"-m           : use array-based MCT, values are coma separated, line by line\n");
 217     fprintf(stdout,"                       no specific separators between lines, no space allowed between values\n");
 218     fprintf(stdout,"-jpip        : write jpip codestream index box in JP2 output file\n");
 219     fprintf(stdout,"               NOTICE: currently supports only RPCL order\n");
 220     fprintf(stdout,"\n");
 221     /* UniPG>> */
 222 #ifdef USE_JPWL
 223     fprintf(stdout,"-W           : adoption of JPWL (Part 11) capabilities (-W params)\n");
 224     fprintf(stdout,"               The parameters can be written and repeated in any order:\n");
 225     fprintf(stdout,"               [h<tilepart><=type>,s<tilepart><=method>,a=<addr>,...\n");
 226     fprintf(stdout,"                ...,z=<size>,g=<range>,p<tilepart:pack><=type>]\n");
 227     fprintf(stdout,"\n");
 228     fprintf(stdout,"                 h selects the header error protection (EPB): 'type' can be\n");
 229     fprintf(stdout,"                   [0=none 1,absent=predefined 16=CRC-16 32=CRC-32 37-128=RS]\n");
 230     fprintf(stdout,"                   if 'tilepart' is absent, it is for main and tile headers\n");
 231     fprintf(stdout,"                   if 'tilepart' is present, it applies from that tile\n");
 232     fprintf(stdout,"                     onwards, up to the next h<> spec, or to the last tilepart\n");
 233     fprintf(stdout,"                     in the codestream (max. %d specs)\n", JPWL_MAX_NO_TILESPECS);
 234     fprintf(stdout,"\n");
 235     fprintf(stdout,"                 p selects the packet error protection (EEP/UEP with EPBs)\n");
 236     fprintf(stdout,"                  to be applied to raw data: 'type' can be\n");
 237     fprintf(stdout,"                   [0=none 1,absent=predefined 16=CRC-16 32=CRC-32 37-128=RS]\n");
 238     fprintf(stdout,"                   if 'tilepart:pack' is absent, it is from tile 0, packet 0\n");
 239     fprintf(stdout,"                   if 'tilepart:pack' is present, it applies from that tile\n");
 240     fprintf(stdout,"                     and that packet onwards, up to the next packet spec\n");
 241     fprintf(stdout,"                     or to the last packet in the last tilepart in the stream\n");
 242     fprintf(stdout,"                     (max. %d specs)\n", JPWL_MAX_NO_PACKSPECS);
 243     fprintf(stdout,"\n");
 244     fprintf(stdout,"                 s enables sensitivity data insertion (ESD): 'method' can be\n");
 245     fprintf(stdout,"                   [-1=NO ESD 0=RELATIVE ERROR 1=MSE 2=MSE REDUCTION 3=PSNR\n");
 246     fprintf(stdout,"                    4=PSNR INCREMENT 5=MAXERR 6=TSE 7=RESERVED]\n");
 247     fprintf(stdout,"                   if 'tilepart' is absent, it is for main header only\n");
 248     fprintf(stdout,"                   if 'tilepart' is present, it applies from that tile\n");
 249     fprintf(stdout,"                     onwards, up to the next s<> spec, or to the last tilepart\n");
 250     fprintf(stdout,"                     in the codestream (max. %d specs)\n", JPWL_MAX_NO_TILESPECS);
 251     fprintf(stdout,"\n");
 252     fprintf(stdout,"                 g determines the addressing mode: <range> can be\n");
 253     fprintf(stdout,"                   [0=PACKET 1=BYTE RANGE 2=PACKET RANGE]\n");
 254     fprintf(stdout,"\n");
 255     fprintf(stdout,"                 a determines the size of data addressing: <addr> can be\n");
 256     fprintf(stdout,"                   2/4 bytes (small/large codestreams). If not set, auto-mode\n");
 257     fprintf(stdout,"\n");
 258     fprintf(stdout,"                 z determines the size of sensitivity values: <size> can be\n");
 259     fprintf(stdout,"                   1/2 bytes, for the transformed pseudo-floating point value\n");
 260     fprintf(stdout,"\n");
 261     fprintf(stdout,"                 ex.:\n");
 262     fprintf(stdout,"                   h,h0=64,h3=16,h5=32,p0=78,p0:24=56,p1,p3:0=0,p3:20=32,s=0,\n");
 263     fprintf(stdout,"                     s0=6,s3=-1,a=0,g=1,z=1\n");
 264     fprintf(stdout,"                 means\n");
 265     fprintf(stdout,"                   predefined EPB in MH, rs(64,32) from TPH 0 to TPH 2,\n");
 266     fprintf(stdout,"                   CRC-16 in TPH 3 and TPH 4, CRC-32 in remaining TPHs,\n");
 267     fprintf(stdout,"                   UEP rs(78,32) for packets 0 to 23 of tile 0,\n");
 268     fprintf(stdout,"                   UEP rs(56,32) for packs. 24 to the last of tilepart 0,\n");
 269     fprintf(stdout,"                   UEP rs default for packets of tilepart 1,\n");
 270     fprintf(stdout,"                   no UEP for packets 0 to 19 of tilepart 3,\n");
 271     fprintf(stdout,"                   UEP CRC-32 for packs. 20 of tilepart 3 to last tilepart,\n");
 272     fprintf(stdout,"                   relative sensitivity ESD for MH,\n");
 273     fprintf(stdout,"                   TSE ESD from TPH 0 to TPH 2, byte range with automatic\n");
 274     fprintf(stdout,"                   size of addresses and 1 byte for each sensitivity value\n");
 275     fprintf(stdout,"\n");
 276     fprintf(stdout,"                 ex.:\n");
 277     fprintf(stdout,"                       h,s,p\n");
 278     fprintf(stdout,"                 means\n");
 279     fprintf(stdout,"                   default protection to headers (MH and TPHs) as well as\n");
 280     fprintf(stdout,"                   data packets, one ESD in MH\n");
 281     fprintf(stdout,"\n");
 282     fprintf(stdout,"                 N.B.: use the following recommendations when specifying\n");
 283     fprintf(stdout,"                       the JPWL parameters list\n");
 284     fprintf(stdout,"                   - when you use UEP, always pair the 'p' option with 'h'\n");
 285     fprintf(stdout,"                 \n");
 286 #endif /* USE_JPWL */
 287     /* <<UniPG */
 288     fprintf(stdout,"IMPORTANT:\n");
 289     fprintf(stdout,"-----------\n");
 290     fprintf(stdout,"\n");
 291     fprintf(stdout,"The index file has the structure below:\n");
 292     fprintf(stdout,"---------------------------------------\n");
 293     fprintf(stdout,"\n");
 294     fprintf(stdout,"Image_height Image_width\n");
 295     fprintf(stdout,"progression order\n");
 296     fprintf(stdout,"Tiles_size_X Tiles_size_Y\n");
 297     fprintf(stdout,"Tiles_nb_X Tiles_nb_Y\n");
 298     fprintf(stdout,"Components_nb\n");
 299     fprintf(stdout,"Layers_nb\n");
 300     fprintf(stdout,"decomposition_levels\n");
 301     fprintf(stdout,"[Precincts_size_X_res_Nr Precincts_size_Y_res_Nr]...\n");
 302     fprintf(stdout,"   [Precincts_size_X_res_0 Precincts_size_Y_res_0]\n");
 303     fprintf(stdout,"Main_header_start_position\n");
 304     fprintf(stdout,"Main_header_end_position\n");
 305     fprintf(stdout,"Codestream_size\n");
 306     fprintf(stdout,"\n");
 307     fprintf(stdout,"INFO ON TILES\n");
 308     fprintf(stdout,"tileno start_pos end_hd end_tile nbparts disto nbpix disto/nbpix\n");
 309     fprintf(stdout,"Tile_0 start_pos end_Theader end_pos NumParts TotalDisto NumPix MaxMSE\n");
 310     fprintf(stdout,"Tile_1   ''           ''        ''        ''       ''    ''      ''\n");
 311     fprintf(stdout,"...\n");
 312     fprintf(stdout,"Tile_Nt   ''           ''        ''        ''       ''    ''     ''\n");
 313     fprintf(stdout,"...\n");
 314     fprintf(stdout,"TILE 0 DETAILS\n");
 315     fprintf(stdout,"part_nb tileno num_packs start_pos end_tph_pos end_pos\n");
 316     fprintf(stdout,"...\n");
 317     fprintf(stdout,"Progression_string\n");
 318     fprintf(stdout,"pack_nb tileno layno resno compno precno start_pos end_ph_pos end_pos disto\n");
 319     fprintf(stdout,"Tpacket_0 Tile layer res. comp. prec. start_pos end_pos disto\n");
 320     fprintf(stdout,"...\n");
 321     fprintf(stdout,"Tpacket_Np ''   ''    ''   ''    ''       ''       ''     ''\n");
 322
 323     fprintf(stdout,"MaxDisto\n");
 324
 325     fprintf(stdout,"TotalDisto\n\n");
 326 }
 327
 328 static OPJ_PROG_ORDER give_progression(const char progression[4]) {
 329     if(strncmp(progression, "LRCP", 4) == 0) {
 330         return OPJ_LRCP;
 331     }
 332     if(strncmp(progression, "RLCP", 4) == 0) {
 333         return OPJ_RLCP;
 334     }
 335     if(strncmp(progression, "RPCL", 4) == 0) {
 336         return OPJ_RPCL;
 337     }
 338     if(strncmp(progression, "PCRL", 4) == 0) {
 339         return OPJ_PCRL;
 340     }
 341     if(strncmp(progression, "CPRL", 4) == 0) {
 342         return OPJ_CPRL;
 343     }
 344
 345     return OPJ_PROG_UNKNOWN;
 346 }
 347
 348 static unsigned int get_num_images(char *imgdirpath){
 349     DIR *dir;
 350     struct dirent* content;
 351     unsigned int num_images = 0;
 352
 353     /*Reading the input images from given input directory*/
 354
 355     dir= opendir(imgdirpath);
 356     if(!dir){
 357         fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
 358         return 0;
 359     }
 360
 361     num_images=0;
 362     while((content=readdir(dir))!=NULL){
 363         if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
 364             continue;
 365         num_images++;
 366     }
 367     return num_images;
 368 }
 369
 370 static int load_images(dircnt_t *dirptr, char *imgdirpath){
 371     DIR *dir;
 372     struct dirent* content;
 373     int i = 0;
 374
 375     /*Reading the input images from given input directory*/
 376
 377     dir= opendir(imgdirpath);
 378     if(!dir){
 379         fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
 380         return 1;
 381     }else       {
 382         fprintf(stderr,"Folder opened successfully\n");
 383     }
 384
 385     while((content=readdir(dir))!=NULL){
 386         if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
 387             continue;
 388
 389         strcpy(dirptr->filename[i],content->d_name);
 390         i++;
 391     }
 392     return 0;
 393 }
 394
 395 static int get_file_format(char *filename) {
 396     unsigned int i;
 397     static const char *extension[] = {
 398         "pgx", "pnm", "pgm", "ppm", "pbm", "pam", "bmp", "tif", "raw", "rawl", "tga", "png", "j2k", "jp2", "j2c", "jpc"
 399     };
 400     static const int format[] = {
 401         PGX_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, BMP_DFMT, TIF_DFMT, RAW_DFMT, RAWL_DFMT, TGA_DFMT, PNG_DFMT, J2K_CFMT, JP2_CFMT, J2K_CFMT, J2K_CFMT
 402     };
 403     char * ext = strrchr(filename, '.');
 404     if (ext == NULL)
 405         return -1;
 406     ext++;
 407     for(i = 0; i < sizeof(format)/sizeof(*format); i++) {
 408         if(strcasecmp(ext, extension[i]) == 0) {
 409             return format[i];
 410         }
 411     }
 412     return -1;
 413 }
 414
 415 static char * get_file_name(char *name){
 416     char *fname;
 417     fname= (char*)malloc(OPJ_PATH_LEN*sizeof(char));
 418     fname= strtok(name,".");
 419     return fname;
 420 }
 421
 422 static char get_next_file(int imageno,dircnt_t *dirptr,img_fol_t *img_fol, opj_cparameters_t *parameters){
 423     char image_filename[OPJ_PATH_LEN], infilename[OPJ_PATH_LEN],outfilename[OPJ_PATH_LEN],temp_ofname[OPJ_PATH_LEN];
 424     char *temp_p, temp1[OPJ_PATH_LEN]="";
 425
 426     strcpy(image_filename,dirptr->filename[imageno]);
 427     fprintf(stderr,"File Number %d \"%s\"\n",imageno,image_filename);
 428     parameters->decod_format = get_file_format(image_filename);
 429     if (parameters->decod_format == -1)
 430         return 1;
 431     sprintf(infilename,"%s/%s",img_fol->imgdirpath,image_filename);
 432     strncpy(parameters->infile, infilename, sizeof(infilename));
 433
 434     /*Set output file*/
 435     strcpy(temp_ofname,get_file_name(image_filename));
 436     while((temp_p = strtok(NULL,".")) != NULL){
 437         strcat(temp_ofname,temp1);
 438         sprintf(temp1,".%s",temp_p);
 439     }
 440     if(img_fol->set_out_format==1){
 441         sprintf(outfilename,"%s/%s.%s",img_fol->imgdirpath,temp_ofname,img_fol->out_format);
 442         strncpy(parameters->outfile, outfilename, sizeof(outfilename));
 443     }
 444     return 0;
 445 }
 446
 447 static int initialise_4K_poc(opj_poc_t *POC, int numres){
 448     POC[0].tile  = 1;
 449     POC[0].resno0  = 0;
 450     POC[0].compno0 = 0;
 451     POC[0].layno1  = 1;
 452     POC[0].resno1  = numres-1;
 453     POC[0].compno1 = 3;
 454     POC[0].prg1 = OPJ_CPRL;
 455     POC[1].tile  = 1;
 456     POC[1].resno0  = numres-1;
 457     POC[1].compno0 = 0;
 458     POC[1].layno1  = 1;
 459     POC[1].resno1  = numres;
 460     POC[1].compno1 = 3;
 461     POC[1].prg1 = OPJ_CPRL;
 462     return 2;
 463 }
 464
 465 static void set_cinema_parameters(opj_cparameters_t *parameters)
 466 {
 467     /* Configure cinema parameters */
 468     float max_rate = 0;
 469     float temp_rate = 0;
 470     int i;
 471
 472     /* profile (Rsiz) */
 473     switch (parameters->cp_cinema){
 474     case OPJ_CINEMA2K_24:
 475     case OPJ_CINEMA2K_48:
 476         parameters->cp_rsiz = OPJ_CINEMA2K;
 477         break;
 478     case OPJ_CINEMA4K_24:
 479         parameters->cp_rsiz = OPJ_CINEMA4K;
 480         break;
 481     }
 482
 483     /* No tiling */
 484     parameters->tile_size_on = OPJ_FALSE;
 485     parameters->cp_tdx=1;
 486     parameters->cp_tdy=1;
 487
 488     /* One tile part for each component */
 489     parameters->tp_flag = 'C';
 490     parameters->tp_on = 1;
 491
 492     /* Tile and Image shall be at (0,0) */
 493     parameters->cp_tx0 = 0;
 494     parameters->cp_ty0 = 0;
 495     parameters->image_offset_x0 = 0;
 496     parameters->image_offset_y0 = 0;
 497
 498     /* Codeblock size= 32*32 */
 499     parameters->cblockw_init = 32;
 500     parameters->cblockh_init = 32;
 501
 502     /* Codeblock style: no mode switch enabled */
 503     parameters->mode = 0;
 504
 505     /* No ROI */
 506     parameters->roi_compno = -1;
 507
 508     /* No subsampling */
 509     parameters->subsampling_dx = 1;
 510     parameters->subsampling_dy = 1;
 511
 512     /* 9-7 transform */
 513     parameters->irreversible = 1;
 514
 515     /* Number of layers */
 516     if (parameters->tcp_numlayers > 1){
 517         fprintf(stdout,"JPEG 2000 Profile-3 and 4 (2k/4k dc profile) requires:\n"
 518                 "1 single quality layer"
 519                 "-> Number of layers forced to 1 (rather than %d)\n",
 520                 parameters->tcp_numlayers);
 521         parameters->tcp_numlayers = 1;
 522     }
 523
 524     /* Resolution levels */
 525     switch (parameters->cp_cinema){
 526     case OPJ_CINEMA2K_24:
 527     case OPJ_CINEMA2K_48:
 528         if(parameters->numresolution > 6){
 529             fprintf(stdout,"JPEG 2000 Profile-3 (2k dc profile) requires:\n"
 530                     "Number of decomposition levels <= 5\n"
 531                     "-> Number of decomposition levels forced to 5 (rather than %d)\n",
 532                     parameters->numresolution+1);
 533             parameters->numresolution = 6;
 534         }
 535         break;
 536     case OPJ_CINEMA4K_24:
 537         if(parameters->numresolution < 2){
 538             fprintf(stdout,"JPEG 2000 Profile-4 (4k dc profile) requires:\n"
 539                     "Number of decomposition levels >= 1 && <= 6\n"
 540                     "-> Number of decomposition levels forced to 1 (rather than %d)\n",
 541                     parameters->numresolution+1);
 542             parameters->numresolution = 1;
 543         }else if(parameters->numresolution > 7){
 544             fprintf(stdout,"JPEG 2000 Profile-4 (4k dc profile) requires:\n"
 545                     "Number of decomposition levels >= 1 && <= 6\n"
 546                     "-> Number of decomposition levels forced to 6 (rather than %d)\n",
 547                     parameters->numresolution+1);
 548             parameters->numresolution = 7;
 549         }
 550         break;
 551     default :
 552         break;
 553     }
 554
 555     /* Precincts */
 556     parameters->csty |= 0x01;
 557     parameters->res_spec = parameters->numresolution-1;
 558     for (i = 0; i<parameters->res_spec; i++) {
 559         parameters->prcw_init[i] = 256;
 560         parameters->prch_init[i] = 256;
 561     }
 562
 563     /* The progression order shall be CPRL */
 564     parameters->prog_order = OPJ_CPRL;
 565
 566     /* Progression order changes for 4K, disallowed for 2K */
 567     if (parameters->cp_cinema == OPJ_CINEMA4K_24) {
 568         parameters->numpocs = initialise_4K_poc(parameters->POC,parameters->numresolution);
 569     } else {
 570         parameters->numpocs = 0;
 571     }
 572 }
 573
 574 static void set_cinema_rate(opj_cparameters_t *parameters, opj_image_t *image)
 575 {
 576     /* Limited bit-rate */
 577     float max_rate = 0;
 578     float temp_rate = 0;
 579     parameters->cp_disto_alloc = 1;
 580     switch (parameters->cp_cinema){
 581     case OPJ_CINEMA2K_24:
 582     case OPJ_CINEMA4K_24:
 583         max_rate = ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/
 584                 (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy);
 585         if (parameters->tcp_rates[0] == 0){
 586             parameters->tcp_rates[0] = max_rate;
 587         }else{
 588             temp_rate =((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/
 589                     (parameters->tcp_rates[0] * 8 * image->comps[0].dx * image->comps[0].dy);
 590             if (temp_rate > CINEMA_24_CS ){
 591                 fprintf(stdout,"JPEG 2000 Profile-3 and 4 (2k/4k dc profile) requires:\n"
 592                         "Maximum 1302083 compressed bytes @ 24fps\n"
 593                         "-> Specified rate (%3.1f) exceeds this limit. Rate will be forced to %3.1f.\n",
 594                         parameters->tcp_rates[0], max_rate);
 595                 parameters->tcp_rates[0]= max_rate;
 596             }else{
 597                 fprintf(stdout,"JPEG 2000 Profile-3 and 4 (2k/4k dc profile):\n"
 598                         "INFO : Specified rate (%3.1f) is below the 2k/4k limit @ 24fps.\n",
 599                         parameters->tcp_rates[0]);
 600             }
 601         }
 602         parameters->max_comp_size = COMP_24_CS;
 603         break;
 604     case OPJ_CINEMA2K_48:
 605         max_rate = ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/
 606                 (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy);
 607         if (parameters->tcp_rates[0] == 0){
 608             parameters->tcp_rates[0] = max_rate;
 609         }else{
 610             temp_rate =((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/
 611                     (parameters->tcp_rates[0] * 8 * image->comps[0].dx * image->comps[0].dy);
 612             if (temp_rate > CINEMA_48_CS ){
 613                 fprintf(stdout,"JPEG 2000 Profile-3 (2k dc profile) requires:\n"
 614                         "Maximum 651041 compressed bytes @ 48fps\n"
 615                         "-> Specified rate (%3.1f) exceeds this limit. Rate will be forced to %3.1f.\n",
 616                         parameters->tcp_rates[0], max_rate);
 617                 parameters->tcp_rates[0]= max_rate;
 618             }else{
 619                 fprintf(stdout,"JPEG 2000 Profile-3 (2k dc profile):\n"
 620                         "INFO : Specified rate (%3.1f) is below the 2k limit @ 48 fps.\n",
 621                         parameters->tcp_rates[0]);
 622             }
 623         }
 624         parameters->max_comp_size = COMP_48_CS;
 625         break;
 626     default:
 627         break;
 628     }
 629 }
 630
 631 static OPJ_BOOL is_cinema_compliant(opj_image_t *image, OPJ_CINEMA_MODE cinema_mode)
 632 {
 633     OPJ_UINT32 i;
 634
 635     /* Number of components */
 636     if (image->numcomps != 3){
 637         fprintf(stdout,"JPEG 2000 Profile-3 (2k dc profile) requires:\n"
 638                 "3 components"
 639                 "-> Number of components of input image (%d) is not compliant\n"
 640                 "-> Non-profile-3 codestream will be generated\n",
 641                 image->numcomps);
 642         return OPJ_FALSE;
 643     }
 644
 645     /* Bitdepth */
 646     for (i = 0; i < image->numcomps; i++) {
 647         if ((image->comps[i].bpp != 12) | (image->comps[i].sgnd)){
 648             char signed_str[] = "signed";
 649             char unsigned_str[] = "unsigned";
 650             char *tmp_str = image->comps[i].sgnd?signed_str:unsigned_str;
 651             fprintf(stdout,"JPEG 2000 Profile-3 (2k dc profile) requires:\n"
 652                     "Precision of each component shall be 12 bits unsigned"
 653                     "-> At least component %d of input image (%d bits, %s) is not compliant\n"
 654                     "-> Non-profile-3 codestream will be generated\n",
 655                     i,image->comps[i].bpp, tmp_str);
 656             return OPJ_FALSE;
 657         }
 658     }
 659
 660     /* Image size */
 661     switch (cinema_mode){
 662     case OPJ_CINEMA2K_24:
 663     case OPJ_CINEMA2K_48:
 664         if (!((image->comps[0].w > 2048) | (image->comps[0].h > 1080))){
 665             fprintf(stdout,"JPEG 2000 Profile-3 (2k dc profile) requires:\n"
 666                     "width <= 2048 and height <= 1080\n"
 667                     "-> Input image size %d x %d is not compliant\n"
 668                     "-> Non-profile-3 codestream will be generated\n",
 669                     image->comps[0].w,image->comps[0].h);
 670             return OPJ_FALSE;
 671         }
 672         break;
 673     case OPJ_CINEMA4K_24:
 674         if (!((image->comps[0].w > 4096) | (image->comps[0].h > 2160))){
 675             fprintf(stdout,"JPEG 2000 Profile-4 (4k dc profile) requires:\n"
 676                     "width <= 4096 and height <= 2160\n"
 677                     "-> Image size %d x %d is not compliant\n"
 678                     "-> Non-profile-4 codestream will be generated\n",
 679                     image->comps[0].w,image->comps[0].h);
 680             return OPJ_FALSE;
 681         }
 682         break;
 683     default :
 684         break;
 685     }
 686
 687     return OPJ_TRUE;
 688 }
 689
 690 /* ------------------------------------------------------------------------------------ */
 691
 692 static int parse_cmdline_encoder(int argc, char **argv, opj_cparameters_t *parameters,
 693                                  img_fol_t *img_fol, raw_cparameters_t *raw_cp, char *indexfilename) {
 694     OPJ_UINT32 i, j;
 695     int totlen, c;
 696     opj_option_t long_option[]={
 697         {"cinema2K",REQ_ARG, NULL ,'w'},
 698         {"cinema4K",NO_ARG, NULL ,'y'},
 699         {"ImgDir",REQ_ARG, NULL ,'z'},
 700         {"TP",REQ_ARG, NULL ,'u'},
 701         {"SOP",NO_ARG, NULL ,'S'},
 702         {"EPH",NO_ARG, NULL ,'E'},
 703         {"OutFor",REQ_ARG, NULL ,'O'},
 704         {"POC",REQ_ARG, NULL ,'P'},
 705         {"ROI",REQ_ARG, NULL ,'R'},
 706         {"jpip",NO_ARG, NULL, 'J'}
 707     };
 708
 709     /* parse the command line */
 710     const char optlist[] = "i:o:r:q:n:b:c:t:p:s:SEM:x:R:d:T:If:P:C:F:u:J"
 711         #ifdef USE_JPWL
 712             "W:"
 713         #endif /* USE_JPWL */
 714             "h";
 715
 716     totlen=sizeof(long_option);
 717     img_fol->set_out_format=0;
 718     raw_cp->rawWidth = 0;
 719
 720     do{
 721         c = opj_getopt_long(argc, argv, optlist,long_option,totlen);
 722         if (c == -1)
 723             break;
 724         switch (c) {
 725         case 'i':                       /* input file */
 726         {
 727             char *infile = opj_optarg;
 728             parameters->decod_format = get_file_format(infile);
 729             switch(parameters->decod_format) {
 730             case PGX_DFMT:
 731             case PXM_DFMT:
 732             case BMP_DFMT:
 733             case TIF_DFMT:
 734             case RAW_DFMT:
 735             case RAWL_DFMT:
 736             case TGA_DFMT:
 737             case PNG_DFMT:
 738                 break;
 739             default:
 740                 fprintf(stderr,
 741                         "!! Unrecognized format for infile : %s "
 742                         "[accept only *.pnm, *.pgm, *.ppm, *.pgx, *png, *.bmp, *.tif, *.raw or *.tga] !!\n\n",
 743                         infile);
 744                 return 1;
 745             }
 746             strncpy(parameters->infile, infile, sizeof(parameters->infile)-1);
 747         }
 748             break;
 749
 750             /* ----------------------------------------------------- */
 751
 752         case 'o':                       /* output file */
 753         {
 754             char *outfile = opj_optarg;
 755             parameters->cod_format = get_file_format(outfile);
 756             switch(parameters->cod_format) {
 757             case J2K_CFMT:
 758             case JP2_CFMT:
 759                 break;
 760             default:
 761                 fprintf(stderr, "Unknown output format image %s [only *.j2k, *.j2c or *.jp2]!! \n", outfile);
 762                 return 1;
 763             }
 764             strncpy(parameters->outfile, outfile, sizeof(parameters->outfile)-1);
 765         }
 766             break;
 767
 768             /* ----------------------------------------------------- */
 769         case 'O':                       /* output format */
 770         {
 771             char outformat[50];
 772             char *of = opj_optarg;
 773             sprintf(outformat,".%s",of);
 774             img_fol->set_out_format = 1;
 775             parameters->cod_format = get_file_format(outformat);
 776             switch(parameters->cod_format) {
 777             case J2K_CFMT:
 778             case JP2_CFMT:
 779                 img_fol->out_format = opj_optarg;
 780                 break;
 781             default:
 782                 fprintf(stderr, "Unknown output format image [only j2k, j2c, jp2]!! \n");
 783                 return 1;
 784             }
 785         }
 786             break;
 787
 788
 789             /* ----------------------------------------------------- */
 790
 791
 792         case 'r':                       /* rates rates/distorsion */
 793         {
 794             char *s = opj_optarg;
 795             parameters->tcp_numlayers = 0;
 796             while (sscanf(s, "%f", &parameters->tcp_rates[parameters->tcp_numlayers]) == 1) {
 797                 parameters->tcp_numlayers++;
 798                 while (*s && *s != ',') {
 799                     s++;
 800                 }
 801                 if (!*s)
 802                     break;
 803                 s++;
 804             }
 805             parameters->cp_disto_alloc = 1;
 806         }
 807             break;
 808
 809             /* ----------------------------------------------------- */
 810
 811
 812         case 'F':                       /* Raw image format parameters */
 813         {
 814             char signo;
 815             char *s = opj_optarg;
 816             if (sscanf(s, "%d,%d,%d,%d,%c", &raw_cp->rawWidth, &raw_cp->rawHeight, &raw_cp->rawComp, &raw_cp->rawBitDepth, &signo) == 5) {
 817                 if (signo == 's') {
 818                     raw_cp->rawSigned = OPJ_TRUE;
 819                     fprintf(stdout,"\nRaw file parameters: %d,%d,%d,%d Signed\n", raw_cp->rawWidth, raw_cp->rawHeight, raw_cp->rawComp, raw_cp->rawBitDepth);
 820                 }
 821                 else if (signo == 'u') {
 822                     raw_cp->rawSigned = OPJ_FALSE;
 823                     fprintf(stdout,"\nRaw file parameters: %d,%d,%d,%d Unsigned\n", raw_cp->rawWidth, raw_cp->rawHeight, raw_cp->rawComp, raw_cp->rawBitDepth);
 824                 }
 825                 else {
 826                     fprintf(stderr,"\nError: invalid raw image parameters: Unknown sign of raw file\n");
 827                     fprintf(stderr,"Please use the Format option -F:\n");
 828                     fprintf(stderr,"-F rawWidth,rawHeight,rawComp,rawBitDepth,s/u (Signed/Unsigned)\n");
 829                     fprintf(stderr,"Example: -i lena.raw -o lena.j2k -F 512,512,3,8,u\n");
 830                     fprintf(stderr,"Aborting\n");
 831                 }
 832             }
 833             else {
 834                 fprintf(stderr,"\nError: invalid raw image parameters\n");
 835                 fprintf(stderr,"Please use the Format option -F:\n");
 836                 fprintf(stderr,"-F rawWidth,rawHeight,rawComp,rawBitDepth,s/u (Signed/Unsigned)\n");
 837                 fprintf(stderr,"Example: -i lena.raw -o lena.j2k -F 512,512,3,8,u\n");
 838                 fprintf(stderr,"Aborting\n");
 839                 return 1;
 840             }
 841         }
 842             break;
 843
 844             /* ----------------------------------------------------- */
 845
 846         case 'q':                       /* add fixed_quality */
 847         {
 848             char *s = opj_optarg;
 849             while (sscanf(s, "%f", &parameters->tcp_distoratio[parameters->tcp_numlayers]) == 1) {
 850                 parameters->tcp_numlayers++;
 851                 while (*s && *s != ',') {
 852                     s++;
 853                 }
 854                 if (!*s)
 855                     break;
 856                 s++;
 857             }
 858             parameters->cp_fixed_quality = 1;
 859         }
 860             break;
 861
 862             /* dda */
 863             /* ----------------------------------------------------- */
 864
 865         case 'f':                       /* mod fixed_quality (before : -q) */
 866         {
 867             int *row = NULL, *col = NULL;
 868             OPJ_UINT32 numlayers = 0, numresolution = 0, matrix_width = 0;
 869
 870             char *s = opj_optarg;
 871             sscanf(s, "%d", &numlayers);
 872             s++;
 873             if (numlayers > 9)
 874                 s++;
 875
 876             parameters->tcp_numlayers = numlayers;
 877             numresolution = parameters->numresolution;
 878             matrix_width = numresolution * 3;
 879             parameters->cp_matrice = (int *) malloc(numlayers * matrix_width * sizeof(int));
 880             s = s + 2;
 881
 882             for (i = 0; i < numlayers; i++) {
 883                 row = &parameters->cp_matrice[i * matrix_width];
 884                 col = row;
 885                 parameters->tcp_rates[i] = 1;
 886                 sscanf(s, "%d,", &col[0]);
 887                 s += 2;
 888                 if (col[0] > 9)
 889                     s++;
 890                 col[1] = 0;
 891                 col[2] = 0;
 892                 for (j = 1; j < numresolution; j++) {
 893                     col += 3;
 894                     sscanf(s, "%d,%d,%d", &col[0], &col[1], &col[2]);
 895                     s += 6;
 896                     if (col[0] > 9)
 897                         s++;
 898                     if (col[1] > 9)
 899                         s++;
 900                     if (col[2] > 9)
 901                         s++;
 902                 }
 903                 if (i < numlayers - 1)
 904                     s++;
 905             }
 906             parameters->cp_fixed_alloc = 1;
 907         }
 908             break;
 909
 910             /* ----------------------------------------------------- */
 911
 912         case 't':                       /* tiles */
 913         {
 914             sscanf(opj_optarg, "%d,%d", &parameters->cp_tdx, &parameters->cp_tdy);
 915             parameters->tile_size_on = OPJ_TRUE;
 916         }
 917             break;
 918
 919             /* ----------------------------------------------------- */
 920
 921         case 'n':                       /* resolution */
 922         {
 923             sscanf(opj_optarg, "%d", &parameters->numresolution);
 924         }
 925             break;
 926
 927             /* ----------------------------------------------------- */
 928         case 'c':                       /* precinct dimension */
 929         {
 930             char sep;
 931             int res_spec = 0;
 932
 933             char *s = opj_optarg;
 934             do {
 935                 sep = 0;
 936                 sscanf(s, "[%d,%d]%c", &parameters->prcw_init[res_spec],
 937                        &parameters->prch_init[res_spec], &sep);
 938                 parameters->csty |= 0x01;
 939                 res_spec++;
 940                 s = strpbrk(s, "]") + 2;
 941             }
 942             while (sep == ',');
 943             parameters->res_spec = res_spec;
 944         }
 945             break;
 946
 947             /* ----------------------------------------------------- */
 948
 949         case 'b':                       /* code-block dimension */
 950         {
 951             int cblockw_init = 0, cblockh_init = 0;
 952             sscanf(opj_optarg, "%d,%d", &cblockw_init, &cblockh_init);
 953             if (cblockw_init * cblockh_init > 4096 || cblockw_init > 1024
 954                     || cblockw_init < 4 || cblockh_init > 1024 || cblockh_init < 4) {
 955                 fprintf(stderr,
 956                         "!! Size of code_block error (option -b) !!\n\nRestriction :\n"
 957                         "    * width*height<=4096\n    * 4<=width,height<= 1024\n\n");
 958                 return 1;
 959             }
 960             parameters->cblockw_init = cblockw_init;
 961             parameters->cblockh_init = cblockh_init;
 962         }
 963             break;
 964
 965             /* ----------------------------------------------------- */
 966
 967         case 'x':                       /* creation of index file */
 968         {
 969             char *index = opj_optarg;
 970             strncpy(indexfilename, index, OPJ_PATH_LEN);
 971         }
 972             break;
 973
 974             /* ----------------------------------------------------- */
 975
 976         case 'p':                       /* progression order */
 977         {
 978             char progression[4];
 979
 980             strncpy(progression, opj_optarg, 4);
 981             parameters->prog_order = give_progression(progression);
 982             if (parameters->prog_order == -1) {
 983                 fprintf(stderr, "Unrecognized progression order "
 984                         "[LRCP, RLCP, RPCL, PCRL, CPRL] !!\n");
 985                 return 1;
 986             }
 987         }
 988             break;
 989
 990             /* ----------------------------------------------------- */
 991
 992         case 's':                       /* subsampling factor */
 993         {
 994             if (sscanf(opj_optarg, "%d,%d", &parameters->subsampling_dx,
 995                        &parameters->subsampling_dy) != 2) {
 996                 fprintf(stderr, "'-s' sub-sampling argument error !  [-s dx,dy]\n");
 997                 return 1;
 998             }
 999         }
1000             break;
1001
1002             /* ----------------------------------------------------- */
1003
1004         case 'd':                       /* coordonnate of the reference grid */
1005         {
1006             if (sscanf(opj_optarg, "%d,%d", &parameters->image_offset_x0,
1007                        &parameters->image_offset_y0) != 2) {
1008                 fprintf(stderr, "-d 'coordonnate of the reference grid' argument "
1009                         "error !! [-d x0,y0]\n");
1010                 return 1;
1011             }
1012         }
1013             break;
1014
1015             /* ----------------------------------------------------- */
1016
1017         case 'h':                       /* display an help description */
1018             encode_help_display();
1019             return 1;
1020
1021             /* ----------------------------------------------------- */
1022
1023         case 'P':                       /* POC */
1024         {
1025             int numpocs = 0;            /* number of progression order change (POC) default 0 */
1026             opj_poc_t *POC = NULL;      /* POC : used in case of Progression order change */
1027
1028             char *s = opj_optarg;
1029             POC = parameters->POC;
1030
1031             while (sscanf(s, "T%d=%d,%d,%d,%d,%d,%4s", &POC[numpocs].tile,
1032                           &POC[numpocs].resno0, &POC[numpocs].compno0,
1033                           &POC[numpocs].layno1, &POC[numpocs].resno1,
1034                           &POC[numpocs].compno1, POC[numpocs].progorder) == 7) {
1035                 POC[numpocs].prg1 = give_progression(POC[numpocs].progorder);
1036                 numpocs++;
1037                 while (*s && *s != '/') {
1038                     s++;
1039                 }
1040                 if (!*s) {
1041                     break;
1042                 }
1043                 s++;
1044             }
1045             parameters->numpocs = numpocs;
1046         }
1047             break;
1048
1049             /* ------------------------------------------------------ */
1050
1051         case 'S':                       /* SOP marker */
1052         {
1053             parameters->csty |= 0x02;
1054         }
1055             break;
1056
1057             /* ------------------------------------------------------ */
1058
1059         case 'E':                       /* EPH marker */
1060         {
1061             parameters->csty |= 0x04;
1062         }
1063             break;
1064
1065             /* ------------------------------------------------------ */
1066
1067         case 'M':                       /* Mode switch pas tous au point !! */
1068         {
1069             int value = 0;
1070             if (sscanf(opj_optarg, "%d", &value) == 1) {
1071                 for (i = 0; i <= 5; i++) {
1072                     int cache = value & (1 << i);
1073                     if (cache)
1074                         parameters->mode |= (1 << i);
1075                 }
1076             }
1077         }
1078             break;
1079
1080             /* ------------------------------------------------------ */
1081
1082         case 'R':                       /* ROI */
1083         {
1084             if (sscanf(opj_optarg, "c=%d,U=%d", &parameters->roi_compno,
1085                        &parameters->roi_shift) != 2) {
1086                 fprintf(stderr, "ROI error !! [-ROI c='compno',U='shift']\n");
1087                 return 1;
1088             }
1089         }
1090             break;
1091
1092             /* ------------------------------------------------------ */
1093
1094         case 'T':                       /* Tile offset */
1095         {
1096             if (sscanf(opj_optarg, "%d,%d", &parameters->cp_tx0, &parameters->cp_ty0) != 2) {
1097                 fprintf(stderr, "-T 'tile offset' argument error !! [-T X0,Y0]");
1098                 return 1;
1099             }
1100         }
1101             break;
1102
1103             /* ------------------------------------------------------ */
1104
1105         case 'C':                       /* add a comment */
1106         {
1107             parameters->cp_comment = (char*)malloc(strlen(opj_optarg) + 1);
1108             if(parameters->cp_comment) {
1109                 strcpy(parameters->cp_comment, opj_optarg);
1110             }
1111         }
1112             break;
1113
1114
1115             /* ------------------------------------------------------ */
1116
1117         case 'I':                       /* reversible or not */
1118         {
1119             parameters->irreversible = 1;
1120         }
1121             break;
1122
1123             /* ------------------------------------------------------ */
1124
1125         case 'u':                       /* Tile part generation*/
1126         {
1127             parameters->tp_flag = opj_optarg[0];
1128             parameters->tp_on = 1;
1129         }
1130             break;
1131
1132             /* ------------------------------------------------------ */
1133
1134         case 'z':                       /* Image Directory path */
1135         {
1136             img_fol->imgdirpath = (char*)malloc(strlen(opj_optarg) + 1);
1137             strcpy(img_fol->imgdirpath,opj_optarg);
1138             img_fol->set_imgdir=1;
1139         }
1140             break;
1141
1142             /* ------------------------------------------------------ */
1143
1144         case 'w':                       /* Digital Cinema 2K profile compliance*/
1145         {
1146             int fps=0;
1147             sscanf(opj_optarg,"%d",&fps);
1148             if(fps == 24){
1149                 parameters->cp_cinema = OPJ_CINEMA2K_24;
1150             }else if(fps == 48 ){
1151                 parameters->cp_cinema = OPJ_CINEMA2K_48;
1152             }else {
1153                 fprintf(stderr,"Incorrect value!! must be 24 or 48\n");
1154                 return 1;
1155             }
1156             fprintf(stdout,"CINEMA 2K profile activated\n"
1157                     "Other options specified could be overriden\n");
1158
1159         }
1160             break;
1161
1162             /* ------------------------------------------------------ */
1163
1164         case 'y':                       /* Digital Cinema 4K profile compliance*/
1165         {
1166             parameters->cp_cinema = OPJ_CINEMA4K_24;
1167             fprintf(stdout,"CINEMA 4K profile activated\n"
1168                     "Other options specified could be overriden\n");
1169         }
1170             break;
1171
1172             /* ------------------------------------------------------ */
1173         case 'm':                       /* mct input file */
1174         {
1175             char *lFilename = opj_optarg;
1176             char *lMatrix;
1177             char *lCurrentPtr ;
1178             float *lCurrentDoublePtr;
1179             float *lSpace;
1180             int *l_int_ptr;
1181             int lNbComp = 0, lTotalComp, lMctComp, i, lStrLen;
1182
1183             /* Open file */
1184             FILE * lFile = fopen(lFilename,"r");
1185             if (lFile == NULL) {
1186                 return 1;
1187             }
1188
1189             /* Set size of file and read its content*/
1190             fseek(lFile,0,SEEK_END);
1191             lStrLen = ftell(lFile);
1192             fseek(lFile,0,SEEK_SET);
1193             lMatrix = (char *) malloc(lStrLen + 1);
1194             fread(lMatrix, lStrLen, 1, lFile);
1195             fclose(lFile);
1196
1197             lMatrix[lStrLen] = 0;
1198             lCurrentPtr = lMatrix;
1199
1200             /* replace ',' by 0 */
1201             while (*lCurrentPtr != 0 ) {
1202                 if (*lCurrentPtr == ' ') {
1203                     *lCurrentPtr = 0;
1204                     ++lNbComp;
1205                 }
1206                 ++lCurrentPtr;
1207             }
1208             ++lNbComp;
1209             lCurrentPtr = lMatrix;
1210
1211             lNbComp = (int) (sqrt(4*lNbComp + 1)/2. - 0.5);
1212             lMctComp = lNbComp * lNbComp;
1213             lTotalComp = lMctComp + lNbComp;
1214             lSpace = (float *) malloc(lTotalComp * sizeof(float));
1215             lCurrentDoublePtr = lSpace;
1216             for (i=0;i<lMctComp;++i) {
1217                 lStrLen = strlen(lCurrentPtr) + 1;
1218                 *lCurrentDoublePtr++ = (float) atof(lCurrentPtr);
1219                 lCurrentPtr += lStrLen;
1220             }
1221
1222             l_int_ptr = (int*) lCurrentDoublePtr;
1223             for (i=0;i<lNbComp;++i) {
1224                 lStrLen = strlen(lCurrentPtr) + 1;
1225                 *l_int_ptr++ = atoi(lCurrentPtr);
1226                 lCurrentPtr += lStrLen;
1227             }
1228
1229             /* TODO should not be here ! */
1230             opj_set_MCT(parameters, lSpace, (int *)(lSpace + lMctComp), lNbComp);
1231
1232             /* Free memory*/
1233             free(lSpace);
1234             free(lMatrix);
1235         }
1236             break;
1237
1238
1239             /* ------------------------------------------------------ */
1240
1241             /* UniPG>> */
1242 #ifdef USE_JPWL
1243             /* ------------------------------------------------------ */
1244
1245         case 'W':                       /* JPWL capabilities switched on */
1246         {
1247             char *token = NULL;
1248             int hprot, pprot, sens, addr, size, range;
1249
1250             /* we need to enable indexing */
1251             if (!indexfilename || !*indexfilename) {
1252                 strncpy(indexfilename, JPWL_PRIVATEINDEX_NAME, OPJ_PATH_LEN);
1253             }
1254
1255             /* search for different protection methods */
1256
1257             /* break the option in comma points and parse the result */
1258             token = strtok(opj_optarg, ",");
1259             while(token != NULL) {
1260
1261                 /* search header error protection method */
1262                 if (*token == 'h') {
1263
1264                     static int tile = 0, tilespec = 0, lasttileno = 0;
1265
1266                     hprot = 1; /* predefined method */
1267
1268                     if(sscanf(token, "h=%d", &hprot) == 1) {
1269                         /* Main header, specified */
1270                         if (!((hprot == 0) || (hprot == 1) || (hprot == 16) || (hprot == 32) ||
1271                               ((hprot >= 37) && (hprot <= 128)))) {
1272                             fprintf(stderr, "ERROR -> invalid main header protection method h = %d\n", hprot);
1273                             return 1;
1274                         }
1275                         parameters->jpwl_hprot_MH = hprot;
1276
1277                     } else if(sscanf(token, "h%d=%d", &tile, &hprot) == 2) {
1278                         /* Tile part header, specified */
1279                         if (!((hprot == 0) || (hprot == 1) || (hprot == 16) || (hprot == 32) ||
1280                               ((hprot >= 37) && (hprot <= 128)))) {
1281                             fprintf(stderr, "ERROR -> invalid tile part header protection method h = %d\n", hprot);
1282                             return 1;
1283                         }
1284                         if (tile < 0) {
1285                             fprintf(stderr, "ERROR -> invalid tile part number on protection method t = %d\n", tile);
1286                             return 1;
1287                         }
1288                         if (tilespec < JPWL_MAX_NO_TILESPECS) {
1289                             parameters->jpwl_hprot_TPH_tileno[tilespec] = lasttileno = tile;
1290                             parameters->jpwl_hprot_TPH[tilespec++] = hprot;
1291                         }
1292
1293                     } else if(sscanf(token, "h%d", &tile) == 1) {
1294                         /* Tile part header, unspecified */
1295                         if (tile < 0) {
1296                             fprintf(stderr, "ERROR -> invalid tile part number on protection method t = %d\n", tile);
1297                             return 1;
1298                         }
1299                         if (tilespec < JPWL_MAX_NO_TILESPECS) {
1300                             parameters->jpwl_hprot_TPH_tileno[tilespec] = lasttileno = tile;
1301                             parameters->jpwl_hprot_TPH[tilespec++] = hprot;
1302                         }
1303
1304
1305                     } else if (!strcmp(token, "h")) {
1306                         /* Main header, unspecified */
1307                         parameters->jpwl_hprot_MH = hprot;
1308
1309                     } else {
1310                         fprintf(stderr, "ERROR -> invalid protection method selection = %s\n", token);
1311                         return 1;
1312                     };
1313
1314                 }
1315
1316                 /* search packet error protection method */
1317                 if (*token == 'p') {
1318
1319                     static int pack = 0, tile = 0, packspec = 0;
1320
1321                     pprot = 1; /* predefined method */
1322
1323                     if (sscanf(token, "p=%d", &pprot) == 1) {
1324                         /* Method for all tiles and all packets */
1325                         if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) ||
1326                               ((pprot >= 37) && (pprot <= 128)))) {
1327                             fprintf(stderr, "ERROR -> invalid default packet protection method p = %d\n", pprot);
1328                             return 1;
1329                         }
1330                         parameters->jpwl_pprot_tileno[0] = 0;
1331                         parameters->jpwl_pprot_packno[0] = 0;
1332                         parameters->jpwl_pprot[0] = pprot;
1333
1334                     } else if (sscanf(token, "p%d=%d", &tile, &pprot) == 2) {
1335                         /* method specified from that tile on */
1336                         if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) ||
1337                               ((pprot >= 37) && (pprot <= 128)))) {
1338                             fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot);
1339                             return 1;
1340                         }
1341                         if (tile < 0) {
1342                             fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile);
1343                             return 1;
1344                         }
1345                         if (packspec < JPWL_MAX_NO_PACKSPECS) {
1346                             parameters->jpwl_pprot_tileno[packspec] = tile;
1347                             parameters->jpwl_pprot_packno[packspec] = 0;
1348                             parameters->jpwl_pprot[packspec++] = pprot;
1349                         }
1350
1351                     } else if (sscanf(token, "p%d:%d=%d", &tile, &pack, &pprot) == 3) {
1352                         /* method fully specified from that tile and that packet on */
1353                         if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) ||
1354                               ((pprot >= 37) && (pprot <= 128)))) {
1355                             fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot);
1356                             return 1;
1357                         }
1358                         if (tile < 0) {
1359                             fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile);
1360                             return 1;
1361                         }
1362                         if (pack < 0) {
1363                             fprintf(stderr, "ERROR -> invalid packet number on protection method p = %d\n", pack);
1364                             return 1;
1365                         }
1366                         if (packspec < JPWL_MAX_NO_PACKSPECS) {
1367                             parameters->jpwl_pprot_tileno[packspec] = tile;
1368                             parameters->jpwl_pprot_packno[packspec] = pack;
1369                             parameters->jpwl_pprot[packspec++] = pprot;
1370                         }
1371
1372                     } else if (sscanf(token, "p%d:%d", &tile, &pack) == 2) {
1373                         /* default method from that tile and that packet on */
1374                         if (!((pprot == 0) || (pprot == 1) || (pprot == 16) || (pprot == 32) ||
1375                               ((pprot >= 37) && (pprot <= 128)))) {
1376                             fprintf(stderr, "ERROR -> invalid packet protection method p = %d\n", pprot);
1377                             return 1;
1378                         }
1379                         if (tile < 0) {
1380                             fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile);
1381                             return 1;
1382                         }
1383                         if (pack < 0) {
1384                             fprintf(stderr, "ERROR -> invalid packet number on protection method p = %d\n", pack);
1385                             return 1;
1386                         }
1387                         if (packspec < JPWL_MAX_NO_PACKSPECS) {
1388                             parameters->jpwl_pprot_tileno[packspec] = tile;
1389                             parameters->jpwl_pprot_packno[packspec] = pack;
1390                             parameters->jpwl_pprot[packspec++] = pprot;
1391                         }
1392
1393                     } else if (sscanf(token, "p%d", &tile) == 1) {
1394                         /* default from a tile on */
1395                         if (tile < 0) {
1396                             fprintf(stderr, "ERROR -> invalid tile part number on protection method p = %d\n", tile);
1397                             return 1;
1398                         }
1399                         if (packspec < JPWL_MAX_NO_PACKSPECS) {
1400                             parameters->jpwl_pprot_tileno[packspec] = tile;
1401                             parameters->jpwl_pprot_packno[packspec] = 0;
1402                             parameters->jpwl_pprot[packspec++] = pprot;
1403                         }
1404
1405
1406                     } else if (!strcmp(token, "p")) {
1407                         /* all default */
1408                         parameters->jpwl_pprot_tileno[0] = 0;
1409                         parameters->jpwl_pprot_packno[0] = 0;
1410                         parameters->jpwl_pprot[0] = pprot;
1411
1412                     } else {
1413                         fprintf(stderr, "ERROR -> invalid protection method selection = %s\n", token);
1414                         return 1;
1415                     };
1416
1417                 }
1418
1419                 /* search sensitivity method */
1420                 if (*token == 's') {
1421
1422                     static int tile = 0, tilespec = 0, lasttileno = 0;
1423
1424                     sens = 0; /* predefined: relative error */
1425
1426                     if(sscanf(token, "s=%d", &sens) == 1) {
1427                         /* Main header, specified */
1428                         if ((sens < -1) || (sens > 7)) {
1429                             fprintf(stderr, "ERROR -> invalid main header sensitivity method s = %d\n", sens);
1430                             return 1;
1431                         }
1432                         parameters->jpwl_sens_MH = sens;
1433
1434                     } else if(sscanf(token, "s%d=%d", &tile, &sens) == 2) {
1435                         /* Tile part header, specified */
1436                         if ((sens < -1) || (sens > 7)) {
1437                             fprintf(stderr, "ERROR -> invalid tile part header sensitivity method s = %d\n", sens);
1438                             return 1;
1439                         }
1440                         if (tile < 0) {
1441                             fprintf(stderr, "ERROR -> invalid tile part number on sensitivity method t = %d\n", tile);
1442                             return 1;
1443                         }
1444                         if (tilespec < JPWL_MAX_NO_TILESPECS) {
1445                             parameters->jpwl_sens_TPH_tileno[tilespec] = lasttileno = tile;
1446                             parameters->jpwl_sens_TPH[tilespec++] = sens;
1447                         }
1448
1449                     } else if(sscanf(token, "s%d", &tile) == 1) {
1450                         /* Tile part header, unspecified */
1451                         if (tile < 0) {
1452                             fprintf(stderr, "ERROR -> invalid tile part number on sensitivity method t = %d\n", tile);
1453                             return 1;
1454                         }
1455                         if (tilespec < JPWL_MAX_NO_TILESPECS) {
1456                             parameters->jpwl_sens_TPH_tileno[tilespec] = lasttileno = tile;
1457                             parameters->jpwl_sens_TPH[tilespec++] = hprot;
1458                         }
1459
1460                     } else if (!strcmp(token, "s")) {
1461                         /* Main header, unspecified */
1462                         parameters->jpwl_sens_MH = sens;
1463
1464                     } else {
1465                         fprintf(stderr, "ERROR -> invalid sensitivity method selection = %s\n", token);
1466                         return 1;
1467                     };
1468
1469                     parameters->jpwl_sens_size = 2; /* 2 bytes for default size */
1470                 }
1471
1472                 /* search addressing size */
1473                 if (*token == 'a') {
1474
1475
1476                     addr = 0; /* predefined: auto */
1477
1478                     if(sscanf(token, "a=%d", &addr) == 1) {
1479                         /* Specified */
1480                         if ((addr != 0) && (addr != 2) && (addr != 4)) {
1481                             fprintf(stderr, "ERROR -> invalid addressing size a = %d\n", addr);
1482                             return 1;
1483                         }
1484                         parameters->jpwl_sens_addr = addr;
1485
1486                     } else if (!strcmp(token, "a")) {
1487                         /* default */
1488                         parameters->jpwl_sens_addr = addr; /* auto for default size */
1489
1490                     } else {
1491                         fprintf(stderr, "ERROR -> invalid addressing selection = %s\n", token);
1492                         return 1;
1493                     };
1494
1495                 }
1496
1497                 /* search sensitivity size */
1498                 if (*token == 'z') {
1499
1500
1501                     size = 1; /* predefined: 1 byte */
1502
1503                     if(sscanf(token, "z=%d", &size) == 1) {
1504                         /* Specified */
1505                         if ((size != 0) && (size != 1) && (size != 2)) {
1506                             fprintf(stderr, "ERROR -> invalid sensitivity size z = %d\n", size);
1507                             return 1;
1508                         }
1509                         parameters->jpwl_sens_size = size;
1510
1511                     } else if (!strcmp(token, "a")) {
1512                         /* default */
1513                         parameters->jpwl_sens_size = size; /* 1 for default size */
1514
1515                     } else {
1516                         fprintf(stderr, "ERROR -> invalid size selection = %s\n", token);
1517                         return 1;
1518                     };
1519
1520                 }
1521
1522                 /* search range method */
1523                 if (*token == 'g') {
1524
1525
1526                     range = 0; /* predefined: 0 (packet) */
1527
1528                     if(sscanf(token, "g=%d", &range) == 1) {
1529                         /* Specified */
1530                         if ((range < 0) || (range > 3)) {
1531                             fprintf(stderr, "ERROR -> invalid sensitivity range method g = %d\n", range);
1532                             return 1;
1533                         }
1534                         parameters->jpwl_sens_range = range;
1535
1536                     } else if (!strcmp(token, "g")) {
1537                         /* default */
1538                         parameters->jpwl_sens_range = range;
1539
1540                     } else {
1541                         fprintf(stderr, "ERROR -> invalid range selection = %s\n", token);
1542                         return 1;
1543                     };
1544
1545                 }
1546
1547                 /* next token or bust */
1548                 token = strtok(NULL, ",");
1549             };
1550
1551
1552             /* some info */
1553             fprintf(stdout, "Info: JPWL capabilities enabled\n");
1554             parameters->jpwl_epc_on = OPJ_TRUE;
1555
1556         }
1557             break;
1558 #endif /* USE_JPWL */
1559             /* <<UniPG */
1560             /* ------------------------------------------------------ */
1561
1562         case 'J':                       /* jpip on */
1563         {
1564             parameters->jpip_on = OPJ_TRUE;
1565         }
1566             break;
1567             /* ------------------------------------------------------ */
1568
1569
1570         default:
1571             fprintf(stderr, "ERROR -> Command line not valid\n");
1572             return 1;
1573         }
1574     }while(c != -1);
1575
1576     if(img_fol->set_imgdir == 1){
1577         if(!(parameters->infile[0] == 0)){
1578             fprintf(stderr, "Error: options -ImgDir and -i cannot be used together !!\n");
1579             return 1;
1580         }
1581         if(img_fol->set_out_format == 0){
1582             fprintf(stderr, "Error: When -ImgDir is used, -OutFor <FORMAT> must be used !!\n");
1583             fprintf(stderr, "Only one format allowed! Valid formats are j2k and jp2!!\n");
1584             return 1;
1585         }
1586         if(!((parameters->outfile[0] == 0))){
1587             fprintf(stderr, "Error: options -ImgDir and -o cannot be used together !!\n");
1588             fprintf(stderr, "Specify OutputFormat using -OutFor<FORMAT> !!\n");
1589             return 1;
1590         }
1591     }else{
1592         if((parameters->infile[0] == 0) || (parameters->outfile[0] == 0)) {
1593             fprintf(stderr, "Example: %s -i image.ppm  -o image.j2k\n",argv[0]);
1594             fprintf(stderr, "    Try: %s -h\n",argv[0]);
1595             return 1;
1596         }
1597     }
1598
1599     if ( (parameters->decod_format == RAW_DFMT && raw_cp->rawWidth == 0)
1600          || (parameters->decod_format == RAWL_DFMT && raw_cp->rawWidth == 0)) {
1601         fprintf(stderr,"\nError: invalid raw image parameters\n");
1602         fprintf(stderr,"Please use the Format option -F:\n");
1603         fprintf(stderr,"-F rawWidth,rawHeight,rawComp,rawBitDepth,s/u (Signed/Unsigned)\n");
1604         fprintf(stderr,"Example: -i lena.raw -o lena.j2k -F 512,512,3,8,u\n");
1605         fprintf(stderr,"Aborting\n");
1606         return 1;
1607     }
1608
1609     if ((parameters->cp_disto_alloc || parameters->cp_fixed_alloc || parameters->cp_fixed_quality)
1610             && (!(parameters->cp_disto_alloc ^ parameters->cp_fixed_alloc ^ parameters->cp_fixed_quality))) {
1611         fprintf(stderr, "Error: options -r -q and -f cannot be used together !!\n");
1612         return 1;
1613     }                           /* mod fixed_quality */
1614
1615     /* if no rate entered, lossless by default */
1616     if (parameters->tcp_numlayers == 0) {
1617         parameters->tcp_rates[0] = 0;   /* MOD antonin : losslessbug */
1618         parameters->tcp_numlayers++;
1619         parameters->cp_disto_alloc = 1;
1620     }
1621
1622     if((parameters->cp_tx0 > parameters->image_offset_x0) || (parameters->cp_ty0 > parameters->image_offset_y0)) {
1623         fprintf(stderr,
1624                 "Error: Tile offset dimension is unnappropriate --> TX0(%d)<=IMG_X0(%d) TYO(%d)<=IMG_Y0(%d) \n",
1625                 parameters->cp_tx0, parameters->image_offset_x0, parameters->cp_ty0, parameters->image_offset_y0);
1626         return 1;
1627     }
1628
1629     for (i = 0; i < parameters->numpocs; i++) {
1630         if (parameters->POC[i].prg == -1) {
1631             fprintf(stderr,
1632                     "Unrecognized progression order in option -P (POC n %d) [LRCP, RLCP, RPCL, PCRL, CPRL] !!\n",
1633                     i + 1);
1634         }
1635     }
1636
1637     return 0;
1638 }
1639
1640 /* -------------------------------------------------------------------------- */
1641
1642 /**
1643 sample error callback expecting a FILE* client object
1644 */
1645 static void error_file_callback(const char *msg, void *client_data) {
1646     FILE *stream = (FILE*)client_data;
1647     fprintf(stream, "[ERROR] %s", msg);
1648 }
1649 /**
1650 sample warning callback expecting a FILE* client object
1651 */
1652 static void warning_file_callback(const char *msg, void *client_data) {
1653     FILE *stream = (FILE*)client_data;
1654     fprintf(stream, "[WARNING] %s", msg);
1655 }
1656 /**
1657 sample debug callback expecting a FILE* client object
1658 */
1659 static void info_file_callback(const char *msg, void *client_data) {
1660     FILE *stream = (FILE*)client_data;
1661     fprintf(stream, "[INFO] %s", msg);
1662 }
1663
1664 /**
1665 sample error debug callback expecting no client object
1666 */
1667 static void error_callback(const char *msg, void *client_data) {
1668     (void)client_data;
1669     fprintf(stdout, "[ERROR] %s", msg);
1670 }
1671 /**
1672 sample warning debug callback expecting no client object
1673 */
1674 static void warning_callback(const char *msg, void *client_data) {
1675     (void)client_data;
1676     fprintf(stdout, "[WARNING] %s", msg);
1677 }
1678 /**
1679 sample debug callback expecting no client object
1680 */
1681 static void info_callback(const char *msg, void *client_data) {
1682     (void)client_data;
1683     fprintf(stdout, "[INFO] %s", msg);
1684 }
1685
1686 /* -------------------------------------------------------------------------- */
1687 /**
1688  * OPJ_COMPRESS MAIN
1689  */
1690 /* -------------------------------------------------------------------------- */
1691 int main(int argc, char **argv) {
1692
1693     opj_cparameters_t parameters;       /* compression parameters */
1694
1695     opj_stream_t *l_stream = 00;
1696     opj_codec_t* l_codec = 00;
1697     opj_image_t *image = NULL;
1698     raw_cparameters_t raw_cp;
1699
1700     char indexfilename[OPJ_PATH_LEN];   /* index file name */
1701
1702     unsigned int i, num_images, imageno;
1703     img_fol_t img_fol;
1704     dircnt_t *dirptr = NULL;
1705
1706     OPJ_BOOL bSuccess;
1707     OPJ_BOOL bUseTiles = OPJ_FALSE; /* OPJ_TRUE */
1708     OPJ_UINT32 l_nb_tiles = 4;
1709
1710     /* set encoding parameters to default values */
1711     opj_set_default_encoder_parameters(&parameters);
1712
1713     /* Initialize indexfilename and img_fol */
1714     *indexfilename = 0;
1715     memset(&img_fol,0,sizeof(img_fol_t));
1716
1717     /* parse input and get user encoding parameters */
1718     if(parse_cmdline_encoder(argc, argv, &parameters,&img_fol, &raw_cp, indexfilename) == 1) {
1719         return 1;
1720     }
1721
1722     if (parameters.cp_cinema){
1723         set_cinema_parameters(&parameters);
1724     }
1725
1726     /* Create comment for codestream */
1727     if(parameters.cp_comment == NULL) {
1728         const char comment[] = "Created by OpenJPEG version ";
1729         const size_t clen = strlen(comment);
1730         const char *version = opj_version();
1731         /* UniPG>> */
1732 #ifdef USE_JPWL
1733         parameters.cp_comment = (char*)malloc(clen+strlen(version)+11);
1734         sprintf(parameters.cp_comment,"%s%s with JPWL", comment, version);
1735 #else
1736         parameters.cp_comment = (char*)malloc(clen+strlen(version)+1);
1737         sprintf(parameters.cp_comment,"%s%s", comment, version);
1738 #endif
1739         /* <<UniPG */
1740     }
1741
1742     /* Read directory if necessary */
1743     if(img_fol.set_imgdir==1){
1744         num_images=get_num_images(img_fol.imgdirpath);
1745         dirptr=(dircnt_t*)malloc(sizeof(dircnt_t));
1746         if(dirptr){
1747             dirptr->filename_buf = (char*)malloc(num_images*OPJ_PATH_LEN*sizeof(char)); /* Stores at max 10 image file names*/
1748             dirptr->filename = (char**) malloc(num_images*sizeof(char*));
1749             if(!dirptr->filename_buf){
1750                 return 0;
1751             }
1752             for(i=0;i<num_images;i++){
1753                 dirptr->filename[i] = dirptr->filename_buf + i*OPJ_PATH_LEN;
1754             }
1755         }
1756         if(load_images(dirptr,img_fol.imgdirpath)==1){
1757             return 0;
1758         }
1759         if (num_images==0){
1760             fprintf(stdout,"Folder is empty\n");
1761             return 0;
1762         }
1763     }else{
1764         num_images=1;
1765     }
1766     /*Encoding image one by one*/
1767     for(imageno=0;imageno<num_images;imageno++) {
1768         image = NULL;
1769         fprintf(stderr,"\n");
1770
1771         if(img_fol.set_imgdir==1){
1772             if (get_next_file(imageno, dirptr,&img_fol, &parameters)) {
1773                 fprintf(stderr,"skipping file...\n");
1774                 continue;
1775             }
1776         }
1777
1778         switch(parameters.decod_format) {
1779         case PGX_DFMT:
1780             break;
1781         case PXM_DFMT:
1782             break;
1783         case BMP_DFMT:
1784             break;
1785         case TIF_DFMT:
1786             break;
1787         case RAW_DFMT:
1788         case RAWL_DFMT:
1789             break;
1790         case TGA_DFMT:
1791             break;
1792         case PNG_DFMT:
1793             break;
1794         default:
1795             fprintf(stderr,"skipping file...\n");
1796             continue;
1797         }
1798
1799         /* decode the source image */
1800         /* ----------------------- */
1801
1802         switch (parameters.decod_format) {
1803         case PGX_DFMT:
1804             image = pgxtoimage(parameters.infile, &parameters);
1805             if (!image) {
1806                 fprintf(stderr, "Unable to load pgx file\n");
1807                 return 1;
1808             }
1809             break;
1810
1811         case PXM_DFMT:
1812             image = pnmtoimage(parameters.infile, &parameters);
1813             if (!image) {
1814                 fprintf(stderr, "Unable to load pnm file\n");
1815                 return 1;
1816             }
1817             break;
1818
1819         case BMP_DFMT:
1820             image = bmptoimage(parameters.infile, &parameters);
1821             if (!image) {
1822                 fprintf(stderr, "Unable to load bmp file\n");
1823                 return 1;
1824             }
1825             break;
1826
1827 #ifdef OPJ_HAVE_LIBTIFF
1828         case TIF_DFMT:
1829             image = tiftoimage(parameters.infile, &parameters);
1830             if (!image) {
1831                 fprintf(stderr, "Unable to load tiff file\n");
1832                 return 1;
1833             }
1834             break;
1835 #endif /* OPJ_HAVE_LIBTIFF */
1836
1837         case RAW_DFMT:
1838             image = rawtoimage(parameters.infile, &parameters, &raw_cp);
1839             if (!image) {
1840                 fprintf(stderr, "Unable to load raw file\n");
1841                 return 1;
1842             }
1843             break;
1844
1845         case RAWL_DFMT:
1846             image = rawltoimage(parameters.infile, &parameters, &raw_cp);
1847             if (!image) {
1848                 fprintf(stderr, "Unable to load raw file\n");
1849                 return 1;
1850             }
1851             break;
1852
1853         case TGA_DFMT:
1854             image = tgatoimage(parameters.infile, &parameters);
1855             if (!image) {
1856                 fprintf(stderr, "Unable to load tga file\n");
1857                 return 1;
1858             }
1859             break;
1860
1861 #ifdef OPJ_HAVE_LIBPNG
1862         case PNG_DFMT:
1863             image = pngtoimage(parameters.infile, &parameters);
1864             if (!image) {
1865                 fprintf(stderr, "Unable to load png file\n");
1866                 return 1;
1867             }
1868             break;
1869 #endif /* OPJ_HAVE_LIBPNG */
1870         }
1871
1872         /* Can happen if input file is TIFF or PNG
1873  * and OPJ_HAVE_LIBTIF or OPJ_HAVE_LIBPNG is undefined
1874 */
1875         if( !image) {
1876             fprintf(stderr, "Unable to load file: got no image\n");
1877             return 1;
1878         }
1879
1880         /* Decide if MCT should be used */
1881         parameters.tcp_mct = image->numcomps == 3 ? 1 : 0;
1882
1883         if(parameters.cp_cinema){
1884             if (!is_cinema_compliant(image,parameters.cp_cinema)) {
1885                 parameters.cp_rsiz = OPJ_STD_RSIZ;
1886             }
1887             /* This cannot be in set_cinema_parameters because we need size and precision of the input image */
1888             set_cinema_rate(&parameters,image);
1889         }
1890
1891         /* encode the destination image */
1892         /* ---------------------------- */
1893
1894         switch(parameters.cod_format) {
1895         case J2K_CFMT:  /* JPEG-2000 codestream */
1896         {
1897             /* Get a decoder handle */
1898             l_codec = opj_create_compress(OPJ_CODEC_J2K);
1899             break;
1900         }
1901         case JP2_CFMT:  /* JPEG 2000 compressed image data */
1902         {
1903             /* Get a decoder handle */
1904             l_codec = opj_create_compress(OPJ_CODEC_JP2);
1905             break;
1906         }
1907         default:
1908             fprintf(stderr, "skipping file..\n");
1909             opj_stream_destroy_v3(l_stream);
1910             continue;
1911         }
1912
1913         /* catch events using our callbacks and give a local context */
1914         opj_set_info_handler(l_codec, info_callback,00);
1915         opj_set_warning_handler(l_codec, warning_callback,00);
1916         opj_set_error_handler(l_codec, error_callback,00);
1917
1918         if( bUseTiles ) {
1919             parameters.cp_tx0 = 0;
1920             parameters.cp_ty0 = 0;
1921             parameters.tile_size_on = OPJ_TRUE;
1922             parameters.cp_tdx = 512;
1923             parameters.cp_tdy = 512;
1924         }
1925         opj_setup_encoder(l_codec, &parameters, image);
1926
1927         /* open a byte stream for writing and allocate memory for all tiles */
1928         l_stream = opj_stream_create_default_file_stream_v3(parameters.outfile,OPJ_FALSE);
1929         if (! l_stream){
1930             return 1;
1931         }
1932
1933         /* encode the image */
1934         bSuccess = opj_start_compress(l_codec,image,l_stream);
1935         if (!bSuccess)  {
1936             fprintf(stderr, "failed to encode image: opj_start_compress\n");
1937         }
1938         if( bUseTiles ) {
1939             OPJ_BYTE *l_data;
1940             OPJ_UINT32 l_data_size = 512*512*3;
1941             l_data = (OPJ_BYTE*) malloc( l_data_size * sizeof(OPJ_BYTE));
1942             memset(l_data, 0, l_data_size );
1943             assert( l_data );
1944             for (i=0;i<l_nb_tiles;++i) {
1945                 if (! opj_write_tile(l_codec,i,l_data,l_data_size,l_stream)) {
1946                     fprintf(stderr, "ERROR -> test_tile_encoder: failed to write the tile %d!\n",i);
1947                     opj_stream_destroy_v3(l_stream);
1948                     opj_destroy_codec(l_codec);
1949                     opj_image_destroy(image);
1950                     return 1;
1951                 }
1952             }
1953             free(l_data);
1954         }
1955         else {
1956             bSuccess = bSuccess && opj_encode(l_codec, l_stream);
1957             if (!bSuccess)  {
1958                 fprintf(stderr, "failed to encode image: opj_encode\n");
1959             }
1960         }
1961         bSuccess = bSuccess && opj_end_compress(l_codec, l_stream);
1962         if (!bSuccess)  {
1963             fprintf(stderr, "failed to encode image: opj_end_compress\n");
1964         }
1965
1966         if (!bSuccess)  {
1967             opj_stream_destroy_v3(l_stream);
1968             opj_destroy_codec(l_codec);
1969             opj_image_destroy(image);
1970             fprintf(stderr, "failed to encode image\n");
1971             return 1;
1972         }
1973
1974         fprintf(stderr,"Generated outfile %s\n",parameters.outfile);
1975         /* close and free the byte stream */
1976         opj_stream_destroy_v3(l_stream);
1977
1978         /* free remaining compression structures */
1979         opj_destroy_codec(l_codec);
1980
1981         /* free image data */
1982         opj_image_destroy(image);
1983
1984     }
1985
1986     /* free user parameters structure */
1987     if(parameters.cp_comment)   free(parameters.cp_comment);
1988     if(parameters.cp_matrice)   free(parameters.cp_matrice);
1989
1990     return 0;
1991 }