Fixed the DCinema filesize allocation. It now includes the SOT marker size
[openjpeg.git] / libopenjpeg / j2k.c
index 1416684fc34ffb23dde07868f5ba7385ec1a3e81..8fabbf96fcbc5cf35075b73ea2795cfbbbd4dcaa 100644 (file)
@@ -5,6 +5,7 @@
  * Copyright (c) 2002-2003, Yannick Verschueren
  * Copyright (c) 2003-2007, Francois-Olivier Devaux and Antonin Descampe
  * Copyright (c) 2005, Herve Drolon, FreeImage Team
+ * Copyright (c) 2006-2007, Parvatha Elangovan
  * All rights reserved.
  *
  * Redistribution and use in source and binary forms, with or without
@@ -175,6 +176,11 @@ Read the PPT marker (packet packet headers, tile-part header)
 */
 static void j2k_read_ppt(opj_j2k_t *j2k);
 /**
+Write the TLM marker (Mainheader)
+@param j2k J2K handle
+*/
+static void j2k_write_tlm(opj_j2k_t *j2k);
+/**
 Write the SOT marker (start of tile-part)
 @param j2k J2K handle
 */
@@ -228,6 +234,66 @@ static void j2k_read_unk(opj_j2k_t *j2k);
 /*@}*/
 
 /* ----------------------------------------------------------------------- */
+typedef struct j2k_prog_order{
+       OPJ_PROG_ORDER enum_prog;
+       char str_prog[4];
+}j2k_prog_order_t;
+
+j2k_prog_order_t j2k_prog_order_list[] = {
+       {CPRL, "CPRL"},
+       {LRCP, "LRCP"},
+       {PCRL, "PCRL"},
+       {RLCP, "RLCP"},
+       {RPCL, "RPCL"},
+       {-1, ""}
+};
+
+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){
+                       break;
+               }
+       }
+       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;
@@ -303,6 +369,61 @@ void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_t * cp) {
 }
 
 /* ----------------------------------------------------------------------- */
+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 = j2k_convert_progression_order(tcp->prg);
+       
+       if(cp->tp_on == 1){
+               for(i=0;i<4;i++){
+                       if(tpend!=1){
+                               if( cp->tp_flag == prog[i] ){
+                                       tpend=1;cp->tp_pos=i;
+                               }
+                               switch(prog[i]){
+                               case 'C':
+                                       tpnum= tpnum * tcp->pocs[pino].compE;
+                                       break;
+                               case 'R':
+                                       tpnum= tpnum * tcp->pocs[pino].resE;
+                                       break;
+                               case 'P':
+                                       tpnum= tpnum * tcp->pocs[pino].prcE;
+                                       break;
+                               case 'L':
+                                       tpnum= tpnum * tcp->pocs[pino].layE;
+                                       break;
+                               }
+                       }
+               }
+       }else{
+               tpnum=1;
+       }
+       return tpnum;
+}
+
+/**    mem allocation for TLM marker*/
+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 = 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;
+       }
+       return totnum_tp;
+}
 
 static void j2k_write_soc(opj_j2k_t *j2k) {
        opj_cio_t *cio = j2k->cio;
@@ -324,7 +445,7 @@ static void j2k_write_siz(opj_j2k_t *j2k) {
        cio_write(cio, J2K_MS_SIZ, 2);  /* SIZ */
        lenp = cio_tell(cio);
        cio_skip(cio, 2);
-       cio_write(cio, 0, 2);                   /* Rsiz (capabilities) */
+       cio_write(cio, cp->rsiz, 2);                    /* Rsiz (capabilities) */
        cio_write(cio, image->x1, 4);   /* Xsiz */
        cio_write(cio, image->y1, 4);   /* Ysiz */
        cio_write(cio, image->x0, 4);   /* X0siz */
@@ -364,7 +485,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
@@ -373,15 +494,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)) {
@@ -399,7 +524,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;
@@ -409,7 +534,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
@@ -418,8 +542,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) {
@@ -436,7 +562,6 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                        
                }
 #endif /* USE_JPWL */
-/* <<UniPG */
 
 
                /* TODO: unused ? */
