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