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