2 * The copyright in this software is being made available under the 2-clauses
3 * BSD License, included below. This software may be subject to other third
4 * party and contributor rights, including patent rights, and no such rights
5 * are granted under this license.
7 * Copyright (c) 2001-2003, David Janssens
8 * Copyright (c) 2002-2003, Yannick Verschueren
9 * Copyright (c) 2003-2005, Francois Devaux and Antonin Descampe
10 * Copyright (c) 2005, Herve Drolon, FreeImage Team
11 * Copyright (c) 2002-2005, Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium
12 * Copyright (c) 2006, M�nica D�ez Garc�a, Image Processing Laboratory, University of Valladolid, Spain
13 * All rights reserved.
15 * Redistribution and use in source and binary forms, with or without
16 * modification, are permitted provided that the following conditions
18 * 1. Redistributions of source code must retain the above copyright
19 * notice, this list of conditions and the following disclaimer.
20 * 2. Redistributions in binary form must reproduce the above copyright
21 * notice, this list of conditions and the following disclaimer in the
22 * documentation and/or other materials provided with the distribution.
24 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
25 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
26 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
27 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
28 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
29 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
30 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
31 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
32 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
33 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
34 * POSSIBILITY OF SUCH DAMAGE.
40 @brief The JPEG-2000 Codestream Reader/Writer (J3D)
42 The functions in J3D.C have for goal to read/write the several parts of the codestream: markers and data.
45 /** @defgroup J3D J3D - JPEG-2000 PART 10 codestream reader/writer */
48 #define J3D_CP_CSTY_PRT 0x01
49 #define J3D_CP_CSTY_SOP 0x02
50 #define J3D_CP_CSTY_EPH 0x04
51 #define J3D_CCP_CSTY_PRT 0x01
53 #define J3D_CCP_CBLKSTY_LAZY 0x01 /* Selective arithmetic coding bypass */
54 #define J3D_CCP_CBLKSTY_RESET 0x02 /* Reset context probabilities on coding pass boundaries */
55 #define J3D_CCP_CBLKSTY_TERMALL 0x04 /* Termination on each coding pass */
56 #define J3D_CCP_CBLKSTY_VSC 0x08 /* Vertically causal context, add also hook for switching off and on 3D context models */
57 #define J3D_CCP_CBLKSTY_PTERM 0x10 /* Predictable termination */
58 #define J3D_CCP_CBLKSTY_SEGSYM 0x20 /* Segmentation symbols are used */
59 #define J3D_CCP_CBLKSTY_3DCTXT 0x40 /* 3D context models (3D-EBCOT) vs 2D context models */
61 #define J3D_CCP_QNTSTY_NOQNT 0 /* Quantization style : no quantization */
62 #define J3D_CCP_QNTSTY_SIQNT 1 /* Quantization style : scalar derived (values signalled only in LLL subband) */
63 #define J3D_CCP_QNTSTY_SEQNT 2 /* Quantization style : scalar expounded (values signalled for each subband) */
65 /* ----------------------------------------------------------------------- */
67 #define J3D_MS_SOC 0xff4f /**< SOC marker value */
68 #define J3D_MS_SOT 0xff90 /**< SOT marker value */
69 #define J3D_MS_SOD 0xff93 /**< SOD marker value */
70 #define J3D_MS_EOC 0xffd9 /**< EOC marker value */
71 #define J3D_MS_CAP 0xff50 /**< CAP marker value */
72 #define J3D_MS_SIZ 0xff51 /**< SIZ marker value */
73 #define J3D_MS_NSI 0xff54 /**< NSI marker value */
74 #define J3D_MS_COD 0xff52 /**< COD marker value */
75 #define J3D_MS_COC 0xff53 /**< COC marker value */
76 #define J3D_MS_RGN 0xff5e /**< RGN marker value */
77 #define J3D_MS_QCD 0xff5c /**< QCD marker value */
78 #define J3D_MS_QCC 0xff5d /**< QCC marker value */
79 #define J3D_MS_POC 0xff5f /**< POC marker value */
80 #define J3D_MS_TLM 0xff55 /**< TLM marker value */
81 #define J3D_MS_PLM 0xff57 /**< PLM marker value */
82 #define J3D_MS_PLT 0xff58 /**< PLT marker value */
83 #define J3D_MS_PPM 0xff60 /**< PPM marker value */
84 #define J3D_MS_PPT 0xff61 /**< PPT marker value */
85 #define J3D_MS_SOP 0xff91 /**< SOP marker value */
86 #define J3D_MS_EPH 0xff92 /**< EPH marker value */
87 #define J3D_MS_CRG 0xff63 /**< CRG marker value */
88 #define J3D_MS_COM 0xff64 /**< COM marker value */
90 #define J3D_MS_DCO 0xff70 /**< DCO marker value */
91 #define J3D_MS_VMS 0xff71 /**< VMS marker value */
92 #define J3D_MS_DFS 0xff72 /**< DFS marker value */
93 #define J3D_MS_ADS 0xff73 /**< ADS marker value */
94 #define J3D_MS_ATK 0xff79 /**< ATK marker value */
95 #define J3D_MS_CBD 0xff78 /**< CBD marker value */
96 #define J3D_MS_MCT 0xff74 /**< MCT marker value */
97 #define J3D_MS_MCC 0xff75 /**< MCC marker value */
98 #define J3D_MS_MCO 0xff77 /**< MCO marker value */
99 #define J3D_MS_NLT 0xff76 /**< NLT marker value */
100 #define J3D_MS_QPD 0xff5a /**< QPD marker value */
101 #define J3D_MS_QPC 0xff5b /**< QPC marker value */
103 /* ----------------------------------------------------------------------- */
104 /* Capability RSIZ parameter, extended */
105 #define J3D_RSIZ_BASIC 0x0000
107 #define J3D_RSIZ_DCO 0x8001 /* Required */
108 #define J3D_RSIZ_VSQNT 0x8002
109 #define J3D_RSIZ_TCQNT 0x8004
110 #define J3D_RSIZ_VMASK 0x8008
111 #define J3D_RSIZ_SSOVL 0x8010
112 #define J3D_RSIZ_ADECS 0x8020
113 #define J3D_RSIZ_ATK 0x8040 /*Required*/
114 #define J3D_RSIZ_SSYMK 0x8080
115 #define J3D_RSIZ_MCT 0x8100 /*Not compatible with DCO*/
116 #define J3D_RSIZ_NLT 0x8200 /*Required*/
117 #define J3D_RSIZ_ASHAP 0x8400
118 #define J3D_RSIZ_PRQNT 0x8800
120 #define J3D_CAP_10 0x00400000
121 /* Arbitrary transformation kernel, 15444-2 */
122 #define J3D_ATK_IRR 0
123 #define J3D_ATK_REV 1
124 #define J3D_ATK_ARB 0
126 #define J3D_ATK_CON 0
127 /* ----------------------------------------------------------------------- */
130 Values that specify the status of the decoding process when decoding the main header.
131 These values may be combined with a | operator.
133 typedef enum J3D_STATUS {
134 /**< a SOC marker is expected */
135 J3D_STATE_MHSOC = 0x0001,
136 /**< a SIZ marker is expected */
137 J3D_STATE_MHSIZ = 0x0002,
138 /**< the decoding process is in the main header */
139 J3D_STATE_MH = 0x0004,
140 /**< the decoding process is in a tile part header and expects a SOT marker */
141 J3D_STATE_TPHSOT = 0x0008,
142 /**< the decoding process is in a tile part header */
143 J3D_STATE_TPH = 0x0010,
144 /**< the EOC marker has just been read */
145 J3D_STATE_MT = 0x0020,
146 /**< the decoding process must not expect a EOC marker because the codestream is truncated */
147 J3D_STATE_NEOC = 0x0040
153 Arbitrary transformation kernel
155 typedef struct opj_atk {
156 /** index of wavelet kernel */
158 /** Numerical type of scaling factor and lifting step parameters */
160 /** Wavelet filter category */
162 /** Wavelet transformation type (REV/IRR) */
164 /** Initial odd/even subsequence */
166 /** Boundary extension method (constant CON / whole-sample symmetric WS) */
168 /** Scaling factor. Only for wt_typ=IRR */
170 /** Number of lifting steps */
172 /** Offset for lifting step s. Only for filt_cat=ARB */
174 /** Base 2 scaling exponent for lifting step s. Only for wt_typ=REV */
176 /** Additive residue for lifting step s. Only for wt_typ=REV */
178 /** Number of lifting coefficients signaled for lifting step s */
180 /** Lifting coefficient k for lifting step s */
181 double Aatk[256][256];
186 Quantization stepsize
188 typedef struct opj_stepsize {
196 Tile-component coding parameters
198 typedef struct opj_tccp {
201 /** number of resolutions of x, y and z-axis */
202 int numresolution[3];
203 /** code-blocks width height & depth*/
205 /** code-block coding style */
207 /** 0: no ATK (only 9-7 or 5-3) 1: ATK defined WT*/
209 /** Arbitrary transformation kernel (15444-2)*/
211 /** DWT identifier for x, y and z-axis (0:WT9-7 1:WT5-3 >1:WT-atk->index) */
213 /** reversible/irreversible wavelet transfomation (0:irrev 1:reversible)*/
215 /** quantisation style */
217 /** stepsizes used for quantization */
218 opj_stepsize_t stepsizes[J3D_MAXBANDS];
219 /** number of guard bits. Table A28 de 15444-1*/
221 /** Region Of Interest shift */
223 /** precinct width height & depth*/
224 int prctsiz[3][J3D_MAXRLVLS];
228 Tile coding parameters : coding/decoding parameters common to all tiles
229 (information like COD, COC in main header)
231 typedef struct opj_tcp {
232 /** 1 : first part-tile of a tile */
236 /** progression order */
238 /** number of layers */
240 /** multi-component transform identifier */
242 /** rates of layers */
244 /** number of progression order changes */
246 /** indicates if a POC marker has been used O:NO, 1:YES */
248 /** progression order changes */
249 opj_poc_t pocs[J3D_MAXRLVLS - 1];
250 /** add fixed_quality */
251 float distoratio[100];
252 /** tile-component coding parameters */
254 /** packet header store there for futur use in t2_decode_packet */
255 unsigned char *ppt_data;
256 /** pointer remaining on the first byte of the first header if ppt is used */
257 unsigned char *ppt_data_first;
258 /** If ppt == 1 --> there was a PPT marker for the present tile */
260 /** used in case of multiple marker PPT (number of info already stored) */
268 typedef struct opj_cp {
269 /** transform format 0: 2DWT, 1: 2D1P, 2: 3DWT, 3: 3RLS */
270 OPJ_TRANSFORM transform_format;
271 /** entropy coding format 0: 2EB, 1: 3EB, 2: 2GR, 3: 3GR, 4: GRI*/
272 OPJ_ENTROPY_CODING encoding_format;
273 /** allocation by rate/distortion */
275 /** allocation by fixed layer */
277 /** add fixed_quality */
279 /** Rsiz: capabilities */
281 /** if != 0, then original dimension divided by 2^(reduce); if == 0 or not used, volume is decoded to the full resolution */
283 /** if != 0, then only the first "layer" layers are decoded; if == 0 or not used, all the quality layers are decoded */
285 /** 0 = no index || 1 = index */
287 /** Big-Endian/Little-endian order */
301 /** comment for coding */
303 /** number of tiles in width, height and depth */
307 /** ID number of the tiles present in the codestream */
309 /** size of the vector tileno */
311 /** tile coding parameters */
316 /** packet header store there for futur use in t2_decode_packet */
317 unsigned char *ppm_data;
318 /** pointer remaining on the first byte of the first header if ppm is used */
319 unsigned char *ppm_data_first;
320 /** if ppm == 1 --> there was a PPM marker for the present tile */
322 /** use in case of multiple marker PPM (number of info already store) */
324 /** use in case of multiple marker PPM (case on non-finished previous info) */
330 Information concerning a packet inside tile
332 typedef struct opj_packet_info {
333 /** start position */
337 /** distorsion introduced */
342 Index structure : information regarding tiles inside volume
344 typedef struct opj_tile_info {
345 /** value of thresh for each layer by tile cfr. Marcela */
347 /** number of tile */
349 /** start position */
351 /** end position of the header */
355 /** precinct number for each resolution level (width, height and depth) */
356 int prctno[3][J3D_MAXRLVLS];
357 /** precinct size (in power of 2), in X for each resolution level */
358 int prctsiz[3][J3D_MAXRLVLS];
359 /** information concerning packets inside tile */
360 opj_packet_info_t *packet;
362 /** add fixed_quality */
364 /** add fixed_quality */
371 typedef struct opj_volume_info {
373 /** transform format 0: 2DWT, 1: 2D1P, 2: 3DWT, 3: 3RLS */
374 OPJ_TRANSFORM transform_format;
375 /** output file format 0: 2EB, 1: 3EB, 2: 2GR, 3: 3GR, 4: GRI*/
376 OPJ_ENTROPY_CODING encoding_format; /** 0 = no index || 1 = index */
378 /** 0 = wt 9-7 || 1 = wt 5-3 || >=2 wt atk defined */
380 /** maximum distortion reduction on the whole volume (add for Marcela) */
384 /** writing the packet in the index with t2_encode_packets */
386 /** volume width, height and depth */
390 /** progression order */
392 /** tile size in x, y and z */
396 /** tile origin in x, y and z */
400 /** number of tiles in X, Y and Z */
404 /** component numbers */
406 /** number of layer */
408 /** number of decomposition in X, Y and Z*/
409 int decomposition[3];
410 /** DC offset (15444-2) */
412 /** main header position */
414 /** codestream's size */
416 /** information regarding tiles inside volume */
417 opj_tile_info_t *tile;
421 JPEG-2000 codestream reader/writer
423 typedef struct opj_j3d {
425 opj_common_ptr cinfo;
426 /** locate in which part of the codestream the decoder is (main header, tile header, end) */
428 /** number of the tile currently concern by coding/decoding */
430 /** locate the position of the end of the tile in the codestream, used to detect a truncated codestream (in j3d_read_sod) */
432 /** locate the start position of the SOT marker of the current coded tile: */
434 /* after encoding the tile, a jump (in j3d_write_sod) is done to the SOT marker to store the value of its length. */
436 /** as the J3D-file is written in several parts during encoding, it enables to make the right correction in position return by cio_tell */
438 /** array used to store the data of each tile */
439 unsigned char **tile_data;
440 /** array used to store the length of each tile */
443 /** decompression only : store decoding parameters common to all tiles */
444 opj_tcp_t *default_tcp;
445 /** pointer to the encoded / decoded volume */
446 opj_volume_t *volume;
447 /** pointer to the coding parameters */
449 /** helper used to write the index file */
450 opj_volume_info_t *volume_info;
451 /** pointer to the byte i/o stream */
455 /** @name Funciones generales */
457 /* ----------------------------------------------------------------------- */
459 Creates a J3D decompression structure
460 @param cinfo Codec context info
461 @return Returns a handle to a J3D decompressor if successful, returns NULL otherwise
463 opj_j3d_t* j3d_create_decompress(opj_common_ptr cinfo);
465 Destroy a J3D decompressor handle
466 @param j3d J3D decompressor handle to destroy
468 void j3d_destroy_decompress(opj_j3d_t *j3d);
470 Setup the decoder decoding parameters using user parameters.
471 Decoding parameters are returned in j3d->cp.
472 @param j3d J3D decompressor handle
473 @param parameters decompression parameters
475 void j3d_setup_decoder(opj_j3d_t *j3d, opj_dparameters_t *parameters);
477 Decode an volume from a JPEG-2000 codestream
478 @param j3d J3D decompressor handle
479 @param cio Input buffer stream
480 @return Returns a decoded volume if successful, returns NULL otherwise
482 opj_volume_t* j3d_decode(opj_j3d_t *j3d, opj_cio_t *cio);
484 Decode an volume form a JPT-stream (JPEG 2000, JPIP)
485 @param j3d J3D decompressor handle
486 @param cio Input buffer stream
487 @return Returns a decoded volume if successful, returns NULL otherwise
489 opj_volume_t* j3d_decode_jpt_stream(opj_j3d_t *j3d, opj_cio_t *cio);
491 Creates a J3D compression structure
492 @param cinfo Codec context info
493 @return Returns a handle to a J3D compressor if successful, returns NULL otherwise
495 opj_j3d_t* j3d_create_compress(opj_common_ptr cinfo);
497 Destroy a J3D compressor handle
498 @param j3d J3D compressor handle to destroy
500 void j3d_destroy_compress(opj_j3d_t *j3d);
502 Setup the encoder parameters using the current volume and using user parameters.
503 Coding parameters are returned in j3d->cp.
504 @param j3d J3D compressor handle
505 @param parameters compression parameters
506 @param volume input filled volume
508 void j3d_setup_encoder(opj_j3d_t *j3d, opj_cparameters_t *parameters,
509 opj_volume_t *volume);
511 Encode an volume into a JPEG-2000 codestream
512 @param j3d J3D compressor handle
513 @param cio Output buffer stream
514 @param volume Volume to encode
515 @param index Name of the index file if required, NULL otherwise
516 @return Returns true if successful, returns false otherwise
518 bool j3d_encode(opj_j3d_t *j3d, opj_cio_t *cio, opj_volume_t *volume,
520 /* ----------------------------------------------------------------------- */