[trunk] WIP: correct some memory leaks in applications functions
[openjpeg.git] / libopenjpeg / j2k.c
index 4bd37eac4ed55fda4fa06f88ed73e52f76747861..e3e4c4987f606c5ae6207d1a5427da35be67f899 100644 (file)
@@ -699,6 +699,75 @@ static opj_bool j2k_read_unk_v2 (  opj_j2k_v2_t *p_j2k,
                                                                        OPJ_UINT32 *output_marker,
                                                                        struct opj_event_mgr * p_manager );
 
+/**
+ * Reads a MCT marker (Multiple Component Transform)
+ *
+ * @param      p_header_data   the data contained in the MCT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the MCT marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_mct ( opj_j2k_v2_t *p_j2k,
+                                                               OPJ_BYTE * p_header_data,
+                                                               OPJ_UINT32 p_header_size,
+                                                               struct opj_event_mgr * p_manager );
+
+/**
+ * Reads a MCC marker (Multiple Component Collection)
+ *
+ * @param      p_header_data   the data contained in the MCC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the MCC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_mcc ( opj_j2k_v2_t *p_j2k,
+                                                       OPJ_BYTE * p_header_data,
+                                                       OPJ_UINT32 p_header_size,
+                                                       struct opj_event_mgr * p_manager );
+
+/**
+ * Reads a MCO marker (Multiple Component Transform Ordering)
+ *
+ * @param      p_header_data   the data contained in the MCO box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the MCO marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_mco ( opj_j2k_v2_t *p_j2k,
+                                                               OPJ_BYTE * p_header_data,
+                                                               OPJ_UINT32 p_header_size,
+                                                               struct opj_event_mgr * p_manager );
+
+static opj_bool j2k_add_mct(opj_tcp_v2_t * p_tcp, opj_image_t * p_image, OPJ_UINT32 p_index);
+
+static void  j2k_read_int16_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_read_int32_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_read_float32_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_read_float64_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+
+static void  j2k_read_int16_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_read_int32_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_read_float32_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_read_float64_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+
+static void  j2k_write_float_to_int16 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_write_float_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_write_float_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+static void  j2k_write_float_to_float64 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+
+/**
+ * Reads a CBD marker (Component bit depth definition)
+ * @param      p_header_data   the data contained in the CBD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the CBD marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_cbd ( opj_j2k_v2_t *p_j2k,
+                                                       OPJ_BYTE * p_header_data,
+                                                       OPJ_UINT32 p_header_size,
+                                                       struct opj_event_mgr * p_manager);
+
+
 /**
  * Copy the image header from the jpeg2000 codec into an external image_header
  *
@@ -707,7 +776,7 @@ static opj_bool j2k_read_unk_v2 (   opj_j2k_v2_t *p_j2k,
  *
  * @return     true                    if the image header is correctly copy.
  */
-static opj_bool j2k_copy_img_header(opj_j2k_v2_t* p_j2k, opj_image_t* p_image);
+static opj_bool j2k_copy_img_header(opj_j2k_v2_t* p_j2k, opj_image_t** p_image);
 
 
 static void j2k_dump_MH_info(opj_j2k_v2_t* p_j2k, FILE* out_stream);
@@ -735,6 +804,43 @@ j2k_prog_order_t j2k_prog_order_list[] = {
        {(OPJ_PROG_ORDER)-1, ""}
 };
 
