X-Git-Url: https://main.carlh.net/gitweb/?a=blobdiff_plain;f=codec%2Fj2k_to_image.c;h=4541badd46820d8acae2fa4cad6b6fa04ada9b0e;hb=355b88c15d2dcafe4074b725e2165511df02f2dd;hp=15edba6901583cab426a545db40cc5b48870cf85;hpb=64c34e0d5260ef8e4dee4638458431aec0e1e614;p=openjpeg.git diff --git a/codec/j2k_to_image.c b/codec/j2k_to_image.c index 15edba69..4541badd 100644 --- a/codec/j2k_to_image.c +++ b/codec/j2k_to_image.c @@ -1,8 +1,11 @@ -/* 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. +/* + * Copyright (c) 2001-2003, David Janssens + * Copyright (c) 2002-2003, Yannick Verschueren + * Copyright (c) 2003-2005, Francois Devaux and Antonin Descampe + * Copyright (c) 2005, HervĂ© Drolon, FreeImage Team + * Copyright (c) 2002-2005, 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: @@ -24,487 +27,365 @@ * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE * POSSIBILITY OF SUCH DAMAGE. */ +#include +#include +#include +#include "openjpeg.h" +#include "compat/getopt.h" +#include "convert.h" + +#ifndef WIN32 +#define stricmp strcasecmp +#define strnicmp strncasecmp +#endif + +/* ----------------------------------------------------------------------- */ + +#define J2K_CFMT 0 +#define JP2_CFMT 1 +#define JPT_CFMT 2 +#define MJ2_CFMT 3 +#define PXM_DFMT 0 +#define PGX_DFMT 1 +#define BMP_DFMT 2 +#define YUV_DFMT 3 + +/* ----------------------------------------------------------------------- */ + +void decode_help_display() { + fprintf(stdout,"HELP\n----\n\n"); + fprintf(stdout,"- the -h option displays this help information on screen\n\n"); + + fprintf(stdout,"List of parameters for the JPEG 2000 encoder:\n"); + fprintf(stdout,"\n"); + fprintf(stdout," -i \n"); + fprintf(stdout," REQUIRED\n"); + fprintf(stdout," Currently accepts J2K-files, JP2-files and JPT-files. The file type\n"); + fprintf(stdout," is identified based on its suffix.\n"); + fprintf(stdout," -o \n"); + fprintf(stdout," REQUIRED\n"); + fprintf(stdout," Currently accepts PGM-files, PPM-files, PNM-files, PGX-files and\n"); + fprintf(stdout," BMP-files. Binary data is written to the file (not ascii). If a PGX\n"); + fprintf(stdout," filename is given, there will be as many output files as there are\n"); + fprintf(stdout," components: an indice starting from 0 will then be appended to the\n"); + fprintf(stdout," output filename, just before the \"pgx\" extension. If a PGM filename\n"); + fprintf(stdout," is given and there are more than one component, only the first component\n"); + fprintf(stdout," will be written to the file.\n"); + fprintf(stdout," -r \n"); + fprintf(stdout," Set the number of highest resolution levels to be discarded. The\n"); + fprintf(stdout," image resolution is effectively divided by 2 to the power of the\n"); + fprintf(stdout," number of discarded levels. The reduce factor is limited by the\n"); + fprintf(stdout," smallest total number of decomposition levels among tiles.\n"); + fprintf(stdout," -l \n"); + fprintf(stdout," Set the maximum number of quality layers to decode. If there are\n"); + fprintf(stdout," less quality layers than the specified number, all the quality layers\n"); + fprintf(stdout," are decoded.\n"); + fprintf(stdout,"\n"); +} +/* -------------------------------------------------------------------------- */ + +int get_file_format(char *filename) { + unsigned int i; + static const char *extension[] = {"pgx", "pnm", "pgm", "ppm", "bmp", "j2k", "jp2", "jpt" }; + static const int format[] = { PGX_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, BMP_DFMT, J2K_CFMT, JP2_CFMT, JPT_CFMT }; + char * ext = strrchr(filename, '.') + 1; + if(ext) { + for(i = 0; i < sizeof(format)/sizeof(*format); i++) { + if(strnicmp(ext, extension[i], 3) == 0) { + return format[i]; + } + } + } -#include -#include -#include -#include -//#include + return -1; +} + +/* -------------------------------------------------------------------------- */ + +int parse_cmdline_decoder(int argc, char **argv, opj_dparameters_t *parameters) { + /* parse the command line */ + + while (1) { + int c = getopt(argc, argv, "i:o:r:q:f:t:n:c:b:x:p:s:d:hP:S:E:M:R:T:C:I"); + if (c == -1) + break; + switch (c) { + case 'i': /* input file */ + { + char *infile = optarg; + parameters->decod_format = get_file_format(infile); + switch(parameters->decod_format) { + case J2K_CFMT: + case JP2_CFMT: + case JPT_CFMT: + break; + default: + fprintf(stderr, + "!! Unrecognized format for infile : %s [accept only *.j2k, *.jp2, *.jpc or *.jpt] !!\n\n", + infile); + return 1; + } + strncpy(parameters->infile, infile, MAX_PATH); + } + break; + + /* ----------------------------------------------------- */ + + case 'o': /* output file */ + { + char *outfile = optarg; + parameters->cod_format = get_file_format(outfile); + switch(parameters->cod_format) { + case PGX_DFMT: + case PXM_DFMT: + case BMP_DFMT: + break; + default: + fprintf(stderr, "Unknown output format image %s [only *.pnm, *.pgm, *.ppm, *.pgx or *.bmp]!! \n", outfile); + return 1; + } + strncpy(parameters->outfile, outfile, MAX_PATH); + } + break; + + /* ----------------------------------------------------- */ + + + case 'r': /* reduce option */ + { + sscanf(optarg, "%d", ¶meters->cp_reduce); + } + break; + + /* ----------------------------------------------------- */ + + + case 'l': /* layering option */ + { + sscanf(optarg, "%d", ¶meters->cp_layer); + } + break; + + /* ----------------------------------------------------- */ + + case 'h': /* display an help description */ + decode_help_display(); + return 1; + + /* ----------------------------------------------------- */ + + default: + fprintf(stderr,"WARNING -> this option is not valid \"-%c %s\"\n",c, optarg); + break; + } + } + + /* check for possible errors */ + + if((parameters->infile[0] == 0) || (parameters->outfile[0] == 0)) { + fprintf(stderr,"ERROR -> At least one required argument is missing\nCheck j2k_to_image -h for usage information\n"); + return 1; + } + + return 0; +} + +/* -------------------------------------------------------------------------- */ -int ceildiv(int a, int b) -{ - return (a + b - 1) / b; +/** +sample error callback expecting a FILE* client object +*/ +void error_callback(const char *msg, void *client_data) { + FILE *stream = (FILE*)client_data; + fprintf(stream, "[ERROR] %s", msg); } +/** +sample warning callback expecting a FILE* client object +*/ +void warning_callback(const char *msg, void *client_data) { + FILE *stream = (FILE*)client_data; + fprintf(stream, "[WARNING] %s", msg); +} +/** +sample debug callback expecting no client object +*/ +void info_callback(const char *msg, void *client_data) { + (void)client_data; + fprintf(stdout, "[INFO] %s", msg); +} + +/* -------------------------------------------------------------------------- */ -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; - j2k_option_t option; - int w, wr, wrr, h, hr, hrr, max; - int i, image_type = -1, compno, pad; - int adjust; - - if (argc < 3) { - fprintf(stderr, - "usage: %s j2k-file image-file -reduce n (<- optional)\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]; - - option.reduce_on = 0; - option.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; - } - option.reduce_on = 1; - sscanf(argv[4], "%d", &option.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; - } - - 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; - - if (S1 == 'j' && S2 == '2' && S3 == 'k') { - if (!j2k_decode(src, len, &img, &cp, option)) { - fprintf(stderr, "j2k_to_image: failed to decode image!\n"); - return 1; - } - } else { - if (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; - } - } else { - fprintf(stderr, - "j2k_to_image : Unknown format image *.%c%c%c [only *.j2k or *.jpt]!! \n", - S1, S2, S3); - return 1; - } - } - - 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); - } - 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]); +int main(int argc, char **argv) { + opj_dparameters_t parameters; /* decompression parameters */ + opj_event_mgr_t event_mgr; /* event manager */ + opj_image_t *image = NULL; + FILE *fsrc = NULL; + unsigned char *src = NULL; + int file_length; + + opj_dinfo_t* dinfo = NULL; /* handle to a decompressor */ + opj_cio_t *cio = NULL; + + /* configure the event callbacks (not required) */ + memset(&event_mgr, 0, sizeof(opj_event_mgr_t)); + event_mgr.error_handler = error_callback; + event_mgr.warning_handler = warning_callback; + event_mgr.info_handler = info_callback; + + /* set decoding parameters to default values */ + opj_set_default_decoder_parameters(¶meters); + + /* parse input and get user decoding parameters */ + if(parse_cmdline_decoder(argc, argv, ¶meters) == 1) { + return 0; } - 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); + + /* read the input file and put it in memory */ + /* ---------------------------------------- */ + fsrc = fopen(parameters.infile, "rb"); + if (!fsrc) { + fprintf(stderr, "ERROR -> failed to open %s for reading\n", parameters.infile); + return 1; + } + fseek(fsrc, 0, SEEK_END); + file_length = ftell(fsrc); + fseek(fsrc, 0, SEEK_SET); + src = (unsigned char *) malloc(file_length); + fread(src, 1, file_length, fsrc); + fclose(fsrc); + + /* decode the code-stream */ + /* ---------------------- */ + + switch(parameters.decod_format) { + case J2K_CFMT: + { + /* JPEG-2000 codestream */ + + /* get a decoder handle */ + dinfo = opj_create_decompress(CODEC_J2K); + + /* catch events using our callbacks and give a local context */ + opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr); + + /* setup the decoder decoding parameters using user parameters */ + opj_setup_decoder(dinfo, ¶meters); + + /* open a byte stream */ + cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length); + + /* decode the stream and fill the image structure */ + image = opj_decode(dinfo, cio); + if(!image) { + fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n"); + opj_destroy_decompress(dinfo); + opj_cio_close(cio); + return 1; + } + + /* close the byte stream */ + opj_cio_close(cio); + } + break; + + case JP2_CFMT: + { + /* JPEG 2000 compressed image data */ + + /* get a decoder handle */ + dinfo = opj_create_decompress(CODEC_JP2); + + /* catch events using our callbacks and give a local context */ + opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr); + + /* setup the decoder decoding parameters using the current image and using user parameters */ + opj_setup_decoder(dinfo, ¶meters); + + /* open a byte stream */ + cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length); + + /* decode the stream and fill the image structure */ + image = opj_decode(dinfo, cio); + if(!image) { + fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n"); + opj_destroy_decompress(dinfo); + opj_cio_close(cio); + return 1; + } + + /* close the byte stream */ + opj_cio_close(cio); + + } + break; + + case JPT_CFMT: + { + /* JPEG 2000, JPIP */ + + /* get a decoder handle */ + dinfo = opj_create_decompress(CODEC_JPT); + + /* catch events using our callbacks and give a local context */ + opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr); + + /* setup the decoder decoding parameters using user parameters */ + opj_setup_decoder(dinfo, ¶meters); + + /* open a byte stream */ + cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length); + + /* decode the stream and fill the image structure */ + image = opj_decode(dinfo, cio); + if(!image) { + fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n"); + opj_destroy_decompress(dinfo); + opj_cio_close(cio); + return 1; + } + + /* close the byte stream */ + opj_cio_close(cio); + } + break; + + default: + fprintf(stderr, "ERROR -> j2k_to_image : Unknown input image format\n"); + return 1; } - fclose(f); - } - } - break; - - /* ------------------------ / */ - /* / / */ - /* / FORMAT : PGX / */ - /* / / */ - /* /----------------------- / */ - case 0: /* PGX */ - for (compno = 0; compno < img->numcomps; compno++) { - j2k_comp_t *comp = &img->comps[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, 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 LM %c %d %d %d\n", comp->sgnd ? '-' : '+', - comp->prec, wr, hr); - for (i = 0; i < wr * hr; i++) { - int v = img->comps[compno].data[i / wr * w + i % wr]; - if (comp->prec <= 8) { - char c = (char) v; - fwrite(&c, 1, 1, f); - } else if (comp->prec <= 16) { - short s = (short) v; - fwrite(&s, 2, 1, f); - } else { - fwrite(&v, 4, 1, f); + + /* free the memory containing the code-stream */ + free(src); + src = NULL; + + /* create output image */ + /* ------------------- */ + + switch (parameters.cod_format) { + case PXM_DFMT: /* PNM PGM PPM */ + imagetopnm(image, parameters.outfile); + break; + + case PGX_DFMT: /* PGX */ + imagetopgx(image, parameters.outfile); + break; + + case BMP_DFMT: /* BMP */ + imagetobmp(image, parameters.outfile); + break; } - } - 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); + + /* free remaining structures */ + if(dinfo) { + opj_destroy_decompress(dinfo); } - } - fclose(f); - } 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); - } - } - break; - default: - break; - } - - return 0; + + /* free image data structure */ + opj_image_destroy(image); + + return 0; } +