ENH: add Java j2kviewer but do not compile it using cmake since cmake has too poor...
[openjpeg.git] / libopenjpeg / j2k.c
index f8c5a3e06a792eec650dbe81571753c537ce2f7a..9a6d7458482a0998fbd019a4a2f9346552ee538a 100644 (file)
@@ -245,10 +245,10 @@ j2k_prog_order_t j2k_prog_order_list[] = {
        {PCRL, "PCRL"},
        {RLCP, "RLCP"},
        {RPCL, "RPCL"},
-       {-1, NULL}
+       {-1, ""}
 };
 
-char *convert_progression_order(OPJ_PROG_ORDER prg_order){
+char *j2k_convert_progression_order(OPJ_PROG_ORDER prg_order){
        j2k_prog_order_t *po;
        for(po = j2k_prog_order_list; po->enum_prog != -1; po++ ){
                if(po->enum_prog == prg_order){
@@ -258,6 +258,43 @@ char *convert_progression_order(OPJ_PROG_ORDER prg_order){
        return po->str_prog;
 }
 
+static void j2k_check_poc_val(opj_cparameters_t *parameters, int numcomps, int numlayers){
+       int index, resno, compno, layno, i;
+       char loss = 0;
+       int step_c = 1;
+       int step_r = numcomps * step_c;
+       int step_l = parameters->numresolution * step_r;
+       int array_size = step_l * numlayers * sizeof(int);
+       int *packet_array = (int *) opj_malloc(array_size);
+       
+       for (i = 0; i < parameters->numpocs ; i++) {
+               int layno0 = 0;
+               if(i > 0)
+                       layno0 = (parameters->POC[i].layno1 > parameters->POC[i-1].layno1 )? parameters->POC[i-1].layno1 : 0;
+               for (resno = parameters->POC[i].resno0 ; resno < parameters->POC[i].resno1 ; resno++) {
+                       for (compno = parameters->POC[i].compno0 ; compno < parameters->POC[i].compno1 ; compno++) {
+                               for (layno = layno0; layno < parameters->POC[i].layno1 ; layno++) {
+                                       index = step_r * resno + step_c * compno + step_l * layno;
+                                       packet_array[index]= 1;
+                               }
+                       }
+               }
+       }
+       for (resno = 0; resno < parameters->numresolution; resno++) {
+               for (compno = 0; compno < numcomps; compno++) {
+                       for (layno = 0; layno < numlayers ; layno++) {
+                               index = step_r * resno + step_c * compno + step_l * layno;
+                               if(!(   packet_array[index]== 1)){
+                                       loss = 1;
+                               }
+                       }
+               }
+       }
+       if(loss == 1)
+               fprintf(stdout,"Missing packets possible loss of data\n");
+       opj_free(packet_array);
+}
+
 void j2k_dump_image(FILE *fd, opj_image_t * img) {
        int compno;
        fprintf(fd, "image {\n");
@@ -332,12 +369,12 @@ void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_t * cp) {
 }
 
 /* ----------------------------------------------------------------------- */
-int j2k_get_num_tp(opj_cp_t *cp,int pino,int tileno){
+static int j2k_get_num_tp(opj_cp_t *cp,int pino,int tileno){
        char *prog;
        int i;
        int tpnum=1,tpend=0;
        opj_tcp_t *tcp = &cp->tcps[tileno];
-       prog = convert_progression_order(tcp->prg);
+       prog = j2k_convert_progression_order(tcp->prg);
        
        if(cp->tp_on == 1){
                for(i=0;i<4;i++){
@@ -368,19 +405,32 @@ int j2k_get_num_tp(opj_cp_t *cp,int pino,int tileno){
 }
 
 /**    mem allocation for TLM marker*/
-int j2k_generate_tlm(opj_cp_t *cp,int img_numcomp,opj_image_t *image ){
+int j2k_calculate_tp(opj_cp_t *cp,int img_numcomp,opj_image_t *image,opj_j2k_t *j2k ){
        int pino,tileno,maxres=0,totnum_tp=0;
+       j2k->cur_totnum_tp = (int *) opj_malloc(cp->tw * cp->th * sizeof(int));
        for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
+               int cur_totnum_tp = 0;
                opj_tcp_t *tcp = &cp->tcps[tileno];
                for(pino = 0; pino <= tcp->numpocs; pino++) {
                        int tp_num=0;
-                       opj_pi_iterator_t *pi = (opj_pi_iterator_t *) opj_malloc(sizeof(opj_pi_iterator_t));
-                       pi = pi_initialise_encode(image, cp, tileno,pino);
+                       opj_pi_iterator_t *pi = pi_initialise_encode(image, cp, tileno,FINAL_PASS);
                        if(!pi) { return -1;}
                        tp_num = j2k_get_num_tp(cp,pino,tileno);
                        totnum_tp = totnum_tp + tp_num;
+                       cur_totnum_tp = cur_totnum_tp + tp_num;
                        pi_destroy(pi, cp, tileno);
                }
+               j2k->cur_totnum_tp[tileno] = cur_totnum_tp;
+/* UniPG>> */
+               /* INDEX >> */
+               if (j2k->cstr_info && j2k->cstr_info->index_on) {
+                       j2k->cstr_info->tile[tileno].num_tps = cur_totnum_tp;
+                       j2k->cstr_info->tile[tileno].tp_start_pos = (int *) opj_malloc(cur_totnum_tp * sizeof(int));
+                       j2k->cstr_info->tile[tileno].tp_end_header = (int *) opj_malloc(cur_totnum_tp * sizeof(int));
+                       j2k->cstr_info->tile[tileno].tp_end_pos = (int *) opj_malloc(cur_totnum_tp      * sizeof(int));
+               }
+               /* << INDEX */
+/* <<UniPG */
        }
        return totnum_tp;
 }
@@ -445,7 +495,7 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
        cp->ty0 = cio_read(cio, 4);             /* YT0siz */
        
        image->numcomps = cio_read(cio, 2);     /* Csiz */
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
                /* if JPWL is on, we check whether TX errors have damaged
@@ -454,15 +504,19 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                        opj_event_msg(j2k->cinfo, EVT_ERROR,
                                "JPWL: bad image size (%d x %d)\n",
                                image->x1, image->y1);
-                       if (!JPWL_ASSUME || JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME || JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                }
                if (image->numcomps != ((len - 38) / 3)) {
                        opj_event_msg(j2k->cinfo, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
                                "JPWL: Csiz is %d => space in SIZ only for %d comps.!!!\n",
                                image->numcomps, ((len - 38) / 3));
-                       if (!JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                        /* we try to correct */
                        opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n");
                        if (image->numcomps < ((len - 38) / 3)) {
@@ -480,7 +534,7 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                cp->exp_comps = image->numcomps;
        }
 #endif /* USE_JPWL */
-/* <<UniPG */
+
        image->comps = (opj_image_comp_t *) opj_malloc(image->numcomps * sizeof(opj_image_comp_t));
        for (i = 0; i < image->numcomps; i++) {
                int tmp, w, h;
@@ -490,7 +544,6 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                image->comps[i].dx = cio_read(cio, 1);  /* XRsiz_i */
                image->comps[i].dy = cio_read(cio, 1);  /* YRsiz_i */
                
-/* UniPG>> */
 #ifdef USE_JPWL
                if (j2k->cp->correct) {
                /* if JPWL is on, we check whether TX errors have damaged
@@ -499,8 +552,10 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                                opj_event_msg(j2k->cinfo, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
                                        "JPWL: bad XRsiz_%d/YRsiz_%d (%d x %d)\n",
                                        i, i, image->comps[i].dx, image->comps[i].dy);
-                               if (!JPWL_ASSUME)
-                                       exit(1);
+                               if (!JPWL_ASSUME) {
+                                       opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                                       return;
+                               }
                                /* we try to correct */
                                opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust them\n");
                                if (!image->comps[i].dx) {
@@ -517,7 +572,6 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                        
                }
 #endif /* USE_JPWL */
-/* <<UniPG */
 
 
                /* TODO: unused ? */
@@ -530,7 +584,7 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
        
        cp->tw = int_ceildiv(image->x1 - cp->tx0, cp->tdx);
        cp->th = int_ceildiv(image->y1 - cp->ty0, cp->tdy);
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
                /* if JPWL is on, we check whether TX errors have damaged
@@ -539,8 +593,10 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                        opj_event_msg(j2k->cinfo, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
                                "JPWL: bad number of tiles (%d x %d)\n",
                                cp->tw, cp->th);
-                       if (!JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                        /* we try to correct */
                        opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust them\n");
                        if (cp->tw < 1) {
@@ -549,11 +605,10 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                                        cp->tw);
                        }
                        if (cp->tw > cp->max_tiles) {
-                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- too large x, increase expectance of %d\n",
-                                       cp->max_tiles);
                                cp->tw= 1;
-                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- setting %d tiles in x => HYPOTHESIS!!!\n",
-                                       cp->tw);
+                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- too large x, increase expectance of %d\n"
+                                       "- setting %d tiles in x => HYPOTHESIS!!!\n",
+                                       cp->max_tiles, cp->tw);
                        }
                        if (cp->th < 1) {
                                cp->th= 1;
@@ -561,32 +616,32 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                                        cp->th);
                        }
                        if (cp->th > cp->max_tiles) {
-                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- too large y, increase expectance of %d to continue\n",
-                                       cp->max_tiles);
                                cp->th= 1;
-                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- setting %d tiles in y => HYPOTHESIS!!!\n",
-                                       cp->th);
+                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- too large y, increase expectance of %d to continue\n",
+                                       "- setting %d tiles in y => HYPOTHESIS!!!\n",
+                                       cp->max_tiles, cp->th);
                        }
                }
        }
 #endif /* USE_JPWL */
