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.c --- support and portability routines: error handling, safe memory
45 * Initial code by Alex Jakulin, Aug. 1995
47 * Modified and optimized: Gadiel Seroussi, October 1995 - ...
49 * Modified and added Restart marker and input tables by:
50 * David Cheng-Hsiu Chu, and Ismail R. Ismail march 1999
64 This program is Copyright (c) University of British Columbia.\n\
65 All rights reserved. It may be freely redistributed in its\n\
66 entirety provided that this copyright notice is not removed.\n\
67 It may not be sold for profit or incorporated in commercial programs\n\
68 without the written permission of the copyright holder.\n\
75 FILE *c_in[MAX_COMPONENTS];
76 FILE *c_out[MAX_COMPONENTS];
77 FILE *msgfile = NULL; /* = stdout;*/
79 /* Context quantization thresholds - initially unset */
86 int verbose = 1; /* verbosity level */
87 int nopause = 0; /* whether to pause the legal notice or not */
88 int nolegal = 0; /* whether to print the legal notice or not */
91 /* parameters for LOSSY images */
92 int quant, /* quantization = 2*NEAR+1 */
93 beta, /* size of extended alphabet */
94 qbeta, /* size of quantized alphabet */
95 ceil_half_qbeta, /* ceil(qbeta/2) */
97 alpha1eps; /* alpha-1+NEAR */
99 int NEAR = DEF_NEAR; /* loss tolerance per symbol, fixed at 0 for lossless */
100 int bpp, /* bits per sample */
101 qbpp, /* bits per sample for quantized prediction errors */
102 limit, /* limit for unary part of Golomb code */
103 limit_reduce; /* reduction on above for EOR states */
106 /* define color mode strings */
107 char *plane_int_string = "plane by plane",
108 *line_int_string = "line intlv",
109 *pixel_int_string = "sample intlv";
112 /* function to print out error messages */
113 void error(char *msg) {
114 fprintf(stderr, msg);
119 /* function to safely call malloc */
120 void *safealloc(size_t size) {
125 error("\nsafealloc: Out of memory. Aborting...\n");
130 /* function to safely call calloc **/
131 void *safecalloc(size_t numels, size_t size) {
134 temp = calloc(numels, size);
136 error("\nsafecalloc: Out of memory. Aborting...\n");
151 return (double)clock()/CLOCKS_PER_SEC;
155 /* Set thresholds to default unless specified by header: */
156 int set_thresholds(int alfa, int NEAR, int *T1p, int *T2p, int *T3p)
168 lambda = (alfa+127)/256;
170 lambda = (4096+127)/256;
175 /* compute lossless default */
177 T1 = lambda*(BASIC_T1 - 2) + 2;
178 else { /* alphabet < 8 bits */
179 T1 = BASIC_T1/ilambda;
180 if ( T1 < 2 ) T1 = 2;
182 /* adjust for lossy */
185 /* check that the default threshold is in bounds */
186 if ( T1 < NEAR+1 || T1 > (alfa-1) )
187 T1 = NEAR+1; /* eliminates the threshold */
190 /* compute lossless default */
192 T2 = lambda*(BASIC_T2 - 3) + 3;
194 T2 = BASIC_T2/ilambda;
195 if ( T2 < 3 ) T2 = 3;
197 /* adjust for lossy */
200 /* check that the default threshold is in bounds */
201 if ( T2 < T1 || T2 > (alfa-1) )
202 T2 = T1; /* eliminates the threshold */
205 /* compute lossless default */
207 T3 = lambda*(BASIC_T3 - 4) + 4;
209 T3 = BASIC_T3/ilambda;
210 if ( T3 < 4 ) T3 = 4;
212 /* adjust for lossy */
215 /* check that the default threshold is in bounds */
216 if ( T3 < T2 || T3 > (alfa-1) )
217 T3 = T2; /* eliminates the threshold */
229 /* We first check compatibility with JPEG-LS, then with this implementation */
231 void check_compatibility(jpeg_ls_header *head_frame, jpeg_ls_header *head_scan, int n_s)
234 int number_of_scans,i;
237 /* Check implemented color modes */
238 if ((head_scan->color_mode>PIXEL_INT)) {
239 fprintf(stderr,"Color mode %d not supported\n",head_scan->color_mode);
243 if (head_scan->color_mode==PLANE_INT)
244 number_of_scans=head_frame->comp;
249 /* Test standard compatibility */
251 if (head_frame->columns<=0 || head_frame->rows <=0) {
252 fprintf(stderr,"Image size must be positive for this implementation.\n");
256 if (head_frame->alp<4) {
257 fprintf(stderr,"Alphabet size must be >= 4, got %d\n",head_frame->alp);
262 if (head_scan->T1>head_scan->T2 || head_scan->T2>head_scan->T3 ||
263 head_scan->T1<head_scan->NEAR+1 || head_scan->T3>=head_scan->alp ) {
264 fprintf(stderr,"Bad thresholds: must be %d <= Ta <= Tb <= Tc <= %d\n",
265 head_scan->NEAR+1,head_scan->alp-1);
269 if (head_frame->comp>255) {
270 fprintf(stderr,"Too many components (must be less than 255)\n");
274 if (head_scan->NEAR>=head_scan->alp) {
275 fprintf(stderr,"Error for near-lossless must be smaller than alphabet (%d), got %d",head_scan->alp,head_scan->NEAR);
280 if (head_scan->RES < MINRESET || head_scan->RES >= head_scan->alp ) {
281 fprintf(stderr,"Reset parameter must be between %d and %d\n",
282 MINRESET, head_scan->alp-1);
287 maxreset = (head_scan->alp >= 256)? (head_scan->alp-1):255;
289 if (head_scan->RES < MINRESET || head_scan->RES > maxreset ) {
290 fprintf(stderr,"Reset parameter must be between %d and %d\n",
291 MINRESET, head_scan->alp-1);
295 for (i=0;i<head_frame->comp;i++)
296 if (head_frame->comp_ids[i] != (i+1)) {
297 fprintf(stderr,"Components id in frame not compatible with this implementation.\n");
301 if (number_of_scans == 1) {
302 if (head_frame->comp != head_scan->comp) {
303 fprintf(stderr,"In this implementation, when single scan, all components must be in the scan.\n");
306 for (i=0;i<head_frame->comp;i++)
307 if (head_scan->comp_ids[i] != (i+1)) {
308 fprintf(stderr,"Components id in single scan not compatible with this implementation.\n");
314 if (head_scan->comp != 1) {
315 fprintf(stderr,"Only 1 component per scan for plane interleaved mode\n");
318 if (head_scan->comp_ids[0] != (n_s+1)) {
319 fprintf(stderr,"Components id in multiple scan not compatible with this implementation.\n");
327 /* for writing disclaimer to command line in DOS */
329 char *ttyfilename = "CON";
333 void fprint_disclaimer(FILE *fp, int nopause)
339 nopause = nopause | !isatty(fileno(fp));
341 if ( !nopause && (ttyf=fopen(ttyfilename,"r"))==NULL ) {
345 for ( i=1, p0=disclaimer; ; i++ ) {
347 if ( !nopause && i%PAUSE==0 ) {
349 fprintf(stderr, "--- (press RETURN to continue) ---");
353 for ( p1=p0; (*p1 != '\n') && (*p1 != 0); p1++ );
355 fprintf(fp,"%s\n",p0);
358 fprintf(fp,"\n"); fflush(fp);
359 if ( !nopause) fclose(ttyf);