applications/codec/j2k_dump.c

   1 /*
   2  * Copyright (c) 2010, Mathieu Malaterre, GDCM
   3  * Copyright (c) 2011, Mickael Savinaud, Communications & Systemes <mickael.savinaud@c-s.fr>
   4  * All rights reserved.
   5  *
   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 #include <stdio.h>
  29 #include <string.h>
  30 #include <stdlib.h>
  31 #include <math.h>
  32
  33 #ifdef _WIN32
  34 #include "windirent.h"
  35 #else
  36 #include <dirent.h>
  37 #endif /* _WIN32 */
  38
  39 #ifdef _WIN32
  40 #include <windows.h>
  41 #else
  42 #include <strings.h>
  43 #define _stricmp strcasecmp
  44 #define _strnicmp strncasecmp
  45 #endif /* _WIN32 */
  46
  47 #include "opj_config.h"
  48 #include "openjpeg.h"
  49 #include "j2k.h"
  50 #include "jp2.h"
  51 #include "opj_getopt.h"
  52 #include "convert.h"
  53 #include "index.h"
  54
  55 #include "format_defs.h"
  56
  57 typedef struct dircnt{
  58         /** Buffer for holding images read from Directory*/
  59         char *filename_buf;
  60         /** Pointer to the buffer*/
  61         char **filename;
  62 }dircnt_t;
  63
  64
  65 typedef struct img_folder{
  66         /** The directory path of the folder containing input images*/
  67         char *imgdirpath;
  68         /** Output format*/
  69         const char *out_format;
  70         /** Enable option*/
  71         char set_imgdir;
  72         /** Enable Cod Format for output*/
  73         char set_out_format;
  74
  75 }img_fol_t;
  76
  77 void decode_help_display(void) {
  78         fprintf(stdout,"HELP for j2k_dump\n----\n\n");
  79         fprintf(stdout,"- the -h option displays this help information on screen\n\n");
  80
  81 /* UniPG>> */
  82         fprintf(stdout,"List of parameters for the JPEG 2000 "
  83 #ifdef USE_JPWL
  84                 "+ JPWL "
  85 #endif /* USE_JPWL */
  86                 "decoder:\n");
  87 /* <<UniPG */
  88         fprintf(stdout,"\n");
  89         fprintf(stdout,"\n");
  90         fprintf(stdout,"  -ImgDir \n");
  91         fprintf(stdout,"        Image file Directory path \n");
  92         fprintf(stdout,"  -i <compressed file>\n");
  93         fprintf(stdout,"    REQUIRED only if an Input image directory not specified\n");
  94         fprintf(stdout,"    Currently accepts J2K-files, JP2-files and JPT-files. The file type\n");
  95         fprintf(stdout,"    is identified based on its suffix.\n");
  96         fprintf(stdout,"  -o <output file>\n");
  97         fprintf(stdout,"    OPTIONAL\n");
  98         fprintf(stdout,"    Output file where file info will be dump.\n");
  99         fprintf(stdout,"    By default it will be in the stdout.\n");
 100         fprintf(stdout,"\n");
 101 }
 102
 103 /* -------------------------------------------------------------------------- */
 104 static void j2k_dump_image(FILE *fd, opj_image_header_t * img);
 105 static void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_v2_t * cp);
 106
 107 int get_num_images(char *imgdirpath){
 108         DIR *dir;
 109         struct dirent* content;
 110         int num_images = 0;
 111
 112         /*Reading the input images from given input directory*/
 113
 114         dir= opendir(imgdirpath);
 115         if(!dir){
 116                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
 117                 return 0;
 118         }
 119
 120         while((content=readdir(dir))!=NULL){
 121                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
 122                         continue;
 123                 num_images++;
 124         }
 125         return num_images;
 126 }
 127
 128 int load_images(dircnt_t *dirptr, char *imgdirpath){
 129         DIR *dir;
 130         struct dirent* content;
 131         int i = 0;
 132
 133         /*Reading the input images from given input directory*/
 134
 135         dir= opendir(imgdirpath);
 136         if(!dir){
 137                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
 138                 return 1;
 139         }else   {
 140                 fprintf(stderr,"Folder opened successfully\n");
 141         }
 142
 143         while((content=readdir(dir))!=NULL){
 144                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
 145                         continue;
 146
 147                 strcpy(dirptr->filename[i],content->d_name);
 148                 i++;
 149         }
 150         return 0;
 151 }
 152
 153 int get_file_format(char *filename) {
 154         unsigned int i;
 155         static const char *extension[] = {"pgx", "pnm", "pgm", "ppm", "bmp","tif", "raw", "tga", "png", "j2k", "jp2", "jpt", "j2c", "jpc"  };
 156         static const int format[] = { PGX_DFMT, PXM_DFMT, PXM_DFMT, PXM_DFMT, BMP_DFMT, TIF_DFMT, RAW_DFMT, TGA_DFMT, PNG_DFMT, J2K_CFMT, JP2_CFMT, JPT_CFMT, J2K_CFMT, J2K_CFMT };
 157         char * ext = strrchr(filename, '.');
 158         if (ext == NULL)
 159                 return -1;
 160         ext++;
 161         if(ext) {
 162                 for(i = 0; i < sizeof(format)/sizeof(*format); i++) {
 163                         if(_strnicmp(ext, extension[i], 3) == 0) {
 164                                 return format[i];
 165                         }
 166                 }
 167         }
 168
 169         return -1;
 170 }
 171
 172 char get_next_file(int imageno,dircnt_t *dirptr,img_fol_t *img_fol, opj_dparameters_t *parameters){
 173         char image_filename[OPJ_PATH_LEN], infilename[OPJ_PATH_LEN],outfilename[OPJ_PATH_LEN],temp_ofname[OPJ_PATH_LEN];
 174         char *temp_p, temp1[OPJ_PATH_LEN]="";
 175
 176         strcpy(image_filename,dirptr->filename[imageno]);
 177         fprintf(stderr,"File Number %d \"%s\"\n",imageno,image_filename);
 178         parameters->decod_format = get_file_format(image_filename);
 179         if (parameters->decod_format == -1)
 180                 return 1;
 181         sprintf(infilename,"%s/%s",img_fol->imgdirpath,image_filename);
 182         strncpy(parameters->infile, infilename, sizeof(infilename));
 183
 184         /* Set output file */
 185         strcpy(temp_ofname,strtok(image_filename,"."));
 186         while((temp_p = strtok(NULL,".")) != NULL){
 187                 strcat(temp_ofname,temp1);
 188                 sprintf(temp1,".%s",temp_p);
 189         }
 190         if(img_fol->set_out_format==1){
 191                 sprintf(outfilename,"%s/%s.%s",img_fol->imgdirpath,temp_ofname,img_fol->out_format);
 192                 strncpy(parameters->outfile, outfilename, sizeof(outfilename));
 193         }
 194         return 0;
 195 }
 196
 197 /* -------------------------------------------------------------------------- */
 198 /* parse the command line */
 199 int parse_cmdline_decoder(int argc, char **argv, opj_dparameters_t *parameters,img_fol_t *img_fol) {
 200         int totlen, c;
 201         opj_option_t long_option[]={
 202                 {"ImgDir",REQ_ARG, NULL ,'y'},
 203         };
 204         const char optlist[] = "i:o:d:h";
 205
 206         totlen=sizeof(long_option);
 207         img_fol->set_out_format = 0;
 208         do {
 209                 c = opj_getopt_long(argc, argv,optlist,long_option,totlen);
 210                 if (c == -1)
 211                         break;
 212                 switch (c) {
 213                         case 'i':                       /* input file */
 214                         {
 215                                 char *infile = opj_optarg;
 216                                 parameters->decod_format = get_file_format(infile);
 217                                 switch(parameters->decod_format) {
 218                                         case J2K_CFMT:
 219                                                 break;
 220                                         case JP2_CFMT:
 221                                                 break;
 222                                         case JPT_CFMT:
 223                                                 break;
 224                                         default:
 225                                                 fprintf(stderr,
 226                                                         "!! Unrecognized format for infile : %s [accept only *.j2k, *.jp2, *.jpc or *.jpt] !!\n\n",
 227                                                         infile);
 228                                                 return 1;
 229                                 }
 230                                 strncpy(parameters->infile, infile, sizeof(parameters->infile)-1);
 231                         }
 232                         break;
 233
 234                                 /* ------------------------------------------------------ */
 235
 236                         case 'o':     /* output file */
 237                         {
 238                           char *outfile = opj_optarg;
 239                           strncpy(parameters->outfile, outfile, sizeof(parameters->outfile)-1);
 240                         }
 241                         break;
 242
 243                                 /* ----------------------------------------------------- */
 244
 245                         case 'h':                       /* display an help description */
 246                                 decode_help_display();
 247                                 return 1;
 248
 249                                 /* ------------------------------------------------------ */
 250
 251                         case 'y':                       /* Image Directory path */
 252                                 {
 253                                         img_fol->imgdirpath = (char*)malloc(strlen(opj_optarg) + 1);
 254                                         strcpy(img_fol->imgdirpath,opj_optarg);
 255                                         img_fol->set_imgdir=1;
 256                                 }
 257                                 break;
 258
 259                                 /* ----------------------------------------------------- */
 260
 261                         case 'd':               /* Input decode ROI */
 262                         {
 263                                 int size_optarg = (int)strlen(optarg) + 1;
 264                                 char *ROI_values = (char*) malloc(size_optarg);
 265                                 ROI_values[0] = '\0';
 266                                 strncpy(ROI_values, optarg, strlen(optarg));
 267                                 ROI_values[strlen(optarg)] = '\0';
 268                                 /*printf("ROI_values = %s [%d / %d]\n", ROI_values, strlen(ROI_values), size_optarg ); */
 269                                 parse_ROI_values( ROI_values, &parameters->ROI_x0, &parameters->ROI_y0, &parameters->ROI_x1, &parameters->ROI_y1);
 270                         }
 271                         break;
 272
 273                                 /* ----------------------------------------------------- */
 274                         default:
 275                                 fprintf(stderr,"WARNING -> this option is not valid \"-%c %s\"\n",c, opj_optarg);
 276                                 break;
 277                 }
 278         }while(c != -1);
 279
 280         /* check for possible errors */
 281         if(img_fol->set_imgdir==1){
 282                 if(!(parameters->infile[0]==0)){
 283                         fprintf(stderr, "Error: options -ImgDir and -i cannot be used together !!\n");
 284                         return 1;
 285                 }
 286                 if(img_fol->set_out_format == 0){
 287                         fprintf(stderr, "Error: When -ImgDir is used, -OutFor <FORMAT> must be used !!\n");
 288                         fprintf(stderr, "Only one format allowed! Valid format PGM, PPM, PNM, PGX, BMP, TIF, RAW and TGA!!\n");
 289                         return 1;
 290                 }
 291                 if(!((parameters->outfile[0] == 0))){
 292                         fprintf(stderr, "Error: options -ImgDir and -o cannot be used together !!\n");
 293                         return 1;
 294                 }
 295         }else{
 296                 if((parameters->infile[0] == 0) ) {
 297                         fprintf(stderr, "Example: %s -i image.j2k\n",argv[0]);
 298                         fprintf(stderr, "    Try: %s -h\n",argv[0]);
 299                         return 1;
 300                 }
 301         }
 302
 303         return 0;
 304 }
 305
 306 /*******************************************************************************
 307  * Parse ROI input values
 308  * separator = ","
 309  *******************************************************************************/
 310 int parse_ROI_values( char* inArg, unsigned int *ROI_x0, unsigned int *ROI_y0, unsigned int *ROI_x1, unsigned int *ROI_y1)
 311 {
 312         int it = 0;
 313         int values[4];
 314         char delims[] = ",";
 315         char *result = NULL;
 316         result = strtok( inArg, delims );
 317
 318         while( (result != NULL) && (it < 4 ) ) {
 319                 values[it] = atoi(result);
 320                 result = strtok( NULL, delims );
 321                 it++;
 322         }
 323
 324         if (it != 4) {
 325                 return EXIT_FAILURE;
 326         }
 327         else{
 328                 *ROI_x0 = values[0]; *ROI_y0 = values[1];
 329                 *ROI_x1 = values[2]; *ROI_y1 = values[3];
 330                 return EXIT_SUCCESS;
 331         }
 332 }
 333
 334 /* -------------------------------------------------------------------------- */
 335
 336 /**
 337 sample error callback expecting a FILE* client object
 338 */
 339 void error_callback(const char *msg, void *client_data) {
 340         FILE *stream = (FILE*)client_data;
 341         fprintf(stream, "[ERROR] %s", msg);
 342 }
 343 /**
 344 sample warning callback expecting a FILE* client object
 345 */
 346 void warning_callback(const char *msg, void *client_data) {
 347         FILE *stream = (FILE*)client_data;
 348         fprintf(stream, "[WARNING] %s", msg);
 349 }
 350 /**
 351 sample debug callback expecting no client object
 352 */
 353 void info_callback(const char *msg, void *client_data) {
 354         (void)client_data;
 355         fprintf(stdout, "[INFO] %s", msg);
 356 }
 357
 358 /* -------------------------------------------------------------------------- */
 359
 360 int main(int argc, char *argv[])
 361 {
 362         int ret;
 363         opj_dparameters_t parameters;   /* decompression parameters */
 364         img_fol_t img_fol;
 365         opj_event_mgr_t event_mgr;              /* event manager */
 366         opj_image_header_t* image = NULL;
 367         opj_file_info_t file_info;
 368         FILE *fsrc = NULL, *fout = NULL;
 369         int num_images;
 370         int i,imageno;
 371         dircnt_t *dirptr = NULL;
 372         opj_codec_t* dinfo = NULL;      /* handle to a decompressor */
 373         opj_stream_t *cio = NULL;
 374         opj_codestream_info_t* cstr_info =NULL;  /* Codestream information structure */
 375         opj_bool l_go_on = OPJ_TRUE;
 376
 377
 378         /* FIXME configure the event callbacks (not required) */
 379         memset(&event_mgr, 0, sizeof(opj_event_mgr_t));
 380         event_mgr.error_handler = error_callback;
 381         event_mgr.warning_handler = warning_callback;
 382         event_mgr.info_handler = info_callback;
 383         event_mgr.client_data = stderr;
 384
 385         /* set decoding parameters to default values */
 386         opj_set_default_decoder_parameters(&parameters);
 387
 388         /* Initialize img_fol */
 389         memset(&img_fol,0,sizeof(img_fol_t));
 390
 391         /* parse input and get user encoding parameters */
 392         if(parse_cmdline_decoder(argc, argv, &parameters,&img_fol) == 1) {
 393                 return EXIT_FAILURE;
 394         }
 395
 396         /* Initialize reading of directory */
 397         if(img_fol.set_imgdir==1){
 398                 num_images=get_num_images(img_fol.imgdirpath);
 399
 400                 dirptr=(dircnt_t*)malloc(sizeof(dircnt_t));
 401                 if(dirptr){
 402                         dirptr->filename_buf = (char*)malloc(num_images*OPJ_PATH_LEN*sizeof(char));     /* Stores at max 10 image file names */
 403                         dirptr->filename = (char**) malloc(num_images*sizeof(char*));
 404
 405                         if(!dirptr->filename_buf){
 406                                 return EXIT_FAILURE;
 407                         }
 408
 409                         for(i=0;i<num_images;i++){
 410                                 dirptr->filename[i] = dirptr->filename_buf + i*OPJ_PATH_LEN;
 411                         }
 412                 }
 413                 if(load_images(dirptr,img_fol.imgdirpath)==1){
 414                         return EXIT_FAILURE;
 415                 }
 416
 417                 if (num_images==0){
 418                         fprintf(stdout,"Folder is empty\n");
 419                         return EXIT_FAILURE;
 420                 }
 421         }else{
 422                 num_images=1;
 423         }
 424
 425         /* Try to open for writing the output file if necessary */
 426         if (parameters.outfile[0] != 0){
 427                 fout = fopen(parameters.outfile,"w");
 428                 if (!fout){
 429                         fprintf(stderr, "ERROR -> failed to open %s for writing\n", parameters.outfile);
 430                         return EXIT_FAILURE;
 431                 }
 432         }
 433         else
 434                 fout = stdout;
 435
 436         /* Read the header of each image one by one */
 437         for(imageno = 0; imageno < num_images ; imageno++){
 438
 439                 fprintf(stderr,"\n");
 440
 441                 if(img_fol.set_imgdir==1){
 442                         if (get_next_file(imageno, dirptr,&img_fol, &parameters)) {
 443                                 fprintf(stderr,"skipping file...\n");
 444                                 continue;
 445                         }
 446                 }
 447
 448                 /* read the input file and put it in memory */
 449                 /* ---------------------------------------- */
 450                 fsrc = fopen(parameters.infile, "rb");
 451                 if (!fsrc) {
 452                         fprintf(stderr, "ERROR -> failed to open %s for reading\n", parameters.infile);
 453                         return EXIT_FAILURE;
 454                 }
 455
 456                 cio = opj_stream_create_default_file_stream(fsrc,1);
 457                 if (!cio){
 458                         fprintf(stderr, "ERROR -> failed to create the stream from the file\n");
 459                         return EXIT_FAILURE;
 460                 }
 461
 462                 /* decode the code-stream */
 463                 /* ---------------------- */
 464
 465                 switch(parameters.decod_format) {
 466                         case J2K_CFMT:
 467                         {       /* JPEG-2000 codestream */
 468
 469                                 /* get a decoder handle */
 470                                 dinfo = opj_create_decompress_v2(CODEC_J2K);
 471                                 break;
 472                         }
 473                         case JP2_CFMT:
 474                         {       /* JPEG 2000 compressed image data */
 475
 476                                 /* get a decoder handle */
 477                                 dinfo = opj_create_decompress_v2(CODEC_JP2);
 478                                 break;
 479                         }
 480                         case JPT_CFMT:
 481                         {       /* JPEG 2000, JPIP */
 482
 483                                 /* get a decoder handle */
 484                                 dinfo = opj_create_decompress_v2(CODEC_JPT);
 485                                 break;
 486                         }
 487                         default:
 488                                 fprintf(stderr, "skipping file..\n");
 489                                 opj_stream_destroy(cio);
 490                                 continue;
 491                 }
 492
 493                 /* setup the decoder decoding parameters using user parameters */
 494                 opj_setup_decoder_v2(dinfo, &parameters, &event_mgr);
 495
 496                 if(! opj_read_header(cio, dinfo, &file_info, OPJ_IMG_INFO | OPJ_J2K_INFO)){
 497                         fprintf(stderr, "ERROR -> j2k_dump: failed to read the header\n");
 498                         opj_stream_destroy(cio);
 499                         fclose(fsrc);
 500                         opj_destroy_codec(dinfo);
 501                         return EXIT_FAILURE;
 502                 }
 503
 504                 printf("Setting decoding area to %d,%d,%d,%d\n",
 505                                 parameters.ROI_x0, parameters.ROI_y0, parameters.ROI_x1, parameters.ROI_y1);
 506
 507                 if (! opj_set_decode_area(      dinfo,
 508                                                                         parameters.ROI_x0, parameters.ROI_y0,
 509                                                                         parameters.ROI_x1, parameters.ROI_y1)){
 510                         fprintf(stderr, "ERROR -> j2k_dump: failed to set the decoded area\n");
 511                         opj_stream_destroy(cio);
 512                         fclose(fsrc);
 513                         opj_destroy_codec(dinfo);
 514                         return EXIT_FAILURE;
 515                 }
 516
 517                 while (l_go_on) {
 518                         OPJ_INT32 l_current_tile_x0,l_current_tile_y0,l_current_tile_x1,l_current_tile_y1;
 519                         OPJ_UINT32 l_nb_comps, l_tile_index, l_data_size;
 520
 521                         if (! opj_read_tile_header(     dinfo,
 522                                                                                 cio,
 523                                                                                 &l_tile_index,
 524                                                                                 &l_data_size,
 525                                                                                 &l_current_tile_x0,
 526                                                                                 &l_current_tile_y0,
 527                                                                                 &l_current_tile_x1,
 528                                                                                 &l_current_tile_y1,
 529                                                                                 &l_nb_comps,
 530                                                                                 &l_go_on
 531                                                                                 )) {
 532                                 fprintf(stderr, "ERROR -> j2k_dump: failed read the tile header\n");
 533                                 opj_stream_destroy(cio);
 534                                 fclose(fsrc);
 535                                 opj_destroy_codec(dinfo);
 536                                 return EXIT_FAILURE;
 537                         }
 538                 }
 539
 540                 /* Dump file informations from header */
 541                 dump_file_info(fout, &file_info);
 542
 543                 /* close the byte stream */
 544                 opj_stream_destroy(cio);
 545                 fclose(fsrc);
 546
 547                 /* free remaining structures */
 548                 if (dinfo) {
 549                         opj_destroy_codec(dinfo);
 550                 }
 551
 552
 553                 opj_image_header_destroy(image);
 554                 //FIXME opj_file_info_destroy(file_info);
 555
 556         }
 557
 558
 559         /* Close the output file */
 560         fclose(fout);
 561
 562   return EXIT_SUCCESS;
 563 }
 564
 565
 566 static void j2k_dump_image(FILE *fd, opj_image_header_t * img) {
 567         int compno;
 568         fprintf(fd, "image {\n");
 569         fprintf(fd, "  x0=%d, y0=%d, x1=%d, y1=%d\n", img->x0, img->y0, img->x1, img->y1);
 570         fprintf(fd, "  numcomps=%d\n", img->numcomps);
 571         for (compno = 0; compno < img->numcomps; compno++) {
 572                 opj_image_comp_header_t *comp = &img->comps[compno];
 573                 fprintf(fd, "  comp %d {\n", compno);
 574                 fprintf(fd, "    dx=%d, dy=%d\n", comp->dx, comp->dy);
 575                 fprintf(fd, "    prec=%d\n", comp->prec);
 576                 /* fprintf(fd, "    bpp=%d\n", comp->bpp); */
 577                 fprintf(fd, "    sgnd=%d\n", comp->sgnd);
 578                 fprintf(fd, "  }\n");
 579         }
 580         //fprintf(fd, "  XTOsiz=%d, YTOsiz=%d, XTsiz=%d, YTsiz=%d\n", img->tile_x0, img->tile_y0, img->tile_width, img->tile_height);
 581         //fprintf(fd, "  Nb of tiles in x direction=%d, Nb of tiles in y direction=%d\n", img->nb_tiles_x, img->nb_tiles_y);
 582         fprintf(fd, "}\n");
 583 }
 584
 585 static void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_v2_t * cp) {
 586         int tileno, compno, layno, bandno, resno, numbands;
 587         fprintf(fd, "coding parameters {\n");
 588         fprintf(fd, "  tx0=%d, ty0=%d\n", cp->tx0, cp->ty0);
 589         fprintf(fd, "  tdx=%d, tdy=%d\n", cp->tdx, cp->tdy);
 590         fprintf(fd, "  tw=%d, th=%d\n", cp->tw, cp->th);
 591         for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
 592                 opj_tcp_v2_t *tcp = &cp->tcps[tileno];
 593                 fprintf(fd, "  tile %d {\n", tileno);
 594                 fprintf(fd, "    csty=%x\n", tcp->csty);
 595                 fprintf(fd, "    prg=%d\n", tcp->prg);
 596                 fprintf(fd, "    numlayers=%d\n", tcp->numlayers);
 597                 fprintf(fd, "    mct=%d\n", tcp->mct);
 598                 fprintf(fd, "    rates=");
 599                 for (layno = 0; layno < tcp->numlayers; layno++) {
 600                         fprintf(fd, "%.1f ", tcp->rates[layno]);
 601                 }
 602                 fprintf(fd, "\n");
 603                 for (compno = 0; compno < img->numcomps; compno++) {
 604                         opj_tccp_t *tccp = &tcp->tccps[compno];
 605                         fprintf(fd, "    comp %d {\n", compno);
 606                         fprintf(fd, "      csty=%x\n", tccp->csty);
 607                         fprintf(fd, "      numresolutions=%d\n", tccp->numresolutions);
 608                         fprintf(fd, "      cblkw=%d\n", tccp->cblkw);
 609                         fprintf(fd, "      cblkh=%d\n", tccp->cblkh);
 610                         fprintf(fd, "      cblksty=%x\n", tccp->cblksty);
 611                         fprintf(fd, "      qmfbid=%d\n", tccp->qmfbid);
 612                         fprintf(fd, "      qntsty=%d\n", tccp->qntsty);
 613                         fprintf(fd, "      numgbits=%d\n", tccp->numgbits);
 614                         fprintf(fd, "      roishift=%d\n", tccp->roishift);
 615                         fprintf(fd, "      stepsizes=");
 616                         numbands = tccp->qntsty == J2K_CCP_QNTSTY_SIQNT ? 1 : tccp->numresolutions * 3 - 2;
 617                         for (bandno = 0; bandno < numbands; bandno++) {
 618                                 fprintf(fd, "(%d,%d) ", tccp->stepsizes[bandno].mant,
 619                                         tccp->stepsizes[bandno].expn);
 620                         }
 621                         fprintf(fd, "\n");
 622
 623                         if (tccp->csty & J2K_CCP_CSTY_PRT) {
 624                                 fprintf(fd, "      prcw=");
 625                                 for (resno = 0; resno < tccp->numresolutions; resno++) {
 626                                         fprintf(fd, "%d ", tccp->prcw[resno]);
 627                                 }
 628                                 fprintf(fd, "\n");
 629                                 fprintf(fd, "      prch=");
 630                                 for (resno = 0; resno < tccp->numresolutions; resno++) {
 631                                         fprintf(fd, "%d ", tccp->prch[resno]);
 632                                 }
 633                                 fprintf(fd, "\n");
 634                         }
 635                         fprintf(fd, "    }\n");
 636                 } /*end of component*/
 637                 fprintf(fd, "  }\n");
 638         } /*end of tile */
 639         fprintf(fd, "}\n");
 640 }
 641