@@ -449,7 +574,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
@@ -458,8 +583,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) {
@@ -468,11 +595,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;
@@ -480,32 +606,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;
@@ -726,7 +852,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) {
 
@@ -735,18 +861,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) {
@@ -825,7 +953,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) {
 
@@ -836,20 +964,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));
 }
 
@@ -863,7 +993,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 */
@@ -901,8 +1031,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 */
@@ -1010,7 +1140,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) {
@@ -1018,12 +1148,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;
                }
@@ -1067,6 +1200,18 @@ static void j2k_read_ppt(opj_j2k_t *j2k) {
        tcp->ppt_store = j;
 }
 
+static void j2k_write_tlm(opj_j2k_t *j2k){
+       int lenp;
+       opj_cio_t *cio = j2k->cio;
+       j2k->tlm_start = cio_tell(cio);
+       cio_write(cio, J2K_MS_TLM, 2);/* TLM */
+       lenp = 4 + (5*j2k->totnum_tp);
+       cio_write(cio,lenp,2);                          /* Ltlm */
+       cio_write(cio, 0,1);                                    /* Ztlm=0*/
+       cio_write(cio,80,1);                                    /* Stlm ST=1(8bits-255 tiles max),SP=1(Ptlm=32bits) */
+       cio_skip(cio,5*j2k->totnum_tp);
+}
+
 static void j2k_write_sot(opj_j2k_t *j2k) {
        int lenp, len;
 
@@ -1078,8 +1223,8 @@ static void j2k_write_sot(opj_j2k_t *j2k) {
        cio_skip(cio, 2);                                       /* Lsot (further) */
        cio_write(cio, j2k->curtileno, 2);      /* Isot */
        cio_skip(cio, 4);                                       /* Psot (further in j2k_write_sod) */
-       cio_write(cio, 0, 1);                           /* TPsot */
-       cio_write(cio, 1, 1);                           /* TNsot */
+       cio_write(cio, j2k->cur_tp_num , 1);    /* TPsot */
+       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 */
@@ -1096,7 +1241,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) {
 
@@ -1107,20 +1252,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;
@@ -1138,7 +1285,7 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
        }
        
        totlen = cio_read(cio, 4);
-/* UniPG>> */
+
 #ifdef USE_JPWL
        if (j2k->cp->correct) {
 
@@ -1147,19 +1294,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;
        
@@ -1197,16 +1345,24 @@ static void j2k_write_sod(opj_j2k_t *j2k, void *tile_coder) {
        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;
        
        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;
+               if (!j2k->cur_tp_num ){
+                       image_info->tile[j2k->curtileno].end_header = cio_tell(cio) + j2k->pos_correction - 1;
+               }else{
+                       if(image_info->tile[j2k->curtileno].packet[image_info->num - 1].end_pos < cio_tell(cio))
+                               image_info->tile[j2k->curtileno].packet[image_info->num].start_pos = cio_tell(cio);
+               }
        }
        /* << INDEX */
        
@@ -1214,7 +1370,7 @@ 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) {
+       if(image_info && (j2k->cur_tp_num == 0)) {
                image_info->num = 0;
        }
        
@@ -1225,6 +1381,13 @@ static void j2k_write_sod(opj_j2k_t *j2k, void *tile_coder) {
        cio_seek(cio, j2k->sot_start + 6);
        cio_write(cio, totlen, 4);
        cio_seek(cio, j2k->sot_start + totlen);
+       /* Writing Ttlm and Ptlm in TLM marker */
+       if(cp->cinema){
+               cio_seek(cio, j2k->tlm_start + 6 + (5*j2k->cur_tp_num));
+               cio_write(cio, j2k->curtileno, 1);
+               cio_write(cio, totlen, 4);
+       }
+       cio_seek(cio, j2k->sot_start + totlen);
 }
 
 static void j2k_read_sod(opj_j2k_t *j2k) {
@@ -1287,7 +1450,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!!! */
@@ -1295,12 +1458,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 */
 }
 
@@ -1313,24 +1478,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;
 }
@@ -1364,20 +1534,20 @@ 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 */
+
   {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;
@@ -1389,8 +1559,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 */
@@ -1434,7 +1604,7 @@ static void j2k_read_unk(opj_j2k_t *j2k) {
 
        };
 #endif /* USE_JPWL */
-/* <<UniPG */
+
 }
 
 /**
@@ -1522,13 +1692,14 @@ void j2k_setup_decoder(opj_j2k_t *j2k, opj_dparameters_t *parameters) {
                opj_cp_t *cp = (opj_cp_t*)opj_malloc(sizeof(opj_cp_t));
                cp->reduce = parameters->cp_reduce;     
                cp->layer = parameters->cp_layer;
-/* UniPG>> */
+               cp->limit_decoding = parameters->cp_limit_decoding;
+
 #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;
@@ -1552,7 +1723,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) {
@@ -1571,34 +1742,44 @@ 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 */
+#endif /* USE_JPWL */
+
                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);
                        return 0;
                }
                e = j2k_dec_mstab_lookup(id);
