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
55 #include "jpegls_config.h"
60 #ifndef CLOCKS_PER_SEC
61 #define CLOCKS_PER_SEC 10e6
70 /* TRUE and FALSE values */
76 #define JPEGLSVERSION "V.2.1"
78 /* Maximal number of components in the implementation*/
79 #define MAX_COMPONENTS 6
80 #define MAX_SCANS MAX_COMPONENTS
83 /* For 1st component of plane interl. mode */
87 /* Different colour modes */
92 #define DEFAULT_COLOR_MODE LINE_INT
94 extern char *plane_int_string,
103 typedef struct jpeg_ls {
105 int columns, /* The number of columns */
106 rows, /* Number of rows */
107 alp, /* alphabet size (Max+1) , 2 bytes*/
108 comp, /* number of components, 1 byte */
109 NEAR, /* near-lossless error, 1 byte */
110 color_mode, /* indicates the color mode , 1 byte */
111 need_lse, /* Indicates non-default parameters */
112 need_table, /* Indicates use of mapping table */
113 need_restart, /* Indicates use of restart markers */
114 restart_interval, /* The number of MCU's between restart markers */
115 shift, /* for sparse images, 1 byte */
116 T1,T2,T3, /* Thresholds, 2 bytes each */
117 RES, /* reset value for counters, 2 bytes */
118 samplingx[MAX_COMPONENTS], /* col. sampling rates 1 byte each*/
119 samplingy[MAX_COMPONENTS], /* row sampling rates */
120 comp_ids[MAX_COMPONENTS], /* component id's */
121 acc_size, /* 1 byte */
122 adds[MAX_COMPONENTS]; /* size given by acc_size */
123 unsigned int TID, /* Table ID, 1 byte */
124 MAXTAB, /* Maximum table index value */
125 Wt, /* Width of each table entry, 1 byte */
126 *TABLE[MAX_COMPONENTS]; /* The table(s) for each component */
130 extern int components;
131 extern int sampling[MAX_COMPONENTS];
133 #define NAME_LENGTH 40
135 /* Output file names */
136 #define OUTFILE "outfile"
137 #define COMPSUFFIX ".out"
140 /* Define max and min macros */
142 #define max(a,b) (((a)>=(b))?(a):(b))
143 #define min(a,b) (((a)<=(b))?(a):(b))
149 /* margins for scan lines */
151 #define RIGHTMARGIN 1
154 extern char *disclaimer;
159 #define MAXA16 (65536)
160 #define LUTMAX8 (256)
161 #define LUTMAX16 (4501)
168 # define highmask (-(alpha))
172 # if (alpha!=2) && (alpha!=4) && (alpha!=8) && (alpha!=16) && (alpha!=32) &&\
173 (alpha!=64) && ( alpha!=128) && (alpha!=256) && (alpha!=512) &&\
174 (alpha!=1024) && ( alpha!=2048) && (alpha!=4096) && (alpha!=8192) &&\
175 (alpha!=16384) && ( alpha!=32768) && (alpha!=65536)
176 # error "Fixed alpha must be a power of 2"
178 # define ceil_half_alpha (alpha/2)
180 extern int alpha; /* alphabet size */
181 extern int ceil_half_alpha; /* ceil(alpha/2) */
182 extern int highmask; /* for powers of 2, a mask for high bits */
187 extern int bpp, /* bits per sample */
188 qbpp, /* bits per sample for quantized prediction errors */
189 limit, /* limit for unary part of Golomb code */
190 limit_reduce; /* reduction on above for EOR states */
207 /* Quantization threshold basic defaults */
208 /* These are the defaults for LOSSLESS, 8 bpp. Defaults for other
209 cases are derived from these basic values */
215 #define CREGIONS (9) /* quantization regions for d-a, a-c, c-b */
217 /* run-length treshold */
225 /* number of different contexts */
226 #define CONTEXTS1 (CREGIONS*CREGIONS*CREGIONS)
228 #define CONTEXTS ((CONTEXTS1+1)/2) /* all regions, with symmetric merging */
231 /* Mandatory for JPEG-LS: */
241 #define MAXCODE (N_R_L_ERROR)
244 /* Number of end-of-run contexts */
245 #define EOR_CONTEXTS 2
248 /* Total number of contexts */
249 #define TOT_CONTEXTS (CONTEXTS + EOR_CONTEXTS)
252 /* index of first end-of-run context */
253 #define EOR_0 (CONTEXTS)
256 /* index of run state */
263 /* The longest code the bit IO can facilitate */
264 #define MAXCODELEN 24
266 /* The stat initialization values */
268 #define MIN_INITABSTAT 2 /* min init value for abstat[] */
269 #define INITABSLACK 6 /* init value for abstat is roughly
270 2^(bpp-INITABSLACK) but not less than above */
271 #define INITBIASTAT 0
273 /* Limit for unary code */
277 #define DEFAULT_RESET 64
282 # define RESET DEFAULT_RESET
288 #define reset RESET /* reset threshold */
293 /****** Type prototypes */
295 /* Portability types */
296 typedef unsigned char byte;
297 typedef unsigned short word;
298 typedef unsigned long dword;
300 typedef unsigned short pixel;
304 /****** Global variables prototypes */
306 extern FILE *in, *out, *msgfile;
307 extern FILE *c_in[MAX_COMPONENTS];
308 extern FILE *c_out[MAX_COMPONENTS];
310 extern int T1, T2, T3, Ta;
319 /* for look-up tables */
321 extern int vLUT[3][2 * LUTMAX16];
323 extern int classmap[CONTEXTS1];
324 extern int *qdiv0, *qdiv, /* quantization table (division via look-up) */
325 *qmul0, *qmul; /* dequantization table */
327 /* statistics tables */
328 extern int N[TOT_CONTEXTS],
334 /****** Function prototypes */
337 void error(char *msg);
338 void *safealloc(size_t size);
339 void *safecalloc(size_t numels, size_t size);
341 int set_thresholds(int alfa, int NEAR, int *T1p, int *T2p, int *T3p);
342 void check_compatibility(jpeg_ls_header *head_frame, jpeg_ls_header *head_scan, int n_s);
346 void init_stats(int);
347 void doscanline(pixel *psl, pixel *sl, int no, int color);
348 int undoscanline(pixel *psl, pixel *sl, int no, int color);
349 void doscanline_pixel(pixel *psl, pixel *sl, int no);
350 int undoscanline_pixel(pixel *psl, pixel *sl, int no);
355 void createzeroLUT();
356 void buffinit(FILE *);
357 void bufiinit(FILE *fil);
361 void init_process_run(int);
362 void close_process_run();
363 int process_run_dec(int,int);
367 void prepare_qtables(int, int);
368 void init_stats(int);
371 int lossless_undoscanline( pixel *psl, pixel *sl, int no, int color);
372 int lossless_undoscanline_pixel(pixel *psl, pixel *sl, int no);
375 int lossy_undoscanline( pixel *psl, pixel *sl, int no, int color);
376 int lossy_undoscanline_pixel( pixel *psl, pixel *sl, int no);
380 # define ENDIAN8(x) (x)
381 # define ENDIAN16(x) (x)
383 # define ENDIAN8(x) (x&0x000000ff)
384 # define ENDIAN16(x) ( ((x>>8)|(x<<8)) & 0x0000ffff)
387 /* ENDIAN function to fix endian of PCs (for 8 bit pixels)
388 #define ENDIAN8(x) (x&0x000000ff)*/
391 /* ENDIAN function to fix endian of PCs (for 16 bit pixels)
392 #define ENDIAN16(x) ( ((x>>8)|(x<<8)) & 0x0000ffff )*/
399 # define clip(x,alpha) \
400 if ( x & highmask ) {\
407 # define clip(x,alpha) \
410 else if (x >= alpha) \
416 /* macro to predict Px */
417 #define predict(Rb, Ra, Rc) \
419 register pixel minx; \
420 register pixel maxx; \
431 else if (Rc <= minx) \