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-1999.
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 /* global.h --- prototypes for functions and global variables
44 * Initial code by Alex Jakulin, Aug. 1995
46 * Modified and optimized: Gadiel Seroussi, October 1995
48 * Modified and added Restart marker and input tables by:
49 * David Cheng-Hsiu Chu, and Ismail R. Ismail march 1999
60 #ifndef CLOCKS_PER_SEC
61 #define CLOCKS_PER_SEC 10e6
70 /* TRUE and FALSE values */
76 #define JPEGLSVERSION "V.2.1"
79 /* Maximal number of components in the implementation*/
80 #define MAX_COMPONENTS 6
81 #define MAX_SCANS MAX_COMPONENTS
84 /* For 1st component of plane interl. mode */
88 /* Different colour modes */
93 #define DEFAULT_COLOR_MODE LINE_INT
95 extern char *plane_int_string,
102 typedef struct jpeg_ls {
104 int columns, /* The number of columns */
105 rows, /* Number of rows */
106 alp, /* alphabet size (Max+1) , 2 bytes*/
107 comp, /* number of components, 1 byte */
108 NEAR, /* near-lossless error, 1 byte */
109 color_mode, /* indicates the color mode , 1 byte */
110 need_lse, /* Indicates non-default parameters */
111 need_table, /* Indicates use of mapping table */
112 need_restart, /* Indicates use of restart markers */
113 restart_interval, /* The number of MCU's between restart markers */
114 shift, /* for sparse images, 1 byte */
115 T1,T2,T3, /* Thresholds, 2 bytes each */
116 RES, /* reset value for counters, 2 bytes */
117 samplingx[MAX_COMPONENTS], /* col. sampling rates 1 byte each*/
118 samplingy[MAX_COMPONENTS], /* row sampling rates */
119 comp_ids[MAX_COMPONENTS], /* component id's */
120 acc_size, /* 1 byte */
121 adds[MAX_COMPONENTS]; /* size given by acc_size */
122 unsigned int TID, /* Table ID, 1 byte */
123 Wt, /* Width of each table entry, 1 byte */
124 *TABLE[MAX_COMPONENTS]; /* The table(s) for each component */
128 extern int components;
129 extern int sampling[MAX_COMPONENTS];
131 #define NAME_LENGTH 40
133 /* Output file names */
134 #define OUTFILE "outfile"
135 #define COMPSUFFIX ".jls"
138 /* Define max and min macros */
140 # define max(a,b) (((a)>=(b))?(a):(b))
141 # define min(a,b) (((a)<=(b))?(a):(b))
147 /* margins for scan lines */
149 #define RIGHTMARGIN 1
152 extern char *disclaimer;
157 #define MAXA16 (65536)
158 #define LUTMAX8 (256)
159 #define LUTMAX16 (4501)
166 # define highmask (-(alpha))
170 # if (alpha!=2) && (alpha!=4) && (alpha!=8) && (alpha!=16) && (alpha!=32) &&\
171 (alpha!=64) && ( alpha!=128) && (alpha!=256) && (alpha!=512) &&\
172 (alpha!=1024) && ( alpha!=2048) && (alpha!=4096) && (alpha!=8192) &&\
173 (alpha!=16384) && ( alpha!=32768) && (alpha!=65536)\
174 # error "Fixed alpha must be a power of 2"
176 # define ceil_half_alpha (alpha/2)
178 extern int alpha; /* alphabet size */
179 extern int ceil_half_alpha; /* ceil(alpha/2) */
180 extern int highmask; /* for powers of 2, a mask for high bits */
185 extern int bpp, /* bits per sample */
186 qbpp, /* bits per sample for quantized prediction errors */
187 limit, /* limit for unary part of Golomb code */
188 limit_reduce; /* reduction on above for EOR states */
205 /* Quantization threshold basic defaults */
206 /* These are the defaults for LOSSLESS, 8 bpp. Defaults for other
207 cases are derived from these basic values */
213 #define CREGIONS (9) /* quantization regions for d-b, b-c, c-a */
215 /* run-length treshold */
223 /* number of different contexts */
224 #define CONTEXTS1 (CREGIONS*CREGIONS*CREGIONS)
226 #define CONTEXTS ((CONTEXTS1+1)/2) /* all regions, with symmetric merging */
229 /* Mandatory for JPEG-LS: */
239 #define MAXCODE (N_R_L_ERROR)
242 /* Number of end-of-run contexts */
243 #define EOR_CONTEXTS 2
246 /* Total number of contexts */
247 #define TOT_CONTEXTS (CONTEXTS + EOR_CONTEXTS)
250 /* index of first end-of-run context */
251 #define EOR_0 (CONTEXTS)
254 /* index of run state */
261 /* The longest code the bit IO can facilitate */
262 #define MAXCODELEN 24
264 /* The stat initialization values */
265 #define INITNSTAT 1 /* init value for N[] */
266 #define MIN_INITABSTAT 2 /* min init value for A[] */
267 #define INITABSLACK 6 /* init value for A is roughly
268 2^(bpp-INITABSLACK) but not less than above */
269 #define INITBIASTAT 0 /* init value for B[] */
271 /* Limit for unary code */
275 #define DEFAULT_RESET 64
280 # define RESET DEFAULT_RESET
286 #define reset RESET /* reset threshold */
291 /****** Type prototypes */
293 /* Portability types */
294 typedef unsigned char byte;
295 typedef unsigned short word;
296 typedef unsigned long dword;
298 typedef unsigned short pixel;
301 /****** Global variables prototypes */
303 extern FILE *in, *out, *msgfile;
304 extern FILE *c_in[MAX_COMPONENTS];
305 extern FILE *c_out[MAX_COMPONENTS];
307 extern int T1, T2, T3, Ta;
312 extern int bpp16; /* Indicates if 16 bits per pixel mode or not */
316 /* for look-up-tables */
318 extern int vLUT[3][2 * LUTMAX16];
320 extern int classmap[CONTEXTS1];
321 extern int *qdiv0, *qdiv, /* quantization table (division via look-up) */
322 *qmul0, *qmul; /* dequantization table */
324 /* statistics tables */
325 extern int N[TOT_CONTEXTS],
331 /*extern byte getk[65][3000];*/
332 /*extern int clipPx[510];*/
335 /****** Function prototypes */
338 void error(char *msg);
339 void *safealloc(size_t size);
340 void *safecalloc(size_t numels, size_t size);
343 void lossless_doscanline_pixel(pixel *psl, pixel *sl, int no);
344 void lossless_doscanline(pixel *psl, pixel *sl, int no, int color);
347 void lossy_doscanline_pixel(pixel *psl, pixel *sl, int no);
348 void lossy_doscanline(pixel *psl, pixel *sl, int no, int color);
353 void buffinit(FILE *);
356 void init_process_run(int);
357 void close_process_run();
358 int process_run(int,int,int);
363 void prepare_qtables(int, int);
364 void init_stats(int);
369 # define ENDIAN8(x) (x)
370 # define ENDIAN16(x) (x)
372 # define ENDIAN8(x) (x&0x000000ff)
373 # define ENDIAN16(x) ( ((x>>8)|(x<<8)) & 0x0000ffff)
376 /* ENDIAN function to fix endian of PCs (for 8 bit pixels)
377 #define ENDIAN8(x) (x&0x000000ff)*/
380 /* ENDIAN function to fix endian of PCs (for 16 bit pixels)
381 #define ENDIAN16(x) ( ((x>>8)|(x<<8)) & 0x0000ffff )*/
387 # define clip(x,alpha) \
388 if ( x & highmask ) {\
395 # define clip(x,alpha) \
398 else if (x >= alpha) \
404 /* macro to predict Px */
405 #define predict(Rb, Ra, Rc) \
407 register pixel minx; \
408 register pixel maxx; \
419 else if (Rc <= minx) \