1 /* Copyright (c) 2001 David Janssens
2 * Copyright (c) 2002-2003 Yannick Verschueren
3 * Copyright (c) 2002-2003 Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
9 * 1. Redistributions of source code must retain the above copyright
10 * notice, this list of conditions and the following disclaimer.
11 * 2. Redistributions in binary form must reproduce the above copyright
12 * notice, this list of conditions and the following disclaimer in the
13 * documentation and/or other materials provided with the distribution.
15 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
16 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
19 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
20 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
21 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
22 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
23 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
24 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
25 * POSSIBILITY OF SUCH DAMAGE.
36 int ceildiv(int a, int b)
38 return (a + b - 1) / b;
41 int main(int argc, char **argv)
45 char *dest, S1, S2, S3;
50 int w, wr, wrr, h, hr, hrr, max;
51 int i, image_type = -1, compno, pad;
56 "usage: %s j2k-file image-file -reduce n (<- optional)\n",
61 f = fopen(argv[1], "rb");
63 fprintf(stderr, "failed to open %s for reading\n", argv[1]);
70 option.reduce_value = 0;
72 /* OPTION REDUCE IS ACTIVE */
74 if (strcmp(argv[3], "-reduce")) {
76 "usage: options " "-reduce n"
77 " where n is the factor of reduction [%s]\n", argv[3]);
81 sscanf(argv[4], "%d", &option.reduce_value);
94 if ((S1 == 'p' && S2 == 'g' && S3 == 'x')
95 || (S1 == 'P' && S2 == 'G' && S3 == 'X')) {
99 if ((S1 == 'p' && S2 == 'n' && S3 == 'm')
100 || (S1 == 'P' && S2 == 'N' && S3 == 'M') || (S1 == 'p' && S2 == 'g'
102 || (S1 == 'P' && S2 == 'G' && S3 == 'M') || (S1 == 'P' && S2 == 'P'
104 || (S1 == 'p' && S2 == 'p' && S3 == 'm')) {
108 if ((S1 == 'b' && S2 == 'm' && S3 == 'p')
109 || (S1 == 'B' && S2 == 'M' && S3 == 'P')) {
113 if (image_type == -1) {
115 "!! Unrecognized format for infile : %c%c%c [accept only *.pnm, *.pgm, *.ppm, *.pgx or *.bmp] !!\n\n",
120 fseek(f, 0, SEEK_END);
122 fseek(f, 0, SEEK_SET);
123 src = (char *) malloc(len);
124 fread(src, 1, len, f);
138 if (S1 == 'j' && S2 == '2' && S3 == 'k') {
139 if (!j2k_decode(src, len, &img, &cp, option)) {
140 fprintf(stderr, "j2k_to_image: failed to decode image!\n");
144 if (S1 == 'j' && S2 == 'p' && S3 == 't') {
145 if (!j2k_decode_jpt_stream(src, len, &img, &cp)) {
146 fprintf(stderr, "j2k_to_image: failed to decode image!\n");
151 "j2k_to_image : Unknown format image *.%c%c%c [only *.j2k or *.jpt]!! \n",
158 /* ------------------ CREATE OUT IMAGE WITH THE RIGHT FORMAT ----------------------- */
160 /* ---------------------------- / */
162 /* / FORMAT : PNM, PGM or PPM / */
164 /* ---------------------------- / */
166 switch (image_type) {
167 case 1: /* PNM PGM PPM */
168 if (img->numcomps == 3 && img->comps[0].dx == img->comps[1].dx
169 && img->comps[1].dx == img->comps[2].dx
170 && img->comps[0].dy == img->comps[1].dy
171 && img->comps[1].dy == img->comps[2].dy
172 && img->comps[0].prec == img->comps[1].prec
173 && img->comps[1].prec == img->comps[2].prec) {
174 f = fopen(argv[2], "wb");
175 w = ceildiv(img->x1 - img->x0, img->comps[0].dx);
176 // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor),img->comps[0].dx);
177 wr = img->comps[0].w;
178 wrr = int_ceildivpow2(img->comps[0].w, img->comps[0].factor);
180 h = ceildiv(img->y1 - img->y0, img->comps[0].dy);
181 // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), img->comps[0].dy);
182 hr = img->comps[0].h;
183 hrr = int_ceildivpow2(img->comps[0].h, img->comps[0].factor);
185 max = img->comps[0].prec > 8 ? 255 : (1 << img->comps[0].prec) - 1;
188 int_ceildivpow2(img->comps[0].x0 -
189 int_ceildiv(img->x0, img->comps[0].dx),
190 img->comps[0].factor);
192 int_ceildivpow2(img->comps[0].y0 -
193 int_ceildiv(img->y0, img->comps[0].dy),
194 img->comps[0].factor);
197 fprintf(f, "P6\n# %d %d %d %d %d\n%d %d\n%d\n",
198 cp->tcps[cp->tileno[0]].tccps[0].numresolutions, w, h,
199 img->comps[0].x0, img->comps[0].y0, wrr, hrr, max);
200 adjust = img->comps[0].prec > 8 ? img->comps[0].prec - 8 : 0;
201 for (i = 0; i < wrr * hrr; i++) {
203 r = img->comps[0].data[i / wrr * wr + i % wrr];
204 r += (img->comps[0].sgnd ? 1 << (img->comps[0].prec - 1) : 0);
207 g = img->comps[1].data[i / wrr * wr + i % wrr];
208 g += (img->comps[1].sgnd ? 1 << (img->comps[1].prec - 1) : 0);
211 b = img->comps[2].data[i / wrr * wr + i % wrr];
212 b += (img->comps[2].sgnd ? 1 << (img->comps[2].prec - 1) : 0);
215 fprintf(f, "%c%c%c", r, g, b);
219 for (compno = 0; compno < img->numcomps; compno++) {
221 if (img->numcomps > 1) {
222 sprintf(name, "%d.%s", compno, argv[2]);
224 sprintf(name, "%s", argv[2]);
226 f = fopen(name, "wb");
227 w = ceildiv(img->x1 - img->x0, img->comps[compno].dx);
228 // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor),img->comps[compno].dx);
229 wr = img->comps[compno].w;
231 int_ceildivpow2(img->comps[compno].w, img->comps[compno].factor);
233 h = ceildiv(img->y1 - img->y0, img->comps[compno].dy);
234 // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), img->comps[compno].dy);
235 hr = img->comps[compno].h;
237 int_ceildivpow2(img->comps[compno].h, img->comps[compno].factor);
240 img->comps[compno].prec >
241 8 ? 255 : (1 << img->comps[compno].prec) - 1;
243 img->comps[compno].x0 =
244 int_ceildivpow2(img->comps[compno].x0 -
246 img->comps[compno].dx),
247 img->comps[compno].factor);
248 img->comps[compno].y0 =
249 int_ceildivpow2(img->comps[compno].y0 -
251 img->comps[compno].dy),
252 img->comps[compno].factor);
254 fprintf(f, "P5\n# %d %d %d %d %d\n%d %d\n%d\n",
255 cp->tcps[cp->tileno[0]].tccps[compno].
256 numresolutions, w, h, img->comps[compno].x0,
257 img->comps[compno].y0, wrr, hrr, max);
259 img->comps[compno].prec > 8 ? img->comps[compno].prec - 8 : 0;
260 for (i = 0; i < wrr * hrr; i++) {
262 l = img->comps[compno].data[i / wrr * wr + i % wrr];
263 l += (img->comps[compno].
264 sgnd ? 1 << (img->comps[compno].prec - 1) : 0);
273 /* ------------------------ / */
275 /* / FORMAT : PGX / */
277 /* /----------------------- / */
279 for (compno = 0; compno < img->numcomps; compno++) {
280 j2k_comp_t *comp = &img->comps[compno];
282 if (img->numcomps > 1)
283 sprintf(name, "%d_%s", compno, argv[2]);
285 sprintf(name, "%s", argv[2]);
287 f = fopen(name, "wb");
288 // w = ceildiv(img->x1 - img->x0, comp->dx);
289 // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor), comp->dx);
290 w = img->comps[compno].w;
291 wr = int_ceildivpow2(img->comps[compno].w,
292 img->comps[compno].factor);
294 // h = ceildiv(img->y1 - img->y0, comp->dy);
295 // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), comp->dy);
296 h = img->comps[compno].h;
297 hr = int_ceildivpow2(img->comps[compno].h,
298 img->comps[compno].factor);
300 fprintf(f, "PG LM %c %d %d %d\n", comp->sgnd ? '-' : '+',
302 for (i = 0; i < wr * hr; i++) {
303 int v = img->comps[compno].data[i / wr * w + i % wr];
304 if (comp->prec <= 8) {
307 } else if (comp->prec <= 16) {
318 /* ------------------------ / */
320 /* / FORMAT : BMP / */
322 /* /----------------------- / */
325 if (img->numcomps == 3 && img->comps[0].dx == img->comps[1].dx
326 && img->comps[1].dx == img->comps[2].dx
327 && img->comps[0].dy == img->comps[1].dy
328 && img->comps[1].dy == img->comps[2].dy
329 && img->comps[0].prec == img->comps[1].prec
330 && img->comps[1].prec == img->comps[2].prec) {
331 /* -->> -->> -->> -->>
335 <<-- <<-- <<-- <<-- */
337 f = fopen(argv[2], "wb");
338 // w = ceildiv(img->x1 - img->x0, img->comps[0].dx);
339 // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor), img->comps[0].dx);
341 wr = int_ceildivpow2(img->comps[0].w, img->comps[0].factor);
343 // h = ceildiv(img->y1 - img->y0, img->comps[0].dy);
344 // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), img->comps[0].dy);
346 hr = int_ceildivpow2(img->comps[0].h, img->comps[0].factor);
352 fprintf(f, "%c%c%c%c",
353 (unsigned char) (hr * wr * 3 + 3 * hr * (wr % 2) +
355 (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
357 (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
359 (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
361 fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
362 ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
363 fprintf(f, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff,
364 ((54) >> 16) & 0xff, ((54) >> 24) & 0xff);
368 fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,
369 ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
370 fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
371 (unsigned char) ((wr) >> 8) & 0xff,
372 (unsigned char) ((wr) >> 16) & 0xff,
373 (unsigned char) ((wr) >> 24) & 0xff);
374 fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
375 (unsigned char) ((hr) >> 8) & 0xff,
376 (unsigned char) ((hr) >> 16) & 0xff,
377 (unsigned char) ((hr) >> 24) & 0xff);
378 fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
379 fprintf(f, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff);
380 fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
381 ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
382 fprintf(f, "%c%c%c%c",
383 (unsigned char) (3 * hr * wr +
384 3 * hr * (wr % 2)) & 0xff,
385 (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
387 (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
389 (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
391 fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
392 ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
393 fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
394 ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
395 fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
396 ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
397 fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
398 ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
400 for (i = 0; i < wr * hr; i++) {
401 unsigned char R, G, B;
403 // R = img->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
404 R = img->comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
405 // G = img->comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
406 G = img->comps[1].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
407 // B = img->comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
408 B = img->comps[2].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
409 fprintf(f, "%c%c%c", B, G, R);
411 if ((i + 1) % wr == 0) {
412 for (pad = (3 * wr) % 4 ? 4 - (3 * wr) % 4 : 0; pad > 0; pad--) /* ADD */
417 } else { /* Gray-scale */
419 /* -->> -->> -->> -->>
421 8 bits non code (Gray scale)
423 <<-- <<-- <<-- <<-- */
424 f = fopen(argv[2], "wb");
425 // w = ceildiv(img->x1 - img->x0, img->comps[0].dx);
426 // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor), img->comps[0].dx);
428 wr = int_ceildivpow2(img->comps[0].w, img->comps[0].factor);
430 // h = ceildiv(img->y1 - img->y0, img->comps[0].dy);
431 // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), img->comps[0].dy);
433 hr = int_ceildivpow2(img->comps[0].h, img->comps[0].factor);
439 fprintf(f, "%c%c%c%c",
440 (unsigned char) (hr * wr + 54 + 1024 +
441 hr * (wr % 2)) & 0xff,
442 (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2))
444 (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2))
446 (unsigned char) ((hr * wr + 54 + 1024 + wr * (wr % 2))
448 fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
449 ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
450 fprintf(f, "%c%c%c%c", (54 + 1024) & 0xff,
451 ((54 + 1024) >> 8) & 0xff, ((54 + 1024) >> 16) & 0xff,
452 ((54 + 1024) >> 24) & 0xff);
456 fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,
457 ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
458 fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
459 (unsigned char) ((wr) >> 8) & 0xff,
460 (unsigned char) ((wr) >> 16) & 0xff,
461 (unsigned char) ((wr) >> 24) & 0xff);
462 fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
463 (unsigned char) ((hr) >> 8) & 0xff,
464 (unsigned char) ((hr) >> 16) & 0xff,
465 (unsigned char) ((hr) >> 24) & 0xff);
466 fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
467 fprintf(f, "%c%c", (8) & 0xff, ((8) >> 8) & 0xff);
468 fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
469 ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
470 fprintf(f, "%c%c%c%c",
471 (unsigned char) (hr * wr + hr * (wr % 2)) & 0xff,
472 (unsigned char) ((hr * wr + hr * (wr % 2)) >> 8) &
474 (unsigned char) ((hr * wr + hr * (wr % 2)) >> 16) &
476 (unsigned char) ((hr * wr + hr * (wr % 2)) >> 24) & 0xff);
477 fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
478 ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
479 fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
480 ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
481 fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
482 ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
483 fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
484 ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
487 for (i = 0; i < 256; i++) {
488 fprintf(f, "%c%c%c%c", i, i, i, 0);
491 for (i = 0; i < wr * hr; i++) {
493 // fprintf(f, "%c", img->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]);
495 img->comps[0].data[w * hr - ((i) / (wr) + 1) * w +
497 /*if (((i + 1) % w == 0 && w % 2))
498 fprintf(f, "%c", 0); */
499 if ((i + 1) % wr == 0) {
500 for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--) /* ADD */