applications/jpip/tools/indexer/ext_libopenjpeg/phix_manager.c

   1 /*
   2  * $Id$
   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 <stdlib.h>
  37 #include <stdio.h>
  38 #include <math.h>
  39 #include "indexbox_manager.h"
  40 #include "cio_ext.h"
  41 #include "cio.h"
  42
  43 /*
  44  * Write faix box of phix
  45  *
  46  * @param[in] compno    component number
  47  * @param[in] cstr_info codestream information
  48  * @param[in] EPHused   true if if EPH option used
  49  * @param[in] j2klen    length of j2k codestream
  50  * @param[in] cio       file output handle
  51  * @return              length of faix box
  52  */
  53 int write_phixfaix( int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio);
  54
  55 int write_phix( opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio)
  56 {
  57   int len, lenp=0, compno, i;
  58   opj_jp2_box_t *box;
  59
  60   box = (opj_jp2_box_t *)calloc( cstr_info.numcomps, sizeof(opj_jp2_box_t));
  61
  62   for( i=0;i<2;i++){
  63     if (i) cio_seek( cio, lenp);
  64
  65     lenp = cio_tell( cio);
  66     cio_skip( cio, 4);              /* L [at the end] */
  67     cio_write( cio, JPIP_PHIX, 4);  /* PHIX           */
  68
  69     write_manf( i, cstr_info.numcomps, box, cio);
  70
  71     for( compno=0; compno<cstr_info.numcomps; compno++){
  72       box[compno].length = write_phixfaix( compno, cstr_info, EPHused, j2klen, cio);
  73       box[compno].type = JPIP_FAIX;
  74     }
  75
  76     len = cio_tell( cio)-lenp;
  77     cio_seek( cio, lenp);
  78     cio_write( cio, len, 4);        /* L              */
  79     cio_seek( cio, lenp+len);
  80   }
  81
  82   free(box);
  83
  84   return len;
  85 }
  86
  87 int write_phixfaix( int compno, opj_codestream_info_t cstr_info, opj_bool EPHused, int j2klen, opj_cio_t *cio)
  88 {
  89   int len, lenp;
  90   int size_of_coding; // 4 or 8
  91   int version;
  92   int tileno, resno, precno, layno, num_packet=0;
  93   int i,nmax=0;
  94
  95   if( j2klen > pow( 2, 32)){
  96     size_of_coding =  8;
  97     version = 1;
  98   }
  99   else{
 100     size_of_coding = 4;
 101     version = 0;
 102   }
 103
 104   lenp = cio_tell( cio);
 105   cio_skip( cio, 4);              /* L [at the end]      */
 106   cio_write( cio, JPIP_FAIX, 4);  /* FAIX                */
 107   cio_write( cio, version,1);     /* Version 0 = 4 bytes */
 108
 109   for( i=0; i<=cstr_info.numdecompos[compno]; i++)
 110     nmax += cstr_info.tile[0].ph[i] * cstr_info.tile[0].pw[i] * cstr_info.numlayers;
 111
 112   cio_ext_write( cio, nmax , size_of_coding); /* NMAX */
 113   cio_ext_write( cio, cstr_info.tw*cstr_info.th, size_of_coding);      /* M    */
 114
 115   for( tileno=0; tileno<cstr_info.tw*cstr_info.th; tileno++){
 116
 117     opj_tile_info_t *tile_Idx = &cstr_info.tile[ tileno];
 118     //    int correction = EPHused ? 3 : 1;
 119     num_packet = 0;
 120
 121     for( resno=0; resno<cstr_info.numdecompos[compno]+1; resno++){
 122       for( precno=0; precno<tile_Idx->pw[resno]*tile_Idx->ph[resno]; precno++){
 123         for( layno=0; layno<cstr_info.numlayers; layno++){
 124           opj_packet_info_t packet = tile_Idx->packet[ num_packet * cstr_info.numcomps + compno];
 125           cio_ext_write( cio, packet.start_pos, size_of_coding);                                   /* start position */
 126           cio_ext_write( cio, packet.end_ph_pos-packet.start_pos+1, size_of_coding); /* length         */
 127
 128           num_packet++;
 129         }
 130       }
 131     }
 132
 133     /* PADDING */
 134     while( num_packet < nmax){
 135       cio_ext_write( cio, 0, size_of_coding); /* start position            */
 136       cio_ext_write( cio, 0, size_of_coding); /* length                    */
 137       num_packet++;
 138     }
 139   }
 140
 141   len = cio_tell( cio)-lenp;
 142   cio_seek( cio, lenp);
 143   cio_write( cio, len, 4);        /* L  */
 144   cio_seek( cio, lenp+len);
 145
 146   return len;
 147 }