applications/codec/j2k_dump.c

   1 /*
   2  * Copyright (c) 20010, Mathieu Malaterre, GDCM
   3  * All rights reserved.
   4  *
   5  * Redistribution and use in source and binary forms, with or without
   6  * modification, are permitted provided that the following conditions
   7  * are met:
   8  * 1. Redistributions of source code must retain the above copyright
   9  *    notice, this list of conditions and the following disclaimer.
  10  * 2. Redistributions in binary form must reproduce the above copyright
  11  *    notice, this list of conditions and the following disclaimer in the
  12  *    documentation and/or other materials provided with the distribution.
  13  *
  14  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
  15  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  16  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  17  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  18  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  19  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  20  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  21  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  22  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  23  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  24  * POSSIBILITY OF SUCH DAMAGE.
  25  */
  26
  27 #include <stdio.h>
  28 #include <string.h>
  29 #include <stdlib.h>
  30 #include <math.h>
  31
  32 #ifdef _WIN32
  33 #include "windirent.h"
  34 #else
  35 #include <dirent.h>
  36 #endif /* _WIN32 */
  37
  38 #ifdef _WIN32
  39 #include <windows.h>
  40 #else
  41 #include <strings.h>
  42 #define _stricmp strcasecmp
  43 #define _strnicmp strncasecmp
  44 #endif /* _WIN32 */
  45
  46 #include "opj_config.h"
  47 #include "openjpeg.h"
  48 #include "j2k.h"
  49 #include "jp2.h"
  50 #include "opj_getopt.h"
  51 #include "convert.h"
  52 #include "index.h"
  53
  54 #include "format_defs.h"
  55
  56 typedef struct dircnt{
  57         /** Buffer for holding images read from Directory*/
  58         char *filename_buf;
  59         /** Pointer to the buffer*/
  60         char **filename;
  61 }dircnt_t;
  62
  63
  64 typedef struct img_folder{
  65         /** The directory path of the folder containing input images*/
  66         char *imgdirpath;
  67         /** Output format*/
  68         const char *out_format;
  69         /** Enable option*/
  70         char set_imgdir;
  71         /** Enable Cod Format for output*/
  72         char set_out_format;
  73
  74 }img_fol_t;
  75
  76 void decode_help_display(void) {
  77         fprintf(stdout,"HELP for j2k_dump\n----\n\n");
  78         fprintf(stdout,"- the -h option displays this help information on screen\n\n");
  79
  80 /* UniPG>> */
  81         fprintf(stdout,"List of parameters for the JPEG 2000 "
  82 #ifdef USE_JPWL
  83                 "+ JPWL "
  84 #endif /* USE_JPWL */
  85                 "decoder:\n");
  86 /* <<UniPG */
  87         fprintf(stdout,"\n");
  88         fprintf(stdout,"\n");
  89         fprintf(stdout,"  -ImgDir \n");
  90         fprintf(stdout,"        Image file Directory path \n");
  91         fprintf(stdout,"  -i <compressed file>\n");
  92         fprintf(stdout,"    REQUIRED only if an Input image directory not specified\n");
  93         fprintf(stdout,"    Currently accepts J2K-files, JP2-files and JPT-files. The file type\n");
  94         fprintf(stdout,"    is identified based on its suffix.\n");
  95         fprintf(stdout,"  -o <output file>\n");
  96   fprintf(stdout,"    OPTIONAL\n");
  97   fprintf(stdout,"    Output file where file info will be dump.\n");
  98   fprintf(stdout,"    By default it will be in the stdout.\n");
  99         fprintf(stdout,"\n");
 100 }
 101
 102 /* -------------------------------------------------------------------------- */
 103 static void j2k_dump_image(FILE *fd, opj_image_t * img);
 104 static void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_t * cp);
 105
 106 int get_num_images(char *imgdirpath){
 107         DIR *dir;
 108         struct dirent* content;
 109         int num_images = 0;
 110
 111         /*Reading the input images from given input directory*/
 112
 113         dir= opendir(imgdirpath);
 114         if(!dir){
 115                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
 116                 return 0;
 117         }
 118
 119         while((content=readdir(dir))!=NULL){
 120                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
 121                         continue;
 122                 num_images++;
 123         }
 124         return num_images;
 125 }
 126
 127 int load_images(dircnt_t *dirptr, char *imgdirpath){
 128         DIR *dir;
 129         struct dirent* content;
 130         int i = 0;
 131
 132         /*Reading the input images from given input directory*/
 133
 134         dir= opendir(imgdirpath);
 135         if(!dir){
 136                 fprintf(stderr,"Could not open Folder %s\n",imgdirpath);
 137                 return 1;
 138         }else   {
 139                 fprintf(stderr,"Folder opened successfully\n");
 140         }
 141
 142         while((content=readdir(dir))!=NULL){
 143                 if(strcmp(".",content->d_name)==0 || strcmp("..",content->d_name)==0 )
 144                         continue;
 145
 146                 strcpy(dirptr->filename[i],content->d_name);
 147                 i++;
 148         }
 149         return 0;
 150 }
 151
 152 int get_file_format(char *filename) {
 153         unsigned int i;
 154         static const char *extension[] = {"pgx", "pnm", "pgm", "ppm", "bmp","tif", "raw", "tga", "png", "j2k", "jp2", "jpt", "j2c", "jpc"  };
 155         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 };
 156         char * ext = strrchr(filename, '.');
 157         if (ext == NULL)
 158                 return -1;
 159         ext++;
 160         if(ext) {
 161                 for(i = 0; i < sizeof(format)/sizeof(*format); i++) {
 162                         if(_strnicmp(ext, extension[i], 3) == 0) {
 163                                 return format[i];
 164                         }
 165                 }
 166         }
 167
 168         return -1;
 169 }
 170
 171 char get_next_file(int imageno,dircnt_t *dirptr,img_fol_t *img_fol, opj_dparameters_t *parameters){
 172         char image_filename[OPJ_PATH_LEN], infilename[OPJ_PATH_LEN],outfilename[OPJ_PATH_LEN],temp_ofname[OPJ_PATH_LEN];
 173         char *temp_p, temp1[OPJ_PATH_LEN]="";
 174
 175         strcpy(image_filename,dirptr->filename[imageno]);
 176         fprintf(stderr,"File Number %d \"%s\"\n",imageno,image_filename);
 177         parameters->decod_format = get_file_format(image_filename);
 178         if (parameters->decod_format == -1)
 179                 return 1;
 180         sprintf(infilename,"%s/%s",img_fol->imgdirpath,image_filename);
 181         strncpy(parameters->infile, infilename, sizeof(infilename));
 182
 183         //Set output file
 184         strcpy(temp_ofname,strtok(image_filename,"."));
 185         while((temp_p = strtok(NULL,".")) != NULL){
 186                 strcat(temp_ofname,temp1);
 187                 sprintf(temp1,".%s",temp_p);
 188         }
 189         if(img_fol->set_out_format==1){
 190                 sprintf(outfilename,"%s/%s.%s",img_fol->imgdirpath,temp_ofname,img_fol->out_format);
 191                 strncpy(parameters->outfile, outfilename, sizeof(outfilename));
 192         }
 193         return 0;
 194 }
 195
 196 /* -------------------------------------------------------------------------- */
 197 int parse_cmdline_decoder(int argc, char **argv, opj_dparameters_t *parameters,img_fol_t *img_fol, char *indexfilename) {
 198         /* parse the command line */
 199         int totlen, c;
 200         opj_option_t long_option[]={
 201                 {"ImgDir",REQ_ARG, NULL ,'y'},
 202         };
 203         const char optlist[] = "i:o:h";
 204
 205         OPJ_ARG_NOT_USED(indexfilename);
 206
 207         totlen=sizeof(long_option);
 208         img_fol->set_out_format = 0;
 209         do {
 210                 c = opj_getopt_long(argc, argv,optlist,long_option,totlen);
 211                 if (c == -1)
 212                         break;
 213                 switch (c) {
 214                         case 'i':                       /* input file */
 215                         {
 216                                 char *infile = opj_optarg;
 217                                 parameters->decod_format = get_file_format(infile);
 218                                 switch(parameters->decod_format) {
 219                                         case J2K_CFMT:
 220                                         case JP2_CFMT:
 221                                         case JPT_CFMT:
 222                                                 break;
 223                                         default:
 224                                                 fprintf(stderr,
 225                                                         "!! Unrecognized format for infile : %s [accept only *.j2k, *.jp2, *.jpc or *.jpt] !!\n\n",
 226                                                         infile);
 227                                                 return 1;
 228                                 }
 229                                 strncpy(parameters->infile, infile, sizeof(parameters->infile)-1);
 230                         }
 231                         break;
 232
 233                         /* ------------------------------------------------------ */
 234
 235                         case 'o':     /* output file */
 236                         {
 237                           char *outfile = opj_optarg;
 238                           strncpy(parameters->outfile, outfile, sizeof(parameters->outfile)-1);
 239                         }
 240                         break;
 241
 242                                 /* ----------------------------------------------------- */
 243
 244                         case 'h':                       /* display an help description */
 245                                 decode_help_display();
 246                                 return 1;
 247
 248                                 /* ------------------------------------------------------ */
 249
 250                         case 'y':                       /* Image Directory path */
 251                                 {
 252                                         img_fol->imgdirpath = (char*)malloc(strlen(opj_optarg) + 1);
 253                                         strcpy(img_fol->imgdirpath,opj_optarg);
 254                                         img_fol->set_imgdir=1;
 255                                 }
 256                                 break;
 257
 258                                 /* ----------------------------------------------------- */
 259
 260                         default:
 261                                 fprintf(stderr,"WARNING -> this option is not valid \"-%c %s\"\n",c, opj_optarg);
 262                                 break;
 263                 }
 264         }while(c != -1);
 265
 266         /* check for possible errors */
 267         if(img_fol->set_imgdir==1){
 268                 if(!(parameters->infile[0]==0)){
 269                         fprintf(stderr, "Error: options -ImgDir and -i cannot be used together !!\n");
 270                         return 1;
 271                 }
 272                 if(img_fol->set_out_format == 0){
 273                         fprintf(stderr, "Error: When -ImgDir is used, -OutFor <FORMAT> must be used !!\n");
 274                         fprintf(stderr, "Only one format allowed! Valid format PGM, PPM, PNM, PGX, BMP, TIF, RAW and TGA!!\n");
 275                         return 1;
 276                 }
 277                 if(!((parameters->outfile[0] == 0))){
 278                         fprintf(stderr, "Error: options -ImgDir and -o cannot be used together !!\n");
 279                         return 1;
 280                 }
 281         }else{
 282                 if((parameters->infile[0] == 0) ) {
 283                         fprintf(stderr, "Example: %s -i image.j2k\n",argv[0]);
 284                         fprintf(stderr, "    Try: %s -h\n",argv[0]);
 285                         return 1;
 286                 }
 287         }
 288
 289         return 0;
 290 }
 291
 292 /* -------------------------------------------------------------------------- */
 293
 294 /**
 295 sample error callback expecting a FILE* client object
 296 */
 297 void error_callback(const char *msg, void *client_data) {
 298         FILE *stream = (FILE*)client_data;
 299         fprintf(stream, "[ERROR] %s", msg);
 300 }
 301 /**
 302 sample warning callback expecting a FILE* client object
 303 */
 304 void warning_callback(const char *msg, void *client_data) {
 305         FILE *stream = (FILE*)client_data;
 306         fprintf(stream, "[WARNING] %s", msg);
 307 }
 308 /**
 309 sample debug callback expecting no client object
 310 */
 311 void info_callback(const char *msg, void *client_data) {
 312         (void)client_data;
 313         fprintf(stdout, "[INFO] %s", msg);
 314 }
 315
 316 /* -------------------------------------------------------------------------- */
 317
 318 int main(int argc, char *argv[])
 319 {
 320         opj_dparameters_t parameters;   /* decompression parameters */
 321         img_fol_t img_fol;
 322         opj_event_mgr_t event_mgr;              /* event manager */
 323         opj_image_t *image = NULL;
 324         FILE *fsrc = NULL, *fout = NULL;
 325         unsigned char *src = NULL;
 326         int file_length;
 327         int num_images;
 328         int i,imageno;
 329         dircnt_t *dirptr = NULL;
 330         opj_dinfo_t* dinfo = NULL;      /* handle to a decompressor */
 331         opj_cio_t *cio = NULL;
 332         opj_codestream_info_t cstr_info;  /* Codestream information structure */
 333         char indexfilename[OPJ_PATH_LEN];       /* index file name */
 334
 335         /* configure the event callbacks (not required) */
 336         memset(&event_mgr, 0, sizeof(opj_event_mgr_t));
 337         event_mgr.error_handler = error_callback;
 338         event_mgr.warning_handler = warning_callback;
 339         event_mgr.info_handler = info_callback;
 340
 341         /* set decoding parameters to default values */
 342         opj_set_default_decoder_parameters(&parameters);
 343
 344         /* Initialize indexfilename and img_fol */
 345         *indexfilename = 0;
 346         memset(&img_fol,0,sizeof(img_fol_t));
 347
 348         /* parse input and get user encoding parameters */
 349         if(parse_cmdline_decoder(argc, argv, &parameters,&img_fol, indexfilename) == 1) {
 350                 return 1;
 351         }
 352
 353         /* Initialize reading of directory */
 354         if(img_fol.set_imgdir==1){
 355                 num_images=get_num_images(img_fol.imgdirpath);
 356
 357                 dirptr=(dircnt_t*)malloc(sizeof(dircnt_t));
 358                 if(dirptr){
 359                         dirptr->filename_buf = (char*)malloc(num_images*OPJ_PATH_LEN*sizeof(char));     // Stores at max 10 image file names
 360                         dirptr->filename = (char**) malloc(num_images*sizeof(char*));
 361
 362                         if(!dirptr->filename_buf){
 363                                 return 1;
 364                         }
 365                         for(i=0;i<num_images;i++){
 366                                 dirptr->filename[i] = dirptr->filename_buf + i*OPJ_PATH_LEN;
 367                         }
 368                 }
 369                 if(load_images(dirptr,img_fol.imgdirpath)==1){
 370                         return 1;
 371                 }
 372                 if (num_images==0){
 373                         fprintf(stdout,"Folder is empty\n");
 374                         return 1;
 375                 }
 376         }else{
 377                 num_images=1;
 378         }
 379
 380         //
 381         if (parameters.outfile[0] != 0)
 382           {
 383           fout = fopen(parameters.outfile,"w");
 384     if (!fout)
 385       {
 386       fprintf(stderr, "ERROR -> failed to open %s for reading\n", parameters.outfile);
 387       return 1;
 388       }
 389           }
 390         else
 391           fout = stdout;
 392
 393         /*Encoding image one by one*/
 394         for(imageno = 0; imageno < num_images ; imageno++)
 395   {
 396                 image = NULL;
 397                 fprintf(stderr,"\n");
 398
 399                 if(img_fol.set_imgdir==1){
 400                         if (get_next_file(imageno, dirptr,&img_fol, &parameters)) {
 401                                 fprintf(stderr,"skipping file...\n");
 402                                 continue;
 403                         }
 404                 }
 405
 406                 /* read the input file and put it in memory */
 407                 /* ---------------------------------------- */
 408                 fsrc = fopen(parameters.infile, "rb");
 409                 if (!fsrc) {
 410                         fprintf(stderr, "ERROR -> failed to open %s for reading\n", parameters.infile);
 411                         return 1;
 412                 }
 413                 fseek(fsrc, 0, SEEK_END);
 414                 file_length = ftell(fsrc);
 415                 fseek(fsrc, 0, SEEK_SET);
 416                 src = (unsigned char *) malloc(file_length);
 417                 if (fread(src, 1, file_length, fsrc) != (size_t)file_length)
 418                 {
 419                         free(src);
 420                         fclose(fsrc);
 421                         fprintf(stderr, "\nERROR: fread return a number of element different from the expected.\n");
 422                         return 1;
 423                 }
 424                 fclose(fsrc);
 425
 426                 /* decode the code-stream */
 427                 /* ---------------------- */
 428
 429                 switch(parameters.decod_format) {
 430                 case J2K_CFMT:
 431                 {
 432                         /* JPEG-2000 codestream */
 433
 434                         /* get a decoder handle */
 435                         dinfo = opj_create_decompress(CODEC_J2K);
 436
 437                         /* catch events using our callbacks and give a local context */
 438                         opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);
 439
 440                         /* setup the decoder decoding parameters using user parameters */
 441                         opj_setup_decoder(dinfo, &parameters);
 442
 443                         /* open a byte stream */
 444                         cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length);
 445
 446                         /* decode the stream and fill the image structure */
 447                         if (*indexfilename)                             // If need to extract codestream information
 448                                 image = opj_decode_with_info(dinfo, cio, &cstr_info);
 449                         else
 450                                 image = opj_decode(dinfo, cio);
 451                         if(!image) {
 452                                 fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n");
 453                                 opj_destroy_decompress(dinfo);
 454                                 opj_cio_close(cio);
 455                                 return 1;
 456                         }
 457                         /* dump image */
 458       j2k_dump_image(fout, image);
 459
 460                         /* dump cp */
 461       j2k_dump_cp(fout, image, ((opj_j2k_t*)dinfo->j2k_handle)->cp);
 462
 463                         /* close the byte stream */
 464                         opj_cio_close(cio);
 465
 466                         /* Write the index to disk */
 467                         if (*indexfilename) {
 468                                 opj_bool bSuccess;
 469                                 bSuccess = write_index_file(&cstr_info, indexfilename);
 470                                 if (bSuccess) {
 471                                         fprintf(stderr, "Failed to output index file\n");
 472                                 }
 473                         }
 474                 }
 475                 break;
 476
 477                 case JP2_CFMT:
 478                 {
 479                         /* JPEG 2000 compressed image data */
 480
 481                         /* get a decoder handle */
 482                         dinfo = opj_create_decompress(CODEC_JP2);
 483
 484                         /* catch events using our callbacks and give a local context */
 485                         opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);
 486
 487                         /* setup the decoder decoding parameters using the current image and user parameters */
 488                         opj_setup_decoder(dinfo, &parameters);
 489
 490                         /* open a byte stream */
 491                         cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length);
 492
 493                         /* decode the stream and fill the image structure */
 494                         if (*indexfilename)                             // If need to extract codestream information
 495                                 image = opj_decode_with_info(dinfo, cio, &cstr_info);
 496                         else
 497                                 image = opj_decode(dinfo, cio);
 498                         if(!image) {
 499                                 fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n");
 500                                 opj_destroy_decompress(dinfo);
 501                                 opj_cio_close(cio);
 502                                 return 1;
 503                         }
 504                         /* dump image */
 505           if(image->icc_profile_buf)
 506          {
 507           free(image->icc_profile_buf); image->icc_profile_buf = NULL;
 508          }
 509       j2k_dump_image(fout, image);
 510
 511                         /* dump cp */
 512       j2k_dump_cp(fout, image, ((opj_jp2_t*)dinfo->jp2_handle)->j2k->cp);
 513
 514                         /* close the byte stream */
 515                         opj_cio_close(cio);
 516
 517                         /* Write the index to disk */
 518                         if (*indexfilename) {
 519                                 opj_bool bSuccess;
 520                                 bSuccess = write_index_file(&cstr_info, indexfilename);
 521                                 if (bSuccess) {
 522                                         fprintf(stderr, "Failed to output index file\n");
 523                                 }
 524                         }
 525                 }
 526                 break;
 527
 528                 case JPT_CFMT:
 529                 {
 530                         /* JPEG 2000, JPIP */
 531
 532                         /* get a decoder handle */
 533                         dinfo = opj_create_decompress(CODEC_JPT);
 534
 535                         /* catch events using our callbacks and give a local context */
 536                         opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);
 537
 538                         /* setup the decoder decoding parameters using user parameters */
 539                         opj_setup_decoder(dinfo, &parameters);
 540
 541                         /* open a byte stream */
 542                         cio = opj_cio_open((opj_common_ptr)dinfo, src, file_length);
 543
 544                         /* decode the stream and fill the image structure */
 545                         if (*indexfilename)                             // If need to extract codestream information
 546                                 image = opj_decode_with_info(dinfo, cio, &cstr_info);
 547                         else
 548                                 image = opj_decode(dinfo, cio);
 549                         if(!image) {
 550                                 fprintf(stderr, "ERROR -> j2k_to_image: failed to decode image!\n");
 551                                 opj_destroy_decompress(dinfo);
 552                                 opj_cio_close(cio);
 553                                 return 1;
 554                         }
 555
 556                         /* close the byte stream */
 557                         opj_cio_close(cio);
 558
 559                         /* Write the index to disk */
 560                         if (*indexfilename) {
 561                                 opj_bool bSuccess;
 562                                 bSuccess = write_index_file(&cstr_info, indexfilename);
 563                                 if (bSuccess) {
 564                                         fprintf(stderr, "Failed to output index file\n");
 565                                 }
 566                         }
 567                 }
 568                 break;
 569
 570                 default:
 571                         fprintf(stderr, "skipping file..\n");
 572                         continue;
 573         }
 574
 575                 /* free the memory containing the code-stream */
 576                 free(src);
 577                 src = NULL;
 578
 579                 /* free remaining structures */
 580                 if(dinfo) {
 581                         opj_destroy_decompress(dinfo);
 582                 }
 583                 /* free codestream information structure */
 584                 if (*indexfilename)
 585                         opj_destroy_cstr_info(&cstr_info);
 586                 /* free image data structure */
 587                 opj_image_destroy(image);
 588
 589         }
 590
 591         fclose(fout);
 592
 593   return EXIT_SUCCESS;
 594 }
 595
 596
 597 static void j2k_dump_image(FILE *fd, opj_image_t * img) {
 598         int compno;
 599         fprintf(fd, "image {\n");
 600         fprintf(fd, "  x0=%d, y0=%d, x1=%d, y1=%d\n", img->x0, img->y0, img->x1, img->y1);
 601         fprintf(fd, "  numcomps=%d\n", img->numcomps);
 602         for (compno = 0; compno < img->numcomps; compno++) {
 603                 opj_image_comp_t *comp = &img->comps[compno];
 604                 fprintf(fd, "  comp %d {\n", compno);
 605                 fprintf(fd, "    dx=%d, dy=%d\n", comp->dx, comp->dy);
 606                 fprintf(fd, "    prec=%d\n", comp->prec);
 607                 //fprintf(fd, "    bpp=%d\n", comp->bpp);
 608                 fprintf(fd, "    sgnd=%d\n", comp->sgnd);
 609                 fprintf(fd, "  }\n");
 610         }
 611         fprintf(fd, "}\n");
 612 }
 613
 614 static void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_t * cp) {
 615         int tileno, compno, layno, bandno, resno, numbands;
 616         fprintf(fd, "coding parameters {\n");
 617         fprintf(fd, "  tx0=%d, ty0=%d\n", cp->tx0, cp->ty0);
 618         fprintf(fd, "  tdx=%d, tdy=%d\n", cp->tdx, cp->tdy);
 619         fprintf(fd, "  tw=%d, th=%d\n", cp->tw, cp->th);
 620         for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
 621                 opj_tcp_t *tcp = &cp->tcps[tileno];
 622                 fprintf(fd, "  tile %d {\n", tileno);
 623                 fprintf(fd, "    csty=%x\n", tcp->csty);
 624                 fprintf(fd, "    prg=%d\n", tcp->prg);
 625                 fprintf(fd, "    numlayers=%d\n", tcp->numlayers);
 626                 fprintf(fd, "    mct=%d\n", tcp->mct);
 627                 fprintf(fd, "    rates=");
 628                 for (layno = 0; layno < tcp->numlayers; layno++) {
 629                         fprintf(fd, "%.1f ", tcp->rates[layno]);
 630                 }
 631                 fprintf(fd, "\n");
 632                 for (compno = 0; compno < img->numcomps; compno++) {
 633                         opj_tccp_t *tccp = &tcp->tccps[compno];
 634                         fprintf(fd, "    comp %d {\n", compno);
 635                         fprintf(fd, "      csty=%x\n", tccp->csty);
 636                         fprintf(fd, "      numresolutions=%d\n", tccp->numresolutions);
 637                         fprintf(fd, "      cblkw=%d\n", tccp->cblkw);
 638                         fprintf(fd, "      cblkh=%d\n", tccp->cblkh);
 639                         fprintf(fd, "      cblksty=%x\n", tccp->cblksty);
 640                         fprintf(fd, "      qmfbid=%d\n", tccp->qmfbid);
 641                         fprintf(fd, "      qntsty=%d\n", tccp->qntsty);
 642                         fprintf(fd, "      numgbits=%d\n", tccp->numgbits);
 643                         fprintf(fd, "      roishift=%d\n", tccp->roishift);
 644                         fprintf(fd, "      stepsizes=");
 645                         numbands = tccp->qntsty == J2K_CCP_QNTSTY_SIQNT ? 1 : tccp->numresolutions * 3 - 2;
 646                         for (bandno = 0; bandno < numbands; bandno++) {
 647                                 fprintf(fd, "(%d,%d) ", tccp->stepsizes[bandno].mant,
 648                                         tccp->stepsizes[bandno].expn);
 649                         }
 650                         fprintf(fd, "\n");
 651
 652                         if (tccp->csty & J2K_CCP_CSTY_PRT) {
 653                                 fprintf(fd, "      prcw=");
 654                                 for (resno = 0; resno < tccp->numresolutions; resno++) {
 655                                         fprintf(fd, "%d ", tccp->prcw[resno]);
 656                                 }
 657                                 fprintf(fd, "\n");
 658                                 fprintf(fd, "      prch=");
 659                                 for (resno = 0; resno < tccp->numresolutions; resno++) {
 660                                         fprintf(fd, "%d ", tccp->prch[resno]);
 661                                 }
 662                                 fprintf(fd, "\n");
 663                         }
 664                         fprintf(fd, "    }\n");
 665                 }
 666                 fprintf(fd, "  }\n");
 667         }
 668         fprintf(fd, "}\n");
 669 }
 670