-/*\r
- * Copyright (c) 2001-2003, David Janssens\r
- * Copyright (c) 2002-2003, Yannick Verschueren\r
- * Copyright (c) 2003-2005, Francois Devaux and Antonin Descampe\r
- * Copyright (c) 2005, Hervé Drolon, FreeImage Team\r
- * Copyright (c) 2002-2005, Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium\r
- * All rights reserved.\r
- *\r
- * Redistribution and use in source and binary forms, with or without\r
- * modification, are permitted provided that the following conditions\r
- * are met:\r
- * 1. Redistributions of source code must retain the above copyright\r
- * notice, this list of conditions and the following disclaimer.\r
- * 2. Redistributions in binary form must reproduce the above copyright\r
- * notice, this list of conditions and the following disclaimer in the\r
- * documentation and/or other materials provided with the distribution.\r
- *\r
- * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'\r
- * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE\r
- * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE\r
- * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE\r
- * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR\r
- * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF\r
- * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS\r
- * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN\r
- * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)\r
- * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE\r
- * POSSIBILITY OF SUCH DAMAGE.\r
- */\r
-#include <stdio.h>\r
-#include <stdlib.h>\r
-#include <string.h>\r
-#include "openjpeg.h"\r
-\r
-/*\r
- * Get logarithm of an integer and round downwards.\r
- *\r
- * log2(a)\r
- */\r
-static int int_floorlog2(int a) {\r
- int l;\r
- for (l = 0; a > 1; l++) {\r
- a >>= 1;\r
- }\r
- return l;\r
-}\r
-\r
-/*\r
- * Divide an integer by a power of 2 and round upwards.\r
- *\r
- * a divided by 2^b\r
- */\r
-static int int_ceildivpow2(int a, int b) {\r
- return (a + (1 << b) - 1) >> b;\r
-}\r
-\r
-/*\r
- * Divide an integer and round upwards.\r
- *\r
- * a divided by b\r
- */\r
-static int int_ceildiv(int a, int b) {\r
- return (a + b - 1) / b;\r
-}\r
-\r
-/* -->> -->> -->> -->>\r
-\r
- BMP IMAGE FORMAT\r
-\r
- <<-- <<-- <<-- <<-- */\r
-\r
-/* WORD defines a two byte word */\r
-typedef unsigned short int WORD;\r
-\r
-/* DWORD defines a four byte word */\r
-typedef unsigned long int DWORD;\r
-\r
-typedef struct {\r
- WORD bfType; /* 'BM' for Bitmap (19776) */\r
- DWORD bfSize; /* Size of the file */\r
- WORD bfReserved1; /* Reserved : 0 */\r
- WORD bfReserved2; /* Reserved : 0 */\r
- DWORD bfOffBits; /* Offset */\r
-} BITMAPFILEHEADER_t;\r
-\r
-typedef struct {\r
- DWORD biSize; /* Size of the structure in bytes */\r
- DWORD biWidth; /* Width of the image in pixels */\r
- DWORD biHeight; /* Heigth of the image in pixels */\r
- WORD biPlanes; /* 1 */\r
- WORD biBitCount; /* Number of color bits by pixels */\r
- DWORD biCompression; /* Type of encoding 0: none 1: RLE8 2: RLE4 */\r
- DWORD biSizeImage; /* Size of the image in bytes */\r
- DWORD biXpelsPerMeter; /* Horizontal (X) resolution in pixels/meter */\r
- DWORD biYpelsPerMeter; /* Vertical (Y) resolution in pixels/meter */\r
- DWORD biClrUsed; /* Number of color used in the image (0: ALL) */\r
- DWORD biClrImportant; /* Number of important color (0: ALL) */\r
-} BITMAPINFOHEADER_t;\r
-\r
-opj_image_t* bmptoimage(char *filename, opj_cparameters_t *parameters) {\r
- int subsampling_dx = parameters->subsampling_dx;\r
- int subsampling_dy = parameters->subsampling_dy;\r
-\r
- int i, numcomps, w, h;\r
- OPJ_COLOR_SPACE color_space;\r
- opj_image_cmptparm_t cmptparm[3]; /* maximum of 3 components */\r
- opj_image_t * image = NULL;\r
-\r
- FILE *IN;\r
- BITMAPFILEHEADER_t File_h;\r
- BITMAPINFOHEADER_t Info_h;\r
- unsigned char *RGB;\r
- unsigned char *table_R, *table_G, *table_B;\r
- unsigned int j, PAD = 0;\r
-\r
- int x, y, index;\r
- int gray_scale = 1, not_end_file = 1; \r
-\r
- unsigned int line = 0, col = 0;\r
- unsigned char v, v2;\r
- DWORD W, H;\r
- \r
- IN = fopen(filename, "rb");\r
- if (!IN) {\r
- fprintf(stderr, "\033[0;33mFailed to open %s for reading !!\033[0;39m\n", filename);\r
- return 0;\r
- }\r
- \r
- File_h.bfType = getc(IN);\r
- File_h.bfType = (getc(IN) << 8) + File_h.bfType;\r
- \r
- if (File_h.bfType != 19778) {\r
- fprintf(stderr,"Error, not a BMP file!\n");\r
- return 0;\r
- } else {\r
- /* FILE HEADER */\r
- /* ------------- */\r
- File_h.bfSize = getc(IN);\r
- File_h.bfSize = (getc(IN) << 8) + File_h.bfSize;\r
- File_h.bfSize = (getc(IN) << 16) + File_h.bfSize;\r
- File_h.bfSize = (getc(IN) << 24) + File_h.bfSize;\r
-\r
- File_h.bfReserved1 = getc(IN);\r
- File_h.bfReserved1 = (getc(IN) << 8) + File_h.bfReserved1;\r
-\r
- File_h.bfReserved2 = getc(IN);\r
- File_h.bfReserved2 = (getc(IN) << 8) + File_h.bfReserved2;\r
-\r
- File_h.bfOffBits = getc(IN);\r
- File_h.bfOffBits = (getc(IN) << 8) + File_h.bfOffBits;\r
- File_h.bfOffBits = (getc(IN) << 16) + File_h.bfOffBits;\r
- File_h.bfOffBits = (getc(IN) << 24) + File_h.bfOffBits;\r
-\r
- /* INFO HEADER */\r
- /* ------------- */\r
-\r
- Info_h.biSize = getc(IN);\r
- Info_h.biSize = (getc(IN) << 8) + Info_h.biSize;\r
- Info_h.biSize = (getc(IN) << 16) + Info_h.biSize;\r
- Info_h.biSize = (getc(IN) << 24) + Info_h.biSize;\r
-\r
- Info_h.biWidth = getc(IN);\r
- Info_h.biWidth = (getc(IN) << 8) + Info_h.biWidth;\r
- Info_h.biWidth = (getc(IN) << 16) + Info_h.biWidth;\r
- Info_h.biWidth = (getc(IN) << 24) + Info_h.biWidth;\r
- w = Info_h.biWidth;\r
-\r
- Info_h.biHeight = getc(IN);\r
- Info_h.biHeight = (getc(IN) << 8) + Info_h.biHeight;\r
- Info_h.biHeight = (getc(IN) << 16) + Info_h.biHeight;\r
- Info_h.biHeight = (getc(IN) << 24) + Info_h.biHeight;\r
- h = Info_h.biHeight;\r
-\r
- Info_h.biPlanes = getc(IN);\r
- Info_h.biPlanes = (getc(IN) << 8) + Info_h.biPlanes;\r
-\r
- Info_h.biBitCount = getc(IN);\r
- Info_h.biBitCount = (getc(IN) << 8) + Info_h.biBitCount;\r
-\r
- Info_h.biCompression = getc(IN);\r
- Info_h.biCompression = (getc(IN) << 8) + Info_h.biCompression;\r
- Info_h.biCompression = (getc(IN) << 16) + Info_h.biCompression;\r
- Info_h.biCompression = (getc(IN) << 24) + Info_h.biCompression;\r
-\r
- Info_h.biSizeImage = getc(IN);\r
- Info_h.biSizeImage = (getc(IN) << 8) + Info_h.biSizeImage;\r
- Info_h.biSizeImage = (getc(IN) << 16) + Info_h.biSizeImage;\r
- Info_h.biSizeImage = (getc(IN) << 24) + Info_h.biSizeImage;\r
-\r
- Info_h.biXpelsPerMeter = getc(IN);\r
- Info_h.biXpelsPerMeter = (getc(IN) << 8) + Info_h.biXpelsPerMeter;\r
- Info_h.biXpelsPerMeter = (getc(IN) << 16) + Info_h.biXpelsPerMeter;\r
- Info_h.biXpelsPerMeter = (getc(IN) << 24) + Info_h.biXpelsPerMeter;\r
-\r
- Info_h.biYpelsPerMeter = getc(IN);\r
- Info_h.biYpelsPerMeter = (getc(IN) << 8) + Info_h.biYpelsPerMeter;\r
- Info_h.biYpelsPerMeter = (getc(IN) << 16) + Info_h.biYpelsPerMeter;\r
- Info_h.biYpelsPerMeter = (getc(IN) << 24) + Info_h.biYpelsPerMeter;\r
-\r
- Info_h.biClrUsed = getc(IN);\r
- Info_h.biClrUsed = (getc(IN) << 8) + Info_h.biClrUsed;\r
- Info_h.biClrUsed = (getc(IN) << 16) + Info_h.biClrUsed;\r
- Info_h.biClrUsed = (getc(IN) << 24) + Info_h.biClrUsed;\r
-\r
- Info_h.biClrImportant = getc(IN);\r
- Info_h.biClrImportant = (getc(IN) << 8) + Info_h.biClrImportant;\r
- Info_h.biClrImportant = (getc(IN) << 16) + Info_h.biClrImportant;\r
- Info_h.biClrImportant = (getc(IN) << 24) + Info_h.biClrImportant;\r
-\r
- /* Read the data and store them in the OUT file */\r
- \r
- if (Info_h.biBitCount == 24) {\r
- numcomps = 3;\r
- color_space = CLRSPC_SRGB;\r
- /* initialize image components */\r
- memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));\r
- for(i = 0; i < numcomps; i++) {\r
- cmptparm[i].prec = 8;\r
- cmptparm[i].bpp = 8;\r
- cmptparm[i].sgnd = 0;\r
- cmptparm[i].dx = subsampling_dx;\r
- cmptparm[i].dy = subsampling_dy;\r
- cmptparm[i].w = w;\r
- cmptparm[i].h = h;\r
- }\r
- /* create the image */\r
- image = opj_image_create(numcomps, &cmptparm[0], color_space);\r
- if(!image) {\r
- fclose(IN);\r
- return NULL;\r
- }\r
-\r
- /* set image offset and reference grid */\r
- image->x0 = parameters->image_offset_x0;\r
- image->y0 = parameters->image_offset_y0;\r
- image->x1 = !image->x0 ? (w - 1) * subsampling_dx + 1 : image->x0 + (w - 1) * subsampling_dx + 1;\r
- image->y1 = !image->y0 ? (h - 1) * subsampling_dy + 1 : image->y0 + (h - 1) * subsampling_dy + 1;\r
-\r
- /* set image data */\r
-\r
- /* Place the cursor at the beginning of the image information */\r
- fseek(IN, 0, SEEK_SET);\r
- fseek(IN, File_h.bfOffBits, SEEK_SET);\r
- \r
- W = Info_h.biWidth;\r
- H = Info_h.biHeight;\r
-\r
- /* PAD = 4 - (3 * W) % 4; */\r
- /* PAD = (PAD == 4) ? 0 : PAD; */\r
- PAD = (3 * W) % 4 ? 4 - (3 * W) % 4 : 0;\r
- \r
- RGB = (unsigned char *) malloc((3 * W + PAD) * H * sizeof(unsigned char));\r
- \r
- fread(RGB, sizeof(unsigned char), (3 * W + PAD) * H, IN);\r
- \r
- index = 0;\r
-\r
- for(y = 0; y < (int)H; y++) {\r
- unsigned char *scanline = RGB + (3 * W + PAD) * (H - 1 - y);\r
- for(x = 0; x < (int)W; x++) {\r
- unsigned char *pixel = &scanline[3 * x];\r
- image->comps[0].data[index] = pixel[2]; /* R */\r
- image->comps[1].data[index] = pixel[1]; /* G */\r
- image->comps[2].data[index] = pixel[0]; /* B */\r
- index++;\r
- }\r
- }\r
-\r
- free(RGB);\r
-\r
- } else if (Info_h.biBitCount == 8 && Info_h.biCompression == 0) {\r
- table_R = (unsigned char *) malloc(256 * sizeof(unsigned char));\r
- table_G = (unsigned char *) malloc(256 * sizeof(unsigned char));\r
- table_B = (unsigned char *) malloc(256 * sizeof(unsigned char));\r
- \r
- for (j = 0; j < Info_h.biClrUsed; j++) {\r
- table_B[j] = getc(IN);\r
- table_G[j] = getc(IN);\r
- table_R[j] = getc(IN);\r
- getc(IN);\r
- if (table_R[j] != table_G[j] && table_R[j] != table_B[j] && table_G[j] != table_B[j])\r
- gray_scale = 0;\r
- }\r
- \r
- /* Place the cursor at the beginning of the image information */\r
- fseek(IN, 0, SEEK_SET);\r
- fseek(IN, File_h.bfOffBits, SEEK_SET);\r
- \r
- W = Info_h.biWidth;\r
- H = Info_h.biHeight;\r
- if (Info_h.biWidth % 2)\r
- W++;\r
- \r
- numcomps = gray_scale ? 1 : 3;\r
- color_space = gray_scale ? CLRSPC_GRAY : CLRSPC_SRGB;\r
- /* initialize image components */\r
- memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));\r
- for(i = 0; i < numcomps; i++) {\r
- cmptparm[i].prec = 8;\r
- cmptparm[i].bpp = 8;\r
- cmptparm[i].sgnd = 0;\r
- cmptparm[i].dx = subsampling_dx;\r
- cmptparm[i].dy = subsampling_dy;\r
- cmptparm[i].w = w;\r
- cmptparm[i].h = h;\r
- }\r
- /* create the image */\r
- image = opj_image_create(numcomps, &cmptparm[0], color_space);\r
- if(!image) {\r
- fclose(IN);\r
- return NULL;\r
- }\r
-\r
- /* set image offset and reference grid */\r
- image->x0 = parameters->image_offset_x0;\r
- image->y0 = parameters->image_offset_y0;\r
- image->x1 = !image->x0 ? (w - 1) * subsampling_dx + 1 : image->x0 + (w - 1) * subsampling_dx + 1;\r
- image->y1 = !image->y0 ? (h - 1) * subsampling_dy + 1 : image->y0 + (h - 1) * subsampling_dy + 1;\r
-\r
- /* set image data */\r
-\r
- RGB = (unsigned char *) malloc(W * H * sizeof(unsigned char));\r
- \r
- fread(RGB, sizeof(unsigned char), W * H, IN);\r
- if (gray_scale) {\r
- index = 0;\r
- for (j = 0; j < W * H; j++) {\r
- if ((j % W < W - 1 && Info_h.biWidth % 2) || !(Info_h.biWidth % 2)) {\r
- image->comps[0].data[index] = table_R[RGB[W * H - ((j) / (W) + 1) * W + (j) % (W)]];\r
- index++;\r
- }\r
- }\r
-\r
- } else { \r
- index = 0;\r
- for (j = 0; j < W * H; j++) {\r
- if ((j % W < W - 1 && Info_h.biWidth % 2) || !(Info_h.biWidth % 2)) {\r
- unsigned char pixel_index = RGB[W * H - ((j) / (W) + 1) * W + (j) % (W)];\r
- image->comps[0].data[index] = table_R[pixel_index];\r
- image->comps[1].data[index] = table_G[pixel_index];\r
- image->comps[2].data[index] = table_B[pixel_index];\r
- index++;\r
- }\r
- }\r
- }\r
- free(RGB);\r
- \r
- } else if (Info_h.biBitCount == 8 && Info_h.biCompression == 1) { \r
- table_R = (unsigned char *) malloc(256 * sizeof(unsigned char));\r
- table_G = (unsigned char *) malloc(256 * sizeof(unsigned char));\r
- table_B = (unsigned char *) malloc(256 * sizeof(unsigned char));\r
- \r
- for (j = 0; j < Info_h.biClrUsed; j++) {\r
- table_B[j] = getc(IN);\r
- table_G[j] = getc(IN);\r
- table_R[j] = getc(IN);\r
- getc(IN);\r
- if (table_R[j] != table_G[j] && table_R[j] != table_B[j] && table_G[j] != table_B[j])\r
- gray_scale = 0;\r
- }\r
-\r
- numcomps = gray_scale ? 1 : 3;\r
- color_space = gray_scale ? CLRSPC_GRAY : CLRSPC_SRGB;\r
- /* initialize image components */\r
- memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));\r
- for(i = 0; i < numcomps; i++) {\r
- cmptparm[i].prec = 8;\r
- cmptparm[i].bpp = 8;\r
- cmptparm[i].sgnd = 0;\r
- cmptparm[i].dx = subsampling_dx;\r
- cmptparm[i].dy = subsampling_dy;\r
- cmptparm[i].w = w;\r
- cmptparm[i].h = h;\r
- }\r
- /* create the image */\r
- image = opj_image_create(numcomps, &cmptparm[0], color_space);\r
- if(!image) {\r
- fclose(IN);\r
- return NULL;\r
- }\r
-\r
- /* set image offset and reference grid */\r
- image->x0 = parameters->image_offset_x0;\r
- image->y0 = parameters->image_offset_y0;\r
- image->x1 = !image->x0 ? (w - 1) * subsampling_dx + 1 : image->x0 + (w - 1) * subsampling_dx + 1;\r
- image->y1 = !image->y0 ? (h - 1) * subsampling_dy + 1 : image->y0 + (h - 1) * subsampling_dy + 1;\r
-\r
- /* set image data */\r
- \r
- /* Place the cursor at the beginning of the image information */\r
- fseek(IN, 0, SEEK_SET);\r
- fseek(IN, File_h.bfOffBits, SEEK_SET);\r
- \r
- RGB = (unsigned char *) malloc(Info_h.biWidth * Info_h.biHeight * sizeof(unsigned char));\r
- \r
- while (not_end_file) {\r
- v = getc(IN);\r
- if (v) {\r
- v2 = getc(IN);\r
- for (i = 0; i < (int) v; i++) {\r
- RGB[line * Info_h.biWidth + col] = v2;\r
- col++;\r
- }\r
- } else {\r
- v = getc(IN);\r
- switch (v) {\r
- case 0:\r
- col = 0;\r
- line++;\r
- break;\r
- case 1:\r
- line++;\r
- not_end_file = 0;\r
- break;\r
- case 2:\r
- fprintf(stderr,"No Delta supported\n");\r
- opj_image_destroy(image);\r
- fclose(IN);\r
- return NULL;\r
- break;\r
- default:\r
- for (i = 0; i < v; i++) {\r
- v2 = getc(IN);\r
- RGB[line * Info_h.biWidth + col] = v2;\r
- col++;\r
- }\r
- if (v % 2)\r
- v2 = getc(IN);\r
- break;\r
- }\r
- }\r
- }\r
- if (gray_scale) {\r
- index = 0;\r
- for (line = 0; line < Info_h.biHeight; line++) {\r
- for (col = 0; col < Info_h.biWidth; col++) {\r
- image->comps[0].data[index] = table_R[(int)RGB[(Info_h.biHeight - line - 1) * Info_h.biWidth + col]];\r
- index++;\r
- }\r
- }\r
- } else {\r
- index = 0;\r
- for (line = 0; line < Info_h.biHeight; line++) {\r
- for (col = 0; col < Info_h.biWidth; col++) {\r
- unsigned char pixel_index = (int)RGB[(Info_h.biHeight - line - 1) * Info_h.biWidth + col];\r
- image->comps[0].data[index] = table_R[pixel_index];\r
- image->comps[1].data[index] = table_G[pixel_index];\r
- image->comps[2].data[index] = table_B[pixel_index];\r
- index++;\r
- }\r
- }\r
- }\r
- free(RGB);\r
- } else {\r
- fprintf(stderr, \r
- "Other system than 24 bits/pixels or 8 bits (no RLE coding) is not yet implemented [%d]\n", Info_h.biBitCount);\r
- }\r
- fclose(IN);\r
- }\r
- \r
- return image;\r
-}\r
-\r
-int imagetobmp(opj_image_t * image, char *outfile) {\r
- int w, wr, h, hr;\r
- int i, pad;\r
- FILE *fdest = NULL;\r
-\r
- if (image->numcomps == 3 && image->comps[0].dx == image->comps[1].dx\r
- && image->comps[1].dx == image->comps[2].dx\r
- && image->comps[0].dy == image->comps[1].dy\r
- && image->comps[1].dy == image->comps[2].dy\r
- && image->comps[0].prec == image->comps[1].prec\r
- && image->comps[1].prec == image->comps[2].prec) {\r
- \r
- /* -->> -->> -->> -->> \r
- 24 bits color \r
- <<-- <<-- <<-- <<-- */\r
- \r
- fdest = fopen(outfile, "wb");\r
- if (!fdest) {\r
- fprintf(stderr, "ERROR -> failed to open %s for writing\n", outfile);\r
- return 1;\r
- }\r
- \r
- /* w = int_ceildiv(image->x1 - image->x0, image->comps[0].dx); */\r
- /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor), image->comps[0].dx); */\r
- w = image->comps[0].w;\r
- wr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);\r
- \r
- /* h = int_ceildiv(image->y1 - image->y0, image->comps[0].dy); */\r
- /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[0].dy); */\r
- h = image->comps[0].h;\r
- hr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);\r
- \r
- fprintf(fdest, "BM");\r
- \r
- /* FILE HEADER */\r
- /* ------------- */\r
- fprintf(fdest, "%c%c%c%c",\r
- (unsigned char) (hr * wr * 3 + 3 * hr * (wr % 2) + 54) & 0xff,\r
- (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) >> 8) & 0xff,\r
- (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) >> 16) & 0xff,\r
- (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff,((54) >> 16) & 0xff, ((54) >> 24) & 0xff);\r
- \r
- /* INFO HEADER */\r
- /* ------------- */\r
- fprintf(fdest, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff, ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (unsigned char) ((wr) & 0xff),\r
- (unsigned char) ((wr) >> 8) & 0xff,\r
- (unsigned char) ((wr) >> 16) & 0xff,\r
- (unsigned char) ((wr) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (unsigned char) ((hr) & 0xff),\r
- (unsigned char) ((hr) >> 8) & 0xff,\r
- (unsigned char) ((hr) >> 16) & 0xff,\r
- (unsigned char) ((hr) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);\r
- fprintf(fdest, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (unsigned char) (3 * hr * wr + 3 * hr * (wr % 2)) & 0xff,\r
- (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> 8) & 0xff,\r
- (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> 16) & 0xff,\r
- (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);\r
- \r
- for (i = 0; i < wr * hr; i++) {\r
- unsigned char R, G, B;\r
- /* a modifier */\r
- /* R = image->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; */\r
- R = image->comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];\r
- /* G = image->comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; */\r
- G = image->comps[1].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];\r
- /* B = image->comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; */\r
- B = image->comps[2].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];\r
- fprintf(fdest, "%c%c%c", B, G, R);\r
- \r
- if ((i + 1) % wr == 0) {\r
- for (pad = (3 * wr) % 4 ? 4 - (3 * wr) % 4 : 0; pad > 0; pad--) /* ADD */\r
- fprintf(fdest, "%c", 0);\r
- }\r
- }\r
- fclose(fdest);\r
- } else { /* Gray-scale */\r
-\r
- /* -->> -->> -->> -->>\r
- 8 bits non code (Gray scale)\r
- <<-- <<-- <<-- <<-- */\r
-\r
- fdest = fopen(outfile, "wb");\r
- /* w = int_ceildiv(image->x1 - image->x0, image->comps[0].dx); */\r
- /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor), image->comps[0].dx); */\r
- w = image->comps[0].w;\r
- wr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);\r
- \r
- /* h = int_ceildiv(image->y1 - image->y0, image->comps[0].dy); */\r
- /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[0].dy); */\r
- h = image->comps[0].h;\r
- hr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);\r
- \r
- fprintf(fdest, "BM");\r
- \r
- /* FILE HEADER */\r
- /* ------------- */\r
- fprintf(fdest, "%c%c%c%c", (unsigned char) (hr * wr + 54 + 1024 + hr * (wr % 2)) & 0xff,\r
- (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2)) >> 8) & 0xff,\r
- (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2)) >> 16) & 0xff,\r
- (unsigned char) ((hr * wr + 54 + 1024 + wr * (wr % 2)) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (54 + 1024) & 0xff, ((54 + 1024) >> 8) & 0xff, \r
- ((54 + 1024) >> 16) & 0xff,\r
- ((54 + 1024) >> 24) & 0xff);\r
- \r
- /* INFO HEADER */\r
- /* ------------- */\r
- fprintf(fdest, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff, ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (unsigned char) ((wr) & 0xff),\r
- (unsigned char) ((wr) >> 8) & 0xff,\r
- (unsigned char) ((wr) >> 16) & 0xff,\r
- (unsigned char) ((wr) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (unsigned char) ((hr) & 0xff),\r
- (unsigned char) ((hr) >> 8) & 0xff,\r
- (unsigned char) ((hr) >> 16) & 0xff,\r
- (unsigned char) ((hr) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);\r
- fprintf(fdest, "%c%c", (8) & 0xff, ((8) >> 8) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (unsigned char) (hr * wr + hr * (wr % 2)) & 0xff,\r
- (unsigned char) ((hr * wr + hr * (wr % 2)) >> 8) & 0xff,\r
- (unsigned char) ((hr * wr + hr * (wr % 2)) >> 16) & 0xff,\r
- (unsigned char) ((hr * wr + hr * (wr % 2)) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff, ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);\r
- fprintf(fdest, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff, ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);\r
- }\r
-\r
- for (i = 0; i < 256; i++) {\r
- fprintf(fdest, "%c%c%c%c", i, i, i, 0);\r
- }\r
-\r
- for (i = 0; i < wr * hr; i++) {\r
- /* a modifier !! */\r
- /* fprintf(fdest, "%c", image->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]); */\r
- fprintf(fdest, "%c", image->comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)]);\r
- /*if (((i + 1) % w == 0 && w % 2))\r
- fprintf(fdest, "%c", 0); */\r
- if ((i + 1) % wr == 0) {\r
- for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--) /* ADD */\r
- fprintf(fdest, "%c", 0);\r
- }\r
- }\r
- fclose(fdest);\r
-\r
- return 0;\r
-}\r
-\r
-/* -->> -->> -->> -->>\r
-\r
-PGX IMAGE FORMAT\r
-\r
-<<-- <<-- <<-- <<-- */\r
-\r
-\r
-unsigned char readuchar(FILE * f)\r
-{\r
- unsigned char c1;\r
- fread(&c1, 1, 1, f);\r
- return c1;\r
-}\r
-\r
-unsigned short readushort(FILE * f, int bigendian)\r
-{\r
- unsigned char c1, c2;\r
- fread(&c1, 1, 1, f);\r
- fread(&c2, 1, 1, f);\r
- if (bigendian)\r
- return (c1 << 8) + c2;\r
- else\r
- return (c2 << 8) + c1;\r
-}\r
-\r
-unsigned int readuint(FILE * f, int bigendian)\r
-{\r
- unsigned char c1, c2, c3, c4;\r
- fread(&c1, 1, 1, f);\r
- fread(&c2, 1, 1, f);\r
- fread(&c3, 1, 1, f);\r
- fread(&c4, 1, 1, f);\r
- if (bigendian)\r
- return (c1 << 24) + (c2 << 16) + (c3 << 8) + c4;\r
- else\r
- return (c4 << 24) + (c3 << 16) + (c2 << 8) + c1;\r
-}\r
-\r
-opj_image_t* pgxtoimage(char *filename, opj_cparameters_t *parameters) {\r
- FILE *f = NULL;\r
- int w, h, prec;\r
- int i, numcomps, max;\r
- OPJ_COLOR_SPACE color_space;\r
- opj_image_cmptparm_t cmptparm; /* maximum of 1 component */\r
- opj_image_t * image = NULL;\r
-\r
- char endian1,endian2,sign;\r
- char signtmp[32];\r
-\r
- char temp[32];\r
- int bigendian;\r
- opj_image_comp_t *comp = NULL;\r
-\r
- numcomps = 1;\r
- color_space = CLRSPC_GRAY;\r
-\r
- memset(&cmptparm, 0, sizeof(opj_image_cmptparm_t));\r
-\r
- max = 0;\r
-\r
- f = fopen(filename, "rb");\r
- if (!f) {\r
- fprintf(stderr, "Failed to open %s for reading !\n", filename);\r
- return NULL;\r
- }\r
-\r
- fseek(f, 0, SEEK_SET);\r
- fscanf(f, "PG%[ \t]%c%c%[ \t+-]%d%[ \t]%d%[ \t]%d",temp,&endian1,&endian2,signtmp,&prec,temp,&w,temp,&h);\r
- \r
- i=0;\r
- sign='+'; \r
- while (signtmp[i]!='\0') {\r
- if (signtmp[i]=='-') sign='-';\r
- i++;\r
- }\r
- \r
- fgetc(f);\r
- if (endian1=='M' && endian2=='L') {\r
- bigendian = 1;\r
- } else if (endian2=='M' && endian1=='L') {\r
- bigendian = 0;\r
- } else {\r
- fprintf(stderr, "Bad pgx header, please check input file\n");\r
- return NULL;\r
- }\r
-\r
- /* initialize image component */\r
-\r
- cmptparm.x0 = parameters->image_offset_x0;\r
- cmptparm.y0 = parameters->image_offset_y0;\r
- cmptparm.w = !cmptparm.x0 ? (w - 1) * parameters->subsampling_dx + 1 : cmptparm.x0 + (w - 1) * parameters->subsampling_dx + 1;\r
- cmptparm.h = !cmptparm.y0 ? (h - 1) * parameters->subsampling_dy + 1 : cmptparm.y0 + (h - 1) * parameters->subsampling_dy + 1;\r
- \r
- if (sign == '-') {\r
- cmptparm.sgnd = 1;\r
- } else {\r
- cmptparm.sgnd = 0;\r
- }\r
- cmptparm.prec = prec;\r
- cmptparm.bpp = prec;\r
- cmptparm.dx = parameters->subsampling_dx;\r
- cmptparm.dy = parameters->subsampling_dy;\r
- \r
- /* create the image */\r
- image = opj_image_create(numcomps, &cmptparm, color_space);\r
- if(!image) {\r
- fclose(f);\r
- return NULL;\r
- }\r
- /* set image offset and reference grid */\r
- image->x0 = cmptparm.x0;\r
- image->y0 = cmptparm.x0;\r
- image->x1 = cmptparm.w;\r
- image->y1 = cmptparm.h;\r
-\r
- /* set image data */\r
-\r
- comp = &image->comps[0];\r
-\r
- for (i = 0; i < w * h; i++) {\r
- int v;\r
- if (comp->prec <= 8) {\r
- if (!comp->sgnd) {\r
- v = readuchar(f);\r
- } else {\r
- v = (char) readuchar(f);\r
- }\r
- } else if (comp->prec <= 16) {\r
- if (!comp->sgnd) {\r
- v = readushort(f, bigendian);\r
- } else {\r
- v = (short) readushort(f, bigendian);\r
- }\r
- } else {\r
- if (!comp->sgnd) {\r
- v = readuint(f, bigendian);\r
- } else {\r
- v = (int) readuint(f, bigendian);\r
- }\r
- }\r
- if (v > max)\r
- max = v;\r
- comp->data[i] = v;\r
- }\r
- fclose(f);\r
- comp->bpp = int_floorlog2(max) + 1;\r
-\r
- return image;\r
-}\r
-\r
-int imagetopgx(opj_image_t * image, char *outfile) {\r
- int w, wr, h, hr;\r
- int i, j, compno;\r
- FILE *fdest = NULL;\r
-\r
- for (compno = 0; compno < image->numcomps; compno++) {\r
- opj_image_comp_t *comp = &image->comps[compno];\r
- char name[256];\r
- int nbytes = 0;\r
- char *tmp = outfile;\r
- while (*tmp) {\r
- tmp++;\r
- }\r
- while (*tmp!='.') {\r
- tmp--;\r
- }\r
- *tmp='\0';\r
-\r
- if (image->numcomps > 1) {\r
- sprintf(name, "%s-%d.pgx", outfile, compno);\r
- } else {\r
- sprintf(name, "%s.pgx", outfile);\r
- }\r
- fdest = fopen(name, "wb");\r
- if (!fdest) {\r
- fprintf(stderr, "ERROR -> failed to open %s for writing\n", name);\r
- return 1;\r
- }\r
- /* w = int_ceildiv(image->x1 - image->x0, comp->dx); */\r
- /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor), comp->dx); */\r
- w = image->comps[compno].w;\r
- wr = int_ceildivpow2(image->comps[compno].w, image->comps[compno].factor);\r
- \r
- /* h = int_ceildiv(image->y1 - image->y0, comp->dy); */\r
- /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), comp->dy); */\r
- h = image->comps[compno].h;\r
- hr = int_ceildivpow2(image->comps[compno].h, image->comps[compno].factor);\r
- \r
- fprintf(fdest, "PG ML %c %d %d %d\n", comp->sgnd ? '-' : '+', comp->prec, wr, hr);\r
- if (comp->prec <= 8) {\r
- nbytes = 1;\r
- } else if (comp->prec <= 16) {\r
- nbytes = 2;\r
- } else {\r
- nbytes = 4;\r
- }\r
- for (i = 0; i < wr * hr; i++) {\r
- int v = image->comps[compno].data[i / wr * w + i % wr];\r
- for (j = nbytes - 1; j >= 0; j--) {\r
- char byte = (char) (v >> (j * 8));\r
- fwrite(&byte, 1, 1, fdest);\r
- }\r
- }\r
- fclose(fdest);\r
- }\r
-\r
- return 0;\r
-}\r
-\r
-/* -->> -->> -->> -->>\r
-\r
-PNM IMAGE FORMAT\r
-\r
-<<-- <<-- <<-- <<-- */\r
-\r
-opj_image_t* pnmtoimage(char *filename, opj_cparameters_t *parameters) {\r
- int subsampling_dx = parameters->subsampling_dx;\r
- int subsampling_dy = parameters->subsampling_dy;\r
-\r
- FILE *f = NULL;\r
- int i, compno, numcomps, w, h;\r
- OPJ_COLOR_SPACE color_space;\r
- opj_image_cmptparm_t cmptparm[3]; /* maximum of 3 components */\r
- opj_image_t * image = NULL;\r
- char value;\r
- char comment[256];\r
-\r
- f = fopen(filename, "rb");\r
- if (!f) {\r
- fprintf(stderr, "\033[0;33mFailed to open %s for reading !!\033[0;39m\n", filename);\r
- return 0;\r
- }\r
-\r
- if (fgetc(f) != 'P')\r
- return 0;\r
- value = fgetc(f);\r
-\r
- switch(value) {\r
- case '2': /* greyscale image type */\r
- case '5':\r
- {\r
- numcomps = 1;\r
- color_space = CLRSPC_GRAY;\r
-\r
- fgetc(f);\r
-\r
- if (fgetc(f) == '#') {\r
- /* skip comments */\r
- fseek(f, 0, SEEK_SET);\r
- if (value == '2') {\r
- fscanf(f, "P2\n");\r
- } else if (value == '5') {\r
- fscanf(f, "P5\n");\r
- }\r
- fgets(comment, 256, f);\r
- fscanf(f, "%d %d\n255", &w, &h);\r
- } else {\r
- fseek(f, 0, SEEK_SET);\r
- if (value == '2') {\r
- fscanf(f, "P2\n%d %d\n255", &w, &h);\r
- } else if (value == '5') {\r
- fscanf(f, "P5\n%d %d\n255", &w, &h);\r
- }\r
- }\r
- \r
- fgetc(f); /* <cr><lf> */\r
- }\r
- break;\r
-\r
- case '3': /* RGB image type */\r
- case '6':\r
- {\r
- numcomps = 3;\r
- color_space = CLRSPC_SRGB;\r
-\r
- fgetc(f);\r
-\r
- if (fgetc(f) == '#') {\r
- /* skip comments */\r
- fseek(f, 0, SEEK_SET);\r
- if (value == '3') {\r
- fscanf(f, "P3\n");\r
- } else if (value == '6') {\r
- fscanf(f, "P6\n");\r
- }\r
- fgets(comment, 256, f);\r
- fscanf(f, "%d %d\n255", &w, &h);\r
- } else {\r
- fseek(f, 0, SEEK_SET);\r
- if (value == '3') {\r
- fscanf(f, "P3\n%d %d\n255", &w, &h);\r
- } else if (value == '6') {\r
- fscanf(f, "P6\n%d %d\n255", &w, &h);\r
- }\r
- }\r
- \r
- fgetc(f); /* <cr><lf> */\r
- }\r
- break;\r
-\r
- default:\r
- fclose(f);\r
- return NULL;\r
- }\r
-\r
- /* initialize image components */\r
- memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));\r
- for(i = 0; i < numcomps; i++) {\r
- cmptparm[i].prec = 8;\r
- cmptparm[i].bpp = 8;\r
- cmptparm[i].sgnd = 0;\r
- cmptparm[i].dx = subsampling_dx;\r
- cmptparm[i].dy = subsampling_dy;\r
- cmptparm[i].w = w;\r
- cmptparm[i].h = h;\r
- }\r
- /* create the image */\r
- image = opj_image_create(numcomps, &cmptparm[0], color_space);\r
- if(!image) {\r
- fclose(f);\r
- return NULL;\r
- }\r
-\r
- /* set image offset and reference grid */\r
- image->x0 = parameters->image_offset_x0;\r
- image->y0 = parameters->image_offset_y0;\r
- image->x1 = parameters->image_offset_x0 + (w - 1) * subsampling_dx + 1;\r
- image->y1 = parameters->image_offset_y0 + (h - 1) * subsampling_dy + 1;\r
-\r
- /* set image data */\r
-\r
- if ((value == '2') || (value == '3')) { /* ASCII */\r
- for (i = 0; i < w * h; i++) {\r
- for(compno = 0; compno < numcomps; compno++) {\r
- unsigned int index = 0;\r
- fscanf(f, "%u", &index);\r
- /* compno : 0 = GREY, (0, 1, 2) = (R, G, B) */\r
- image->comps[compno].data[i] = index;\r
- }\r
- }\r
- } else if ((value == '5') || (value == '6')) { /* BINARY */\r
- for (i = 0; i < w * h; i++) {\r
- for(compno = 0; compno < numcomps; compno++) {\r
- unsigned char index = 0;\r
- fread(&index, 1, 1, f);\r
- /* compno : 0 = GREY, (0, 1, 2) = (R, G, B) */\r
- image->comps[compno].data[i] = index;\r
- }\r
- }\r
- }\r
-\r
- fclose(f);\r
-\r
- return image;\r
-}\r
-\r
-int imagetopnm(opj_image_t * image, char *outfile) {\r
- int w, wr, wrr, h, hr, hrr, max;\r
- int i, compno;\r
- int adjust;\r
- FILE *fdest = NULL;\r
- char S2;\r
- char *tmp = outfile;\r
-\r
- while (*tmp) {\r
- tmp++;\r
- }\r
- tmp--;\r
- tmp--;\r
- S2 = *tmp;\r
-\r
- if (image->numcomps == 3 && image->comps[0].dx == image->comps[1].dx\r
- && image->comps[1].dx == image->comps[2].dx\r
- && image->comps[0].dy == image->comps[1].dy\r
- && image->comps[1].dy == image->comps[2].dy\r
- && image->comps[0].prec == image->comps[1].prec\r
- && image->comps[1].prec == image->comps[2].prec\r
- && S2 !='g' && S2 !='G') {\r
-\r
- fdest = fopen(outfile, "wb");\r
- if (!fdest) {\r
- fprintf(stderr, "ERROR -> failed to open %s for writing\n", outfile);\r
- return 1;\r
- }\r
-\r
- w = int_ceildiv(image->x1 - image->x0, image->comps[0].dx);\r
- /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor),image->comps[0].dx); */\r
- wr = image->comps[0].w;\r
- wrr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);\r
- \r
- h = int_ceildiv(image->y1 - image->y0, image->comps[0].dy);\r
- /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[0].dy); */\r
- hr = image->comps[0].h;\r
- hrr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);\r
- \r
- max = image->comps[0].prec > 8 ? 255 : (1 << image->comps[0].prec) - 1;\r
- \r
- image->comps[0].x0 = int_ceildivpow2(image->comps[0].x0 - int_ceildiv(image->x0, image->comps[0].dx), image->comps[0].factor);\r
- image->comps[0].y0 = int_ceildivpow2(image->comps[0].y0 - int_ceildiv(image->y0, image->comps[0].dy), image->comps[0].factor);\r
-\r
- fprintf(fdest, "P6\n%d %d\n%d\n", wrr, hrr, max);\r
- adjust = image->comps[0].prec > 8 ? image->comps[0].prec - 8 : 0;\r
- for (i = 0; i < wrr * hrr; i++) {\r
- int r, g, b;\r
- unsigned char rc,gc,bc;\r
- r = image->comps[0].data[i / wrr * wr + i % wrr];\r
- r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);\r
- rc = (unsigned char) ((r >> adjust)+((r >> (adjust-1))%2));\r
-\r
- g = image->comps[1].data[i / wrr * wr + i % wrr];\r
- g += (image->comps[1].sgnd ? 1 << (image->comps[1].prec - 1) : 0);\r
- gc = (unsigned char) ((g >> adjust)+((g >> (adjust-1))%2));\r
- \r
- b = image->comps[2].data[i / wrr * wr + i % wrr];\r
- b += (image->comps[2].sgnd ? 1 << (image->comps[2].prec - 1) : 0);\r
- bc = (unsigned char) ((b >> adjust)+((b >> (adjust-1))%2));\r
- \r
- fprintf(fdest, "%c%c%c", rc, gc, bc);\r
- }\r
- fclose(fdest);\r
- } else {\r
- int ncomp=(S2=='g' || S2=='G')?1:image->numcomps;\r
- if (image->numcomps>ncomp) {\r
- fprintf(stderr,"WARNING -> [PGM files] Only the first component\n");\r
- fprintf(stderr," is written to the file\n");\r
- }\r
- for (compno = 0; compno < ncomp; compno++) {\r
- char name[256];\r
- if (ncomp > 1) {\r
- sprintf(name, "%d.%s", compno, outfile);\r
- } else {\r
- sprintf(name, "%s", outfile);\r
- }\r
- \r
- fdest = fopen(name, "wb");\r
- if (!fdest) {\r
- fprintf(stderr, "ERROR -> failed to open %s for writing\n", name);\r
- return 1;\r
- }\r
- \r
- w = int_ceildiv(image->x1 - image->x0, image->comps[compno].dx);\r
- /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor),image->comps[compno].dx); */\r
- wr = image->comps[compno].w;\r
- wrr = int_ceildivpow2(image->comps[compno].w, image->comps[compno].factor);\r
- \r
- h = int_ceildiv(image->y1 - image->y0, image->comps[compno].dy);\r
- /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[compno].dy); */\r
- hr = image->comps[compno].h;\r
- hrr = int_ceildivpow2(image->comps[compno].h, image->comps[compno].factor);\r
- \r
- max = image->comps[compno].prec > 8 ? 255 : (1 << image->comps[compno].prec) - 1;\r
- \r
- image->comps[compno].x0 = int_ceildivpow2(image->comps[compno].x0 - int_ceildiv(image->x0, image->comps[compno].dx), image->comps[compno].factor);\r
- image->comps[compno].y0 = int_ceildivpow2(image->comps[compno].y0 - int_ceildiv(image->y0, image->comps[compno].dy), image->comps[compno].factor);\r
- \r
- fprintf(fdest, "P5\n%d %d\n%d\n", wrr, hrr, max);\r
- adjust = image->comps[compno].prec > 8 ? image->comps[compno].prec - 8 : 0;\r
-\r
- for (i = 0; i < wrr * hrr; i++) {\r
- int l;\r
- unsigned char lc;\r
- l = image->comps[compno].data[i / wrr * wr + i % wrr];\r
- l += (image->comps[compno].sgnd ? 1 << (image->comps[compno].prec - 1) : 0);\r
- lc = (unsigned char) ((l >> adjust)+((l >> (adjust-1))%2));\r
- fprintf(fdest, "%c", lc);\r
- }\r
- fclose(fdest);\r
- }\r
- }\r
-\r
- return 0;\r
-}\r
+/*
+ * 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:
+ * 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 <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include "openjpeg.h"
+
+/*
+ * Get logarithm of an integer and round downwards.
+ *
+ * log2(a)
+ */
+static int int_floorlog2(int a) {
+ int l;
+ for (l = 0; a > 1; l++) {
+ a >>= 1;
+ }
+ return l;
+}
+
+/*
+ * Divide an integer by a power of 2 and round upwards.
+ *
+ * a divided by 2^b
+ */
+static int int_ceildivpow2(int a, int b) {
+ return (a + (1 << b) - 1) >> b;
+}
+
+/*
+ * Divide an integer and round upwards.
+ *
+ * a divided by b
+ */
+static int int_ceildiv(int a, int b) {
+ return (a + b - 1) / b;
+}
+
+/* -->> -->> -->> -->>
+
+ BMP IMAGE FORMAT
+
+ <<-- <<-- <<-- <<-- */
+
+/* WORD defines a two byte word */
+typedef unsigned short int WORD;
+
+/* DWORD defines a four byte word */
+typedef unsigned long int DWORD;
+
+typedef struct {
+ WORD bfType; /* 'BM' for Bitmap (19776) */
+ DWORD bfSize; /* Size of the file */
+ WORD bfReserved1; /* Reserved : 0 */
+ WORD bfReserved2; /* Reserved : 0 */
+ DWORD bfOffBits; /* Offset */
+} BITMAPFILEHEADER_t;
+
+typedef struct {
+ DWORD biSize; /* Size of the structure in bytes */
+ DWORD biWidth; /* Width of the image in pixels */
+ DWORD biHeight; /* Heigth of the image in pixels */
+ WORD biPlanes; /* 1 */
+ WORD biBitCount; /* Number of color bits by pixels */
+ DWORD biCompression; /* Type of encoding 0: none 1: RLE8 2: RLE4 */
+ DWORD biSizeImage; /* Size of the image in bytes */
+ DWORD biXpelsPerMeter; /* Horizontal (X) resolution in pixels/meter */
+ DWORD biYpelsPerMeter; /* Vertical (Y) resolution in pixels/meter */
+ DWORD biClrUsed; /* Number of color used in the image (0: ALL) */
+ DWORD biClrImportant; /* Number of important color (0: ALL) */
+} BITMAPINFOHEADER_t;
+
+opj_image_t* bmptoimage(char *filename, opj_cparameters_t *parameters) {
+ int subsampling_dx = parameters->subsampling_dx;
+ int subsampling_dy = parameters->subsampling_dy;
+
+ int i, numcomps, w, h;
+ OPJ_COLOR_SPACE color_space;
+ opj_image_cmptparm_t cmptparm[3]; /* maximum of 3 components */
+ opj_image_t * image = NULL;
+
+ FILE *IN;
+ BITMAPFILEHEADER_t File_h;
+ BITMAPINFOHEADER_t Info_h;
+ unsigned char *RGB;
+ unsigned char *table_R, *table_G, *table_B;
+ unsigned int j, PAD = 0;
+
+ int x, y, index;
+ int gray_scale = 1, not_end_file = 1;
+
+ unsigned int line = 0, col = 0;
+ unsigned char v, v2;
+ DWORD W, H;
+
+ IN = fopen(filename, "rb");
+ if (!IN) {
+ fprintf(stderr, "\033[0;33mFailed to open %s for reading !!\033[0;39m\n", filename);
+ return 0;
+ }
+
+ File_h.bfType = getc(IN);
+ File_h.bfType = (getc(IN) << 8) + File_h.bfType;
+
+ if (File_h.bfType != 19778) {
+ fprintf(stderr,"Error, not a BMP file!\n");
+ return 0;
+ } else {
+ /* FILE HEADER */
+ /* ------------- */
+ File_h.bfSize = getc(IN);
+ File_h.bfSize = (getc(IN) << 8) + File_h.bfSize;
+ File_h.bfSize = (getc(IN) << 16) + File_h.bfSize;
+ File_h.bfSize = (getc(IN) << 24) + File_h.bfSize;
+
+ File_h.bfReserved1 = getc(IN);
+ File_h.bfReserved1 = (getc(IN) << 8) + File_h.bfReserved1;
+
+ File_h.bfReserved2 = getc(IN);
+ File_h.bfReserved2 = (getc(IN) << 8) + File_h.bfReserved2;
+
+ File_h.bfOffBits = getc(IN);
+ File_h.bfOffBits = (getc(IN) << 8) + File_h.bfOffBits;
+ File_h.bfOffBits = (getc(IN) << 16) + File_h.bfOffBits;
+ File_h.bfOffBits = (getc(IN) << 24) + File_h.bfOffBits;
+
+ /* INFO HEADER */
+ /* ------------- */
+
+ Info_h.biSize = getc(IN);
+ Info_h.biSize = (getc(IN) << 8) + Info_h.biSize;
+ Info_h.biSize = (getc(IN) << 16) + Info_h.biSize;
+ Info_h.biSize = (getc(IN) << 24) + Info_h.biSize;
+
+ Info_h.biWidth = getc(IN);
+ Info_h.biWidth = (getc(IN) << 8) + Info_h.biWidth;
+ Info_h.biWidth = (getc(IN) << 16) + Info_h.biWidth;
+ Info_h.biWidth = (getc(IN) << 24) + Info_h.biWidth;
+ w = Info_h.biWidth;
+
+ Info_h.biHeight = getc(IN);
+ Info_h.biHeight = (getc(IN) << 8) + Info_h.biHeight;
+ Info_h.biHeight = (getc(IN) << 16) + Info_h.biHeight;
+ Info_h.biHeight = (getc(IN) << 24) + Info_h.biHeight;
+ h = Info_h.biHeight;
+
+ Info_h.biPlanes = getc(IN);
+ Info_h.biPlanes = (getc(IN) << 8) + Info_h.biPlanes;
+
+ Info_h.biBitCount = getc(IN);
+ Info_h.biBitCount = (getc(IN) << 8) + Info_h.biBitCount;
+
+ Info_h.biCompression = getc(IN);
+ Info_h.biCompression = (getc(IN) << 8) + Info_h.biCompression;
+ Info_h.biCompression = (getc(IN) << 16) + Info_h.biCompression;
+ Info_h.biCompression = (getc(IN) << 24) + Info_h.biCompression;
+
+ Info_h.biSizeImage = getc(IN);
+ Info_h.biSizeImage = (getc(IN) << 8) + Info_h.biSizeImage;
+ Info_h.biSizeImage = (getc(IN) << 16) + Info_h.biSizeImage;
+ Info_h.biSizeImage = (getc(IN) << 24) + Info_h.biSizeImage;
+
+ Info_h.biXpelsPerMeter = getc(IN);
+ Info_h.biXpelsPerMeter = (getc(IN) << 8) + Info_h.biXpelsPerMeter;
+ Info_h.biXpelsPerMeter = (getc(IN) << 16) + Info_h.biXpelsPerMeter;
+ Info_h.biXpelsPerMeter = (getc(IN) << 24) + Info_h.biXpelsPerMeter;
+
+ Info_h.biYpelsPerMeter = getc(IN);
+ Info_h.biYpelsPerMeter = (getc(IN) << 8) + Info_h.biYpelsPerMeter;
+ Info_h.biYpelsPerMeter = (getc(IN) << 16) + Info_h.biYpelsPerMeter;
+ Info_h.biYpelsPerMeter = (getc(IN) << 24) + Info_h.biYpelsPerMeter;
+
+ Info_h.biClrUsed = getc(IN);
+ Info_h.biClrUsed = (getc(IN) << 8) + Info_h.biClrUsed;
+ Info_h.biClrUsed = (getc(IN) << 16) + Info_h.biClrUsed;
+ Info_h.biClrUsed = (getc(IN) << 24) + Info_h.biClrUsed;
+
+ Info_h.biClrImportant = getc(IN);
+ Info_h.biClrImportant = (getc(IN) << 8) + Info_h.biClrImportant;
+ Info_h.biClrImportant = (getc(IN) << 16) + Info_h.biClrImportant;
+ Info_h.biClrImportant = (getc(IN) << 24) + Info_h.biClrImportant;
+
+ /* Read the data and store them in the OUT file */
+
+ if (Info_h.biBitCount == 24) {
+ numcomps = 3;
+ color_space = CLRSPC_SRGB;
+ /* initialize image components */
+ memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));
+ for(i = 0; i < numcomps; i++) {
+ cmptparm[i].prec = 8;
+ cmptparm[i].bpp = 8;
+ cmptparm[i].sgnd = 0;
+ cmptparm[i].dx = subsampling_dx;
+ cmptparm[i].dy = subsampling_dy;
+ cmptparm[i].w = w;
+ cmptparm[i].h = h;
+ }
+ /* create the image */
+ image = opj_image_create(numcomps, &cmptparm[0], color_space);
+ if(!image) {
+ fclose(IN);
+ return NULL;
+ }
+
+ /* set image offset and reference grid */
+ image->x0 = parameters->image_offset_x0;
+ image->y0 = parameters->image_offset_y0;
+ image->x1 = !image->x0 ? (w - 1) * subsampling_dx + 1 : image->x0 + (w - 1) * subsampling_dx + 1;
+ image->y1 = !image->y0 ? (h - 1) * subsampling_dy + 1 : image->y0 + (h - 1) * subsampling_dy + 1;
+
+ /* set image data */
+
+ /* Place the cursor at the beginning of the image information */
+ fseek(IN, 0, SEEK_SET);
+ fseek(IN, File_h.bfOffBits, SEEK_SET);
+
+ W = Info_h.biWidth;
+ H = Info_h.biHeight;
+
+ /* PAD = 4 - (3 * W) % 4; */
+ /* PAD = (PAD == 4) ? 0 : PAD; */
+ PAD = (3 * W) % 4 ? 4 - (3 * W) % 4 : 0;
+
+ RGB = (unsigned char *) malloc((3 * W + PAD) * H * sizeof(unsigned char));
+
+ fread(RGB, sizeof(unsigned char), (3 * W + PAD) * H, IN);
+
+ index = 0;
+
+ for(y = 0; y < (int)H; y++) {
+ unsigned char *scanline = RGB + (3 * W + PAD) * (H - 1 - y);
+ for(x = 0; x < (int)W; x++) {
+ unsigned char *pixel = &scanline[3 * x];
+ image->comps[0].data[index] = pixel[2]; /* R */
+ image->comps[1].data[index] = pixel[1]; /* G */
+ image->comps[2].data[index] = pixel[0]; /* B */
+ index++;
+ }
+ }
+
+ free(RGB);
+
+ } else if (Info_h.biBitCount == 8 && Info_h.biCompression == 0) {
+ table_R = (unsigned char *) malloc(256 * sizeof(unsigned char));
+ table_G = (unsigned char *) malloc(256 * sizeof(unsigned char));
+ table_B = (unsigned char *) malloc(256 * sizeof(unsigned char));
+
+ for (j = 0; j < Info_h.biClrUsed; j++) {
+ table_B[j] = getc(IN);
+ table_G[j] = getc(IN);
+ table_R[j] = getc(IN);
+ getc(IN);
+ if (table_R[j] != table_G[j] && table_R[j] != table_B[j] && table_G[j] != table_B[j])
+ gray_scale = 0;
+ }
+
+ /* Place the cursor at the beginning of the image information */
+ fseek(IN, 0, SEEK_SET);
+ fseek(IN, File_h.bfOffBits, SEEK_SET);
+
+ W = Info_h.biWidth;
+ H = Info_h.biHeight;
+ if (Info_h.biWidth % 2)
+ W++;
+
+ numcomps = gray_scale ? 1 : 3;
+ color_space = gray_scale ? CLRSPC_GRAY : CLRSPC_SRGB;
+ /* initialize image components */
+ memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));
+ for(i = 0; i < numcomps; i++) {
+ cmptparm[i].prec = 8;
+ cmptparm[i].bpp = 8;
+ cmptparm[i].sgnd = 0;
+ cmptparm[i].dx = subsampling_dx;
+ cmptparm[i].dy = subsampling_dy;
+ cmptparm[i].w = w;
+ cmptparm[i].h = h;
+ }
+ /* create the image */
+ image = opj_image_create(numcomps, &cmptparm[0], color_space);
+ if(!image) {
+ fclose(IN);
+ return NULL;
+ }
+
+ /* set image offset and reference grid */
+ image->x0 = parameters->image_offset_x0;
+ image->y0 = parameters->image_offset_y0;
+ image->x1 = !image->x0 ? (w - 1) * subsampling_dx + 1 : image->x0 + (w - 1) * subsampling_dx + 1;
+ image->y1 = !image->y0 ? (h - 1) * subsampling_dy + 1 : image->y0 + (h - 1) * subsampling_dy + 1;
+
+ /* set image data */
+
+ RGB = (unsigned char *) malloc(W * H * sizeof(unsigned char));
+
+ fread(RGB, sizeof(unsigned char), W * H, IN);
+ if (gray_scale) {
+ index = 0;
+ for (j = 0; j < W * H; j++) {
+ if ((j % W < W - 1 && Info_h.biWidth % 2) || !(Info_h.biWidth % 2)) {
+ image->comps[0].data[index] = table_R[RGB[W * H - ((j) / (W) + 1) * W + (j) % (W)]];
+ index++;
+ }
+ }
+
+ } else {
+ index = 0;
+ for (j = 0; j < W * H; j++) {
+ if ((j % W < W - 1 && Info_h.biWidth % 2) || !(Info_h.biWidth % 2)) {
+ unsigned char pixel_index = RGB[W * H - ((j) / (W) + 1) * W + (j) % (W)];
+ image->comps[0].data[index] = table_R[pixel_index];
+ image->comps[1].data[index] = table_G[pixel_index];
+ image->comps[2].data[index] = table_B[pixel_index];
+ index++;
+ }
+ }
+ }
+ free(RGB);
+
+ } else if (Info_h.biBitCount == 8 && Info_h.biCompression == 1) {
+ table_R = (unsigned char *) malloc(256 * sizeof(unsigned char));
+ table_G = (unsigned char *) malloc(256 * sizeof(unsigned char));
+ table_B = (unsigned char *) malloc(256 * sizeof(unsigned char));
+
+ for (j = 0; j < Info_h.biClrUsed; j++) {
+ table_B[j] = getc(IN);
+ table_G[j] = getc(IN);
+ table_R[j] = getc(IN);
+ getc(IN);
+ if (table_R[j] != table_G[j] && table_R[j] != table_B[j] && table_G[j] != table_B[j])
+ gray_scale = 0;
+ }
+
+ numcomps = gray_scale ? 1 : 3;
+ color_space = gray_scale ? CLRSPC_GRAY : CLRSPC_SRGB;
+ /* initialize image components */
+ memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));
+ for(i = 0; i < numcomps; i++) {
+ cmptparm[i].prec = 8;
+ cmptparm[i].bpp = 8;
+ cmptparm[i].sgnd = 0;
+ cmptparm[i].dx = subsampling_dx;
+ cmptparm[i].dy = subsampling_dy;
+ cmptparm[i].w = w;
+ cmptparm[i].h = h;
+ }
+ /* create the image */
+ image = opj_image_create(numcomps, &cmptparm[0], color_space);
+ if(!image) {
+ fclose(IN);
+ return NULL;
+ }
+
+ /* set image offset and reference grid */
+ image->x0 = parameters->image_offset_x0;
+ image->y0 = parameters->image_offset_y0;
+ image->x1 = !image->x0 ? (w - 1) * subsampling_dx + 1 : image->x0 + (w - 1) * subsampling_dx + 1;
+ image->y1 = !image->y0 ? (h - 1) * subsampling_dy + 1 : image->y0 + (h - 1) * subsampling_dy + 1;
+
+ /* set image data */
+
+ /* Place the cursor at the beginning of the image information */
+ fseek(IN, 0, SEEK_SET);
+ fseek(IN, File_h.bfOffBits, SEEK_SET);
+
+ RGB = (unsigned char *) malloc(Info_h.biWidth * Info_h.biHeight * sizeof(unsigned char));
+
+ while (not_end_file) {
+ v = getc(IN);
+ if (v) {
+ v2 = getc(IN);
+ for (i = 0; i < (int) v; i++) {
+ RGB[line * Info_h.biWidth + col] = v2;
+ col++;
+ }
+ } else {
+ v = getc(IN);
+ switch (v) {
+ case 0:
+ col = 0;
+ line++;
+ break;
+ case 1:
+ line++;
+ not_end_file = 0;
+ break;
+ case 2:
+ fprintf(stderr,"No Delta supported\n");
+ opj_image_destroy(image);
+ fclose(IN);
+ return NULL;
+ break;
+ default:
+ for (i = 0; i < v; i++) {
+ v2 = getc(IN);
+ RGB[line * Info_h.biWidth + col] = v2;
+ col++;
+ }
+ if (v % 2)
+ v2 = getc(IN);
+ break;
+ }
+ }
+ }
+ if (gray_scale) {
+ index = 0;
+ for (line = 0; line < Info_h.biHeight; line++) {
+ for (col = 0; col < Info_h.biWidth; col++) {
+ image->comps[0].data[index] = table_R[(int)RGB[(Info_h.biHeight - line - 1) * Info_h.biWidth + col]];
+ index++;
+ }
+ }
+ } else {
+ index = 0;
+ for (line = 0; line < Info_h.biHeight; line++) {
+ for (col = 0; col < Info_h.biWidth; col++) {
+ unsigned char pixel_index = (int)RGB[(Info_h.biHeight - line - 1) * Info_h.biWidth + col];
+ image->comps[0].data[index] = table_R[pixel_index];
+ image->comps[1].data[index] = table_G[pixel_index];
+ image->comps[2].data[index] = table_B[pixel_index];
+ index++;
+ }
+ }
+ }
+ free(RGB);
+ } else {
+ fprintf(stderr,
+ "Other system than 24 bits/pixels or 8 bits (no RLE coding) is not yet implemented [%d]\n", Info_h.biBitCount);
+ }
+ fclose(IN);
+ }
+
+ return image;
+}
+
+int imagetobmp(opj_image_t * image, char *outfile) {
+ int w, wr, h, hr;
+ int i, pad;
+ FILE *fdest = NULL;
+
+ if (image->numcomps == 3 && image->comps[0].dx == image->comps[1].dx
+ && image->comps[1].dx == image->comps[2].dx
+ && image->comps[0].dy == image->comps[1].dy
+ && image->comps[1].dy == image->comps[2].dy
+ && image->comps[0].prec == image->comps[1].prec
+ && image->comps[1].prec == image->comps[2].prec) {
+
+ /* -->> -->> -->> -->>
+ 24 bits color
+ <<-- <<-- <<-- <<-- */
+
+ fdest = fopen(outfile, "wb");
+ if (!fdest) {
+ fprintf(stderr, "ERROR -> failed to open %s for writing\n", outfile);
+ return 1;
+ }
+
+ /* w = int_ceildiv(image->x1 - image->x0, image->comps[0].dx); */
+ /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor), image->comps[0].dx); */
+ w = image->comps[0].w;
+ wr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);
+
+ /* h = int_ceildiv(image->y1 - image->y0, image->comps[0].dy); */
+ /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[0].dy); */
+ h = image->comps[0].h;
+ hr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);
+
+ fprintf(fdest, "BM");
+
+ /* FILE HEADER */
+ /* ------------- */
+ fprintf(fdest, "%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(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff,((54) >> 16) & 0xff, ((54) >> 24) & 0xff);
+
+ /* INFO HEADER */
+ /* ------------- */
+ fprintf(fdest, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff, ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
+ (unsigned char) ((wr) >> 8) & 0xff,
+ (unsigned char) ((wr) >> 16) & 0xff,
+ (unsigned char) ((wr) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
+ (unsigned char) ((hr) >> 8) & 0xff,
+ (unsigned char) ((hr) >> 16) & 0xff,
+ (unsigned char) ((hr) >> 24) & 0xff);
+ fprintf(fdest, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
+ fprintf(fdest, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
+ fprintf(fdest, "%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(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
+ fprintf(fdest, "%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 = image->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; */
+ R = image->comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
+ /* G = image->comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; */
+ G = image->comps[1].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
+ /* B = image->comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]; */
+ B = image->comps[2].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
+ fprintf(fdest, "%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(fdest, "%c", 0);
+ }
+ }
+ fclose(fdest);
+ } else { /* Gray-scale */
+
+ /* -->> -->> -->> -->>
+ 8 bits non code (Gray scale)
+ <<-- <<-- <<-- <<-- */
+
+ fdest = fopen(outfile, "wb");
+ /* w = int_ceildiv(image->x1 - image->x0, image->comps[0].dx); */
+ /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor), image->comps[0].dx); */
+ w = image->comps[0].w;
+ wr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);
+
+ /* h = int_ceildiv(image->y1 - image->y0, image->comps[0].dy); */
+ /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[0].dy); */
+ h = image->comps[0].h;
+ hr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);
+
+ fprintf(fdest, "BM");
+
+ /* FILE HEADER */
+ /* ------------- */
+ fprintf(fdest, "%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(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (54 + 1024) & 0xff, ((54 + 1024) >> 8) & 0xff,
+ ((54 + 1024) >> 16) & 0xff,
+ ((54 + 1024) >> 24) & 0xff);
+
+ /* INFO HEADER */
+ /* ------------- */
+ fprintf(fdest, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff, ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
+ (unsigned char) ((wr) >> 8) & 0xff,
+ (unsigned char) ((wr) >> 16) & 0xff,
+ (unsigned char) ((wr) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
+ (unsigned char) ((hr) >> 8) & 0xff,
+ (unsigned char) ((hr) >> 16) & 0xff,
+ (unsigned char) ((hr) >> 24) & 0xff);
+ fprintf(fdest, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
+ fprintf(fdest, "%c%c", (8) & 0xff, ((8) >> 8) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff, ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
+ fprintf(fdest, "%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(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff, ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff, ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
+ fprintf(fdest, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff, ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
+ }
+
+ for (i = 0; i < 256; i++) {
+ fprintf(fdest, "%c%c%c%c", i, i, i, 0);
+ }
+
+ for (i = 0; i < wr * hr; i++) {
+ /* a modifier !! */
+ /* fprintf(fdest, "%c", image->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]); */
+ fprintf(fdest, "%c", image->comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)]);
+ /*if (((i + 1) % w == 0 && w % 2))
+ fprintf(fdest, "%c", 0); */
+ if ((i + 1) % wr == 0) {
+ for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--) /* ADD */
+ fprintf(fdest, "%c", 0);
+ }
+ }
+ fclose(fdest);
+
+ return 0;
+}
+
+/* -->> -->> -->> -->>
+
+PGX IMAGE FORMAT
+
+<<-- <<-- <<-- <<-- */
+
+
+unsigned char readuchar(FILE * f)
+{
+ unsigned char c1;
+ fread(&c1, 1, 1, f);
+ return c1;
+}
+
+unsigned short readushort(FILE * f, int bigendian)
+{
+ unsigned char c1, c2;
+ fread(&c1, 1, 1, f);
+ fread(&c2, 1, 1, f);
+ if (bigendian)
+ return (c1 << 8) + c2;
+ else
+ return (c2 << 8) + c1;
+}
+
+unsigned int readuint(FILE * f, int bigendian)
+{
+ unsigned char c1, c2, c3, c4;
+ fread(&c1, 1, 1, f);
+ fread(&c2, 1, 1, f);
+ fread(&c3, 1, 1, f);
+ fread(&c4, 1, 1, f);
+ if (bigendian)
+ return (c1 << 24) + (c2 << 16) + (c3 << 8) + c4;
+ else
+ return (c4 << 24) + (c3 << 16) + (c2 << 8) + c1;
+}
+
+opj_image_t* pgxtoimage(char *filename, opj_cparameters_t *parameters) {
+ FILE *f = NULL;
+ int w, h, prec;
+ int i, numcomps, max;
+ OPJ_COLOR_SPACE color_space;
+ opj_image_cmptparm_t cmptparm; /* maximum of 1 component */
+ opj_image_t * image = NULL;
+
+ char endian1,endian2,sign;
+ char signtmp[32];
+
+ char temp[32];
+ int bigendian;
+ opj_image_comp_t *comp = NULL;
+
+ numcomps = 1;
+ color_space = CLRSPC_GRAY;
+
+ memset(&cmptparm, 0, sizeof(opj_image_cmptparm_t));
+
+ max = 0;
+
+ f = fopen(filename, "rb");
+ if (!f) {
+ fprintf(stderr, "Failed to open %s for reading !\n", filename);
+ return NULL;
+ }
+
+ fseek(f, 0, SEEK_SET);
+ fscanf(f, "PG%[ \t]%c%c%[ \t+-]%d%[ \t]%d%[ \t]%d",temp,&endian1,&endian2,signtmp,&prec,temp,&w,temp,&h);
+
+ i=0;
+ sign='+';
+ while (signtmp[i]!='\0') {
+ if (signtmp[i]=='-') sign='-';
+ i++;
+ }
+
+ fgetc(f);
+ if (endian1=='M' && endian2=='L') {
+ bigendian = 1;
+ } else if (endian2=='M' && endian1=='L') {
+ bigendian = 0;
+ } else {
+ fprintf(stderr, "Bad pgx header, please check input file\n");
+ return NULL;
+ }
+
+ /* initialize image component */
+
+ cmptparm.x0 = parameters->image_offset_x0;
+ cmptparm.y0 = parameters->image_offset_y0;
+ cmptparm.w = !cmptparm.x0 ? (w - 1) * parameters->subsampling_dx + 1 : cmptparm.x0 + (w - 1) * parameters->subsampling_dx + 1;
+ cmptparm.h = !cmptparm.y0 ? (h - 1) * parameters->subsampling_dy + 1 : cmptparm.y0 + (h - 1) * parameters->subsampling_dy + 1;
+
+ if (sign == '-') {
+ cmptparm.sgnd = 1;
+ } else {
+ cmptparm.sgnd = 0;
+ }
+ cmptparm.prec = prec;
+ cmptparm.bpp = prec;
+ cmptparm.dx = parameters->subsampling_dx;
+ cmptparm.dy = parameters->subsampling_dy;
+
+ /* create the image */
+ image = opj_image_create(numcomps, &cmptparm, color_space);
+ if(!image) {
+ fclose(f);
+ return NULL;
+ }
+ /* set image offset and reference grid */
+ image->x0 = cmptparm.x0;
+ image->y0 = cmptparm.x0;
+ image->x1 = cmptparm.w;
+ image->y1 = cmptparm.h;
+
+ /* set image data */
+
+ comp = &image->comps[0];
+
+ for (i = 0; i < w * h; i++) {
+ int v;
+ if (comp->prec <= 8) {
+ if (!comp->sgnd) {
+ v = readuchar(f);
+ } else {
+ v = (char) readuchar(f);
+ }
+ } else if (comp->prec <= 16) {
+ if (!comp->sgnd) {
+ v = readushort(f, bigendian);
+ } else {
+ v = (short) readushort(f, bigendian);
+ }
+ } else {
+ if (!comp->sgnd) {
+ v = readuint(f, bigendian);
+ } else {
+ v = (int) readuint(f, bigendian);
+ }
+ }
+ if (v > max)
+ max = v;
+ comp->data[i] = v;
+ }
+ fclose(f);
+ comp->bpp = int_floorlog2(max) + 1;
+
+ return image;
+}
+
+int imagetopgx(opj_image_t * image, char *outfile) {
+ int w, wr, h, hr;
+ int i, j, compno;
+ FILE *fdest = NULL;
+
+ for (compno = 0; compno < image->numcomps; compno++) {
+ opj_image_comp_t *comp = &image->comps[compno];
+ char name[256];
+ int nbytes = 0;
+ char *tmp = outfile;
+ while (*tmp) {
+ tmp++;
+ }
+ while (*tmp!='.') {
+ tmp--;
+ }
+ *tmp='\0';
+
+ if (image->numcomps > 1) {
+ sprintf(name, "%s-%d.pgx", outfile, compno);
+ } else {
+ sprintf(name, "%s.pgx", outfile);
+ }
+ fdest = fopen(name, "wb");
+ if (!fdest) {
+ fprintf(stderr, "ERROR -> failed to open %s for writing\n", name);
+ return 1;
+ }
+ /* w = int_ceildiv(image->x1 - image->x0, comp->dx); */
+ /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor), comp->dx); */
+ w = image->comps[compno].w;
+ wr = int_ceildivpow2(image->comps[compno].w, image->comps[compno].factor);
+
+ /* h = int_ceildiv(image->y1 - image->y0, comp->dy); */
+ /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), comp->dy); */
+ h = image->comps[compno].h;
+ hr = int_ceildivpow2(image->comps[compno].h, image->comps[compno].factor);
+
+ fprintf(fdest, "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 = image->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, fdest);
+ }
+ }
+ fclose(fdest);
+ }
+
+ return 0;
+}
+
+/* -->> -->> -->> -->>
+
+PNM IMAGE FORMAT
+
+<<-- <<-- <<-- <<-- */
+
+opj_image_t* pnmtoimage(char *filename, opj_cparameters_t *parameters) {
+ int subsampling_dx = parameters->subsampling_dx;
+ int subsampling_dy = parameters->subsampling_dy;
+
+ FILE *f = NULL;
+ int i, compno, numcomps, w, h;
+ OPJ_COLOR_SPACE color_space;
+ opj_image_cmptparm_t cmptparm[3]; /* maximum of 3 components */
+ opj_image_t * image = NULL;
+ char value;
+ char comment[256];
+
+ f = fopen(filename, "rb");
+ if (!f) {
+ fprintf(stderr, "\033[0;33mFailed to open %s for reading !!\033[0;39m\n", filename);
+ return 0;
+ }
+
+ if (fgetc(f) != 'P')
+ return 0;
+ value = fgetc(f);
+
+ switch(value) {
+ case '2': /* greyscale image type */
+ case '5':
+ {
+ numcomps = 1;
+ color_space = CLRSPC_GRAY;
+
+ fgetc(f);
+
+ if (fgetc(f) == '#') {
+ /* skip comments */
+ fseek(f, 0, SEEK_SET);
+ if (value == '2') {
+ fscanf(f, "P2\n");
+ } else if (value == '5') {
+ fscanf(f, "P5\n");
+ }
+ fgets(comment, 256, f);
+ fscanf(f, "%d %d\n255", &w, &h);
+ } else {
+ fseek(f, 0, SEEK_SET);
+ if (value == '2') {
+ fscanf(f, "P2\n%d %d\n255", &w, &h);
+ } else if (value == '5') {
+ fscanf(f, "P5\n%d %d\n255", &w, &h);
+ }
+ }
+
+ fgetc(f); /* <cr><lf> */
+ }
+ break;
+
+ case '3': /* RGB image type */
+ case '6':
+ {
+ numcomps = 3;
+ color_space = CLRSPC_SRGB;
+
+ fgetc(f);
+
+ if (fgetc(f) == '#') {
+ /* skip comments */
+ fseek(f, 0, SEEK_SET);
+ if (value == '3') {
+ fscanf(f, "P3\n");
+ } else if (value == '6') {
+ fscanf(f, "P6\n");
+ }
+ fgets(comment, 256, f);
+ fscanf(f, "%d %d\n255", &w, &h);
+ } else {
+ fseek(f, 0, SEEK_SET);
+ if (value == '3') {
+ fscanf(f, "P3\n%d %d\n255", &w, &h);
+ } else if (value == '6') {
+ fscanf(f, "P6\n%d %d\n255", &w, &h);
+ }
+ }
+
+ fgetc(f); /* <cr><lf> */
+ }
+ break;
+
+ default:
+ fclose(f);
+ return NULL;
+ }
+
+ /* initialize image components */
+ memset(&cmptparm[0], 0, 3 * sizeof(opj_image_cmptparm_t));
+ for(i = 0; i < numcomps; i++) {
+ cmptparm[i].prec = 8;
+ cmptparm[i].bpp = 8;
+ cmptparm[i].sgnd = 0;
+ cmptparm[i].dx = subsampling_dx;
+ cmptparm[i].dy = subsampling_dy;
+ cmptparm[i].w = w;
+ cmptparm[i].h = h;
+ }
+ /* create the image */
+ image = opj_image_create(numcomps, &cmptparm[0], color_space);
+ if(!image) {
+ fclose(f);
+ return NULL;
+ }
+
+ /* set image offset and reference grid */
+ image->x0 = parameters->image_offset_x0;
+ image->y0 = parameters->image_offset_y0;
+ image->x1 = parameters->image_offset_x0 + (w - 1) * subsampling_dx + 1;
+ image->y1 = parameters->image_offset_y0 + (h - 1) * subsampling_dy + 1;
+
+ /* set image data */
+
+ if ((value == '2') || (value == '3')) { /* ASCII */
+ for (i = 0; i < w * h; i++) {
+ for(compno = 0; compno < numcomps; compno++) {
+ unsigned int index = 0;
+ fscanf(f, "%u", &index);
+ /* compno : 0 = GREY, (0, 1, 2) = (R, G, B) */
+ image->comps[compno].data[i] = index;
+ }
+ }
+ } else if ((value == '5') || (value == '6')) { /* BINARY */
+ for (i = 0; i < w * h; i++) {
+ for(compno = 0; compno < numcomps; compno++) {
+ unsigned char index = 0;
+ fread(&index, 1, 1, f);
+ /* compno : 0 = GREY, (0, 1, 2) = (R, G, B) */
+ image->comps[compno].data[i] = index;
+ }
+ }
+ }
+
+ fclose(f);
+
+ return image;
+}
+
+int imagetopnm(opj_image_t * image, char *outfile) {
+ int w, wr, wrr, h, hr, hrr, max;
+ int i, compno;
+ int adjust;
+ FILE *fdest = NULL;
+ char S2;
+ char *tmp = outfile;
+
+ while (*tmp) {
+ tmp++;
+ }
+ tmp--;
+ tmp--;
+ S2 = *tmp;
+
+ if (image->numcomps == 3 && image->comps[0].dx == image->comps[1].dx
+ && image->comps[1].dx == image->comps[2].dx
+ && image->comps[0].dy == image->comps[1].dy
+ && image->comps[1].dy == image->comps[2].dy
+ && image->comps[0].prec == image->comps[1].prec
+ && image->comps[1].prec == image->comps[2].prec
+ && S2 !='g' && S2 !='G') {
+
+ fdest = fopen(outfile, "wb");
+ if (!fdest) {
+ fprintf(stderr, "ERROR -> failed to open %s for writing\n", outfile);
+ return 1;
+ }
+
+ w = int_ceildiv(image->x1 - image->x0, image->comps[0].dx);
+ /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor),image->comps[0].dx); */
+ wr = image->comps[0].w;
+ wrr = int_ceildivpow2(image->comps[0].w, image->comps[0].factor);
+
+ h = int_ceildiv(image->y1 - image->y0, image->comps[0].dy);
+ /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[0].dy); */
+ hr = image->comps[0].h;
+ hrr = int_ceildivpow2(image->comps[0].h, image->comps[0].factor);
+
+ max = image->comps[0].prec > 8 ? 255 : (1 << image->comps[0].prec) - 1;
+
+ image->comps[0].x0 = int_ceildivpow2(image->comps[0].x0 - int_ceildiv(image->x0, image->comps[0].dx), image->comps[0].factor);
+ image->comps[0].y0 = int_ceildivpow2(image->comps[0].y0 - int_ceildiv(image->y0, image->comps[0].dy), image->comps[0].factor);
+
+ fprintf(fdest, "P6\n%d %d\n%d\n", wrr, hrr, max);
+ adjust = image->comps[0].prec > 8 ? image->comps[0].prec - 8 : 0;
+ for (i = 0; i < wrr * hrr; i++) {
+ int r, g, b;
+ unsigned char rc,gc,bc;
+ r = image->comps[0].data[i / wrr * wr + i % wrr];
+ r += (image->comps[0].sgnd ? 1 << (image->comps[0].prec - 1) : 0);
+ rc = (unsigned char) ((r >> adjust)+((r >> (adjust-1))%2));
+
+ g = image->comps[1].data[i / wrr * wr + i % wrr];
+ g += (image->comps[1].sgnd ? 1 << (image->comps[1].prec - 1) : 0);
+ gc = (unsigned char) ((g >> adjust)+((g >> (adjust-1))%2));
+
+ b = image->comps[2].data[i / wrr * wr + i % wrr];
+ b += (image->comps[2].sgnd ? 1 << (image->comps[2].prec - 1) : 0);
+ bc = (unsigned char) ((b >> adjust)+((b >> (adjust-1))%2));
+
+ fprintf(fdest, "%c%c%c", rc, gc, bc);
+ }
+ fclose(fdest);
+ } else {
+ int ncomp=(S2=='g' || S2=='G')?1:image->numcomps;
+ if (image->numcomps>ncomp) {
+ fprintf(stderr,"WARNING -> [PGM files] Only the first component\n");
+ fprintf(stderr," is written to the file\n");
+ }
+ for (compno = 0; compno < ncomp; compno++) {
+ char name[256];
+ if (ncomp > 1) {
+ sprintf(name, "%d.%s", compno, outfile);
+ } else {
+ sprintf(name, "%s", outfile);
+ }
+
+ fdest = fopen(name, "wb");
+ if (!fdest) {
+ fprintf(stderr, "ERROR -> failed to open %s for writing\n", name);
+ return 1;
+ }
+
+ w = int_ceildiv(image->x1 - image->x0, image->comps[compno].dx);
+ /* wr = int_ceildiv(int_ceildivpow2(image->x1 - image->x0,image->factor),image->comps[compno].dx); */
+ wr = image->comps[compno].w;
+ wrr = int_ceildivpow2(image->comps[compno].w, image->comps[compno].factor);
+
+ h = int_ceildiv(image->y1 - image->y0, image->comps[compno].dy);
+ /* hr = int_ceildiv(int_ceildivpow2(image->y1 - image->y0,image->factor), image->comps[compno].dy); */
+ hr = image->comps[compno].h;
+ hrr = int_ceildivpow2(image->comps[compno].h, image->comps[compno].factor);
+
+ max = image->comps[compno].prec > 8 ? 255 : (1 << image->comps[compno].prec) - 1;
+
+ image->comps[compno].x0 = int_ceildivpow2(image->comps[compno].x0 - int_ceildiv(image->x0, image->comps[compno].dx), image->comps[compno].factor);
+ image->comps[compno].y0 = int_ceildivpow2(image->comps[compno].y0 - int_ceildiv(image->y0, image->comps[compno].dy), image->comps[compno].factor);
+
+ fprintf(fdest, "P5\n%d %d\n%d\n", wrr, hrr, max);
+ adjust = image->comps[compno].prec > 8 ? image->comps[compno].prec - 8 : 0;
+
+ for (i = 0; i < wrr * hrr; i++) {
+ int l;
+ unsigned char lc;
+ l = image->comps[compno].data[i / wrr * wr + i % wrr];
+ l += (image->comps[compno].sgnd ? 1 << (image->comps[compno].prec - 1) : 0);
+ lc = (unsigned char) ((l >> adjust)+((l >> (adjust-1))%2));
+ fprintf(fdest, "%c", lc);
+ }
+ fclose(fdest);
+ }
+ }
+
+ return 0;
+}