-/* <<UniPG */
+
        cp->tcps = (opj_tcp_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tcp_t));
        cp->tileno = (int *) opj_malloc(cp->tw * cp->th * sizeof(int));
        cp->tileno_size = 0;
        
-/* UniPG>> */
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
                if (!cp->tcps) {
                        opj_event_msg(j2k->cinfo, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
                                "JPWL: could not alloc tcps field of cp\n");
-                       if (!JPWL_ASSUME || JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME || JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                }
        }
 #endif /* USE_JPWL */
-/* <<UniPG */  
+
        for (i = 0; i < cp->tw * cp->th; i++) {
                cp->tcps[i].POC = 0;
                cp->tcps[i].numpocs = 0;
@@ -807,7 +862,7 @@ static void j2k_read_qcx(opj_j2k_t *j2k, int compno, int len) {
        tccp->numgbits = tmp >> 5;
        numbands = (tccp->qntsty == J2K_CCP_QNTSTY_SIQNT) ? 
                1 : ((tccp->qntsty == J2K_CCP_QNTSTY_NOQNT) ? len - 1 : (len - 1) / 2);
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
 
@@ -816,18 +871,20 @@ static void j2k_read_qcx(opj_j2k_t *j2k, int compno, int len) {
                        opj_event_msg(j2k->cinfo, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
                                "JPWL: bad number of subbands in Sqcx (%d)\n",
                                numbands);
-                       if (!JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                        /* we try to correct */
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust them\n");
                        numbands = 1;
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- setting number of bands to %d => HYPOTHESIS!!!\n",
+                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust them\n"
+                               "- setting number of bands to %d => HYPOTHESIS!!!\n",
                                numbands);
                };
 
        };
 #endif /* USE_JPWL */
-/* <<UniPG */
+
        for (bandno = 0; bandno < numbands; bandno++) {
                int expn, mant;
                if (tccp->qntsty == J2K_CCP_QNTSTY_NOQNT) {
@@ -906,7 +963,7 @@ static void j2k_read_qcc(opj_j2k_t *j2k) {
        
        len = cio_read(cio, 2); /* Lqcc */
        compno = cio_read(cio, numcomp <= 256 ? 1 : 2); /* Cqcc */
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
 
@@ -917,20 +974,22 @@ static void j2k_read_qcc(opj_j2k_t *j2k) {
                        opj_event_msg(j2k->cinfo, EVT_ERROR,
                                "JPWL: bad component number in QCC (%d out of a maximum of %d)\n",
                                compno, numcomp);
-                       if (!JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                        /* we try to correct */
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n");
                        compno = backup_compno % numcomp;
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- setting component number to %d\n",
-                               compno);                                
+                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n"
+                               "- setting component number to %d\n",
+                               compno);
                }
 
                /* keep your private count of tiles */
                backup_compno++;
        };
 #endif /* USE_JPWL */
-/* <<UniPG */
+
        j2k_read_qcx(j2k, compno, len - 2 - (numcomp <= 256 ? 1 : 2));
 }
 
@@ -944,7 +1003,7 @@ static void j2k_write_poc(opj_j2k_t *j2k) {
        opj_tccp_t *tccp = &tcp->tccps[0];
        opj_cio_t *cio = j2k->cio;
 
-       numpchgs = tcp->numpocs;
+       numpchgs = 1 + tcp->numpocs;
        cio_write(cio, J2K_MS_POC, 2);  /* POC  */
        len = 2 + (5 + 2 * (numcomps <= 256 ? 1 : 2)) * numpchgs;
        cio_write(cio, len, 2);         /* Lpoc */
@@ -982,8 +1041,8 @@ static void j2k_read_poc(opj_j2k_t *j2k) {
                poc = &tcp->pocs[i];
                poc->resno0 = cio_read(cio, 1); /* RSpoc_i */
                poc->compno0 = cio_read(cio, numcomps <= 256 ? 1 : 2);  /* CSpoc_i */
-               poc->layno1 = int_min(cio_read(cio, 2), (unsigned int) tcp->numlayers); /* LYEpoc_i */
-               poc->resno1 = int_min(cio_read(cio, 1), (unsigned int) tccp->numresolutions);   /* REpoc_i */
+               poc->layno1 = cio_read(cio, 2);    /* LYEpoc_i */
+               poc->resno1 = cio_read(cio, 1);    /* REpoc_i */
                poc->compno1 = int_min(
                        cio_read(cio, numcomps <= 256 ? 1 : 2), (unsigned int) numcomps);       /* CEpoc_i */
                poc->prg = (OPJ_PROG_ORDER)cio_read(cio, 1);    /* Ppoc_i */
@@ -1091,7 +1150,7 @@ static void j2k_read_ppm(opj_j2k_t *j2k) {
                        cp->ppm_len = N_ppm;
                } else {                        /* NON-first PPM marker */
                        cp->ppm_data = (unsigned char *) opj_realloc(cp->ppm_data, (N_ppm +     cp->ppm_store) * sizeof(unsigned char));
-/* UniPG>> */
+
 #ifdef USE_JPWL
                        /* this memory allocation check could be done even in non-JPWL cases */
                        if (cp->correct) {
@@ -1099,12 +1158,15 @@ static void j2k_read_ppm(opj_j2k_t *j2k) {
                                        opj_event_msg(j2k->cinfo, EVT_ERROR,
                                                "JPWL: failed memory allocation during PPM marker parsing (pos. %x)\n",
                                                cio_tell(cio));
-                                       if (!JPWL_ASSUME || JPWL_ASSUME)
-                                               exit(1);
+                                       if (!JPWL_ASSUME || JPWL_ASSUME) {
+                                               free(cp->ppm_data);
+                                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                                               return;
+                                       }
                                }
                        }
 #endif
-/* <<UniPG */
+
                        cp->ppm_data_first = cp->ppm_data;
                        cp->ppm_len = N_ppm + cp->ppm_store;
                }
@@ -1172,7 +1234,7 @@ static void j2k_write_sot(opj_j2k_t *j2k) {
        cio_write(cio, j2k->curtileno, 2);      /* Isot */
        cio_skip(cio, 4);                                       /* Psot (further in j2k_write_sod) */
        cio_write(cio, j2k->cur_tp_num , 1);    /* TPsot */
-       cio_write(cio, j2k->cur_totnum_tp, 1);          /* TNsot */
+       cio_write(cio, j2k->cur_totnum_tp[j2k->curtileno], 1);          /* TNsot */
        len = cio_tell(cio) - lenp;
        cio_seek(cio, lenp);
        cio_write(cio, len, 2);                         /* Lsot */
@@ -1189,7 +1251,7 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
        
        len = cio_read(cio, 2);
        tileno = cio_read(cio, 2);
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
 
@@ -1200,20 +1262,22 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
                        opj_event_msg(j2k->cinfo, EVT_ERROR,
                                "JPWL: bad tile number (%d out of a maximum of %d)\n",
                                tileno, (cp->tw * cp->th));
-                       if (!JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                        /* we try to correct */
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n");
                        tileno = backup_tileno;
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- setting tile number to %d\n",
-                               tileno);                                
+                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n"
+                               "- setting tile number to %d\n",
+                               tileno);
                }
 
                /* keep your private count of tiles */
                backup_tileno++;
        };
 #endif /* USE_JPWL */
-/* <<UniPG */
+
        
        if (cp->tileno_size == 0) {
                cp->tileno[cp->tileno_size] = tileno;
@@ -1231,7 +1295,7 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
        }
        
        totlen = cio_read(cio, 4);
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
 
@@ -1240,19 +1304,20 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
                        opj_event_msg(j2k->cinfo, EVT_ERROR,
                                "JPWL: bad tile byte size (%d bytes against %d bytes left)\n",
                                totlen, cio_numbytesleft(cio) + 8);
-                       if (!JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                        /* we try to correct */
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n");
                        totlen = 0;
-                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- setting Psot to %d => "
-                               "assuming it is the last tile\n",
-                               totlen);                                
+                       opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n"
+                               "- setting Psot to %d => assuming it is the last tile\n",
+                               totlen);
                }
 
        };
 #endif /* USE_JPWL */
