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
  50
  51 int opj_write_ppix( int coff, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_stream_private_t *cio,
  52               opj_event_mgr_t * p_manager )
  53 {
  54   OPJ_BYTE l_data_header [4];
  55   int compno, i;
  56   opj_jp2_box_t *box;
  57   OPJ_OFF_T lenp;
  58   OPJ_UINT32 len;
  59
  60   /*  printf("cstr_info.packno %d\n", cstr_info.packno); //NMAX? */
  61
  62   lenp = -1;
  63   box = (opj_jp2_box_t *)opj_calloc( cstr_info.numcomps, sizeof(opj_jp2_box_t));
  64
  65   for (i=0;i<2;i++){
  66     if (i)
  67
  68       opj_stream_seek( cio, lenp, p_manager);
  69
  70     lenp = (OPJ_UINT32)(opj_stream_tell(cio));
  71     opj_stream_skip( cio, 4, p_manager);       /* L [at the end] */
  72     opj_write_bytes(l_data_header,JPIP_PPIX,4);/* PPIX           */
  73     opj_stream_write_data(cio,l_data_header,4,p_manager);
  74
  75     opj_write_manf( i, cstr_info.numcomps, box, cio, p_manager);
  76
  77     for (compno=0; compno<cstr_info.numcomps; compno++){
  78       box[compno].length = opj_write_ppixfaix( coff, compno, cstr_info, EPHused, j2klen, cio,p_manager);
  79       box[compno].type = JPIP_FAIX;
  80     }
  81
  82
  83   len = (OPJ_UINT32)(opj_stream_tell(cio)-lenp);
  84   opj_stream_seek(cio, lenp,p_manager);
  85   opj_write_bytes(l_data_header,len,4);/* L              */
  86   opj_stream_write_data(cio,l_data_header,4,p_manager);
  87   opj_stream_seek(cio, lenp+len,p_manager);
  88   }
  89
  90   opj_free(box);
  91
  92   return len;
  93 }
  94
  95
  96
  97 int opj_write_ppixfaix( int coff, int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_stream_private_t *cio,
  98               opj_event_mgr_t * p_manager )
  99 {
 100   OPJ_BYTE l_data_header [8];
 101   int tileno, version, i, nmax, size_of_coding; /* 4 or 8*/
 102   OPJ_UINT32 len;
 103   OPJ_OFF_T lenp;
 104   opj_tile_info_t *tile_Idx;
 105   opj_packet_info_t packet;
 106   int resno, precno, layno, num_packet;
 107   int numOfres, numOfprec, numOflayers;
 108   packet.end_pos = packet.end_ph_pos = packet.start_pos = -1;
 109   (void)EPHused; /* unused ? */
 110
 111   if( j2klen > pow( 2, 32)){
 112     size_of_coding =  8;
 113     version = 1;
 114   }
 115   else{
 116     size_of_coding = 4;
 117     version = 0;
 118   }
 119
 120   lenp = opj_stream_tell(cio);
 121   opj_stream_skip(cio, 4, p_manager);         /* L [at the end]      */
 122   opj_write_bytes(l_data_header,JPIP_FAIX,4); /* FAIX */
 123   opj_write_bytes(l_data_header,version,1);
 124   opj_stream_write_data(cio,l_data_header,1,p_manager);/* Version 0 = 4 bytes */
 125
 126   nmax = 0;
 127   for( i=0; i<=cstr_info.numdecompos[compno]; i++)
 128     nmax += cstr_info.tile[0].ph[i] * cstr_info.tile[0].pw[i] * cstr_info.numlayers;
 129
 130   opj_write_bytes(l_data_header,nmax,size_of_coding);         /* NMAX           */
 131   opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 132   opj_write_bytes(l_data_header,cstr_info.tw*cstr_info.th,size_of_coding);  /* M              */
 133   opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 134
 135   for( tileno=0; tileno<cstr_info.tw*cstr_info.th; tileno++){
 136     tile_Idx = &cstr_info.tile[ tileno];
 137
 138     num_packet=0;
 139     numOfres = cstr_info.numdecompos[compno] + 1;
 140
 141     for( resno=0; resno<numOfres ; resno++){
 142       numOfprec = tile_Idx->pw[resno]*tile_Idx->ph[resno];
 143       for( precno=0; precno<numOfprec; precno++){
 144         numOflayers = cstr_info.numlayers;
 145         for( layno=0; layno<numOflayers; layno++){
 146
 147           switch ( cstr_info.prog){
 148           case LRCP:
 149             packet = tile_Idx->packet[ ((layno*numOfres+resno)*cstr_info.numcomps+compno)*numOfprec+precno];
 150             break;
 151           case RLCP:
 152             packet = tile_Idx->packet[ ((resno*numOflayers+layno)*cstr_info.numcomps+compno)*numOfprec+precno];
 153             break;
 154           case RPCL:
 155             packet = tile_Idx->packet[ ((resno*numOfprec+precno)*cstr_info.numcomps+compno)*numOflayers+layno];
 156             break;
 157           case PCRL:
 158             packet = tile_Idx->packet[ ((precno*cstr_info.numcomps+compno)*numOfres+resno)*numOflayers + layno];
 159             break;
 160           case CPRL:
 161             packet = tile_Idx->packet[ ((compno*numOfprec+precno)*numOfres+resno)*numOflayers + layno];
 162             break;
 163           default:
 164             fprintf( stderr, "failed to ppix indexing\n");
 165           }
 166
 167     opj_write_bytes(l_data_header,(OPJ_UINT32)(packet.start_pos-coff),size_of_coding);            /* start position */
 168     opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 169     opj_write_bytes(l_data_header,(OPJ_UINT32)(packet.end_pos-packet.start_pos+1),size_of_coding); /* length         */
 170     opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 171
 172           num_packet++;
 173         }
 174       }
 175     }
 176
 177     while( num_packet < nmax){     /* PADDING */
 178       opj_write_bytes(l_data_header,0,size_of_coding);/* start position            */
 179       opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 180       opj_write_bytes(l_data_header,0,size_of_coding);/* length                    */
 181       opj_stream_write_data(cio,l_data_header,size_of_coding,p_manager);
 182       num_packet++;
 183     }
 184   }
 185
 186   len = (OPJ_UINT32)(opj_stream_tell(cio)-lenp);
 187   opj_stream_seek(cio, lenp,p_manager);
 188   opj_write_bytes(l_data_header,len,4);/* L  */
 189   opj_stream_write_data(cio,l_data_header,4,p_manager);
 190   opj_stream_seek(cio, lenp+len,p_manager);
 191
 192   return len;
 193 }