]> Creatis software - gdcm.git/blob - src/gdcmjasper/src/libjasper/pgx/pgx_enc.c
ENH: Ok since OJ warnings are not going to be fixed anytime soon remove them
[gdcm.git] / src / gdcmjasper / src / libjasper / pgx / pgx_enc.c
1 /*
2  * Copyright (c) 2001-2003 Michael David Adams.
3  * All rights reserved.
4  */
5
6 /* __START_OF_JASPER_LICENSE__
7  * 
8  * JasPer License Version 2.0
9  * 
10  * Copyright (c) 1999-2000 Image Power, Inc.
11  * Copyright (c) 1999-2000 The University of British Columbia
12  * Copyright (c) 2001-2003 Michael David Adams
13  * 
14  * All rights reserved.
15  * 
16  * Permission is hereby granted, free of charge, to any person (the
17  * "User") obtaining a copy of this software and associated documentation
18  * files (the "Software"), to deal in the Software without restriction,
19  * including without limitation the rights to use, copy, modify, merge,
20  * publish, distribute, and/or sell copies of the Software, and to permit
21  * persons to whom the Software is furnished to do so, subject to the
22  * following conditions:
23  * 
24  * 1.  The above copyright notices and this permission notice (which
25  * includes the disclaimer below) shall be included in all copies or
26  * substantial portions of the Software.
27  * 
28  * 2.  The name of a copyright holder shall not be used to endorse or
29  * promote products derived from the Software without specific prior
30  * written permission.
31  * 
32  * THIS DISCLAIMER OF WARRANTY CONSTITUTES AN ESSENTIAL PART OF THIS
33  * LICENSE.  NO USE OF THE SOFTWARE IS AUTHORIZED HEREUNDER EXCEPT UNDER
34  * THIS DISCLAIMER.  THE SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS
35  * "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING
36  * BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
37  * PARTICULAR PURPOSE AND NONINFRINGEMENT OF THIRD PARTY RIGHTS.  IN NO
38  * EVENT SHALL THE COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, OR ANY SPECIAL
39  * INDIRECT OR CONSEQUENTIAL DAMAGES, OR ANY DAMAGES WHATSOEVER RESULTING
40  * FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT,
41  * NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
42  * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.  NO ASSURANCES ARE
43  * PROVIDED BY THE COPYRIGHT HOLDERS THAT THE SOFTWARE DOES NOT INFRINGE
44  * THE PATENT OR OTHER INTELLECTUAL PROPERTY RIGHTS OF ANY OTHER ENTITY.
45  * EACH COPYRIGHT HOLDER DISCLAIMS ANY LIABILITY TO THE USER FOR CLAIMS
46  * BROUGHT BY ANY OTHER ENTITY BASED ON INFRINGEMENT OF INTELLECTUAL
47  * PROPERTY RIGHTS OR OTHERWISE.  AS A CONDITION TO EXERCISING THE RIGHTS
48  * GRANTED HEREUNDER, EACH USER HEREBY ASSUMES SOLE RESPONSIBILITY TO SECURE
49  * ANY OTHER INTELLECTUAL PROPERTY RIGHTS NEEDED, IF ANY.  THE SOFTWARE
50  * IS NOT FAULT-TOLERANT AND IS NOT INTENDED FOR USE IN MISSION-CRITICAL
51  * SYSTEMS, SUCH AS THOSE USED IN THE OPERATION OF NUCLEAR FACILITIES,
52  * AIRCRAFT NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL
53  * SYSTEMS, DIRECT LIFE SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH
54  * THE FAILURE OF THE SOFTWARE OR SYSTEM COULD LEAD DIRECTLY TO DEATH,
55  * PERSONAL INJURY, OR SEVERE PHYSICAL OR ENVIRONMENTAL DAMAGE ("HIGH
56  * RISK ACTIVITIES").  THE COPYRIGHT HOLDERS SPECIFICALLY DISCLAIM ANY
57  * EXPRESS OR IMPLIED WARRANTY OF FITNESS FOR HIGH RISK ACTIVITIES.
58  * 
59  * __END_OF_JASPER_LICENSE__
60  */
61
62 /******************************************************************************\
63 * Includes.
64 \******************************************************************************/
65
66 #include <assert.h>
67
68 #include "jasper/jas_tvp.h"
69 #include "jasper/jas_stream.h"
70 #include "jasper/jas_image.h"
71 #include "jasper/jas_string.h"
72 #include "jasper/jas_debug.h"
73
74 #include "pgx_cod.h"
75 #include "pgx_enc.h"
76
77 /******************************************************************************\
78 * Local functions.
79 \******************************************************************************/
80
81 static int pgx_puthdr(jas_stream_t *out, pgx_hdr_t *hdr);
82 static int pgx_putdata(jas_stream_t *out, pgx_hdr_t *hdr, jas_image_t *image, int cmpt);
83 static int pgx_putword(jas_stream_t *out, bool bigendian, int prec,
84   uint_fast32_t val);
85 static uint_fast32_t pgx_inttoword(int_fast32_t val, int prec, bool sgnd);
86
87 /******************************************************************************\
88 * Code for save operation.
89 \******************************************************************************/
90
91 /* Save an image to a stream in the the PGX format. */
92
93 int pgx_encode(jas_image_t *image, jas_stream_t *out, char *optstr)
94 {
95   pgx_hdr_t hdr;
96   uint_fast32_t width;
97   uint_fast32_t height;
98   bool sgnd;
99   int prec;
100   pgx_enc_t encbuf;
101   pgx_enc_t *enc = &encbuf;
102
103   /* Avoid compiler warnings about unused parameters. */
104   optstr = 0;
105
106   switch (jas_clrspc_fam(jas_image_clrspc(image))) {
107   case JAS_CLRSPC_FAM_GRAY:
108     if ((enc->cmpt = jas_image_getcmptbytype(image,
109       JAS_IMAGE_CT_COLOR(JAS_CLRSPC_CHANIND_GRAY_Y))) < 0) {
110       jas_eprintf("error: missing color component\n");
111       return -1;
112     }
113     break;
114   default:
115     jas_eprintf("error: BMP format does not support color space\n");
116     return -1;
117     break;
118   }
119
120   width = jas_image_cmptwidth(image, enc->cmpt);
121   height = jas_image_cmptheight(image, enc->cmpt);
122   prec = jas_image_cmptprec(image, enc->cmpt);
123   sgnd = jas_image_cmptsgnd(image, enc->cmpt);
124
125   /* The PGX format is quite limited in the set of image geometries
126     that it can handle.  Here, we check to ensure that the image to
127     be saved can actually be represented reasonably accurately using the
128     PGX format. */
129   /* There must be exactly one component. */
130   if (jas_image_numcmpts(image) > 1 || prec > 16) {
131     fprintf(stderr, "The PNM format cannot be used to represent an image with this geometry.\n");
132     return -1;
133   }
134
135   hdr.magic = PGX_MAGIC;
136   hdr.bigendian = true;
137   hdr.sgnd = sgnd;
138   hdr.prec = prec;
139   hdr.width = width;
140   hdr.height = height;
141
142 #ifdef PGX_DEBUG
143   pgx_dumphdr(stderr, &hdr);
144 #endif
145
146   if (pgx_puthdr(out, &hdr)) {
147     return -1;
148   }
149
150   if (pgx_putdata(out, &hdr, image, enc->cmpt)) {
151     return -1;
152   }
153
154   return 0;
155 }
156
157 /******************************************************************************\
158 \******************************************************************************/
159
160 static int pgx_puthdr(jas_stream_t *out, pgx_hdr_t *hdr)
161 {
162   jas_stream_printf(out, "%c%c", hdr->magic >> 8, hdr->magic & 0xff);
163   jas_stream_printf(out, " %s %s %d %ld %ld\n", hdr->bigendian ? "ML" : "LM",
164     hdr->sgnd ? "-" : "+", hdr->prec, (long) hdr->width, (long) hdr->height);
165   if (jas_stream_error(out)) {
166     return -1;
167   }
168   return 0;
169 }
170
171 static int pgx_putdata(jas_stream_t *out, pgx_hdr_t *hdr, jas_image_t *image, int cmpt)
172 {
173   jas_matrix_t *data;
174   uint_fast32_t x;
175   uint_fast32_t y;
176   int_fast32_t v;
177   uint_fast32_t word;
178
179   data = 0;
180
181   if (!(data = jas_matrix_create(1, hdr->width))) {
182     goto error;
183   }
184   for (y = 0; y < hdr->height; ++y) {
185     if (jas_image_readcmpt(image, cmpt, 0, y, hdr->width, 1, data)) {
186       goto error;
187     }
188     for (x = 0; x < hdr->width; ++x) {
189       v = jas_matrix_get(data, 0, x);
190       word = pgx_inttoword(v, hdr->prec, hdr->sgnd);
191       if (pgx_putword(out, hdr->bigendian, hdr->prec, word)) {
192         goto error;
193       }
194     }
195   }
196   jas_matrix_destroy(data);
197   data = 0;
198   return 0;
199
200 error:
201   if (data) {
202     jas_matrix_destroy(data);
203   }
204   return -1;
205 }
206
207 static int pgx_putword(jas_stream_t *out, bool bigendian, int prec,
208   uint_fast32_t val)
209 {
210   int i;
211   int j;
212   int wordsize;
213
214   val &= (1 << prec) - 1;
215   wordsize = (prec + 7) /8;
216   for (i = 0; i < wordsize; ++i) {
217     j = bigendian ? (wordsize - 1 - i) : i;
218     if (jas_stream_putc(out, (val >> (8 * j)) & 0xff) == EOF) {
219       return -1;
220     }
221   }
222   return 0;
223 }
224
225 static uint_fast32_t pgx_inttoword(jas_seqent_t v, int prec, bool sgnd)
226 {
227   uint_fast32_t ret;
228   ret = ((sgnd && v < 0) ? ((1 << prec) + v) : v) & ((1 << prec) - 1);
229   return ret;
230 }