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",
  58             argv[0]);
  59     return 1;
  60   }
  61
  62   f = fopen(argv[1], "rb");
  63   if (!f) {
  64     fprintf(stderr, "failed to open %s for reading\n", argv[1]);
  65     return 1;
  66   }
  67
  68   dest = argv[2];
  69
  70   cp.reduce_on = 0;
  71   cp.reduce_value = 0;
  72
  73   /* OPTION REDUCE IS ACTIVE */
  74   if (argc == 5) {
  75     if (strcmp(argv[3], "-reduce")) {
  76       fprintf(stderr,
  77               "usage: options " "-reduce n"
  78               " where n is the factor of reduction [%s]\n", argv[3]);
  79       return 1;
  80     }
  81     cp.reduce_on = 1;
  82     sscanf(argv[4], "%d", &cp.reduce_value);
  83   }
  84
  85   while (*dest) {
  86     dest++;
  87   }
  88   dest--;
  89   S3 = *dest;
  90   dest--;
  91   S2 = *dest;
  92   dest--;
  93   S1 = *dest;
  94
  95   if ((S1 == 'p' && S2 == 'g' && S3 == 'x')
  96       || (S1 == 'P' && S2 == 'G' && S3 == 'X')) {
  97     image_type = 0;
  98
  99     dest--;
 100
 101     *dest = '\0';
 102   }
 103
 104   if ((S1 == 'p' && S2 == 'n' && S3 == 'm')
 105       || (S1 == 'P' && S2 == 'N' && S3 == 'M') || (S1 == 'p' && S2 == 'g'
 106                                                    && S3 == 'm')
 107       || (S1 == 'P' && S2 == 'G' && S3 == 'M') || (S1 == 'P' && S2 == 'P'
 108                                                    && S3 == 'M')
 109       || (S1 == 'p' && S2 == 'p' && S3 == 'm')) {
 110     image_type = 1;
 111   }
 112
 113   if ((S1 == 'b' && S2 == 'm' && S3 == 'p')
 114       || (S1 == 'B' && S2 == 'M' && S3 == 'P')) {
 115     image_type = 2;
 116   }
 117
 118   if (image_type == -1) {
 119     fprintf(stderr,
 120             "!! Unrecognized format for infile : %c%c%c [accept only *.pnm, *.pgm, *.ppm, *.pgx or *.bmp] !!\n\n",
 121             S1, S2, S3);
 122     return 1;
 123   }
 124
 125   fseek(f, 0, SEEK_END);
 126   len = ftell(f);
 127   fseek(f, 0, SEEK_SET);
 128   src = (char *) malloc(len);
 129   fread(src, 1, len, f);
 130   fclose(f);
 131
 132   src_name = argv[1];
 133   while (*src_name) {
 134     src_name++;
 135   }
 136   src_name--;
 137   S3 = *src_name;
 138   src_name--;
 139   S2 = *src_name;
 140   src_name--;
 141   S1 = *src_name;
 142
 143   /* J2K format */
 144   if ((S1 == 'j' && S2 == '2' && S3 == 'k') || (S1 == 'J' && S2 == '2' && S3 == 'K')) {
 145
 146     if (!j2k_decode(src, len, &img, &cp)) {
 147       fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 148       return 1;
 149     }
 150   }
 151   /* JP2 format */
 152   else if ((S1 == 'j' && S2 == 'p' && S3 == '2') || (S1 == 'J' && S2 == 'P' && S3 == '2')) {
 153     jp2_struct = (jp2_struct_t *) malloc(sizeof(jp2_struct_t));
 154
 155     jp2_struct->image=&img;
 156
 157     if (jp2_decode(src,len,jp2_struct,&cp)) {
 158       fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 159       return 1;
 160     }
 161     /* Insert code here if you want to create actions on jp2_struct before deleting it */
 162     free(jp2_struct);
 163   }
 164   /* JPT format */
 165   else if ((S1 == 'j' && S2 == 'p' && S3 == 't') || (S1 == 'J' && S2 == 'P' && S3 == 'T')){
 166       if (!j2k_decode_jpt_stream(src, len, &img, &cp)) {
 167         fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 168         return 1;
 169       }
 170     else {
 171       fprintf(stderr,
 172               "j2k_to_image : Unknown format image *.%c%c%c [only *.j2k, *.jp2 or *.jpt]!! \n",
 173               S1, S2, S3);
 174       return 1;
 175     }
 176   }
 177
 178   free(src);
 179   /* ------------------  CREATE OUT IMAGE WITH THE RIGHT FORMAT ----------------------- */
 180
 181   /* ---------------------------- / */
 182   /* /                            / */
 183   /* /  FORMAT : PNM, PGM or PPM  / */
 184   /* /                            / */
 185   /* ---------------------------- / */
 186
 187   switch (image_type) {
 188   case 1:                       /* PNM PGM PPM */
 189     if (img.numcomps == 3 && img.comps[0].dx == img.comps[1].dx
 190         && img.comps[1].dx == img.comps[2].dx
 191         && img.comps[0].dy == img.comps[1].dy
 192         && img.comps[1].dy == img.comps[2].dy
 193         && img.comps[0].prec == img.comps[1].prec
 194         && img.comps[1].prec == img.comps[2].prec) {
 195       f = fopen(argv[2], "wb");
 196       w = ceildiv(img.x1 - img.x0, img.comps[0].dx);
 197       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor),img.comps[0].dx);
 198       wr = img.comps[0].w;
 199       wrr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor);
 200
 201       h = ceildiv(img.y1 - img.y0, img.comps[0].dy);
 202       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy);
 203       hr = img.comps[0].h;
 204       hrr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor);
 205
 206       max = img.comps[0].prec > 8 ? 255 : (1 << img.comps[0].prec) - 1;
 207
 208       img.comps[0].x0 =
 209         int_ceildivpow2(img.comps[0].x0 -
 210                         int_ceildiv(img.x0, img.comps[0].dx),
 211                         img.comps[0].factor);
 212       img.comps[0].y0 =
 213         int_ceildivpow2(img.comps[0].y0 -
 214                         int_ceildiv(img.y0, img.comps[0].dy),
 215                         img.comps[0].factor);
 216
 217
 218       fprintf(f, "P6\n# %d %d %d %d %d\n%d %d\n%d\n",
 219               cp.tcps[cp.tileno[0]].tccps[0].numresolutions, w, h,
 220               img.comps[0].x0, img.comps[0].y0, wrr, hrr, max);
 221       adjust = img.comps[0].prec > 8 ? img.comps[0].prec - 8 : 0;
 222       for (i = 0; i < wrr * hrr; i++) {
 223         char r, g, b;
 224         r = img.comps[0].data[i / wrr * wr + i % wrr];
 225         r += (img.comps[0].sgnd ? 1 << (img.comps[0].prec - 1) : 0);
 226         r = r >> adjust;
 227
 228         g = img.comps[1].data[i / wrr * wr + i % wrr];
 229         g += (img.comps[1].sgnd ? 1 << (img.comps[1].prec - 1) : 0);
 230         g = g >> adjust;
 231
 232         b = img.comps[2].data[i / wrr * wr + i % wrr];
 233         b += (img.comps[2].sgnd ? 1 << (img.comps[2].prec - 1) : 0);
 234         b = b >> adjust;
 235
 236         fprintf(f, "%c%c%c", r, g, b);
 237       }
 238       free(img.comps[0].data);
 239       free(img.comps[1].data);
 240       free(img.comps[2].data);
 241       fclose(f);
 242     } else {
 243       for (compno = 0; compno < img.numcomps; compno++) {
 244         char name[256];
 245         if (img.numcomps > 1) {
 246           sprintf(name, "%d.%s", compno, argv[2]);
 247         } else {
 248           sprintf(name, "%s", argv[2]);
 249         }
 250         f = fopen(name, "wb");
 251         w = ceildiv(img.x1 - img.x0, img.comps[compno].dx);
 252         // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor),img.comps[compno].dx);
 253         wr = img.comps[compno].w;
 254         wrr =
 255           int_ceildivpow2(img.comps[compno].w, img.comps[compno].factor);
 256
 257         h = ceildiv(img.y1 - img.y0, img.comps[compno].dy);
 258         // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[compno].dy);
 259         hr = img.comps[compno].h;
 260         hrr =
 261           int_ceildivpow2(img.comps[compno].h, img.comps[compno].factor);
 262
 263         max =
 264           img.comps[compno].prec >
 265           8 ? 255 : (1 << img.comps[compno].prec) - 1;
 266
 267         img.comps[compno].x0 =
 268           int_ceildivpow2(img.comps[compno].x0 -
 269                           int_ceildiv(img.x0,
 270                                       img.comps[compno].dx),
 271                           img.comps[compno].factor);
 272         img.comps[compno].y0 =
 273           int_ceildivpow2(img.comps[compno].y0 -
 274                           int_ceildiv(img.y0,
 275                                       img.comps[compno].dy),
 276                           img.comps[compno].factor);
 277
 278         fprintf(f, "P5\n# %d %d %d %d %d\n%d %d\n%d\n",
 279                 cp.tcps[cp.tileno[0]].tccps[compno].
 280                 numresolutions, w, h, img.comps[compno].x0,
 281                 img.comps[compno].y0, wrr, hrr, max);
 282         adjust =
 283           img.comps[compno].prec > 8 ? img.comps[compno].prec - 8 : 0;
 284         for (i = 0; i < wrr * hrr; i++) {
 285           char l;
 286           l = img.comps[compno].data[i / wrr * wr + i % wrr];
 287           l += (img.comps[compno].
 288                 sgnd ? 1 << (img.comps[compno].prec - 1) : 0);
 289           l = l >> adjust;
 290           fprintf(f, "%c", l);
 291         }
 292         fclose(f);
 293         free(img.comps[compno].data);
 294       }
 295     }
 296     break;
 297
 298     /* ------------------------ / */
 299     /* /                        / */
 300     /* /     FORMAT : PGX       / */
 301     /* /                        / */
 302     /* /----------------------- / */
 303   case 0:                       /* PGX */
 304     for (compno = 0; compno < img.numcomps; compno++) {
 305       j2k_comp_t *comp = &img.comps[compno];
 306       char name[256];
 307
 308       int nbytes=0;
 309       //if (img.numcomps > 1)
 310         sprintf(name, "%s-%d.pgx", argv[2], compno);
 311
 312       //else
 313
 314         //sprintf(name, "%s.pgx", argv[2]);
 315
 316       f = fopen(name, "wb");
 317       // w = ceildiv(img.x1 - img.x0, comp->dx);
 318       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), comp->dx);
 319       w = img.comps[compno].w;
 320       wr = int_ceildivpow2(img.comps[compno].w,
 321                            img.comps[compno].factor);
 322
 323       // h = ceildiv(img.y1 - img.y0, comp->dy);
 324       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), comp->dy);
 325       h = img.comps[compno].h;
 326       hr = int_ceildivpow2(img.comps[compno].h,
 327                            img.comps[compno].factor);
 328
 329       fprintf(f, "PG ML %c %d %d %d\n", comp->sgnd ? '-' : '+',
 330               comp->prec, wr, hr);
 331
 332       if (comp->prec <= 8) nbytes=1;
 333
 334       else if (comp->prec <= 16) nbytes=2;
 335
 336       else nbytes=4;
 337       for (i = 0; i < wr * hr; i++) {
 338         int v = img.comps[compno].data[i / wr * w + i % wr];
 339
 340         for (j=nbytes-1 ; j>=0 ; j--) {
 341
 342           char byte=(char)(v>>(j*8));
 343
 344           fwrite(&byte, 1, 1, f);
 345
 346         }
 347       }
 348       free(img.comps[compno].data);
 349       fclose(f);
 350     }
 351     break;
 352
 353     /* ------------------------ / */
 354     /* /                        / */
 355     /* /     FORMAT : BMP       / */
 356     /* /                        / */
 357     /* /----------------------- / */
 358
 359   case 2:                       /* BMP */
 360     if (img.numcomps == 3 && img.comps[0].dx == img.comps[1].dx
 361         && img.comps[1].dx == img.comps[2].dx
 362         && img.comps[0].dy == img.comps[1].dy
 363         && img.comps[1].dy == img.comps[2].dy
 364         && img.comps[0].prec == img.comps[1].prec
 365         && img.comps[1].prec == img.comps[2].prec) {
 366       /* -->> -->> -->> -->>
 367
 368          24 bits color
 369
 370          <<-- <<-- <<-- <<-- */
 371
 372       f = fopen(argv[2], "wb");
 373       // w = ceildiv(img.x1 - img.x0, img.comps[0].dx);
 374       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), img.comps[0].dx);
 375       w = img.comps[0].w;
 376       wr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor);
 377
 378       // h = ceildiv(img.y1 - img.y0, img.comps[0].dy);
 379       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy);
 380       h = img.comps[0].h;
 381       hr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor);
 382
 383       fprintf(f, "BM");
 384
 385       /* FILE HEADER */
 386       /* ------------- */
 387       fprintf(f, "%c%c%c%c",
 388               (unsigned char) (hr * wr * 3 + 3 * hr * (wr % 2) +
 389                                54) & 0xff,
 390               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 391                                >> 8) & 0xff,
 392               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 393                                >> 16) & 0xff,
 394               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 395                                >> 24) & 0xff);
 396       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 397               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 398       fprintf(f, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff,
 399               ((54) >> 16) & 0xff, ((54) >> 24) & 0xff);
 400
 401       /* INFO HEADER   */
 402       /* ------------- */
 403       fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,
 404               ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
 405       fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
 406               (unsigned char) ((wr) >> 8) & 0xff,
 407               (unsigned char) ((wr) >> 16) & 0xff,
 408               (unsigned char) ((wr) >> 24) & 0xff);
 409       fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
 410               (unsigned char) ((hr) >> 8) & 0xff,
 411               (unsigned char) ((hr) >> 16) & 0xff,
 412               (unsigned char) ((hr) >> 24) & 0xff);
 413       fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
 414       fprintf(f, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff);
 415       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 416               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 417       fprintf(f, "%c%c%c%c",
 418               (unsigned char) (3 * hr * wr +
 419                                3 * hr * (wr % 2)) & 0xff,
 420               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 421                                8) & 0xff,
 422               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 423                                16) & 0xff,
 424               (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 425                                24) & 0xff);
 426       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 427               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 428       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 429               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 430       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 431               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 432       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 433               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 434
 435       for (i = 0; i < wr * hr; i++) {
 436         unsigned char R, G, B;
 437         /* a modifier */
 438         // R = img.comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 439         R = img.comps[0].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
 440         // G = img.comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 441         G = img.comps[1].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
 442         // B = img.comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 443         B = img.comps[2].data[w * hr - ((i) / (wr) + 1) * w + (i) % (wr)];
 444         fprintf(f, "%c%c%c", B, G, R);
 445
 446         if ((i + 1) % wr == 0) {
 447           for (pad = (3 * wr) % 4 ? 4 - (3 * wr) % 4 : 0; pad > 0; pad--)       /* ADD */
 448             fprintf(f, "%c", 0);
 449         }
 450       }
 451       fclose(f);
 452       free(img.comps[1].data);
 453       free(img.comps[2].data);
 454     } else {                    /* Gray-scale */
 455
 456       /* -->> -->> -->> -->>
 457
 458          8 bits non code (Gray scale)
 459
 460          <<-- <<-- <<-- <<-- */
 461       f = fopen(argv[2], "wb");
 462       // w = ceildiv(img.x1 - img.x0, img.comps[0].dx);
 463       // wr = ceildiv(int_ceildivpow2(img.x1 - img.x0,img.factor), img.comps[0].dx);
 464       w = img.comps[0].w;
 465       wr = int_ceildivpow2(img.comps[0].w, img.comps[0].factor);
 466
 467       // h = ceildiv(img.y1 - img.y0, img.comps[0].dy);
 468       // hr = ceildiv(int_ceildivpow2(img.y1 - img.y0,img.factor), img.comps[0].dy);
 469       h = img.comps[0].h;
 470       hr = int_ceildivpow2(img.comps[0].h, img.comps[0].factor);
 471
 472       fprintf(f, "BM");
 473
 474       /* FILE HEADER */
 475       /* ------------- */
 476       fprintf(f, "%c%c%c%c",
 477               (unsigned char) (hr * wr + 54 + 1024 +
 478                                hr * (wr % 2)) & 0xff,
 479               (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2))
 480                                >> 8) & 0xff,
 481               (unsigned char) ((hr * wr + 54 + 1024 + hr * (wr % 2))
 482                                >> 16) & 0xff,
 483               (unsigned char) ((hr * wr + 54 + 1024 + wr * (wr % 2))
 484                                >> 24) & 0xff);
 485       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 486               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 487       fprintf(f, "%c%c%c%c", (54 + 1024) & 0xff,
 488               ((54 + 1024) >> 8) & 0xff, ((54 + 1024) >> 16) & 0xff,
 489               ((54 + 1024) >> 24) & 0xff);
 490
 491       /* INFO HEADER */
 492       /* ------------- */
 493       fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,
 494               ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
 495       fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
 496               (unsigned char) ((wr) >> 8) & 0xff,
 497               (unsigned char) ((wr) >> 16) & 0xff,
 498               (unsigned char) ((wr) >> 24) & 0xff);
 499       fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
 500               (unsigned char) ((hr) >> 8) & 0xff,
 501               (unsigned char) ((hr) >> 16) & 0xff,
 502               (unsigned char) ((hr) >> 24) & 0xff);
 503       fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
 504       fprintf(f, "%c%c", (8) & 0xff, ((8) >> 8) & 0xff);
 505       fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 506               ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 507       fprintf(f, "%c%c%c%c",
 508               (unsigned char) (hr * wr + hr * (wr % 2)) & 0xff,
 509               (unsigned char) ((hr * wr + hr * (wr % 2)) >> 8) &
 510               0xff,
 511               (unsigned char) ((hr * wr + hr * (wr % 2)) >> 16) &
 512               0xff,
 513               (unsigned char) ((hr * wr + hr * (wr % 2)) >> 24) & 0xff);
 514       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 515               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 516       fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 517               ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 518       fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
 519               ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
 520       fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
 521               ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
 522     }
 523
 524     for (i = 0; i < 256; i++) {
 525       fprintf(f, "%c%c%c%c", i, i, i, 0);
 526     }
 527
 528     for (i = 0; i < wr * hr; i++) {
 529       /* a modifier !! */
 530       // fprintf(f, "%c", img.comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]);
 531       fprintf(f, "%c",
 532               img.comps[0].data[w * hr - ((i) / (wr) + 1) * w +
 533                                  (i) % (wr)]);
 534       /*if (((i + 1) % w == 0 && w % 2))
 535          fprintf(f, "%c", 0); */
 536       if ((i + 1) % wr == 0) {
 537         for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--)     /* ADD */
 538           fprintf(f, "%c", 0);
 539       }
 540     }
 541     fclose(f);
 542     free(img.comps[0].data);
 543     break;
 544   default:
 545     break;
 546   }
 547
 548   return 0;
 549 }