+               // Check if the marker is known
                if (!(j2k->state & e->states)) {
                        opj_image_destroy(image);
                        opj_event_msg(cinfo, EVT_ERROR, "%.8x: unexpected marker %x\n", cio_tell(cio) - 2, id);
                        return 0;
                }
+               // Check if the decoding is limited to the main header
+               if (e->id == J2K_MS_SOT && j2k->cp->limit_decoding == LIMIT_TO_MAIN_HEADER) {
+                       opj_event_msg(cinfo, EVT_INFO, "Main Header decoded.\n");
+                       return image;
+               }               
+
                if (e->handler) {
                        (*e->handler)(j2k);
                }
@@ -1772,7 +1953,9 @@ 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;
        cp->fixed_quality = parameters->cp_fixed_quality;
@@ -1818,7 +2001,17 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                cp->tdy = image->y1 - cp->ty0;
        }
 
-/* UniPG>> */
+       if(parameters->tp_on){
+               cp->tp_flag = parameters->tp_flag;
+               cp->tp_on = 1;
+       }
+       
+       cp->img_size = 0;
+       for(i=0;i<image->numcomps ;i++){
+       cp->img_size += (image->comps[i].w *image->comps[i].h * image->comps[i].prec);
+       }
+
+
 #ifdef USE_JPWL
        /*
        calculate JPWL encoding parameters
@@ -1874,7 +2067,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 */
        /* ---------------------------- */
@@ -1884,21 +2077,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];
@@ -1907,13 +2108,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));
                
@@ -1932,42 +2135,60 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                        } else {
                                tccp->roishift = 0;
                        }
