086918c77a0fe80743b6e332193ca8152a746182
[openjpeg.git] / src / bin / mj2 / opj_mj2_decompress.c
1 /*
2 * Copyright (c) 2003-2004, Francois-Olivier Devaux
3 * Copyright (c) 2002-2004,  Communications and remote sensing Laboratory, Universite catholique de Louvain, Belgium
4 * All rights reserved.
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 * 1. Redistributions of source code must retain the above copyright
10 *    notice, this list of conditions and the following disclaimer.
11 * 2. Redistributions in binary form must reproduce the above copyright
12 *    notice, this list of conditions and the following disclaimer in the
13 *    documentation and/or other materials provided with the distribution.
14 *
15 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
16 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18 * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
19 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
20 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
21 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
22 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
23 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
24 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
25 * POSSIBILITY OF SUCH DAMAGE.
26 */
27
28 #include <stdio.h>
29 #include <stdlib.h>
30 #include <string.h>
31
32 #include "opj_apps_config.h"
33 #include "openjpeg.h"
34 #include "j2k_lib.h"
35 #include "cio.h"
36 #include "j2k.h"
37 #include "jp2.h"
38 #include "mj2.h"
39 #include "mj2_convert.h"
40
41 #ifdef OPJ_HAVE_LIBLCMS2
42 #include <lcms2.h>
43 #endif
44 #ifdef OPJ_HAVE_LIBLCMS1
45 #include <lcms.h>
46 #endif
47 #include "color.h"
48 /* -------------------------------------------------------------------------- */
49
50 /**
51 sample error callback expecting a FILE* client object
52 */
53 static void error_callback(const char *msg, void *client_data) {
54         FILE *stream = (FILE*)client_data;
55         fprintf(stream, "[ERROR] %s", msg);
56 }
57 /**
58 sample warning callback expecting a FILE* client object
59 */
60 static void warning_callback(const char *msg, void *client_data) {
61         FILE *stream = (FILE*)client_data;
62         fprintf(stream, "[WARNING] %s", msg);
63 }
64
65 /* -------------------------------------------------------------------------- */
66
67
68 int main(int argc, char *argv[]) {
69         mj2_dparameters_t mj2_parameters;                       /* decompression parameters */
70         opj_dinfo_t* dinfo; 
71         opj_event_mgr_t event_mgr;              /* event manager */     
72         opj_cio_t *cio = NULL;
73   unsigned int tnum, snum;
74   opj_mj2_t *movie;
75   mj2_tk_t *track;
76   mj2_sample_t *sample;
77   unsigned char* frame_codestream;
78   FILE *file, *outfile;
79   char outfilename[50];
80   opj_image_t *img = NULL;
81         unsigned int max_codstrm_size = 0;
82         double total_time = 0;
83         unsigned int numframes = 0;
84                         
85   if (argc != 3) {
86     printf("Usage: %s inputfile.mj2 outputfile.yuv\n",argv[0]); 
87     return 1;
88   }
89   
90   file = fopen(argv[1], "rb");
91   
92   if (!file) {
93     fprintf(stderr, "failed to open %s for reading\n", argv[1]);
94     return 1;
95   }
96         
97   /* Checking output file */
98   outfile = fopen(argv[2], "w");
99   if (!file) {
100     fprintf(stderr, "failed to open %s for writing\n", argv[2]);
101     return 1;
102   }
103   fclose(outfile);
104         
105         /*
106         configure the event callbacks (not required)
107         setting of each callback is optionnal
108         */
109         memset(&event_mgr, 0, sizeof(opj_event_mgr_t));
110         event_mgr.error_handler = error_callback;
111         event_mgr.warning_handler = warning_callback;
112         event_mgr.info_handler = NULL;
113         
114         /* get a MJ2 decompressor handle */
115         dinfo = mj2_create_decompress();
116         movie = (opj_mj2_t*)dinfo->mj2_handle;
117         
118         /* catch events using our callbacks and give a local context */
119         opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);           
120
121         memset(&mj2_parameters, 0, sizeof(mj2_dparameters_t));
122         /* set J2K decoding parameters to default values */
123         opj_set_default_decoder_parameters(&mj2_parameters.j2k_parameters);
124         
125         /* setup the decoder decoding parameters using user parameters */
126         mj2_setup_decoder(movie, &mj2_parameters);
127                         
128   if (mj2_read_struct(file, movie)) /* Creating the movie structure */
129     return 1;   
130         
131   /* Decode first video track */
132         for (tnum=0; tnum < (unsigned int)(movie->num_htk + movie->num_stk + movie->num_vtk); tnum++) {
133                 if (movie->tk[tnum].track_type == 0) 
134                         break;
135         }
136         
137         if (movie->tk[tnum].track_type != 0) {
138                 printf("Error. Movie does not contain any video track\n");
139                 return 1;
140         }
141         
142   track = &movie->tk[tnum];
143         
144   /* Output info on first video tracl */
145   fprintf(stdout,"The first video track contains %d frames.\nWidth: %d, Height: %d \n\n",
146     track->num_samples, track->w, track->h);
147         
148         max_codstrm_size = track->sample[0].sample_size-8;
149         frame_codestream = (unsigned char*) malloc(max_codstrm_size * sizeof(unsigned char)); 
150
151         numframes = track->num_samples;
152         
153   for (snum=0; snum < numframes; snum++)
154   {
155                 double init_time = opj_clock();
156                 double elapsed_time;
157
158     sample = &track->sample[snum];
159                 if (sample->sample_size-8 > max_codstrm_size) {
160                         max_codstrm_size =  sample->sample_size-8;
161                         if ((frame_codestream = (unsigned char*)
162                                 realloc(frame_codestream, max_codstrm_size)) == NULL) {
163                                 printf("Error reallocation memory\n");
164                                 return 1;
165                         };              
166                 }
167     fseek(file,sample->offset+8,SEEK_SET);
168     fread(frame_codestream, sample->sample_size-8, 1, file);  /* Assuming that jp and ftyp markers size do */
169                 
170                 /* open a byte stream */
171                 cio = opj_cio_open((opj_common_ptr)dinfo, frame_codestream, sample->sample_size-8);
172                 
173                 img = opj_decode(dinfo, cio); /* Decode J2K to image */
174
175 #ifdef WANT_SYCC_TO_RGB
176         if(img->color_space == CLRSPC_SYCC)
177   {
178         color_sycc_to_rgb(img);
179   }
180 #endif
181
182         if(img->icc_profile_buf)
183   {
184 #if defined(OPJ_HAVE_LIBLCMS1) || defined(OPJ_HAVE_LIBLCMS2)
185         color_apply_icc_profile(img);
186 #endif
187
188         free(img->icc_profile_buf);
189         img->icc_profile_buf = NULL; img->icc_profile_len = 0;
190   }
191
192     if (((img->numcomps == 3) && (img->comps[0].dx == img->comps[1].dx / 2) 
193       && (img->comps[0].dx == img->comps[2].dx / 2 ) && (img->comps[0].dx == 1)) 
194       || (img->numcomps == 1)) {
195       
196       if (!imagetoyuv(img, argv[2]))    /* Convert image to YUV */
197                                 return 1;
198     }
199     else if ((img->numcomps == 3) && 
200       (img->comps[0].dx == 1) && (img->comps[1].dx == 1)&&
201       (img->comps[2].dx == 1))/* If YUV 4:4:4 input --> to bmp */
202     {
203       fprintf(stdout,"The frames will be output in a bmp format (output_1.bmp, ...)\n");
204       sprintf(outfilename,"output_%d.bmp",snum);
205       if (imagetobmp(img, outfilename)) /* Convert image to BMP */
206                                 return 1;
207       
208     }
209     else {
210       fprintf(stdout,"Image component dimensions are unknown. Unable to output image\n");
211       fprintf(stdout,"The frames will be output in a j2k file (output_1.j2k, ...)\n");
212                         
213       sprintf(outfilename,"output_%d.j2k",snum);
214       outfile = fopen(outfilename, "wb");
215       if (!outfile) {
216                                 fprintf(stderr, "failed to open %s for writing\n",outfilename);
217                                 return 1;
218       }
219       fwrite(frame_codestream,sample->sample_size-8,1,outfile);
220       fclose(outfile);
221     }
222                 /* close the byte stream */
223                 opj_cio_close(cio);     
224                 /* free image data structure */
225                 opj_image_destroy(img);
226                 elapsed_time = opj_clock()-init_time;
227                 fprintf(stderr, "Frame number %d/%d decoded in %.2f mseconds\n", snum + 1, numframes, elapsed_time*1000);
228                 total_time += elapsed_time;
229
230   }
231         
232         free(frame_codestream); 
233   fclose(file); 
234
235         /* free remaining structures */
236         if(dinfo) {
237                 mj2_destroy_decompress((opj_mj2_t*)dinfo->mj2_handle);
238         }
239         free(dinfo);
240         
241         fprintf(stdout, "%d frame(s) correctly decompressed\n", snum);
242         fprintf(stdout,"Total decoding time: %.2f seconds (%.1f fps)\n", total_time, (float)numframes/total_time);
243                 
244   return 0;
245 }