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