src/lib/openjpip/phix_manager.c

   1 /*
   2  * $Id: phix_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 /*
  40  * Write faix box of phix
  41  *
  42  * @param[in] coff      offset of j2k codestream
  43  * @param[in] compno    component number
  44  * @param[in] cstr_info codestream information
  45  * @param[in] EPHused   true if if EPH option used
  46  * @param[in] j2klen    length of j2k codestream
  47  * @param[in] cio       file output handle
  48  * @return              length of faix box
  49  */
  50 int write_phixfaix( int coff, int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio);
  51
  52 int write_phix( int coff, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio)
  53 {
  54   int len, lenp=0, compno, i;
  55   opj_jp2_box_t *box;
  56
  57   box = (opj_jp2_box_t *)opj_calloc( cstr_info.numcomps, sizeof(opj_jp2_box_t));
  58
  59   for( i=0;i<2;i++){
  60     if (i) cio_seek( cio, lenp);
  61
  62     lenp = cio_tell( cio);
  63     cio_skip( cio, 4);              /* L [at the end] */
  64     cio_write( cio, JPIP_PHIX, 4);  /* PHIX           */
  65
  66     write_manf( i, cstr_info.numcomps, box, cio);
  67
  68     for( compno=0; compno<cstr_info.numcomps; compno++){
  69       box[compno].length = write_phixfaix( coff, compno, cstr_info, EPHused, j2klen, cio);
  70       box[compno].type = JPIP_FAIX;
  71     }
  72
  73     len = cio_tell( cio)-lenp;
  74     cio_seek( cio, lenp);
  75     cio_write( cio, len, 4);        /* L              */
  76     cio_seek( cio, lenp+len);
  77   }
  78
  79   opj_free(box);
  80
  81   return len;
  82 }
  83
  84 int write_phix_v2( int coff, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_stream_private_t *cio,
  85               opj_event_mgr_t * p_manager )
  86 {
  87   OPJ_BYTE l_data_header [8];
  88   int len, lenp=0, compno, i;
  89   opj_jp2_box_t *box;
  90
  91   box = (opj_jp2_box_t *)opj_calloc( cstr_info.numcomps, sizeof(opj_jp2_box_t));
  92
  93   for( i=0;i<2;i++){
  94     if (i)
  95 #if 0
  96       cio_seek( cio, lenp);
  97 #else
  98       opj_stream_seek( cio, lenp, p_manager);
  99 #endif
 100
 101 #if 0
 102     lenp = cio_tell( cio);
 103     cio_skip( cio, 4);              /* L [at the end] */
 104     cio_write( cio, JPIP_PHIX, 4);  /* PHIX           */
 105 #else
 106     lenp = opj_stream_tell(cio);
 107     opj_stream_skip(cio, 4, p_manager);         /* L [at the end]      */
 108     opj_write_bytes(l_data_header,JPIP_PHIX,4); /* PHIX */
 109     opj_stream_write_data(cio,l_data_header,4,p_manager);
 110 #endif
 111
 112     write_manf_v2( i, cstr_info.numcomps, box, cio);
 113
 114     for( compno=0; compno<cstr_info.numcomps; compno++){
 115       box[compno].length = write_phixfaix_v2( coff, compno, cstr_info, EPHused, j2klen, cio);
 116       box[compno].type = JPIP_FAIX;
 117     }
 118
 119 #if 0
 120     len = cio_tell( cio)-lenp;
 121     cio_seek( cio, lenp);
 122     cio_write( cio, len, 4);        /* L              */
 123     cio_seek( cio, lenp+len);
 124 #else
 125     lenp = opj_stream_tell(cio);
 126     opj_stream_seek(cio, 4, p_manager);
 127     opj_write_bytes(l_data_header,len,4);/* L              */
 128     opj_stream_write_data(cio,l_data_header,4,p_manager);
 129     opj_stream_seek( cio, lenp+len,p_manager);
 130 #endif
 131   }
 132
 133   opj_free(box);
 134
 135   return len;
 136 }
 137
 138 int write_phixfaix( int coff, int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio)
 139 {
 140   int len, lenp, tileno, version, i, nmax, size_of_coding; /* 4 or 8 */
 141   opj_tile_info_t *tile_Idx;
 142   opj_packet_info_t packet;
 143   int resno, precno, layno, num_packet;
 144   int numOfres, numOfprec, numOflayers;
 145   packet.end_ph_pos = packet.start_pos = -1;
 146   (void)EPHused; /* unused ? */
 147
 148   if( j2klen > pow( 2, 32)){
 149     size_of_coding =  8;
 150     version = 1;
 151   }
 152   else{
 153     size_of_coding = 4;
 154     version = 0;
 155   }
 156
 157   lenp = cio_tell( cio);
 158   cio_skip( cio, 4);              /* L [at the end]      */
 159   cio_write( cio, JPIP_FAIX, 4);  /* FAIX                */
 160   cio_write( cio, version,1);     /* Version 0 = 4 bytes */
 161
 162   nmax = 0;
 163   for( i=0; i<=cstr_info.numdecompos[compno]; i++)
 164     nmax += cstr_info.tile[0].ph[i] * cstr_info.tile[0].pw[i] * cstr_info.numlayers;
 165
 166   cio_write( cio, nmax, size_of_coding); /* NMAX */
 167   cio_write( cio, cstr_info.tw*cstr_info.th, size_of_coding);      /* M    */
 168
 169   for( tileno=0; tileno<cstr_info.tw*cstr_info.th; tileno++){
 170     tile_Idx = &cstr_info.tile[ tileno];
 171
 172     num_packet = 0;
 173     numOfres = cstr_info.numdecompos[compno] + 1;
 174
 175     for( resno=0; resno<numOfres ; resno++){
 176       numOfprec = tile_Idx->pw[resno]*tile_Idx->ph[resno];
 177       for( precno=0; precno<numOfprec; precno++){
 178         numOflayers = cstr_info.numlayers;
 179         for( layno=0; layno<numOflayers; layno++){
 180
 181           switch ( cstr_info.prog){
 182           case LRCP:
 183             packet = tile_Idx->packet[ ((layno*numOfres+resno)*cstr_info.numcomps+compno)*numOfprec+precno];
 184             break;
 185           case RLCP:
 186             packet = tile_Idx->packet[ ((resno*numOflayers+layno)*cstr_info.numcomps+compno)*numOfprec+precno];
 187             break;
 188           case RPCL:
 189             packet = tile_Idx->packet[ ((resno*numOfprec+precno)*cstr_info.numcomps+compno)*numOflayers+layno];
 190             break;
 191           case PCRL:
 192             packet = tile_Idx->packet[ ((precno*cstr_info.numcomps+compno)*numOfres+resno)*numOflayers + layno];
 193             break;
 194           case CPRL:
 195             packet = tile_Idx->packet[ ((compno*numOfprec+precno)*numOfres+resno)*numOflayers + layno];
 196             break;
 197           default:
 198             fprintf( stderr, "failed to ppix indexing\n");
 199           }
 200
 201           cio_write( cio, packet.start_pos-coff, size_of_coding);                /* start position */
 202           cio_write( cio, packet.end_ph_pos-packet.start_pos+1, size_of_coding); /* length         */
 203
 204           num_packet++;
 205         }
 206       }
 207     }
 208
 209     /* PADDING */
 210     while( num_packet < nmax){
 211       cio_write( cio, 0, size_of_coding); /* start position            */
 212       cio_write( cio, 0, size_of_coding); /* length                    */
 213       num_packet++;
 214     }
 215   }
 216
 217   len = cio_tell( cio)-lenp;
 218   cio_seek( cio, lenp);
 219   cio_write( cio, len, 4);        /* L  */
 220   cio_seek( cio, lenp+len);
 221
 222   return len;
 223 }
 224
 225 int write_phixfaix_v2( int coff, int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_stream_private_t *cio,
 226               opj_event_mgr_t * p_manager )
 227 {
 228   int len, lenp, tileno, version, i, nmax, size_of_coding; /* 4 or 8 */
 229   opj_tile_info_t *tile_Idx;
 230   opj_packet_info_t packet;
 231   int resno, precno, layno, num_packet;
 232   int numOfres, numOfprec, numOflayers;
 233   OPJ_BYTE l_data_header [8];
 234   packet.end_ph_pos = packet.start_pos = -1;
 235   (void)EPHused; /* unused ? */
 236
 237   if( j2klen > pow( 2, 32)){
 238     size_of_coding =  8;
 239     version = 1;
 240   }
 241   else{
 242     size_of_coding = 4;
 243     version = 0;
 244   }
 245
 246 #if 0
 247   lenp = cio_tell( cio);
 248   cio_skip( cio, 4);              /* L [at the end]      */
 249   cio_write( cio, JPIP_FAIX, 4);  /* FAIX                */
 250   cio_write( cio, version,1);     /* Version 0 = 4 bytes */
 251 #else
 252   lenp = opj_stream_tell(cio);
 253   opj_stream_skip(cio, 4, p_manager);         /* L [at the end]      */
 254   opj_write_bytes(l_data_header,JPIP_FAIX,4); /* FAIX */
 255   opj_stream_write_data(cio,l_data_header,4,p_manager);
 256   opj_write_bytes(l_data_header,version,1);   /* Version 0 = 4 bytes */
 257   opj_stream_write_data(cio,l_data_header,1,p_manager);
 258 #endif
 259
 260   nmax = 0;
 261   for( i=0; i<=cstr_info.numdecompos[compno]; i++)
 262     nmax += cstr_info.tile[0].ph[i] * cstr_info.tile[0].pw[i] * cstr_info.numlayers;
 263
 264 #if 0
 265   cio_write( cio, nmax, size_of_coding); /* NMAX */
 266   cio_write( cio, cstr_info.tw*cstr_info.th, size_of_coding);      /* M    */
 267 #else
 268   opj_write_bytes(l_data_header,nmax,size_of_coding);         /* NMAX           */
 269   opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 270   opj_write_bytes(l_data_header,cstr_info.tw*cstr_info.th,size_of_coding);  /* M              */
 271   opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 272 #endif
 273
 274   for( tileno=0; tileno<cstr_info.tw*cstr_info.th; tileno++){
 275     tile_Idx = &cstr_info.tile[ tileno];
 276
 277     num_packet = 0;
 278     numOfres = cstr_info.numdecompos[compno] + 1;
 279
 280     for( resno=0; resno<numOfres ; resno++){
 281       numOfprec = tile_Idx->pw[resno]*tile_Idx->ph[resno];
 282       for( precno=0; precno<numOfprec; precno++){
 283         numOflayers = cstr_info.numlayers;
 284         for( layno=0; layno<numOflayers; layno++){
 285
 286           switch ( cstr_info.prog){
 287           case LRCP:
 288             packet = tile_Idx->packet[ ((layno*numOfres+resno)*cstr_info.numcomps+compno)*numOfprec+precno];
 289             break;
 290           case RLCP:
 291             packet = tile_Idx->packet[ ((resno*numOflayers+layno)*cstr_info.numcomps+compno)*numOfprec+precno];
 292             break;
 293           case RPCL:
 294             packet = tile_Idx->packet[ ((resno*numOfprec+precno)*cstr_info.numcomps+compno)*numOflayers+layno];
 295             break;
 296           case PCRL:
 297             packet = tile_Idx->packet[ ((precno*cstr_info.numcomps+compno)*numOfres+resno)*numOflayers + layno];
 298             break;
 299           case CPRL:
 300             packet = tile_Idx->packet[ ((compno*numOfprec+precno)*numOfres+resno)*numOflayers + layno];
 301             break;
 302           default:
 303             fprintf( stderr, "failed to ppix indexing\n");
 304           }
 305
 306 #if 0
 307           cio_write( cio, packet.start_pos-coff, size_of_coding);                /* start position */
 308           cio_write( cio, packet.end_ph_pos-packet.start_pos+1, size_of_coding); /* length         */
 309 #else
 310     opj_write_bytes(l_data_header,packet.start_pos-coff,size_of_coding);            /* start position */
 311     opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 312     opj_write_bytes(l_data_header,packet.end_ph_pos-packet.start_pos+1,size_of_coding); /* length         */
 313     opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 314 #endif
 315
 316           num_packet++;
 317         }
 318       }
 319     }
 320
 321     /* PADDING */
 322     while( num_packet < nmax){
 323 #if 0
 324       cio_write( cio, 0, size_of_coding); /* start position            */
 325       cio_write( cio, 0, size_of_coding); /* length                    */
 326 #else
 327       opj_write_bytes(l_data_header,0,size_of_coding);/* start position            */
 328       opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 329       opj_write_bytes(l_data_header,0,size_of_coding);/* length                    */
 330       opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 331 #endif
 332       num_packet++;
 333     }
 334   }
 335
 336 #if 0
 337   len = cio_tell( cio)-lenp;
 338   cio_seek( cio, lenp);
 339   cio_write( cio, len, 4);        /* L  */
 340   cio_seek( cio, lenp+len);
 341 #else
 342   len = opj_stream_tell(cio)-lenp;
 343   opj_stream_seek(cio, lenp,p_manager);
 344   opj_write_bytes(l_data_header,len,4);/* L  */
 345   opj_stream_write_data(cio,l_data_header,4,p_manager);
 346   opj_stream_seek(cio, lenp+len,p_manager);
 347 #endif
 348
 349   return len;
 350 }