1 /* SPMG/JPEG-LS IMPLEMENTATION V.2.1
2 =====================================
3 These programs are Copyright (c) University of British Columbia. All rights reserved.
4 They may be freely redistributed in their entirety provided that this copyright
5 notice is not removed. They may not be sold for profit or incorporated in
6 commercial programs without the written permission of the copyright holder.
7 Each program is provided as is, without any express or implied warranty,
8 without even the warranty of fitness for a particular purpose.
10 =========================================================
11 THIS SOFTWARE IS BASED ON HP's implementation of jpeg-ls:
12 =========================================================
14 LOCO-I/JPEG-LS IMPLEMENTATION V.0.90
15 -------------------------------------------------------------------------------
16 (c) COPYRIGHT HEWLETT-PACKARD COMPANY, 1995-1997.
17 HEWLETT-PACKARD COMPANY ("HP") DOES NOT WARRANT THE ACCURACY OR
18 COMPLETENESS OF THE INFORMATION GIVEN HERE. ANY USE MADE OF, OR
19 RELIANCE ON, SUCH INFORMATION IS ENTIRELY AT USER'S OWN RISK.
20 BY DOWNLOADING THE LOCO-I/JPEG-LS COMPRESSORS/DECOMPRESSORS
21 ("THE SOFTWARE") YOU AGREE TO BE BOUND BY THE TERMS AND CONDITIONS
22 OF THIS LICENSING AGREEMENT.
23 YOU MAY DOWNLOAD AND USE THE SOFTWARE FOR NON-COMMERCIAL PURPOSES
24 FREE OF CHARGE OR FURTHER OBLIGATION. YOU MAY NOT, DIRECTLY OR
25 INDIRECTLY, DISTRIBUTE THE SOFTWARE FOR A FEE, INCORPORATE THIS
26 SOFTWARE INTO ANY PRODUCT OFFERED FOR SALE, OR USE THE SOFTWARE
27 TO PROVIDE A SERVICE FOR WHICH A FEE IS CHARGED.
28 YOU MAY MAKE COPIES OF THE SOFTWARE AND DISTRIBUTE SUCH COPIES TO
29 OTHER PERSONS PROVIDED THAT SUCH COPIES ARE ACCOMPANIED BY
30 HEWLETT-PACKARD'S COPYRIGHT NOTICE AND THIS AGREEMENT AND THAT
31 SUCH OTHER PERSONS AGREE TO BE BOUND BY THE TERMS OF THIS AGREEMENT.
32 THE SOFTWARE IS NOT OF PRODUCT QUALITY AND MAY HAVE ERRORS OR DEFECTS.
33 THE JPEG-LS STANDARD IS STILL UNDER DEVELOPMENT. THE SOFTWARE IS NOT A
34 FINAL OR FULL IMPLEMENTATION OF THE STANDARD. HP GIVES NO EXPRESS OR
35 IMPLIED WARRANTY OF ANY KIND AND ANY IMPLIED WARRANTIES OF
36 MERCHANTABILITY AND FITNESS FOR PURPOSE ARE DISCLAIMED.
37 HP SHALL NOT BE LIABLE FOR ANY DIRECT, INDIRECT, SPECIAL, INCIDENTAL,
38 OR CONSEQUENTIAL DAMAGES ARISING OUT OF ANY USE OF THE SOFTWARE.
39 -------------------------------------------------------------------------------
42 /* melcode.c --- for processing in run mode
44 * Initial code by Alex Jakulin, Aug. 1995
46 * Modified and optimized: Gadiel Seroussi, October 1995
48 * Color Enhancement: Guillermo Sapiro, August 1996
50 * Modified and added Restart marker and input tables by:
51 * David Cheng-Hsiu Chu, and Ismail R. Ismail march 1999
58 #define MELCSTATES 32 /* number of melcode states */
60 static int J[MELCSTATES] = {
61 0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,5,5,6,6,7,
62 7,8,9,10,11,12,13,14,15
66 static int melcstate[MAX_COMPONENTS], /* index to the state array */
67 melclen[MAX_COMPONENTS]; /* contents of the state array location
68 indexed by melcstate: the "expected"
69 run length is 2^melclen, shorter runs are
70 encoded by a 1 followed by the run length
71 in binary representation, wit a fixed length
74 static unsigned long melcorder[MAX_COMPONENTS]; /* 2^ melclen */
77 void init_process_run(int maxrun) /* maxrun is ignoreed when using MELCODE,
78 kept here for function compatibility */
83 for (n_c=0;n_c<components;n_c++)
87 melcorder[n_c] = 1<<melclen[n_c];
94 /* decoding routine: reads bits from the input and returns a run length. */
95 /* argument is the number of pixels left to end-of-line (bound on run length) */
96 int process_run_dec(int lineleft, int color)
101 register int temp, hits;
102 temp = zeroLUT[(byte)(~(reg >> 24))]; /* number of leading ones in the
103 input stream, up to 8 */
104 for ( hits = 1; hits<=temp; hits++ )
106 runlen += melcorder[color];
107 if ( runlen >= lineleft )
108 { /* reached end-of-line */
109 if ( runlen==lineleft && melcstate[color] < MELCSTATES )
111 melclen[color] = J[++melcstate[color]];
112 melcorder[color] = (1L<<melclen[color]);
114 FILLBUFFER(hits); /* actual # of 1's consumed */
117 if ( melcstate[color] < MELCSTATES )
119 melclen[color] = J[++melcstate[color]];
120 melcorder[color] = (1L<<melclen[color]);
125 FILLBUFFER(temp + 1); /* consume the leading
126 0 of the remainder encoding */
132 /* read the length of the remainder */
133 if ( melclen[color] )
136 GETBITS(temp, melclen[color]); /*** GETBITS is a macro, not a function */
139 limit_reduce = melclen[color]+1;
141 /* adjust melcoder parameters */
142 if ( melcstate[color] )
144 melclen[color] = J[--melcstate[color]];
145 melcorder[color] = (1L<<melclen[color]);
156 /* retained for compatibility with ranked runs */