-/* <<UniPG */
+
        if (!totlen)
                totlen = cio_numbytesleft(cio) + 8;
        
@@ -1285,24 +1350,29 @@ static void j2k_write_sod(opj_j2k_t *j2k, void *tile_coder) {
        int l, layno;
        int totlen;
        opj_tcp_t *tcp = NULL;
-       opj_image_info_t *image_info = NULL;
+       opj_codestream_info_t *cstr_info = NULL;
        
        opj_tcd_t *tcd = (opj_tcd_t*)tile_coder;        /* cast is needed because of conflicts in header inclusions */
        opj_cp_t *cp = j2k->cp;
        opj_cio_t *cio = j2k->cio;
 
+       tcd->tp_num = j2k->tp_num ;
        tcd->cur_tp_num = j2k->cur_tp_num;
-       tcd->cur_totnum_tp = j2k->cur_totnum_tp;
        
        cio_write(cio, J2K_MS_SOD, 2);
        if (j2k->curtileno == 0) {
                j2k->sod_start = cio_tell(cio) + j2k->pos_correction;
        }
-       
+
        /* INDEX >> */
-       image_info = j2k->image_info;
-       if (image_info && image_info->index_on) {
-               image_info->tile[j2k->curtileno].end_header = cio_tell(cio) + j2k->pos_correction - 1;
+       cstr_info = j2k->cstr_info;
+       if (cstr_info && cstr_info->index_on) {
+               if (!j2k->cur_tp_num ){
+                       cstr_info->tile[j2k->curtileno].end_header = cio_tell(cio) + j2k->pos_correction - 1;
+               }else{
+                       if(cstr_info->tile[j2k->curtileno].packet[cstr_info->num - 1].end_pos < cio_tell(cio))
+                               cstr_info->tile[j2k->curtileno].packet[cstr_info->num].start_pos = cio_tell(cio);
+               }
        }
        /* << INDEX */
        
@@ -1310,11 +1380,11 @@ static void j2k_write_sod(opj_j2k_t *j2k, void *tile_coder) {
        for (layno = 0; layno < tcp->numlayers; layno++) {
                tcp->rates[layno] -= tcp->rates[layno] ? (j2k->sod_start / (cp->th * cp->tw)) : 0;
        }
-       if(image_info) {
-               image_info->num = 0;
+       if(cstr_info && (j2k->cur_tp_num == 0)) {
+               cstr_info->num = 0;
        }
        
-       l = tcd_encode_tile(tcd, j2k->curtileno, cio_getbp(cio), cio_numbytesleft(cio) - 2, image_info);
+       l = tcd_encode_tile(tcd, j2k->curtileno, cio_getbp(cio), cio_numbytesleft(cio) - 2, cstr_info);
        
        /* Writing Psot in SOT marker */
        totlen = cio_tell(cio) + l - j2k->sot_start;
@@ -1390,7 +1460,7 @@ static void j2k_read_rgn(opj_j2k_t *j2k) {
        len = cio_read(cio, 2);                                                                         /* Lrgn */
        compno = cio_read(cio, numcomps <= 256 ? 1 : 2);                        /* Crgn */
        roisty = cio_read(cio, 1);                                                                      /* Srgn */
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
                /* totlen is negative or larger than the bytes left!!! */
@@ -1398,12 +1468,14 @@ static void j2k_read_rgn(opj_j2k_t *j2k) {
                        opj_event_msg(j2k->cinfo, EVT_ERROR,
                                "JPWL: bad component number in RGN (%d when there are only %d)\n",
                                compno, numcomps);
-                       if (!JPWL_ASSUME || JPWL_ASSUME)
-                               exit(1);
+                       if (!JPWL_ASSUME || JPWL_ASSUME) {
+                               opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
+                               return;
+                       }
                }
        };
 #endif /* USE_JPWL */
-/* <<UniPG */
+
        tcp->tccps[compno].roishift = cio_read(cio, 1);                         /* SPrgn */
 }
 
@@ -1416,24 +1488,29 @@ static void j2k_write_eoc(opj_j2k_t *j2k) {
 static void j2k_read_eoc(opj_j2k_t *j2k) {
        int i, tileno;
 
-#ifndef NO_PACKETS_DECODING  
-       opj_tcd_t *tcd = tcd_create(j2k->cinfo);
-       tcd_malloc_decode(tcd, j2k->image, j2k->cp);
-       for (i = 0; i < j2k->cp->tileno_size; i++) {
-               tileno = j2k->cp->tileno[i];
-               tcd_decode_tile(tcd, j2k->tile_data[tileno], j2k->tile_len[tileno], tileno);
-               opj_free(j2k->tile_data[tileno]);
-               j2k->tile_data[tileno] = NULL;
+       /* if packets should be decoded */
+       if (j2k->cp->limit_decoding != DECODE_ALL_BUT_PACKETS) {
+               opj_tcd_t *tcd = tcd_create(j2k->cinfo);
+               tcd_malloc_decode(tcd, j2k->image, j2k->cp);
+               for (i = 0; i < j2k->cp->tileno_size; i++) {
+                       tcd_malloc_decode_tile(tcd, j2k->image, j2k->cp, i);
+                       tileno = j2k->cp->tileno[i];
+                       tcd_decode_tile(tcd, j2k->tile_data[tileno], j2k->tile_len[tileno], tileno);
+                       opj_free(j2k->tile_data[tileno]);
+                       j2k->tile_data[tileno] = NULL;
+                       tcd_free_decode_tile(tcd, i);
+               }
+               tcd_free_decode(tcd);
+               tcd_destroy(tcd);
        }
-       tcd_free_decode(tcd);
-       tcd_destroy(tcd);
-#else 
-       for (i = 0; i < j2k->cp->tileno_size; i++) {
-               tileno = j2k->cp->tileno[i];
-               opj_free(j2k->tile_data[tileno]);
-               j2k->tile_data[tileno] = NULL;
+       /* if packets should not be decoded  */
+       else {
+               for (i = 0; i < j2k->cp->tileno_size; i++) {
+                       tileno = j2k->cp->tileno[i];
+                       opj_free(j2k->tile_data[tileno]);
+                       j2k->tile_data[tileno] = NULL;
+               }
        }
-#endif
        
        j2k->state = J2K_STATE_MT;
 }
@@ -1467,20 +1544,24 @@ opj_dec_mstabent_t j2k_dec_mstab[] = {
   {J2K_MS_SOP, 0, 0},
   {J2K_MS_CRG, J2K_STATE_MH, j2k_read_crg},
   {J2K_MS_COM, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_com},
-/* UniPG>> */
+
 #ifdef USE_JPWL
   {J2K_MS_EPC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_epc},
   {J2K_MS_EPB, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_epb},
   {J2K_MS_ESD, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_esd},
   {J2K_MS_RED, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_red},
 #endif /* USE_JPWL */
-/* <<UniPG */
+#ifdef USE_JPSEC
+  {J2K_MS_SEC, J2K_STATE_MH, j2k_read_sec},
+  {J2K_MS_INSEC, 0, j2k_read_insec},
+#endif /* USE_JPSEC */
+
   {0, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_unk}
 };
 
 static void j2k_read_unk(opj_j2k_t *j2k) {
        opj_event_msg(j2k->cinfo, EVT_WARNING, "Unknown marker\n");
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
                int m = 0, id, i;
@@ -1492,8 +1573,8 @@ static void j2k_read_unk(opj_j2k_t *j2k) {
                        id);
                if (!JPWL_ASSUME) {
                        opj_event_msg(j2k->cinfo, EVT_ERROR,
-                               "- possible synch loss due to uncorrectable channel errors => Exiting\n");
-                       exit(1);
+                               "- possible synch loss due to uncorrectable codestream errors => giving up\n");
+                       return;
                }
                /* OK, activate this at your own risk!!! */
                /* we look for the marker at the minimum hamming distance from this */
@@ -1537,7 +1618,7 @@ static void j2k_read_unk(opj_j2k_t *j2k) {
 
        };
 #endif /* USE_JPWL */
-/* <<UniPG */
+
 }
 
 /**
@@ -1626,13 +1707,13 @@ void j2k_setup_decoder(opj_j2k_t *j2k, opj_dparameters_t *parameters) {
                cp->reduce = parameters->cp_reduce;     
                cp->layer = parameters->cp_layer;
                cp->limit_decoding = parameters->cp_limit_decoding;
-/* UniPG>> */
+
 #ifdef USE_JPWL
                cp->correct = parameters->jpwl_correct;
                cp->exp_comps = parameters->jpwl_exp_comps;
                cp->max_tiles = parameters->jpwl_max_tiles;
 #endif /* USE_JPWL */
-/* <<UniPG */
+
 
                /* keep a link to cp so that we can destroy it later in j2k_destroy_decompress */
                j2k->cp = cp;
@@ -1656,7 +1737,7 @@ opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio) {
                opj_dec_mstabent_t *e;
                int id = cio_read(cio, 2);
 
-/* UniPG>> */
+
 #ifdef USE_JPWL
                /* we try to honor JPWL correction power */
                if (j2k->cp->correct) {
@@ -1675,23 +1756,26 @@ opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio) {
 
                        /* check whether it begins with ff */
                        if (id >> 8 != 0xff) {
-                               opj_event_msg(j2k->cinfo, EVT_ERROR,
+                               opj_event_msg(cinfo, EVT_ERROR,
                                        "JPWL: possible bad marker %x at %d\n",
                                        id, cio_tell(cio) - 2);
-                               if (!JPWL_ASSUME)
-                                       exit(1);
+                               if (!JPWL_ASSUME) {
+                                       opj_image_destroy(image);
+                                       opj_event_msg(cinfo, EVT_ERROR, "JPWL: giving up\n");
+                                       return 0;
+                               }
                                /* we try to correct */
-                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- trying to adjust this\n");
                                id = id | 0xff00;
                                cio_seek(cio, cio_tell(cio) - 2);
                                cio_write(cio, id, 2);
-                               opj_event_msg(j2k->cinfo, EVT_WARNING, "- setting marker to %x\n",
-                                       id);                            
+                               opj_event_msg(cinfo, EVT_WARNING, "- trying to adjust this\n"
+                                       "- setting marker to %x\n",
+                                       id);
                        }
 
                }
 #endif /* USE_JPWL */
-/* <<UniPG */
+
                if (id >> 8 != 0xff) {
                        opj_image_destroy(image);
                        opj_event_msg(cinfo, EVT_ERROR, "%.8x: expected a marker instead of %x\n", cio_tell(cio) - 2, id);
@@ -1746,6 +1830,7 @@ opj_image_t* j2k_decode_jpt_stream(opj_j2k_t *j2k, opj_cio_t *cio) {
 
        /* create an empty image */
        image = opj_image_create0();
+       j2k->image = image;
 
        j2k->state = J2K_STATE_MHSOC;
        
@@ -1830,18 +1915,24 @@ void j2k_destroy_compress(opj_j2k_t *j2k) {
 
        if(!j2k) return;
 
-       if(j2k->image_info != NULL) {
-               opj_image_info_t *image_info = j2k->image_info;
-               if (image_info->index_on && j2k->cp) {
+       if(j2k->cstr_info != NULL) {
+               opj_codestream_info_t *cstr_info = j2k->cstr_info;
+               if (cstr_info->index_on && j2k->cp) {
                        opj_cp_t *cp = j2k->cp;
                        for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
-                               opj_tile_info_t *tile_info = &image_info->tile[tileno];
+                               opj_tile_info_t *tile_info = &cstr_info->tile[tileno];
                                opj_free(tile_info->thresh);
                                opj_free(tile_info->packet);
+/* UniPG>> */
+                /* INDEX >> */
+                opj_free(tile_info->tp_start_pos);
+                opj_free(tile_info->tp_end_header);
+                opj_free(tile_info->tp_end_pos);
+                /* << INDEX */
+/* <<UniPG */
                        }
-                       opj_free(image_info->tile);
+                       opj_free(cstr_info->tile);
                }
-               opj_free(image_info);
        }
        if(j2k->cp != NULL) {
                opj_cp_t *cp = j2k->cp;
@@ -1884,6 +1975,7 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
        copy user encoding parameters 
        */
        cp->cinema = parameters->cp_cinema;
+       cp->max_comp_size =     parameters->max_comp_size;
        cp->rsiz   = parameters->cp_rsiz;
        cp->disto_alloc = parameters->cp_disto_alloc;
        cp->fixed_alloc = parameters->cp_fixed_alloc;
@@ -1898,9 +1990,6 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
 
        /* creation of an index file ? */
        cp->index_on = parameters->index_on;
-       if(cp->index_on) {
-               j2k->image_info = (opj_image_info_t*)opj_malloc(sizeof(opj_image_info_t));
-       }
 
        /* tiles */
        cp->tdx = parameters->cp_tdx;
@@ -1940,7 +2029,7 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
        cp->img_size += (image->comps[i].w *image->comps[i].h * image->comps[i].prec);
        }
 
-/* UniPG>> */
+
 #ifdef USE_JPWL
        /*
        calculate JPWL encoding parameters
@@ -1996,7 +2085,7 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                cp->epc_on = false;
        }
 #endif /* USE_JPWL */
-/* <<UniPG */
+
 
        /* initialize the mutiple tiles */
        /* ---------------------------- */
@@ -2006,21 +2095,29 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                opj_tcp_t *tcp = &cp->tcps[tileno];
                tcp->numlayers = parameters->tcp_numlayers;
                for (j = 0; j < tcp->numlayers; j++) {
-                       if (cp->fixed_quality) {        /* add fixed_quality */
-                               tcp->distoratio[j] = parameters->tcp_distoratio[j];
-                       } else {
+                       if(cp->cinema){
+                               if (cp->fixed_quality) {
+                                       tcp->distoratio[j] = parameters->tcp_distoratio[j];
+                               }
                                tcp->rates[j] = parameters->tcp_rates[j];
+                       }else{
+                               if (cp->fixed_quality) {        /* add fixed_quality */
+                                       tcp->distoratio[j] = parameters->tcp_distoratio[j];
+                               } else {
+                                       tcp->rates[j] = parameters->tcp_rates[j];
+                               }
                        }
                }
                tcp->csty = parameters->csty;
                tcp->prg = parameters->prog_order;
-               tcp->mct = image->numcomps == 3 ? 1 : 0;
+               tcp->mct = parameters->tcp_mct; 
 
                numpocs_tile = 0;
                tcp->POC = 0;
                if (parameters->numpocs) {
                        /* initialisation of POC */
                        tcp->POC = 1;
+                       j2k_check_poc_val(parameters, image->numcomps, tcp->numlayers);
                        for (i = 0; i < parameters->numpocs; i++) {
                                if((tileno == parameters->POC[i].tile - 1) || (parameters->POC[i].tile == -1)) {
                                        opj_poc_t *tcp_poc = &tcp->pocs[numpocs_tile];
@@ -2029,13 +2126,15 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                                        tcp_poc->layno1         = parameters->POC[numpocs_tile].layno1;
                                        tcp_poc->resno1         = parameters->POC[numpocs_tile].resno1;
                                        tcp_poc->compno1        = parameters->POC[numpocs_tile].compno1;
-                                       tcp_poc->prg            = parameters->POC[numpocs_tile].prg;
+                                       tcp_poc->prg1           = parameters->POC[numpocs_tile].prg1;
                                        tcp_poc->tile           = parameters->POC[numpocs_tile].tile;
                                        numpocs_tile++;
                                }
                        }
+                       tcp->numpocs = numpocs_tile -1 ;
+               }else{ 
+                       tcp->numpocs = 0;
                }
-               tcp->numpocs = numpocs_tile;
 
                tcp->tccps = (opj_tccp_t *) opj_malloc(image->numcomps * sizeof(opj_tccp_t));
                
@@ -2116,244 +2215,8 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
        }
 }
 
-/**
-Create an index file
-@param j2k
-@param cio
-@param image_info
-@param index Index filename
-@return Returns 1 if successful, returns 0 otherwise
-*/
-static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *image_info, char *index) {
-       int tileno, compno, layno, resno, precno, pack_nb, x, y;
-       FILE *stream = NULL;
-       double total_disto = 0;
-
-       image_info->codestream_size = cio_tell(cio) + j2k->pos_correction;      /* Correction 14/4/03 suite rmq de Patrick */
-
-/* UniPG>> */
-#ifdef USE_JPWL
-       /* if JPWL is enabled and the name coincides with our own set
-          then discard the creation of the file: this was just done to
-          enable indexing, we do not want an index file
-       */
-       if (j2k->cp->epc_on && !strcmp(index, JPWL_PRIVATEINDEX_NAME))
-               return 1;
-#endif /* USE_JPWL */
-/* <<UniPG */
-
-       stream = fopen(index, "w");
-       if (!stream) {
-               opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to open %s for writing\n", index);
-               return 0;
-       }
-       
-       fprintf(stream, "%d %d\n", image_info->image_w, image_info->image_h);
-       fprintf(stream, "%d\n", image_info->prog);
-       fprintf(stream, "%d %d\n", image_info->tile_x, image_info->tile_y);
-       fprintf(stream, "%d %d\n", image_info->tw, image_info->th);
-       fprintf(stream, "%d\n", image_info->comp);
-       fprintf(stream, "%d\n", image_info->layer);
-       fprintf(stream, "%d\n", image_info->decomposition);
-       
-       for (resno = image_info->decomposition; resno >= 0; resno--) {
-               fprintf(stream, "[%d,%d] ", 
-                       (1 << image_info->tile[0].pdx[resno]), (1 << image_info->tile[0].pdx[resno]));  /* based on tile 0 */
-       }
-       fprintf(stream, "\n");
-       fprintf(stream, "%d\n", image_info->main_head_end);
-       fprintf(stream, "%d\n", image_info->codestream_size);
-       
-       for (tileno = 0; tileno < image_info->tw * image_info->th; tileno++) {
-               fprintf(stream, "%4d %9d %9d %9d %9e %9d %9e\n",
-                       image_info->tile[tileno].num_tile,
-                       image_info->tile[tileno].start_pos,
-                       image_info->tile[tileno].end_header,
-                       image_info->tile[tileno].end_pos,
-                       image_info->tile[tileno].distotile, image_info->tile[tileno].nbpix,
-                       image_info->tile[tileno].distotile / image_info->tile[tileno].nbpix);
-       }
-       
-       for (tileno = 0; tileno < image_info->tw * image_info->th; tileno++) {
-               int start_pos, end_pos;
-               double disto = 0;
-               pack_nb = 0;
-               
-               /*
-               fprintf(stream, "pkno tileno layerno resno compno precno start_pos   end_pos       deltaSE        \n");
-               */
-               
-               if (image_info->prog == LRCP) { /* LRCP */
-                       /*
-                       fprintf(stream, "pack_nb tileno layno resno compno precno start_pos  end_pos   disto");
-                       */
-                       for (layno = 0; layno < image_info->layer; layno++) {
-                               for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                                       for (compno = 0; compno < image_info->comp; compno++) {
-                                               int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                               for (precno = 0; precno < prec_max; precno++) {
-                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                       fprintf(stream, "%4d %6d %7d %5d %6d %6d %9d %9d %8e\n",
-                                                               pack_nb, tileno, layno, resno, compno, precno, start_pos, end_pos, disto);
-                                                       total_disto += disto;
-                                                       pack_nb++;
-                                               }
-                                       }
-                               }
-                       }
-               } /* LRCP */
-               else if (image_info->prog == RLCP) {    /* RLCP */
-                       /*
-                       fprintf(stream, "pack_nb tileno resno layno compno precno start_pos  end_pos   disto");
-                       */
-                       for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                               for (layno = 0; layno < image_info->layer; layno++) {
-                                       for (compno = 0; compno < image_info->comp; compno++) {
-                                               int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                               for (precno = 0; precno < prec_max; precno++) {
-                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                       fprintf(stream, "%4d %6d %5d %7d %6d %6d %9d %9d %8e\n",
-                                                               pack_nb, tileno, resno, layno, compno, precno, start_pos, end_pos, disto);
-                                                       total_disto += disto;
-                                                       pack_nb++;
-                                               }
-                                       }
-                               }
-                       }
-               } /* RLCP */
-               else if (image_info->prog == RPCL) {    /* RPCL */
-                       /*
-                       fprintf(stream, "\npack_nb tileno resno precno compno layno start_pos  end_pos   disto\n"); 
-                       */
-                       for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                               /* I suppose components have same XRsiz, YRsiz */
-                               int x0 = image_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tw * image_info->tile_x;
-                               int y0 = image_info->tile_Ox + (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tile_y;
-                               int x1 = x0 + image_info->tile_x;
-                               int y1 = y0 + image_info->tile_y;
-                               for(y = y0; y < y1; y++) {
-                                       for(x = x0; x < x1; x++) {
-                                               for (compno = 0; compno < image_info->comp; compno++) {
-                                                       int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                                       for (precno = 0; precno < prec_max; precno++) {
-                                                               int pcnx = image_info->tile[tileno].pw[resno];
-                                                               int pcx = (int) pow( 2, image_info->tile[tileno].pdx[resno] + image_info->decomposition - resno );
-                                                               int pcy = (int) pow( 2, image_info->tile[tileno].pdy[resno] + image_info->decomposition - resno );
-                                                               int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
-                                                               int precno_y = (int) floor( (float)precno/(float)pcnx );
-                                                               if (precno_y*pcy == y ) {
-                                                                       if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < image_info->layer; layno++) {
-                                                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                                                       fprintf(stream, "%4d %6d %5d %6d %6d %7d %9d %9d %8e\n",
-                                                                                               pack_nb, tileno, resno, precno, compno, layno, start_pos, end_pos, disto); 
-                                                                                       total_disto += disto;
-                                                                                       pack_nb++; 
-                                                                               }
-                                                                       }
-                                                               }
-                                                       } /* precno */
-                                               } /* compno */
-                                       } /* x = x0..x1 */
-                               } /* y = y0..y1 */
-                       } /* resno */
-               } /* RPCL */
-               else if (image_info->prog == PCRL) {    /* PCRL */
-                       /* I suppose components have same XRsiz, YRsiz */
-                       int x0 = image_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tw * image_info->tile_x;
-                       int y0 = image_info->tile_Ox + (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tile_y;
-                       int x1 = x0 + image_info->tile_x;
-                       int y1 = y0 + image_info->tile_y;
-                       /*
-                       fprintf(stream, "\npack_nb tileno precno compno resno layno start_pos  end_pos   disto\n"); 
-                       */
-                       for(y = y0; y < y1; y++) {
-                               for(x = x0; x < x1; x++) {
-                                       for (compno = 0; compno < image_info->comp; compno++) {
-                                               for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                                                       int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                                       for (precno = 0; precno < prec_max; precno++) {
-                                                               int pcnx = image_info->tile[tileno].pw[resno];
-                                                               int pcx = (int) pow( 2, image_info->tile[tileno].pdx[resno] + image_info->decomposition - resno );
-                                                               int pcy = (int) pow( 2, image_info->tile[tileno].pdy[resno] + image_info->decomposition - resno );
-                                                               int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
-                                                               int precno_y = (int) floor( (float)precno/(float)pcnx );
-                                                               if (precno_y*pcy == y ) {
-                                                                       if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < image_info->layer; layno++) {
-                                                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                                                       fprintf(stream, "%4d %6d %6d %6d %5d %7d %9d %9d %8e\n",
-                                                                                               pack_nb, tileno, precno, compno, resno, layno, start_pos, end_pos, disto); 
-                                                                                       total_disto += disto;
-                                                                                       pack_nb++; 
-                                                                               }
-                                                                       }
-                                                               }
-                                                       } /* precno */
-                                               } /* resno */
-                                       } /* compno */
-                               } /* x = x0..x1 */
-                       } /* y = y0..y1 */
-               } /* PCRL */
-               else {  /* CPRL */
-                       /*
-                       fprintf(stream, "\npack_nb tileno compno precno resno layno start_pos  end_pos   disto\n"); 
-                       */
-                       for (compno = 0; compno < image_info->comp; compno++) {
-                               /* I suppose components have same XRsiz, YRsiz */
-                               int x0 = image_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tw * image_info->tile_x;
-                               int y0 = image_info->tile_Ox + (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tile_y;
-                               int x1 = x0 + image_info->tile_x;
-                               int y1 = y0 + image_info->tile_y;
-                               for(y = y0; y < y1; y++) {
-                                       for(x = x0; x < x1; x++) {
-                                               for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                                                       int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                                       for (precno = 0; precno < prec_max; precno++) {
-                                                               int pcnx = image_info->tile[tileno].pw[resno];
-                                                               int pcx = (int) pow( 2, image_info->tile[tileno].pdx[resno] + image_info->decomposition - resno );
-                                                               int pcy = (int) pow( 2, image_info->tile[tileno].pdy[resno] + image_info->decomposition - resno );
-                                                               int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
-                                                               int precno_y = (int) floor( (float)precno/(float)pcnx );
-                                                               if (precno_y*pcy == y ) {
-                                                                       if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < image_info->layer; layno++) {
-                                                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                                                       fprintf(stream, "%4d %6d %6d %6d %5d %7d %9d %9d %8e\n",
-                                                                                               pack_nb, tileno, compno, precno, resno, layno, start_pos, end_pos, disto); 
-                                                                                       total_disto += disto;
-                                                                                       pack_nb++; 
-                                                                               }
-                                                                       }
-                                                               }
-                                                       } /* precno */
-                                               } /* resno */
-                                       } /* x = x0..x1 */
-                               } /* y = y0..y1 */
-                       } /* comno */
-               } /* CPRL */   
-       } /* tileno */
-       
-       fprintf(stream, "%8e\n", image_info->D_max); /* SE max */
-       fprintf(stream, "%.8e\n", total_disto); /* SE totale */
-       fclose(stream);
-
-       return 1;
-}
-
-bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index) {
+bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, opj_codestream_info_t *cstr_info) {
        int tileno, compno;
-       opj_image_info_t *image_info = NULL;
        opj_cp_t *cp = NULL;
 
        opj_tcd_t *tcd = NULL;  /* TCD component */
@@ -2366,35 +2229,34 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
        /* j2k_dump_cp(stdout, image, cp); */
 
        /* INDEX >> */
-       image_info = j2k->image_info;
-       if (image_info && cp->index_on) {
-               image_info->index_on = cp->index_on;
-               image_info->tile = (opj_tile_info_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tile_info_t));
-               image_info->image_w = image->x1 - image->x0;
-               image_info->image_h = image->y1 - image->y0;
-               image_info->prog = (&cp->tcps[0])->prg;
-               image_info->tw = cp->tw;
-               image_info->th = cp->th;
-               image_info->tile_x = cp->tdx;   /* new version parser */
-               image_info->tile_y = cp->tdy;   /* new version parser */
-               image_info->tile_Ox = cp->tx0;  /* new version parser */
-               image_info->tile_Oy = cp->ty0;  /* new version parser */
-               image_info->comp = image->numcomps;
-               image_info->layer = (&cp->tcps[0])->numlayers;
-               image_info->decomposition = (&cp->tcps[0])->tccps->numresolutions - 1;
-               image_info->D_max = 0;          /* ADD Marcela */
+       j2k->cstr_info = cstr_info;
+       if (cstr_info && cp->index_on) {
+               cstr_info->index_on = cp->index_on;
+               cstr_info->tile = (opj_tile_info_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tile_info_t));
+               cstr_info->image_w = image->x1 - image->x0;
+               cstr_info->image_h = image->y1 - image->y0;
+               cstr_info->prog = (&cp->tcps[0])->prg;
+               cstr_info->tw = cp->tw;
+               cstr_info->th = cp->th;
+               cstr_info->tile_x = cp->tdx;    /* new version parser */
+               cstr_info->tile_y = cp->tdy;    /* new version parser */
+               cstr_info->tile_Ox = cp->tx0;   /* new version parser */
+               cstr_info->tile_Oy = cp->ty0;   /* new version parser */
+               cstr_info->comp = image->numcomps;
+               cstr_info->layer = (&cp->tcps[0])->numlayers;
+               cstr_info->decomposition = (&cp->tcps[0])->tccps->numresolutions - 1;
+               cstr_info->D_max = 0;           /* ADD Marcela */
+/* UniPG>> */
+               cstr_info->main_head_start = cio_tell(cio); /* position of SOC */
+/* <<UniPG */
+       }
+       else if (cstr_info) {
+               cstr_info->index_on = 0;
        }
        /* << INDEX */
        
        j2k_write_soc(j2k);
        j2k_write_siz(j2k);