-                       if (parameters->csty & J2K_CCP_CSTY_PRT) {
-                               int p = 0;
-                               for (j = tccp->numresolutions - 1; j >= 0; j--) {
-                                       if (p < parameters->res_spec) {
-                                               if (parameters->prcw_init[p] < 1) {
-                                                       tccp->prcw[j] = 1;
-                                               } else {
-                                                       tccp->prcw[j] = int_floorlog2(parameters->prcw_init[p]);
-                                               }
-                                               if (parameters->prch_init[p] < 1) {
-                                                       tccp->prch[j] = 1;
-                                               } else {
-                                                       tccp->prch[j] = int_floorlog2(parameters->prch_init[p]);
-                                               }
-                                       } else {
-                                               int res_spec = parameters->res_spec;
-                                               int size_prcw = parameters->prcw_init[res_spec - 1] >> (p - (res_spec - 1));
-                                               int size_prch = parameters->prch_init[res_spec - 1] >> (p - (res_spec - 1));
-                                               if (size_prcw < 1) {
-                                                       tccp->prcw[j] = 1;
-                                               } else {
-                                                       tccp->prcw[j] = int_floorlog2(size_prcw);
-                                               }
-                                               if (size_prch < 1) {
-                                                       tccp->prch[j] = 1;
+
+                       if(parameters->cp_cinema)
+                       {
+                               //Precinct size for lowest frequency subband=128
+                               tccp->prcw[0] = 7;
+                               tccp->prch[0] = 7;
+                               //Precinct size at all other resolutions = 256
+                               for (j = 1; j < tccp->numresolutions; j++) {
+                                       tccp->prcw[j] = 8;
+                                       tccp->prch[j] = 8;
+                               }
+                       }else{
+                               if (parameters->csty & J2K_CCP_CSTY_PRT) {
+                                       int p = 0;
+                                       for (j = tccp->numresolutions - 1; j >= 0; j--) {
+                                               if (p < parameters->res_spec) {
+                                                       
+                                                       if (parameters->prcw_init[p] < 1) {
+                                                               tccp->prcw[j] = 1;
+                                                       } else {
+                                                               tccp->prcw[j] = int_floorlog2(parameters->prcw_init[p]);
+                                                       }
+                                                       
+                                                       if (parameters->prch_init[p] < 1) {
+                                                               tccp->prch[j] = 1;
+                                                       }else {
+                                                               tccp->prch[j] = int_floorlog2(parameters->prch_init[p]);
+                                                       }
+
                                                } else {
-                                                       tccp->prch[j] = int_floorlog2(size_prch);
+                                                       int res_spec = parameters->res_spec;
+                                                       int size_prcw = parameters->prcw_init[res_spec - 1] >> (p - (res_spec - 1));
+                                                       int size_prch = parameters->prch_init[res_spec - 1] >> (p - (res_spec - 1));
+                                                       
+                                                       if (size_prcw < 1) {
+                                                               tccp->prcw[j] = 1;
+                                                       } else {
+                                                               tccp->prcw[j] = int_floorlog2(size_prcw);
+                                                       }
+                                                       
+                                                       if (size_prch < 1) {
+                                                               tccp->prch[j] = 1;
+                                                       } else {
+                                                               tccp->prch[j] = int_floorlog2(size_prch);
+                                                       }
                                                }
+                                               p++;
+                                               /*printf("\nsize precinct for level %d : %d,%d\n", j,tccp->prcw[j], tccp->prch[j]); */
+                                       }       //end for
+                               } else {
+                                       for (j = 0; j < tccp->numresolutions; j++) {
+                                               tccp->prcw[j] = 15;
+                                               tccp->prch[j] = 15;
                                        }
-                                       p++;
-                                       /*printf("\nsize precinct for level %d : %d,%d\n", j,tccp->prcw[j], tccp->prch[j]); */
-                               }
-                       } else {
-                               for (j = 0; j < tccp->numresolutions; j++) {
-                                       tccp->prcw[j] = 15;
-                                       tccp->prch[j] = 15;
                                }
                        }
 
@@ -1991,7 +2212,7 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
 
        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
@@ -2000,7 +2221,7 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
        if (j2k->cp->epc_on && !strcmp(index, JPWL_PRIVATEINDEX_NAME))
                return 1;
 #endif /* USE_JPWL */
-/* <<UniPG */
+
 
        stream = fopen(index, "w");
        if (!stream) {
@@ -2038,15 +2259,11 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
                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");
-                       */
+
+                       fprintf(stream, "pack_nb tileno layno resno compno precno start_pos  end_pos disto\n");
+
                        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++) {
@@ -2065,9 +2282,9 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
                        }
                } /* LRCP */
                else if (image_info->prog == RLCP) {    /* RLCP */
-                       /*
-                       fprintf(stream, "pack_nb tileno resno layno compno precno start_pos  end_pos   disto");
-                       */
+
+                       fprintf(stream, "pack_nb tileno resno layno compno precno start_pos  end_pos disto\n");
+
                        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++) {
@@ -2086,26 +2303,26 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
                        }
                } /* RLCP */
                else if (image_info->prog == RPCL) {    /* RPCL */
-                       /*
-                       fprintf(stream, "\npack_nb tileno resno precno compno layno start_pos  end_pos   disto\n"); 
-                       */
+
+                       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 ) {
+                               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 );
+                                               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 < image_info->layer; layno++) {
                                                                                        start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
@@ -2117,11 +2334,11 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
                                                                                        pack_nb++; 
                                                                                }
                                                                        }
-                                                               }
-                                                       } /* precno */
-                                               } /* compno */
-                                       } /* x = x0..x1 */
-                               } /* y = y0..y1 */
+                                                               }/* x = x0..x1 */
+                                                       } 
+                                               }  /* y = y0..y1 */
+                                       } /* precno */
+                               } /* compno */
                        } /* resno */
                } /* RPCL */
                else if (image_info->prog == PCRL) {    /* PCRL */
@@ -2130,21 +2347,21 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
                        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 ) {
+
+                       fprintf(stream, "\npack_nb tileno precno compno resno layno start_pos  end_pos disto\n"); 
+
+                       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 );
+                                               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 < image_info->layer; layno++) {
                                                                                        start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
@@ -2156,34 +2373,35 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
                                                                                        pack_nb++; 
                                                                                }
                                                                        }
