diff options
Diffstat (limited to 'jpwl/decoder/libopenjpeg/j2k.c')
| -rw-r--r-- | jpwl/decoder/libopenjpeg/j2k.c | 1769 |
1 files changed, 0 insertions, 1769 deletions
diff --git a/jpwl/decoder/libopenjpeg/j2k.c b/jpwl/decoder/libopenjpeg/j2k.c deleted file mode 100644 index 7360a108..00000000 --- a/jpwl/decoder/libopenjpeg/j2k.c +++ /dev/null @@ -1,1769 +0,0 @@ - // 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 |