-/* UniPG>> */
-#ifdef USE_JPWL
-       /** THIS CODE IS NOT USED */
-       //if(image_info && image_info->index_on && cp->epc_on)
-       //      j2k_write_epc(j2k);
-#endif /* USE_JPWL */
-/* <<UniPG */
        j2k_write_cod(j2k);
        j2k_write_qcd(j2k);
        
@@ -2413,24 +2275,30 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
        if (cp->comment != NULL) {
                j2k_write_com(j2k);
        }
-       /* INDEX >> */
-       if(image_info && image_info->index_on) {
-               image_info->main_head_end = cio_tell(cio) - 1;
-       }
-       /* << INDEX */
-
-       /*       TLM Marker*/
+       
+       j2k->totnum_tp = j2k_calculate_tp(cp,image->numcomps,image,j2k);
+       /* TLM Marker*/
        if(cp->cinema){
-               j2k->totnum_tp = j2k_generate_tlm(cp,image->numcomps,image);
                j2k_write_tlm(j2k);
+               if (cp->cinema == CINEMA4K_24) {
+                       j2k_write_poc(j2k);
+               }
        }
+
+       /* uncomment only for testing JPSEC marker writing */
+       /* j2k_write_sec(j2k); */
+
+       /* INDEX >> */
+       if(cstr_info && cstr_info->index_on) {
+               cstr_info->main_head_end = cio_tell(cio) - 1;
+       }
+       /* << INDEX */
        /**** Main Header ENDS here ***/
        
        /* create the tile encoder */
        tcd = tcd_create(j2k->cinfo);
 
        /* encode each tile */
