src/lib/openjp2/ppix_manager.c

   1 /*
   2  * $Id: ppix_manager.c 897 2011-08-28 21:43:57Z Kaori.Hagihara@gmail.com $
   3  *
   4  * Copyright (c) 2002-2011, Communications and Remote Sensing Laboratory, Universite catholique de Louvain (UCL), Belgium
   5  * Copyright (c) 2002-2011, Professor Benoit Macq
   6  * Copyright (c) 2003-2004, Yannick Verschueren
   7  * Copyright (c) 2010-2011, Kaori Hagihara
   8  * All rights reserved.
   9  *
  10  * Redistribution and use in source and binary forms, with or without
  11  * modification, are permitted provided that the following conditions
  12  * are met:
  13  * 1. Redistributions of source code must retain the above copyright
  14  *    notice, this list of conditions and the following disclaimer.
  15  * 2. Redistributions in binary form must reproduce the above copyright
  16  *    notice, this list of conditions and the following disclaimer in the
  17  *    documentation and/or other materials provided with the distribution.
  18  *
  19  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
  20  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  21  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  22  * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  23  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  24  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  25  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  26  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  27  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  28  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  29  * POSSIBILITY OF SUCH DAMAGE.
  30  */
  31
  32 /*! \file
  33  *  \brief Modification of jpip.c from 2KAN indexer
  34  */
  35
  36 #include "opj_includes.h"
  37
  38 /*
  39  * Write faix box of ppix
  40  *
  41  * @param[in] coff offset of j2k codestream
  42  * @param[in] compno    component number
  43  * @param[in] cstr_info codestream information
  44  * @param[in] EPHused   true if if EPH option used
  45  * @param[in] j2klen    length of j2k codestream
  46  * @param[in] cio       file output handle
  47  * @return              length of faix box
  48  */
  49 int write_ppixfaix( int coff, int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio);
  50
  51 int write_ppix( int coff, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio)
  52 {
  53   int len, lenp, compno, i;
  54   opj_jp2_box_t *box;
  55
  56   /*  printf("cstr_info.packno %d\n", cstr_info.packno); //NMAX? */
  57
  58   lenp = -1;
  59   box = (opj_jp2_box_t *)opj_calloc( cstr_info.numcomps, sizeof(opj_jp2_box_t));
  60
  61   for (i=0;i<2;i++){
  62     if (i) cio_seek( cio, lenp);
  63
  64     lenp = cio_tell( cio);
  65     cio_skip( cio, 4);              /* L [at the end] */
  66     cio_write( cio, JPIP_PPIX, 4);  /* PPIX           */
  67
  68     write_manf( i, cstr_info.numcomps, box, cio);
  69
  70     for (compno=0; compno<cstr_info.numcomps; compno++){
  71       box[compno].length = write_ppixfaix( coff, compno, cstr_info, EPHused, j2klen, cio);
  72       box[compno].type = JPIP_FAIX;
  73     }
  74
  75     len = cio_tell( cio)-lenp;
  76     cio_seek( cio, lenp);
  77     cio_write( cio, len, 4);        /* L              */
  78     cio_seek( cio, lenp+len);
  79   }
  80
  81   opj_free(box);
  82
  83   return len;
  84 }
  85
  86 int write_ppix_v2( int coff, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_stream_private_t *cio,
  87               opj_event_mgr_t * p_manager )
  88 {
  89   OPJ_BYTE l_data_header [4];
  90   int len, lenp, compno, i;
  91   opj_jp2_box_t *box;
  92
  93   /*  printf("cstr_info.packno %d\n", cstr_info.packno); //NMAX? */
  94
  95   lenp = -1;
  96   box = (opj_jp2_box_t *)opj_calloc( cstr_info.numcomps, sizeof(opj_jp2_box_t));
  97
  98   for (i=0;i<2;i++){
  99     if (i)
 100 #if 0
 101       cio_seek( cio, lenp);
 102 #else
 103       opj_stream_seek( cio, lenp, p_manager);
 104 #endif
 105
 106 #if 0
 107     lenp = cio_tell( cio);
 108     cio_skip( cio, 4);              /* L [at the end] */
 109     cio_write( cio, JPIP_PPIX, 4);  /* PPIX           */
 110 #else
 111     lenp = opj_stream_tell(cio);
 112     opj_stream_skip( cio, 4, p_manager);       /* L [at the end] */
 113     opj_write_bytes(l_data_header,JPIP_PPIX,4);/* PPIX           */
 114     opj_stream_write_data(cio,l_data_header,4,p_manager);
 115 #endif
 116
 117     write_manf_v2( i, cstr_info.numcomps, box, cio);
 118
 119     for (compno=0; compno<cstr_info.numcomps; compno++){
 120       box[compno].length = write_ppixfaix_v2( coff, compno, cstr_info, EPHused, j2klen, cio,p_manager);
 121       box[compno].type = JPIP_FAIX;
 122     }
 123
 124 #if 0
 125     len = cio_tell( cio)-lenp;
 126     cio_seek( cio, lenp);
 127     cio_write( cio, len, 4);        /* L              */
 128     cio_seek( cio, lenp+len);
 129 #else
 130   len = opj_stream_tell(cio)-lenp;
 131   opj_stream_seek(cio, lenp,p_manager);
 132   opj_write_bytes(l_data_header,len,4);/* L              */
 133   opj_stream_write_data(cio,l_data_header,4,p_manager);
 134   opj_stream_seek(cio, lenp+len,p_manager);
 135 #endif
 136   }
 137
 138   opj_free(box);
 139
 140   return len;
 141 }
 142
 143 int write_ppixfaix( int coff, int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio)
 144 {
 145   int len, lenp, tileno, version, i, nmax, size_of_coding; /* 4 or 8*/
 146   opj_tile_info_t *tile_Idx;
 147   opj_packet_info_t packet;
 148   int resno, precno, layno, num_packet;
 149   int numOfres, numOfprec, numOflayers;
 150   packet.end_pos = packet.end_ph_pos = packet.start_pos = -1;
 151   (void)EPHused; /* unused ? */
 152
 153   if( j2klen > pow( 2, 32)){
 154     size_of_coding =  8;
 155     version = 1;
 156   }
 157   else{
 158     size_of_coding = 4;
 159     version = 0;
 160   }
 161
 162   lenp = cio_tell( cio);
 163   cio_skip( cio, 4);              /* L [at the end]      */
 164   cio_write( cio, JPIP_FAIX, 4);  /* FAIX                */
 165   cio_write( cio, version, 1);     /* Version 0 = 4 bytes */
 166
 167   nmax = 0;
 168   for( i=0; i<=cstr_info.numdecompos[compno]; i++)
 169     nmax += cstr_info.tile[0].ph[i] * cstr_info.tile[0].pw[i] * cstr_info.numlayers;
 170
 171   cio_write( cio, nmax, size_of_coding); /* NMAX */
 172   cio_write( cio, cstr_info.tw*cstr_info.th, size_of_coding);      /* M    */
 173
 174   for( tileno=0; tileno<cstr_info.tw*cstr_info.th; tileno++){
 175     tile_Idx = &cstr_info.tile[ tileno];
 176
 177     num_packet=0;
 178     numOfres = cstr_info.numdecompos[compno] + 1;
 179
 180     for( resno=0; resno<numOfres ; resno++){
 181       numOfprec = tile_Idx->pw[resno]*tile_Idx->ph[resno];
 182       for( precno=0; precno<numOfprec; precno++){
 183         numOflayers = cstr_info.numlayers;
 184         for( layno=0; layno<numOflayers; layno++){
 185
 186           switch ( cstr_info.prog){
 187           case LRCP:
 188             packet = tile_Idx->packet[ ((layno*numOfres+resno)*cstr_info.numcomps+compno)*numOfprec+precno];
 189             break;
 190           case RLCP:
 191             packet = tile_Idx->packet[ ((resno*numOflayers+layno)*cstr_info.numcomps+compno)*numOfprec+precno];
 192             break;
 193           case RPCL:
 194             packet = tile_Idx->packet[ ((resno*numOfprec+precno)*cstr_info.numcomps+compno)*numOflayers+layno];
 195             break;
 196           case PCRL:
 197             packet = tile_Idx->packet[ ((precno*cstr_info.numcomps+compno)*numOfres+resno)*numOflayers + layno];
 198             break;
 199           case CPRL:
 200             packet = tile_Idx->packet[ ((compno*numOfprec+precno)*numOfres+resno)*numOflayers + layno];
 201             break;
 202           default:
 203             fprintf( stderr, "failed to ppix indexing\n");
 204           }
 205
 206           cio_write( cio, packet.start_pos-coff, size_of_coding);             /* start position */
 207           cio_write( cio, packet.end_pos-packet.start_pos+1, size_of_coding); /* length         */
 208
 209           num_packet++;
 210         }
 211       }
 212     }
 213
 214     while( num_packet < nmax){     /* PADDING */
 215       cio_write( cio, 0, size_of_coding); /* start position            */
 216       cio_write( cio, 0, size_of_coding); /* length                    */
 217       num_packet++;
 218     }
 219   }
 220
 221   len = cio_tell( cio)-lenp;
 222   cio_seek( cio, lenp);
 223   cio_write( cio, len, 4);        /* L  */
 224   cio_seek( cio, lenp+len);
 225
 226   return len;
 227 }
 228
 229 int write_ppixfaix_v2( int coff, int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_stream_private_t *cio,
 230               opj_event_mgr_t * p_manager )
 231 {
 232   OPJ_BYTE l_data_header [8];
 233   int len, lenp, tileno, version, i, nmax, size_of_coding; /* 4 or 8*/
 234   opj_tile_info_t *tile_Idx;
 235   opj_packet_info_t packet;
 236   int resno, precno, layno, num_packet;
 237   int numOfres, numOfprec, numOflayers;
 238   packet.end_pos = packet.end_ph_pos = packet.start_pos = -1;
 239   (void)EPHused; /* unused ? */
 240
 241   if( j2klen > pow( 2, 32)){
 242     size_of_coding =  8;
 243     version = 1;
 244   }
 245   else{
 246     size_of_coding = 4;
 247     version = 0;
 248   }
 249
 250 #if 0
 251   lenp = cio_tell( cio);
 252   cio_skip( cio, 4);              /* L [at the end]      */
 253   cio_write( cio, JPIP_FAIX, 4);  /* FAIX                */
 254   cio_write( cio, version, 1);     /* Version 0 = 4 bytes */
 255 #else
 256   lenp = opj_stream_tell(cio);
 257   opj_stream_skip(cio, 4, p_manager);         /* L [at the end]      */
 258   opj_write_bytes(l_data_header,JPIP_FAIX,4); /* FAIX */
 259   opj_write_bytes(l_data_header,version,1);
 260   opj_stream_write_data(cio,l_data_header,1,p_manager);/* Version 0 = 4 bytes */
 261 #endif
 262
 263   nmax = 0;
 264   for( i=0; i<=cstr_info.numdecompos[compno]; i++)
 265     nmax += cstr_info.tile[0].ph[i] * cstr_info.tile[0].pw[i] * cstr_info.numlayers;
 266
 267 #if 0
 268   cio_write( cio, nmax, size_of_coding); /* NMAX */
 269   cio_write( cio, cstr_info.tw*cstr_info.th, size_of_coding);      /* M    */
 270 #else
 271   opj_write_bytes(l_data_header,nmax,size_of_coding);         /* NMAX           */
 272   opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 273   opj_write_bytes(l_data_header,cstr_info.tw*cstr_info.th,size_of_coding);  /* M              */
 274   opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 275 #endif
 276
 277   for( tileno=0; tileno<cstr_info.tw*cstr_info.th; tileno++){
 278     tile_Idx = &cstr_info.tile[ tileno];
 279
 280     num_packet=0;
 281     numOfres = cstr_info.numdecompos[compno] + 1;
 282
 283     for( resno=0; resno<numOfres ; resno++){
 284       numOfprec = tile_Idx->pw[resno]*tile_Idx->ph[resno];
 285       for( precno=0; precno<numOfprec; precno++){
 286         numOflayers = cstr_info.numlayers;
 287         for( layno=0; layno<numOflayers; layno++){
 288
 289           switch ( cstr_info.prog){
 290           case LRCP:
 291             packet = tile_Idx->packet[ ((layno*numOfres+resno)*cstr_info.numcomps+compno)*numOfprec+precno];
 292             break;
 293           case RLCP:
 294             packet = tile_Idx->packet[ ((resno*numOflayers+layno)*cstr_info.numcomps+compno)*numOfprec+precno];
 295             break;
 296           case RPCL:
 297             packet = tile_Idx->packet[ ((resno*numOfprec+precno)*cstr_info.numcomps+compno)*numOflayers+layno];
 298             break;
 299           case PCRL:
 300             packet = tile_Idx->packet[ ((precno*cstr_info.numcomps+compno)*numOfres+resno)*numOflayers + layno];
 301             break;
 302           case CPRL:
 303             packet = tile_Idx->packet[ ((compno*numOfprec+precno)*numOfres+resno)*numOflayers + layno];
 304             break;
 305           default:
 306             fprintf( stderr, "failed to ppix indexing\n");
 307           }
 308
 309 #if 0
 310           cio_write( cio, packet.start_pos-coff, size_of_coding);             /* start position */
 311           cio_write( cio, packet.end_pos-packet.start_pos+1, size_of_coding); /* length         */
 312 #else
 313     opj_write_bytes(l_data_header,packet.start_pos-coff,size_of_coding);            /* start position */
 314     opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 315     opj_write_bytes(l_data_header,packet.end_pos-packet.start_pos+1,size_of_coding); /* length         */
 316     opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 317 #endif
 318
 319           num_packet++;
 320         }
 321       }
 322     }
 323
 324     while( num_packet < nmax){     /* PADDING */
 325 #if 0
 326       cio_write( cio, 0, size_of_coding); /* start position            */
 327       cio_write( cio, 0, size_of_coding); /* length                    */
 328 #else
 329       opj_write_bytes(l_data_header,0,size_of_coding);/* start position            */
 330       opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 331       opj_write_bytes(l_data_header,0,size_of_coding);/* length                    */
 332       opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 333 #endif
 334       num_packet++;
 335     }
 336   }
 337
 338 #if 0
 339   len = cio_tell( cio)-lenp;
 340   cio_seek( cio, lenp);
 341   cio_write( cio, len, 4);        /* L  */
 342   cio_seek( cio, lenp+len);
 343 #else
 344   len = opj_stream_tell(cio)-lenp;
 345   opj_stream_seek(cio, lenp,p_manager);
 346   opj_write_bytes(l_data_header,len,4);/* L  */
 347   opj_stream_write_data(cio,l_data_header,4,p_manager);
 348   opj_stream_seek(cio, lenp+len,p_manager);
 349 #endif
 350
 351   return len;
 352 }