-                                                               }
-                                                       } /* precno */
-                                               } /* resno */
-                                       } /* compno */
-                               } /* x = x0..x1 */
-                       } /* y = y0..y1 */
+                                                               }/* x = x0..x1 */
+                                                       } 
+                                               }  /* y = y0..y1 */
+                                       } /* precno */
+                               } /* resno */
+                       } /* compno */
                } /* PCRL */
                else {  /* CPRL */
-                       /*
-                       fprintf(stream, "\npack_nb tileno compno precno resno layno start_pos  end_pos   disto\n"); 
-                       */
+
+                       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 ) {
+                               
+                               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 );
+                                               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 < image_info->layer; layno++) {
                                                                                        start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
@@ -2195,12 +2413,12 @@ static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *im
                                                                                        pack_nb++; 
                                                                                }
                                                                        }
-                                                               }
-                                                       } /* precno */
-                                               } /* resno */
-                                       } /* x = x0..x1 */
-                               } /* y = y0..y1 */
-                       } /* comno */
+                                                               }/* x = x0..x1 */
+                                                       }
+                                               } /* y = y0..y1 */
+                                       } /* precno */
+                               } /* resno */
+                       } /* compno */
                } /* CPRL */   
        } /* tileno */
        
@@ -2248,15 +2466,16 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *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);
+       
+       if(cp->cinema){
+               for (compno = 1; compno < image->numcomps; compno++) {
+                       j2k_write_coc(j2k, compno);
+                       j2k_write_qcc(j2k, compno);
+               }
+       }
+
        for (compno = 0; compno < image->numcomps; compno++) {
                opj_tcp_t *tcp = &cp->tcps[0];
                if (tcp->tccps[compno].roishift)
@@ -2271,16 +2490,31 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
        }
        /* << INDEX */
 
+       j2k->totnum_tp = j2k_calculate_tp(cp,image->numcomps,image,j2k);
+       /*       TLM Marker*/
+       if(cp->cinema){
+               j2k_write_tlm(j2k);
+               if (cp->cinema == CINEMA4K_24) {
+                       j2k_write_poc(j2k);
+               }
+       }
+       /**** 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;
+
+               opj_tcp_t *tcp = &cp->tcps[tileno];
                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);
@@ -2295,18 +2529,33 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
                }
                /* << INDEX */
 
-               j2k_write_sot(j2k);
-               
-               for (compno = 1; compno < image->numcomps; compno++) {
-                       j2k_write_coc(j2k, compno);
-                       j2k_write_qcc(j2k, compno);
-               }
-               if (cp->tcps[tileno].numpocs) {
-                       j2k_write_poc(j2k);
-               }
+               for(pino = 0; pino <= tcp->numpocs; pino++) {
+                       int tot_num_tp;
+                       tcd->cur_pino=pino;
+                       
+                       /*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->tp_num = tilepartno;
+                               j2k_write_sot(j2k);
+
+                               if(j2k->cur_tp_num == 0 && cp->cinema == 0){
+                                       for (compno = 1; compno < image->numcomps; compno++) {
+                                               j2k_write_coc(j2k, compno);
+                                               j2k_write_qcc(j2k, compno);
+                                       }
+                                       if (cp->tcps[tileno].numpocs) {
+                                               j2k_write_poc(j2k);
+                                       }
+                               }
 
-               j2k_write_sod(j2k, tcd);
-               
+                               j2k_write_sod(j2k, tcd);
+                               j2k->cur_tp_num ++;
+                       }
+                       
+               }
                /* INDEX >> */
                if(image_info && image_info->index_on) {
                        image_info->tile[j2k->curtileno].end_pos = cio_tell(cio) + j2k->pos_correction - 1;
@@ -2346,7 +2595,7 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
                }
        }
 
-/* UniPG>> */
+
 #ifdef USE_JPWL
        /*
        preparation of JPWL marker segments: can be finalized only when the whole
@@ -2377,8 +2626,9 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
 
        }
 #endif /* USE_JPWL */
-/* <<UniPG */
-         
+
+       
        return true;
 }
 
+