-
        for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
                int pino;
                int tilepartno=0;
@@ -2439,7 +2307,8 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
                opj_event_msg(j2k->cinfo, EVT_INFO, "tile number %d / %d\n", tileno + 1, cp->tw * cp->th);
                
                j2k->curtileno = tileno;
-
+               j2k->cur_tp_num = 0;
+               tcd->cur_totnum_tp = j2k->cur_totnum_tp[j2k->curtileno];
                /* initialisation before tile encoding  */
                if (tileno == 0) {
                        tcd_malloc_encode(tcd, image, cp, j2k->curtileno);
@@ -2448,26 +2317,29 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
                }
                
                /* INDEX >> */
-               if(image_info && image_info->index_on) {
-                       image_info->tile[j2k->curtileno].num_tile = j2k->curtileno;
-                       image_info->tile[j2k->curtileno].start_pos = cio_tell(cio) + j2k->pos_correction;
+               if(cstr_info && cstr_info->index_on) {
+                       cstr_info->tile[j2k->curtileno].num_tile = j2k->curtileno;
+                       cstr_info->tile[j2k->curtileno].start_pos = cio_tell(cio) + j2k->pos_correction;
                }
                /* << INDEX */
 
                for(pino = 0; pino <= tcp->numpocs; pino++) {
                        int tot_num_tp;
-                       opj_pi_iterator_t *pi = NULL;
                        tcd->cur_pino=pino;
-                       tcd->pi = pi_initialise_encode(image, cp, tileno,pino);
-                       pi=tcd->pi;
                        
                        /*Get number of tile parts*/
                        tot_num_tp = j2k_get_num_tp(cp,pino,tileno);
                        tcd->tp_pos = cp->tp_pos;
 
                        for(tilepartno = 0; tilepartno < tot_num_tp ; tilepartno++){
-                               j2k->cur_tp_num = tilepartno;
-                               j2k->cur_totnum_tp = tot_num_tp;
+                               j2k->tp_num = tilepartno;
+/* UniPG>> */
+                               /* INDEX >> */
+                               if(cstr_info && cstr_info->index_on)
+                                       cstr_info->tile[j2k->curtileno].tp_start_pos[j2k->cur_tp_num] =
+                                               cio_tell(cio) + j2k->pos_correction;
+                               /* << INDEX */
+/* <<UniPG */                          
                                j2k_write_sot(j2k);
 
                                if(j2k->cur_tp_num == 0 && cp->cinema == 0){
@@ -2480,13 +2352,28 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
                                        }
                                }
 
+/* UniPG>> */
+                               /* INDEX >> */
+                               if(cstr_info && cstr_info->index_on)
+                                       cstr_info->tile[j2k->curtileno].tp_end_header[j2k->cur_tp_num] =
+                                               cio_tell(cio) + j2k->pos_correction + 1;
+                               /* << INDEX */
+/* <<UniPG */
                                j2k_write_sod(j2k, tcd);
+/* UniPG>> */
+                               /* INDEX >> */
+                               if(cstr_info && cstr_info->index_on)
+                                       cstr_info->tile[j2k->curtileno].tp_end_pos[j2k->cur_tp_num] =
+                                               cio_tell(cio) + j2k->pos_correction - 1;
+                               /* << INDEX */
+/* <<UniPG */
+                               j2k->cur_tp_num ++;
                        }
-                       pi_destroy(pi, cp, tileno);
+                       
                }
                /* INDEX >> */
