Initial revision

1 files changed, 1769 insertions, 0 deletions

diff --git a/jpwl/decoder/libopenjpeg/j2k.c b/jpwl/decoder/libopenjpeg/j2k.c
new file mode 100644
index 00000000..7360a108
--- /dev/null
+++ b/jpwl/decoder/libopenjpeg/j2k.c
@@ -0,0 +1,1769 @@
+	// Questo file per la versione jpwldec1.6 è modificato con l'aggiunta della funzione per 
+	// lettura del marker segment ESD!!!
+/*
+ * Copyright (c) 2001-2002, David Janssens
+ * Copyright (c) 2002-2004, Yannick Verschueren
+ * Copyright (c) 2002-2004, Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium
+ * All rights reserved.
+ *
+ * Redistribution and use in source and binary forms, with or without
+ * modification, are permitted provided that the following conditions
+ * are met:
+ * 1. Redistributions of source code must retain the above copyright
+ *    notice, this list of conditions and the following disclaimer.
+ * 2. Redistributions in binary form must reproduce the above copyright
+ *    notice, this list of conditions and the following disclaimer in the
+ *    documentation and/or other materials provided with the distribution.
+ *
+ * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
+ * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+ * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+ * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+ * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+ * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+ * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+ * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+ * POSSIBILITY OF SUCH DAMAGE.
+ */
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <setjmp.h>
+#include <math.h>
+
+#include "j2k.h"
+#include "cio.h"
+#include "tcd.h"
+#include "dwt.h"
+#include "int.h"
+#include "jpt.h"
+#include "jpw.h"
+
+#define J2K_MS_SOC 0xff4f
+#define J2K_MS_SOT 0xff90
+#define J2K_MS_SOD 0xff93
+#define J2K_MS_EOC 0xffd9
+#define J2K_MS_SIZ 0xff51
+#define J2K_MS_COD 0xff52
+#define J2K_MS_COC 0xff53
+#define J2K_MS_RGN 0xff5e
+#define J2K_MS_QCD 0xff5c
+#define J2K_MS_QCC 0xff5d
+#define J2K_MS_POC 0xff5f
+#define J2K_MS_TLM 0xff55
+#define J2K_MS_PLM 0xff57
+#define J2K_MS_PLT 0xff58
+#define J2K_MS_PPM 0xff60
+#define J2K_MS_PPT 0xff61
+#define J2K_MS_SOP 0xff91
+#define J2K_MS_EPH 0xff92
+#define J2K_MS_CRG 0xff63
+#define J2K_MS_COM 0xff64
+
+#define JPWL_MS_EPB 0xff96
+#define JPWL_MS_EPC 0xff97
+#define JPWL_MS_ESD 0xff98 // marker segment ESD: aggiunto per JPWL
+#define JPWL_MS_RED 0xff99
+
+#define J2K_STATE_MHSOC 0x0001
+#define J2K_STATE_MHSIZ 0x0002
+#define J2K_STATE_MH 0x0004
+#define J2K_STATE_TPHSOT 0x0008
+#define J2K_STATE_TPH 0x0010
+#define J2K_STATE_MT 0x0020
+#define J2K_STATE_NEOC 0x0040
+
+
+jmp_buf j2k_error;
+
+static int j2k_state;
+static int j2k_curtileno;
+static j2k_tcp_t j2k_default_tcp;
+static unsigned char *j2k_eot;
+static int j2k_sot_start;
+static int pos_correction;
+
+static j2k_image_t *j2k_img;
+static j2k_cp_t *j2k_cp;
+
+static unsigned char **j2k_tile_data;
+static int *j2k_tile_len;
+
+static info_image info_IM;
+
+static int nesd; // Aggiunto per JPWL (gestione ESD!)
+extern int cslen; // Aggiunto per JPWL (gestione RED!)
+
+/* Add Patrick */
+void j2k_clean()
+{
+  int tileno = 0;
+  tcd_free_encode(j2k_img, j2k_cp, j2k_curtileno);
+
+  if (info_IM.index_on) {
+    for (tileno = 0; tileno < j2k_cp->tw * j2k_cp->th; tileno++) {
+      free(info_IM.tile[tileno].packet);
+    }
+    free(info_IM.tile);
+  }
+}
+
+/* \Add Patrick */
+
+void j2k_dump_image(j2k_image_t * img)
+{
+  int compno;
+  fprintf(stderr, "image {\n");
+  fprintf(stderr, "  x0=%d, y0=%d, x1=%d, y1=%d\n", img->x0, img->y0,
+	  img->x1, img->y1);
+  fprintf(stderr, "  numcomps=%d\n", img->numcomps);
+  for (compno = 0; compno < img->numcomps; compno++) {
+    j2k_comp_t *comp = &img->comps[compno];
+    fprintf(stderr, "  comp %d {\n", compno);
+    fprintf(stderr, "    dx=%d, dy=%d\n", comp->dx, comp->dy);
+    fprintf(stderr, "    prec=%d\n", comp->prec);
+    fprintf(stderr, "    sgnd=%d\n", comp->sgnd);
+    fprintf(stderr, "  }\n");
+  }
+  fprintf(stderr, "}\n");
+}
+
+void j2k_dump_cp(j2k_image_t * img, j2k_cp_t * cp)
+{
+  int tileno, compno, layno, bandno, resno, numbands;
+  fprintf(stderr, "coding parameters {\n");
+  fprintf(stderr, "  tx0=%d, ty0=%d\n", cp->tx0, cp->ty0);
+  fprintf(stderr, "  tdx=%d, tdy=%d\n", cp->tdx, cp->tdy);
+  fprintf(stderr, "  tw=%d, th=%d\n", cp->tw, cp->th);
+  for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
+    j2k_tcp_t *tcp = &cp->tcps[tileno];
+    fprintf(stderr, "  tile %d {\n", tileno);
+    fprintf(stderr, "    csty=%x\n", tcp->csty);
+    fprintf(stderr, "    prg=%d\n", tcp->prg);
+    fprintf(stderr, "    numlayers=%d\n", tcp->numlayers);
+    fprintf(stderr, "    mct=%d\n", tcp->mct);
+    fprintf(stderr, "    rates=");
+    for (layno = 0; layno < tcp->numlayers; layno++) {
+      fprintf(stderr, "%d ", tcp->rates[layno]);
+    }
+    fprintf(stderr, "\n");
+    for (compno = 0; compno < img->numcomps; compno++) {
+      j2k_tccp_t *tccp = &tcp->tccps[compno];
+      fprintf(stderr, "    comp %d {\n", compno);
+      fprintf(stderr, "      csty=%x\n", tccp->csty);
+      fprintf(stderr, "      numresolutions=%d\n", tccp->numresolutions);
+      fprintf(stderr, "      cblkw=%d\n", tccp->cblkw);
+      fprintf(stderr, "      cblkh=%d\n", tccp->cblkh);
+      fprintf(stderr, "      cblksty=%x\n", tccp->cblksty);
+      fprintf(stderr, "      qmfbid=%d\n", tccp->qmfbid);
+      fprintf(stderr, "      qntsty=%d\n", tccp->qntsty);
+      fprintf(stderr, "      numgbits=%d\n", tccp->numgbits);
+      fprintf(stderr, "      roishift=%d\n", tccp->roishift);
+      fprintf(stderr, "      stepsizes=");
+      numbands =
+	tccp->qntsty ==
+	J2K_CCP_QNTSTY_SIQNT ? 1 : tccp->numresolutions * 3 - 2;
+      for (bandno = 0; bandno < numbands; bandno++) {
+	fprintf(stderr, "(%d,%d) ", tccp->stepsizes[bandno].mant,
+		tccp->stepsizes[bandno].expn);
+      }
+      fprintf(stderr, "\n");
+
+      if (tccp->csty & J2K_CCP_CSTY_PRT) {
+	fprintf(stderr, "      prcw=");
+	for (resno = 0; resno < tccp->numresolutions; resno++) {
+	  fprintf(stderr, "%d ", tccp->prcw[resno]);
+	}
+	fprintf(stderr, "\n");
+	fprintf(stderr, "      prch=");
+	for (resno = 0; resno < tccp->numresolutions; resno++) {
+	  fprintf(stderr, "%d ", tccp->prch[resno]);
+	}
+	fprintf(stderr, "\n");
+      }
+      fprintf(stderr, "    }\n");
+    }
+    fprintf(stderr, "  }\n");
+  }
+  fprintf(stderr, "}\n");
+}
+
+void j2k_write_soc()
+{
+  cio_write(J2K_MS_SOC, 2);
+}
+
+void j2k_read_soc()
+{
+  j2k_state = J2K_STATE_MHSIZ;
+}
+
+void j2k_write_siz()
+{
+  int i;
+  int lenp, len;
+
+  cio_write(J2K_MS_SIZ, 2);	/* SIZ                 */
+  lenp = cio_tell();
+  cio_skip(2);
+  cio_write(0, 2);		/* Rsiz (capabilities) */
+  cio_write(j2k_img->x1, 4);	/* Xsiz                */
+  cio_write(j2k_img->y1, 4);	/* Ysiz                */
+  cio_write(j2k_img->x0, 4);	/* X0siz               */
+  cio_write(j2k_img->y0, 4);	/* Y0siz               */
+  cio_write(j2k_cp->tdx, 4);	/* XTsiz               */
+  cio_write(j2k_cp->tdy, 4);	/* YTsiz               */
+  cio_write(j2k_cp->tx0, 4);	/* XT0siz              */
+  cio_write(j2k_cp->ty0, 4);	/* YT0siz              */
+  cio_write(j2k_img->numcomps, 2);	/* Csiz                */
+  for (i = 0; i < j2k_img->numcomps; i++) {
+    cio_write(j2k_img->comps[i].prec - 1 + (j2k_img->comps[i].sgnd << 7), 1);	/* Ssiz_i */
+    cio_write(j2k_img->comps[i].dx, 1);	/* XRsiz_i             */
+    cio_write(j2k_img->comps[i].dy, 1);	/* YRsiz_i             */
+  }
+  len = cio_tell() - lenp;
+  cio_seek(lenp);
+  cio_write(len, 2);		/* Lsiz                */
+  cio_seek(lenp + len);
+
+}
+
+void j2k_read_siz()
+{
+  int len, i;
+
+  len = cio_read(2);		/* Lsiz                */
+  cio_read(2);			/* Rsiz (capabilities) */
+  j2k_img->x1 = cio_read(4);	/* Xsiz                */
+  j2k_img->y1 = cio_read(4);	/* Ysiz                */
+  j2k_img->x0 = cio_read(4);	/* X0siz               */
+  j2k_img->y0 = cio_read(4);	/* Y0siz               */
+  j2k_cp->tdx = cio_read(4);	/* XTsiz               */
+  j2k_cp->tdy = cio_read(4);	/* YTsiz               */
+  j2k_cp->tx0 = cio_read(4);	/* XT0siz              */
+  j2k_cp->ty0 = cio_read(4);	/* YT0siz              */
+
+  j2k_img->numcomps = cio_read(2);	/* Csiz                */
+  j2k_img->comps =
+    (j2k_comp_t *) malloc(j2k_img->numcomps * sizeof(j2k_comp_t));
+  for (i = 0; i < j2k_img->numcomps; i++) {
+    int tmp, w, h;
+    tmp = cio_read(1);		/* Ssiz_i          */
+    j2k_img->comps[i].prec = (tmp & 0x7f) + 1;
+    j2k_img->comps[i].sgnd = tmp >> 7;
+    j2k_img->comps[i].dx = cio_read(1);	/* XRsiz_i         */
+    j2k_img->comps[i].dy = cio_read(1);	/* YRsiz_i         */
+    w = int_ceildiv(j2k_img->x1 - j2k_img->x0, j2k_img->comps[i].dx);
+    h = int_ceildiv(j2k_img->y1 - j2k_img->y0, j2k_img->comps[i].dy);
+    j2k_img->comps[i].resno_decoded = 0;	/* number of resolution decoded */
+    j2k_img->comps[i].factor = 0;	/* reducing factor by component */
+  }
+
+  j2k_cp->tw = int_ceildiv(j2k_img->x1 - j2k_cp->tx0, j2k_cp->tdx);
+  j2k_cp->th = int_ceildiv(j2k_img->y1 - j2k_cp->ty0, j2k_cp->tdy);
+  j2k_cp->tcps =
+    (j2k_tcp_t *) calloc(j2k_cp->tw * j2k_cp->th, sizeof(j2k_tcp_t));
+  j2k_cp->tileno = (int *) calloc(j2k_cp->tw * j2k_cp->th, sizeof(int));
+  j2k_cp->tileno_size = 0;
+
+  for (i = 0; i < j2k_cp->tw * j2k_cp->th; i++) {
+    j2k_cp->tcps[i].POC = 0;
+    j2k_cp->tcps[i].numpocs = 0;
+    j2k_cp->tcps[i].first = 1;
+  }
+
+  /* Initialization for PPM marker */
+  j2k_cp->ppm = 0;
+  j2k_cp->ppm_data = NULL;
+  j2k_cp->ppm_previous = 0;
+  j2k_cp->ppm_store = 0;
+
+  j2k_default_tcp.tccps =
+    (j2k_tccp_t *) calloc(sizeof(j2k_tccp_t), j2k_img->numcomps);
+  for (i = 0; i < j2k_cp->tw * j2k_cp->th; i++) {
+    j2k_cp->tcps[i].tccps =
+      (j2k_tccp_t *) calloc(sizeof(j2k_tccp_t), j2k_img->numcomps);
+  }
+  j2k_tile_data =
+    (unsigned char **) calloc(j2k_cp->tw * j2k_cp->th, sizeof(char *));
+  j2k_tile_len = (int *) calloc(j2k_cp->tw * j2k_cp->th, sizeof(int));
+  j2k_state = J2K_STATE_MH;
+
+
+}
+
+void j2k_write_com()
+{
+  unsigned int i;
+  int lenp, len;
+  char str[256];
+  sprintf(str, "%s", j2k_cp->comment);
+
+  cio_write(J2K_MS_COM, 2);
+  lenp = cio_tell();
+  cio_skip(2);
+  cio_write(0, 2);
+  for (i = 0; i < strlen(str); i++) {
+    cio_write(str[i], 1);
+  }
+  len = cio_tell() - lenp;
+  cio_seek(lenp);
+  cio_write(len, 2);
+  cio_seek(lenp + len);
+
+}
+
+void j2k_read_com()
+{
+  int len;
+
+  len = cio_read(2);
+  cio_skip(len - 2);
+
+}
+
+void j2k_write_cox(int compno)
+{
+  int i;
+  j2k_tcp_t *tcp;
+  j2k_tccp_t *tccp;
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+  tccp = &tcp->tccps[compno];
+
+  cio_write(tccp->numresolutions - 1, 1);	/* SPcox (D) */
+  cio_write(tccp->cblkw - 2, 1);	/* SPcox (E) */
+  cio_write(tccp->cblkh - 2, 1);	/* SPcox (F) */
+  cio_write(tccp->cblksty, 1);	/* SPcox (G) */
+  cio_write(tccp->qmfbid, 1);	/* SPcox (H) */
+
+  if (tccp->csty & J2K_CCP_CSTY_PRT) {
+    for (i = 0; i < tccp->numresolutions; i++) {
+      cio_write(tccp->prcw[i] + (tccp->prch[i] << 4), 1);	/* SPcox (I_i) */
+    }
+  }
+}
+
+void j2k_read_cox(int compno)
+{
+  int i;
+  j2k_tcp_t *tcp;
+  j2k_tccp_t *tccp;
+  tcp =
+    j2k_state ==
+    J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
+  tccp = &tcp->tccps[compno];
+  tccp->numresolutions = cio_read(1) + 1;	/* SPcox (D) */
+  tccp->cblkw = cio_read(1) + 2;	/* SPcox (E) */
+  tccp->cblkh = cio_read(1) + 2;	/* SPcox (F) */
+  tccp->cblksty = cio_read(1);	/* SPcox (G) */
+  tccp->qmfbid = cio_read(1);	/* SPcox (H) */
+  if (tccp->csty & J2K_CP_CSTY_PRT) {
+    for (i = 0; i < tccp->numresolutions; i++) {
+      int tmp = cio_read(1);	/* SPcox (I_i) */
+      tccp->prcw[i] = tmp & 0xf;
+      tccp->prch[i] = tmp >> 4;
+    }
+  }
+}
+
+void j2k_write_cod()
+{
+  j2k_tcp_t *tcp;
+  int lenp, len;
+
+  cio_write(J2K_MS_COD, 2);	/* COD */
+
+  lenp = cio_tell();
+  cio_skip(2);
+
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+  cio_write(tcp->csty, 1);	/* Scod */
+  cio_write(tcp->prg, 1);	/* SGcod (A) */
+  cio_write(tcp->numlayers, 2);	/* SGcod (B) */
+  cio_write(tcp->mct, 1);	/* SGcod (C) */
+
+  j2k_write_cox(0);
+  len = cio_tell() - lenp;
+  cio_seek(lenp);
+  cio_write(len, 2);		/* Lcod */
+  cio_seek(lenp + len);
+}
+
+void j2k_read_cod()
+{
+  int len, i, pos;
+  j2k_tcp_t *tcp;
+
+  tcp =
+    j2k_state ==
+    J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
+  len = cio_read(2);		/* Lcod */
+  tcp->csty = cio_read(1);	/* Scod */
+  tcp->prg = cio_read(1);	/* SGcod (A) */
+  tcp->numlayers = cio_read(2);	/* SGcod (B) */
+  tcp->mct = cio_read(1);	/* SGcod (C) */
+
+  pos = cio_tell();
+  for (i = 0; i < j2k_img->numcomps; i++) {
+    tcp->tccps[i].csty = tcp->csty & J2K_CP_CSTY_PRT;
+    cio_seek(pos);
+    j2k_read_cox(i);
+  }
+}
+
+void j2k_write_coc(int compno)
+{
+  j2k_tcp_t *tcp;
+  int lenp, len;
+
+  cio_write(J2K_MS_COC, 2);	/* COC */
+  lenp = cio_tell();
+  cio_skip(2);
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+  cio_write(compno, j2k_img->numcomps <= 256 ? 1 : 2);	/* Ccoc */
+  cio_write(tcp->tccps[compno].csty, 1);	/* Scoc */
+  j2k_write_cox(compno);
+  len = cio_tell() - lenp;
+  cio_seek(lenp);
+  cio_write(len, 2);		/* Lcoc */
+  cio_seek(lenp + len);
+}
+
+void j2k_read_coc()
+{
+  int len, compno;
+  j2k_tcp_t *tcp;
+
+  tcp =
+    j2k_state ==
+    J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
+  len = cio_read(2);		/* Lcoc */
+  compno = cio_read(j2k_img->numcomps <= 256 ? 1 : 2);	/* Ccoc */
+  tcp->tccps[compno].csty = cio_read(1);	/* Scoc */
+  j2k_read_cox(compno);
+}
+
+void j2k_write_qcx(int compno)
+{
+  j2k_tcp_t *tcp;
+  j2k_tccp_t *tccp;
+  int bandno, numbands;
+  int expn, mant;
+
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+  tccp = &tcp->tccps[compno];
+
+  cio_write(tccp->qntsty + (tccp->numgbits << 5), 1);	/* Sqcx */
+  numbands =
+    tccp->qntsty ==
+    J2K_CCP_QNTSTY_SIQNT ? 1 : tccp->numresolutions * 3 - 2;
+
+  for (bandno = 0; bandno < numbands; bandno++) {
+    expn = tccp->stepsizes[bandno].expn;
+    mant = tccp->stepsizes[bandno].mant;
+
+    if (tccp->qntsty == J2K_CCP_QNTSTY_NOQNT) {
+      cio_write(expn << 3, 1);	/* SPqcx_i */
+    } else {
+      cio_write((expn << 11) + mant, 2);	/* SPqcx_i */
+    }
+  }
+
+}
+
+void j2k_read_qcx(int compno, int len)
+{
+  int tmp;
+  j2k_tcp_t *tcp;
+  j2k_tccp_t *tccp;
+  int bandno, numbands;
+
+  tcp =
+    j2k_state ==
+    J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
+  tccp = &tcp->tccps[compno];
+  tmp = cio_read(1);		/* Sqcx */
+  tccp->qntsty = tmp & 0x1f;
+  tccp->numgbits = tmp >> 5;
+  numbands =
+    tccp->qntsty == J2K_CCP_QNTSTY_SIQNT ? 1 : (tccp->qntsty ==
+						J2K_CCP_QNTSTY_NOQNT ?
+						len - 1 : (len - 1) / 2);
+  for (bandno = 0; bandno < numbands; bandno++) {
+    int expn, mant;
+    if (tccp->qntsty == J2K_CCP_QNTSTY_NOQNT) {
+      expn = cio_read(1) >> 3;	/* SPqcx_i */
+      mant = 0;
+    } else {
+      tmp = cio_read(2);	/* SPqcx_i */
+      expn = tmp >> 11;
+      mant = tmp & 0x7ff;
+    }
+    tccp->stepsizes[bandno].expn = expn;
+    tccp->stepsizes[bandno].mant = mant;
+  }
+
+
+
+  /* Add Antonin : if scalar_derived -> compute other stepsizes */
+
+
+
+  if (tccp->qntsty==J2K_CCP_QNTSTY_SIQNT) {
+
+    for (bandno=1 ; bandno<J2K_MAXBANDS ; bandno++) {
+
+      tccp->stepsizes[bandno].expn = ((tccp->stepsizes[0].expn)-((bandno-1)/3)>0)?(tccp->stepsizes[0].expn)-((bandno-1)/3):0;
+
+      tccp->stepsizes[bandno].mant = tccp->stepsizes[0].mant;
+
+    }
+
+  }
+
+
+
+  /* ddA */
+}
+
+void j2k_write_qcd()
+{
+  int lenp, len;
+
+  cio_write(J2K_MS_QCD, 2);	/* QCD */
+  lenp = cio_tell();
+  cio_skip(2);
+  j2k_write_qcx(0);
+  len = cio_tell() - lenp;
+  cio_seek(lenp);
+  cio_write(len, 2);		/* Lqcd */
+  cio_seek(lenp + len);
+}
+
+void j2k_read_qcd()
+{
+  int len, i, pos;
+
+  len = cio_read(2);		/* Lqcd */
+  pos = cio_tell();
+  for (i = 0; i < j2k_img->numcomps; i++) {
+    cio_seek(pos);
+    j2k_read_qcx(i, len - 2);
+  }
+}
+
+void j2k_write_qcc(int compno)
+{
+  int lenp, len;
+
+  cio_write(J2K_MS_QCC, 2);	/* QCC */
+  lenp = cio_tell();
+  cio_skip(2);
+  cio_write(compno, j2k_img->numcomps <= 256 ? 1 : 2);	/* Cqcc */
+  j2k_write_qcx(compno);
+  len = cio_tell() - lenp;
+  cio_seek(lenp);
+  cio_write(len, 2);		/* Lqcc */
+  cio_seek(lenp + len);
+}
+
+void j2k_read_qcc()
+{
+  int len, compno;
+
+  len = cio_read(2);		/* Lqcc */
+  compno = cio_read(j2k_img->numcomps <= 256 ? 1 : 2);	/* Cqcc */
+  j2k_read_qcx(compno, len - 2 - (j2k_img->numcomps <= 256 ? 1 : 2));
+}
+
+void j2k_write_poc()
+{
+  int len, numpchgs, i;
+  j2k_tcp_t *tcp;
+  j2k_tccp_t *tccp;
+
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+  tccp = &tcp->tccps[0];
+  numpchgs = tcp->numpocs;
+  cio_write(J2K_MS_POC, 2);	/* POC  */
+  len = 2 + (5 + 2 * (j2k_img->numcomps <= 256 ? 1 : 2)) * numpchgs;
+  cio_write(len, 2);		/* Lpoc */
+  for (i = 0; i < numpchgs; i++) {
+    // MODIF
+    j2k_poc_t *poc;
+    poc = &tcp->pocs[i];
+    cio_write(poc->resno0, 1);	/* RSpoc_i */
+    cio_write(poc->compno0, (j2k_img->numcomps <= 256 ? 1 : 2));	/* CSpoc_i */
+    cio_write(poc->layno1, 2);	/* LYEpoc_i */
+    poc->layno1 = int_min(poc->layno1, tcp->numlayers);
+    cio_write(poc->resno1, 1);	/* REpoc_i */
+    poc->resno1 = int_min(poc->resno1, tccp->numresolutions);
+    cio_write(poc->compno1, (j2k_img->numcomps <= 256 ? 1 : 2));	/* CEpoc_i */
+    poc->compno1 = int_min(poc->compno1, j2k_img->numcomps);
+    cio_write(poc->prg, 1);	/* Ppoc_i */
+  }
+}
+
+void j2k_read_poc()
+{
+  int len, numpchgs, i, old_poc;
+  j2k_tcp_t *tcp;
+  j2k_tccp_t *tccp;
+
+  tcp =
+    j2k_state ==
+    J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
+
+  old_poc = tcp->POC ? tcp->numpocs + 1 : 0;
+  tcp->POC = 1;
+  tccp = &tcp->tccps[0];
+  len = cio_read(2);		/* Lpoc */
+  numpchgs = (len - 2) / (5 + 2 * (j2k_img->numcomps <= 256 ? 1 : 2));
+
+  for (i = old_poc; i < numpchgs + old_poc; i++) {
+    j2k_poc_t *poc;
+    poc = &tcp->pocs[i];
+    poc->resno0 = cio_read(1);	/* RSpoc_i */
+    poc->compno0 = cio_read(j2k_img->numcomps <= 256 ? 1 : 2);	/* CSpoc_i */
+    poc->layno1 = int_min(cio_read(2), tcp->numlayers);	/* LYEpoc_i */
+    poc->resno1 = int_min(cio_read(1), tccp->numresolutions);	/* REpoc_i */
+    poc->compno1 = int_min(cio_read(j2k_img->numcomps <= 256 ? 1 : 2), j2k_img->numcomps);	/* CEpoc_i */
+    poc->prg = cio_read(1);	/* Ppoc_i */
+  }
+
+  tcp->numpocs = numpchgs + old_poc - 1;
+}
+
+void j2k_read_crg()
+{
+  int len, i, Xcrg_i, Ycrg_i;
+
+  len = cio_read(2);		/* Lcrg */
+  for (i = 0; i < j2k_img->numcomps; i++) {
+    Xcrg_i = cio_read(2);	/* Xcrg_i */
+    Ycrg_i = cio_read(2);	/* Ycrg_i */
+  }
+}
+
+void j2k_read_tlm()
+{
+  int len, Ztlm, Stlm, ST, SP, tile_tlm, i;
+  long int Ttlm_i, Ptlm_i;
+
+  len = cio_read(2);		/* Ltlm */
+  Ztlm = cio_read(1);		/* Ztlm */
+  Stlm = cio_read(1);		/* Stlm */
+  ST = ((Stlm >> 4) & 0x01) + ((Stlm >> 4) & 0x02);
+  SP = (Stlm >> 6) & 0x01;
+  tile_tlm = (len - 4) / ((SP + 1) * 2 + ST);
+  for (i = 0; i < tile_tlm; i++) {
+    Ttlm_i = cio_read(ST);	/* Ttlm_i */
+    Ptlm_i = cio_read(SP ? 4 : 2);	/* Ptlm_i */
+  }
+}
+
+void j2k_read_plm()
+{
+  int len, i, Zplm, Nplm, add, packet_len = 0;
+
+  len = cio_read(2);		/* Lplm */
+  Zplm = cio_read(1);		/* Zplm */
+  len -= 3;
+  while (len > 0) {
+    Nplm = cio_read(4);		/* Nplm */
+    len -= 4;
+    for (i = Nplm; i > 0; i--) {
+      add = cio_read(1);
+      len--;
+      packet_len = (packet_len << 7) + add;	/* Iplm_ij */
+      if ((add & 0x80) == 0) {
+	/* New packet */
+	packet_len = 0;
+      }
+      if (len <= 0)
+	break;
+    }
+  }
+}
+
+void j2k_read_plt()
+{
+  int len, i, Zplt, packet_len = 0, add;
+
+  len = cio_read(2);		/* Lplt */
+  Zplt = cio_read(1);		/* Zplt */
+  for (i = len - 3; i > 0; i--) {
+    add = cio_read(1);
+    packet_len = (packet_len << 7) + add;	/* Iplt_i */
+    if ((add & 0x80) == 0) {
+      /* New packet */
+      packet_len = 0;
+    }
+  }
+}
+
+void j2k_read_ppm()
+{
+  int len, Z_ppm, i, j;
+  int N_ppm;
+
+  len = cio_read(2);
+  j2k_cp->ppm = 1;
+
+  Z_ppm = cio_read(1);		/* Z_ppm */
+  len -= 3;
+  while (len > 0) {
+    if (j2k_cp->ppm_previous == 0) {
+      N_ppm = cio_read(4);	/* N_ppm */
+      len -= 4;
+    } else {
+      N_ppm = j2k_cp->ppm_previous;
+    }
+
+    j = j2k_cp->ppm_store;
+    if (Z_ppm == 0) {		/* First PPM marker */
+      j2k_cp->ppm_data =
+	(unsigned char *) calloc(N_ppm, sizeof(unsigned char));
+
+      j2k_cp->ppm_len = N_ppm;	      //Add antonin : ppmbug1
+
+    } else {			/* NON-first PPM marker */
+      j2k_cp->ppm_data =
+	(unsigned char *) realloc(j2k_cp->ppm_data,
+				  (N_ppm +
+				   j2k_cp->ppm_store) *
+				  sizeof(unsigned char));
+
+      j2k_cp->ppm_len = N_ppm + j2k_cp->ppm_store; //Add antonin : ppmbug1
+
+    }
+
+    for (i = N_ppm; i > 0; i--) {	/* Read packet header */
+      j2k_cp->ppm_data[j] = cio_read(1);
+      j++;
+      len--;
+      if (len == 0)
+	break;			/* Case of non-finished packet header in present marker but finished in next one */
+    }
+
+    j2k_cp->ppm_previous = i - 1;
+    j2k_cp->ppm_store = j;
+  }
+}
+
+void j2k_read_ppt()
+{
+  int len, Z_ppt, i, j = 0;
+  j2k_tcp_t *tcp;
+
+  len = cio_read(2);
+  Z_ppt = cio_read(1);
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+  tcp->ppt = 1;
+  if (Z_ppt == 0) {		/* First PPT marker */
+    tcp->ppt_data =
+      (unsigned char *) calloc(len - 3, sizeof(unsigned char));
+    tcp->ppt_store = 0;
+
+    tcp->ppt_len = len-3;	//Add antonin : ppmbug1
+  } else {			/* NON-first PPT marker */
+    tcp->ppt_data =
+      (unsigned char *) realloc(tcp->ppt_data,
+				(len - 3 +
+				 tcp->ppt_store) * sizeof(unsigned char));
+
+    tcp->ppt_len=len - 3 + tcp->ppt_store;    //Add antonin : ppmbug1
+
+  }
+
+  j = tcp->ppt_store;
+  for (i = len - 3; i > 0; i--) {
+    tcp->ppt_data[j] = cio_read(1);
+    j++;
+  }
+  tcp->ppt_store = j;
+}
+
+void j2k_write_sot()
+{
+  int lenp, len;
+
+  j2k_sot_start = cio_tell();
+  cio_write(J2K_MS_SOT, 2);	/* SOT */
+  lenp = cio_tell();
+  cio_skip(2);			/* Lsot (further) */
+  cio_write(j2k_curtileno, 2);	/* Isot */
+  cio_skip(4);			/* Psot (further in j2k_write_sod) */
+  cio_write(0, 1);		/* TPsot */
+  cio_write(1, 1);		/* TNsot */
+  len = cio_tell() - lenp;
+  cio_seek(lenp);
+  cio_write(len, 2);		/* Lsot */
+  cio_seek(lenp + len);
+}
+
+void j2k_read_sot()
+{
+  int len, tileno, totlen, partno, numparts, i;
+  j2k_tcp_t *tcp;
+  j2k_tccp_t *tmp;
+  char status = 0;
+
+  len = cio_read(2);
+  tileno = cio_read(2);
+
+  if (j2k_cp->tileno_size == 0) {
+    j2k_cp->tileno[j2k_cp->tileno_size] = tileno;
+    j2k_cp->tileno_size++;
+  } else {
+    i = 0;
+    while (i < j2k_cp->tileno_size && status == 0) {
+      status = j2k_cp->tileno[i] == tileno ? 1 : 0;
+      i++;
+    }
+    if (status == 0) {
+      j2k_cp->tileno[j2k_cp->tileno_size] = tileno;
+      j2k_cp->tileno_size++;
+    }
+  }
+
+  totlen = cio_read(4);
+  if (!totlen)
+    totlen = cio_numbytesleft() + 8;
+
+  partno = cio_read(1);
+  numparts = cio_read(1);
+
+  j2k_curtileno = tileno;
+  j2k_eot = cio_getbp() - 12 + totlen;
+  j2k_state = J2K_STATE_TPH;
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+
+  if (tcp->first == 1) {
+    tmp = tcp->tccps;
+    *tcp = j2k_default_tcp;
+
+    /* Initialization PPT */
+    tcp->ppt = 0;
+    tcp->ppt_data = NULL;
+
+    tcp->tccps = tmp;
+    for (i = 0; i < j2k_img->numcomps; i++) {
+      tcp->tccps[i] = j2k_default_tcp.tccps[i];
+    }
+    j2k_cp->tcps[j2k_curtileno].first = 0;
+  }
+}
+
+void j2k_write_sod()
+{
+  int l, layno;
+  int totlen;
+  j2k_tcp_t *tcp;
+  static int j2k_sod_start;
+
+  cio_write(J2K_MS_SOD, 2);
+  if (j2k_curtileno == 0) {
+    j2k_sod_start = cio_tell() + pos_correction;
+  }
+
+  /* INDEX >> */
+  if (info_IM.index_on) {
+    info_IM.tile[j2k_curtileno].end_header =
+      cio_tell() + pos_correction - 1;
+    info_IM.tile[j2k_curtileno].packet =
+      (info_packet *) calloc(info_IM.Comp * info_IM.Layer *
+			     (info_IM.Decomposition + 1) * 100,
+			     sizeof(info_packet));
+  }
+  /* << INDEX */
+
+  tcp = &j2k_cp->tcps[j2k_curtileno];
+  for (layno = 0; layno < tcp->numlayers; layno++) {
+    tcp->rates[layno] -= (j2k_sod_start / (j2k_cp->th * j2k_cp->tw));
+  }
+
+  info_IM.num = 0;
+  if (j2k_cp->image_type)
+    l = tcd_encode_tile_pxm(j2k_curtileno, cio_getbp(),
+			    cio_numbytesleft() - 2, &info_IM);
+  else
+    l = tcd_encode_tile_pgx(j2k_curtileno, cio_getbp(),
+			    cio_numbytesleft() - 2, &info_IM);
+
+  /* Writing Psot in SOT marker */
+  totlen = cio_tell() + l - j2k_sot_start;
+  cio_seek(j2k_sot_start + 6);
+  cio_write(totlen, 4);
+  cio_seek(j2k_sot_start + totlen);
+}
+
+void j2k_read_sod()
+{
+  int len, truncate = 0, i;
+  unsigned char *data;
+
+  len = int_min(j2k_eot - cio_getbp(), cio_numbytesleft() + 1);
+  if (len == cio_numbytesleft() + 1)
+    truncate = 1;		/* Case of a truncate codestream */
+
+  data =
+    (unsigned char *) malloc((j2k_tile_len[j2k_curtileno] + len) *
+			     sizeof(unsigned char));
+  for (i = 0; i < j2k_tile_len[j2k_curtileno]; i++)
+    data[i] = j2k_tile_data[j2k_curtileno][i];
+  for (i = 0; i < len; i++)
+    data[i + j2k_tile_len[j2k_curtileno]] = cio_read(1);
+
+  j2k_tile_len[j2k_curtileno] += len;
+  free(j2k_tile_data[j2k_curtileno]);
+  j2k_tile_data[j2k_curtileno] = data;
+  data = NULL;
+
+  if (!truncate)
+    j2k_state = J2K_STATE_TPHSOT;
+  else
+    j2k_state = J2K_STATE_NEOC;	/* RAJOUTE !! */
+}
+
+void j2k_write_rgn(int compno, int tileno)
+{
+  j2k_tcp_t *tcp = &j2k_cp->tcps[tileno];
+
+  cio_write(J2K_MS_RGN, 2);	/* RGN  */
+  cio_write(j2k_img->numcomps <= 256 ? 5 : 6, 2);	/* Lrgn */
+  cio_write(compno, j2k_img->numcomps <= 256 ? 1 : 2);	/* Crgn */
+  cio_write(0, 1);		/* Srgn */
+  cio_write(tcp->tccps[compno].roishift, 1);	/* SPrgn */
+}
+
+void j2k_read_rgn()
+{
+  int len, compno, roisty;
+  j2k_tcp_t *tcp;
+
+  tcp =
+    j2k_state ==
+    J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
+  len = cio_read(2);		/* Lrgn */
+  compno = cio_read(j2k_img->numcomps <= 256 ? 1 : 2);	/* Crgn */
+  roisty = cio_read(1);		/* Srgn */
+  tcp->tccps[compno].roishift = cio_read(1);	/* SPrgn */
+}
+
+void j2k_write_eoc()
+{
+  /* fprintf(stderr, "%.8x: EOC\n", cio_tell() + pos_correction); */
+  cio_write(J2K_MS_EOC, 2);
+}
+
+void j2k_read_eoc()
+{
+  int i, tileno;
+
+  tcd_init(j2k_img, j2k_cp);
+
+  for (i = 0; i < j2k_cp->tileno_size; i++) {
+    tileno = j2k_cp->tileno[i];
+    tcd_decode_tile(j2k_tile_data[tileno], j2k_tile_len[tileno], tileno);
+    free(j2k_tile_data[tileno]);
+  }
+
+  j2k_state = J2K_STATE_MT;
+  longjmp(j2k_error, 1);
+}
+
+void j2k_read_unk()
+{
+  //c'era un bug in questa funzione occorre leggere la lunghezza di un marker sconosciuto per proseguire il parsing 
+	int len;
+	fprintf(stderr, "warning: unknown marker at %x\n",cio_tell());
+	len = cio_read(2);
+	cio_skip(len - 2);			//bug corretto
+}
+
+void j2k_read_jpwlms()
+{
+  //c'era un bug in questa funzione occorre leggere la lunghezza di un marker sconosciuto per proseguire il parsing 
+	int len;
+	len = cio_read(2);
+	cio_skip(len - 2);			//bug corretto
+}
+
+LIBJ2K_API int j2k_encode(j2k_image_t * img, j2k_cp_t * cp, char *output,
+			  int len, char *index)
+{
+  int tileno, compno, layno, resno, precno, pack_nb;
+  char *dest = NULL;
+  FILE *INDEX = NULL;
+  FILE *f = NULL;
+
+  if (setjmp(j2k_error)) {
+    return 0;
+  }
+
+  if (cp->intermed_file==1) {
+    f = fopen(output, "wb");
+    if (!f) {
+      fprintf(stderr, "failed to open %s for writing\n", output);
+      return 1;
+    }
+    dest = (char *) malloc(len);
+    cio_init(dest, len);
+  }
+
+  j2k_img = img;
+  j2k_cp = cp;
+  /* j2k_dump_cp(j2k_img, j2k_cp); */
+
+  /* INDEX >> */
+  info_IM.index_on = j2k_cp->index_on;
+  if (info_IM.index_on) {
+    info_IM.tile =
+      (info_tile *) malloc(j2k_cp->tw * j2k_cp->th * sizeof(info_tile));
+    info_IM.Im_w = j2k_img->x1 - j2k_img->x0;
+    info_IM.Im_h = j2k_img->y1 - j2k_img->y0;
+    info_IM.Prog = (&j2k_cp->tcps[0])->prg;
+    info_IM.tw=j2k_cp->tw; 
+    info_IM.th=j2k_cp->th; 
+    info_IM.Tile_x = j2k_cp->tdx;	/* new version parser */
+    info_IM.Tile_y = j2k_cp->tdy;	/* new version parser */
+    info_IM.Comp = j2k_img->numcomps;
+    info_IM.Layer = (&j2k_cp->tcps[0])->numlayers;
+    info_IM.Decomposition = (&j2k_cp->tcps[0])->tccps->numresolutions - 1;
+    info_IM.D_max = 0;		/* ADD Marcela */
+  }
+  /* << INDEX */
+
+  j2k_write_soc();
+  j2k_write_siz();
+  j2k_write_cod();
+  j2k_write_qcd();
+  for (compno = 0; compno < j2k_img->numcomps; compno++) {
+    j2k_tcp_t *tcp = &j2k_cp->tcps[0];
+    if (tcp->tccps[compno].roishift)
+      j2k_write_rgn(compno, 0);
+  }
+  if (j2k_cp->comment != NULL)
+    j2k_write_com();
+
+  if (cp->intermed_file==1) {
+    /* Writing the main header */
+    pos_correction = cio_tell();
+    fwrite(dest, 1, cio_tell(), f);
+  }
+
+  /* INDEX >> */
+  if (info_IM.index_on) {
+    info_IM.Main_head_end = cio_tell() - 1;
+  }
+  /* << INDEX */
+
+
+  for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
+    fprintf(stderr, "Tile number %d / %d ", tileno + 1,
+      cp->tw * cp->th);
+    
+    if (cp->intermed_file==1) {
+      /* new dest for each tile  */
+      free(dest);
+      dest = (char *) malloc(len);
+      cio_init(dest, len);
+    }
+    j2k_curtileno = tileno;
+    /* initialisation before tile encoding  */
+
+    if (tileno == 0) {
+      tcd_malloc_encode(j2k_img, j2k_cp, j2k_curtileno);
+    } else {
+      tcd_init_encode(j2k_img, j2k_cp, j2k_curtileno);
+    }
+
+    /* INDEX >> */
+    if (info_IM.index_on) {
+      info_IM.tile[j2k_curtileno].num_tile = j2k_curtileno;
+      info_IM.tile[j2k_curtileno].start_pos = cio_tell() + pos_correction;
+    }
+    /* << INDEX */
+    j2k_write_sot();
+
+    for (compno = 1; compno < img->numcomps; compno++) {
+      j2k_write_coc(compno);
+      j2k_write_qcc(compno);
+    }
+
+    if (cp->tcps[tileno].numpocs)
+      j2k_write_poc();
+    j2k_write_sod();
+
+    /* INDEX >> */
+    if (info_IM.index_on) {
+      info_IM.tile[j2k_curtileno].end_pos =
+	cio_tell() + pos_correction - 1;
+    }
+    /* << INDEX */
+
+    /*
+       if (tile->PPT)  BAD PPT !!!
+       {
+       FILE *PPT_file;
+
+       int i;
+       PPT_file=fopen("PPT","rb");
+       fprintf(stderr,"%c%c%c%c",255,97,tile->len_ppt/256,tile->len_ppt%256);
+       for (i=0;i<tile->len_ppt;i++)
+       {
+       unsigned char elmt;
+       fread(&elmt, 1, 1, PPT_file);
+       fwrite(&elmt,1,1,f);
+       }
+       fclose(PPT_file);
+       unlink("PPT");
+       }
+     */
+    if (cp->intermed_file==1) {
+      fwrite(dest, 1, cio_tell(), f);
+      pos_correction = cio_tell() + pos_correction;
+    }
+  }
+
+  if (cp->intermed_file==1) {
+    free(dest);
+    dest = (char *) malloc(len);
+    cio_init(dest, len);
+  }
+
+  j2k_write_eoc();
+
+  if (cp->intermed_file==1) {
+    fwrite(dest, 1, 2, f);
+    free(dest);
+    /* closing file *.j2k */
+    fclose(f);
+  }
+
+  /* Creation of the index file     */
+
+  if (info_IM.index_on) {
+
+    double DistoTotal = 0;
+
+    info_IM.codestream_size = cio_tell() + pos_correction;	/* Correction 14/4/03 suite rmq de Patrick */
+
+    INDEX = fopen(index, "w");
+
+
+
+    if (!INDEX) {
+
+      fprintf(stderr, "failed to open %s for writing\n", index);
+
+      return 1;
+
+    }
+
+
+
+    fprintf(INDEX, "%d %d\n", info_IM.Im_w, info_IM.Im_h);
+
+    fprintf(INDEX, "%d\n", info_IM.Prog);
+
+    fprintf(INDEX, "%d %d\n", info_IM.Tile_x, info_IM.Tile_y);
+
+    fprintf(INDEX, "%d %d\n", info_IM.tw, info_IM.th);
+
+    fprintf(INDEX, "%d\n", info_IM.Comp);
+
+    fprintf(INDEX, "%d\n", info_IM.Layer);
+
+    fprintf(INDEX, "%d\n", info_IM.Decomposition);
+
+    for (resno=info_IM.Decomposition;resno>=0;resno--) {
+
+      fprintf(INDEX, "[%d,%d] ", (1<<info_IM.tile[0].pdx[resno]), (1<<info_IM.tile[0].pdx[resno])); //based on tile 0
+
+    }
+
+    fprintf(INDEX,"\n");
+
+    fprintf(INDEX, "%d\n", info_IM.Main_head_end);
+
+    fprintf(INDEX, "%d\n", info_IM.codestream_size);
+
+    for (tileno = 0; tileno < info_IM.tw * info_IM.th; tileno++) {
+
+      fprintf(INDEX, "%4d %9d %9d %9d %9e %9d %9e\n",
+
+	      info_IM.tile[tileno].num_tile,
+
+	      info_IM.tile[tileno].start_pos,
+
+	      info_IM.tile[tileno].end_header,
+
+	      info_IM.tile[tileno].end_pos, info_IM.tile[tileno].distotile, info_IM.tile[tileno].nbpix,
+
+	      info_IM.tile[tileno].distotile / info_IM.tile[tileno].nbpix);
+
+    }
+
+    for (tileno = 0; tileno < info_IM.tw * info_IM.th; tileno++) {
+
+      int start_pos, end_pos;
+
+      double disto = 0;
+
+      pack_nb = 0;
+
+      /* fprintf(INDEX,
+
+	      "pkno tileno layerno resno compno precno start_pos   end_pos       deltaSE        \n");*/
+
+      if (info_IM.Prog == 0) {	/* LRCP */
+
+	for (layno = 0; layno < info_IM.Layer; layno++) {
+
+	  for (resno = 0; resno < info_IM.Decomposition + 1; resno++) {
+
+	    for (compno = 0; compno < info_IM.Comp; compno++) {
+
+	      for (precno = 0;
+
+		   precno <
+
+		   info_IM.tile[tileno].pw[resno] * info_IM.tile[tileno].ph[resno];
+
+		   precno++) {
+
+		start_pos = info_IM.tile[tileno].packet[pack_nb].start_pos;
+
+		end_pos = info_IM.tile[tileno].packet[pack_nb].end_pos;
+
+		disto = info_IM.tile[tileno].packet[pack_nb].disto;
+
+		fprintf(INDEX, "%4d %6d %7d %5d %6d %6d %9d %9d %8e\n",
+
+			pack_nb, tileno, layno, resno, compno, precno,
+
+			start_pos, end_pos, disto);
+
+		DistoTotal += disto;
+
+		pack_nb++;
+
+	      }
+
+	    }
+
+	  }
+
+	}
+
+      } else if (info_IM.Prog == 1) {	/* RLCP */
+
+	for (resno = 0; resno < info_IM.Decomposition + 1; resno++) {
+
+	  for (layno = 0; layno < info_IM.Layer; layno++) {
+
+	    for (compno = 0; compno < info_IM.Comp; compno++) {
+
+	      for (precno = 0; precno < info_IM.tile[tileno].pw[resno] * info_IM.tile[tileno].ph[resno]; precno++) {
+
+		start_pos = info_IM.tile[tileno].packet[pack_nb].start_pos;
+
+		end_pos = info_IM.tile[tileno].packet[pack_nb].end_pos;
+
+		disto = info_IM.tile[tileno].packet[pack_nb].disto;
+
+		fprintf(INDEX, "%4d %6d %7d %5d %6d %6d %9d %9d %8e\n",
+
+			pack_nb, tileno, layno, resno, compno, precno,
+
+			start_pos, end_pos, disto);
+
+		DistoTotal += disto;
+
+		pack_nb++;
+
+	      }
+
+	    }
+
+	  }
+
+	}
+
+      } else if (info_IM.Prog == 2) {	/* RPCL */
+
+	for (resno = 0; resno < info_IM.Decomposition + 1; resno++) {
+
+	  for (precno = 0; precno < info_IM.tile[tileno].pw[resno] * info_IM.tile[tileno].ph[resno]; precno++) {
+
+	    for (compno = 0; compno < info_IM.Comp; compno++) {
+
+	      for (layno = 0; layno < info_IM.Layer; layno++) {
+
+		start_pos = info_IM.tile[tileno].packet[pack_nb].start_pos;
+
+		end_pos = info_IM.tile[tileno].packet[pack_nb].end_pos;
+
+		disto = info_IM.tile[tileno].packet[pack_nb].disto;
+
+		fprintf(INDEX, "%4d %6d %7d %5d %6d %6d %9d %9d %8e\n",
+
+			pack_nb, tileno, layno, resno, compno, precno,
+
+			start_pos, end_pos, disto);
+
+		DistoTotal += disto;
+
+		pack_nb++;
+
+	      }
+
+	    }
+
+	  }
+
+	}
+
+      } else if (info_IM.Prog == 3) {	/* PCRL */
+
+	for (precno = 0; precno < info_IM.tile[tileno].pw[resno] * info_IM.tile[tileno].ph[resno]; precno++) {
+
+	  for (compno = 0; compno < info_IM.Comp; compno++) {
+
+	    for (resno = 0; resno < info_IM.Decomposition + 1; resno++) {
+
+	      for (layno = 0; layno < info_IM.Layer; layno++) {
+
+		start_pos = info_IM.tile[tileno].packet[pack_nb].start_pos;
+
+		end_pos = info_IM.tile[tileno].packet[pack_nb].end_pos;
+
+		disto = info_IM.tile[tileno].packet[pack_nb].disto;
+
+		fprintf(INDEX, "%4d %6d %7d %5d %6d %6d %9d %9d %8e\n",
+
+			pack_nb, tileno, layno, resno, compno, precno,
+
+			start_pos, end_pos, disto);
+
+		DistoTotal += disto;
+
+		pack_nb++;
+
+	      }
+
+	    }
+
+	  }
+
+	}
+
+      } else {			/* CPRL */
+
+
+
+	for (compno = 0; compno < info_IM.Comp; compno++) {
+
+	  for (precno = 0; precno < info_IM.tile[tileno].pw[resno] * info_IM.tile[tileno].ph[resno]; precno++) {
+
+	    for (resno = 0; resno < info_IM.Decomposition + 1; resno++) {
+
+	      for (layno = 0; layno < info_IM.Layer; layno++) {
+
+		start_pos = info_IM.tile[tileno].packet[pack_nb].start_pos;
+
+		end_pos = info_IM.tile[tileno].packet[pack_nb].end_pos;
+
+		disto = info_IM.tile[tileno].packet[pack_nb].disto;
+
+		fprintf(INDEX, "%4d %6d %7d %5d %6d %6d %9d %9d %8e\n",
+
+			pack_nb, tileno, layno, resno, compno, precno,
+
+			start_pos, end_pos, disto);
+
+		DistoTotal += disto;
+
+		pack_nb++;
+
+	      }
+
+	    }
+
+	  }
+
+	}
+
+      }
+
+    }
+
+    fprintf(INDEX, "SE max : %8e\n", info_IM.D_max);
+
+    fprintf(INDEX, "SE total : %.8e\n", DistoTotal);
+
+    fclose(INDEX);
+
+  }
+
+  j2k_clean();
+
+  return cio_tell();
+}
+
+typedef struct {
+  int id;
+  int states;
+  void (*handler) ();
+} j2k_dec_mstabent_t;
+
+j2k_dec_mstabent_t j2k_dec_mstab[] = {
+  {J2K_MS_SOC, J2K_STATE_MHSOC, j2k_read_soc},
+  {J2K_MS_SOT, J2K_STATE_MH | J2K_STATE_TPHSOT, j2k_read_sot},
+  {J2K_MS_SOD, J2K_STATE_TPH, j2k_read_sod},
+  {J2K_MS_EOC, J2K_STATE_TPHSOT, j2k_read_eoc},
+  {J2K_MS_SIZ, J2K_STATE_MHSIZ, j2k_read_siz},
+  {J2K_MS_COD, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_cod},
+  {J2K_MS_COC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_coc},
+  {J2K_MS_RGN, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_rgn},
+  {J2K_MS_QCD, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_qcd},
+  {J2K_MS_QCC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_qcc},
+  {J2K_MS_POC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_poc},
+  {J2K_MS_TLM, J2K_STATE_MH, j2k_read_tlm},
+  {J2K_MS_PLM, J2K_STATE_MH, j2k_read_plm},
+  {J2K_MS_PLT, J2K_STATE_TPH, j2k_read_plt},
+  {J2K_MS_PPM, J2K_STATE_MH, j2k_read_ppm},
+  {J2K_MS_PPT, J2K_STATE_TPH, j2k_read_ppt},
+  {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},
+  {JPWL_MS_EPB, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_jpwlms},
+  {JPWL_MS_EPC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_jpwlms},	
+  {JPWL_MS_ESD, J2K_STATE_MH | J2K_STATE_TPH, jpwl_read_esd},
+  {JPWL_MS_RED, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_jpwlms},
+  {0, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_unk}
+};
+
+j2k_dec_mstabent_t *j2k_dec_mstab_lookup(int id)
+{
+  j2k_dec_mstabent_t *e;
+  for (e = j2k_dec_mstab; e->id != 0; e++) {
+    if (e->id == id) {
+      break;
+    }
+  }
+  return e;
+}
+
+
+
+LIBJ2K_API int j2k_decode(unsigned char *src, int len, j2k_image_t * img,
+
+			  j2k_cp_t * cp)
+{
+	/*int temp;
+	printf("ciao!\n");
+	temp=setjmp(j2k_error);
+	 printf("temp: %d\n");*/
+  if (setjmp(j2k_error)) {
+	  //printf("temp: %d\n");
+	  if (j2k_state != J2K_STATE_MT) {
+      fprintf(stderr, "WARNING: incomplete bitstream\n");
+      return 0;
+    }
+    j2k_clean();
+    return cio_numbytes();    /* Correct way of ending j2k_decode */
+	//return len;
+  }
+
+  j2k_img = img;
+
+  j2k_cp = cp;
+
+  j2k_state = J2K_STATE_MHSOC;
+  cio_init(src, len);
+
+  
+  if (decode_JPWL(src, len)) {
+	  //printf("ciao!\n");
+	  _exit(-4);
+	  //printf("Sto leggendo una codestream JPWL!\n");
+	  nesd = 0;
+	  len = cslen;
+	  cio_init(src, len);
+  }
+  //printf("ciao!\n");
+  //else
+	  //printf("Non ho riconosciuto la codestream come JPWL!\n");
+  //cio_init(src, cslen);
+  // len = cslen;
+  exit(-5);
+  //cio_init(src, len);
+  for (;;) {
+    j2k_dec_mstabent_t *e;
+    int id = cio_read(2);
+    if (id >> 8 != 0xff) {
+      fprintf(stderr, "%.8x: expected a marker instead of %x\n",
+	      cio_tell() - 2, id);
+      return 0;
+    }
+    e = j2k_dec_mstab_lookup(id);
+    if (!(j2k_state & e->states)) {
+      fprintf(stderr, "%.8x: unexpected marker %x\n", cio_tell() - 2, id);
+      return 0;
+    }
+    if (e->handler) {
+      (*e->handler) ();
+    }
+    if (j2k_state == J2K_STATE_NEOC)
+      break;			/* RAJOUTE */
+  }
+  if (j2k_state == J2K_STATE_NEOC)
+    j2k_read_eoc();		/* RAJOUTE */
+
+  return 0;
+}
+void jpwl_read_esd()
+{
+	unsigned int lesd, cesd2;
+	unsigned char pesd, cesd1;
+	unsigned long esdata;
+	int mode, metric, i, temp;
+	short int ex, mant;
+	double value;
+	FILE *f;
+
+	
+	//printf("Eccomi qua!\n");
+	lesd = cio_read(2);
+	//esd[nesd].lesd = cio_read(2);
+	if (j2k_img->numcomps<257)
+	{
+		cesd1 = cio_read(1);
+		//esd[nesd].cesd1 = cio_read(1);
+		esdata = lesd - 4;
+	}
+	else
+	{
+		cesd2 = cio_read(2);
+		//esd[nesd].cesd2 = cio_read(2);
+		esdata = lesd - 5;
+	}
+	pesd = cio_read(1);
+	//esd[nesd].pesd = cio_read(1);
+	if (j2k_state==J2K_STATE_MH)
+	{
+			if ((f = fopen("esdmap","w"))==NULL)
+		{
+			printf("Unable to open file 'esdmap'!\n");
+			exit(1);
+		}
+	}
+	else
+	if ((f = fopen("esdmap","a"))==NULL)
+	{
+		printf("Unable to open file 'esdmap'!\n");
+		exit(1);
+	}
+	mode = (pesd >> 6)&3; // 0:pm, 1:brm, 2:prm, 3:reserved !
+	//mode = (esd[nesd].pesd >> 6)&3; // 0:pm, 1:brm, 2:prm, 3:reserved !
+	metric = (pesd >> 3) & 7; // 0:relative, 1:MSE, 2:MSEred, 3:PSNR, 4:PSNRincr, 5:MAXERR !
+	//metric = (esd[nesd].pesd >> 3) & 7; // 0:relative, 1:MSE, 2:MSEred, 3:PSNR, 4:PSNRincr, 5:MAXERR !
+	fprintf(f,"Lesd: %x\n",lesd);
+	if (j2k_img->numcomps<257)
+		fprintf(f,"Cesd: %x\n",cesd1);
+	else
+		fprintf(f,"Cesd: %x\n",cesd2);
+	fprintf(f,"Pesd: %x\n",pesd);
+	//fprintf(f,"mode: %d\n",mode);
+	//fprintf(f,"metric: %d\n",metric);
+	switch (mode) {
+		case 0: // packed mode
+			//esdata = esdata/(1+((pesd>>2)&1)+2*(1+((pesd>>1)&1))); //numero livelli
+			esdata = esdata/(1+((pesd>>2)&1)); //numero livelli
+			fprintf(f,"LEVEL\tVALUE\n");
+			for (i=0; i<esdata; i++)
+			{
+				//fprintf(f,"%d\t%d\n",i,cio_read(1+((pesd>>2)&1)));
+				temp = cio_read(1+((pesd>>2)&1));
+				if (metric!=0)
+				{
+					mant = temp & 0x07FF;
+					ex = (temp >> 11) & 0x001F;
+					value = pow(2,ex-15)*(1+(mant/pow(2,11)));
+					//fprintf(f,"%x\n",cio_read(1));
+					fprintf(f,"%d\t%f\n",i,value);
+				}
+				else	// relative error sensivity
+				{
+					fprintf(f,"%d\t%d\n",i,temp);
+				}
+			}
+			break;
+		case 1: // byte range mode
+			esdata = esdata/(1+((pesd>>2)&1)+2*(2+2*((pesd>>1)&1))); //numero livelli
+			fprintf(f,"LEVEL\tSTART BYTE\tEND BYTE\tVALUE\n");
+			for (i=0; i<esdata; i++)
+			{
+				//fprintf(f,"%d\t%d\n",i,cio_read(1+((pesd>>2)&1)));
+				//fprintf(f,"%d\t%d\t%d\t",i,cio_read(2+2*((pesd>>1)&1)),
+				//	cio_read(2+2*((pesd>>1)&1)));
+				fprintf(f,"%d\t%d\t\t",i,cio_read(2+2*((pesd>>1)&1)));
+				fprintf(f,"%d\t\t",cio_read(2+2*((pesd>>1)&1)));
+				temp = cio_read(1+((pesd>>2)&1));
+				if (metric!=0)
+				{
+					mant = temp & 0x07FF;
+					ex = (temp >> 11) & 0x001F;
+					value = pow(2,ex-15)*(1+(mant/pow(2,11)));
+					//fprintf(f,"%x\n",cio_read(1));
+					fprintf(f,"%f\n",value);
+				}
+				else	// relative error sensivity
+				{
+					fprintf(f,"%d\n",temp);
+				}
+			}
+			break;
+		case 2: // packed range mode
+			esdata = esdata/(1+((pesd>>2)&1)+2*(2+2*((pesd>>1)&1))); //numero livelli
+			//printf("esdata: %d\n",esdata);
+			fprintf(f,"LEVEL\tSTART PCK\tEND PCK\t\tVALUE\n");
+			for (i=0; i<esdata; i++)
+			{
+				//fprintf(f,"%d\t%d\t\t%d\t\t",i,cio_read(2+2*((pesd>>1)&1)),
+				//	cio_read(2+2*((pesd>>1)&1))); // Cosi' non funziona bene...mah!!!
+				fprintf(f,"%d\t%d\t\t",i,cio_read(2+2*((pesd>>1)&1)));
+				fprintf(f,"%d\t\t",cio_read(2+2*((pesd>>1)&1)));
+				temp = cio_read(1+((pesd>>2)&1));
+				if (metric!=0)
+				{
+					mant = temp & 0x07FF;
+					ex = (temp >> 11) & 0x001F;
+					value = pow(2,ex-15)*(1+(mant/pow(2,11)));
+					//fprintf(f,"%x\n",cio_read(1));
+					fprintf(f,"%f\n",value);
+				}
+				else	// relative error sensivity
+				{
+					fprintf(f,"%d\n",temp);
+				}
+			}
+			//printf("Ho letto ESD!\n");
+			break;
+	}
+	
+	//esdata = esdata/2;
+	//for (i=0; i<esdata; i++)
+	//{
+	//	temp = cio_read(2);
+	//	mant = temp & 0x07FF;
+	//	ex = (temp >> 11) & 0x001F;
+	//	value = pow(2,ex-15)*(1+(mant/pow(2,11)));
+		//fprintf(f,"%x\n",cio_read(1));
+	//	fprintf(f,"level %d: %f\n",i,value);
+	//}
+	//printf("temp: %x\n",cio_read(2));
+	//cio_skip(-2);
+	fclose(f);
+}
+/*
+ * Read a JPT-stream and decode file
+ *
+ */
+int j2k_decode_jpt_stream(unsigned char *src, int len, j2k_image_t * img,
+			  j2k_cp_t * cp)
+{
+  jpt_msg_header_struct_t header;
+  int position;
+
+  if (setjmp(j2k_error)) {
+    if (j2k_state != J2K_STATE_MT) {
+      fprintf(stderr, "WARNING: incomplete bitstream\n");
+      return 0;
+    }
+    return cio_numbytes();
+  }
+
+  j2k_img = img;
+
+  j2k_cp = cp;
+
+  j2k_state = J2K_STATE_MHSOC;
+  cio_init(src, len);
+
+  /* Initialize the header */
+  jpt_init_Msg_Header(&header);
+  /* Read the first header of the message */
+  jpt_read_Msg_Header(&header);
+
+  position = cio_tell();
+  if (header.Class_Id != 6) {	/* 6 : Main header data-bin message */
+    fprintf(stderr,
+	    "[JPT-stream] : Expecting Main header first [class_Id %d] !\n",
+	    header.Class_Id);
+    return 0;
+  }
+
+  for (;;) {
+    j2k_dec_mstabent_t *e;
+    int id;
+
+    if (!cio_numbytesleft()) {
+      j2k_read_eoc();
+      return 0;
+    }
+    /* data-bin read -> need to read a new header */
+    if ((unsigned int) (cio_tell() - position) == header.Msg_length) {
+      jpt_read_Msg_Header(&header);
+      position = cio_tell();
+      if (header.Class_Id != 4) {	/* 4 : Tile data-bin message */
+	fprintf(stderr, "[JPT-stream] : Expecting Tile info !\n");
+	return 0;
+      }
+    }
+
+    id = cio_read(2);
+    if (id >> 8 != 0xff) {
+      fprintf(stderr, "%.8x: expected a marker instead of %x\n",
+	      cio_tell() - 2, id);
+      return 0;
+    }
+    e = j2k_dec_mstab_lookup(id);
+    if (!(j2k_state & e->states)) {
+      fprintf(stderr, "%.8x: unexpected marker %x\n", cio_tell() - 2, id);
+      return 0;
+    }
+    if (e->handler) {
+      (*e->handler) ();
+    }
+    if (j2k_state == J2K_STATE_NEOC)
+      break;			/* RAJOUTE */
+  }
+  if (j2k_state == J2K_STATE_NEOC)
+    j2k_read_eoc();		/* RAJOUTE */
+
+  return 0;
+}
+
+#ifdef WIN32
+#include <windows.h>
+
+BOOL APIENTRY DllMain(HANDLE hModule, DWORD ul_reason_for_call,
+		      LPVOID lpReserved)
+{
+  switch (ul_reason_for_call) {
+  case DLL_PROCESS_ATTACH:
+  case DLL_THREAD_ATTACH:
+  case DLL_THREAD_DETACH:
+  case DLL_PROCESS_DETACH:
+    break;
+  }
+  return TRUE;
+}
+#endif