Avoided ABI breakage
[openjpeg.git] / codec / image_to_j2k.c
index cbb9d0bc62684fec4b1e0e8543c870d32e6d736d..54ca14a1a5c47bc5bd42a691b662d4cf60653e61 100644 (file)
@@ -80,6 +80,8 @@ typedef struct img_folder{
        char set_imgdir;
        /** Enable Cod Format for output*/
        char set_out_format;
+       /** User specified rate stored in case of cinema option*/
+       float *rates;
 }img_fol_t;
 
 void encode_help_display() {
@@ -197,7 +199,7 @@ void encode_help_display() {
        fprintf(stdout,"\n");
        fprintf(stdout,"-ROI         : c=%%d,U=%%d : quantization indices upshifted \n");
        fprintf(stdout,"               for component c=%%d [%%d = 0,1,2]\n");
-       fprintf(stdout,"               with a value of U=%%d [0 <= %%d <= 37] (i.e. -ROI:c=0,U=25) \n");
+       fprintf(stdout,"               with a value of U=%%d [0 <= %%d <= 37] (i.e. -ROI c=0,U=25) \n");
        fprintf(stdout,"\n");
        fprintf(stdout,"-d           : offset of the origin of the image (-d 150,300) \n");
        fprintf(stdout,"\n");
@@ -287,12 +289,22 @@ void encode_help_display() {
        fprintf(stdout,"decomposition_levels\n");
        fprintf(stdout,"[Precincts_size_X_res_Nr Precincts_size_Y_res_Nr]...\n");
        fprintf(stdout,"   [Precincts_size_X_res_0 Precincts_size_Y_res_0]\n");
+       fprintf(stdout,"Main_header_start_position\n");
        fprintf(stdout,"Main_header_end_position\n");
        fprintf(stdout,"Codestream_size\n");
-       fprintf(stdout,"Tile_0 start_pos end_Theader end_pos TotalDisto NumPix MaxMSE\n");
-       fprintf(stdout,"Tile_1   ''           ''        ''        ''       ''    ''\n");
+       fprintf(stdout,"\n");
+       fprintf(stdout,"INFO ON TILES\n");
+       fprintf(stdout,"tileno start_pos end_hd end_tile nbparts disto nbpix disto/nbpix\n");
+       fprintf(stdout,"Tile_0 start_pos end_Theader end_pos NumParts TotalDisto NumPix MaxMSE\n");
+       fprintf(stdout,"Tile_1   ''           ''        ''        ''       ''    ''      ''\n");
+       fprintf(stdout,"...\n");
+       fprintf(stdout,"Tile_Nt   ''           ''        ''        ''       ''    ''     ''\n");
        fprintf(stdout,"...\n");
-       fprintf(stdout,"Tile_Nt   ''           ''        ''        ''       ''    ''\n");
+       fprintf(stdout,"TILE 0 DETAILS\n");
+       fprintf(stdout,"part_nb tileno num_packs start_pos end_tph_pos end_pos\n");
+       fprintf(stdout,"...\n");
+       fprintf(stdout,"Progression_string\n");
+       fprintf(stdout,"pack_nb tileno layno resno compno precno start_pos end_ph_pos end_pos disto\n");
        fprintf(stdout,"Tpacket_0 Tile layer res. comp. prec. start_pos end_pos disto\n");
        fprintf(stdout,"...\n");
        fprintf(stdout,"Tpacket_Np ''   ''    ''   ''    ''       ''       ''     ''\n");
@@ -472,7 +484,7 @@ void cinema_parameters(opj_cparameters_t *parameters){
 
 }
 
-void cinema_setup_encoder(opj_cparameters_t *parameters,opj_image_t *image){
+void cinema_setup_encoder(opj_cparameters_t *parameters,opj_image_t *image, img_fol_t *img_fol){
        int i;
        float temp_rate;
        opj_poc_t *POC = NULL;
@@ -510,17 +522,19 @@ void cinema_setup_encoder(opj_cparameters_t *parameters,opj_image_t *image){
        switch (parameters->cp_cinema){
                case CINEMA2K_24:
                case CINEMA4K_24:
-                       for(i=0;i<parameters->tcp_numlayers;i++){
+                       for(i=0 ; i<parameters->tcp_numlayers ; i++){
                                temp_rate = 0 ;
-                               if (parameters->tcp_rates[i]== 0){
+                               if (img_fol->rates[i]== 0){
                                        parameters->tcp_rates[0]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ 
                                        (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy);
                                }else{
                                        temp_rate =((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ 
-                                               (parameters->tcp_rates[i] * 8 * image->comps[0].dx * image->comps[0].dy);
+                                               (img_fol->rates[i] * 8 * image->comps[0].dx * image->comps[0].dy);
                                        if (temp_rate > CINEMA_24_CS ){
                                                parameters->tcp_rates[i]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ 
-                                                                                                                                                               (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy);
+                                               (CINEMA_24_CS * 8 * image->comps[0].dx * image->comps[0].dy);
+                                       }else{
+                                               parameters->tcp_rates[i]= img_fol->rates[i];
                                        }
                                }
                        }
@@ -528,17 +542,19 @@ void cinema_setup_encoder(opj_cparameters_t *parameters,opj_image_t *image){
                        break;
                
                case CINEMA2K_48:
-       for(i=0;i<parameters->tcp_numlayers;i++){
+                       for(i=0 ; i<parameters->tcp_numlayers ; i++){
                                temp_rate = 0 ;
-                               if (parameters->tcp_rates[i]== 0){
+                               if (img_fol->rates[i]== 0){
                                        parameters->tcp_rates[0]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ 
                                        (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy);
                                }else{
                                        temp_rate =((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ 
-                                               (parameters->tcp_rates[i] * 8 * image->comps[0].dx * image->comps[0].dy);
+                                               (img_fol->rates[i] * 8 * image->comps[0].dx * image->comps[0].dy);
                                        if (temp_rate > CINEMA_48_CS ){
                                                parameters->tcp_rates[0]= ((float) (image->numcomps * image->comps[0].w * image->comps[0].h * image->comps[0].prec))/ 
-                                                                                                                                                               (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy);
+                                               (CINEMA_48_CS * 8 * image->comps[0].dx * image->comps[0].dy);
+                                       }else{
+                                               parameters->tcp_rates[i]= img_fol->rates[i];
                                        }
                                }
                        }
@@ -582,14 +598,15 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
        fprintf(stream, "%d\n", cstr_info->prog);
        fprintf(stream, "%d %d\n", cstr_info->tile_x, cstr_info->tile_y);
        fprintf(stream, "%d %d\n", cstr_info->tw, cstr_info->th);
-       fprintf(stream, "%d\n", cstr_info->comp);
-       fprintf(stream, "%d\n", cstr_info->layer);
-       fprintf(stream, "%d\n", cstr_info->decomposition);
-       
-       for (resno = cstr_info->decomposition; resno >= 0; resno--) {
+       fprintf(stream, "%d\n", cstr_info->numcomps);
+       fprintf(stream, "%d\n", cstr_info->numlayers);
+       fprintf(stream, "%d\n", cstr_info->numdecompos);
+
+       for (resno = cstr_info->numdecompos[0]; resno >= 0; resno--) {
                fprintf(stream, "[%d,%d] ", 
-                       (1 << cstr_info->tile[0].pdx[resno]), (1 << cstr_info->tile[0].pdx[resno]));    /* based on tile 0 */
+                       (1 << cstr_info->tile[0].pdx[resno]), (1 << cstr_info->tile[0].pdx[resno]));    /* based on tile 0 and component 0 */
        }
+
        fprintf(stream, "\n");
 /* UniPG>> */
        fprintf(stream, "%d\n", cstr_info->main_head_start);
@@ -598,53 +615,50 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
        fprintf(stream, "%d\n", cstr_info->codestream_size);
        
        fprintf(stream, "\nINFO ON TILES\n");
-       fprintf(stream, "tileno start_pos  end_hd  end_tile"
-/* UniPG>> */
-               "   nbparts"
-/* <<UniPG */
-               "         disto     nbpix   disto/nbpix\n");
+       fprintf(stream, "tileno start_pos  end_hd  end_tile   nbparts         disto     nbpix   disto/nbpix\n");
        for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) {
-               fprintf(stream, "%4d %9d %9d %9d "
-/* UniPG>> */
-                       "%9d "
-/* <<UniPG */
-                       "%9e %9d %9e\n",
-                       cstr_info->tile[tileno].num_tile,
+               fprintf(stream, "%4d %9d %9d %9d %9d %9e %9d %9e\n",
+                       cstr_info->tile[tileno].tileno,
                        cstr_info->tile[tileno].start_pos,
                        cstr_info->tile[tileno].end_header,
                        cstr_info->tile[tileno].end_pos,
-/* UniPG>> */
                        cstr_info->tile[tileno].num_tps,
-/* <<UniPG */
-                       cstr_info->tile[tileno].distotile, cstr_info->tile[tileno].nbpix,
-                       cstr_info->tile[tileno].distotile / cstr_info->tile[tileno].nbpix);
+                       cstr_info->tile[tileno].distotile, cstr_info->tile[tileno].numpix,
+                       cstr_info->tile[tileno].distotile / cstr_info->tile[tileno].numpix);
        }
                
        for (tileno = 0; tileno < cstr_info->tw * cstr_info->th; tileno++) {
                int start_pos, end_ph_pos, end_pos;
                double disto = 0;
+               int max_numdecompos = 0;
                pack_nb = 0;
 
+               for (compno = 0; compno < cstr_info->numcomps; compno++) {
+                       if (max_numdecompos < cstr_info->numdecompos[compno])
+                               max_numdecompos = cstr_info->numdecompos[compno];
+               }       
+
                fprintf(stream, "\nTILE %d DETAILS\n", tileno); 
-/* UniPG>> */
-               fprintf(stream, "part_nb tileno   pack_nb start_pos end_tph_pos   end_pos\n");
+               fprintf(stream, "part_nb tileno  start_pack num_packs  start_pos end_tph_pos   end_pos\n");
                for (tilepartno = 0; tilepartno < cstr_info->tile[tileno].num_tps; tilepartno++)
-                       fprintf(stream, "%4d %9d %9d %9d %11d %9d\n",
+                       fprintf(stream, "%4d %9d   %9d %9d  %9d %11d %9d\n",
                                tilepartno, tileno,
-                               cstr_info->tile[tileno].tp_num[tilepartno],
-                               cstr_info->tile[tileno].tp_start_pos[tilepartno],
-                               cstr_info->tile[tileno].tp_end_header[tilepartno],
-                               cstr_info->tile[tileno].tp_end_pos[tilepartno]
+                               cstr_info->tile[tileno].tp[tilepartno].tp_start_pack,
+                               cstr_info->tile[tileno].tp[tilepartno].tp_numpacks,
+                               cstr_info->tile[tileno].tp[tilepartno].tp_start_pos,
+                               cstr_info->tile[tileno].tp[tilepartno].tp_end_header,
+                               cstr_info->tile[tileno].tp[tilepartno].tp_end_pos
                                );
-/* <<UniPG */
                if (cstr_info->prog == LRCP) {  /* LRCP */
-
                        fprintf(stream, "LRCP\npack_nb tileno layno resno compno precno start_pos end_ph_pos end_pos disto\n");
 
-                       for (layno = 0; layno < cstr_info->layer; layno++) {
-                               for (resno = 0; resno < cstr_info->decomposition + 1; resno++) {
-                                       for (compno = 0; compno < cstr_info->comp; compno++) {
-                                               int prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno];
+                       for (layno = 0; layno < cstr_info->numlayers; layno++) {
+                               for (resno = 0; resno < max_numdecompos + 1; resno++) {
+                                       for (compno = 0; compno < cstr_info->numcomps; compno++) {
+                                               int prec_max;
+                                               if (resno > cstr_info->numdecompos[compno])
+                                                       break;
+                                               prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno];
                                                for (precno = 0; precno < prec_max; precno++) {
                                                        start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos;
                                                        end_ph_pos = cstr_info->tile[tileno].packet[pack_nb].end_ph_pos;
@@ -663,10 +677,13 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
 
                        fprintf(stream, "RLCP\npack_nb tileno resno layno compno precno start_pos end_ph_pos end_pos disto\n");
 
-                       for (resno = 0; resno < cstr_info->decomposition + 1; resno++) {
-                               for (layno = 0; layno < cstr_info->layer; layno++) {
-                                       for (compno = 0; compno < cstr_info->comp; compno++) {
-                                               int prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno];
+                       for (resno = 0; resno < max_numdecompos + 1; resno++) {
+                               for (layno = 0; layno < cstr_info->numlayers; layno++) {
+                                       for (compno = 0; compno < cstr_info->numcomps; compno++) {
+                                               int prec_max; 
+                                               if (resno > cstr_info->numdecompos[compno])
+                                                       break;
+                                               prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno];
                                                for (precno = 0; precno < prec_max; precno++) {
                                                        start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos;
                                                        end_ph_pos = cstr_info->tile[tileno].packet[pack_nb].end_ph_pos;
@@ -685,25 +702,27 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
 
                        fprintf(stream, "RPCL\npack_nb tileno resno precno compno layno start_pos end_ph_pos end_pos disto\n"); 
 
-                       for (resno = 0; resno < cstr_info->decomposition + 1; resno++) {
+                       for (resno = 0; resno < max_numdecompos + 1; resno++) {
                                /* I suppose components have same XRsiz, YRsiz */
                                int x0 = cstr_info->tile_Ox + tileno - (int)floor((float)tileno/(float)cstr_info->tw ) * cstr_info->tw * cstr_info->tile_x;
                                int y0 = cstr_info->tile_Ox + (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tile_y;
                                int x1 = x0 + cstr_info->tile_x;
                                int y1 = y0 + cstr_info->tile_y;
-                               for (compno = 0; compno < cstr_info->comp; compno++) {
+                               for (compno = 0; compno < cstr_info->numcomps; compno++) {
                                        int prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno];
+                                       if (resno > cstr_info->numdecompos[compno])
+                                                       break;
                                        for (precno = 0; precno < prec_max; precno++) {
                                                int pcnx = cstr_info->tile[tileno].pw[resno];
-                                               int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->decomposition - resno );
-                                               int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->decomposition - resno );
+                                               int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->numdecompos[compno] - resno );
+                                               int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->numdecompos[compno] - resno );
                                                int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
                                                int precno_y = (int) floor( (float)precno/(float)pcnx );
                                                for(y = y0; y < y1; y++) {                                                      
                                                        if (precno_y*pcy == y ) {
                                                                for (x = x0; x < x1; x++) {                                                                     
                                                                        if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < cstr_info->layer; layno++) {
+                                                                               for (layno = 0; layno < cstr_info->numlayers; layno++) {
                                                                                        start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos;
                                                                                        end_ph_pos = cstr_info->tile[tileno].packet[pack_nb].end_ph_pos;
                                                                                        end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos;
@@ -730,20 +749,20 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
 
                        fprintf(stream, "PCRL\npack_nb tileno precno compno resno layno start_pos end_ph_pos end_pos disto\n"); 
 
-                       for (compno = 0; compno < cstr_info->comp; compno++) {
-                               for (resno = 0; resno < cstr_info->decomposition + 1; resno++) {
+                       for (compno = 0; compno < cstr_info->numcomps; compno++) {
+                               for (resno = 0; resno < cstr_info->numdecompos[compno] + 1; resno++) {
                                        int prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno];
                                        for (precno = 0; precno < prec_max; precno++) {
                                                int pcnx = cstr_info->tile[tileno].pw[resno];
-                                               int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->decomposition - resno );
-                                               int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->decomposition - resno );
+                                               int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->numdecompos[compno] - resno );
+                                               int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->numdecompos[compno] - resno );
                                                int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
                                                int precno_y = (int) floor( (float)precno/(float)pcnx );
                                                for(y = y0; y < y1; y++) {                                                      
                                                        if (precno_y*pcy == y ) {
                                                                for (x = x0; x < x1; x++) {                                                                     
                                                                        if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < cstr_info->layer; layno++) {
+                                                                               for (layno = 0; layno < cstr_info->numlayers; layno++) {
                                                                                        start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos;
                                                                                        end_ph_pos = cstr_info->tile[tileno].packet[pack_nb].end_ph_pos;
                                                                                        end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos;
@@ -765,26 +784,26 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
 
                        fprintf(stream, "CPRL\npack_nb tileno compno precno resno layno start_pos end_ph_pos end_pos disto\n"); 
 
-                       for (compno = 0; compno < cstr_info->comp; compno++) {
+                       for (compno = 0; compno < cstr_info->numcomps; compno++) {
                                /* I suppose components have same XRsiz, YRsiz */
                                int x0 = cstr_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tw * cstr_info->tile_x;
                                int y0 = cstr_info->tile_Ox + (int)floor( (float)tileno/(float)cstr_info->tw ) * cstr_info->tile_y;
                                int x1 = x0 + cstr_info->tile_x;
                                int y1 = y0 + cstr_info->tile_y;
                                
-                               for (resno = 0; resno < cstr_info->decomposition + 1; resno++) {
+                               for (resno = 0; resno < cstr_info->numdecompos[compno] + 1; resno++) {
                                        int prec_max = cstr_info->tile[tileno].pw[resno] * cstr_info->tile[tileno].ph[resno];
                                        for (precno = 0; precno < prec_max; precno++) {
                                                int pcnx = cstr_info->tile[tileno].pw[resno];
-                                               int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->decomposition - resno );
-                                               int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->decomposition - resno );
+                                               int pcx = (int) pow( 2, cstr_info->tile[tileno].pdx[resno] + cstr_info->numdecompos[compno] - resno );
+                                               int pcy = (int) pow( 2, cstr_info->tile[tileno].pdy[resno] + cstr_info->numdecompos[compno] - resno );
                                                int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
                                                int precno_y = (int) floor( (float)precno/(float)pcnx );
                                                for(y = y0; y < y1; y++) {
                                                        if (precno_y*pcy == y ) {
                                                                for (x = x0; x < x1; x++) {
                                                                        if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < cstr_info->layer; layno++) {
+                                                                               for (layno = 0; layno < cstr_info->numlayers; layno++) {
                                                                                        start_pos = cstr_info->tile[tileno].packet[pack_nb].start_pos;
                                                                                        end_ph_pos = cstr_info->tile[tileno].packet[pack_nb].end_ph_pos;
                                                                                        end_pos = cstr_info->tile[tileno].packet[pack_nb].end_pos;
@@ -806,6 +825,14 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
        
        fprintf(stream, "%8e\n", cstr_info->D_max); /* SE max */
        fprintf(stream, "%.8e\n", total_disto); /* SE totale */
+/* UniPG>> */
+       /* print the markers' list */
+       fprintf(stream, "\nMARKER LIST\n");
+       fprintf(stream, "%d\n", cstr_info->marknum);
+       fprintf(stream, "type\tstart_pos    length\n");
+       for (x = 0; x < cstr_info->marknum; x++)
+               fprintf(stream, "%X\t%9d %9d\n", cstr_info->marker[x].type, cstr_info->marker[x].pos, cstr_info->marker[x].len);
+/* <<UniPG */
        fclose(stream);
 
        fprintf(stderr,"Generated index file %s\n", index);
@@ -816,7 +843,7 @@ int write_index_file(opj_codestream_info_t *cstr_info, char *index) {
 /* ------------------------------------------------------------------------------------ */
 
 int parse_cmdline_encoder(int argc, char **argv, opj_cparameters_t *parameters,
-                                                                                                       img_fol_t *img_fol, raw_cparameters_t *raw_cp) {
+                                                                                                       img_fol_t *img_fol, raw_cparameters_t *raw_cp, char *indexfilename) {
        int i, j,totlen;
        option_t long_option[]={
                {"cinema2K",REQ_ARG, NULL ,'w'},
@@ -827,6 +854,7 @@ int parse_cmdline_encoder(int argc, char **argv, opj_cparameters_t *parameters,
                {"EPH",NO_ARG, NULL ,'E'},
                {"OutFor",REQ_ARG, NULL ,'O'},
                {"POC",REQ_ARG, NULL ,'P'},
+               {"ROI",REQ_ARG, NULL ,'R'},
        };
 
        /* parse the command line */
@@ -1087,8 +1115,7 @@ int parse_cmdline_encoder(int argc, char **argv, opj_cparameters_t *parameters,
                        case 'x':                       /* creation of index file */
                        {
                                char *index = optarg;
-                               strncpy(parameters->index, index, sizeof(parameters->index)-1);
-                               parameters->index_on = 1;
+                               strncpy(indexfilename, index, OPJ_PATH_LEN);
                        }
                        break;
 
@@ -1202,9 +1229,9 @@ int parse_cmdline_encoder(int argc, char **argv, opj_cparameters_t *parameters,
 
                        case 'R':                       /* ROI */
                        {
-                               if (sscanf(optarg, "OI:c=%d,U=%d", &parameters->roi_compno,
+                               if (sscanf(optarg, "c=%d,U=%d", &parameters->roi_compno,
                                            &parameters->roi_shift) != 2) {
-                                       fprintf(stderr, "ROI error !! [-ROI:c='compno',U='shift']\n");
+                                       fprintf(stderr, "ROI error !! [-ROI c='compno',U='shift']\n");
                                        return 1;
                                }
                        }
@@ -1302,9 +1329,8 @@ int parse_cmdline_encoder(int argc, char **argv, opj_cparameters_t *parameters,
                                int hprot, pprot, sens, addr, size, range;
 
                                /* we need to enable indexing */
-                               if (!parameters->index_on) {
-                                       strncpy(parameters->index, JPWL_PRIVATEINDEX_NAME, sizeof(parameters->index)-1);
-                                       parameters->index_on = 1;
+                               if (!indexfilename) {
+                                       strncpy(indexfilename, JPWL_PRIVATEINDEX_NAME, OPJ_PATH_LEN);
                                }
 
                                /* search for different protection methods */
@@ -1731,7 +1757,8 @@ int main(int argc, char **argv) {
        int imageno;
        dircnt_t *dirptr;
        raw_cparameters_t raw_cp;
-       opj_codestream_info_t cstr_info;
+       opj_codestream_info_t cstr_info;                /* Codestream information structure */
+       char indexfilename[OPJ_PATH_LEN];       /* index file name */
 
        /*
        configure the event callbacks (not required)
@@ -1745,15 +1772,20 @@ int main(int argc, char **argv) {
        /* set encoding parameters to default values */
        opj_set_default_encoder_parameters(&parameters);
 
-       /* need to initialize img_fol since parameters will be read in parse_cmdline_decoder */
+       /* Initialize indexfilename and img_fol */
+       *indexfilename = 0;
        memset(&img_fol,0,sizeof(img_fol_t));
 
        /* parse input and get user encoding parameters */
-       if(parse_cmdline_encoder(argc, argv, &parameters,&img_fol, &raw_cp) == 1) {
-               return 0;
+       if(parse_cmdline_encoder(argc, argv, &parameters,&img_fol, &raw_cp, indexfilename) == 1) {
+               return 1;
        }
        
        if (parameters.cp_cinema){
+               img_fol.rates = (float*)malloc(parameters.tcp_numlayers * sizeof(float));
+               for(i=0; i< parameters.tcp_numlayers; i++){
+                       img_fol.rates[i] = parameters.tcp_rates[i];
+               }
                cinema_parameters(&parameters);
        }                               
 
@@ -1882,7 +1914,7 @@ int main(int argc, char **argv) {
                        parameters.tcp_mct = image->numcomps == 3 ? 1 : 0;
 
                        if(parameters.cp_cinema){
-                               cinema_setup_encoder(&parameters,image);
+                               cinema_setup_encoder(&parameters,image,&img_fol);
                        }
 
                        /* encode the destination image */
@@ -1907,7 +1939,10 @@ int main(int argc, char **argv) {
                                cio = opj_cio_open((opj_common_ptr)cinfo, NULL, 0);
 
                                /* encode the image */
-                               bSuccess = opj_encode(cinfo, cio, image, &cstr_info);
+                               if (*indexfilename)                                     // If need to extract codestream information
+                                       bSuccess = opj_encode_with_info(cinfo, cio, image, &cstr_info);
+                               else
+                                       bSuccess = opj_encode(cinfo, cio, image, NULL);
                                if (!bSuccess) {
                                        opj_cio_close(cio);
                                        fprintf(stderr, "failed to encode image\n");
@@ -1929,8 +1964,8 @@ int main(int argc, char **argv) {
                                opj_cio_close(cio);
 
                                /* Write the index to disk */
-                               if (parameters.index_on) {
-                                       bSuccess = write_index_file(&cstr_info, parameters.index);
+                               if (*indexfilename) {
+                                       bSuccess = write_index_file(&cstr_info, indexfilename);
                                        if (bSuccess) {
                                                fprintf(stderr, "Failed to output index file\n");
                                        }
@@ -1938,6 +1973,8 @@ int main(int argc, char **argv) {
 
                                /* free remaining compression structures */
                                opj_destroy_compress(cinfo);
+                               if (*indexfilename)
+                                       opj_destroy_cstr_info(&cstr_info);
                        } else {                        /* JP2 format output */
                                int codestream_length;
                                opj_cio_t *cio = NULL;
@@ -1957,7 +1994,10 @@ int main(int argc, char **argv) {
                                cio = opj_cio_open((opj_common_ptr)cinfo, NULL, 0);
 
                                /* encode the image */
-                               bSuccess = opj_encode(cinfo, cio, image, &cstr_info);
+                               if (*indexfilename)                                     // If need to extract codestream information
+                                       bSuccess = opj_encode_with_info(cinfo, cio, image, &cstr_info);
+                               else
+                                       bSuccess = opj_encode(cinfo, cio, image, NULL);
                                if (!bSuccess) {
                                        opj_cio_close(cio);
                                        fprintf(stderr, "failed to encode image\n");
@@ -1978,8 +2018,8 @@ int main(int argc, char **argv) {
                                opj_cio_close(cio);
                                
                                /* Write the index to disk */
-                               if (parameters.index_on) {
-                                       bSuccess = write_index_file(&cstr_info, parameters.index);
+                               if (*indexfilename) {
+                                       bSuccess = write_index_file(&cstr_info, indexfilename);
                                        if (bSuccess) {
                                                fprintf(stderr, "Failed to output index file\n");
                                        }
@@ -1987,7 +2027,8 @@ int main(int argc, char **argv) {
 
                                /* free remaining compression structures */
                                opj_destroy_compress(cinfo);
-
+                               if (*indexfilename)
+                                       opj_destroy_cstr_info(&cstr_info);
                        }
        
                        /* free image data */