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 {
  55     FILE *stream = (FILE*)client_data;
  56     fprintf(stream, "[ERROR] %s", msg);
  57 }
  58 /**
  59 sample warning callback expecting a FILE* client object
  60 */
  61 static void warning_callback(const char *msg, void *client_data)
  62 {
  63     FILE *stream = (FILE*)client_data;
  64     fprintf(stream, "[WARNING] %s", msg);
  65 }
  66
  67 /* -------------------------------------------------------------------------- */
  68
  69
  70 int main(int argc, char *argv[])
  71 {
  72     mj2_dparameters_t mj2_parameters;           /* decompression parameters */
  73     opj_dinfo_t* dinfo;
  74     opj_event_mgr_t event_mgr;      /* event manager */
  75     opj_cio_t *cio = NULL;
  76     unsigned int tnum, snum;
  77     opj_mj2_t *movie;
  78     mj2_tk_t *track;
  79     mj2_sample_t *sample;
  80     unsigned char* frame_codestream;
  81     FILE *file, *outfile;
  82     char outfilename[50];
  83     opj_image_t *img = NULL;
  84     unsigned int max_codstrm_size = 0;
  85     double total_time = 0;
  86     unsigned int numframes = 0;
  87
  88     if (argc != 3) {
  89         printf("Usage: %s inputfile.mj2 outputfile.yuv\n", argv[0]);
  90         return 1;
  91     }
  92
  93     file = fopen(argv[1], "rb");
  94
  95     if (!file) {
  96         fprintf(stderr, "failed to open %s for reading\n", argv[1]);
  97         return 1;
  98     }
  99
 100     /* Checking output file */
 101     outfile = fopen(argv[2], "w");
 102     if (!file) {
 103         fprintf(stderr, "failed to open %s for writing\n", argv[2]);
 104         return 1;
 105     }
 106     fclose(outfile);
 107
 108     /*
 109     configure the event callbacks (not required)
 110     setting of each callback is optionnal
 111     */
 112     memset(&event_mgr, 0, sizeof(opj_event_mgr_t));
 113     event_mgr.error_handler = error_callback;
 114     event_mgr.warning_handler = warning_callback;
 115     event_mgr.info_handler = NULL;
 116
 117     /* get a MJ2 decompressor handle */
 118     dinfo = mj2_create_decompress();
 119     movie = (opj_mj2_t*)dinfo->mj2_handle;
 120
 121     /* catch events using our callbacks and give a local context */
 122     opj_set_event_mgr((opj_common_ptr)dinfo, &event_mgr, stderr);
 123
 124     memset(&mj2_parameters, 0, sizeof(mj2_dparameters_t));
 125     /* set J2K decoding parameters to default values */
 126     opj_set_default_decoder_parameters(&mj2_parameters.j2k_parameters);
 127
 128     /* setup the decoder decoding parameters using user parameters */
 129     mj2_setup_decoder(movie, &mj2_parameters);
 130
 131     if (mj2_read_struct(file, movie)) { /* Creating the movie structure */
 132         return 1;
 133     }
 134
 135     /* Decode first video track */
 136     for (tnum = 0;
 137             tnum < (unsigned int)(movie->num_htk + movie->num_stk + movie->num_vtk);
 138             tnum++) {
 139         if (movie->tk[tnum].track_type == 0) {
 140             break;
 141         }
 142     }
 143
 144     if (movie->tk[tnum].track_type != 0) {
 145         printf("Error. Movie does not contain any video track\n");
 146         return 1;
 147     }
 148
 149     track = &movie->tk[tnum];
 150
 151     /* Output info on first video tracl */
 152     fprintf(stdout,
 153             "The first video track contains %d frames.\nWidth: %d, Height: %d \n\n",
 154             track->num_samples, track->w, track->h);
 155
 156     max_codstrm_size = track->sample[0].sample_size - 8;
 157     frame_codestream = (unsigned char*) malloc(max_codstrm_size * sizeof(
 158                            unsigned char));
 159
 160     numframes = track->num_samples;
 161
 162     for (snum = 0; snum < numframes; snum++) {
 163         double init_time = opj_clock();
 164         double elapsed_time;
 165
 166         sample = &track->sample[snum];
 167         if (sample->sample_size - 8 > max_codstrm_size) {
 168             max_codstrm_size =  sample->sample_size - 8;
 169             if ((frame_codestream = (unsigned char*)
 170                                     realloc(frame_codestream, max_codstrm_size)) == NULL) {
 171                 printf("Error reallocation memory\n");
 172                 return 1;
 173             };
 174         }
 175         fseek(file, sample->offset + 8, SEEK_SET);
 176         fread(frame_codestream, sample->sample_size - 8, 1,
 177               file); /* Assuming that jp and ftyp markers size do */
 178
 179         /* open a byte stream */
 180         cio = opj_cio_open((opj_common_ptr)dinfo, frame_codestream,
 181                            sample->sample_size - 8);
 182
 183         img = opj_decode(dinfo, cio); /* Decode J2K to image */
 184
 185 #ifdef WANT_SYCC_TO_RGB
 186         if (img->color_space == CLRSPC_SYCC) {
 187             color_sycc_to_rgb(img);
 188         }
 189 #endif
 190
 191         if (img->icc_profile_buf) {
 192 #if defined(OPJ_HAVE_LIBLCMS1) || defined(OPJ_HAVE_LIBLCMS2)
 193             color_apply_icc_profile(img);
 194 #endif
 195
 196             free(img->icc_profile_buf);
 197             img->icc_profile_buf = NULL;
 198             img->icc_profile_len = 0;
 199         }
 200
 201         if (((img->numcomps == 3) && (img->comps[0].dx == img->comps[1].dx / 2)
 202                 && (img->comps[0].dx == img->comps[2].dx / 2) && (img->comps[0].dx == 1))
 203                 || (img->numcomps == 1)) {
 204
 205             if (!imagetoyuv(img, argv[2])) {  /* Convert image to YUV */
 206                 return 1;
 207             }
 208         } else if ((img->numcomps == 3) &&
 209                    (img->comps[0].dx == 1) && (img->comps[1].dx == 1) &&
 210                    (img->comps[2].dx == 1)) { /* If YUV 4:4:4 input --> to bmp */
 211             fprintf(stdout,
 212                     "The frames will be output in a bmp format (output_1.bmp, ...)\n");
 213             sprintf(outfilename, "output_%d.bmp", snum);
 214             if (imagetobmp(img, outfilename)) { /* Convert image to BMP */
 215                 return 1;
 216             }
 217
 218         } else {
 219             fprintf(stdout,
 220                     "Image component dimensions are unknown. Unable to output image\n");
 221             fprintf(stdout,
 222                     "The frames will be output in a j2k file (output_1.j2k, ...)\n");
 223
 224             sprintf(outfilename, "output_%d.j2k", snum);
 225             outfile = fopen(outfilename, "wb");
 226             if (!outfile) {
 227                 fprintf(stderr, "failed to open %s for writing\n", outfilename);
 228                 return 1;
 229             }
 230             fwrite(frame_codestream, sample->sample_size - 8, 1, outfile);
 231             fclose(outfile);
 232         }
 233         /* close the byte stream */
 234         opj_cio_close(cio);
 235         /* free image data structure */
 236         opj_image_destroy(img);
 237         elapsed_time = opj_clock() - init_time;
 238         fprintf(stderr, "Frame number %d/%d decoded in %.2f mseconds\n", snum + 1,
 239                 numframes, elapsed_time * 1000);
 240         total_time += elapsed_time;
 241
 242     }
 243
 244     free(frame_codestream);
 245     fclose(file);
 246
 247     /* free remaining structures */
 248     if (dinfo) {
 249         mj2_destroy_decompress((opj_mj2_t*)dinfo->mj2_handle);
 250     }
 251     free(dinfo);
 252
 253     fprintf(stdout, "%d frame(s) correctly decompressed\n", snum);
 254     fprintf(stdout, "Total decoding time: %.2f seconds (%.1f fps)\n", total_time,
 255             (float)numframes / total_time);
 256
 257     return 0;
 258 }