2
* Copyright (c) 2002-2007, Communications and Remote Sensing Laboratory, Universite catholique de Louvain (UCL), Belgium
3
* Copyright (c) 2002-2007, Professor Benoit Macq
4
* Copyright (c) 2001-2003, David Janssens
5
* Copyright (c) 2002-2003, Yannick Verschueren
6
* Copyright (c) 2003-2007, Francois-Olivier Devaux and Antonin Descampe
7
* Copyright (c) 2005, Herve Drolon, FreeImage Team
8
* Copyright (c) 2007, Callum Lerwick <seg@haxxed.com>
11
* Redistribution and use in source and binary forms, with or without
12
* modification, are permitted provided that the following conditions
14
* 1. Redistributions of source code must retain the above copyright
15
* notice, this list of conditions and the following disclaimer.
16
* 2. Redistributions in binary form must reproduce the above copyright
17
* notice, this list of conditions and the following disclaimer in the
18
* documentation and/or other materials provided with the distribution.
20
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
21
* AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
22
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
23
* ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
24
* LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
25
* CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
26
* SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
27
* INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
28
* CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
29
* ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
30
* POSSIBILITY OF SUCH DAMAGE.
33
#include "opj_includes.h"
36
/** @defgroup T1 T1 - Implementation of the tier-1 coding */
39
/** @name Local static functions */
42
static INLINE char t1_getctxno_zc(int f, int orient);
43
static char t1_getctxno_sc(int f);
44
static INLINE int t1_getctxno_mag(int f);
45
static char t1_getspb(int f);
46
static short t1_getnmsedec_sig(int x, int bitpos);
47
static short t1_getnmsedec_ref(int x, int bitpos);
48
static void t1_updateflags(flag_t *flagsp, int s, int stride);
50
Encode significant pass
52
static void t1_enc_sigpass_step(
63
Decode significant pass
65
static void t1_dec_sigpass_step(
74
Encode significant pass
76
static void t1_enc_sigpass(
84
Decode significant pass
86
static void t1_dec_sigpass(
93
Encode refinement pass
95
static void t1_enc_refpass_step(
105
Decode refinement pass
107
static void t1_dec_refpass_step(
116
Encode refinement pass
118
static void t1_enc_refpass(
125
Decode refinement pass
127
static void t1_dec_refpass(
135
static void t1_enc_clnpass_step(
148
static void t1_dec_clnpass_step(
159
static void t1_enc_clnpass(
168
static void t1_dec_clnpass(
173
static double t1_getwmsedec(
185
@param cblk Code-block coding parameters
187
@param compno Component number
191
@param cblksty Code-block style
195
static void t1_encode_cblk(
197
opj_tcd_cblk_enc_t* cblk,
205
opj_tcd_tile_t * tile);
209
@param cblk Code-block coding parameters
211
@param roishift Region of interest shifting value
212
@param cblksty Code-block style
214
static void t1_decode_cblk(
216
opj_tcd_cblk_dec_t* cblk,
225
/* ----------------------------------------------------------------------- */
227
static char t1_getctxno_zc(int f, int orient) {
228
return lut_ctxno_zc[(orient << 8) | (f & T1_SIG_OTH)];
231
static char t1_getctxno_sc(int f) {
232
return lut_ctxno_sc[(f & (T1_SIG_PRIM | T1_SGN)) >> 4];
235
static int t1_getctxno_mag(int f) {
236
int tmp1 = (f & T1_SIG_OTH) ? T1_CTXNO_MAG + 1 : T1_CTXNO_MAG;
237
int tmp2 = (f & T1_REFINE) ? T1_CTXNO_MAG + 2 : tmp1;
241
static char t1_getspb(int f) {
242
return lut_spb[(f & (T1_SIG_PRIM | T1_SGN)) >> 4];
245
static short t1_getnmsedec_sig(int x, int bitpos) {
246
if (bitpos > T1_NMSEDEC_FRACBITS) {
247
return lut_nmsedec_sig[(x >> (bitpos - T1_NMSEDEC_FRACBITS)) & ((1 << T1_NMSEDEC_BITS) - 1)];
250
return lut_nmsedec_sig0[x & ((1 << T1_NMSEDEC_BITS) - 1)];
253
static short t1_getnmsedec_ref(int x, int bitpos) {
254
if (bitpos > T1_NMSEDEC_FRACBITS) {
255
return lut_nmsedec_ref[(x >> (bitpos - T1_NMSEDEC_FRACBITS)) & ((1 << T1_NMSEDEC_BITS) - 1)];
258
return lut_nmsedec_ref0[x & ((1 << T1_NMSEDEC_BITS) - 1)];
261
static void t1_updateflags(flag_t *flagsp, int s, int stride) {
262
flag_t *np = flagsp - stride;
263
flag_t *sp = flagsp + stride;
265
static const flag_t mod[] = {
266
T1_SIG_S, T1_SIG_S|T1_SGN_S,
267
T1_SIG_E, T1_SIG_E|T1_SGN_E,
268
T1_SIG_W, T1_SIG_W|T1_SGN_W,
269
T1_SIG_N, T1_SIG_N|T1_SGN_N
276
flagsp[-1] |= mod[s+2];
278
flagsp[1] |= mod[s+4];
285
static void t1_enc_sigpass_step(
298
opj_mqc_t *mqc = t1->mqc; /* MQC component */
300
flag = vsc ? ((*flagsp) & (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) : (*flagsp);
301
if ((flag & T1_SIG_OTH) && !(flag & (T1_SIG | T1_VISIT))) {
302
v = int_abs(*datap) & one ? 1 : 0;
303
mqc_setcurctx(mqc, t1_getctxno_zc(flag, orient)); /* ESSAI */
304
if (type == T1_TYPE_RAW) { /* BYPASS/LAZY MODE */
305
mqc_bypass_enc(mqc, v);
310
v = *datap < 0 ? 1 : 0;
311
*nmsedec += t1_getnmsedec_sig(int_abs(*datap), bpno + T1_NMSEDEC_FRACBITS);
312
mqc_setcurctx(mqc, t1_getctxno_sc(flag)); /* ESSAI */
313
if (type == T1_TYPE_RAW) { /* BYPASS/LAZY MODE */
314
mqc_bypass_enc(mqc, v);
316
mqc_encode(mqc, v ^ t1_getspb(flag));
318
t1_updateflags(flagsp, v, t1->flags_stride);
324
static void t1_dec_sigpass_step(
335
opj_raw_t *raw = t1->raw; /* RAW component */
336
opj_mqc_t *mqc = t1->mqc; /* MQC component */
338
flag = vsc ? ((*flagsp) & (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) : (*flagsp);
339
if ((flag & T1_SIG_OTH) && !(flag & (T1_SIG | T1_VISIT))) {
340
if (type == T1_TYPE_RAW) {
341
if (raw_decode(raw)) {
342
v = raw_decode(raw); /* ESSAI */
343
*datap = v ? -oneplushalf : oneplushalf;
344
t1_updateflags(flagsp, v, t1->flags_stride);
347
mqc_setcurctx(mqc, t1_getctxno_zc(flag, orient));
348
if (mqc_decode(mqc)) {
349
mqc_setcurctx(mqc, t1_getctxno_sc(flag));
350
v = mqc_decode(mqc) ^ t1_getspb(flag);
351
*datap = v ? -oneplushalf : oneplushalf;
352
t1_updateflags(flagsp, v, t1->flags_stride);
357
} /* VSC and BYPASS by Antonin */
359
static void t1_enc_sigpass(
367
int i, j, k, one, vsc;
369
one = 1 << (bpno + T1_NMSEDEC_FRACBITS);
370
for (k = 0; k < t1->h; k += 4) {
371
for (i = 0; i < t1->w; ++i) {
372
for (j = k; j < k + 4 && j < t1->h; ++j) {
373
vsc = ((cblksty & J2K_CCP_CBLKSTY_VSC) && (j == k + 3 || j == t1->h - 1)) ? 1 : 0;
376
&t1->flags[((j+1) * t1->flags_stride) + i + 1],
377
&t1->data[(j * t1->w) + i],
389
static void t1_dec_sigpass(
396
int i, j, k, one, half, oneplushalf, vsc;
399
oneplushalf = one | half;
400
for (k = 0; k < t1->h; k += 4) {
401
for (i = 0; i < t1->w; ++i) {
402
for (j = k; j < k + 4 && j < t1->h; ++j) {
403
vsc = ((cblksty & J2K_CCP_CBLKSTY_VSC) && (j == k + 3 || j == t1->h - 1)) ? 1 : 0;
406
&t1->flags[((j+1) * t1->flags_stride) + i + 1],
407
&t1->data[(j * t1->w) + i],
415
} /* VSC and BYPASS by Antonin */
417
static void t1_enc_refpass_step(
429
opj_mqc_t *mqc = t1->mqc; /* MQC component */
431
flag = vsc ? ((*flagsp) & (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) : (*flagsp);
432
if ((flag & (T1_SIG | T1_VISIT)) == T1_SIG) {
433
*nmsedec += t1_getnmsedec_ref(int_abs(*datap), bpno + T1_NMSEDEC_FRACBITS);
434
v = int_abs(*datap) & one ? 1 : 0;
435
mqc_setcurctx(mqc, t1_getctxno_mag(flag)); /* ESSAI */
436
if (type == T1_TYPE_RAW) { /* BYPASS/LAZY MODE */
437
mqc_bypass_enc(mqc, v);
441
*flagsp |= T1_REFINE;
445
static void t1_dec_refpass_step(
456
opj_mqc_t *mqc = t1->mqc; /* MQC component */
457
opj_raw_t *raw = t1->raw; /* RAW component */
459
flag = vsc ? ((*flagsp) & (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) : (*flagsp);
460
if ((flag & (T1_SIG | T1_VISIT)) == T1_SIG) {
461
mqc_setcurctx(mqc, t1_getctxno_mag(flag)); /* ESSAI */
462
if (type == T1_TYPE_RAW) {
467
t = v ? poshalf : neghalf;
468
*datap += *datap < 0 ? -t : t;
469
*flagsp |= T1_REFINE;
471
} /* VSC and BYPASS by Antonin */
473
static void t1_enc_refpass(
480
int i, j, k, one, vsc;
482
one = 1 << (bpno + T1_NMSEDEC_FRACBITS);
483
for (k = 0; k < t1->h; k += 4) {
484
for (i = 0; i < t1->w; ++i) {
485
for (j = k; j < k + 4 && j < t1->h; ++j) {
486
vsc = ((cblksty & J2K_CCP_CBLKSTY_VSC) && (j == k + 3 || j == t1->h - 1)) ? 1 : 0;
489
&t1->flags[((j+1) * t1->flags_stride) + i + 1],
490
&t1->data[(j * t1->w) + i],
501
static void t1_dec_refpass(
507
int i, j, k, one, poshalf, neghalf;
511
neghalf = bpno > 0 ? -poshalf : -1;
512
for (k = 0; k < t1->h; k += 4) {
513
for (i = 0; i < t1->w; ++i) {
514
for (j = k; j < k + 4 && j < t1->h; ++j) {
515
vsc = ((cblksty & J2K_CCP_CBLKSTY_VSC) && (j == k + 3 || j == t1->h - 1)) ? 1 : 0;
518
&t1->flags[((j+1) * t1->flags_stride) + i + 1],
519
&t1->data[(j * t1->w) + i],
527
} /* VSC and BYPASS by Antonin */
529
static void t1_enc_clnpass_step(
542
opj_mqc_t *mqc = t1->mqc; /* MQC component */
544
flag = vsc ? ((*flagsp) & (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) : (*flagsp);
548
if (!(*flagsp & (T1_SIG | T1_VISIT))) {
549
mqc_setcurctx(mqc, t1_getctxno_zc(flag, orient));
550
v = int_abs(*datap) & one ? 1 : 0;
554
*nmsedec += t1_getnmsedec_sig(int_abs(*datap), bpno + T1_NMSEDEC_FRACBITS);
555
mqc_setcurctx(mqc, t1_getctxno_sc(flag));
556
v = *datap < 0 ? 1 : 0;
557
mqc_encode(mqc, v ^ t1_getspb(flag));
558
t1_updateflags(flagsp, v, t1->flags_stride);
561
*flagsp &= ~T1_VISIT;
564
static void t1_dec_clnpass_step(
575
opj_mqc_t *mqc = t1->mqc; /* MQC component */
577
flag = vsc ? ((*flagsp) & (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) : (*flagsp);
581
if (!(flag & (T1_SIG | T1_VISIT))) {
582
mqc_setcurctx(mqc, t1_getctxno_zc(flag, orient));
583
if (mqc_decode(mqc)) {
585
mqc_setcurctx(mqc, t1_getctxno_sc(flag));
586
v = mqc_decode(mqc) ^ t1_getspb(flag);
587
*datap = v ? -oneplushalf : oneplushalf;
588
t1_updateflags(flagsp, v, t1->flags_stride);
591
*flagsp &= ~T1_VISIT;
592
} /* VSC and BYPASS by Antonin */
594
static void t1_enc_clnpass(
601
int i, j, k, one, agg, runlen, vsc;
603
opj_mqc_t *mqc = t1->mqc; /* MQC component */
606
one = 1 << (bpno + T1_NMSEDEC_FRACBITS);
607
for (k = 0; k < t1->h; k += 4) {
608
for (i = 0; i < t1->w; ++i) {
610
if (cblksty & J2K_CCP_CBLKSTY_VSC) {
611
agg = !(MACRO_t1_flags(1 + k,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
612
|| MACRO_t1_flags(1 + k + 1,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
613
|| MACRO_t1_flags(1 + k + 2,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
614
|| (MACRO_t1_flags(1 + k + 3,1 + i)
615
& (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) & (T1_SIG | T1_VISIT | T1_SIG_OTH));
617
agg = !(MACRO_t1_flags(1 + k,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
618
|| MACRO_t1_flags(1 + k + 1,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
619
|| MACRO_t1_flags(1 + k + 2,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
620
|| MACRO_t1_flags(1 + k + 3,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH));
626
for (runlen = 0; runlen < 4; ++runlen) {
627
if (int_abs(t1->data[((k + runlen)*t1->w) + i]) & one)
630
mqc_setcurctx(mqc, T1_CTXNO_AGG);
631
mqc_encode(mqc, runlen != 4);
635
mqc_setcurctx(mqc, T1_CTXNO_UNI);
636
mqc_encode(mqc, runlen >> 1);
637
mqc_encode(mqc, runlen & 1);
641
for (j = k + runlen; j < k + 4 && j < t1->h; ++j) {
642
vsc = ((cblksty & J2K_CCP_CBLKSTY_VSC) && (j == k + 3 || j == t1->h - 1)) ? 1 : 0;
645
&t1->flags[((j+1) * t1->flags_stride) + i + 1],
646
&t1->data[(j * t1->w) + i],
651
agg && (j == k + runlen),
658
static void t1_dec_clnpass(
664
int i, j, k, one, half, oneplushalf, agg, runlen, vsc;
665
int segsym = cblksty & J2K_CCP_CBLKSTY_SEGSYM;
667
opj_mqc_t *mqc = t1->mqc; /* MQC component */
671
oneplushalf = one | half;
672
for (k = 0; k < t1->h; k += 4) {
673
for (i = 0; i < t1->w; ++i) {
675
if (cblksty & J2K_CCP_CBLKSTY_VSC) {
676
agg = !(MACRO_t1_flags(1 + k,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
677
|| MACRO_t1_flags(1 + k + 1,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
678
|| MACRO_t1_flags(1 + k + 2,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
679
|| (MACRO_t1_flags(1 + k + 3,1 + i)
680
& (~(T1_SIG_S | T1_SIG_SE | T1_SIG_SW | T1_SGN_S))) & (T1_SIG | T1_VISIT | T1_SIG_OTH));
682
agg = !(MACRO_t1_flags(1 + k,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
683
|| MACRO_t1_flags(1 + k + 1,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
684
|| MACRO_t1_flags(1 + k + 2,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH)
685
|| MACRO_t1_flags(1 + k + 3,1 + i) & (T1_SIG | T1_VISIT | T1_SIG_OTH));
691
mqc_setcurctx(mqc, T1_CTXNO_AGG);
692
if (!mqc_decode(mqc)) {
695
mqc_setcurctx(mqc, T1_CTXNO_UNI);
696
runlen = mqc_decode(mqc);
697
runlen = (runlen << 1) | mqc_decode(mqc);
701
for (j = k + runlen; j < k + 4 && j < t1->h; ++j) {
702
vsc = ((cblksty & J2K_CCP_CBLKSTY_VSC) && (j == k + 3 || j == t1->h - 1)) ? 1 : 0;
705
&t1->flags[((j+1) * t1->flags_stride) + i + 1],
706
&t1->data[(j * t1->w) + i],
709
agg && (j == k + runlen),
716
mqc_setcurctx(mqc, T1_CTXNO_UNI);
718
v = (v << 1) | mqc_decode(mqc);
719
v = (v << 1) | mqc_decode(mqc);
720
v = (v << 1) | mqc_decode(mqc);
723
opj_event_msg(t1->cinfo, EVT_WARNING, "Bad segmentation symbol %x\n", v);
727
} /* VSC and BYPASS by Antonin */
730
/** mod fixed_quality */
731
static double t1_getwmsedec(
741
double w1, w2, wmsedec;
743
w1 = (numcomps > 1) ? mct_getnorm(compno) : 1.0;
744
w2 = dwt_getnorm(level, orient);
745
} else { /* if (qmfbid == 0) */
746
w1 = (numcomps > 1) ? mct_getnorm_real(compno) : 1.0;
747
w2 = dwt_getnorm_real(level, orient);
749
wmsedec = w1 * w2 * stepsize * (1 << bpno);
750
wmsedec *= wmsedec * nmsedec / 8192.0;
755
static bool allocate_buffers(
763
if(datasize > t1->datasize){
764
opj_aligned_free(t1->data);
765
t1->data = (int*) opj_aligned_malloc(datasize * sizeof(int));
769
t1->datasize=datasize;
771
memset(t1->data,0,datasize * sizeof(int));
773
t1->flags_stride=w+2;
774
flagssize=t1->flags_stride * (h+2);
776
if(flagssize > t1->flagssize){
777
opj_aligned_free(t1->flags);
778
t1->flags = (flag_t*) opj_aligned_malloc(flagssize * sizeof(flag_t));
782
t1->flagssize=flagssize;
784
memset(t1->flags,0,flagssize * sizeof(flag_t));
792
/** mod fixed_quality */
793
static void t1_encode_cblk(
795
opj_tcd_cblk_enc_t* cblk,
803
opj_tcd_tile_t * tile)
805
double cumwmsedec = 0.0;
807
opj_mqc_t *mqc = t1->mqc; /* MQC component */
809
int passno, bpno, passtype;
812
char type = T1_TYPE_MQ;
816
for (i = 0; i < t1->w * t1->h; ++i) {
817
int tmp = abs(t1->data[i]);
818
max = int_max(max, tmp);
821
cblk->numbps = max ? (int_floorlog2(max) + 1) - T1_NMSEDEC_FRACBITS : 0;
823
bpno = cblk->numbps - 1;
826
mqc_resetstates(mqc);
827
mqc_setstate(mqc, T1_CTXNO_UNI, 0, 46);
828
mqc_setstate(mqc, T1_CTXNO_AGG, 0, 3);
829
mqc_setstate(mqc, T1_CTXNO_ZC, 0, 4);
830
mqc_init_enc(mqc, cblk->data);
832
for (passno = 0; bpno >= 0; ++passno) {
833
opj_tcd_pass_t *pass = &cblk->passes[passno];
835
type = ((bpno < (cblk->numbps - 4)) && (passtype < 2) && (cblksty & J2K_CCP_CBLKSTY_LAZY)) ? T1_TYPE_RAW : T1_TYPE_MQ;
839
t1_enc_sigpass(t1, bpno, orient, &nmsedec, type, cblksty);
842
t1_enc_refpass(t1, bpno, &nmsedec, type, cblksty);
845
t1_enc_clnpass(t1, bpno, orient, &nmsedec, cblksty);
846
/* code switch SEGMARK (i.e. SEGSYM) */
847
if (cblksty & J2K_CCP_CBLKSTY_SEGSYM)
848
mqc_segmark_enc(mqc);
853
tempwmsedec = t1_getwmsedec(nmsedec, compno, level, orient, bpno, qmfbid, stepsize, numcomps);
854
cumwmsedec += tempwmsedec;
855
tile->distotile += tempwmsedec;
857
/* Code switch "RESTART" (i.e. TERMALL) */
858
if ((cblksty & J2K_CCP_CBLKSTY_TERMALL) && !((passtype == 2) && (bpno - 1 < 0))) {
859
if (type == T1_TYPE_RAW) {
862
/* correction = mqc_bypass_flush_enc(); */
863
} else { /* correction = mqc_restart_enc(); */
869
if (((bpno < (cblk->numbps - 4) && (passtype > 0))
870
|| ((bpno == (cblk->numbps - 4)) && (passtype == 2))) && (cblksty & J2K_CCP_CBLKSTY_LAZY)) {
871
if (type == T1_TYPE_RAW) {
874
/* correction = mqc_bypass_flush_enc(); */
875
} else { /* correction = mqc_restart_enc(); */
885
if (++passtype == 3) {
890
if (pass->term && bpno > 0) {
891
type = ((bpno < (cblk->numbps - 4)) && (passtype < 2) && (cblksty & J2K_CCP_CBLKSTY_LAZY)) ? T1_TYPE_RAW : T1_TYPE_MQ;
892
if (type == T1_TYPE_RAW)
893
mqc_bypass_init_enc(mqc);
895
mqc_restart_init_enc(mqc);
898
pass->distortiondec = cumwmsedec;
899
pass->rate = mqc_numbytes(mqc) + correction; /* FIXME */
901
/* Code-switch "RESET" */
902
if (cblksty & J2K_CCP_CBLKSTY_RESET)
906
/* Code switch "ERTERM" (i.e. PTERM) */
907
if (cblksty & J2K_CCP_CBLKSTY_PTERM)
909
else /* Default coding */ if (!(cblksty & J2K_CCP_CBLKSTY_LAZY))
912
cblk->totalpasses = passno;
914
for (passno = 0; passno<cblk->totalpasses; passno++) {
915
opj_tcd_pass_t *pass = &cblk->passes[passno];
916
if (pass->rate > mqc_numbytes(mqc))
917
pass->rate = mqc_numbytes(mqc);
918
/*Preventing generation of FF as last data byte of a pass*/
919
if((pass->rate>1) && (cblk->data[pass->rate - 1] == 0xFF)){
922
pass->len = pass->rate - (passno == 0 ? 0 : cblk->passes[passno - 1].rate);
926
static void t1_decode_cblk(
928
opj_tcd_cblk_dec_t* cblk,
933
opj_raw_t *raw = t1->raw; /* RAW component */
934
opj_mqc_t *mqc = t1->mqc; /* MQC component */
938
char type = T1_TYPE_MQ; /* BYPASS mode */
940
if(!allocate_buffers(
943
cblk->y1 - cblk->y0))
948
bpno = roishift + cblk->numbps - 1;
951
mqc_resetstates(mqc);
952
mqc_setstate(mqc, T1_CTXNO_UNI, 0, 46);
953
mqc_setstate(mqc, T1_CTXNO_AGG, 0, 3);
954
mqc_setstate(mqc, T1_CTXNO_ZC, 0, 4);
956
for (segno = 0; segno < cblk->numsegs; ++segno) {
957
opj_tcd_seg_t *seg = &cblk->segs[segno];
960
type = ((bpno <= (cblk->numbps - 1) - 4) && (passtype < 2) && (cblksty & J2K_CCP_CBLKSTY_LAZY)) ? T1_TYPE_RAW : T1_TYPE_MQ;
961
/* FIXME: slviewer gets here with a null pointer. Why? Partially downloaded and/or corrupt textures? */
962
if(seg->data == NULL){
965
if (type == T1_TYPE_RAW) {
966
raw_init_dec(raw, (*seg->data) + seg->dataindex, seg->len);
968
mqc_init_dec(mqc, (*seg->data) + seg->dataindex, seg->len);
971
for (passno = 0; passno < seg->numpasses; ++passno) {
974
t1_dec_sigpass(t1, bpno+1, orient, type, cblksty);
977
t1_dec_refpass(t1, bpno+1, type, cblksty);
980
t1_dec_clnpass(t1, bpno+1, orient, cblksty);
984
if ((cblksty & J2K_CCP_CBLKSTY_RESET) && type == T1_TYPE_MQ) {
985
mqc_resetstates(mqc);
986
mqc_setstate(mqc, T1_CTXNO_UNI, 0, 46);
987
mqc_setstate(mqc, T1_CTXNO_AGG, 0, 3);
988
mqc_setstate(mqc, T1_CTXNO_ZC, 0, 4);
990
if (++passtype == 3) {
998
/* ----------------------------------------------------------------------- */
1000
opj_t1_t* t1_create(opj_common_ptr cinfo) {
1001
opj_t1_t *t1 = (opj_t1_t*) opj_malloc(sizeof(opj_t1_t));
1006
/* create MQC and RAW handles */
1007
t1->mqc = mqc_create();
1008
t1->raw = raw_create();
1018
void t1_destroy(opj_t1_t *t1) {
1020
/* destroy MQC and RAW handles */
1021
mqc_destroy(t1->mqc);
1022
raw_destroy(t1->raw);
1023
opj_aligned_free(t1->data);
1024
opj_aligned_free(t1->flags);
1029
void t1_encode_cblks(
1031
opj_tcd_tile_t *tile,
1034
int compno, resno, bandno, precno, cblkno;
1036
tile->distotile = 0; /* fixed_quality */
1038
for (compno = 0; compno < tile->numcomps; ++compno) {
1039
opj_tcd_tilecomp_t* tilec = &tile->comps[compno];
1040
opj_tccp_t* tccp = &tcp->tccps[compno];
1041
int tile_w = tilec->x1 - tilec->x0;
1043
for (resno = 0; resno < tilec->numresolutions; ++resno) {
1044
opj_tcd_resolution_t *res = &tilec->resolutions[resno];
1046
for (bandno = 0; bandno < res->numbands; ++bandno) {
1047
opj_tcd_band_t* restrict band = &res->bands[bandno];
1049
for (precno = 0; precno < res->pw * res->ph; ++precno) {
1050
opj_tcd_precinct_t *prc = &band->precincts[precno];
1052
for (cblkno = 0; cblkno < prc->cw * prc->ch; ++cblkno) {
1053
opj_tcd_cblk_enc_t* cblk = &prc->cblks.enc[cblkno];
1054
int* restrict datap;
1055
int* restrict tiledp;
1060
int x = cblk->x0 - band->x0;
1061
int y = cblk->y0 - band->y0;
1062
if (band->bandno & 1) {
1063
opj_tcd_resolution_t *pres = &tilec->resolutions[resno - 1];
1064
x += pres->x1 - pres->x0;
1066
if (band->bandno & 2) {
1067
opj_tcd_resolution_t *pres = &tilec->resolutions[resno - 1];
1068
y += pres->y1 - pres->y0;
1071
if(!allocate_buffers(
1073
cblk->x1 - cblk->x0,
1074
cblk->y1 - cblk->y0))
1083
tiledp=&tilec->data[(y * tile_w) + x];
1084
if (tccp->qmfbid == 1) {
1085
for (j = 0; j < cblk_h; ++j) {
1086
for (i = 0; i < cblk_w; ++i) {
1087
int tmp = tiledp[(j * tile_w) + i];
1088
datap[(j * cblk_w) + i] = tmp << T1_NMSEDEC_FRACBITS;
1091
} else { /* if (tccp->qmfbid == 0) */
1092
for (j = 0; j < cblk_h; ++j) {
1093
for (i = 0; i < cblk_w; ++i) {
1094
int tmp = tiledp[(j * tile_w) + i];
1095
datap[(j * cblk_w) + i] =
1098
8192 * 8192 / ((int) floor(band->stepsize * 8192))) >> (11 - T1_NMSEDEC_FRACBITS);
1108
tilec->numresolutions - 1 - resno,
1122
void t1_decode_cblks(
1124
opj_tcd_tilecomp_t* tilec,
1127
int resno, bandno, precno, cblkno;
1129
int tile_w = tilec->x1 - tilec->x0;
1131
for (resno = 0; resno < tilec->numresolutions; ++resno) {
1132
opj_tcd_resolution_t* res = &tilec->resolutions[resno];
1134
for (bandno = 0; bandno < res->numbands; ++bandno) {
1135
opj_tcd_band_t* restrict band = &res->bands[bandno];
1137
for (precno = 0; precno < res->pw * res->ph; ++precno) {
1138
opj_tcd_precinct_t* precinct = &band->precincts[precno];
1140
for (cblkno = 0; cblkno < precinct->cw * precinct->ch; ++cblkno) {
1141
opj_tcd_cblk_dec_t* cblk = &precinct->cblks.dec[cblkno];
1142
int* restrict datap;
1143
void* restrict tiledp;
1155
x = cblk->x0 - band->x0;
1156
y = cblk->y0 - band->y0;
1157
if (band->bandno & 1) {
1158
opj_tcd_resolution_t* pres = &tilec->resolutions[resno - 1];
1159
x += pres->x1 - pres->x0;
1161
if (band->bandno & 2) {
1162
opj_tcd_resolution_t* pres = &tilec->resolutions[resno - 1];
1163
y += pres->y1 - pres->y0;
1170
if (tccp->roishift) {
1171
int thresh = 1 << tccp->roishift;
1172
for (j = 0; j < cblk_h; ++j) {
1173
for (i = 0; i < cblk_w; ++i) {
1174
int val = datap[(j * cblk_w) + i];
1176
if (mag >= thresh) {
1177
mag >>= tccp->roishift;
1178
datap[(j * cblk_w) + i] = val < 0 ? -mag : mag;
1184
tiledp=(void*)&tilec->data[(y * tile_w) + x];
1185
if (tccp->qmfbid == 1) {
1186
for (j = 0; j < cblk_h; ++j) {
1187
for (i = 0; i < cblk_w; ++i) {
1188
int tmp = datap[(j * cblk_w) + i];
1189
((int*)tiledp)[(j * tile_w) + i] = tmp / 2;
1192
} else { /* if (tccp->qmfbid == 0) */
1193
for (j = 0; j < cblk_h; ++j) {
1194
for (i = 0; i < cblk_w; ++i) {
1195
float tmp = datap[(j * cblk_w) + i] * band->stepsize;
1196
((float*)tiledp)[(j * tile_w) + i] = tmp;
1200
opj_free(cblk->data);
1201
opj_free(cblk->segs);
1203
opj_free(precinct->cblks.dec);