codec/j2k_to_image.c

   1 /* Copyright (c) 2001 David Janssens
   2  * Copyright (c) 2002-2003 Yannick Verschueren
   3  * Copyright (c) 2002-2003 Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium
   4  *
   5  * All rights reserved.
   6  * Redistribution and use in source and binary forms, with or without
   7  * modification, are permitted provided that the following conditions
   8  * are met:
   9  * 1. Redistributions of source code must retain the above copyright
  10  *    notice, this list of conditions and the following disclaimer.
  11  * 2. Redistributions in binary form must reproduce the above copyright
  12  *    notice, this list of conditions and the following disclaimer in the
  13  *    documentation and/or other materials provided with the distribution.
  14  *
  15  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
  16  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  17  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  18  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  19  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  20  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  21  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  22  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  23  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  24  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  25  * POSSIBILITY OF SUCH DAMAGE.
  26  */
  27
  28
  29
  30 #include <openjpeg.h>
  31 #include <stdio.h>
  32 #include <stdlib.h>
  33 #include <string.h>
  34
  35 int ceildiv(int a, int b)
  36 {
  37   return (a + b - 1) / b;
  38 }
  39
  40 int main(int argc, char **argv)
  41 {
  42   FILE *f;
  43   char *src, *src_name;
  44   char *dest, S1, S2, S3;
  45   int len;
  46
  47   j2k_image_t img;
  48
  49   j2k_cp_t cp;
  50   int w, wr, wrr, h, hr, hrr, max;
  51   int i, image_type = -1, compno, pad, j;
  52   int adjust;
  53   jp2_struct_t *jp2_struct;
  54
  55   if (argc < 3) {
  56     fprintf(stderr,
  57             "usage: %s j2k-file image-file [-reduce n]\n", argv[0]);
  58     return 1;
  59   }
  60
  61   f = fopen(argv[1], "rb");
  62   if (!f) {
  63     fprintf(stderr, "failed to open %s for reading\n", argv[1]);
  64     return 1;
  65   }
  66
  67   dest = argv[2];
  68
  69   cp.reduce_on = 0;
  70   cp.reduce_value = 0;
  71
  72   /* OPTION REDUCE IS ACTIVE */
  73   if (argc == 5) {
  74     if (strcmp(argv[3], "-reduce")) {
  75       fprintf(stderr,
  76               "usage: options " "-reduce n"
  77               " where n is the factor of reduction [%s]\n", argv[3]);
  78       return 1;
  79     }
  80     cp.reduce_on = 1;
  81     sscanf(argv[4], "%d", &cp.reduce_value);
  82   }
  83
  84   while (*dest) {
  85     dest++;
  86   }
  87   dest--;
  88   S3 = *dest;
  89   dest--;
  90   S2 = *dest;
  91   dest--;
  92   S1 = *dest;
  93
  94   if ((S1 == 'p' && S2 == 'g' && S3 == 'x')
  95       || (S1 == 'P' && S2 == 'G' && S3 == 'X')) {
  96     image_type = 0;
  97
  98     dest--;
  99
 100     *dest = '\0';
 101   }
 102
 103   if ((S1 == 'p' && S2 == 'n' && S3 == 'm')
 104       || (S1 == 'P' && S2 == 'N' && S3 == 'M') || (S1 == 'p' && S2 == 'g'
 105                                                    && S3 == 'm')
 106       || (S1 == 'P' && S2 == 'G' && S3 == 'M') || (S1 == 'P' && S2 == 'P'
 107                                                    && S3 == 'M')
 108       || (S1 == 'p' && S2 == 'p' && S3 == 'm')) {
 109     image_type = 1;
 110   }
 111
 112   if ((S1 == 'b' && S2 == 'm' && S3 == 'p')
 113       || (S1 == 'B' && S2 == 'M' && S3 == 'P')) {
 114     image_type = 2;
 115   }
 116
 117   if (image_type == -1) {
 118     fprintf(stderr,
 119             "!! Unrecognized format for infile : %c%c%c [accept only *.pnm, *.pgm, *.ppm, *.pgx or *.bmp] !!\n\n",
 120             S1, S2, S3);
 121     return 1;
 122   }
 123
 124   fseek(f, 0, SEEK_END);
 125   len = ftell(f);
 126   fseek(f, 0, SEEK_SET);
 127   src = (char *) malloc(len);
 128   fread(src, 1, len, f);
 129   fclose(f);
 130
 131   src_name = argv[1];
 132   while (*src_name) {
 133     src_name++;
 134   }
 135   src_name--;
 136   S3 = *src_name;
 137   src_name--;
 138   S2 = *src_name;
 139   src_name--;
 140   S1 = *src_name;
 141
 142   /* J2K format */
 143   if ((S1 == 'j' && S2 == '2' && S3 == 'k')
 144       || (S1 == 'J' && S2 == '2' && S3 == 'K') || (S1 == 'j' && S2 == '2'
 145                                                    && S3 == 'c')
 146       || (S1 == 'J' && S2 == '2' && S3 == 'C')) {
 147     if (!j2k_decode(src, len, &img, &cp)) {
 148       fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 149       return 1;
 150     }
 151   }
 152
 153   /* JP2 format */
 154   else if ((S1 == 'j' && S2 == 'p' && S3 == '2')
 155            || (S1 == 'J' && S2 == 'P' && S3 == '2')) {
 156     jp2_struct = (jp2_struct_t *) malloc(sizeof(jp2_struct_t));
 157
 158     jp2_struct->image = &img;
 159
 160     if (jp2_decode(src, len, jp2_struct, &cp)) {
 161       fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 162       return 1;
 163     }
 164     /* Insert code here if you want to create actions on jp2_struct before deleting it */
 165     free(jp2_struct);
 166   }
 167
 168   /* JPT format */
 169   else if ((S1 == 'j' && S2 == 'p' && S3 == 't')
 170            || (S1 == 'J' && S2 == 'P' && S3 == 'T')) {
 171     if (!j2k_decode_jpt_stream(src, len, &img, &cp)) {
 172       fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 173       return 1;
 174     }
 175   }
 176
 177   /* otherwise : error */
 178   else {
 179     fprintf(stderr,
 180             "j2k_to_image : Unknown format image *.%c%c%c [only *.j2k, *.jp2, *.jpc or *.jpt]!! \n",
 181             S1, S2, S3);
 182     return 1;
 183   }
 184
 185   free(src);
 186   /* ------------------  CREATE OUT IMAGE WITH THE RIGHT FORMAT ----------------------- */
 187
 188   /* ---------------------------- / */
 189   /* /                            / */
 190   /* /  FORMAT : PNM, PGM or PPM  / */
 191   /* /                            / */
 192   /* ---------------------------- / */
 193
 194   switch (image_type) {
 195   case 1:                       /* PNM PGM PPM */
 196     if (img.numcomps == 3 && img.comps[0].dx == img.comps[1].dx
 197         && img.comps[1].dx == img.comps[2].dx
 198         && img.comps[0].dy == img.comps[1].dy
 199         && img.comps[1].dy == img.comps[2].dy
 200         && img.comps[0].prec == img.comps[1].prec
 201         && img.comps[1].prec == img.comps[2].prec) {
 202       f = fopen(argv[2], "wb");
 203       w = ceildiv(img.x1 - img.x0, img.comps[0].dx);
 204       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor),img.comps[0].dx);
 205       wr = img.comps[0].w;
 206       wrr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor);
 207
 208       h = ceildiv(img.y1 - img.y0, img.comps[0].dy);
 209       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy);
 210       hr = img.comps[0].h;
 211       hrr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor);
 212
 213       max = img.comps[0].prec > 8 ? 255 : (1 << img.comps[0].prec) - 1;
 214
 215       img.comps[0].x0 =
 216         int_ceildivpow2(img.comps[0].x0 -
 217                         int_ceildiv(img.x0, img.comps[0].dx),
 218                         img.comps[0].factor);
 219       img.comps[0].y0 =
 220         int_ceildivpow2(img.comps[0].y0 -
 221                         int_ceildiv(img.y0, img.comps[0].dy),
 222                         img.comps[0].factor);
 223
 224
 225       fprintf(f, "P6\n# %d %d %d %d %d\n%d %d\n%d\n",
 226               cp.tcps[cp.tileno[0]].tccps[0].numresolutions, w, h,
 227               img.comps[0].x0, img.comps[0].y0, wrr, hrr, max);
 228       adjust = img.comps[0].prec > 8 ? img.comps[0].prec - 8 : 0;
 229       for (i = 0; i < wrr * hrr; i++) {
 230         char r, g, b;
 231         r = img.comps[0].data[i / wrr * wr + i % wrr];
 232         r += (img.comps[0].sgnd ? 1 << (img.comps[0].prec - 1) : 0);
 233         r = r >> adjust;
 234
 235         g = img.comps[1].data[i / wrr * wr + i % wrr];
 236         g += (img.comps[1].sgnd ? 1 << (img.comps[1].prec - 1) : 0);
 237         g = g >> adjust;
 238
 239         b = img.comps[2].data[i / wrr * wr + i % wrr];
 240         b += (img.comps[2].sgnd ? 1 << (img.comps[2].prec - 1) : 0);
 241         b = b >> adjust;
 242
 243         fprintf(f, "%c%c%c", r, g, b);
 244       }
 245       free(img.comps[0].data);
 246       free(img.comps[1].data);
 247       free(img.comps[2].data);
 248       fclose(f);
 249     } else {
 250       for (compno = 0; compno < img.numcomps; compno++) {
 251         char name[256];
 252         if (img.numcomps > 1) {
 253           sprintf(name, "%d.%s", compno, argv[2]);
 254         } else {
 255           sprintf(name, "%s", argv[2]);
 256         }
 257         f = fopen(name, "wb");
 258         w = ceildiv(img.x1 - img.x0, img.comps[compno].dx);
 259         // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor),img.comps[compno].dx);
 260         wr = img.comps[compno].w;
 261         wrr =
 262           int_ceildivpow2(img.comps[compno].w, img.comps[compno].factor);
 263
 264         h = ceildiv(img.y1 - img.y0, img.comps[compno].dy);
 265         // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[compno].dy);
 266         hr = img.comps[compno].h;
 267         hrr =
 268           int_ceildivpow2(img.comps[compno].h, img.comps[compno].factor);
 269
 270         max =
 271           img.comps[compno].prec >
 272           8 ? 255 : (1 << img.comps[compno].prec) - 1;
 273
 274         img.comps[compno].x0 =
 275           int_ceildivpow2(img.comps[compno].x0 -
 276                           int_ceildiv(img.x0,
 277                                       img.comps[compno].dx),
 278                           img.comps[compno].factor);
 279         img.comps[compno].y0 =
 280           int_ceildivpow2(img.comps[compno].y0 -
 281                           int_ceildiv(img.y0,
 282                                       img.comps[compno].dy),
 283                           img.comps[compno].factor);
 284
 285         fprintf(f, "P5\n# %d %d %d %d %d\n%d %d\n%d\n",
 286                 cp.tcps[cp.tileno[0]].tccps[compno].
 287                 numresolutions, w, h, img.comps[compno].x0,
 288                 img.comps[compno].y0, wrr, hrr, max);
 289         adjust =
 290           img.comps[compno].prec > 8 ? img.comps[compno].prec - 8 : 0;
 291         for (i = 0; i < wrr * hrr; i++) {
 292           char l;
 293           l = img.comps[compno].data[i / wrr * wr + i % wrr];
 294           l += (img.comps[compno].
 295                 sgnd ? 1 << (img.comps[compno].prec - 1) : 0);
 296           l = l >> adjust;
 297           fprintf(f, "%c", l);
 298         }
 299         fclose(f);
 300         free(img.comps[compno].data);
 301       }
 302     }
 303     break;
 304
 305     /* ------------------------ / */
 306     /* /                        / */
 307     /* /     FORMAT : PGX       / */
 308     /* /                        / */
 309     /* /----------------------- / */
 310   case 0:                       /* PGX */
 311     for (compno = 0; compno < img.numcomps; compno++) {
 312       j2k_comp_t *comp = &img.comps[compno];
 313       char name[256];
 314
 315       int nbytes = 0;
 316       //if (img.numcomps > 1)
 317       sprintf(name, "%s-%d.pgx", argv[2], compno);
 318
 319       //else
 320
 321       //sprintf(name, "%s.pgx", argv[2]);
 322
 323       f = fopen(name, "wb");
 324       // w = ceildiv(img.x1 - img.x0, comp->dx);
 325       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), comp->dx);
 326       w = img.comps[compno].w;
 327       wr = int_ceildivpow2(img.comps[compno].w, img.comps[compno].factor);
 328
 329       // h = ceildiv(img.y1 - img.y0, comp->dy);
 330       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), comp->dy);
 331       h = img.comps[compno].h;
 332       hr = int_ceildivpow2(img.comps[compno].h, img.comps[compno].factor);
 333
 334       fprintf(f, "PG ML %c %d %d %d\n", comp->sgnd ? '-' : '+',
 335               comp->prec, wr, hr);
 336
 337       if (comp->prec <= 8)
 338         nbytes = 1;
 339
 340       else if (comp->prec <= 16)
 341         nbytes = 2;
 342
 343       else
 344         nbytes = 4;
 345       for (i = 0; i < wr * hr; i++) {
 346         int v = img.comps[compno].data[i / wr * w + i % wr];
 347
 348         for (j = nbytes - 1; j >= 0; j--) {
 349
 350           char byte = (char) (v >> (j * 8));
 351
 352           fwrite(&byte, 1, 1, f);
 353
 354         }
 355       }
 356       free(img.comps[compno].data);
 357       fclose(f);
 358     }
 359     break;
 360
 361     /* ------------------------ / */
 362     /* /                        / */
 363     /* /     FORMAT : BMP       / */
 364     /* /                        / */
 365     /* /----------------------- / */
 366
 367   case 2:                       /* BMP */
 368     if (img.numcomps == 3 && img.comps[0].dx == img.comps[1].dx
 369         && img.comps[1].dx == img.comps[2].dx
 370         && img.comps[0].dy == img.comps[1].dy
 371         && img.comps[1].dy == img.comps[2].dy
 372         && img.comps[0].prec == img.comps[1].prec
 373         && img.comps[1].prec == img.comps[2].prec) {
 374       /* -->> -->> -->> -->>
 375
 376          24 bits color
 377
 378          <<-- <<-- <<-- <<-- */
 379
 380       f = fopen(argv[2], "wb");
 381       // w = ceildiv(img.x1 - img.x0, img.comps[0].dx);
 382       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), img.comps[0].dx);
 383       w = img.comps[0].w;
 384       wr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor);
 385
 386       // h = ceildiv(img.y1 - img.y0, img.comps[0].dy);
 387       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy);
 388       h = img.comps[0].h;
 389       hr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor);
 390
 391       fprintf(f, "BM");
 392
 393       /* FILE HEADER */
 394       /* ------------- */
 395       fprintf(f, "%c%c%c%c",
 396               (unsigned char) (hr * wr * 3 + 3 * hr * (wr % 2) +
 397                                54) & 0xff,
 398               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 399                                >> 8) & 0xff,
 400               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 401                                >> 16) & 0xff,
 402               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 403                                >> 24) & 0xff);
 404       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 405               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 406       fprintf(f, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff,
 407               ((54) >> 16) & 0xff, ((54) >> 24) & 0xff);
 408
 409       /* INFO HEADER   */
 410       /* ------------- */
 411       fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,
 412               ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
 413       fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
 414               (unsigned char) ((wr) >> 8) & 0xff,
 415               (unsigned char) ((wr) >> 16) & 0xff,
 416               (unsigned char) ((wr) >> 24) & 0xff);
 417       fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
 418               (unsigned char) ((hr) >> 8) & 0xff,
 419               (unsigned char) ((hr) >> 16) & 0xff,
 420               (unsigned char) ((hr) >> 24) & 0xff);
 421       fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
 422       fprintf(f, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff);
 423       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 424               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 425       fprintf(f, "%c%c%c%c",
 426               (unsigned char) (3 * hr * wr +
 427                                3 * hr * (wr % 2)) & 0xff,
 428               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 429                                8) & 0xff,
 430               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 431                                16) & 0xff,
 432               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 433                                24) & 0xff);
 434       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 435               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 436       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 437               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 438       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 439               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 440       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 441               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 442
 443       for (i = 0; i < wr * hr; i++) {
 444         unsigned char R, G, B;
 445         /* a modifier */
 446         // R = img.comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 447         R = img.comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
 448         // G = img.comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 449         G = img.comps[1].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
 450         // B = img.comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 451         B = img.comps[2].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
 452         fprintf(f, "%c%c%c", B, G, R);
 453
 454         if ((i + 1) % wr == 0) {
 455           for (pad = (3 * wr) % 4 ? 4 - (3 * wr) % 4 : 0; pad > 0; pad--)       /* ADD */
 456             fprintf(f, "%c", 0);
 457         }
 458       }
 459       fclose(f);
 460       free(img.comps[1].data);
 461       free(img.comps[2].data);
 462     } else {                    /* Gray-scale */
 463
 464       /* -->> -->> -->> -->>
 465
 466          8 bits non code (Gray scale)
 467
 468          <<-- <<-- <<-- <<-- */
 469       f = fopen(argv[2], "wb");
 470       // w = ceildiv(img.x1 - img.x0, img.comps[0].dx);
 471       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), img.comps[0].dx);
 472       w = img.comps[0].w;
 473       wr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor);
 474
 475       // h = ceildiv(img.y1 - img.y0, img.comps[0].dy);
 476       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy);
 477       h = img.comps[0].h;
 478       hr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor);
 479
 480       fprintf(f, "BM");
 481
 482       /* FILE HEADER */
 483       /* ------------- */
 484       fprintf(f, "%c%c%c%c",
 485               (unsigned char) (hr * wr + 54 + 1024 +
 486                                hr * (wr % 2)) & 0xff,
 487               (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2))
 488                                >> 8) & 0xff,
 489               (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2))
 490                                >> 16) & 0xff,
 491               (unsigned char) ((hr * wr + 54 + 1024 + wr * (wr % 2))
 492                                >> 24) & 0xff);
 493       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 494               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 495       fprintf(f, "%c%c%c%c", (54 + 1024) & 0xff,
 496               ((54 + 1024) >> 8) & 0xff, ((54 + 1024) >> 16) & 0xff,
 497               ((54 + 1024) >> 24) & 0xff);
 498
 499       /* INFO HEADER */
 500       /* ------------- */
 501       fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,
 502               ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
 503       fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
 504               (unsigned char) ((wr) >> 8) & 0xff,
 505               (unsigned char) ((wr) >> 16) & 0xff,
 506               (unsigned char) ((wr) >> 24) & 0xff);
 507       fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
 508               (unsigned char) ((hr) >> 8) & 0xff,
 509               (unsigned char) ((hr) >> 16) & 0xff,
 510               (unsigned char) ((hr) >> 24) & 0xff);
 511       fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
 512       fprintf(f, "%c%c", (8) & 0xff, ((8) >> 8) & 0xff);
 513       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 514               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 515       fprintf(f, "%c%c%c%c",
 516               (unsigned char) (hr * wr + hr * (wr % 2)) & 0xff,
 517               (unsigned char) ((hr * wr + hr * (wr % 2)) >> 8) &
 518               0xff,
 519               (unsigned char) ((hr * wr + hr * (wr % 2)) >> 16) &
 520               0xff,
 521               (unsigned char) ((hr * wr + hr * (wr % 2)) >> 24) & 0xff);
 522       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 523               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 524       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 525               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 526       fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
 527               ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
 528       fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
 529               ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
 530     }
 531
 532     for (i = 0; i < 256; i++) {
 533       fprintf(f, "%c%c%c%c", i, i, i, 0);
 534     }
 535
 536     for (i = 0; i < wr * hr; i++) {
 537       /* a modifier !! */
 538       // fprintf(f, "%c", img.comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]);
 539       fprintf(f, "%c",
 540               img.comps[0].data[w * hr - ((i) / (wr) + 1) * w +
 541                                 (i) % (wr)]);
 542       /*if (((i + 1) % w == 0 && w % 2))
 543          fprintf(f, "%c", 0); */
 544       if ((i + 1) % wr == 0) {
 545         for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--)     /* ADD */
 546           fprintf(f, "%c", 0);
 547       }
 548     }
 549     fclose(f);
 550     free(img.comps[0].data);
 551     break;
 552   default:
 553     break;
 554   }
 555
 556   return 0;
 557 }