1228 lines
35 KiB
C
1228 lines
35 KiB
C
/*
|
|
* Copyright (c) 2001-2002, David Janssens
|
|
* Copyright (c) 2003-2004, Yannick Verschueren
|
|
* Copyright (c) 2003-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 <unistd.h>
|
|
|
|
#include "j2k.h"
|
|
#include "cio.h"
|
|
#include "tcd.h"
|
|
#include "int.h"
|
|
#include "jpip.h"
|
|
#include "jp2.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 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 START_NB 10
|
|
#define INCREMENT 5
|
|
|
|
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 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_t img;
|
|
|
|
|
|
void j2k_clean() {
|
|
int tileno = 0;
|
|
int compno=0, resno=0, precno=0;
|
|
|
|
tcd_free(j2k_img, j2k_cp);
|
|
for (tileno = 0; tileno < j2k_cp->tw * j2k_cp->th; tileno++) {
|
|
info_tile_t *tile_Idx = &img.tile[tileno];
|
|
|
|
for (compno = 0; compno < img.Comp; compno++)
|
|
{
|
|
info_compo_t *compo_Idx = &tile_Idx->compo[compno];
|
|
for(resno = 0; resno < img.Decomposition + 1; resno++)
|
|
{
|
|
info_reso_t *reso_Idx = &compo_Idx->reso[resno];
|
|
for (precno = 0; precno < img.tile[tileno].pw * img.tile[tileno].ph; precno++)
|
|
{
|
|
info_prec_t *prec_Idx = &reso_Idx->prec[precno];
|
|
free(prec_Idx->layer);
|
|
}
|
|
free(reso_Idx->prec);
|
|
}
|
|
free(compo_Idx->reso);
|
|
}
|
|
free(tile_Idx->compo);
|
|
free(tile_Idx->marker);
|
|
free(tile_Idx->tile_parts);
|
|
free(tile_Idx->marker_mul.COC);
|
|
free(tile_Idx->marker_mul.RGN);
|
|
free(tile_Idx->marker_mul.QCC);
|
|
free(tile_Idx->marker_mul.PLT);
|
|
free(tile_Idx->marker_mul.PPT);
|
|
free(tile_Idx->marker_mul.COM);
|
|
}
|
|
free(img.tile);
|
|
free(img.marker);
|
|
free(img.marker_mul.COC);
|
|
free(img.marker_mul.RGN);
|
|
free(img.marker_mul.QCC);
|
|
free(img.marker_mul.PLM);
|
|
free(img.marker_mul.PPM);
|
|
free(img.marker_mul.COM);
|
|
}
|
|
|
|
|
|
|
|
void j2k_read_soc() {
|
|
j2k_state=J2K_STATE_MHSIZ;
|
|
}
|
|
|
|
|
|
|
|
void j2k_read_siz() {
|
|
int len, i;
|
|
info_tile_t *tile;
|
|
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
img.marker[img.num_marker].type = J2K_MS_SIZ;
|
|
img.marker[img.num_marker].start_pos = cio_tell()-2;
|
|
img.marker[img.num_marker].len = len;
|
|
img.num_marker++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
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); /* YTOsiz */
|
|
|
|
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);
|
|
j2k_img->comps[i].prec = (tmp & 0x7f) + 1;
|
|
j2k_img->comps[i].sgnd = tmp >> 7;
|
|
j2k_img->comps[i].dx = cio_read(1);
|
|
j2k_img->comps[i].dy = cio_read(1);
|
|
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].data = (int*)malloc(sizeof(int) * w * h);
|
|
}
|
|
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));
|
|
|
|
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*)malloc(j2k_img->numcomps * sizeof(j2k_tccp_t));
|
|
for (i = 0; i < j2k_cp->tw * j2k_cp->th; i++) {
|
|
j2k_cp->tcps[i].tccps = (j2k_tccp_t*)malloc(j2k_img->numcomps * sizeof(j2k_tccp_t));
|
|
}
|
|
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;
|
|
|
|
/* <INDEX> */
|
|
img.Im_w = j2k_img->x1 - j2k_img->x0;
|
|
img.Im_h = j2k_img->y1 - j2k_img->y0;
|
|
img.Tile_x = j2k_cp->tdx;
|
|
img.Tile_y = j2k_cp->tdy;
|
|
img.Comp = j2k_img->numcomps;
|
|
img.tw = j2k_cp->tw;
|
|
img.th = j2k_cp->th;
|
|
img.tile = (info_tile_t*)malloc(img.tw * img.th * sizeof(info_tile_t));
|
|
|
|
for (i = 0; i < img.tw * img.th; i++)
|
|
{
|
|
tile = &img.tile[i];
|
|
tile->marker = (info_marker_t*)malloc(32 * sizeof(info_marker_t));
|
|
tile->num_marker = 0;
|
|
tile->marker_mul.num_COC = 0;
|
|
tile->marker_mul.CzCOC = START_NB;
|
|
tile->marker_mul.num_RGN = 0;
|
|
tile->marker_mul.CzRGN = START_NB;
|
|
tile->marker_mul.num_QCC = 0;
|
|
tile->marker_mul.CzQCC = START_NB;
|
|
tile->marker_mul.num_PLT = 0;
|
|
tile->marker_mul.CzPLT = START_NB;
|
|
tile->marker_mul.num_PPT = 0;
|
|
tile->marker_mul.CzPPT = START_NB;
|
|
tile->marker_mul.num_COM = 0;
|
|
tile->marker_mul.CzCOM = START_NB;
|
|
}
|
|
/* </INDEX> */
|
|
|
|
|
|
}
|
|
|
|
void j2k_read_com() {
|
|
int len;
|
|
info_tile_t *tile;
|
|
info_marker_t *tmp;
|
|
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (j2k_state == J2K_STATE_MH)
|
|
{
|
|
if (!img.marker_mul.num_COM)
|
|
img.marker_mul.COM = (info_marker_t*)malloc(img.marker_mul.CzCOM * sizeof(info_marker_t));
|
|
if (img.marker_mul.num_COM >= img.marker_mul.CzCOM)
|
|
{
|
|
tmp = (info_marker_t*)malloc(2 * img.marker_mul.CzCOM * sizeof(info_marker_t));
|
|
memcpy(tmp,img.marker_mul.COM,img.marker_mul.CzCOM);
|
|
img.marker_mul.CzCOM *= 2;
|
|
free(img.marker_mul.COM);
|
|
img.marker_mul.COM = tmp;
|
|
}
|
|
|
|
img.marker_mul.COM[img.marker_mul.num_COM].type = J2K_MS_COM;
|
|
img.marker_mul.COM[img.marker_mul.num_COM].start_pos = cio_tell()-2;
|
|
img.marker_mul.COM[img.marker_mul.num_COM].len = len;
|
|
img.marker_mul.num_COM++;
|
|
} else
|
|
{
|
|
tile = &img.tile[j2k_curtileno];
|
|
if (!tile->marker_mul.num_COM)
|
|
tile->marker_mul.COM = (info_marker_t*)calloc(START_NB, sizeof(info_marker_t));
|
|
if (tile->marker_mul.num_COM >= tile->marker_mul.CzCOM)
|
|
{
|
|
tmp = (info_marker_t*)malloc(2 * tile->marker_mul.CzCOM * sizeof(info_marker_t));
|
|
memcpy(tmp,tile->marker_mul.COM,tile->marker_mul.CzCOM);
|
|
tile->marker_mul.CzCOM *= 2;
|
|
free(tile->marker_mul.COM);
|
|
tile->marker_mul.COM = tmp;
|
|
}
|
|
tile->marker_mul.COM[tile->marker_mul.num_COM].type = J2K_MS_COM;
|
|
tile->marker_mul.COM[tile->marker_mul.num_COM].start_pos = cio_tell()-2;
|
|
tile->marker_mul.COM[tile->marker_mul.num_COM].len = len;
|
|
tile->marker_mul.num_COM++;
|
|
}
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
cio_skip(len - 2);
|
|
}
|
|
|
|
|
|
|
|
|
|
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;
|
|
|
|
img.Decomposition = tccp->numresolutions - 1; /* <INDEX> */
|
|
|
|
tccp->cblkw = cio_read(1) + 2;
|
|
tccp->cblkh = cio_read(1) + 2;
|
|
tccp->cblksty = cio_read(1);
|
|
tccp->qmfbid = cio_read(1);
|
|
if (tccp->csty&J2K_CP_CSTY_PRT) {
|
|
for (i = 0; i < tccp->numresolutions; i++) {
|
|
int tmp = cio_read(1);
|
|
tccp->prcw[i] = tmp&0xf;
|
|
tccp->prch[i] = tmp>>4;
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_cod() {
|
|
int len, i, pos;
|
|
j2k_tcp_t *tcp;
|
|
info_tile_t *tile;
|
|
|
|
tcp = j2k_state == J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (j2k_state == J2K_STATE_MH)
|
|
{
|
|
img.marker[img.num_marker].type = J2K_MS_COD;
|
|
img.marker[img.num_marker].start_pos = cio_tell()-2;
|
|
img.marker[img.num_marker].len = len;
|
|
img.num_marker++;
|
|
}
|
|
else
|
|
{
|
|
tile = &img.tile[j2k_curtileno];
|
|
tile->marker[tile->num_marker].type = J2K_MS_COD;
|
|
tile->marker[tile->num_marker].start_pos = cio_tell()-2;
|
|
tile->marker[tile->num_marker].len = len;
|
|
tile->num_marker++;
|
|
}
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
tcp->csty = cio_read(1);
|
|
tcp->prg = cio_read(1);
|
|
tcp->numlayers = cio_read(2);
|
|
tcp->mct = cio_read(1);
|
|
|
|
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);
|
|
}
|
|
|
|
/* <INDEX> */
|
|
img.Prog = tcp->prg;
|
|
img.Layer = tcp->numlayers;
|
|
/* </INDEX> */
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_coc() {
|
|
int len, compno;
|
|
j2k_tcp_t *tcp;
|
|
info_tile_t *tile;
|
|
info_marker_t *tmp;
|
|
|
|
tcp = j2k_state == J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (j2k_state == J2K_STATE_MH)
|
|
{
|
|
if (!img.marker_mul.num_COC)
|
|
img.marker_mul.COC = (info_marker_t*)malloc(img.marker_mul.CzCOC * sizeof(info_marker_t));
|
|
if (img.marker_mul.num_COC >= img.marker_mul.CzCOC)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + img.marker_mul.CzCOC) * sizeof(info_marker_t));
|
|
memcpy(tmp,img.marker_mul.COC,img.marker_mul.CzCOC);
|
|
img.marker_mul.CzCOC += INCREMENT;
|
|
free(img.marker_mul.COC);
|
|
img.marker_mul.COC = tmp;
|
|
}
|
|
img.marker_mul.COC[img.marker_mul.num_COC].type = J2K_MS_COC;
|
|
img.marker_mul.COC[img.marker_mul.num_COC].start_pos = cio_tell()-2;
|
|
img.marker_mul.COC[img.marker_mul.num_COC].len = len;
|
|
img.marker_mul.num_COC++;
|
|
} else
|
|
{
|
|
tile = &img.tile[j2k_curtileno];
|
|
if (!tile->marker_mul.num_COC)
|
|
tile->marker_mul.COC = (info_marker_t*)malloc(tile->marker_mul.CzCOC * sizeof(info_marker_t));
|
|
if (tile->marker_mul.num_COC >= tile->marker_mul.CzCOC)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + tile->marker_mul.CzCOC) * sizeof(info_marker_t));
|
|
memcpy(tmp,tile->marker_mul.COC,tile->marker_mul.CzCOC);
|
|
tile->marker_mul.CzCOC += INCREMENT;
|
|
free(tile->marker_mul.COC);
|
|
tile->marker_mul.COC = tmp;
|
|
}
|
|
tile->marker_mul.COC[tile->marker_mul.num_COC].type = J2K_MS_COC;
|
|
tile->marker_mul.COC[tile->marker_mul.num_COC].start_pos = cio_tell() - 2;
|
|
tile->marker_mul.COC[tile->marker_mul.num_COC].len = len;
|
|
tile->marker_mul.num_COC++;
|
|
}
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
compno =cio_read(j2k_img->numcomps <= 256 ? 1 : 2);
|
|
|
|
tcp->tccps[compno].csty = cio_read(1);
|
|
j2k_read_cox(compno);
|
|
}
|
|
|
|
|
|
|
|
|
|
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);
|
|
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) { /* WHY STEPSIZES WHEN NOQNT ? */
|
|
expn = cio_read(1) >> 3;
|
|
mant = 0;
|
|
} else {
|
|
tmp = cio_read(2);
|
|
expn = tmp >> 11;
|
|
mant = tmp & 0x7ff;
|
|
}
|
|
tccp->stepsizes[bandno].expn = expn;
|
|
tccp->stepsizes[bandno].mant = mant;
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_qcd() {
|
|
int len, i, pos;
|
|
info_tile_t *tile;
|
|
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (j2k_state == J2K_STATE_MH)
|
|
{
|
|
img.marker[img.num_marker].type = J2K_MS_QCD;
|
|
img.marker[img.num_marker].start_pos = cio_tell()-2;
|
|
img.marker[img.num_marker].len = len;
|
|
img.num_marker++;
|
|
} else
|
|
{
|
|
tile = &img.tile[j2k_curtileno];
|
|
tile->marker[tile->num_marker].type = J2K_MS_QCD;
|
|
tile->marker[tile->num_marker].start_pos = cio_tell()-2;
|
|
tile->marker[tile->num_marker].len = len;
|
|
tile->num_marker++;
|
|
}
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
|
|
pos=cio_tell();
|
|
for (i = 0; i < j2k_img->numcomps; i++) {
|
|
cio_seek(pos);
|
|
j2k_read_qcx(i, len - 2);
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_qcc() {
|
|
int len, compno;
|
|
info_tile_t *tile;
|
|
info_marker_t *tmp;
|
|
|
|
len = cio_read(2);
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (j2k_state == J2K_STATE_MH)
|
|
{
|
|
if (!img.marker_mul.num_QCC)
|
|
img.marker_mul.QCC = (info_marker_t*)malloc(img.marker_mul.CzQCC * sizeof(info_marker_t));
|
|
if (img.marker_mul.num_QCC >= img.marker_mul.CzQCC)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + img.marker_mul.CzQCC) * sizeof(info_marker_t));
|
|
memcpy(tmp,img.marker_mul.QCC,img.marker_mul.CzQCC);
|
|
img.marker_mul.CzQCC += INCREMENT;
|
|
free(img.marker_mul.QCC);
|
|
img.marker_mul.QCC = tmp;
|
|
}
|
|
img.marker_mul.QCC[img.marker_mul.num_QCC].type = J2K_MS_QCC;
|
|
img.marker_mul.QCC[img.marker_mul.num_QCC].start_pos = cio_tell() - 2;
|
|
img.marker_mul.QCC[img.marker_mul.num_QCC].len = len;
|
|
img.marker_mul.num_QCC++;
|
|
} else
|
|
{
|
|
tile = &img.tile[j2k_curtileno];
|
|
if (!tile->marker_mul.num_QCC)
|
|
tile->marker_mul.QCC = (info_marker_t*)malloc(tile->marker_mul.CzQCC * sizeof(info_marker_t));
|
|
if (tile->marker_mul.num_QCC >= tile->marker_mul.CzQCC)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + tile->marker_mul.CzQCC) * sizeof(info_marker_t));
|
|
memcpy(tmp,tile->marker_mul.QCC,tile->marker_mul.CzQCC);
|
|
tile->marker_mul.CzQCC += INCREMENT;
|
|
free(tile->marker_mul.QCC);
|
|
tile->marker_mul.QCC = tmp;
|
|
}
|
|
tile->marker_mul.QCC[tile->marker_mul.num_QCC].type = J2K_MS_QCC;
|
|
tile->marker_mul.QCC[tile->marker_mul.num_QCC].start_pos = cio_tell()-2;
|
|
tile->marker_mul.QCC[tile->marker_mul.num_QCC].len = len;
|
|
tile->marker_mul.num_QCC++;
|
|
}
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
compno = cio_read(j2k_img->numcomps <= 256 ? 1 : 2);
|
|
j2k_read_qcx(compno, len - 2 - (j2k_img->numcomps <= 256 ? 1 : 2));
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_poc() {
|
|
int len, numpchgs, i, old_poc;
|
|
j2k_tcp_t *tcp;
|
|
j2k_tccp_t *tccp;
|
|
info_tile_t *tile;
|
|
|
|
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);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (j2k_state == J2K_STATE_MH)
|
|
{
|
|
img.marker[img.num_marker].type = J2K_MS_POC;
|
|
img.marker[img.num_marker].start_pos = cio_tell()-2;
|
|
img.marker[img.num_marker].len = len;
|
|
img.num_marker++;
|
|
} else
|
|
{
|
|
tile = &img.tile[j2k_curtileno];
|
|
tile->marker[tile->num_marker].type = J2K_MS_POC;
|
|
tile->marker[tile->num_marker].start_pos = cio_tell()-2;
|
|
tile->marker[tile->num_marker].len = len;
|
|
tile->num_marker++;
|
|
}
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
numpchgs = (len - 2) / (5 + 2 * (j2k_img->numcomps <= 256 ? 1 : 2));
|
|
for (i = 0; i < numpchgs; i++) {
|
|
j2k_poc_t *poc;
|
|
poc = &tcp->pocs[i];
|
|
poc->resno0 = cio_read(1);
|
|
poc->compno0 = cio_read(j2k_img->numcomps <= 256 ? 1 : 2);
|
|
poc->layno1 = int_min(cio_read(2), tcp->numlayers);
|
|
poc->resno1 = int_min(cio_read(1), tccp->numresolutions);
|
|
poc->compno1 = int_min(cio_read(j2k_img->numcomps <= 256 ? 1 : 2), j2k_img->numcomps);
|
|
poc->prg = cio_read(1);
|
|
}
|
|
|
|
tcp->numpocs = numpchgs + old_poc - 1;
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_crg() {
|
|
int len, i, Xcrg_i, Ycrg_i;
|
|
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
img.marker[img.num_marker].type = J2K_MS_CRG;
|
|
img.marker[img.num_marker].start_pos = cio_tell()-2;
|
|
img.marker[img.num_marker].len = len;
|
|
img.num_marker++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
for (i = 0; i < j2k_img->numcomps; i++)
|
|
{
|
|
Xcrg_i = cio_read(2);
|
|
Ycrg_i = cio_read(2);
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_tlm() {
|
|
int len, Ztlm, Stlm, ST, SP, tile_tlm, i;
|
|
long int Ttlm_i, Ptlm_i;
|
|
info_marker_t *tmp;
|
|
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (!img.marker_mul.num_TLM)
|
|
img.marker_mul.TLM = (info_marker_t*)malloc(img.marker_mul.CzTLM * sizeof(info_marker_t));
|
|
if (img.marker_mul.num_TLM >= img.marker_mul.CzTLM)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + img.marker_mul.CzTLM) * sizeof(info_marker_t));
|
|
memcpy(tmp,img.marker_mul.TLM,img.marker_mul.CzTLM);
|
|
img.marker_mul.CzTLM += INCREMENT;
|
|
free(img.marker_mul.TLM);
|
|
img.marker_mul.TLM = tmp;
|
|
}
|
|
img.marker_mul.TLM[img.marker_mul.num_TLM].type = J2K_MS_TLM;
|
|
img.marker_mul.TLM[img.marker_mul.num_TLM].start_pos = cio_tell()-2;
|
|
img.marker_mul.TLM[img.marker_mul.num_TLM].len = len;
|
|
img.marker_mul.num_TLM++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
Ztlm = cio_read(1);
|
|
Stlm = cio_read(1);
|
|
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);
|
|
Ptlm_i = cio_read(SP ? 4 : 2);
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_plm() {
|
|
int len, i, Z_plm, N_plm, add, packet_len=0;
|
|
info_marker_t *tmp;
|
|
|
|
len=cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (!img.marker_mul.num_PLM)
|
|
img.marker_mul.PLM = (info_marker_t*)malloc(img.marker_mul.CzPLM * sizeof(info_marker_t));
|
|
if (img.marker_mul.num_PLM >= img.marker_mul.CzPLM)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + img.marker_mul.CzPLM) * sizeof(info_marker_t));
|
|
memcpy(tmp,img.marker_mul.PLM,img.marker_mul.CzPLM);
|
|
img.marker_mul.CzPLM += INCREMENT;
|
|
free(img.marker_mul.PLM);
|
|
img.marker_mul.PLM = tmp;
|
|
}
|
|
img.marker_mul.PLM[img.marker_mul.num_PLM].type = J2K_MS_PLM;
|
|
img.marker_mul.PLM[img.marker_mul.num_PLM].start_pos = cio_tell()-2;
|
|
img.marker_mul.PLM[img.marker_mul.num_PLM].len = len;
|
|
img.marker_mul.num_PLM++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
Z_plm = cio_read(1);
|
|
len -= 3;
|
|
while (len > 0)
|
|
{
|
|
N_plm = cio_read(4);
|
|
len -= 4;
|
|
for (i = N_plm ; i > 0 ; i--)
|
|
{
|
|
add = cio_read(1);
|
|
len--;
|
|
packet_len = (packet_len << 7) + add;
|
|
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;
|
|
info_tile_t *tile;
|
|
info_marker_t *tmp;
|
|
;
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
tile = &img.tile[j2k_curtileno];
|
|
if (!tile->marker_mul.num_PLT)
|
|
tile->marker_mul.PLT = (info_marker_t*)malloc(tile->marker_mul.CzPLT * sizeof(info_marker_t));
|
|
if (tile->marker_mul.num_PLT >= tile->marker_mul.CzPLT)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + tile->marker_mul.CzPLT) * sizeof(info_marker_t));
|
|
memcpy(tmp,tile->marker_mul.PLT,tile->marker_mul.CzPLT);
|
|
tile->marker_mul.CzPLT += INCREMENT;
|
|
free(tile->marker_mul.PLT);
|
|
tile->marker_mul.PLT = tmp;
|
|
}
|
|
|
|
tile->marker_mul.PLT[tile->marker_mul.num_PLT].type = J2K_MS_PLT;
|
|
tile->marker_mul.PLT[tile->marker_mul.num_PLT].start_pos = cio_tell()-2;
|
|
tile->marker_mul.PLT[tile->marker_mul.num_PLT].len = len;
|
|
tile->marker_mul.num_PLT++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
Zplt = cio_read(1);
|
|
for (i = len-3; i > 0; i--)
|
|
{
|
|
add = cio_read(1);
|
|
packet_len = (packet_len << 7) + add;
|
|
if ((add & 0x80) == 0)
|
|
{
|
|
/* New packet */
|
|
packet_len = 0;
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_ppm() {
|
|
int len, Z_ppm, i, j;
|
|
int N_ppm;
|
|
info_marker_t *tmp;
|
|
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
if (!img.marker_mul.num_PPM)
|
|
img.marker_mul.PPM = (info_marker_t*)malloc(img.marker_mul.CzPPM * sizeof(info_marker_t));
|
|
if (img.marker_mul.num_PPM >= img.marker_mul.CzPPM)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + img.marker_mul.CzPPM) * sizeof(info_marker_t));
|
|
memcpy(tmp,img.marker_mul.PPM,img.marker_mul.CzPPM);
|
|
img.marker_mul.CzPPM += INCREMENT;
|
|
free(img.marker_mul.PPM);
|
|
img.marker_mul.PPM = tmp;
|
|
}
|
|
img.marker_mul.PLM[img.marker_mul.num_PPM].type = J2K_MS_PPM;
|
|
img.marker_mul.PLM[img.marker_mul.num_PPM].start_pos = cio_tell()-2;
|
|
img.marker_mul.PLM[img.marker_mul.num_PPM].len = len;
|
|
img.marker_mul.num_PPM++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
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));
|
|
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));
|
|
|
|
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;
|
|
info_tile_t *tile;
|
|
len=cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX] */
|
|
tile = & img.tile[j2k_curtileno];
|
|
tile->marker[tile->num_marker].type = J2K_MS_PPT;
|
|
tile->marker[tile->num_marker].start_pos = cio_tell()-2;
|
|
tile->marker[tile->num_marker].len = len;
|
|
tile->num_marker++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
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;
|
|
}
|
|
else /* NON-first PPT marker */
|
|
tcp->ppt_data = (unsigned char*)realloc(tcp->ppt_data, (len - 3 + tcp->ppt_store) * sizeof(unsigned char));
|
|
|
|
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_read_sot() {
|
|
int len, tileno, totlen, partno, numparts, i;
|
|
j2k_tcp_t *tcp;
|
|
j2k_tccp_t *tmp;
|
|
info_tile_t *tile;
|
|
info_tile_part_t *tilepart_tmp;
|
|
|
|
|
|
//fprintf(stderr,"SOT\n");
|
|
len = cio_read(2);
|
|
tileno = cio_read(2);
|
|
/* <INDEX> [MHIX BOX] */
|
|
tile = & img.tile[tileno];
|
|
tile->marker[tile->num_marker].type = J2K_MS_SOT;
|
|
tile->marker[tile->num_marker].start_pos = cio_tell() - 4;
|
|
tile->marker[tile->num_marker].len = len;
|
|
tile->num_marker++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
totlen = cio_read(4);
|
|
if (!totlen) totlen = cio_numbytesleft() + 8;
|
|
partno = cio_read(1);
|
|
numparts = cio_read(1);
|
|
|
|
/* <INDEX> */
|
|
if (tileno == 0 && partno == 0 )
|
|
img.Main_head_end = cio_tell() - 12; /* Correction End = First byte of first SOT */
|
|
|
|
img.tile[tileno].num_tile = tileno;
|
|
/* </INDEX> */
|
|
|
|
tile->numparts = partno + 1; /* INDEX : Number of tile_parts for the tile */
|
|
img.num_max_tile_parts = int_max(tile->numparts, img.num_max_tile_parts); /* INDEX : Maximum number of tile_part per tile */
|
|
|
|
if (partno == 0)
|
|
{
|
|
tile->tile_parts = (info_tile_part_t*)malloc(START_NB * sizeof(info_tile_part_t*));
|
|
tile->Cztile_parts = START_NB;
|
|
}
|
|
if (partno >= tile->Cztile_parts)
|
|
{
|
|
tilepart_tmp = (info_tile_part_t*)malloc((INCREMENT + tile->Cztile_parts) * sizeof(info_tile_part_t));
|
|
memcpy(tilepart_tmp, tile->tile_parts, tile->Cztile_parts); /*add tilepart_ a tmp*/
|
|
tile->Cztile_parts += INCREMENT;
|
|
free(tile->tile_parts);
|
|
tile->tile_parts = tilepart_tmp;
|
|
}
|
|
|
|
tile->tile_parts[partno].start_pos = cio_tell() - 12; /* INDEX : start_pos of the tile_part */
|
|
tile->tile_parts[partno].length = totlen; /* INDEX : length of the tile_part */
|
|
tile->tile_parts[partno].end_pos = totlen + cio_tell() - 12; /* INDEX : end position of the tile_part */
|
|
//tile->tile_parts[partno].num_reso_AUX = j2k_default_tcp.tccps[0].numresolutions;
|
|
|
|
j2k_curtileno = tileno;
|
|
j2k_eot = cio_getbp() - 12 + totlen;
|
|
j2k_state = J2K_STATE_TPH;
|
|
tcp = &j2k_cp->tcps[j2k_curtileno];
|
|
|
|
if (partno == 0)
|
|
// 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;
|
|
}
|
|
// tile->tile_parts[numparts].num_reso_AUX = tcp->tccps[0].numresolutions; /* INDEX : AUX value for TPIX */
|
|
}
|
|
|
|
|
|
|
|
void j2k_read_rgn() {
|
|
int len, compno, roisty;
|
|
j2k_tcp_t *tcp;
|
|
info_tile_t *tile;
|
|
info_marker_t *tmp;
|
|
// fprintf(stderr,"RGN\n");
|
|
tcp = j2k_state == J2K_STATE_TPH ? &j2k_cp->tcps[j2k_curtileno] : &j2k_default_tcp;
|
|
len = cio_read(2);
|
|
|
|
/* <INDEX> [MHIX BOX]*/
|
|
if (j2k_state == J2K_STATE_MH)
|
|
{
|
|
if (!img.marker_mul.num_RGN)
|
|
img.marker_mul.RGN = (info_marker_t*)malloc(img.marker_mul.CzRGN * sizeof(info_marker_t));
|
|
if (img.marker_mul.num_RGN >= img.marker_mul.CzRGN)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + img.marker_mul.CzRGN) * sizeof(info_marker_t));
|
|
memcpy(tmp,img.marker_mul.RGN, img.marker_mul.CzRGN);
|
|
img.marker_mul.CzRGN += INCREMENT;
|
|
free(img.marker_mul.RGN);
|
|
img.marker_mul.RGN = tmp;
|
|
}
|
|
img.marker_mul.RGN[img.marker_mul.num_RGN].type = J2K_MS_RGN;
|
|
img.marker_mul.RGN[img.marker_mul.num_RGN].start_pos = cio_tell() - 2;
|
|
img.marker_mul.RGN[img.marker_mul.num_RGN].len = len;
|
|
img.marker_mul.num_RGN++;
|
|
} else
|
|
{
|
|
tile = &img.tile[j2k_curtileno];
|
|
if (!tile->marker_mul.num_RGN)
|
|
tile->marker_mul.RGN = (info_marker_t*)malloc(tile->marker_mul.CzRGN * sizeof(info_marker_t));
|
|
if (tile->marker_mul.num_RGN >= tile->marker_mul.CzRGN)
|
|
{
|
|
tmp = (info_marker_t*)malloc((INCREMENT + tile->marker_mul.CzRGN) * sizeof(info_marker_t));
|
|
memcpy(tmp,tile->marker_mul.RGN,tile->marker_mul.CzRGN);
|
|
tile->marker_mul.CzRGN += INCREMENT;
|
|
free(tile->marker_mul.RGN);
|
|
tile->marker_mul.RGN = tmp;
|
|
}
|
|
|
|
tile->marker_mul.RGN[tile->marker_mul.num_RGN].type = J2K_MS_RGN;
|
|
tile->marker_mul.RGN[tile->marker_mul.num_RGN].start_pos = cio_tell() - 2;
|
|
tile->marker_mul.RGN[tile->marker_mul.num_RGN].len = len;
|
|
tile->marker_mul.num_RGN++;
|
|
}
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
compno = cio_read(j2k_img->numcomps <= 256 ? 1 : 2);
|
|
roisty = cio_read(1);
|
|
tcp->tccps[compno].roishift = cio_read(1);
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
void j2k_read_sod() {
|
|
int len, i;
|
|
unsigned char *data;
|
|
info_tile_t *tile;
|
|
info_tile_part_t *tile_part;
|
|
/* <INDEX> [MHIX BOX] */
|
|
tile = &img.tile[j2k_curtileno];
|
|
tile->marker[tile->num_marker].type = J2K_MS_SOD;
|
|
tile->marker[tile->num_marker].start_pos = cio_tell();
|
|
tile->marker[tile->num_marker].len = 0;
|
|
tile->num_marker++;
|
|
/* </INDEX> [MHIX BOX] */
|
|
|
|
tile_part = &tile->tile_parts[tile->numparts - 1]; /* INDEX : Current tilepart of a tile */
|
|
tile_part->length_header = cio_tell() - 1 - tile_part->start_pos; /* INDEX : length of the tile-part header */
|
|
tile_part->end_header = cio_tell() - 1; /* INDEX : end header position of the tile-part header */
|
|
|
|
len = int_min(j2k_eot - cio_getbp(), cio_numbytesleft());
|
|
|
|
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;
|
|
|
|
j2k_state = J2K_STATE_TPHSOT;
|
|
}
|
|
|
|
void j2k_read_eoc() {
|
|
int tileno;
|
|
|
|
tcd_init(j2k_img, j2k_cp, &img);
|
|
for (tileno = 0; tileno<j2k_cp->tw * j2k_cp->th; tileno++) {
|
|
tcd_decode_tile(j2k_tile_data[tileno], j2k_tile_len[tileno], tileno, &img);
|
|
}
|
|
|
|
j2k_state = J2K_STATE_MT;
|
|
longjmp(j2k_error, 1);
|
|
}
|
|
|
|
|
|
|
|
|
|
void j2k_read_unk() {
|
|
fprintf(stderr, "warning: unknown marker\n");
|
|
}
|
|
|
|
|
|
|
|
|
|
int j2k_index_JPIP(char *Idx_file, char *J2K_file, int len, int version){
|
|
FILE *dest;
|
|
unsigned char *index;
|
|
int pos_iptr, end_pos;
|
|
int len_cidx, pos_cidx;
|
|
int len_jp2c, pos_jp2c;
|
|
int len_fidx, pos_fidx;
|
|
|
|
dest=fopen(Idx_file, "wb");
|
|
if (!dest) {
|
|
fprintf(stderr, "Failed to open %s for reading !!\n", Idx_file);
|
|
return 0;
|
|
}
|
|
|
|
/* INDEX MODE JPIP */
|
|
index = (unsigned char*)malloc(len);
|
|
cio_init(index, len);
|
|
jp2_write_jp();
|
|
jp2_write_ftyp();
|
|
|
|
jp2_write_jp2h(j2k_img);
|
|
jp2_write_dbtl(Idx_file);
|
|
|
|
pos_iptr=cio_tell();
|
|
cio_skip(24); /* IPTR further ! */
|
|
|
|
pos_jp2c = cio_tell();
|
|
len_jp2c = jp2_write_jp2c(J2K_file);
|
|
|
|
pos_cidx = cio_tell();
|
|
len_cidx = jpip_write_cidx(pos_jp2c + 8,img, j2k_cp, version); /* Correction len_jp2C --> pos_jp2c + 8 */
|
|
|
|
|
|
pos_fidx = cio_tell();
|
|
len_fidx = jpip_write_fidx(pos_jp2c, len_jp2c, pos_cidx, len_cidx);
|
|
|
|
end_pos = cio_tell();
|
|
|
|
cio_seek(pos_iptr);
|
|
jpip_write_iptr(pos_fidx,len_fidx);
|
|
cio_seek(end_pos);
|
|
|
|
fwrite(index, 1, cio_tell(), dest);
|
|
free(index);
|
|
|
|
fclose(dest);
|
|
return 1;
|
|
}
|
|
|
|
|
|
|
|
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},
|
|
{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;
|
|
}
|
|
|
|
int j2k_decode(unsigned char *src, int len, j2k_image_t **image, j2k_cp_t **cp) {
|
|
if (setjmp(j2k_error)) {
|
|
if (j2k_state != J2K_STATE_MT) {
|
|
fprintf(stderr, "WARNING: incomplete bitstream\n");
|
|
return 0;
|
|
}
|
|
return cio_numbytes();
|
|
}
|
|
j2k_img = (j2k_image_t*)calloc(1, sizeof(j2k_image_t));
|
|
j2k_cp = (j2k_cp_t*)calloc(1, sizeof(j2k_cp_t));
|
|
*image = j2k_img;
|
|
*cp = j2k_cp;
|
|
j2k_state = J2K_STATE_MHSOC;
|
|
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)();
|
|
}
|
|
}
|
|
|
|
}
|
|
|
|
|
|
#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
|
|
|
|
int main(int argc, char **argv)
|
|
{
|
|
FILE *src;
|
|
int totlen;
|
|
unsigned char *j2kfile;
|
|
j2k_image_t *imgg;
|
|
j2k_cp_t *cp;
|
|
int version;
|
|
|
|
if (argc != 4)
|
|
{
|
|
fprintf(stderr,"\nERROR in entry : index_create J2K-file JP2-file version\n\nVersion : 0, 1, 2 or 3\n0 : [faix] 4-byte and no AUX fields\n1 : [faix] 8-byte and no AUX fields\n2 : [faix] 4-byte and AUX fields\n3 : [faix] 8-byte and AUX fields\n");
|
|
return 1;
|
|
}
|
|
|
|
src=fopen(argv[1], "rb");
|
|
if (!src) {
|
|
fprintf(stderr, "Failed to open %s for reading !!\n", argv[1]);
|
|
return 1;
|
|
}
|
|
|
|
/* length of the codestream */
|
|
fseek(src, 0, SEEK_END);
|
|
totlen = ftell(src);
|
|
fseek(src, 0, SEEK_SET);
|
|
|
|
j2kfile = (unsigned char*)malloc(totlen);
|
|
fread(j2kfile, 1, totlen, src);
|
|
fclose(src);
|
|
|
|
img.marker = (info_marker_t*)malloc(32 * sizeof(info_marker_t));
|
|
img.num_marker = 0;
|
|
img.num_max_tile_parts = 0;
|
|
img.marker_mul.num_COC = 0;
|
|
img.marker_mul.CzCOC = START_NB;
|
|
img.marker_mul.num_RGN = 0;
|
|
img.marker_mul.CzRGN = START_NB;
|
|
img.marker_mul.num_QCC = 0;
|
|
img.marker_mul.CzQCC = START_NB;
|
|
img.marker_mul.num_TLM = 0;
|
|
img.marker_mul.CzTLM = START_NB;
|
|
img.marker_mul.num_PLM = 0;
|
|
img.marker_mul.CzPLM = START_NB;
|
|
img.marker_mul.num_PPM = 0;
|
|
img.marker_mul.CzPPM = START_NB;
|
|
img.marker_mul.num_COM = 0;
|
|
img.marker_mul.CzCOM = START_NB;
|
|
|
|
/* decode */
|
|
|
|
if (!j2k_decode(j2kfile, totlen, &imgg, &cp)) {
|
|
fprintf(stderr, "Index_creator: failed to decode image!\n");
|
|
return 1;
|
|
}
|
|
free(j2kfile);
|
|
|
|
// fseek(src, 0, SEEK_SET);
|
|
img.codestream_size = totlen;
|
|
sscanf(argv[3], "%d", &version);
|
|
if (version > 3)
|
|
{
|
|
fprintf(stderr,"Error : value of version unauthorized !! Value accepted : 0, 1, 2 or 3 !!\n");
|
|
return 0;
|
|
}
|
|
|
|
j2k_index_JPIP(argv[2], argv[1], totlen * 2 > 60000 ? totlen * 2 : 60000, version);
|
|
|
|
j2k_clean();
|
|
return 0;
|
|
}
|
|
|