+/**
+ * FIXME DOC
+ */
+const OPJ_UINT32 MCT_ELEMENT_SIZE [] =
+{
+       2,
+       4,
+       4,
+       8
+};
+
+typedef void (* j2k_mct_function) (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem);
+
+const j2k_mct_function j2k_mct_read_functions_to_float [] =
+{
+       j2k_read_int16_to_float,
+       j2k_read_int32_to_float,
+       j2k_read_float32_to_float,
+       j2k_read_float64_to_float
+};
+
+const j2k_mct_function j2k_mct_read_functions_to_int32 [] =
+{
+       j2k_read_int16_to_int32,
+       j2k_read_int32_to_int32,
+       j2k_read_float32_to_int32,
+       j2k_read_float64_to_int32
+};
+
+const j2k_mct_function j2k_mct_write_functions_from_float [] =
+{
+       j2k_write_float_to_int16,
+       j2k_write_float_to_int32,
+       j2k_write_float_to_float,
+       j2k_write_float_to_float64
+};
+
 typedef struct opj_dec_memory_marker_handler
 {
        /** marker value */
@@ -791,12 +897,10 @@ const opj_dec_memory_marker_handler_t j2k_memory_marker_handler_tab [] =
   {J2K_MS_SOP, 0, 0},
   {J2K_MS_CRG, J2K_STATE_MH, j2k_read_crg_v2},
   {J2K_MS_COM, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_com_v2},
-#ifdef TODO_MS /* FIXME */
   {J2K_MS_MCT, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_mct},
   {J2K_MS_CBD, J2K_STATE_MH , j2k_read_cbd},
   {J2K_MS_MCC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_mcc},
   {J2K_MS_MCO, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_mco},
-#endif
 #ifdef USE_JPWL
 #ifdef TODO_MS /* FIXME */
   {J2K_MS_EPC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_epc},
@@ -814,11 +918,211 @@ const opj_dec_memory_marker_handler_t j2k_memory_marker_handler_tab [] =
 
 
 
+void  j2k_read_int16_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_FLOAT32 * l_dest_data = (OPJ_FLOAT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_UINT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_bytes(l_src_data,&l_temp,2);
+
+               l_src_data+=sizeof(OPJ_INT16);
+
+               *(l_dest_data++) = (OPJ_FLOAT32) l_temp;
+       }
+}
+
+void  j2k_read_int32_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_FLOAT32 * l_dest_data = (OPJ_FLOAT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_UINT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_bytes(l_src_data,&l_temp,4);
+
+               l_src_data+=sizeof(OPJ_INT32);
+
+               *(l_dest_data++) = (OPJ_FLOAT32) l_temp;
+       }
+}
+
+void  j2k_read_float32_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_FLOAT32 * l_dest_data = (OPJ_FLOAT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_FLOAT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_float(l_src_data,&l_temp);
+
+               l_src_data+=sizeof(OPJ_FLOAT32);
+
+               *(l_dest_data++) = l_temp;
+       }
+}
+
+void  j2k_read_float64_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_FLOAT32 * l_dest_data = (OPJ_FLOAT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_FLOAT64 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_double(l_src_data,&l_temp);
+
+               l_src_data+=sizeof(OPJ_FLOAT64);
+
+               *(l_dest_data++) = (OPJ_FLOAT32) l_temp;
+       }
+}
+
+void  j2k_read_int16_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_INT32 * l_dest_data = (OPJ_INT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_UINT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_bytes(l_src_data,&l_temp,2);
+
+               l_src_data+=sizeof(OPJ_INT16);
+
+               *(l_dest_data++) = (OPJ_INT32) l_temp;
+       }
+}
+
+void  j2k_read_int32_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_INT32 * l_dest_data = (OPJ_INT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_UINT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_bytes(l_src_data,&l_temp,4);
+
+               l_src_data+=sizeof(OPJ_INT32);
+
+               *(l_dest_data++) = (OPJ_INT32) l_temp;
+       }
+}
+
+void  j2k_read_float32_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_INT32 * l_dest_data = (OPJ_INT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_FLOAT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_float(l_src_data,&l_temp);
+
+               l_src_data+=sizeof(OPJ_FLOAT32);
+
+               *(l_dest_data++) = (OPJ_INT32) l_temp;
+       }
+}
+
+void  j2k_read_float64_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_src_data = (OPJ_BYTE *) p_src_data;
+       OPJ_INT32 * l_dest_data = (OPJ_INT32 *) p_dest_data;
+       OPJ_UINT32 i;
+       OPJ_FLOAT64 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               opj_read_double(l_src_data,&l_temp);
+
+               l_src_data+=sizeof(OPJ_FLOAT64);
+
+               *(l_dest_data++) = (OPJ_INT32) l_temp;
+       }
+}
+
+void  j2k_write_float_to_int16 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_dest_data = (OPJ_BYTE *) p_dest_data;
+       OPJ_FLOAT32 * l_src_data = (OPJ_FLOAT32 *) p_src_data;
+       OPJ_UINT32 i;
+       OPJ_UINT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               l_temp = (OPJ_UINT32) *(l_src_data++);
+
+               opj_write_bytes(l_dest_data,l_temp,sizeof(OPJ_INT16));
+
+               l_dest_data+=sizeof(OPJ_INT16);
+       }
+}
+
+void  j2k_write_float_to_int32 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_dest_data = (OPJ_BYTE *) p_dest_data;
+       OPJ_FLOAT32 * l_src_data = (OPJ_FLOAT32 *) p_src_data;
+       OPJ_UINT32 i;
+       OPJ_UINT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               l_temp = (OPJ_UINT32) *(l_src_data++);
+
+               opj_write_bytes(l_dest_data,l_temp,sizeof(OPJ_INT32));
+
+               l_dest_data+=sizeof(OPJ_INT32);
+       }
+}
+
+void  j2k_write_float_to_float (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_dest_data = (OPJ_BYTE *) p_dest_data;
+       OPJ_FLOAT32 * l_src_data = (OPJ_FLOAT32 *) p_src_data;
+       OPJ_UINT32 i;
+       OPJ_FLOAT32 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               l_temp = (OPJ_FLOAT32) *(l_src_data++);
+
+               opj_write_float(l_dest_data,l_temp);
+
+               l_dest_data+=sizeof(OPJ_FLOAT32);
+       }
+}
+
+void  j2k_write_float_to_float64 (const void * p_src_data, void * p_dest_data, OPJ_UINT32 p_nb_elem)
+{
+       OPJ_BYTE * l_dest_data = (OPJ_BYTE *) p_dest_data;
+       OPJ_FLOAT32 * l_src_data = (OPJ_FLOAT32 *) p_src_data;
+       OPJ_UINT32 i;
+       OPJ_FLOAT64 l_temp;
+
+       for (i=0;i<p_nb_elem;++i) {
+               l_temp = (OPJ_FLOAT64) *(l_src_data++);
+
+               opj_write_double(l_dest_data,l_temp);
+
+               l_dest_data+=sizeof(OPJ_FLOAT64);
+       }
+}
+
+
+/**
+ * Converts an enum type progression order to string type.
+ *
+ * @param prg_order            the progression order to get.
+ *
+ * @return     the string representation of the given progression order.
+ */
 char *j2k_convert_progression_order(OPJ_PROG_ORDER prg_order){
        j2k_prog_order_t *po;
        for(po = j2k_prog_order_list; po->enum_prog != -1; po++ ){
                if(po->enum_prog == prg_order){
-                       break;
+                       return po->str_prog;
                }
        }
        return po->str_prog;
@@ -3352,10 +3656,9 @@ opj_bool j2k_read_sot_v2 (
 
        /* Ref A.4.2: Psot could be equal zero if it is the last tile-part of the codestream.*/
        if (!l_tot_len) {
-               opj_event_msg_v2(p_manager, EVT_ERROR, "Psot value of the current tile-part is equal to zero, "
-                               "for the moment we couldn't manage this case (need to compute the number of byte left"
-                               " in the codestream).\n");
-               return OPJ_FALSE;
+               opj_event_msg_v2(p_manager, EVT_INFO, "Psot value of the current tile-part is equal to zero, "
+                               "we assuming it is the last tile-part of the codestream.\n");
+               p_j2k->m_specific_param.m_decoder.m_last_tile_part = 1;
        }
 
        opj_read_bytes(p_header_data,&l_current_part ,1);       /* TPsot */
@@ -3575,7 +3878,11 @@ opj_bool j2k_read_sod_v2 (
        assert(p_stream != 00);
 
        l_tcp = &(p_j2k->m_cp.tcps[p_j2k->m_current_tile_number]);
-       p_j2k->m_specific_param.m_decoder.m_sot_length -= 2;
+
+       if (p_j2k->m_specific_param.m_decoder.m_last_tile_part)
+               p_j2k->m_specific_param.m_decoder.m_sot_length = opj_stream_get_number_byte_left(p_stream) - 2;
+       else
+               p_j2k->m_specific_param.m_decoder.m_sot_length -= 2;
 
        l_current_data = &(l_tcp->m_data);
        l_tile_len = &l_tcp->m_data_size;
@@ -4029,6 +4336,538 @@ static opj_dec_mstabent_t *j2k_dec_mstab_lookup(int id) {
        return e;
 }
 
+
+/**
+ * Reads a MCT marker (Multiple Component Transform)
+ *
+ * @param      p_header_data   the data contained in the MCT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the MCT marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_mct (        opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager )
+{
+       OPJ_UINT32 i;
+       opj_tcp_v2_t *l_tcp = 00;
+       OPJ_UINT32 l_tmp;
+       OPJ_UINT32 l_indix;
+       opj_mct_data_t * l_mct_data;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+
+       l_tcp = p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH ?
+                       &p_j2k->m_cp.tcps[p_j2k->m_current_tile_number] :
+                       p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       if (p_header_size < 2) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCT marker\n");
+               return OPJ_FALSE;
+       }
+
+       /* first marker */
+       opj_read_bytes(p_header_data,&l_tmp,2);                         /* Zmct */
+       p_header_data += 2;
+       if (l_tmp != 0) {
+               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge mct data within multiple MCT records\n");
+               return OPJ_TRUE;
+       }
+
+       if(p_header_size <= 6) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCT marker\n");
+               return OPJ_FALSE;
+       }
+
+       /* Imct -> no need for other values, take the first, type is double with decorrelation x0000 1101 0000 0000*/
+       opj_read_bytes(p_header_data,&l_tmp,2);                         /* Imct */
+       p_header_data += 2;
+
+       l_indix = l_tmp & 0xff;
+       l_mct_data = l_tcp->m_mct_records;
+
+       for (i=0;i<l_tcp->m_nb_mct_records;++i) {
+               if (l_mct_data->m_index == l_indix) {
+                       break;
+               }
+               ++l_mct_data;
+       }
+
+       /* NOT FOUND */
+       if (i == l_tcp->m_nb_mct_records) {
+               if (l_tcp->m_nb_mct_records == l_tcp->m_nb_max_mct_records) {
+                       l_tcp->m_nb_max_mct_records += J2K_MCT_DEFAULT_NB_RECORDS;
+
+                       l_tcp->m_mct_records = (opj_mct_data_t*)opj_realloc(l_tcp->m_mct_records,l_tcp->m_nb_max_mct_records * sizeof(opj_mct_data_t));
+                       if(! l_tcp->m_mct_records) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCT marker\n");
+                               return OPJ_FALSE;
+                       }
+
+                       l_mct_data = l_tcp->m_mct_records + l_tcp->m_nb_mct_records;
+                       memset(l_mct_data ,0,(l_tcp->m_nb_max_mct_records - l_tcp->m_nb_mct_records) * sizeof(opj_mct_data_t));
+               }
+
+               l_mct_data = l_tcp->m_mct_records + l_tcp->m_nb_mct_records;
+       }
+
+       if (l_mct_data->m_data) {
+               opj_free(l_mct_data->m_data);
+               l_mct_data->m_data = 00;
+       }
+
+       l_mct_data->m_index = l_indix;
+       l_mct_data->m_array_type = (J2K_MCT_ARRAY_TYPE)((l_tmp  >> 8) & 3);
+       l_mct_data->m_element_type = (J2K_MCT_ELEMENT_TYPE)((l_tmp  >> 10) & 3);
+
+       opj_read_bytes(p_header_data,&l_tmp,2);                         /* Ymct */
+       p_header_data+=2;
+       if (l_tmp != 0) {
+               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge multiple MCT markers\n");
+               return OPJ_TRUE;
+       }
+
+       p_header_size -= 6;
+
+       l_mct_data->m_data = (OPJ_BYTE*)opj_malloc(p_header_size);
+       if (! l_mct_data->m_data) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCT marker\n");
+               return OPJ_FALSE;
+       }
+       memcpy(l_mct_data->m_data,p_header_data,p_header_size);
+
+       l_mct_data->m_data_size = p_header_size;
+       ++l_tcp->m_nb_mct_records;
+
+       return OPJ_TRUE;
+}
+
+/**
+ * Reads a MCC marker (Multiple Component Collection)
+ *
+ * @param      p_header_data   the data contained in the MCC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the MCC marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_mcc (        opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager )
+{
+       OPJ_UINT32 i,j;
+       OPJ_UINT32 l_tmp;
+       OPJ_UINT32 l_indix;
+       opj_tcp_v2_t * l_tcp;
+       opj_simple_mcc_decorrelation_data_t * l_mcc_record;
+       opj_mct_data_t * l_mct_data;
+       OPJ_UINT32 l_nb_collections;
+       OPJ_UINT32 l_nb_comps;
+       OPJ_UINT32 l_nb_bytes_by_comp;
+
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_tcp = p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH ?
+                       &p_j2k->m_cp.tcps[p_j2k->m_current_tile_number] :
+                       p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       if (p_header_size < 2) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+               return OPJ_FALSE;
+       }
+
+       /* first marker */
+       opj_read_bytes(p_header_data,&l_tmp,2);                         /* Zmcc */
+       p_header_data += 2;
+       if (l_tmp != 0) {
+               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge multiple data spanning\n");
+               return OPJ_TRUE;
+       }
+
+       if (p_header_size < 7) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(p_header_data,&l_indix,1); /* Imcc -> no need for other values, take the first */
+       ++p_header_data;
+
+       l_mcc_record = l_tcp->m_mcc_records;
+
+       for(i=0;i<l_tcp->m_nb_mcc_records;++i) {
+               if (l_mcc_record->m_index == l_indix) {
+                       break;
+               }
+               ++l_mcc_record;
+       }
+
+       /** NOT FOUND */
+       if (i == l_tcp->m_nb_mcc_records) {
+               if (l_tcp->m_nb_mcc_records == l_tcp->m_nb_max_mcc_records) {
+                       l_tcp->m_nb_max_mcc_records += J2K_MCC_DEFAULT_NB_RECORDS;
+
+                       l_tcp->m_mcc_records = (opj_simple_mcc_decorrelation_data_t*)
+                                       opj_realloc(l_tcp->m_mcc_records,l_tcp->m_nb_max_mcc_records * sizeof(opj_simple_mcc_decorrelation_data_t));
+                       if (! l_tcp->m_mcc_records) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+                               return OPJ_FALSE;
+                       }
+                       l_mcc_record = l_tcp->m_mcc_records + l_tcp->m_nb_mcc_records;
+                       memset(l_mcc_record,0,(l_tcp->m_nb_max_mcc_records-l_tcp->m_nb_mcc_records) * sizeof(opj_simple_mcc_decorrelation_data_t));
+               }
+               l_mcc_record = l_tcp->m_mcc_records + l_tcp->m_nb_mcc_records;
+       }
+       l_mcc_record->m_index = l_indix;
+
+       /* only one marker atm */
+       opj_read_bytes(p_header_data,&l_tmp,2);                         /* Ymcc */
+       p_header_data+=2;
+       if (l_tmp != 0) {
+               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge multiple data spanning\n");
+               return OPJ_TRUE;
+       }
+
+       opj_read_bytes(p_header_data,&l_nb_collections,2);                              /* Qmcc -> number of collections -> 1 */
+       p_header_data+=2;
+
+       if (l_nb_collections > 1) {
+               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge multiple collections\n");
+               return OPJ_TRUE;
+       }
+
+       p_header_size -= 7;
+
+       for (i=0;i<l_nb_collections;++i) {
+               if (p_header_size < 3) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+                       return OPJ_FALSE;
+               }
+
+               opj_read_bytes(p_header_data,&l_tmp,1); /* Xmcci type of component transformation -> array based decorrelation */
+               ++p_header_data;
+
+               if (l_tmp != 1) {
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge collections other than array decorrelation\n");
+                       return OPJ_TRUE;
+               }
+
+               opj_read_bytes(p_header_data,&l_nb_comps,2);
+
+               p_header_data+=2;
+               p_header_size-=3;
+
+               l_nb_bytes_by_comp = 1 + (l_nb_comps>>15);
+               l_mcc_record->m_nb_comps = l_nb_comps & 0x7fff;
+
+               if (p_header_size < (l_nb_bytes_by_comp * l_mcc_record->m_nb_comps + 2)) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+                       return OPJ_FALSE;
+               }
+
+               p_header_size -= (l_nb_bytes_by_comp * l_mcc_record->m_nb_comps + 2);
+
+               for (j=0;j<l_mcc_record->m_nb_comps;++j) {
+                       opj_read_bytes(p_header_data,&l_tmp,l_nb_bytes_by_comp);        /* Cmccij Component offset*/
+                       p_header_data+=l_nb_bytes_by_comp;
+
+                       if (l_tmp != j) {
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge collections with indix shuffle\n");
+                               return OPJ_TRUE;
+                       }
+               }
+
+               opj_read_bytes(p_header_data,&l_nb_comps,2);
+               p_header_data+=2;
+
+               l_nb_bytes_by_comp = 1 + (l_nb_comps>>15);
+               l_nb_comps &= 0x7fff;
+
+               if (l_nb_comps != l_mcc_record->m_nb_comps) {
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge collections without same number of indixes\n");
+                       return OPJ_TRUE;
+               }
+
+               if (p_header_size < (l_nb_bytes_by_comp * l_mcc_record->m_nb_comps + 3)) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+                       return OPJ_FALSE;
+               }
+
+               p_header_size -= (l_nb_bytes_by_comp * l_mcc_record->m_nb_comps + 3);
+
+               for (j=0;j<l_mcc_record->m_nb_comps;++j) {
+                       opj_read_bytes(p_header_data,&l_tmp,l_nb_bytes_by_comp);        /* Wmccij Component offset*/
+                       p_header_data+=l_nb_bytes_by_comp;
+
+                       if (l_tmp != j) {
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge collections with indix shuffle\n");
+                               return OPJ_TRUE;
+                       }
+               }
+
+               opj_read_bytes(p_header_data,&l_tmp,3); /* Wmccij Component offset*/
+               p_header_data += 3;
+
+               l_mcc_record->m_is_irreversible = ! ((l_tmp>>16) & 1);
+               l_mcc_record->m_decorrelation_array = 00;
+               l_mcc_record->m_offset_array = 00;
+
+               l_indix = l_tmp & 0xff;
+               if (l_indix != 0) {
+                       l_mct_data = l_tcp->m_mct_records;
+                       for (j=0;j<l_tcp->m_nb_mct_records;++j) {
+                               if (l_mct_data->m_index == l_indix) {
+                                       l_mcc_record->m_decorrelation_array = l_mct_data;
+                                       break;
+                               }
+                               ++l_mct_data;
+                       }
+
+                       if (l_mcc_record->m_decorrelation_array == 00) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+                               return OPJ_FALSE;
+                       }
+               }
+
+               l_indix = (l_tmp >> 8) & 0xff;
+               if (l_indix != 0) {
+                       l_mct_data = l_tcp->m_mct_records;
+                       for (j=0;j<l_tcp->m_nb_mct_records;++j) {
+                               if (l_mct_data->m_index == l_indix) {
+                                       l_mcc_record->m_offset_array = l_mct_data;
+                                       break;
+                               }
+                               ++l_mct_data;
+                       }
+
+                       if (l_mcc_record->m_offset_array == 00) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+                               return OPJ_FALSE;
+                       }
+               }
+       }
+
+       if (p_header_size != 0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCC marker\n");
+               return OPJ_FALSE;
+       }
+
+       ++l_tcp->m_nb_mcc_records;
+
+       return OPJ_TRUE;
+}
+
+/**
+ * Reads a MCO marker (Multiple Component Transform Ordering)
+ *
+ * @param      p_header_data   the data contained in the MCO box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the MCO marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_mco (        opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager )
+{
+       OPJ_UINT32 l_tmp, i;
+       OPJ_UINT32 l_nb_stages;
+       opj_tcp_v2_t * l_tcp;
+       opj_tccp_t * l_tccp;
+       opj_image_t * l_image;
+       opj_image_comp_t * l_img_comp;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_image = p_j2k->m_image;
+       l_tcp = p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH ?
+                       &p_j2k->m_cp.tcps[p_j2k->m_current_tile_number] :
+                       p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       if (p_header_size < 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading MCO marker\n");
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(p_header_data,&l_nb_stages,1);                           /* Nmco : only one tranform stage*/
+       ++p_header_data;
+
+       if (l_nb_stages > 1) {
+               opj_event_msg_v2(p_manager, EVT_WARNING, "Cannot take in charge multiple transformation stages.\n");
+               return OPJ_TRUE;
+       }
+
+       if (p_header_size != l_nb_stages + 1) {
+               opj_event_msg_v2(p_manager, EVT_WARNING, "Error reading MCO marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_tccp = l_tcp->tccps;
+       l_img_comp = l_image->comps;
+
+       for (i=0;i<l_image->numcomps;++i) {
+               l_tccp->m_dc_level_shift = 0;
+               ++l_tccp;
+       }
+
+       if (l_tcp->m_mct_decoding_matrix) {
+               opj_free(l_tcp->m_mct_decoding_matrix);
+               l_tcp->m_mct_decoding_matrix = 00;
+       }
+
+       for (i=0;i<l_nb_stages;++i) {
+               opj_read_bytes(p_header_data,&l_tmp,1);
+               ++p_header_data;
+
+               if (! j2k_add_mct(l_tcp,p_j2k->m_image,l_tmp)) {
+                       return OPJ_FALSE;
+               }
+       }
+
+       return OPJ_TRUE;
+}
+
+opj_bool j2k_add_mct(opj_tcp_v2_t * p_tcp, opj_image_t * p_image, OPJ_UINT32 p_index)
+{
+       OPJ_UINT32 i;
+       opj_simple_mcc_decorrelation_data_t * l_mcc_record;
+       opj_mct_data_t * l_deco_array, * l_offset_array;
+       OPJ_UINT32 l_data_size,l_mct_size, l_offset_size;
+       OPJ_UINT32 l_nb_elem;
+       OPJ_UINT32 * l_offset_data, * l_current_offset_data;
+       opj_tccp_t * l_tccp;
+
+       /* preconditions */
+       assert(p_tcp != 00);
+
+       l_mcc_record = p_tcp->m_mcc_records;
+
+       for (i=0;i<p_tcp->m_nb_mcc_records;++i) {
+               if (l_mcc_record->m_index == p_index) {
+                       break;
+               }
+       }
+
+       if (i==p_tcp->m_nb_mcc_records) {
+               /** element discarded **/
+               return OPJ_TRUE;
+       }
+
+       if (l_mcc_record->m_nb_comps != p_image->numcomps) {
+               /** do not support number of comps != image */
+               return OPJ_TRUE;
+       }
+
+       l_deco_array = l_mcc_record->m_decorrelation_array;
+
+       if (l_deco_array) {
+               l_data_size = MCT_ELEMENT_SIZE[l_deco_array->m_element_type] * p_image->numcomps * p_image->numcomps;
+               if (l_deco_array->m_data_size != l_data_size) {
+                       return OPJ_FALSE;
+               }
+
+               l_nb_elem = p_image->numcomps * p_image->numcomps;
+               l_mct_size = l_nb_elem * sizeof(OPJ_FLOAT32);
+               p_tcp->m_mct_decoding_matrix = (OPJ_FLOAT32*)opj_malloc(l_mct_size);
+
+               if (! p_tcp->m_mct_decoding_matrix ) {
+                       return OPJ_FALSE;
+               }
+
+               j2k_mct_read_functions_to_float[l_deco_array->m_element_type](l_deco_array->m_data,p_tcp->m_mct_decoding_matrix,l_nb_elem);
+       }
+
+       l_offset_array = l_mcc_record->m_offset_array;
+
+       if (l_offset_array) {
+               l_data_size = MCT_ELEMENT_SIZE[l_offset_array->m_element_type] * p_image->numcomps;
+               if (l_offset_array->m_data_size != l_data_size) {
+                       return OPJ_FALSE;
+               }
+
+               l_nb_elem = p_image->numcomps;
+               l_offset_size = l_nb_elem * sizeof(OPJ_UINT32);
+               l_offset_data = (OPJ_UINT32*)opj_malloc(l_offset_size);
+
+               if (! l_offset_data ) {
+                       return OPJ_FALSE;
+               }
+
+               j2k_mct_read_functions_to_int32[l_offset_array->m_element_type](l_offset_array->m_data,l_offset_data,l_nb_elem);
+
+               l_tccp = p_tcp->tccps;
+               l_current_offset_data = l_offset_data;
+
+               for (i=0;i<p_image->numcomps;++i) {
+                       l_tccp->m_dc_level_shift = *(l_current_offset_data++);
+                       ++l_tccp;
+               }
+
+               opj_free(l_offset_data);
+       }
+
+       return OPJ_TRUE;
+}
+
+/**
+ * Reads a CBD marker (Component bit depth definition)
+ * @param      p_header_data   the data contained in the CBD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the CBD marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_cbd (        opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager)
+{
+       OPJ_UINT32 l_nb_comp,l_num_comp;
+       OPJ_UINT32 l_comp_def;
+       OPJ_UINT32 i;
+       opj_image_comp_t * l_comp = 00;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_num_comp = p_j2k->m_image->numcomps;
+
+       if (p_header_size != (p_j2k->m_image->numcomps + 2)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Crror reading CBD marker\n");
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(p_header_data,&l_nb_comp,2);                             /* Ncbd */
+       p_header_data+=2;
+
+       if (l_nb_comp != l_num_comp) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Crror reading CBD marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_comp = p_j2k->m_image->comps;
+       for (i=0;i<l_num_comp;++i) {
+               opj_read_bytes(p_header_data,&l_comp_def,1);                    /* Component bit depth */
+               ++p_header_data;
+        l_comp->sgnd = (l_comp_def>>7) & 1;
+               l_comp->prec = (l_comp_def&0x7f) + 1;
+               ++l_comp;
+       }
+
+       return OPJ_TRUE;
+}
+
+
 /* ----------------------------------------------------------------------- */
 /* J2K / JPT decoder interface                                             */
 /* ----------------------------------------------------------------------- */