-               if(image_info && image_info->index_on) {
-                       image_info->tile[j2k->curtileno].end_pos = cio_tell(cio) + j2k->pos_correction - 1;
+               if(cstr_info && cstr_info->index_on) {
+                       cstr_info->tile[j2k->curtileno].end_pos = cio_tell(cio) + j2k->pos_correction - 1;
                }
                /* << INDEX */
                
@@ -2513,23 +2400,26 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
        tcd_free_encode(tcd);
        tcd_destroy(tcd);
 
+       free(j2k->cur_totnum_tp);
+
        j2k_write_eoc(j2k);
-       
-       /* Creation of the index file */
-       if(image_info && image_info->index_on) {
-               if(!j2k_create_index(j2k, cio, image_info, index)) {
-                       opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to create index file %s\n", index);
-                       return false;
-               }
-       }
 
+       if(cstr_info && cstr_info->index_on) {
+               cstr_info->codestream_size = cio_tell(cio) + j2k->pos_correction;
 /* UniPG>> */
+               /* The following adjustment is done to adjust the codestream size */
+               /* if SOD is not at 0 in the buffer. Useful in case of JP2, where */
+               /* the first unch of bytes is not in the codestream               */
+               cstr_info->codestream_size -= cstr_info->main_head_start;
+/* <<UniPG */
+       }
+
 #ifdef USE_JPWL
        /*
        preparation of JPWL marker segments: can be finalized only when the whole
        codestream is known
        */
-       if(image_info && image_info->index_on && cp->epc_on) {
+       if(cstr_info && cstr_info->index_on && cp->epc_on) {
 
                /* let's begin creating a marker list, according to user wishes */
                jpwl_prepare_marks(j2k, cio, image);
@@ -2540,22 +2430,11 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
                /* do not know exactly what is this for,
                but it gets called during index creation */
                j2k->pos_correction = 0;
-
-               /* Re-creation of the index file, with updated info */
-               if(image_info && image_info->index_on) {
-                       if(!j2k_create_index(j2k, cio, image_info, index)) {
-                               opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to re-create index file %s\n", index);
-                               return false;
-                       }
-               }
-
-               /* now we finalize the marker contents */
-               /*jpwl_finalize_marks(j2k, cio, image);*/
-
        }
 #endif /* USE_JPWL */
-/* <<UniPG */
-         
+
        return true;
 }
 
+
+