diff options
Diffstat (limited to 'jpwl/decoder/jpwldec/j2k_to_image.c')
| -rw-r--r-- | jpwl/decoder/jpwldec/j2k_to_image.c | 561 |
1 files changed, 561 insertions, 0 deletions
diff --git a/jpwl/decoder/jpwldec/j2k_to_image.c b/jpwl/decoder/jpwldec/j2k_to_image.c new file mode 100644 index 00000000..0c61f1cf --- /dev/null +++ b/jpwl/decoder/jpwldec/j2k_to_image.c @@ -0,0 +1,561 @@ +/* Copyright (c) 2001 David Janssens + * Copyright (c) 2002-2003 Yannick Verschueren + * Copyright (c) 2002-2003 Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium + * + * All rights reserved. + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in the + * documentation and/or other materials provided with the distribution. + * + * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS' + * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE + * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE + * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE + * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR + * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF + * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS + * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN + * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) + * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE + * POSSIBILITY OF SUCH DAMAGE. + */ + + + +#include <openjpeg.h> +#include <stdio.h> +#include <stdlib.h> +#include <string.h> + +extern int cslen; // Aggiunta in JPWL (per gestire l'inserimento di RED nella CS!) + +int ceildiv(int a, int b) +{ + return (a + b - 1) / b; +} + +int main(int argc, char **argv) +{ + FILE *f; + char *src, *src_name; + char *dest, S1, S2, S3; + int len; + + j2k_image_t img; + + j2k_cp_t cp; + int w, wr, wrr, h, hr, hrr, max; + int i, image_type = -1, compno, pad, j; + int adjust; + jp2_struct_t *jp2_struct; + + if (argc < 3) { + fprintf(stderr, + "usage: %s j2k-file image-file [-reduce n]\n", argv[0]); + return 1; + } + + f = fopen(argv[1], "rb"); + if (!f) { + fprintf(stderr, "failed to open %s for reading\n", argv[1]); + return 1; + } + + dest = argv[2]; + + cp.reduce_on = 0; + cp.reduce_value = 0; + + /* OPTION REDUCE IS ACTIVE */ + if (argc == 5) { + if (strcmp(argv[3], "-reduce")) { + fprintf(stderr, + "usage: options " "-reduce n" + " where n is the factor of reduction [%s]\n", argv[3]); + return 1; + } + cp.reduce_on = 1; + sscanf(argv[4], "%d", &cp.reduce_value); + } + + while (*dest) { + dest++; + } + dest--; + S3 = *dest; + dest--; + S2 = *dest; + dest--; + S1 = *dest; + + if ((S1 == 'p' && S2 == 'g' && S3 == 'x') + || (S1 == 'P' && S2 == 'G' && S3 == 'X')) { + image_type = 0; + + dest--; + + *dest = '\0'; + } + + if ((S1 == 'p' && S2 == 'n' && S3 == 'm') + || (S1 == 'P' && S2 == 'N' && S3 == 'M') || (S1 == 'p' && S2 == 'g' + && S3 == 'm') + || (S1 == 'P' && S2 == 'G' && S3 == 'M') || (S1 == 'P' && S2 == 'P' + && S3 == 'M') + || (S1 == 'p' && S2 == 'p' && S3 == 'm')) { + image_type = 1; + } + + if ((S1 == 'b' && S2 == 'm' && S3 == 'p') + || (S1 == 'B' && S2 == 'M' && S3 == 'P')) { + image_type = 2; + } + + if (image_type == -1) { + fprintf(stderr, + "!! Unrecognized format for infile : %c%c%c [accept only *.pnm, *.pgm, *.ppm, *.pgx or *.bmp] !!\n\n", + S1, S2, S3); + return 1; + } + + fseek(f, 0, SEEK_END); + len = ftell(f); + fseek(f, 0, SEEK_SET); + src = (char *) malloc(len); + fread(src, 1, len, f); + fclose(f); + + src_name = argv[1]; + while (*src_name) { + src_name++; + } + src_name--; + S3 = *src_name; + src_name--; + S2 = *src_name; + src_name--; + S1 = *src_name; + + /* J2K format */ + if ((S1 == 'j' && S2 == '2' && S3 == 'k') + || (S1 == 'J' && S2 == '2' && S3 == 'K') || (S1 == 'j' && S2 == '2' + && S3 == 'c') + || (S1 == 'J' && S2 == '2' && S3 == 'C')) { + if (!j2k_decode(src, len, &img, &cp)) { + fprintf(stderr, "j2k_to_image: failed to decode image!\n"); + return 1; + } + } + + /* JP2 format */ + else if ((S1 == 'j' && S2 == 'p' && S3 == '2') + || (S1 == 'J' && S2 == 'P' && S3 == '2')) { + jp2_struct = (jp2_struct_t *) malloc(sizeof(jp2_struct_t)); + + jp2_struct->image = &img; + + if (jp2_decode(src, len, jp2_struct, &cp)) { + fprintf(stderr, "j2k_to_image: failed to decode image!\n"); + return 1; + } + /* Insert code here if you want to create actions on jp2_struct before deleting it */ + free(jp2_struct); + } + + /* JPT format */ + else if ((S1 == 'j' && S2 == 'p' && S3 == 't') + || (S1 == 'J' && S2 == 'P' && S3 == 'T')) { + if (!j2k_decode_jpt_stream(src, len, &img, &cp)) { + fprintf(stderr, "j2k_to_image: failed to decode image!\n"); + return 1; + } + } + + /* otherwise : error */ + else { + fprintf(stderr, + "j2k_to_image : Unknown format image *.%c%c%c [only *.j2k, *.jp2, *.jpc or *.jpt]!! \n", + S1, S2, S3); + return 1; + } + + src = realloc(src, cslen); // Aggiunta in JPWL per evitare errori di free a causa + // della modifica di CS in caso di aggiunta RED! + free(src); + /* ------------------ CREATE OUT IMAGE WITH THE RIGHT FORMAT ----------------------- */ + + /* ---------------------------- / */ + /* / / */ + /* / FORMAT : PNM, PGM or PPM / */ + /* / / */ + /* ---------------------------- / */ + + switch (image_type) { + case 1: /* PNM PGM PPM */ + if (img.numcomps == 3 && img.comps[0].dx == img.comps[1].dx + && img.comps[1].dx == img.comps[2].dx + && img.comps[0].dy == img.comps[1].dy + && img.comps[1].dy == img.comps[2].dy + && img.comps[0].prec == img.comps[1].prec + && img.comps[1].prec == img.comps[2].prec) { + f = fopen(argv[2], "wb"); + w = ceildiv(img.x1 - img.x0, img.comps[0].dx); + // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor),img.comps[0].dx); + wr = img.comps[0].w; + wrr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor); + + h = ceildiv(img.y1 - img.y0, img.comps[0].dy); + // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy); + hr = img.comps[0].h; + hrr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor); + + max = img.comps[0].prec > 8 ? 255 : (1 << img.comps[0].prec) - 1; + + img.comps[0].x0 = + int_ceildivpow2(img.comps[0].x0 - + int_ceildiv(img.x0, img.comps[0].dx), + img.comps[0].factor); + img.comps[0].y0 = + int_ceildivpow2(img.comps[0].y0 - + int_ceildiv(img.y0, img.comps[0].dy), + img.comps[0].factor); + + + fprintf(f, "P6\n# %d %d %d %d %d\n%d %d\n%d\n", + cp.tcps[cp.tileno[0]].tccps[0].numresolutions, w, h, + img.comps[0].x0, img.comps[0].y0, wrr, hrr, max); + adjust = img.comps[0].prec > 8 ? img.comps[0].prec - 8 : 0; + for (i = 0; i < wrr * hrr; i++) { + char r, g, b; + r = img.comps[0].data[i / wrr * wr + i % wrr]; + r += (img.comps[0].sgnd ? 1 << (img.comps[0].prec - 1) : 0); + r = r >> adjust; + + g = img.comps[1].data[i / wrr * wr + i % wrr]; + g += (img.comps[1].sgnd ? 1 << (img.comps[1].prec - 1) : 0); + g = g >> adjust; + + b = img.comps[2].data[i / wrr * wr + i % wrr]; + b += (img.comps[2].sgnd ? 1 << (img.comps[2].prec - 1) : 0); + b = b >> adjust; + + fprintf(f, "%c%c%c", r, g, b); + } + free(img.comps[0].data); + free(img.comps[1].data); + free(img.comps[2].data); + fclose(f); + } else { + for (compno = 0; compno < img.numcomps; compno++) { + char name[256]; + if (img.numcomps > 1) { + sprintf(name, "%d.%s", compno, argv[2]); + } else { + sprintf(name, "%s", argv[2]); + } + f = fopen(name, "wb"); + w = ceildiv(img.x1 - img.x0, img.comps[compno].dx); + // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor),img.comps[compno].dx); + wr = img.comps[compno].w; + wrr = + int_ceildivpow2(img.comps[compno].w, img.comps[compno].factor); + + h = ceildiv(img.y1 - img.y0, img.comps[compno].dy); + // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[compno].dy); + hr = img.comps[compno].h; + hrr = + int_ceildivpow2(img.comps[compno].h, img.comps[compno].factor); + + max = + img.comps[compno].prec > + 8 ? 255 : (1 << img.comps[compno].prec) - 1; + + img.comps[compno].x0 = + int_ceildivpow2(img.comps[compno].x0 - + int_ceildiv(img.x0, + img.comps[compno].dx), + img.comps[compno].factor); + img.comps[compno].y0 = + int_ceildivpow2(img.comps[compno].y0 - + int_ceildiv(img.y0, + img.comps[compno].dy), + img.comps[compno].factor); + + fprintf(f, "P5\n# %d %d %d %d %d\n%d %d\n%d\n", + cp.tcps[cp.tileno[0]].tccps[compno]. + numresolutions, w, h, img.comps[compno].x0, + img.comps[compno].y0, wrr, hrr, max); + adjust = + img.comps[compno].prec > 8 ? img.comps[compno].prec - 8 : 0; + for (i = 0; i < wrr * hrr; i++) { + char l; + l = img.comps[compno].data[i / wrr * wr + i % wrr]; + l += (img.comps[compno]. + sgnd ? 1 << (img.comps[compno].prec - 1) : 0); + l = l >> adjust; + fprintf(f, "%c", l); + } + fclose(f); + free(img.comps[compno].data); + } + } + break; + + /* ------------------------ / */ + /* / / */ + /* / FORMAT : PGX / */ + /* / / */ + /* /----------------------- / */ + case 0: /* PGX */ + for (compno = 0; compno < img.numcomps; compno++) { + j2k_comp_t *comp = &img.comps[compno]; + char name[256]; + + int nbytes = 0; + //if (img.numcomps > 1) + sprintf(name, "%s-%d.pgx", argv[2], compno); + + //else + + //sprintf(name, "%s.pgx", argv[2]); + + f = fopen(name, "wb"); + // w = ceildiv(img.x1 - img.x0, comp->dx); + // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), comp->dx); + w = img.comps[compno].w; + wr = int_ceildivpow2(img.comps[compno].w, img.comps[compno].factor); + + // h = ceildiv(img.y1 - img.y0, comp->dy); + // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), comp->dy); + h = img.comps[compno].h; + hr = int_ceildivpow2(img.comps[compno].h, img.comps[compno].factor); + + fprintf(f, "PG ML %c %d %d %d\n", comp->sgnd ? '-' : '+', + comp->prec, wr, hr); + + if (comp->prec <= 8) + nbytes = 1; + + else if (comp->prec <= 16) + nbytes = 2; + + else + nbytes = 4; + for (i = 0; i < wr * hr; i++) { + int v = img.comps[compno].data[i / wr * w + i % wr]; + + for (j = nbytes - 1; j >= 0; j--) { + + char byte = (char) (v >> (j * 8)); + + fwrite(&byte, 1, 1, f); + + } + } + free(img.comps[compno].data); + fclose(f); + } + break; + + /* ------------------------ / */ + /* / / */ + /* / FORMAT : BMP / */ + /* / / */ + /* /----------------------- / */ + + case 2: /* BMP */ + if (img.numcomps == 3 && img.comps[0].dx == img.comps[1].dx + && img.comps[1].dx == img.comps[2].dx + && img.comps[0].dy == img.comps[1].dy + && img.comps[1].dy == img.comps[2].dy + && img.comps[0].prec == img.comps[1].prec + && img.comps[1].prec == img.comps[2].prec) { + /* -->> -->> -->> -->> + + 24 bits color + + <<-- <<-- <<-- <<-- */ + + f = fopen(argv[2], "wb"); + // w = ceildiv(img.x1 - img.x0, img.comps[0].dx); + // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), img.comps[0].dx); + w = img.comps[0].w; + wr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor); + + // h = ceildiv(img.y1 - img.y0, img.comps[0].dy); + // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy); + h = img.comps[0].h; + hr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor); + + fprintf(f, "BM"); + + /* FILE HEADER */ + /* ------------- */ + fprintf(f, "%c%c%c%c", + (unsigned char) (hr * wr * 3 + 3 * hr * (wr % 2) + + 54) & 0xff, + (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) + >> 8) & 0xff, + (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) + >> 16) & 0xff, + (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) + >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, + ((0) >> 16) & 0xff, ((0) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff, + ((54) >> 16) & 0xff, ((54) >> 24) & 0xff); + + /* INFO HEADER */ + /* ------------- */ + fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff, + ((40) >> 16) & 0xff, ((40) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff), + (unsigned char) ((wr) >> 8) & 0xff, + (unsigned char) ((wr) >> 16) & 0xff, + (unsigned char) ((wr) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff), + (unsigned char) ((hr) >> 8) & 0xff, + (unsigned char) ((hr) >> 16) & 0xff, + (unsigned char) ((hr) >> 24) & 0xff); + fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff); + fprintf(f, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff); + fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, + ((0) >> 16) & 0xff, ((0) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", + (unsigned char) (3 * hr * wr + + 3 * hr * (wr % 2)) & 0xff, + (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> + 8) & 0xff, + (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> + 16) & 0xff, + (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> + 24) & 0xff); + fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, + ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, + ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, + ((0) >> 16) & 0xff, ((0) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, + ((0) >> 16) & 0xff, ((0) >> 24) & 0xff); + + for (i = 0; i < wr * hr; i++) { + unsigned char R, G, B; + /* a modifier */ + // R = img.comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; + R = img.comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)]; + // G = img.comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; + G = img.comps[1].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)]; + // B = img.comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; + B = img.comps[2].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)]; + fprintf(f, "%c%c%c", B, G, R); + + if ((i + 1) % wr == 0) { + for (pad = (3 * wr) % 4 ? 4 - (3 * wr) % 4 : 0; pad > 0; pad--) /* ADD */ + fprintf(f, "%c", 0); + } + } + fclose(f); + free(img.comps[1].data); + free(img.comps[2].data); + } else { /* Gray-scale */ + + /* -->> -->> -->> -->> + + 8 bits non code (Gray scale) + + <<-- <<-- <<-- <<-- */ + f = fopen(argv[2], "wb"); + // w = ceildiv(img.x1 - img.x0, img.comps[0].dx); + // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), img.comps[0].dx); + w = img.comps[0].w; + wr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor); + + // h = ceildiv(img.y1 - img.y0, img.comps[0].dy); + // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy); + h = img.comps[0].h; + hr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor); + + fprintf(f, "BM"); + + /* FILE HEADER */ + /* ------------- */ + fprintf(f, "%c%c%c%c", + (unsigned char) (hr * wr + 54 + 1024 + + hr * (wr % 2)) & 0xff, + (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2)) + >> 8) & 0xff, + (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2)) + >> 16) & 0xff, + (unsigned char) ((hr * wr + 54 + 1024 + wr * (wr % 2)) + >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, + ((0) >> 16) & 0xff, ((0) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (54 + 1024) & 0xff, + ((54 + 1024) >> 8) & 0xff, ((54 + 1024) >> 16) & 0xff, + ((54 + 1024) >> 24) & 0xff); + + /* INFO HEADER */ + /* ------------- */ + fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff, + ((40) >> 16) & 0xff, ((40) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff), + (unsigned char) ((wr) >> 8) & 0xff, + (unsigned char) ((wr) >> 16) & 0xff, + (unsigned char) ((wr) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff), + (unsigned char) ((hr) >> 8) & 0xff, + (unsigned char) ((hr) >> 16) & 0xff, + (unsigned char) ((hr) >> 24) & 0xff); + fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff); + fprintf(f, "%c%c", (8) & 0xff, ((8) >> 8) & 0xff); + fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, + ((0) >> 16) & 0xff, ((0) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", + (unsigned char) (hr * wr + hr * (wr % 2)) & 0xff, + (unsigned char) ((hr * wr + hr * (wr % 2)) >> 8) & + 0xff, + (unsigned char) ((hr * wr + hr * (wr % 2)) >> 16) & + 0xff, + (unsigned char) ((hr * wr + hr * (wr % 2)) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, + ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, + ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff, + ((256) >> 16) & 0xff, ((256) >> 24) & 0xff); + fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff, + ((256) >> 16) & 0xff, ((256) >> 24) & 0xff); + } + + for (i = 0; i < 256; i++) { + fprintf(f, "%c%c%c%c", i, i, i, 0); + } + + for (i = 0; i < wr * hr; i++) { + /* a modifier !! */ + // fprintf(f, "%c", img.comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]); + fprintf(f, "%c", + img.comps[0].data[w * hr - ((i) / (wr) + 1) * w + + (i) % (wr)]); + /*if (((i + 1) % w == 0 && w % 2)) + fprintf(f, "%c", 0); */ + if ((i + 1) % wr == 0) { + for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--) /* ADD */ + fprintf(f, "%c", 0); + } + } + fclose(f); + free(img.comps[0].data); + break; + default: + break; + } + + return 0; +} |
