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 //#include <unistd.h>
  35
  36 int ceildiv(int a, int b)
  37 {
  38     return (a + b - 1) / b;
  39 }
  40
  41 int main(int argc, char **argv)
  42 {
  43     FILE *f;
  44     char *src, *src_name;
  45     char *dest, S1, S2, S3;
  46     int len;
  47     j2k_image_t *img;
  48     j2k_cp_t *cp;
  49     j2k_option_t option;
  50     int w, wr, wrr, h, hr, hrr, max;
  51     int i, image_type = -1, compno, pad;
  52     int adjust;
  53
  54     if (argc < 3) {
  55         fprintf(stderr,
  56                 "usage: %s j2k-file image-file -reduce n (<- optional)\n",
  57                 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     option.reduce_on = 0;
  70     option.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         option.reduce_on = 1;
  81         sscanf(argv[4], "%d", &option.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
  99     if ((S1 == 'p' && S2 == 'n' && S3 == 'm')
 100         || (S1 == 'P' && S2 == 'N' && S3 == 'M') || (S1 == 'p' && S2 == 'g'
 101                                                      && S3 == 'm')
 102         || (S1 == 'P' && S2 == 'G' && S3 == 'M') || (S1 == 'P' && S2 == 'P'
 103                                                      && S3 == 'M')
 104         || (S1 == 'p' && S2 == 'p' && S3 == 'm')) {
 105         image_type = 1;
 106     }
 107
 108     if ((S1 == 'b' && S2 == 'm' && S3 == 'p')
 109         || (S1 == 'B' && S2 == 'M' && S3 == 'P')) {
 110         image_type = 2;
 111     }
 112
 113     if (image_type == -1) {
 114         fprintf(stderr,
 115                 "!! Unrecognized format for infile : %c%c%c [accept only *.pnm, *.pgm, *.ppm, *.pgx or *.bmp] !!\n\n",
 116                 S1, S2, S3);
 117         return 1;
 118     }
 119
 120     fseek(f, 0, SEEK_END);
 121     len = ftell(f);
 122     fseek(f, 0, SEEK_SET);
 123     src = (char *) malloc(len);
 124     fread(src, 1, len, f);
 125     fclose(f);
 126
 127     src_name = argv[1];
 128     while (*src_name) {
 129         src_name++;
 130     }
 131     src_name--;
 132     S3 = *src_name;
 133     src_name--;
 134     S2 = *src_name;
 135     src_name--;
 136     S1 = *src_name;
 137
 138     if (S1 == 'j' && S2 == '2' && S3 == 'k') {
 139         if (!j2k_decode(src, len, &img, &cp, option)) {
 140             fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 141             return 1;
 142         }
 143     } else {
 144         if (S1 == 'j' && S2 == 'p' && S3 == 't') {
 145             if (!j2k_decode_jpt_stream(src, len, &img, &cp)) {
 146                 fprintf(stderr, "j2k_to_image: failed to decode image!\n");
 147                 return 1;
 148             }
 149         } else {
 150             fprintf(stderr,
 151                     "j2k_to_image : Unknown format image *.%c%c%c [only *.j2k or *.jpt]!! \n",
 152                     S1, S2, S3);
 153             return 1;
 154         }
 155     }
 156
 157     free(src);
 158     /* ------------------  CREATE OUT IMAGE WITH THE RIGHT FORMAT ----------------------- */
 159
 160     /* ---------------------------- / */
 161     /* /                            / */
 162     /* /  FORMAT : PNM, PGM or PPM  / */
 163     /* /                            / */
 164     /* ---------------------------- / */
 165
 166     switch (image_type) {
 167     case 1:                     /* PNM PGM PPM */
 168         if (img->numcomps == 3 && img->comps[0].dx == img->comps[1].dx
 169             && img->comps[1].dx == img->comps[2].dx
 170             && img->comps[0].dy == img->comps[1].dy
 171             && img->comps[1].dy == img->comps[2].dy
 172             && img->comps[0].prec == img->comps[1].prec
 173             && img->comps[1].prec == img->comps[2].prec) {
 174             f = fopen(argv[2], "wb");
 175             w = ceildiv(img->x1 - img->x0, img->comps[0].dx);
 176             // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor),img->comps[0].dx);
 177             wr = img->comps[0].w;
 178             wrr = int_ceildivpow2(img->comps[0].w, img->comps[0].factor);
 179
 180             h = ceildiv(img->y1 - img->y0, img->comps[0].dy);
 181             // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), img->comps[0].dy);
 182             hr = img->comps[0].h;
 183             hrr = int_ceildivpow2(img->comps[0].h, img->comps[0].factor);
 184
 185             max =
 186                 img->comps[0].prec >
 187                 8 ? 255 : (1 << img->comps[0].prec) - 1;
 188
 189             img->comps[0].x0 =
 190                 int_ceildivpow2(img->comps[0].x0 -
 191                                 int_ceildiv(img->x0, img->comps[0].dx),
 192                                 img->comps[0].factor);
 193             img->comps[0].y0 =
 194                 int_ceildivpow2(img->comps[0].y0 -
 195                                 int_ceildiv(img->y0, img->comps[0].dy),
 196                                 img->comps[0].factor);
 197
 198
 199             fprintf(f, "P6\n# %d %d %d %d %d\n%d %d\n%d\n",
 200                     cp->tcps[cp->tileno[0]].tccps[0].numresolutions, w, h,
 201                     img->comps[0].x0, img->comps[0].y0, wrr, hrr, max);
 202             adjust = img->comps[0].prec > 8 ? img->comps[0].prec - 8 : 0;
 203             for (i = 0; i < wrr * hrr; i++) {
 204                 char r, g, b;
 205                 r = img->comps[0].data[i / wrr * wr + i % wrr];
 206                 r += (img->comps[0].
 207                       sgnd ? 1 << (img->comps[0].prec - 1) : 0);
 208                 r = r >> adjust;
 209
 210                 g = img->comps[1].data[i / wrr * wr + i % wrr];
 211                 g += (img->comps[1].
 212                       sgnd ? 1 << (img->comps[1].prec - 1) : 0);
 213                 g = g >> adjust;
 214
 215                 b = img->comps[2].data[i / wrr * wr + i % wrr];
 216                 b += (img->comps[2].
 217                       sgnd ? 1 << (img->comps[2].prec - 1) : 0);
 218                 b = b >> adjust;
 219
 220                 fprintf(f, "%c%c%c", r, g, b);
 221             }
 222             fclose(f);
 223         } else {
 224             for (compno = 0; compno < img->numcomps; compno++) {
 225                 char name[256];
 226                 if (img->numcomps > 1) {
 227                     sprintf(name, "%d.%s", compno, argv[2]);
 228                 } else {
 229                     sprintf(name, "%s", argv[2]);
 230                 }
 231                 f = fopen(name, "wb");
 232                 w = ceildiv(img->x1 - img->x0, img->comps[compno].dx);
 233                 // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor),img->comps[compno].dx);
 234                 wr = img->comps[compno].w;
 235                 wrr =
 236                     int_ceildivpow2(img->comps[compno].w,
 237                                     img->comps[compno].factor);
 238
 239                 h = ceildiv(img->y1 - img->y0, img->comps[compno].dy);
 240                 // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), img->comps[compno].dy);
 241                 hr = img->comps[compno].h;
 242                 hrr =
 243                     int_ceildivpow2(img->comps[compno].h,
 244                                     img->comps[compno].factor);
 245
 246                 max =
 247                     img->comps[compno].prec >
 248                     8 ? 255 : (1 << img->comps[compno].prec) - 1;
 249
 250                 img->comps[compno].x0 =
 251                     int_ceildivpow2(img->comps[compno].x0 -
 252                                     int_ceildiv(img->x0,
 253                                                 img->comps[compno].dx),
 254                                     img->comps[compno].factor);
 255                 img->comps[compno].y0 =
 256                     int_ceildivpow2(img->comps[compno].y0 -
 257                                     int_ceildiv(img->y0,
 258                                                 img->comps[compno].dy),
 259                                     img->comps[compno].factor);
 260
 261                 fprintf(f, "P5\n# %d %d %d %d %d\n%d %d\n%d\n",
 262                         cp->tcps[cp->tileno[0]].tccps[compno].
 263                         numresolutions, w, h, img->comps[compno].x0,
 264                         img->comps[compno].y0, wrr, hrr, max);
 265                 adjust =
 266                     img->comps[compno].prec >
 267                     8 ? img->comps[compno].prec - 8 : 0;
 268                 for (i = 0; i < wrr * hrr; i++) {
 269                     char l;
 270                     l = img->comps[compno].data[i / wrr * wr + i % wrr];
 271                     l += (img->comps[compno].
 272                           sgnd ? 1 << (img->comps[compno].prec - 1) : 0);
 273                     l = l >> adjust;
 274                     fprintf(f, "%c", l);
 275                 }
 276                 fclose(f);
 277             }
 278         }
 279         break;
 280
 281         /* ------------------------ / */
 282         /* /                        / */
 283         /* /     FORMAT : PGX       / */
 284         /* /                        / */
 285         /* /----------------------- / */
 286     case 0:                     /* PGX */
 287         for (compno = 0; compno < img->numcomps; compno++) {
 288             j2k_comp_t *comp = &img->comps[compno];
 289             char name[256];
 290             if (img->numcomps > 1)
 291                 sprintf(name, "%d_%s", compno, argv[2]);
 292             else
 293                 sprintf(name, "%s", argv[2]);
 294
 295             f = fopen(name, "wb");
 296             // w = ceildiv(img->x1 - img->x0, comp->dx);
 297             // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor), comp->dx);
 298             w = img->comps[compno].w;
 299             wr = int_ceildivpow2(img->comps[compno].w,
 300                                  img->comps[compno].factor);
 301
 302             // h = ceildiv(img->y1 - img->y0, comp->dy);
 303             // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), comp->dy);
 304             h = img->comps[compno].h;
 305             hr = int_ceildivpow2(img->comps[compno].h,
 306                                  img->comps[compno].factor);
 307
 308             fprintf(f, "PG LM %c %d %d %d\n", comp->sgnd ? '-' : '+',
 309                     comp->prec, wr, hr);
 310             for (i = 0; i < wr * hr; i++) {
 311                 int v = img->comps[compno].data[i / wr * w + i % wr];
 312                 if (comp->prec <= 8) {
 313                     char c = (char) v;
 314                     fwrite(&c, 1, 1, f);
 315                 } else if (comp->prec <= 16) {
 316                     short s = (short) v;
 317                     fwrite(&s, 2, 1, f);
 318                 } else {
 319                     fwrite(&v, 4, 1, f);
 320                 }
 321             }
 322             fclose(f);
 323         }
 324         break;
 325
 326         /* ------------------------ / */
 327         /* /                        / */
 328         /* /     FORMAT : BMP       / */
 329         /* /                        / */
 330         /* /----------------------- / */
 331
 332     case 2:                     /* BMP */
 333         if (img->numcomps == 3 && img->comps[0].dx == img->comps[1].dx
 334             && img->comps[1].dx == img->comps[2].dx
 335             && img->comps[0].dy == img->comps[1].dy
 336             && img->comps[1].dy == img->comps[2].dy
 337             && img->comps[0].prec == img->comps[1].prec
 338             && img->comps[1].prec == img->comps[2].prec) {
 339             /* -->> -->> -->> -->>
 340
 341                24 bits color
 342
 343                <<-- <<-- <<-- <<-- */
 344
 345             f = fopen(argv[2], "wb");
 346             // w = ceildiv(img->x1 - img->x0, img->comps[0].dx);
 347             // wr = ceildiv(int_ceildivpow2(img->x1 - img->x0,img->factor), img->comps[0].dx);
 348             w = img->comps[0].w;
 349             wr = int_ceildivpow2(img->comps[0].w, img->comps[0].factor);
 350
 351             // h = ceildiv(img->y1 - img->y0, img->comps[0].dy);
 352             // hr = ceildiv(int_ceildivpow2(img->y1 - img->y0,img->factor), img->comps[0].dy);
 353             h = img->comps[0].h;
 354             hr = int_ceildivpow2(img->comps[0].h, img->comps[0].factor);
 355
 356             fprintf(f, "BM");
 357
 358             /* FILE HEADER */
 359             /* ------------- */
 360             fprintf(f, "%c%c%c%c",
 361                     (unsigned char) (hr * wr * 3 + 3 * hr * (wr % 2) +
 362                                      54) & 0xff,
 363                     (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 364                                      >> 8) & 0xff,
 365                     (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 366                                      >> 16) & 0xff,
 367                     (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2) + 54)
 368                                      >> 24) & 0xff);
 369             fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 370                     ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 371             fprintf(f, "%c%c%c%c", (54) & 0xff, ((54) >> 8) & 0xff,
 372                     ((54) >> 16) & 0xff, ((54) >> 24) & 0xff);
 373
 374             /* INFO HEADER   */
 375             /* ------------- */
 376             fprintf(f, "%c%c%c%c", (40) & 0xff, ((40) >> 8) & 0xff,
 377                     ((40) >> 16) & 0xff, ((40) >> 24) & 0xff);
 378             fprintf(f, "%c%c%c%c", (unsigned char) ((wr) & 0xff),
 379                     (unsigned char) ((wr) >> 8) & 0xff,
 380                     (unsigned char) ((wr) >> 16) & 0xff,
 381                     (unsigned char) ((wr) >> 24) & 0xff);
 382             fprintf(f, "%c%c%c%c", (unsigned char) ((hr) & 0xff),
 383                     (unsigned char) ((hr) >> 8) & 0xff,
 384                     (unsigned char) ((hr) >> 16) & 0xff,
 385                     (unsigned char) ((hr) >> 24) & 0xff);
 386             fprintf(f, "%c%c", (1) & 0xff, ((1) >> 8) & 0xff);
 387             fprintf(f, "%c%c", (24) & 0xff, ((24) >> 8) & 0xff);
 388             fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 389                     ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 390             fprintf(f, "%c%c%c%c",
 391                     (unsigned char) (3 * hr * wr +
 392                                      3 * hr * (wr % 2)) & 0xff,
 393                     (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 394                                      8) & 0xff,
 395                     (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 396                                      16) & 0xff,
 397                     (unsigned char) ((hr * wr * 3 + 3 * hr * (wr % 2)) >>
 398                                      24) & 0xff);
 399             fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 400                     ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 401             fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 402                     ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 403             fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 404                     ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 405             fprintf(f, "%c%c%c%c", (0) & 0xff, ((0) >> 8) & 0xff,
 406                     ((0) >> 16) & 0xff, ((0) >> 24) & 0xff);
 407
 408             for (i = 0; i < wr * hr; i++) {
 409                 unsigned char R, G, B;
 410                 /* a modifier */
 411                 // R = img->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 412                 R = img->comps[0].data[w * hr - ((i) / (wr) + 1) * w +
 413                                        (i) % (wr)];
 414                 // G = img->comps[1].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 415                 G = img->comps[1].data[w * hr - ((i) / (wr) + 1) * w +
 416                                        (i) % (wr)];
 417                 // B = img->comps[2].data[w * h - ((i) / (w) + 1) * w + (i) % (w)];
 418                 B = img->comps[2].data[w * hr - ((i) / (wr) + 1) * w +
 419                                        (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) &
 488                     0xff);
 489             fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 490                     ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 491             fprintf(f, "%c%c%c%c", (7834) & 0xff, ((7834) >> 8) & 0xff,
 492                     ((7834) >> 16) & 0xff, ((7834) >> 24) & 0xff);
 493             fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
 494                     ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
 495             fprintf(f, "%c%c%c%c", (256) & 0xff, ((256) >> 8) & 0xff,
 496                     ((256) >> 16) & 0xff, ((256) >> 24) & 0xff);
 497         }
 498
 499         for (i = 0; i < 256; i++) {
 500             fprintf(f, "%c%c%c%c", i, i, i, 0);
 501         }
 502
 503         for (i = 0; i < wr * hr; i++) {
 504             /* a modifier !! */
 505             // fprintf(f, "%c", img->comps[0].data[w * h - ((i) / (w) + 1) * w + (i) % (w)]);
 506             fprintf(f, "%c",
 507                     img->comps[0].data[w * hr - ((i) / (wr) + 1) * w +
 508                                        (i) % (wr)]);
 509             /*if (((i + 1) % w == 0 && w % 2))
 510                fprintf(f, "%c", 0); */
 511             if ((i + 1) % wr == 0) {
 512                 for (pad = wr % 4 ? 4 - wr % 4 : 0; pad > 0; pad--)     /* ADD */
 513                     fprintf(f, "%c", 0);
 514             }
 515         }
 516         break;
 517     default:
 518         break;
 519     }
 520
 521     return 0;
 522 }