@@ -4957,13 +5796,12 @@ opj_bool j2k_end_decompress(
  */
 opj_bool j2k_read_header(      struct opj_stream_private *p_stream,
                                                        opj_j2k_v2_t* p_j2k,
-                                                       opj_image_t* p_image,
+                                                       opj_image_t** p_image,
                                                        struct opj_event_mgr* p_manager )
 {
        /* preconditions */
        assert(p_j2k != 00);
        assert(p_stream != 00);
-       assert(p_image != 00);
        assert(p_manager != 00);
 
        /* create an empty image header */
@@ -4992,6 +5830,11 @@ opj_bool j2k_read_header(        struct opj_stream_private *p_stream,
                return OPJ_FALSE;
        }
 
+       *p_image = opj_image_create0();
+       if (! (*p_image)) {
+               return OPJ_FALSE;
+       }
+
        if (! j2k_copy_img_header(p_j2k, p_image)){
                opj_image_destroy(p_j2k->m_image);
                p_j2k->m_image = NULL;
@@ -5002,42 +5845,41 @@ opj_bool j2k_read_header(       struct opj_stream_private *p_stream,
        return OPJ_TRUE;
 }
 
-opj_bool j2k_copy_img_header(opj_j2k_v2_t* p_j2k, opj_image_t* p_image)
+opj_bool j2k_copy_img_header(opj_j2k_v2_t* p_j2k, opj_image_t** p_image)
 {
        OPJ_UINT16 compno;
 
        /* preconditions */
        assert(p_j2k != 00);
-       assert(p_image != 00);
 
-       p_image->x0 = p_j2k->m_image->x0;
-       p_image->y0 = p_j2k->m_image->y0;
-       p_image->x1 = p_j2k->m_image->x1;
-       p_image->y1 = p_j2k->m_image->y1;
+       (*p_image)->x0 = p_j2k->m_image->x0;
+       (*p_image)->y0 = p_j2k->m_image->y0;
+       (*p_image)->x1 = p_j2k->m_image->x1;
+       (*p_image)->y1 = p_j2k->m_image->y1;
 
-       p_image->numcomps = p_j2k->m_image->numcomps;
-       p_image->comps = (opj_image_comp_t*)opj_malloc(p_image->numcomps * sizeof(opj_image_comp_t));
-       if (!p_image->comps)
+       (*p_image)->numcomps = p_j2k->m_image->numcomps;
+       (*p_image)->comps = (opj_image_comp_t*)opj_malloc((*p_image)->numcomps * sizeof(opj_image_comp_t));
+       if (!(*p_image)->comps)
                return OPJ_FALSE;
-       for (compno=0; compno < p_image->numcomps; compno++){
-               memcpy( &(p_image->comps[compno]),
+       for (compno=0; compno < (*p_image)->numcomps; compno++){
+               memcpy( &((*p_image)->comps[compno]),
                                &(p_j2k->m_image->comps[compno]),
                                sizeof(opj_image_comp_t));
        }
 
-       p_image->color_space = p_j2k->m_image->color_space;
-       p_image->icc_profile_len = p_j2k->m_image->icc_profile_len;
-       if (!p_image->icc_profile_len) {
+       (*p_image)->color_space = p_j2k->m_image->color_space;
+       (*p_image)->icc_profile_len = p_j2k->m_image->icc_profile_len;
+       if (!(*p_image)->icc_profile_len) {
 
-               p_image->icc_profile_buf = (unsigned char*)opj_malloc(p_image->icc_profile_len);
-               if (!p_image->icc_profile_buf)
+               (*p_image)->icc_profile_buf = (unsigned char*)opj_malloc((*p_image)->icc_profile_len);
+               if (!(*p_image)->icc_profile_buf)
                        return OPJ_FALSE;
-               memcpy( p_image->icc_profile_buf,
+               memcpy( (*p_image)->icc_profile_buf,
                                p_j2k->m_image->icc_profile_buf,
                                p_j2k->m_image->icc_profile_len);
        }
        else
-               p_image->icc_profile_buf = NULL;
+               (*p_image)->icc_profile_buf = NULL;
 
        return OPJ_TRUE;
 }
@@ -5488,7 +6330,7 @@ void j2k_destroy (opj_j2k_v2_t *p_j2k)
        tcd_destroy_v2(p_j2k->m_tcd);
 
        j2k_cp_destroy(&(p_j2k->m_cp));
-       memset(&(p_j2k->m_cp),0,sizeof(opj_cp_t));
+       memset(&(p_j2k->m_cp),0,sizeof(opj_cp_v2_t));
 
        opj_procedure_list_destroy(p_j2k->m_procedure_list);
        p_j2k->m_procedure_list = 00;
@@ -5499,22 +6341,26 @@ void j2k_destroy (opj_j2k_v2_t *p_j2k)
        j2k_destroy_cstr_index(p_j2k->cstr_index);
        p_j2k->cstr_index = NULL;
 
+       opj_image_destroy(p_j2k->m_image);
+       p_j2k->m_image = NULL;
+
        opj_free(p_j2k);
 }
 
 void j2k_destroy_cstr_index (opj_codestream_index_t *p_cstr_ind)
 {
-       if (!p_cstr_ind) {
-               return;
-       }
+       if (p_cstr_ind) {
 
-       if (p_cstr_ind->marker) {
-               opj_free(p_cstr_ind->marker);
-               p_cstr_ind->marker = NULL;
-       }
+               if (p_cstr_ind->marker) {
+                       opj_free(p_cstr_ind->marker);
+                       p_cstr_ind->marker = NULL;
+               }
+
+               if (p_cstr_ind->tile_index) {
+                       // FIXME not used for the moment
+               }
 
-       if (p_cstr_ind->tile_index) {
-               // FIXME
+               opj_free(p_cstr_ind);
        }
 }
 
@@ -6758,6 +7604,8 @@ opj_codestream_info_v2_t* j2k_get_cstr_info(opj_j2k_v2_t* p_j2k)
        cstr_info->tw = p_j2k->m_cp.tw;
        cstr_info->th = p_j2k->m_cp.th;
 
+       cstr_info->tile_info = NULL; /* Not fill from the main header*/
+
        l_default_tile = p_j2k->m_specific_param.m_decoder.m_default_tcp;
 
        cstr_info->m_default_tile_info.csty = l_default_tile->csty;
@@ -6765,7 +7613,7 @@ opj_codestream_info_v2_t* j2k_get_cstr_info(opj_j2k_v2_t* p_j2k)
        cstr_info->m_default_tile_info.numlayers = l_default_tile->numlayers;
        cstr_info->m_default_tile_info.mct = l_default_tile->mct;
 
-       cstr_info->m_default_tile_info.tccp_info = (opj_tccp_info_t*) opj_calloc(1,sizeof(opj_tccp_info_t));
+       cstr_info->m_default_tile_info.tccp_info = (opj_tccp_info_t*) opj_calloc(cstr_info->nbcomps, sizeof(opj_tccp_info_t));
 
        for (compno = 0; compno < numcomps; compno++) {
                opj_tccp_t *l_tccp = &(l_default_tile->tccps[compno]);
@@ -6831,6 +7679,7 @@ opj_codestream_index_t* j2k_get_cstr_index(opj_j2k_v2_t* p_j2k)
 
        memcpy(l_cstr_index->marker, p_j2k->cstr_index->marker, l_cstr_index->marknum * sizeof(opj_marker_info_t) );
 
+       l_cstr_index->tile_index = NULL; /* FIXME not used for the moment*/
 
        return l_cstr_index;
 }
@@ -6926,20 +7775,24 @@ opj_bool j2k_decode_v2( opj_j2k_v2_t * p_j2k,
 {
        OPJ_UINT32 compno;
 
+       if (!p_image)
+               return OPJ_FALSE;
+
        /* customization of the decoding */
        j2k_setup_decoding(p_j2k);
 
-       /* write header */
+       /* Decode the codestream */
        if (! j2k_exec (p_j2k,p_j2k->m_procedure_list,p_stream,p_manager)) {
                opj_image_destroy(p_j2k->m_image);
                p_j2k->m_image = NULL;
                return OPJ_FALSE;
        }
 
+       /* Copy data from codec to output image*/
        for (compno = 0; compno < p_image->numcomps; compno++) {
                p_image->comps[compno].data = p_j2k->m_image->comps[compno].data;
                p_j2k->m_image->comps[compno].data = NULL;
        }
 
-       return OPJ_TRUE /*p_j2k->m_image*/;
+       return OPJ_TRUE;
 }