[trunk] WIP: add output elements about decoding of jp2 files with last tile part...
[openjpeg.git] / libopenjpeg / j2k.c
index 556833c63e611172de12b45074f5ca4716c7a5b9..4bd37eac4ed55fda4fa06f88ed73e52f76747861 100644 (file)
@@ -6,6 +6,7 @@
  * Copyright (c) 2003-2007, Francois-Olivier Devaux and Antonin Descampe
  * Copyright (c) 2005, Herve Drolon, FreeImage Team
  * Copyright (c) 2006-2007, Parvatha Elangovan
+ * Copyright (c) 2010-2011, Kaori Hagihara
  * All rights reserved.
  *
  * Redistribution and use in source and binary forms, with or without
 /** @name Local static functions */
 /*@{*/
 
+/**
+ * Sets up the procedures to do on reading header. Developpers wanting to extend the library can add their own reading procedures.
+ */
+void j2k_setup_header_reading (opj_j2k_v2_t *p_j2k);
+
+/**
+ * The read header procedure.
+ */
+opj_bool j2k_read_header_procedure(
+                                                           opj_j2k_v2_t *p_j2k,
+                                                               struct opj_stream_private *p_stream,
+                                                               struct opj_event_mgr * p_manager);
+
+/**
+ * The default decoding validation procedure without any extension.
+ *
+ * @param      p_j2k                   the jpeg2000 codec to validate.
+ * @param      p_stream                                the input stream to validate.
+ * @param      p_manager               the user event manager.
+ *
+ * @return true if the parameters are correct.
+ */
+opj_bool j2k_decoding_validation (
+                                                               opj_j2k_v2_t * p_j2k,
+                                                               opj_stream_private_t *p_stream,
+                                                               opj_event_mgr_t * p_manager
+                                                       );
+
+/**
+ * Sets up the validation ,i.e. adds the procedures to lauch to make sure the codec parameters
+ * are valid. Developpers wanting to extend the library can add their own validation procedures.
+ */
+static void j2k_setup_decoding_validation (opj_j2k_v2_t *p_j2k);
+
+/**
+ * Builds the tcd decoder to use to decode tile.
+ */
+opj_bool j2k_build_decoder (
+                                               opj_j2k_v2_t * p_j2k,
+                                               opj_stream_private_t *p_stream,
+                                               opj_event_mgr_t * p_manager
+                                               );
+
+/**
+ * Excutes the given procedures on the given codec.
+ *
+ * @param      p_procedure_list        the list of procedures to execute
+ * @param      p_j2k                                   the jpeg2000 codec to execute the procedures on.
+ * @param      p_stream                                        the stream to execute the procedures on.
+ * @param      p_manager                       the user manager.
+ *
+ * @return     true                            if all the procedures were successfully executed.
+ */
+static opj_bool j2k_exec (
+                                       opj_j2k_v2_t * p_j2k,
+                                       opj_procedure_list_t * p_procedure_list,
+                                       opj_stream_private_t *p_stream,
+                                       opj_event_mgr_t * p_manager
+                                 );
+
+/**
+ * Copies the decoding tile parameters onto all the tile parameters.
+ * Creates also the tile decoder.
+ */
+opj_bool j2k_copy_default_tcp_and_create_tcd ( opj_j2k_v2_t * p_j2k,
+                                                                                               opj_stream_private_t *p_stream,
+                                                                                               opj_event_mgr_t * p_manager );
+
+/**
+ * Reads the lookup table containing all the marker, status and action, and returns the handler associated
+ * with the marker value.
+ * @param      p_id            Marker value to look up
+ *
+ * @return     the handler associated with the id.
+*/
+static const struct opj_dec_memory_marker_handler * j2k_get_marker_handler (OPJ_UINT32 p_id);
+
+/**
+ * Destroys a tile coding parameter structure.
+ *
+ * @param      p_tcp           the tile coding parameter to destroy.
+ */
+static void j2k_tcp_destroy (opj_tcp_v2_t *p_tcp);
+
+/**
+ * Destroys a codestream index structure.
+ *
+ * @param      p_cstr_ind      the codestream index parameter to destroy.
+ */
+static void j2k_destroy_cstr_index (opj_codestream_index_t* p_cstr_ind);
+
+/**
+ * Destroys a coding parameter structure.
+ *
+ * @param      p_cp            the coding parameter to destroy.
+ */
+static void j2k_cp_destroy (opj_cp_v2_t *p_cp);
+
+
+/**
+ * Reads a SPCod or SPCoc element, i.e. the coding style of a given component of a tile.
+ * @param      p_header_data   the data contained in the COM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the COM marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_SPCod_SPCoc(
+                                                       opj_j2k_v2_t *p_j2k,
+                                                       OPJ_UINT32 compno,
+                                                       OPJ_BYTE * p_header_data,
+                                                       OPJ_UINT32 * p_header_size,
+                                                       struct opj_event_mgr * p_manager
+                                                       );
+
+/**
+ * Reads a SQcd or SQcc element, i.e. the quantization values of a band in the QCD or QCC.
+ *
+ * @param      p_tile_no               the tile to output.
+ * @param      p_comp_no               the component number to output.
+ * @param      p_data                  the data buffer.
+ * @param      p_header_size   pointer to the size of the data buffer, it is changed by the function.
+ * @param      p_j2k                           J2K codec.
+ * @param      p_manager               the user event manager.
+ *
+*/
+static opj_bool j2k_read_SQcd_SQcc(
+                                                       opj_j2k_v2_t *p_j2k,
+                                                       OPJ_UINT32 compno,
+                                                       OPJ_BYTE * p_header_data,
+                                                       OPJ_UINT32 * p_header_size,
+                                                       struct opj_event_mgr * p_manager
+                                       );
+
+/**
+ * Copies the tile component parameters of all the component from the first tile component.
+ *
+ * @param              p_j2k           the J2k codec.
+ */
+static void j2k_copy_tile_component_parameters(
+                                                       opj_j2k_v2_t *p_j2k
+                                                       );
+
+/**
+ * Copies the tile quantization parameters of all the component from the first tile component.
+ *
+ * @param              p_j2k           the J2k codec.
+ */
+static void j2k_copy_tile_quantization_parameters(
+                                                       opj_j2k_v2_t *p_j2k
+                                                       );
+
+/**
+ * Reads the tiles.
+ */
+opj_bool j2k_decode_tiles (            opj_j2k_v2_t *p_j2k,
+                                                               opj_stream_private_t *p_stream,
+                                                               opj_event_mgr_t * p_manager);
+
+static opj_bool j2k_update_image_data (opj_tcd_v2_t * p_tcd, OPJ_BYTE * p_data);
+
+
+/*
+ * -----------------------------------------------------------------------
+ * -----------------------------------------------------------------------
+ * -----------------------------------------------------------------------
+ */
+
 /**
 Write the SOC marker (Start Of Codestream)
 @param j2k J2K handle
@@ -48,6 +216,20 @@ Read the SOC marker (Start of Codestream)
 @param j2k J2K handle
 */
 static void j2k_read_soc(opj_j2k_t *j2k);
+
+/**
+ * Reads a SOC marker (Start of Codestream)
+ * @param      p_header_data   the data contained in the SOC box.
+ * @param      jp2                             the jpeg2000 file codec.
+ * @param      p_header_size   the size of the data contained in the SOC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_soc_v2(
+                                       opj_j2k_v2_t *p_j2k,
+                                       struct opj_stream_private *p_stream,
+                                       struct opj_event_mgr * p_manager
+                                );
+
 /**
 Write the SIZ marker (image and tile size)
 @param j2k J2K handle
@@ -58,6 +240,21 @@ Read the SIZ marker (image and tile size)
 @param j2k J2K handle
 */
 static void j2k_read_siz(opj_j2k_t *j2k);
+
+/**
+ * Reads a SIZ marker (image and tile size)
+ * @param      p_header_data   the data contained in the SIZ box.
+ * @param      jp2                             the jpeg2000 file codec.
+ * @param      p_header_size   the size of the data contained in the SIZ marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_siz_v2 (
+                                                 opj_j2k_v2_t *p_j2k,
+                                                 OPJ_BYTE * p_header_data,
+                                                 OPJ_UINT32 p_header_size,
+                                                 struct opj_event_mgr * p_manager
+                                       );
+
 /**
 Write the COM marker (comment)
 @param j2k J2K handle
@@ -68,6 +265,19 @@ Read the COM marker (comment)
 @param j2k J2K handle
 */
 static void j2k_read_com(opj_j2k_t *j2k);
+/**
+ * Reads a COM marker (comments)
+ * @param      p_header_data   the data contained in the COM box.
+ * @param      jp2                             the jpeg2000 file codec.
+ * @param      p_header_size   the size of the data contained in the COM marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_com_v2 (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       );
 /**
 Write the value concerning the specified component in the marker COD and COC
 @param j2k J2K handle
@@ -90,6 +300,21 @@ Read the COD marker (coding style default)
 @param j2k J2K handle
 */
 static void j2k_read_cod(opj_j2k_t *j2k);
+
+/**
+ * Reads a COD marker (Coding Styke defaults)
+ * @param      p_header_data   the data contained in the COD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the COD marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_cod_v2 (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       );
+
 /**
 Write the COC marker (coding style component)
 @param j2k J2K handle
@@ -101,6 +326,21 @@ Read the COC marker (coding style component)
 @param j2k J2K handle
 */
 static void j2k_read_coc(opj_j2k_t *j2k);
+
+/**
+ * Reads a COC marker (Coding Style Component)
+ * @param      p_header_data   the data contained in the COC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the COC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_coc_v2 (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       );
+
 /**
 Write the value concerning the specified component in the marker QCD and QCC
 @param j2k J2K handle
@@ -124,6 +364,21 @@ Read the QCD marker (quantization default)
 @param j2k J2K handle
 */
 static void j2k_read_qcd(opj_j2k_t *j2k);
+
+/**
+ * Reads a QCD marker (Quantization defaults)
+ * @param      p_header_data   the data contained in the QCD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the QCD marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_qcd_v2 (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       );
+
 /**
 Write the QCC marker (quantization component)
 @param j2k J2K handle
@@ -135,6 +390,19 @@ Read the QCC marker (quantization component)
 @param j2k J2K handle
 */
 static void j2k_read_qcc(opj_j2k_t *j2k);
+/**
+ * Reads a QCC marker (Quantization component)
+ * @param      p_header_data   the data contained in the QCC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the QCC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_qcc_v2(
+                                                       opj_j2k_v2_t *p_j2k,
+                                                       OPJ_BYTE * p_header_data,
+                                                       OPJ_UINT32 p_header_size,
+                                                       struct opj_event_mgr * p_manager);
+
 /**
 Write the POC marker (progression order change)
 @param j2k J2K handle
@@ -145,36 +413,143 @@ Read the POC marker (progression order change)
 @param j2k J2K handle
 */
 static void j2k_read_poc(opj_j2k_t *j2k);
+/**
+ * Reads a POC marker (Progression Order Change)
+ *
+ * @param      p_header_data   the data contained in the POC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the POC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_poc_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
 /**
 Read the CRG marker (component registration)
 @param j2k J2K handle
 */
 static void j2k_read_crg(opj_j2k_t *j2k);
+/**
+ * Reads a CRG marker (Component registration)
+ *
+ * @param      p_header_data   the data contained in the TLM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the TLM marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_crg_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
 /**
 Read the TLM marker (tile-part lengths)
 @param j2k J2K handle
 */
 static void j2k_read_tlm(opj_j2k_t *j2k);
+/**
+ * Reads a TLM marker (Tile Length Marker)
+ *
+ * @param      p_header_data   the data contained in the TLM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the TLM marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_tlm_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
 /**
 Read the PLM marker (packet length, main header)
 @param j2k J2K handle
 */
 static void j2k_read_plm(opj_j2k_t *j2k);
+
+/**
+ * Reads a PLM marker (Packet length, main header marker)
+ *
+ * @param      p_header_data   the data contained in the TLM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the TLM marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_plm_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
 /**
 Read the PLT marker (packet length, tile-part header)
 @param j2k J2K handle
 */
 static void j2k_read_plt(opj_j2k_t *j2k);
+/**
+ * Reads a PLT marker (Packet length, tile-part header)
+ *
+ * @param      p_header_data   the data contained in the PLT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the PLT marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_plt_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
 /**
 Read the PPM marker (packet packet headers, main header)
 @param j2k J2K handle
 */
 static void j2k_read_ppm(opj_j2k_t *j2k);
+/**
+ * Reads a PPM marker (Packed packet headers, main header)
+ *
+ * @param      p_header_data   the data contained in the POC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the POC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_ppm_v2 (
+                                               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_read_ppm_v3 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
+
 /**
 Read the PPT marker (packet packet headers, tile-part header)
 @param j2k J2K handle
 */
 static void j2k_read_ppt(opj_j2k_t *j2k);
+/**
+ * Reads a PPT marker (Packed packet headers, tile-part header)
+ *
+ * @param      p_header_data   the data contained in the PPT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the PPT marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_ppt_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
 /**
 Write the TLM marker (Mainheader)
 @param j2k J2K handle
@@ -190,6 +565,21 @@ Read the SOT marker (start of tile-part)
 @param j2k J2K handle
 */
 static void j2k_read_sot(opj_j2k_t *j2k);
+
+/**
+ * Reads a PPT marker (Packed packet headers, tile-part header)
+ *
+ * @param      p_header_data   the data contained in the PPT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the PPT marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_sot_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       );
 /**
 Write the SOD marker (start of data)
 @param j2k J2K handle
@@ -201,6 +591,21 @@ Read the SOD marker (start of data)
 @param j2k J2K handle
 */
 static void j2k_read_sod(opj_j2k_t *j2k);
+
+/**
+ * Reads a SOD marker (Start Of Data)
+ *
+ * @param      p_header_data   the data contained in the SOD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the SOD marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_sod_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               struct opj_stream_private *p_stream,
+                                               struct opj_event_mgr * p_manager
+                                       );
+
 /**
 Write the RGN marker (region-of-interest)
 @param j2k J2K handle
@@ -213,6 +618,22 @@ Read the RGN marker (region-of-interest)
 @param j2k J2K handle
 */
 static void j2k_read_rgn(opj_j2k_t *j2k);
+
+/**
+ * Reads a RGN marker (Region Of Interest)
+ *
+ * @param      p_header_data   the data contained in the POC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the POC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_rgn_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       ) ;
+
 /**
 Write the EOC marker (end of codestream)
 @param j2k J2K handle
@@ -223,11 +644,77 @@ Read the EOC marker (end of codestream)
 @param j2k J2K handle
 */
 static void j2k_read_eoc(opj_j2k_t *j2k);
+
+/**
+ * Reads a EOC marker (End Of Codestream)
+ *
+ * @param      p_header_data   the data contained in the SOD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the SOD marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_eoc_v2 (
+                                           opj_j2k_v2_t *p_j2k,
+                                               struct opj_stream_private *p_stream,
+                                               struct opj_event_mgr * p_manager
+                                       ) ;
+
 /**
 Read an unknown marker
 @param j2k J2K handle
 */
 static void j2k_read_unk(opj_j2k_t *j2k);
+/**
+Add main header marker information
+@param cstr_info Codestream information structure
+@param type marker type
+@param pos byte offset of marker segment
+@param len length of marker segment
+ */
+static void j2k_add_mhmarker(opj_codestream_info_t *cstr_info, unsigned short int type, int pos, int len);
+
+static void j2k_add_mhmarker_v2(opj_codestream_index_t *cstr_index, OPJ_UINT32 type, OPJ_UINT32 pos, OPJ_UINT32 len) ;
+/**
+Add tile header marker information
+@param tileno tile index number
+@param cstr_info Codestream information structure
+@param type marker type
+@param pos byte offset of marker segment
+@param len length of marker segment
+ */
+static void j2k_add_tlmarker( int tileno, opj_codestream_info_t *cstr_info, unsigned short int type, int pos, int len);
+
+
+/**
+ * Reads an unknown marker
+ *
+ * @param      p_stream                the stream object to read from.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_manager               the user event manager.
+ *
+ * @return     true                    if the marker could be deduced.
+*/
+static opj_bool j2k_read_unk_v2 (      opj_j2k_v2_t *p_j2k,
+                                                                       struct opj_stream_private *p_stream,
+                                                                       OPJ_UINT32 *output_marker,
+                                                                       struct opj_event_mgr * p_manager );
+
+/**
+ * Copy the image header from the jpeg2000 codec into an external image_header
+ *
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_image_header  the external empty image header
+ *
+ * @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 void j2k_dump_MH_info(opj_j2k_v2_t* p_j2k, FILE* out_stream);
+
+static void j2k_dump_MH_index(opj_j2k_v2_t* p_j2k, FILE* out_stream);
+
+static opj_codestream_index_t* j2k_create_cstr_index(void);
 
 /*@}*/
 
@@ -236,7 +723,7 @@ static void j2k_read_unk(opj_j2k_t *j2k);
 /* ----------------------------------------------------------------------- */
 typedef struct j2k_prog_order{
        OPJ_PROG_ORDER enum_prog;
-       char str_prog[4];
+       char str_prog[5];
 }j2k_prog_order_t;
 
 j2k_prog_order_t j2k_prog_order_list[] = {
@@ -245,9 +732,88 @@ j2k_prog_order_t j2k_prog_order_list[] = {
        {PCRL, "PCRL"},
        {RLCP, "RLCP"},
        {RPCL, "RPCL"},
-       {-1, ""}
+       {(OPJ_PROG_ORDER)-1, ""}
+};
+
+typedef struct opj_dec_memory_marker_handler
+{
+       /** marker value */
+       OPJ_UINT32 id;
+       /** value of the state when the marker can appear */
+       OPJ_UINT32 states;
+       /** action linked to the marker */
+       opj_bool (*handler) (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                               );
+}
+opj_dec_memory_marker_handler_t;
+
+const opj_dec_memory_marker_handler_t j2k_memory_marker_handler_tab [] =
+{
+#ifdef TODO_MS
+  {J2K_MS_SOT, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPHSOT, j2k_read_sot},
+  {J2K_MS_COD, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_cod},
+  {J2K_MS_COC, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_coc},
+  {J2K_MS_RGN, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_rgn},
+  {J2K_MS_QCD, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_qcd},
+  {J2K_MS_QCC, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_qcc},
+  {J2K_MS_POC, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_poc},
+  {J2K_MS_SIZ, J2K_DEC_STATE_MHSIZ , j2k_read_siz},
+  {J2K_MS_TLM, J2K_DEC_STATE_MH, j2k_read_tlm},
+  {J2K_MS_PLM, J2K_DEC_STATE_MH, j2k_read_plm},
+  {J2K_MS_PLT, J2K_DEC_STATE_TPH, j2k_read_plt},
+  {J2K_MS_PPM, J2K_DEC_STATE_MH, j2k_read_ppm},
+  {J2K_MS_PPT, J2K_DEC_STATE_TPH, j2k_read_ppt},
+  {J2K_MS_SOP, 0, 0},
+  {J2K_MS_CRG, J2K_DEC_STATE_MH, j2k_read_crg},
+  {J2K_MS_COM, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_com},
+  {J2K_MS_MCT, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_mct},
+  {J2K_MS_CBD, J2K_DEC_STATE_MH , j2k_read_cbd},
+  {J2K_MS_MCC, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_mcc},
+  {J2K_MS_MCO, J2K_DEC_STATE_MH | J2K_DEC_STATE_TPH, j2k_read_mco},
+#endif
+  {J2K_MS_SOT, J2K_STATE_MH | J2K_STATE_TPHSOT, j2k_read_sot_v2},
+  {J2K_MS_COD, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_cod_v2},
+  {J2K_MS_COC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_coc_v2},
+  {J2K_MS_RGN, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_rgn_v2},
+  {J2K_MS_QCD, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_qcd_v2},
+  {J2K_MS_QCC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_qcc_v2},
+  {J2K_MS_POC, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_poc_v2},
+  {J2K_MS_SIZ, J2K_STATE_MHSIZ , j2k_read_siz_v2},
+  {J2K_MS_TLM, J2K_STATE_MH, j2k_read_tlm_v2},
+  {J2K_MS_PLM, J2K_STATE_MH, j2k_read_plm_v2},
+  {J2K_MS_PLT, J2K_STATE_TPH, j2k_read_plt_v2},
+  {J2K_MS_PPM, J2K_STATE_MH, j2k_read_ppm_v3},
+  {J2K_MS_PPT, J2K_STATE_TPH, j2k_read_ppt_v2},
+  {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},
+  {J2K_MS_EPB, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_epb},
+  {J2K_MS_ESD, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_esd},
+  {J2K_MS_RED, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_red},
+#endif
+#endif /* USE_JPWL */
+#ifdef USE_JPSEC
+  {J2K_MS_SEC, J2K_DEC_STATE_MH, j2k_read_sec},
+  {J2K_MS_INSEC, 0, j2k_read_insec}
+#endif /* USE_JPSEC */
+  {J2K_MS_UNK, J2K_STATE_MH | J2K_STATE_TPH, 0}//j2k_read_unk_v2}
 };
 
+
+
 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++ ){
@@ -258,116 +824,6 @@ char *j2k_convert_progression_order(OPJ_PROG_ORDER prg_order){
        return po->str_prog;
 }
 
-static void j2k_check_poc_val(opj_cparameters_t *parameters, int numcomps, int numlayers){
-       int index, resno, compno, layno, i;
-       char loss = 0;
-       int step_c = 1;
-       int step_r = numcomps * step_c;
-       int step_l = parameters->numresolution * step_r;
-       int array_size = step_l * numlayers * sizeof(int);
-       int *packet_array = (int *) opj_malloc(array_size);
-       
-       for (i = 0; i < parameters->numpocs ; i++) {
-               int layno0 = 0;
-               if(i > 0)
-                       layno0 = (parameters->POC[i].layno1 > parameters->POC[i-1].layno1 )? parameters->POC[i-1].layno1 : 0;
-               for (resno = parameters->POC[i].resno0 ; resno < parameters->POC[i].resno1 ; resno++) {
-                       for (compno = parameters->POC[i].compno0 ; compno < parameters->POC[i].compno1 ; compno++) {
-                               for (layno = layno0; layno < parameters->POC[i].layno1 ; layno++) {
-                                       index = step_r * resno + step_c * compno + step_l * layno;
-                                       packet_array[index]= 1;
-                               }
-                       }
-               }
-       }
-       for (resno = 0; resno < parameters->numresolution; resno++) {
-               for (compno = 0; compno < numcomps; compno++) {
-                       for (layno = 0; layno < numlayers ; layno++) {
-                               index = step_r * resno + step_c * compno + step_l * layno;
-                               if(!(   packet_array[index]== 1)){
-                                       loss = 1;
-                               }
-                       }
-               }
-       }
-       if(loss == 1)
-               fprintf(stdout,"Missing packets possible loss of data\n");
-       opj_free(packet_array);
-}
-
-void j2k_dump_image(FILE *fd, opj_image_t * img) {
-       int compno;
-       fprintf(fd, "image {\n");
-       fprintf(fd, "  x0=%d, y0=%d, x1=%d, y1=%d\n", img->x0, img->y0, img->x1, img->y1);
-       fprintf(fd, "  numcomps=%d\n", img->numcomps);
-       for (compno = 0; compno < img->numcomps; compno++) {
-               opj_image_comp_t *comp = &img->comps[compno];
-               fprintf(fd, "  comp %d {\n", compno);
-               fprintf(fd, "    dx=%d, dy=%d\n", comp->dx, comp->dy);
-               fprintf(fd, "    prec=%d\n", comp->prec);
-               fprintf(fd, "    sgnd=%d\n", comp->sgnd);
-               fprintf(fd, "  }\n");
-       }
-       fprintf(fd, "}\n");
-}
-
-void j2k_dump_cp(FILE *fd, opj_image_t * img, opj_cp_t * cp) {
-       int tileno, compno, layno, bandno, resno, numbands;
-       fprintf(fd, "coding parameters {\n");
-       fprintf(fd, "  tx0=%d, ty0=%d\n", cp->tx0, cp->ty0);
-       fprintf(fd, "  tdx=%d, tdy=%d\n", cp->tdx, cp->tdy);
-       fprintf(fd, "  tw=%d, th=%d\n", cp->tw, cp->th);
-       for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
-               opj_tcp_t *tcp = &cp->tcps[tileno];
-               fprintf(fd, "  tile %d {\n", tileno);
-               fprintf(fd, "    csty=%x\n", tcp->csty);
-               fprintf(fd, "    prg=%d\n", tcp->prg);
-               fprintf(fd, "    numlayers=%d\n", tcp->numlayers);
-               fprintf(fd, "    mct=%d\n", tcp->mct);
-               fprintf(fd, "    rates=");
-               for (layno = 0; layno < tcp->numlayers; layno++) {
-                       fprintf(fd, "%.1f ", tcp->rates[layno]);
-               }
-               fprintf(fd, "\n");
-               for (compno = 0; compno < img->numcomps; compno++) {
-                       opj_tccp_t *tccp = &tcp->tccps[compno];
-                       fprintf(fd, "    comp %d {\n", compno);
-                       fprintf(fd, "      csty=%x\n", tccp->csty);
-                       fprintf(fd, "      numresolutions=%d\n", tccp->numresolutions);
-                       fprintf(fd, "      cblkw=%d\n", tccp->cblkw);
-                       fprintf(fd, "      cblkh=%d\n", tccp->cblkh);
-                       fprintf(fd, "      cblksty=%x\n", tccp->cblksty);
-                       fprintf(fd, "      qmfbid=%d\n", tccp->qmfbid);
-                       fprintf(fd, "      qntsty=%d\n", tccp->qntsty);
-                       fprintf(fd, "      numgbits=%d\n", tccp->numgbits);
-                       fprintf(fd, "      roishift=%d\n", tccp->roishift);
-                       fprintf(fd, "      stepsizes=");
-                       numbands = tccp->qntsty == J2K_CCP_QNTSTY_SIQNT ? 1 : tccp->numresolutions * 3 - 2;
-                       for (bandno = 0; bandno < numbands; bandno++) {
-                               fprintf(fd, "(%d,%d) ", tccp->stepsizes[bandno].mant,
-                                       tccp->stepsizes[bandno].expn);
-                       }
-                       fprintf(fd, "\n");
-                       
-                       if (tccp->csty & J2K_CCP_CSTY_PRT) {
-                               fprintf(fd, "      prcw=");
-                               for (resno = 0; resno < tccp->numresolutions; resno++) {
-                                       fprintf(fd, "%d ", tccp->prcw[resno]);
-                               }
-                               fprintf(fd, "\n");
-                               fprintf(fd, "      prch=");
-                               for (resno = 0; resno < tccp->numresolutions; resno++) {
-                                       fprintf(fd, "%d ", tccp->prch[resno]);
-                               }
-                               fprintf(fd, "\n");
-                       }
-                       fprintf(fd, "    }\n");
-               }
-               fprintf(fd, "  }\n");
-       }
-       fprintf(fd, "}\n");
-}
-
 /* ----------------------------------------------------------------------- */
 static int j2k_get_num_tp(opj_cp_t *cp,int pino,int tileno){
        char *prog;
@@ -406,7 +862,10 @@ static int j2k_get_num_tp(opj_cp_t *cp,int pino,int tileno){
 
 /**    mem allocation for TLM marker*/
 int j2k_calculate_tp(opj_cp_t *cp,int img_numcomp,opj_image_t *image,opj_j2k_t *j2k ){
-       int pino,tileno,maxres=0,totnum_tp=0;
+       int pino,tileno,totnum_tp=0;
+
+       OPJ_ARG_NOT_USED(img_numcomp);
+
        j2k->cur_totnum_tp = (int *) opj_malloc(cp->tw * cp->th * sizeof(int));
        for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
                int cur_totnum_tp = 0;
@@ -421,6 +880,12 @@ int j2k_calculate_tp(opj_cp_t *cp,int img_numcomp,opj_image_t *image,opj_j2k_t *
                        pi_destroy(pi, cp, tileno);
                }
                j2k->cur_totnum_tp[tileno] = cur_totnum_tp;
+               /* INDEX >> */
+               if (j2k->cstr_info) {
+                       j2k->cstr_info->tile[tileno].num_tps = cur_totnum_tp;
+                       j2k->cstr_info->tile[tileno].tp = (opj_tp_info_t *) opj_malloc(cur_totnum_tp * sizeof(opj_tp_info_t));
+               }
+               /* << INDEX */
        }
        return totnum_tp;
 }
@@ -428,10 +893,68 @@ int j2k_calculate_tp(opj_cp_t *cp,int img_numcomp,opj_image_t *image,opj_j2k_t *
 static void j2k_write_soc(opj_j2k_t *j2k) {
        opj_cio_t *cio = j2k->cio;
        cio_write(cio, J2K_MS_SOC, 2);
+
+       if(j2k->cstr_info)
+         j2k_add_mhmarker(j2k->cstr_info, J2K_MS_SOC, cio_tell(cio), 0);
+
+/* UniPG>> */
+#ifdef USE_JPWL
+
+       /* update markers struct */
+       j2k_add_marker(j2k->cstr_info, J2K_MS_SOC, cio_tell(cio) - 2, 2);
+#endif /* USE_JPWL */
+/* <<UniPG */
 }
 
-static void j2k_read_soc(opj_j2k_t *j2k) {
+static void j2k_read_soc(opj_j2k_t *j2k) {     
        j2k->state = J2K_STATE_MHSIZ;
+       /* Index */
+       if (j2k->cstr_info) {
+               j2k->cstr_info->main_head_start = cio_tell(j2k->cio) - 2;
+               j2k->cstr_info->codestream_size = cio_numbytesleft(j2k->cio) + 2 - j2k->cstr_info->main_head_start;
+       }
+}
+
+/**
+ * Reads a SOC marker (Start of Codestream)
+ * @param      p_header_data   the data contained in the SOC box.
+ * @param      jp2                             the jpeg2000 file codec.
+ * @param      p_header_size   the size of the data contained in the SOC marker.
+ * @param      p_manager               the user event manager.
+*/
+static opj_bool j2k_read_soc_v2(       opj_j2k_v2_t *p_j2k,
+                                                                       struct opj_stream_private *p_stream,
+                                                                       struct opj_event_mgr * p_manager )
+{
+       OPJ_BYTE l_data [2];
+       OPJ_UINT32 l_marker;
+
+       /* preconditions */
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       assert(p_stream != 00);
+
+       if (opj_stream_read_data(p_stream,l_data,2,p_manager) != 2) {
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(l_data,&l_marker,2);
+       if (l_marker != J2K_MS_SOC) {
+               return OPJ_FALSE;
+       }
+
+       /* Next marker should be a SIZ marker in the main header */
+       p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_MHSIZ;
+
+       /* FIXME move it in a index structure included in p_j2k*/
+       p_j2k->cstr_index->main_head_start = (OPJ_UINT32) opj_stream_tell(p_stream) - 2;
+
+       opj_event_msg_v2(p_manager, EVT_INFO, "Start to read j2k main header (%d).\n", p_j2k->cstr_index->main_head_start);
+
+       /* Add the marker to the codestream index*/
+       j2k_add_mhmarker_v2(p_j2k->cstr_index, J2K_MS_SOC, p_j2k->cstr_index->main_head_start, 2);
+
+       return OPJ_TRUE;
 }
 
 static void j2k_write_siz(opj_j2k_t *j2k) {
@@ -441,7 +964,7 @@ static void j2k_write_siz(opj_j2k_t *j2k) {
        opj_cio_t *cio = j2k->cio;
        opj_image_t *image = j2k->image;
        opj_cp_t *cp = j2k->cp;
-       
+
        cio_write(cio, J2K_MS_SIZ, 2);  /* SIZ */
        lenp = cio_tell(cio);
        cio_skip(cio, 2);
@@ -464,6 +987,9 @@ static void j2k_write_siz(opj_j2k_t *j2k) {
        cio_seek(cio, lenp);
        cio_write(cio, len, 2);         /* Lsiz */
        cio_seek(cio, lenp + len);
+       
+       if(j2k->cstr_info)
+         j2k_add_mhmarker(j2k->cstr_info, J2K_MS_SIZ, lenp, len);
 }
 
 static void j2k_read_siz(opj_j2k_t *j2k) {
@@ -484,6 +1010,13 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
        cp->tx0 = cio_read(cio, 4);             /* XT0siz */
        cp->ty0 = cio_read(cio, 4);             /* YT0siz */
        
+       if ((image->x0<0)||(image->x1<0)||(image->y0<0)||(image->y1<0)) {
+               opj_event_msg(j2k->cinfo, EVT_ERROR,
+                                                                       "%s: invalid image size (x0:%d, x1:%d, y0:%d, y1:%d)\n",
+                                                                       image->x0,image->x1,image->y0,image->y1);
+               return;
+       }
+       
        image->numcomps = cio_read(cio, 2);     /* Csiz */
 
 #ifdef USE_JPWL
@@ -525,7 +1058,7 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
        }
 #endif /* USE_JPWL */
 
-       image->comps = (opj_image_comp_t *) opj_malloc(image->numcomps * sizeof(opj_image_comp_t));
+       image->comps = (opj_image_comp_t*) opj_calloc(image->numcomps, sizeof(opj_image_comp_t));
        for (i = 0; i < image->numcomps; i++) {
                int tmp, w, h;
                tmp = cio_read(cio, 1);         /* Ssiz_i */
@@ -563,13 +1096,12 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
                }
 #endif /* USE_JPWL */
 
-
                /* TODO: unused ? */
                w = int_ceildiv(image->x1 - image->x0, image->comps[i].dx);
                h = int_ceildiv(image->y1 - image->y0, image->comps[i].dy);
 
                image->comps[i].resno_decoded = 0;      /* number of resolution decoded */
-               image->comps[i].factor = 0;                     /* reducing factor per component */
+               image->comps[i].factor = cp->reduce; /* reducing factor per component */
        }
        
        cp->tw = int_ceildiv(image->x1 - cp->tx0, cp->tdx);
@@ -615,8 +1147,8 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
        }
 #endif /* USE_JPWL */
 
-       cp->tcps = (opj_tcp_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tcp_t));
-       cp->tileno = (int *) opj_malloc(cp->tw * cp->th * sizeof(int));
+       cp->tcps = (opj_tcp_t*) opj_calloc(cp->tw * cp->th, sizeof(opj_tcp_t));
+       cp->tileno = (int*) opj_malloc(cp->tw * cp->th * sizeof(int));
        cp->tileno_size = 0;
        
 #ifdef USE_JPWL
@@ -644,39 +1176,377 @@ static void j2k_read_siz(opj_j2k_t *j2k) {
        cp->ppm_data_first = NULL;
        cp->ppm_previous = 0;
        cp->ppm_store = 0;
-       
-       j2k->default_tcp->tccps = (opj_tccp_t *) opj_malloc(sizeof(opj_tccp_t) * image->numcomps);
+
+       j2k->default_tcp->tccps = (opj_tccp_t*) opj_calloc(image->numcomps, sizeof(opj_tccp_t));
        for (i = 0; i < cp->tw * cp->th; i++) {
-               cp->tcps[i].tccps = (opj_tccp_t *) opj_malloc(sizeof(opj_tccp_t) * image->numcomps);
-       }
-       j2k->tile_data = (unsigned char **) opj_malloc(cp->tw * cp->th * sizeof(unsigned char *));
-       j2k->tile_len = (int *) opj_malloc(cp->tw * cp->th * sizeof(int));
+               cp->tcps[i].tccps = (opj_tccp_t*) opj_malloc(image->numcomps * sizeof(opj_tccp_t));
+       }       
+       j2k->tile_data = (unsigned char**) opj_calloc(cp->tw * cp->th, sizeof(unsigned char*));
+       j2k->tile_len = (int*) opj_calloc(cp->tw * cp->th, sizeof(int));
        j2k->state = J2K_STATE_MH;
+
+       /* Index */
+       if (j2k->cstr_info) {
+               opj_codestream_info_t *cstr_info = j2k->cstr_info;
+               cstr_info->image_w = image->x1 - image->x0;
+               cstr_info->image_h = image->y1 - image->y0;
+               cstr_info->numcomps = image->numcomps;
+               cstr_info->tw = cp->tw;
+               cstr_info->th = cp->th;
+               cstr_info->tile_x = cp->tdx;    
+               cstr_info->tile_y = cp->tdy;    
+               cstr_info->tile_Ox = cp->tx0;   
+               cstr_info->tile_Oy = cp->ty0;                   
+               cstr_info->tile = (opj_tile_info_t*) opj_calloc(cp->tw * cp->th, sizeof(opj_tile_info_t));              
+       }
 }
 
-static void j2k_write_com(opj_j2k_t *j2k) {
-       unsigned int i;
-       int lenp, len;
 
-       if(j2k->cp->comment) {
-               opj_cio_t *cio = j2k->cio;
-               char *comment = j2k->cp->comment;
+/**
+ * Reads a SIZ marker (image and tile size)
+ * @param      p_header_data   the data contained in the SIZ box.
+ * @param      jp2                             the jpeg2000 file codec.
+ * @param      p_header_size   the size of the data contained in the SIZ marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_siz_v2 (
+                                   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_size, i;
+       OPJ_UINT32 l_nb_comp;
+       OPJ_UINT32 l_nb_comp_remain;
+       OPJ_UINT32 l_remaining_size;
+       OPJ_UINT32 l_nb_tiles;
+       OPJ_UINT32 l_tmp;
+       opj_image_t *l_image = 00;
+       opj_cp_v2_t *l_cp = 00;
+       opj_image_comp_t * l_img_comp = 00;
+       opj_tcp_v2_t * l_current_tile_param = 00;
+
+       // preconditions
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       assert(p_header_data != 00);
+
+       l_image = p_j2k->m_image;
+       l_cp = &(p_j2k->m_cp);
+
+       // minimum size == 39 - 3 (= minimum component parameter)
+       if (p_header_size < 36) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error with SIZ marker size\n");
+               return OPJ_FALSE;
+       }
 
-               cio_write(cio, J2K_MS_COM, 2);
-               lenp = cio_tell(cio);
-               cio_skip(cio, 2);
-               cio_write(cio, 0, 2);
-               for (i = 0; i < strlen(comment); i++) {
-                       cio_write(cio, comment[i], 1);
-               }
-               len = cio_tell(cio) - lenp;
-               cio_seek(cio, lenp);
-               cio_write(cio, len, 2);
-               cio_seek(cio, lenp + len);
+       l_remaining_size = p_header_size - 36;
+       l_nb_comp = l_remaining_size / 3;
+       l_nb_comp_remain = l_remaining_size % 3;
+       if (l_nb_comp_remain != 0){
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error with SIZ marker size\n");
+               return OPJ_FALSE;
        }
-}
 
-static void j2k_read_com(opj_j2k_t *j2k) {
+       l_size = p_header_size + 2;                                                                             /* Lsiz */
+
+       opj_read_bytes(p_header_data,&l_tmp ,2);                                                /* Rsiz (capabilities) */
+       p_header_data+=2;
+       l_cp->rsiz = (OPJ_RSIZ_CAPABILITIES) l_tmp;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_image->x1, 4);   /* Xsiz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_image->y1, 4);   /* Ysiz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_image->x0, 4);   /* X0siz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_image->y0, 4);   /* Y0siz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_cp->tdx, 4);             /* XTsiz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_cp->tdy, 4);             /* YTsiz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_cp->tx0, 4);             /* XT0siz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_cp->ty0, 4);             /* YT0siz */
+       p_header_data+=4;
+       opj_read_bytes(p_header_data, (OPJ_UINT32*) &l_tmp, 2);                 /* Csiz */
+       p_header_data+=2;
+       if (l_tmp < 16385)
+               l_image->numcomps = (OPJ_UINT16) l_tmp;
+       else {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error with SIZ marker: number of component is illegal -> %d\n", l_tmp);
+               return OPJ_FALSE;
+       }
+
+       if (l_image->numcomps != l_nb_comp) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error with SIZ marker: number of component is not compatible with the remaining number of parameters ( %d vs %d)\n", l_image->numcomps, l_nb_comp);
+               return OPJ_FALSE;
+       }
+
+#ifdef USE_JPWL
+       if (l_cp->correct) {
+               /* if JPWL is on, we check whether TX errors have damaged
+                 too much the SIZ parameters */
+               if (!(l_image->x1 * l_image->y1)) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR,
+                               "JPWL: bad image size (%d x %d)\n",
+                               l_image->x1, l_image->y1);
+                       if (!JPWL_ASSUME || JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+               }
+
+       /* FIXME check previously in the function so why keep this piece of code ? Need by the norm ?
+               if (l_image->numcomps != ((len - 38) / 3)) {
+                       opj_event_msg_v2(p_manager, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
+                               "JPWL: Csiz is %d => space in SIZ only for %d comps.!!!\n",
+                               l_image->numcomps, ((len - 38) / 3));
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+       */              /* we try to correct */
+       /*              opj_event_msg_v2(p_manager, EVT_WARNING, "- trying to adjust this\n");
+                       if (l_image->numcomps < ((len - 38) / 3)) {
+                               len = 38 + 3 * l_image->numcomps;
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "- setting Lsiz to %d => HYPOTHESIS!!!\n",
+                                       len);
+                       } else {
+                               l_image->numcomps = ((len - 38) / 3);
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "- setting Csiz to %d => HYPOTHESIS!!!\n",
+                                       l_image->numcomps);
+                       }
+               }
+       */
+
+               /* update components number in the jpwl_exp_comps filed */
+               l_cp->exp_comps = l_image->numcomps;
+       }
+#endif /* USE_JPWL */
+
+       // Allocate the resulting image components
+       l_image->comps = (opj_image_comp_t*) opj_calloc(l_image->numcomps, sizeof(opj_image_comp_t));
+       if (l_image->comps == 00){
+               l_image->numcomps = 0;
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory to take in charge SIZ marker\n");
+               return OPJ_FALSE;
+       }
+
+       memset(l_image->comps,0,l_image->numcomps * sizeof(opj_image_comp_t));
+       l_img_comp = l_image->comps;
+
+       // Read the component information
+       for (i = 0; i < l_image->numcomps; ++i){
+               OPJ_UINT32 tmp;
+               opj_read_bytes(p_header_data,&tmp,1);   /* Ssiz_i */
+               ++p_header_data;
+               l_img_comp->prec = (tmp & 0x7f) + 1;
+               l_img_comp->sgnd = tmp >> 7;
+               opj_read_bytes(p_header_data,&tmp,1);   /* XRsiz_i */
+               ++p_header_data;
+               l_img_comp->dx = (OPJ_INT32)tmp; // should be between 1 and 255
+               opj_read_bytes(p_header_data,&tmp,1);   /* YRsiz_i */
+               ++p_header_data;
+               l_img_comp->dy = (OPJ_INT32)tmp; // should be between 1 and 255
+
+#ifdef USE_JPWL
+               if (l_cp->correct) {
+               /* if JPWL is on, we check whether TX errors have damaged
+                       too much the SIZ parameters, again */
+                       if (!(l_image->comps[i].dx * l_image->comps[i].dy)) {
+                               opj_event_msg_v2(p_manager, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
+                                       "JPWL: bad XRsiz_%d/YRsiz_%d (%d x %d)\n",
+                                       i, i, l_image->comps[i].dx, l_image->comps[i].dy);
+                               if (!JPWL_ASSUME) {
+                                       opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                                       return OPJ_FALSE;
+                               }
+                               /* we try to correct */
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "- trying to adjust them\n");
+                               if (!l_image->comps[i].dx) {
+                                       l_image->comps[i].dx = 1;
+                                       opj_event_msg_v2(p_manager, EVT_WARNING, "- setting XRsiz_%d to %d => HYPOTHESIS!!!\n",
+                                               i, l_image->comps[i].dx);
+                               }
+                               if (!l_image->comps[i].dy) {
+                                       l_image->comps[i].dy = 1;
+                                       opj_event_msg_v2(p_manager, EVT_WARNING, "- setting YRsiz_%d to %d => HYPOTHESIS!!!\n",
+                                               i, l_image->comps[i].dy);
+                               }
+                       }
+               }
+#endif /* USE_JPWL */
+               l_img_comp->resno_decoded = 0;                                                          /* number of resolution decoded */
+               l_img_comp->factor = l_cp->m_specific_param.m_dec.m_reduce; /* reducing factor per component */
+               ++l_img_comp;
+       }
+
+       // Compute the number of tiles
+       l_cp->tw = int_ceildiv(l_image->x1 - l_cp->tx0, l_cp->tdx);
+       l_cp->th = int_ceildiv(l_image->y1 - l_cp->ty0, l_cp->tdy);
+       l_nb_tiles = l_cp->tw * l_cp->th;
+
+       // Define the tiles which will be decoded
+       if (p_j2k->m_specific_param.m_decoder.m_discard_tiles) {
+               p_j2k->m_specific_param.m_decoder.m_start_tile_x = (p_j2k->m_specific_param.m_decoder.m_start_tile_x - l_cp->tx0) / l_cp->tdx;
+               p_j2k->m_specific_param.m_decoder.m_start_tile_y = (p_j2k->m_specific_param.m_decoder.m_start_tile_y - l_cp->ty0) / l_cp->tdy;
+               p_j2k->m_specific_param.m_decoder.m_end_tile_x = int_ceildiv((p_j2k->m_specific_param.m_decoder.m_end_tile_x - l_cp->tx0), l_cp->tdx);
+               p_j2k->m_specific_param.m_decoder.m_end_tile_y = int_ceildiv((p_j2k->m_specific_param.m_decoder.m_end_tile_y - l_cp->ty0), l_cp->tdy);
+       }
+       else {
+               p_j2k->m_specific_param.m_decoder.m_start_tile_x = 0;
+               p_j2k->m_specific_param.m_decoder.m_start_tile_y = 0;
+               p_j2k->m_specific_param.m_decoder.m_end_tile_x = l_cp->tw;
+               p_j2k->m_specific_param.m_decoder.m_end_tile_y = l_cp->th;
+       }
+
+#ifdef USE_JPWL
+       if (l_cp->correct) {
+               /* if JPWL is on, we check whether TX errors have damaged
+                 too much the SIZ parameters */
+               if ((l_cp->tw < 1) || (l_cp->th < 1) || (l_cp->tw > l_cp->max_tiles) || (l_cp->th > l_cp->max_tiles)) {
+                       opj_event_msg_v2(p_manager, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
+                               "JPWL: bad number of tiles (%d x %d)\n",
+                               l_cp->tw, l_cp->th);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+                       /* we try to correct */
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "- trying to adjust them\n");
+                       if (l_cp->tw < 1) {
+                               l_cp->tw= 1;
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "- setting %d tiles in x => HYPOTHESIS!!!\n",
+                                               l_cp->tw);
+                       }
+                       if (l_cp->tw > l_cp->max_tiles) {
+                               l_cp->tw= 1;
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "- too large x, increase expectance of %d\n"
+                                       "- setting %d tiles in x => HYPOTHESIS!!!\n",
+                                       l_cp->max_tiles, l_cp->tw);
+                       }
+                       if (l_cp->th < 1) {
+                               l_cp->th= 1;
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "- setting %d tiles in y => HYPOTHESIS!!!\n",
+                                               l_cp->th);
+                       }
+                       if (l_cp->th > l_cp->max_tiles) {
+                               l_cp->th= 1;
+                               opj_event_msg_v2(p_manager, EVT_WARNING, "- too large y, increase expectance of %d to continue\n",
+                                       "- setting %d tiles in y => HYPOTHESIS!!!\n",
+                                       l_cp->max_tiles, l_cp->th);
+                       }
+               }
+       }
+#endif /* USE_JPWL */
+
+       /* memory allocations */
+       l_cp->tcps = (opj_tcp_v2_t*) opj_calloc(l_nb_tiles, sizeof(opj_tcp_v2_t));
+       if (l_cp->tcps == 00) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory to take in charge SIZ marker\n");
+               return OPJ_FALSE;
+       }
+       memset(l_cp->tcps,0,l_nb_tiles*sizeof(opj_tcp_t));
+
+#ifdef USE_JPWL
+       if (l_cp->correct) {
+               if (!l_cp->tcps) {
+                       opj_event_msg_v2(p_manager, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
+                               "JPWL: could not alloc tcps field of cp\n");
+                       if (!JPWL_ASSUME || JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+               }
+       }
+#endif /* USE_JPWL */
+
+       p_j2k->m_specific_param.m_decoder.m_default_tcp->tccps =
+                       (opj_tccp_t*) opj_calloc(l_image->numcomps, sizeof(opj_tccp_t));
+       if(p_j2k->m_specific_param.m_decoder.m_default_tcp->tccps  == 00) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory to take in charge SIZ marker\n");
+               return OPJ_FALSE;
+       }
+       memset(p_j2k->m_specific_param.m_decoder.m_default_tcp->tccps ,0,l_image->numcomps*sizeof(opj_tccp_t));
+
+       p_j2k->m_specific_param.m_decoder.m_default_tcp->m_mct_records =
+                       (opj_mct_data_t*)opj_malloc(J2K_MCT_DEFAULT_NB_RECORDS * sizeof(opj_mct_data_t));
+
+       if (! p_j2k->m_specific_param.m_decoder.m_default_tcp->m_mct_records) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory to take in charge SIZ marker\n");
+               return OPJ_FALSE;
+       }
+       memset(p_j2k->m_specific_param.m_decoder.m_default_tcp->m_mct_records,0,J2K_MCT_DEFAULT_NB_RECORDS * sizeof(opj_mct_data_t));
+       p_j2k->m_specific_param.m_decoder.m_default_tcp->m_nb_max_mct_records = J2K_MCT_DEFAULT_NB_RECORDS;
+
+       p_j2k->m_specific_param.m_decoder.m_default_tcp->m_mcc_records =
+                       (opj_simple_mcc_decorrelation_data_t*)
+                       opj_malloc(J2K_MCC_DEFAULT_NB_RECORDS * sizeof(opj_simple_mcc_decorrelation_data_t));
+
+       if (! p_j2k->m_specific_param.m_decoder.m_default_tcp->m_mcc_records) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory to take in charge SIZ marker\n");
+               return OPJ_FALSE;
+       }
+       memset(p_j2k->m_specific_param.m_decoder.m_default_tcp->m_mcc_records,0,J2K_MCC_DEFAULT_NB_RECORDS * sizeof(opj_simple_mcc_decorrelation_data_t));
+       p_j2k->m_specific_param.m_decoder.m_default_tcp->m_nb_max_mcc_records = J2K_MCC_DEFAULT_NB_RECORDS;
+
+       /* set up default dc level shift */
+       for (i=0;i<l_image->numcomps;++i) {
+               if (! l_image->comps[i].sgnd) {
+                       p_j2k->m_specific_param.m_decoder.m_default_tcp->tccps[i].m_dc_level_shift = 1 << (l_image->comps[i].prec - 1);
+               }
+       }
+
+       l_current_tile_param = l_cp->tcps;
+       for     (i = 0; i < l_nb_tiles; ++i) {
+               l_current_tile_param->tccps = (opj_tccp_t*) opj_malloc(l_image->numcomps * sizeof(opj_tccp_t));
+               if (l_current_tile_param->tccps == 00) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory to take in charge SIZ marker\n");
+                       return OPJ_FALSE;
+               }
+               memset(l_current_tile_param->tccps,0,l_image->numcomps * sizeof(opj_tccp_t));
+
+               ++l_current_tile_param;
+       }
+
+       p_j2k->m_specific_param.m_decoder.m_state =  J2K_STATE_MH; // FIXME J2K_DEC_STATE_MH;
+       opj_image_comp_header_update(l_image,l_cp);
+
+       return OPJ_TRUE;
+}
+
+
+
+static void j2k_write_com(opj_j2k_t *j2k) {
+       unsigned int i;
+       int lenp, len;
+
+       if(j2k->cp->comment) {
+               opj_cio_t *cio = j2k->cio;
+               char *comment = j2k->cp->comment;
+
+               cio_write(cio, J2K_MS_COM, 2);
+               lenp = cio_tell(cio);
+               cio_skip(cio, 2);
+               cio_write(cio, 1, 2);           /* General use (IS 8859-15:1999 (Latin) values) */
+               for (i = 0; i < strlen(comment); i++) {
+                       cio_write(cio, comment[i], 1);
+               }
+               len = cio_tell(cio) - lenp;
+               cio_seek(cio, lenp);
+               cio_write(cio, len, 2);
+               cio_seek(cio, lenp + len);
+
+               
+               if(j2k->cstr_info)
+                 j2k_add_mhmarker(j2k->cstr_info, J2K_MS_COM, lenp, len);
+
+       }
+}
+
+static void j2k_read_com(opj_j2k_t *j2k) {
        int len;
        
        opj_cio_t *cio = j2k->cio;
@@ -684,6 +1554,27 @@ static void j2k_read_com(opj_j2k_t *j2k) {
        len = cio_read(cio, 2);
        cio_skip(cio, len - 2);  
 }
+/**
+ * Reads a COM marker (comments)
+ * @param      p_header_data   the data contained in the COM box.
+ * @param      jp2                             the jpeg2000 file codec.
+ * @param      p_header_size   the size of the data contained in the COM marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_com_v2 (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       )
+{
+       // preconditions
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       assert(p_header_data != 00);
+
+       return OPJ_TRUE;
+}
 
 static void j2k_write_cox(opj_j2k_t *j2k, int compno) {
        int i;
@@ -716,8 +1607,13 @@ static void j2k_read_cox(opj_j2k_t *j2k, int compno) {
 
        tccp->numresolutions = cio_read(cio, 1) + 1;    /* SPcox (D) */
 
-       /* check the reduce value */
-       cp->reduce = int_min((tccp->numresolutions)-1, cp->reduce);
+       // If user wants to remove more resolutions than the codestream contains, return error
+       if (cp->reduce >= tccp->numresolutions) {
+               opj_event_msg(j2k->cinfo, EVT_ERROR, "Error decoding component %d.\nThe number of resolutions to remove is higher than the number "
+                                       "of resolutions of this component\nModify the cp_reduce parameter.\n\n", compno);
+               j2k->state |= J2K_STATE_ERR;
+       }
+
        tccp->cblkw = cio_read(cio, 1) + 2;     /* SPcox (E) */
        tccp->cblkh = cio_read(cio, 1) + 2;     /* SPcox (F) */
        tccp->cblksty = cio_read(cio, 1);       /* SPcox (G) */
@@ -729,6 +1625,21 @@ static void j2k_read_cox(opj_j2k_t *j2k, int compno) {
                        tccp->prch[i] = tmp >> 4;
                }
        }
+
+       /* INDEX >> */
+       if(j2k->cstr_info && compno == 0) {
+               for (i = 0; i < tccp->numresolutions; i++) {
+                       if (tccp->csty & J2K_CP_CSTY_PRT) {
+                               j2k->cstr_info->tile[j2k->curtileno].pdx[i] = tccp->prcw[i];
+                               j2k->cstr_info->tile[j2k->curtileno].pdy[i] = tccp->prch[i];
+                       }
+                       else {
+                               j2k->cstr_info->tile[j2k->curtileno].pdx[i] = 15;
+                               j2k->cstr_info->tile[j2k->curtileno].pdx[i] = 15;
+                       }
+               }
+       }
+       /* << INDEX */
 }
 
 static void j2k_write_cod(opj_j2k_t *j2k) {
@@ -756,6 +1667,10 @@ static void j2k_write_cod(opj_j2k_t *j2k) {
        cio_seek(cio, lenp);
        cio_write(cio, len, 2);         /* Lcod */
        cio_seek(cio, lenp + len);
+
+       if(j2k->cstr_info)
+         j2k_add_mhmarker(j2k->cstr_info, J2K_MS_COD, lenp, len);
+
 }
 
 static void j2k_read_cod(opj_j2k_t *j2k) {
@@ -778,6 +1693,110 @@ static void j2k_read_cod(opj_j2k_t *j2k) {
                cio_seek(cio, pos);
                j2k_read_cox(j2k, i);
        }
+
+       /* Index */
+       if (j2k->cstr_info) {
+               opj_codestream_info_t *cstr_info = j2k->cstr_info;
+               cstr_info->prog = tcp->prg;
+               cstr_info->numlayers = tcp->numlayers;
+               cstr_info->numdecompos = (int*) opj_malloc(image->numcomps * sizeof(int));
+               for (i = 0; i < image->numcomps; i++) {
+                       cstr_info->numdecompos[i] = tcp->tccps[i].numresolutions - 1;
+               }
+       }
+}
+
+/**
+ * Reads a COD marker (Coding Styke defaults)
+ * @param      p_header_data   the data contained in the COD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the COD marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_cod_v2 (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       )
+{
+       // loop
+       OPJ_UINT32 i;
+       OPJ_UINT32 l_tmp;
+       opj_cp_v2_t *l_cp = 00;
+       opj_tcp_v2_t *l_tcp = 00;
+       opj_image_t *l_image = 00;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_image = p_j2k->m_image;
+       l_cp = &(p_j2k->m_cp);
+
+       /* If we are in the first tile-part header of the current tile */
+       l_tcp = (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH) ?
+                               &l_cp->tcps[p_j2k->m_current_tile_number] :
+                               p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       /* Make sure room is sufficient */
+       if (p_header_size < 5) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading COD marker\n");
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(p_header_data,&l_tcp->csty,1);           /* Scod */
+       ++p_header_data;
+       opj_read_bytes(p_header_data,&l_tmp,1);                         /* SGcod (A) */
+       ++p_header_data;
+       l_tcp->prg = (OPJ_PROG_ORDER) l_tmp;
+       opj_read_bytes(p_header_data,&l_tcp->numlayers,2);      /* SGcod (B) */
+       p_header_data+=2;
+
+       // If user didn't set a number layer to decode take the max specify in the codestream.
+       if      (l_cp->m_specific_param.m_dec.m_layer) {
+               l_tcp->num_layers_to_decode = l_cp->m_specific_param.m_dec.m_layer;
+       }
+       else {
+               l_tcp->num_layers_to_decode = l_tcp->numlayers;
+       }
+
+       opj_read_bytes(p_header_data,&l_tcp->mct,1);            /* SGcod (C) */
+       ++p_header_data;
+
+       p_header_size -= 5;
+       for     (i = 0; i < l_image->numcomps; ++i) {
+               l_tcp->tccps[i].csty = l_tcp->csty & J2K_CCP_CSTY_PRT;
+       }
+
+       if (! j2k_read_SPCod_SPCoc(p_j2k,0,p_header_data,&p_header_size,p_manager)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading COD marker\n");
+               return OPJ_FALSE;
+       }
+
+       if (p_header_size != 0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading COD marker\n");
+               return OPJ_FALSE;
+       }
+
+       /* Apply the coding style to other components of the current tile or the m_default_tcp*/
+       j2k_copy_tile_component_parameters(p_j2k);
+
+       /* Index */
+#ifdef WIP_REMOVE_MSD
+       if (p_j2k->cstr_info) {
+               //opj_codestream_info_t *l_cstr_info = p_j2k->cstr_info;
+               p_j2k->cstr_info->prog = l_tcp->prg;
+               p_j2k->cstr_info->numlayers = l_tcp->numlayers;
+               p_j2k->cstr_info->numdecompos = (OPJ_INT32*) opj_malloc(l_image->numcomps * sizeof(OPJ_UINT32));
+               for     (i = 0; i < l_image->numcomps; ++i) {
+                       p_j2k->cstr_info->numdecompos[i] = l_tcp->tccps[i].numresolutions - 1;
+               }
+       }
+#endif
+
+       return OPJ_TRUE;
 }
 
 static void j2k_write_coc(opj_j2k_t *j2k, int compno) {
@@ -814,6 +1833,68 @@ static void j2k_read_coc(opj_j2k_t *j2k) {
        j2k_read_cox(j2k, compno);
 }
 
+/**
+ * Reads a COC marker (Coding Style Component)
+ * @param      p_header_data   the data contained in the COC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the COC marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_coc_v2 (
+                                       opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       )
+{
+       opj_cp_v2_t *l_cp = NULL;
+       opj_tcp_v2_t *l_tcp = NULL;
+       opj_image_t *l_image = NULL;
+       OPJ_UINT32 l_comp_room;
+       OPJ_UINT32 l_comp_no;
+
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_cp = &(p_j2k->m_cp);
+       l_tcp = (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH ) ? /*FIXME J2K_DEC_STATE_TPH*/
+                               &l_cp->tcps[p_j2k->m_current_tile_number] :
+                               p_j2k->m_specific_param.m_decoder.m_default_tcp;
+       l_image = p_j2k->m_image;
+
+       l_comp_room = l_image->numcomps <= 256 ? 1 : 2;
+
+       // make sure room is sufficient
+       if (p_header_size < l_comp_room + 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading COC marker\n");
+               return OPJ_FALSE;
+       }
+       p_header_size -= l_comp_room + 1;
+
+       opj_read_bytes(p_header_data,&l_comp_no,l_comp_room);                   /* Ccoc */
+       p_header_data += l_comp_room;
+       if (l_comp_no >= l_image->numcomps) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading COC marker (bad number of components)\n");
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(p_header_data,&l_tcp->tccps[l_comp_no].csty,1);                  /* Scoc */
+       ++p_header_data ;
+
+       if (! j2k_read_SPCod_SPCoc(p_j2k,l_comp_no,p_header_data,&p_header_size,p_manager)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading COC marker\n");
+               return OPJ_FALSE;
+       }
+
+       if (p_header_size != 0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading COC marker\n");
+               return OPJ_FALSE;
+       }
+       return OPJ_TRUE;
+}
+
 static void j2k_write_qcx(opj_j2k_t *j2k, int compno) {
        int bandno, numbands;
        int expn, mant;
@@ -873,6 +1954,15 @@ static void j2k_read_qcx(opj_j2k_t *j2k, int compno, int len) {
                };
 
        };
+
+#else
+       /* We check whether there are too many subbands */
+       if ((numbands < 0) || (numbands >= J2K_MAXBANDS)) {
+               opj_event_msg(j2k->cinfo, EVT_WARNING ,
+                                       "bad number of subbands in Sqcx (%d) regarding to J2K_MAXBANDS (%d) \n"
+                                   "- limiting number of bands to J2K_MAXBANDS and try to move to the next markers\n", numbands, J2K_MAXBANDS);
+       }
+
 #endif /* USE_JPWL */
 
        for (bandno = 0; bandno < numbands; bandno++) {
@@ -885,8 +1975,10 @@ static void j2k_read_qcx(opj_j2k_t *j2k, int compno, int len) {
                        expn = tmp >> 11;
                        mant = tmp & 0x7ff;
                }
-               tccp->stepsizes[bandno].expn = expn;
-               tccp->stepsizes[bandno].mant = mant;
+               if (bandno < J2K_MAXBANDS){
+                       tccp->stepsizes[bandno].expn = expn;
+                       tccp->stepsizes[bandno].mant = mant;
+               }
        }
        
        /* Add Antonin : if scalar_derived -> compute other stepsizes */
@@ -914,6 +2006,9 @@ static void j2k_write_qcd(opj_j2k_t *j2k) {
        cio_seek(cio, lenp);
        cio_write(cio, len, 2);                 /* Lqcd */
        cio_seek(cio, lenp + len);
+
+       if(j2k->cstr_info)
+         j2k_add_mhmarker(j2k->cstr_info, J2K_MS_QCD, lenp, len);
 }
 
 static void j2k_read_qcd(opj_j2k_t *j2k) {
@@ -930,6 +2025,41 @@ static void j2k_read_qcd(opj_j2k_t *j2k) {
        }
 }
 
+/**
+ * Reads a QCD marker (Quantization defaults)
+ * @param      p_header_data   the data contained in the QCD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the QCD marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_qcd_v2 (
+                                   opj_j2k_v2_t *p_j2k,
+                                       OPJ_BYTE * p_header_data,
+                                       OPJ_UINT32 p_header_size,
+                                       struct opj_event_mgr * p_manager
+                                       )
+{
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       if (! j2k_read_SQcd_SQcc(p_j2k,0,p_header_data,&p_header_size,p_manager)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading QCD marker\n");
+               return OPJ_FALSE;
+       }
+
+       if (p_header_size != 0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading QCD marker\n");
+               return OPJ_FALSE;
+       }
+
+       /* Apply the quantization parameters to other components of the current tile or the m_default_tcp */
+       j2k_copy_tile_quantization_parameters(p_j2k);
+
+       return OPJ_TRUE;
+}
+
 static void j2k_write_qcc(opj_j2k_t *j2k, int compno) {
        int lenp, len;
 
@@ -950,7 +2080,7 @@ static void j2k_read_qcc(opj_j2k_t *j2k) {
        int len, compno;
        int numcomp = j2k->image->numcomps;
        opj_cio_t *cio = j2k->cio;
-       
+
        len = cio_read(cio, 2); /* Lqcc */
        compno = cio_read(cio, numcomp <= 256 ? 1 : 2); /* Cqcc */
 
@@ -983,6 +2113,86 @@ static void j2k_read_qcc(opj_j2k_t *j2k) {
        j2k_read_qcx(j2k, compno, len - 2 - (numcomp <= 256 ? 1 : 2));
 }
 
+/**
+ * Reads a QCC marker (Quantization component)
+ * @param      p_header_data   the data contained in the QCC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the QCC marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_qcc_v2(      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_num_comp,l_comp_no;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_num_comp = p_j2k->m_image->numcomps;
+
+       if (l_num_comp <= 256) {
+               if (p_header_size < 1) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading QCC marker\n");
+                       return OPJ_FALSE;
+               }
+               opj_read_bytes(p_header_data,&l_comp_no,1);
+               ++p_header_data;
+               --p_header_size;
+       }
+       else {
+               if (p_header_size < 2) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading QCC marker\n");
+                       return OPJ_FALSE;
+               }
+               opj_read_bytes(p_header_data,&l_comp_no,2);
+               p_header_data+=2;
+               p_header_size-=2;
+       }
+
+#ifdef USE_JPWL
+       if (p_j2k->m_cp.correct) {
+
+               static OPJ_UINT32 backup_compno = 0;
+
+               /* compno is negative or larger than the number of components!!! */
+               if ((l_comp_no < 0) || (l_comp_no >= l_num_comp)) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR,
+                               "JPWL: bad component number in QCC (%d out of a maximum of %d)\n",
+                               l_comp_no, l_num_comp);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+                       /* we try to correct */
+                       l_comp_no = backup_compno % l_num_comp;
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "- trying to adjust this\n"
+                               "- setting component number to %d\n",
+                               l_comp_no);
+               }
+
+               /* keep your private count of tiles */
+               backup_compno++;
+       };
+#endif /* USE_JPWL */
+
+       if (! j2k_read_SQcd_SQcc(p_j2k,l_comp_no,p_header_data,&p_header_size,p_manager)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading QCC marker\n");
+               return OPJ_FALSE;
+       }
+
+       if (p_header_size != 0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading QCC marker\n");
+               return OPJ_FALSE;
+       }
+
+       return OPJ_TRUE;
+}
+
+
 static void j2k_write_poc(opj_j2k_t *j2k) {
        int len, numpchgs, i;
 
@@ -1018,7 +2228,6 @@ static void j2k_read_poc(opj_j2k_t *j2k) {
        
        opj_cp_t *cp = j2k->cp;
        opj_tcp_t *tcp = j2k->state == J2K_STATE_TPH ? &cp->tcps[j2k->curtileno] : j2k->default_tcp;
-       opj_tccp_t *tccp = &tcp->tccps[0];
        opj_cio_t *cio = j2k->cio;
        
        old_poc = tcp->POC ? tcp->numpocs + 1 : 0;
@@ -1041,6 +2250,87 @@ static void j2k_read_poc(opj_j2k_t *j2k) {
        tcp->numpocs = numpchgs + old_poc - 1;
 }
 
+/**
+ * Reads a POC marker (Progression Order Change)
+ *
+ * @param      p_header_data   the data contained in the POC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the POC marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_poc_v2 (
+                                               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, l_nb_comp, l_tmp;
+       opj_image_t * l_image = 00;
+       OPJ_UINT32 l_old_poc_nb, l_current_poc_nb, l_current_poc_remaining;
+       OPJ_UINT32 l_chunk_size, l_comp_room;
+
+       opj_cp_v2_t *l_cp = 00;
+       opj_tcp_v2_t *l_tcp = 00;
+       opj_poc_t *l_current_poc = 00;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_image = p_j2k->m_image;
+       l_nb_comp = l_image->numcomps;
+       if (l_nb_comp <= 256) {
+               l_comp_room = 1;
+       }
+       else {
+               l_comp_room = 2;
+       }
+       l_chunk_size = 5 + 2 * l_comp_room;
+       l_current_poc_nb = p_header_size / l_chunk_size;
+       l_current_poc_remaining = p_header_size % l_chunk_size;
+
+       if ((l_current_poc_nb <= 0) || (l_current_poc_remaining != 0)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading POC marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_cp = &(p_j2k->m_cp);
+       l_tcp = (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH) ?
+                               &l_cp->tcps[p_j2k->m_current_tile_number] :
+                               p_j2k->m_specific_param.m_decoder.m_default_tcp;
+       l_old_poc_nb = l_tcp->POC ? l_tcp->numpocs + 1 : 0;
+       l_current_poc_nb += l_old_poc_nb;
+
+       assert(l_current_poc_nb < 32);
+
+       /* now poc is in use.*/
+       l_tcp->POC = 1;
+
+       l_current_poc = &l_tcp->pocs[l_old_poc_nb];
+       for     (i = l_old_poc_nb; i < l_current_poc_nb; ++i) {
+               opj_read_bytes(p_header_data,&(l_current_poc->resno0),1);                               /* RSpoc_i */
+               ++p_header_data;
+               opj_read_bytes(p_header_data,&(l_current_poc->compno0),l_comp_room);    /* CSpoc_i */
+               p_header_data+=l_comp_room;
+               opj_read_bytes(p_header_data,&(l_current_poc->layno1),2);                               /* LYEpoc_i */
+               p_header_data+=2;
+               opj_read_bytes(p_header_data,&(l_current_poc->resno1),1);                               /* REpoc_i */
+               ++p_header_data;
+               opj_read_bytes(p_header_data,&(l_current_poc->compno1),l_comp_room);    /* CEpoc_i */
+               p_header_data+=l_comp_room;
+               opj_read_bytes(p_header_data,&l_tmp,1);                                                                 /* Ppoc_i */
+               ++p_header_data;
+               l_current_poc->prg = (OPJ_PROG_ORDER) l_tmp;
+               /* make sure comp is in acceptable bounds */
+               l_current_poc->compno1 = uint_min(l_current_poc->compno1, l_nb_comp);
+               ++l_current_poc;
+       }
+
+       l_tcp->numpocs = l_current_poc_nb - 1;
+       return OPJ_TRUE;
+}
+
 static void j2k_read_crg(opj_j2k_t *j2k) {
        int len, i, Xcrg_i, Ycrg_i;
        
@@ -1054,6 +2344,47 @@ static void j2k_read_crg(opj_j2k_t *j2k) {
        }
 }
 
+/**
+ * Reads a CRG marker (Component registration)
+ *
+ * @param      p_header_data   the data contained in the TLM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the TLM marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_crg_v2 (
+                                               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;
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_nb_comp = p_j2k->m_image->numcomps;
+
+       if (p_header_size != l_nb_comp *4) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading CRG marker\n");
+               return OPJ_FALSE;
+       }
+       /* Do not care of this at the moment since only local variables are set here */
+       /*
+       for
+               (i = 0; i < l_nb_comp; ++i)
+       {
+               opj_read_bytes(p_header_data,&l_Xcrg_i,2);                              // Xcrg_i
+               p_header_data+=2;
+               opj_read_bytes(p_header_data,&l_Ycrg_i,2);                              // Xcrg_i
+               p_header_data+=2;
+       }
+       */
+       return OPJ_TRUE;
+}
+
 static void j2k_read_tlm(opj_j2k_t *j2k) {
        int len, Ztlm, Stlm, ST, SP, tile_tlm, i;
        long int Ttlm_i, Ptlm_i;
@@ -1072,6 +2403,64 @@ static void j2k_read_tlm(opj_j2k_t *j2k) {
        }
 }
 
+/**
+ * Reads a TLM marker (Tile Length Marker)
+ *
+ * @param      p_header_data   the data contained in the TLM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the TLM marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_tlm_v2 (
+                                               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_Ztlm, l_Stlm, l_ST, l_SP, l_tot_num_tp, l_tot_num_tp_remaining, l_quotient, l_Ptlm_size;
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       if (p_header_size < 2) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading TLM marker\n");
+               return OPJ_FALSE;
+       }
+       p_header_size -= 2;
+
+       opj_read_bytes(p_header_data,&l_Ztlm,1);                                /* Ztlm */
+       ++p_header_data;
+       opj_read_bytes(p_header_data,&l_Stlm,1);                                /* Stlm */
+       ++p_header_data;
+
+       l_ST = ((l_Stlm >> 4) & 0x3);
+       l_SP = (l_Stlm >> 6) & 0x1;
+
+       l_Ptlm_size = (l_SP + 1) * 2;
+       l_quotient = l_Ptlm_size + l_ST;
+
+       l_tot_num_tp = p_header_size / l_quotient;
+       l_tot_num_tp_remaining = p_header_size % l_quotient;
+
+       if (l_tot_num_tp_remaining != 0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading TLM marker\n");
+               return OPJ_FALSE;
+       }
+       /* FIXME Do not care of this at the moment since only local variables are set here */
+       /*
+       for
+               (i = 0; i < l_tot_num_tp; ++i)
+       {
+               opj_read_bytes(p_header_data,&l_Ttlm_i,l_ST);                           // Ttlm_i
+               p_header_data += l_ST;
+               opj_read_bytes(p_header_data,&l_Ptlm_i,l_Ptlm_size);            // Ptlm_i
+               p_header_data += l_Ptlm_size;
+       }*/
+       return OPJ_TRUE;
+}
+
 static void j2k_read_plm(opj_j2k_t *j2k) {
        int len, i, Zplm, Nplm, add, packet_len = 0;
        
@@ -1097,9 +2486,80 @@ static void j2k_read_plm(opj_j2k_t *j2k) {
        }
 }
 
-static void j2k_read_plt(opj_j2k_t *j2k) {
-       int len, i, Zplt, packet_len = 0, add;
-       
+/**
+ * Reads a PLM marker (Packet length, main header marker)
+ *
+ * @param      p_header_data   the data contained in the TLM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the TLM marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_plm_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       )
+{
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       if (p_header_size < 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PLM marker\n");
+               return OPJ_FALSE;
+       }
+       /* Do not care of this at the moment since only local variables are set here */
+       /*
+       opj_read_bytes(p_header_data,&l_Zplm,1);                                        // Zplm
+       ++p_header_data;
+       --p_header_size;
+
+       while
+               (p_header_size > 0)
+       {
+               opj_read_bytes(p_header_data,&l_Nplm,1);                                // Nplm
+               ++p_header_data;
+               p_header_size -= (1+l_Nplm);
+               if
+                       (p_header_size < 0)
+               {
+                       opj_event_msg(p_manager, EVT_ERROR, "Error reading PLM marker\n");
+                       return false;
+               }
+               for
+                       (i = 0; i < l_Nplm; ++i)
+               {
+                       opj_read_bytes(p_header_data,&l_tmp,1);                         // Iplm_ij
+                       ++p_header_data;
+                       // take only the last seven bytes
+                       l_packet_len |= (l_tmp & 0x7f);
+                       if
+                               (l_tmp & 0x80)
+                       {
+                               l_packet_len <<= 7;
+                       }
+                       else
+                       {
+                // store packet length and proceed to next packet
+                               l_packet_len = 0;
+                       }
+               }
+               if
+                       (l_packet_len != 0)
+               {
+                       opj_event_msg(p_manager, EVT_ERROR, "Error reading PLM marker\n");
+                       return false;
+               }
+       }
+       */
+       return OPJ_TRUE;
+}
+
+static void j2k_read_plt(opj_j2k_t *j2k) {
+       int len, i, Zplt, packet_len = 0, add;
+       
        opj_cio_t *cio = j2k->cio;
        
        len = cio_read(cio, 2);         /* Lplt */
@@ -1114,6 +2574,59 @@ static void j2k_read_plt(opj_j2k_t *j2k) {
        }
 }
 
+/**
+ * Reads a PLT marker (Packet length, tile-part header)
+ *
+ * @param      p_header_data   the data contained in the PLT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the PLT marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_plt_v2 (
+                                               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_Zplt, l_tmp, l_packet_len = 0, i;
+
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       if (p_header_size < 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PLM marker\n");
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(p_header_data,&l_Zplt,1);                /* Zplt */
+       ++p_header_data;
+       --p_header_size;
+
+       for (i = 0; i < p_header_size; ++i) {
+               opj_read_bytes(p_header_data,&l_tmp,1);         /* Iplt_ij */
+               ++p_header_data;
+               // take only the last seven bytes
+               l_packet_len |= (l_tmp & 0x7f);
+               if (l_tmp & 0x80) {
+                       l_packet_len <<= 7;
+               }
+               else {
+            // store packet length and proceed to next packet
+                       l_packet_len = 0;
+               }
+       }
+
+       if (l_packet_len != 0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PLM marker\n");
+               return OPJ_FALSE;
+       }
+
+       return OPJ_TRUE;
+}
+
 static void j2k_read_ppm(opj_j2k_t *j2k) {
        int len, Z_ppm, i, j;
        int N_ppm;
@@ -1149,7 +2662,7 @@ static void j2k_read_ppm(opj_j2k_t *j2k) {
                                                "JPWL: failed memory allocation during PPM marker parsing (pos. %x)\n",
                                                cio_tell(cio));
                                        if (!JPWL_ASSUME || JPWL_ASSUME) {
-                                               free(cp->ppm_data);
+                                               opj_free(cp->ppm_data);
                                                opj_event_msg(j2k->cinfo, EVT_ERROR, "JPWL: giving up\n");
                                                return;
                                        }
@@ -1171,6 +2684,298 @@ static void j2k_read_ppm(opj_j2k_t *j2k) {
                cp->ppm_store = j;
        }
 }
+/**
+ * Reads a PPM marker (Packed packet headers, main header)
+ *
+ * @param      p_header_data   the data contained in the POC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the POC marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_ppm_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       )
+{
+
+       opj_cp_v2_t *l_cp = 00;
+       OPJ_UINT32 l_remaining_data, l_Z_ppm, l_N_ppm;
+
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       if (p_header_size < 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PPM marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_cp = &(p_j2k->m_cp);
+       l_cp->ppm = 1;
+
+       opj_read_bytes(p_header_data,&l_Z_ppm,1);               /* Z_ppm */
+       ++p_header_data;
+       --p_header_size;
+
+       // First PPM marker
+       if (l_Z_ppm == 0) {
+               if (p_header_size < 4) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PPM marker\n");
+                       return OPJ_FALSE;
+               }
+
+               opj_read_bytes(p_header_data,&l_N_ppm,4);               /* N_ppm */
+               p_header_data+=4;
+               p_header_size-=4;
+
+               /* First PPM marker: Initialization */
+               l_cp->ppm_len = l_N_ppm;
+               l_cp->ppm_data_size = 0;
+
+               l_cp->ppm_buffer = (OPJ_BYTE *) opj_malloc(l_cp->ppm_len);
+               if (l_cp->ppm_buffer == 00) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory reading ppm marker\n");
+                       return OPJ_FALSE;
+               }
+               memset(l_cp->ppm_buffer,0,l_cp->ppm_len);
+
+               l_cp->ppm_data = l_cp->ppm_buffer;
+       }
+
+       while (1) {
+               if (l_cp->ppm_data_size == l_cp->ppm_len) {
+                       if (p_header_size >= 4) {
+                               // read a N_ppm
+                               opj_read_bytes(p_header_data,&l_N_ppm,4);               /* N_ppm */
+                               p_header_data+=4;
+                               p_header_size-=4;
+                               l_cp->ppm_len += l_N_ppm ;
+
+                               l_cp->ppm_buffer = (OPJ_BYTE *) opj_realloc(l_cp->ppm_buffer, l_cp->ppm_len);
+                               if (l_cp->ppm_buffer == 00) {
+                                       opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory reading ppm marker\n");
+                                       return OPJ_FALSE;
+                               }
+                               memset(l_cp->ppm_buffer+l_cp->ppm_data_size,0,l_N_ppm);
+
+                               l_cp->ppm_data = l_cp->ppm_buffer;
+                       }
+                       else {
+                               return OPJ_FALSE;
+                       }
+               }
+
+               l_remaining_data = l_cp->ppm_len - l_cp->ppm_data_size;
+
+               if (l_remaining_data <= p_header_size) {
+                       /* we must store less information than available in the packet */
+                       memcpy(l_cp->ppm_buffer + l_cp->ppm_data_size , p_header_data , l_remaining_data);
+                       l_cp->ppm_data_size = l_cp->ppm_len;
+                       p_header_size -= l_remaining_data;
+                       p_header_data += l_remaining_data;
+               }
+               else {
+                       memcpy(l_cp->ppm_buffer + l_cp->ppm_data_size , p_header_data , p_header_size);
+                       l_cp->ppm_data_size += p_header_size;
+                       p_header_data += p_header_size;
+                       p_header_size = 0;
+                       break;
+               }
+       }
+
+       return OPJ_TRUE;
+}
+
+
+
+/**
+ * Reads a PPM marker (Packed packet headers, main header)
+ *
+ * @param      p_header_data   the data contained in the POC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the POC marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_ppm_v3 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       )
+{
+       opj_cp_v2_t *l_cp = 00;
+       OPJ_UINT32 l_remaining_data, l_Z_ppm, l_N_ppm;
+
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       // Minimum size of PPM marker is equal to the size of Zppm element
+       if (p_header_size < 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PPM marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_cp = &(p_j2k->m_cp);
+       l_cp->ppm = 1;
+
+       opj_read_bytes(p_header_data,&l_Z_ppm,1);               /* Z_ppm */
+       ++p_header_data;
+       --p_header_size;
+
+       // First PPM marker
+       if (l_Z_ppm == 0) {
+               // We need now at least the Nppm^0 element
+               if (p_header_size < 4) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PPM marker\n");
+                       return OPJ_FALSE;
+               }
+
+               opj_read_bytes(p_header_data,&l_N_ppm,4);               /* First N_ppm */
+               p_header_data+=4;
+               p_header_size-=4;
+
+               /* First PPM marker: Initialization */
+               l_cp->ppm_len = l_N_ppm;
+               l_cp->ppm_data_read = 0;
+
+               l_cp->ppm_data = (OPJ_BYTE *) opj_malloc(l_cp->ppm_len);
+               if (l_cp->ppm_data == 00) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory reading ppm marker\n");
+                       return OPJ_FALSE;
+               }
+               memset(l_cp->ppm_data,0,l_cp->ppm_len);
+
+               l_cp->ppm_data_current = l_cp->ppm_data;
+
+               //l_cp->ppm_data = l_cp->ppm_buffer;
+       }
+       else {
+               if (p_header_size < 4) {
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "Empty PPM marker\n");
+                       return OPJ_TRUE;
+               }
+               else {
+                       // Uncompleted Ippm series in the previous PPM marker?
+                       if (l_cp->ppm_data_read < l_cp->ppm_len) {
+                               // Get the place where add the remaining Ippm series
+                               l_cp->ppm_data_current = &(l_cp->ppm_data[l_cp->ppm_data_read]);
+                               l_N_ppm = l_cp->ppm_len - l_cp->ppm_data_read;
+                       }
+                       else {
+                               opj_read_bytes(p_header_data,&l_N_ppm,4);               /* First N_ppm */
+                               p_header_data+=4;
+                               p_header_size-=4;
+
+                               // Increase the size of ppm_data to add the new Ippm series
+                               l_cp->ppm_data = (OPJ_BYTE *) opj_realloc(l_cp->ppm_data, l_cp->ppm_len + l_N_ppm);
+
+                               // Keep the position of the place where concatenate the new series
+                               l_cp->ppm_data_current = &(l_cp->ppm_data[l_cp->ppm_len]);
+                               l_cp->ppm_len += l_N_ppm;
+                       }
+               }
+       }
+
+       l_remaining_data = p_header_size;
+
+       while (l_remaining_data >= l_N_ppm) {
+               // read a complete Ippm series
+               memcpy(l_cp->ppm_data_current, p_header_data, l_N_ppm);
+               p_header_size -= l_N_ppm;
+               p_header_data += l_N_ppm;
+
+               l_cp->ppm_data_read += l_N_ppm; // Increase the number of data read
+
+               if (p_header_size)
+               {
+                       opj_read_bytes(p_header_data,&l_N_ppm,4);               /* N_ppm^i */
+                       p_header_data+=4;
+                       p_header_size-=4;
+               }
+               else {
+                       l_remaining_data = p_header_size;
+                       break;
+               }
+
+               l_remaining_data = p_header_size;
+
+               // Next Ippm series is a complete series ?
+               if (l_remaining_data > l_N_ppm) {
+                       // Increase the size of ppm_data to add the new Ippm series
+                       l_cp->ppm_data = (OPJ_BYTE *) opj_realloc(l_cp->ppm_data, l_cp->ppm_len + l_N_ppm);
+
+                       // Keep the position of the place where concatenate the new series
+                       l_cp->ppm_data_current = &(l_cp->ppm_data[l_cp->ppm_len]);
+                       l_cp->ppm_len += l_N_ppm;
+               }
+
+       }
+
+       // Need to read an incomplete Ippm series
+       if (l_remaining_data) {
+               l_cp->ppm_data = (OPJ_BYTE *) opj_realloc(l_cp->ppm_data, l_cp->ppm_len + l_N_ppm);
+
+               // Keep the position of the place where concatenate the new series
+               l_cp->ppm_data_current = &(l_cp->ppm_data[l_cp->ppm_len]);
+               l_cp->ppm_len += l_N_ppm;
+
+               // Read incomplete Ippm series
+               memcpy(l_cp->ppm_data_current, p_header_data, l_remaining_data);
+               p_header_size -= l_remaining_data;
+               p_header_data += l_remaining_data;
+
+               l_cp->ppm_data_read += l_remaining_data; // Increase the number of data read
+       }
+
+#ifdef CLEAN_MSD
+
+               if (l_cp->ppm_data_size == l_cp->ppm_len) {
+                       if (p_header_size >= 4) {
+                               // read a N_ppm
+                               opj_read_bytes(p_header_data,&l_N_ppm,4);               /* N_ppm */
+                               p_header_data+=4;
+                               p_header_size-=4;
+                               l_cp->ppm_len += l_N_ppm ;
+
+                               l_cp->ppm_buffer = (OPJ_BYTE *) opj_realloc(l_cp->ppm_buffer, l_cp->ppm_len);
+                               if (l_cp->ppm_buffer == 00) {
+                                       opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory reading ppm marker\n");
+                                       return OPJ_FALSE;
+                               }
+                               memset(l_cp->ppm_buffer+l_cp->ppm_data_size,0,l_N_ppm);
+
+                               l_cp->ppm_data = l_cp->ppm_buffer;
+                       }
+                       else {
+                               return OPJ_FALSE;
+                       }
+               }
+
+               l_remaining_data = l_cp->ppm_len - l_cp->ppm_data_size;
+
+               if (l_remaining_data <= p_header_size) {
+                       /* we must store less information than available in the packet */
+                       memcpy(l_cp->ppm_buffer + l_cp->ppm_data_size , p_header_data , l_remaining_data);
+                       l_cp->ppm_data_size = l_cp->ppm_len;
+                       p_header_size -= l_remaining_data;
+                       p_header_data += l_remaining_data;
+               }
+               else {
+                       memcpy(l_cp->ppm_buffer + l_cp->ppm_data_size , p_header_data , p_header_size);
+                       l_cp->ppm_data_size += p_header_size;
+                       p_header_data += p_header_size;
+                       p_header_size = 0;
+                       break;
+               }
+       }
+#endif
+       return OPJ_TRUE;
+}
 
 static void j2k_read_ppt(opj_j2k_t *j2k) {
        int len, Z_ppt, i, j = 0;
@@ -1200,6 +3005,83 @@ static void j2k_read_ppt(opj_j2k_t *j2k) {
        tcp->ppt_store = j;
 }
 
+/**
+ * Reads a PPT marker (Packed packet headers, tile-part header)
+ *
+ * @param      p_header_data   the data contained in the PPT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the PPT marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_ppt_v2 (     opj_j2k_v2_t *p_j2k,
+                                                       OPJ_BYTE * p_header_data,
+                                                       OPJ_UINT32 p_header_size,
+                                                       struct opj_event_mgr * p_manager )
+{
+       opj_cp_v2_t *l_cp = 00;
+       opj_tcp_v2_t *l_tcp = 00;
+       OPJ_UINT32 l_Z_ppt;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       /* We need to have the Z_ppt element at minimum */
+       if (p_header_size < 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PPT marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_cp = &(p_j2k->m_cp);
+       if (l_cp->ppm){
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading PPT marker: packet header have been previously found in the main header (PPM marker).\n");
+               return OPJ_FALSE;
+       }
+
+       l_tcp = &(l_cp->tcps[p_j2k->m_current_tile_number]);
+       l_tcp->ppt = 1;
+
+       opj_read_bytes(p_header_data,&l_Z_ppt,1);               /* Z_ppt */
+       ++p_header_data;
+       --p_header_size;
+
+       /* Allocate buffer to read the packet header */
+       if (l_Z_ppt == 0) {
+               /* First PPT marker */
+               l_tcp->ppt_data_size = 0;
+               l_tcp->ppt_len = p_header_size;
+
+               l_tcp->ppt_buffer = (OPJ_BYTE *) opj_calloc(l_tcp->ppt_len, sizeof(OPJ_BYTE) );
+               if (l_tcp->ppt_buffer == 00) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory reading PPT marker\n");
+                       return OPJ_FALSE;
+               }
+               l_tcp->ppt_data = l_tcp->ppt_buffer;
+
+               /* memset(l_tcp->ppt_buffer,0,l_tcp->ppt_len); */
+       }
+       else {
+               l_tcp->ppt_len += p_header_size;
+
+               l_tcp->ppt_buffer = (OPJ_BYTE *) opj_realloc(l_tcp->ppt_buffer,l_tcp->ppt_len);
+               if (l_tcp->ppt_buffer == 00) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Not enough memory reading PPT marker\n");
+                       return OPJ_FALSE;
+               }
+               l_tcp->ppt_data = l_tcp->ppt_buffer;
+
+               memset(l_tcp->ppt_buffer+l_tcp->ppt_data_size,0,p_header_size);
+       }
+
+       /* Read packet header from buffer */
+       memcpy(l_tcp->ppt_buffer+l_tcp->ppt_data_size,p_header_data,p_header_size);
+
+       l_tcp->ppt_data_size += p_header_size;
+
+       return OPJ_TRUE;
+}
+
 static void j2k_write_tlm(opj_j2k_t *j2k){
        int lenp;
        opj_cio_t *cio = j2k->cio;
@@ -1229,6 +3111,17 @@ static void j2k_write_sot(opj_j2k_t *j2k) {
        cio_seek(cio, lenp);
        cio_write(cio, len, 2);                         /* Lsot */
        cio_seek(cio, lenp + len);
+
+       /* UniPG>> */
+#ifdef USE_JPWL
+       /* update markers struct */
+       j2k_add_marker(j2k->cstr_info, J2K_MS_SOT, j2k->sot_start, len + 2);
+#endif /* USE_JPWL */
+       /* <<UniPG */
+
+       if( j2k->cstr_info && j2k->cur_tp_num==0){
+         j2k_add_tlmarker( j2k->curtileno, j2k->cstr_info, J2K_MS_SOT, lenp, len);
+       }
 }
 
 static void j2k_read_sot(opj_j2k_t *j2k) {
@@ -1238,7 +3131,7 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
 
        opj_cp_t *cp = j2k->cp;
        opj_cio_t *cio = j2k->cio;
-       
+
        len = cio_read(cio, 2);
        tileno = cio_read(cio, 2);
 
@@ -1267,7 +3160,6 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
                backup_tileno++;
        };
 #endif /* USE_JPWL */
-
        
        if (cp->tileno_size == 0) {
                cp->tileno[cp->tileno_size] = tileno;
@@ -1313,14 +3205,40 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
        
        partno = cio_read(cio, 1);
        numparts = cio_read(cio, 1);
+  
+  if (partno >= numparts) {
+    opj_event_msg(j2k->cinfo, EVT_WARNING, "SOT marker inconsistency in tile %d: tile-part index greater (%d) than number of tile-parts (%d)\n", tileno, partno, numparts);
+    numparts = partno+1;
+  }
        
        j2k->curtileno = tileno;
+       j2k->cur_tp_num = partno;
        j2k->eot = cio_getbp(cio) - 12 + totlen;
        j2k->state = J2K_STATE_TPH;
        tcp = &cp->tcps[j2k->curtileno];
+
+       /* Index */
+       if (j2k->cstr_info) {
+               if (tcp->first) {
+                       if (tileno == 0) 
+                               j2k->cstr_info->main_head_end = cio_tell(cio) - 13;
+                       j2k->cstr_info->tile[tileno].tileno = tileno;
+                       j2k->cstr_info->tile[tileno].start_pos = cio_tell(cio) - 12;
+                       j2k->cstr_info->tile[tileno].end_pos = j2k->cstr_info->tile[tileno].start_pos + totlen - 1;                             
+    } else {
+                       j2k->cstr_info->tile[tileno].end_pos += totlen;
+               }
+    j2k->cstr_info->tile[tileno].num_tps = numparts;
+    if (numparts)
+      j2k->cstr_info->tile[tileno].tp = (opj_tp_info_t *) opj_realloc(j2k->cstr_info->tile[tileno].tp, numparts * sizeof(opj_tp_info_t));
+    else
+      j2k->cstr_info->tile[tileno].tp = (opj_tp_info_t *) opj_realloc(j2k->cstr_info->tile[tileno].tp, 10 * sizeof(opj_tp_info_t)); // Fixme (10)
+               j2k->cstr_info->tile[tileno].tp[partno].tp_start_pos = cio_tell(cio) - 12;
+               j2k->cstr_info->tile[tileno].tp[partno].tp_end_pos = 
+                       j2k->cstr_info->tile[tileno].tp[partno].tp_start_pos + totlen - 1;
+       }
        
-       if (tcp->first == 1) {
-               
+       if (tcp->first == 1) {          
                /* Initialization PPT */
                opj_tccp_t *tmp = tcp->tccps;
                memcpy(tcp, j2k->default_tcp, sizeof(opj_tcp_t));
@@ -1336,11 +3254,192 @@ static void j2k_read_sot(opj_j2k_t *j2k) {
        }
 }
 
+/**
+ * Reads a PPT marker (Packed packet headers, tile-part header)
+ *
+ * @param      p_header_data   the data contained in the PPT box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the PPT marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_sot_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               OPJ_BYTE * p_header_data,
+                                               OPJ_UINT32 p_header_size,
+                                               struct opj_event_mgr * p_manager
+                                       )
+{
+
+       opj_cp_v2_t *l_cp = 00;
+       opj_tcp_v2_t *l_tcp = 00;
+       OPJ_UINT32 l_tot_len, l_num_parts = 0;
+       OPJ_UINT32 l_current_part;
+       OPJ_UINT32 l_tile_x,l_tile_y;
+
+       /* preconditions */
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       /* Size of this marker is fixed = 12 (we have already read marker and its size)*/
+       if (p_header_size != 8) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading SOT marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_cp = &(p_j2k->m_cp);
+       opj_read_bytes(p_header_data,&(p_j2k->m_current_tile_number),2);                /* Isot */
+       p_header_data+=2;
+
+       l_tcp = &l_cp->tcps[p_j2k->m_current_tile_number];
+       l_tile_x = p_j2k->m_current_tile_number % l_cp->tw;
+       l_tile_y = p_j2k->m_current_tile_number / l_cp->tw;
+
+#ifdef USE_JPWL
+       if (l_cp->correct) {
+
+               int tileno = p_j2k->m_current_tile_number;
+               static int backup_tileno = 0;
+
+               /* tileno is negative or larger than the number of tiles!!! */
+               if ((tileno < 0) || (tileno > (l_cp->tw * l_cp->th))) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR,
+                               "JPWL: bad tile number (%d out of a maximum of %d)\n",
+                               tileno, (l_cp->tw * l_cp->th));
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+                       /* we try to correct */
+                       tileno = backup_tileno;
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "- trying to adjust this\n"
+                               "- setting tile number to %d\n",
+                               tileno);
+               }
+
+               /* keep your private count of tiles */
+               backup_tileno++;
+       };
+#endif /* USE_JPWL */
+
+       /* look for the tile in the list of already processed tile (in parts). */
+       /* Optimization possible here with a more complex data structure and with the removing of tiles */
+       /* since the time taken by this function can only grow at the time */
+
+       opj_read_bytes(p_header_data,&l_tot_len,4);             /* Psot */
+       p_header_data+=4;
+
+#ifdef USE_JPWL
+       if (l_cp->correct) {
+
+               /* totlen is negative or larger than the bytes left!!! */
+               if ((l_tot_len < 0) || (l_tot_len > p_header_size ) ) { /* FIXME it seems correct; for info in V1 -> (p_stream_numbytesleft(p_stream) + 8))) { */
+                       opj_event_msg_v2(p_manager, EVT_ERROR,
+                               "JPWL: bad tile byte size (%d bytes against %d bytes left)\n",
+                               l_tot_len, p_header_size ); /* FIXME it seems correct; for info in V1 -> p_stream_numbytesleft(p_stream) + 8); */
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+                       /* we try to correct */
+                       l_tot_len = 0;
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "- trying to adjust this\n"
+                               "- setting Psot to %d => assuming it is the last tile\n",
+                               l_tot_len);
+               }
+       };
+#endif /* USE_JPWL */
+
+       /* 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_read_bytes(p_header_data,&l_current_part ,1);       /* TPsot */
+       ++p_header_data;
+
+       opj_read_bytes(p_header_data,&l_num_parts ,1);          /* TNsot */
+       ++p_header_data;
+
+       if (l_num_parts != 0) { /* Number of tile-part header is provided by this tile-part header */
+               /* Useful to manage the case of textGBR.jp2 file because two values of TNSot are allowed: the correct numbers of
+                * tile-parts for that tile and zero (A.4.2 of 15444-1 : 2002). */
+               if (l_tcp->m_nb_tile_parts) {
+                       if (l_current_part >= l_tcp->m_nb_tile_parts){
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "In SOT marker, TPSot (%d) is not valid regards to the current "
+                                               "number of tile-part (%d), giving up\n", l_current_part, l_tcp->m_nb_tile_parts );
+                               p_j2k->m_specific_param.m_decoder.m_last_tile_part = 1;
+                               return OPJ_FALSE;
+                       }
+               }
+               l_tcp->m_nb_tile_parts = l_num_parts;
+       }
+
+       /* If know the number of tile part header we will check if we didn't read the last*/
+       if (l_tcp->m_nb_tile_parts) {
+               if (l_tcp->m_nb_tile_parts == (l_current_part + 1)) {
+                       p_j2k->m_specific_param.m_decoder.m_can_decode = 1; /* Process the last tile-part header*/
+               }
+       }
+
+       if (!p_j2k->m_specific_param.m_decoder.m_last_tile_part){
+               /* Keep the size of data to skip after this marker */
+               p_j2k->m_specific_param.m_decoder.m_sot_length = l_tot_len - 12; /* SOT_marker_size = 12 */
+       }
+       else {
+               /* FIXME: need to be computed from the number of bytes remaining in the codestream */
+               p_j2k->m_specific_param.m_decoder.m_sot_length = 0;
+       }
+
+       p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_TPH;
+
+       /* Check if the current tile is outside the area we want decode (in tile index)*/
+       p_j2k->m_specific_param.m_decoder.m_skip_data =
+                       (l_tile_x < p_j2k->m_specific_param.m_decoder.m_start_tile_x)
+               ||      (l_tile_x >= p_j2k->m_specific_param.m_decoder.m_end_tile_x)
+               ||  (l_tile_y < p_j2k->m_specific_param.m_decoder.m_start_tile_y)
+               ||      (l_tile_y >= p_j2k->m_specific_param.m_decoder.m_end_tile_y);
+
+       /* Index */
+
+       /* FIXME move this onto a separate method to call before reading any SOT, remove part about main_end header, use a index struct inside p_j2k */
+       /* if (p_j2k->cstr_info) {
+               if (l_tcp->first) {
+                       if (tileno == 0) {
+                               p_j2k->cstr_info->main_head_end = p_stream_tell(p_stream) - 13;
+                       }
+
+                       p_j2k->cstr_info->tile[tileno].tileno = tileno;
+                       p_j2k->cstr_info->tile[tileno].start_pos = p_stream_tell(p_stream) - 12;
+                       p_j2k->cstr_info->tile[tileno].end_pos = p_j2k->cstr_info->tile[tileno].start_pos + totlen - 1;
+                       p_j2k->cstr_info->tile[tileno].num_tps = numparts;
+
+                       if (numparts) {
+                               p_j2k->cstr_info->tile[tileno].tp = (opj_tp_info_t *) opj_malloc(numparts * sizeof(opj_tp_info_t));
+                       }
+                       else {
+                               p_j2k->cstr_info->tile[tileno].tp = (opj_tp_info_t *) opj_malloc(10 * sizeof(opj_tp_info_t)); // Fixme (10)
+                       }
+               }
+               else {
+                       p_j2k->cstr_info->tile[tileno].end_pos += totlen;
+               }
+
+               p_j2k->cstr_info->tile[tileno].tp[partno].tp_start_pos = p_stream_tell(p_stream) - 12;
+               p_j2k->cstr_info->tile[tileno].tp[partno].tp_end_pos =
+               p_j2k->cstr_info->tile[tileno].tp[partno].tp_start_pos + totlen - 1;
+       }*/
+       return OPJ_TRUE;
+}
+
 static void j2k_write_sod(opj_j2k_t *j2k, void *tile_coder) {
        int l, layno;
        int totlen;
        opj_tcp_t *tcp = NULL;
-       opj_image_info_t *image_info = NULL;
+       opj_codestream_info_t *cstr_info = NULL;
        
        opj_tcd_t *tcd = (opj_tcd_t*)tile_coder;        /* cast is needed because of conflicts in header inclusions */
        opj_cp_t *cp = j2k->cp;
@@ -1348,34 +3447,52 @@ static void j2k_write_sod(opj_j2k_t *j2k, void *tile_coder) {
 
        tcd->tp_num = j2k->tp_num ;
        tcd->cur_tp_num = j2k->cur_tp_num;
-       tcd->cur_totnum_tp = j2k->cur_totnum_tp[j2k->curtileno];
        
        cio_write(cio, J2K_MS_SOD, 2);
+
+       if( j2k->cstr_info && j2k->cur_tp_num==0){
+         j2k_add_tlmarker( j2k->curtileno, j2k->cstr_info, J2K_MS_SOD, cio_tell(cio), 0);
+       }
+
        if (j2k->curtileno == 0) {
                j2k->sod_start = cio_tell(cio) + j2k->pos_correction;
        }
 
        /* INDEX >> */
-       image_info = j2k->image_info;
-       if (image_info && image_info->index_on) {
-               if (!j2k->cur_tp_num ){
-                       image_info->tile[j2k->curtileno].end_header = cio_tell(cio) + j2k->pos_correction - 1;
-               }else{
-                       if(image_info->tile[j2k->curtileno].packet[image_info->num - 1].end_pos < cio_tell(cio))
-                               image_info->tile[j2k->curtileno].packet[image_info->num].start_pos = cio_tell(cio);
+       cstr_info = j2k->cstr_info;
+       if (cstr_info) {
+               if (!j2k->cur_tp_num ) {
+                       cstr_info->tile[j2k->curtileno].end_header = cio_tell(cio) + j2k->pos_correction - 1;
+                       j2k->cstr_info->tile[j2k->curtileno].tileno = j2k->curtileno;
                }
+               else{
+                       if(cstr_info->tile[j2k->curtileno].packet[cstr_info->packno - 1].end_pos < cio_tell(cio))
+                               cstr_info->tile[j2k->curtileno].packet[cstr_info->packno].start_pos = cio_tell(cio);
+               }
+               /* UniPG>> */
+#ifdef USE_JPWL
+               /* update markers struct */
+               j2k_add_marker(j2k->cstr_info, J2K_MS_SOD, j2k->sod_start, 2);
+#endif /* USE_JPWL */
+               /* <<UniPG */
        }
        /* << INDEX */
        
        tcp = &cp->tcps[j2k->curtileno];
        for (layno = 0; layno < tcp->numlayers; layno++) {
-               tcp->rates[layno] -= tcp->rates[layno] ? (j2k->sod_start / (cp->th * cp->tw)) : 0;
+               if (tcp->rates[layno]>(j2k->sod_start / (cp->th * cp->tw))) {
+                       tcp->rates[layno]-=(j2k->sod_start / (cp->th * cp->tw));
+               } else if (tcp->rates[layno]) {
+                       tcp->rates[layno]=1;
+               }
        }
-       if(image_info && (j2k->cur_tp_num == 0)) {
-               image_info->num = 0;
+       if(j2k->cur_tp_num == 0){
+               tcd->tcd_image->tiles->packno = 0;
+               if(cstr_info)
+                       cstr_info->packno = 0;
        }
        
-       l = tcd_encode_tile(tcd, j2k->curtileno, cio_getbp(cio), cio_numbytesleft(cio) - 2, image_info);
+       l = tcd_encode_tile(tcd, j2k->curtileno, cio_getbp(cio), cio_numbytesleft(cio) - 2, cstr_info);
        
        /* Writing Psot in SOT marker */
        totlen = cio_tell(cio) + l - j2k->sot_start;
@@ -1397,26 +3514,31 @@ static void j2k_read_sod(opj_j2k_t *j2k) {
 
        opj_cio_t *cio = j2k->cio;
        int curtileno = j2k->curtileno;
+
+       /* Index */
+       if (j2k->cstr_info) {
+               j2k->cstr_info->tile[j2k->curtileno].tp[j2k->cur_tp_num].tp_end_header =
+                       cio_tell(cio) + j2k->pos_correction - 1;
+               if (j2k->cur_tp_num == 0)
+                       j2k->cstr_info->tile[j2k->curtileno].end_header = cio_tell(cio) + j2k->pos_correction - 1;
+               j2k->cstr_info->packno = 0;
+       }
        
        len = int_min(j2k->eot - cio_getbp(cio), cio_numbytesleft(cio) + 1);
-       
+
        if (len == cio_numbytesleft(cio) + 1) {
                truncate = 1;           /* Case of a truncate codestream */
-       }
-       
-       data = (unsigned char *) opj_malloc((j2k->tile_len[curtileno] + len) * sizeof(unsigned char));
+       }       
 
-       for (i = 0; i < j2k->tile_len[curtileno]; i++) {
-               data[i] = j2k->tile_data[curtileno][i];
-       }
+       data = j2k->tile_data[curtileno];
+       data = (unsigned char*) opj_realloc(data, (j2k->tile_len[curtileno] + len) * sizeof(unsigned char));
 
        data_ptr = data + j2k->tile_len[curtileno];
        for (i = 0; i < len; i++) {
                data_ptr[i] = cio_read(cio, 1);
        }
-       
+
        j2k->tile_len[curtileno] += len;
-       opj_free(j2k->tile_data[curtileno]);
        j2k->tile_data[curtileno] = data;
        
        if (!truncate) {
@@ -1424,25 +3546,98 @@ static void j2k_read_sod(opj_j2k_t *j2k) {
        } else {
                j2k->state = J2K_STATE_NEOC;    /* RAJOUTE !! */
        }
+       j2k->cur_tp_num++;
 }
 
-static void j2k_write_rgn(opj_j2k_t *j2k, int compno, int tileno) {
-       
-       opj_cp_t *cp = j2k->cp;
-       opj_tcp_t *tcp = &cp->tcps[tileno];
-       opj_cio_t *cio = j2k->cio;
-       int numcomps = j2k->image->numcomps;
-       
-       cio_write(cio, J2K_MS_RGN, 2);                                          /* RGN  */
-       cio_write(cio, numcomps <= 256 ? 5 : 6, 2);                     /* Lrgn */
-       cio_write(cio, compno, numcomps <= 256 ? 1 : 2);        /* Crgn */
-       cio_write(cio, 0, 1);                                                           /* Srgn */
-       cio_write(cio, tcp->tccps[compno].roishift, 1);         /* SPrgn */
-}
-
-static void j2k_read_rgn(opj_j2k_t *j2k) {
-       int len, compno, roisty;
-
+/**
+ * Reads a SOD marker (Start Of Data)
+ *
+ * @param      p_header_data   the data contained in the SOD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the SOD marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_sod_v2 (
+                                               opj_j2k_v2_t *p_j2k,
+                                               struct opj_stream_private *p_stream,
+                                               struct opj_event_mgr * p_manager
+                                       )
+{
+       OPJ_UINT32 l_current_read_size;
+       opj_codestream_index_t * l_cstr_index = 00;
+       OPJ_BYTE ** l_current_data = 00;
+       opj_tcp_v2_t * l_tcp = 00;
+       OPJ_UINT32 * l_tile_len = 00;
+
+       /* preconditions */
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       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;
+
+       l_current_data = &(l_tcp->m_data);
+       l_tile_len = &l_tcp->m_data_size;
+
+       if (! *l_current_data) {
+               *l_current_data = (OPJ_BYTE*) opj_malloc/*FIXME V2 -> my_opj_malloc*/(p_j2k->m_specific_param.m_decoder.m_sot_length);
+       }
+       else {
+               *l_current_data = (OPJ_BYTE*) opj_realloc/*FIXME V2 -> my_opj_realloc*/(*l_current_data, *l_tile_len + p_j2k->m_specific_param.m_decoder.m_sot_length);
+       }
+
+       if (*l_current_data == 00) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Cannot decode tile\n");
+               return OPJ_FALSE;
+       }
+
+#ifdef TODO_MSD /* FIXME */
+       l_cstr_index = p_j2k->cstr_index;
+
+       /* Index */
+       if (l_cstr_index) {
+               OPJ_SIZE_T l_current_pos = opj_stream_tell(p_stream)-1;
+               l_cstr_index->tile_index[p_j2k->m_current_tile_number].tp_index[p_j2k->m_specific_param.m_encoder.m_current_tile_part_number].tp_end_header = l_current_pos;
+
+               l_cstr_index->packno = 0;
+       }
+#endif
+
+       l_current_read_size = opj_stream_read_data(     p_stream,
+                                                                                               *l_current_data + *l_tile_len,
+                                                                                               p_j2k->m_specific_param.m_decoder.m_sot_length,
+                                                                                               p_manager);
+
+       if (l_current_read_size != p_j2k->m_specific_param.m_decoder.m_sot_length) {
+               p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_NEOC;
+       }
+       else {
+               p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_TPHSOT;
+       }
+
+       *l_tile_len +=  l_current_read_size;
+
+       return OPJ_TRUE;
+}
+
+
+static void j2k_write_rgn(opj_j2k_t *j2k, int compno, int tileno) {
+       opj_cp_t *cp = j2k->cp;
+       opj_tcp_t *tcp = &cp->tcps[tileno];
+       opj_cio_t *cio = j2k->cio;
+       int numcomps = j2k->image->numcomps;
+       
+       cio_write(cio, J2K_MS_RGN, 2);                                          /* RGN  */
+       cio_write(cio, numcomps <= 256 ? 5 : 6, 2);                     /* Lrgn */
+       cio_write(cio, compno, numcomps <= 256 ? 1 : 2);        /* Crgn */
+       cio_write(cio, 0, 1);                                                           /* Srgn */
+       cio_write(cio, tcp->tccps[compno].roishift, 1);         /* SPrgn */
+}
+
+static void j2k_read_rgn(opj_j2k_t *j2k) {
+       int len, compno, roisty;
+
        opj_cp_t *cp = j2k->cp;
        opj_tcp_t *tcp = j2k->state == J2K_STATE_TPH ? &cp->tcps[j2k->curtileno] : j2k->default_tcp;
        opj_cio_t *cio = j2k->cio;
@@ -1474,20 +3669,106 @@ static void j2k_write_eoc(opj_j2k_t *j2k) {
        opj_cio_t *cio = j2k->cio;
        /* opj_event_msg(j2k->cinfo, "%.8x: EOC\n", cio_tell(cio) + j2k->pos_correction); */
        cio_write(cio, J2K_MS_EOC, 2);
+
+/* UniPG>> */
+#ifdef USE_JPWL
+       /* update markers struct */
+       j2k_add_marker(j2k->cstr_info, J2K_MS_EOC, cio_tell(cio) - 2, 2);
+#endif /* USE_JPWL */
+/* <<UniPG */
+}
+
+/**
+ * Reads a RGN marker (Region Of Interest)
+ *
+ * @param      p_header_data   the data contained in the POC box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the POC marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_rgn_v2 (
+                                               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;
+       opj_image_t * l_image = 00;
+
+       opj_cp_v2_t *l_cp = 00;
+       opj_tcp_v2_t *l_tcp = 00;
+       OPJ_UINT32 l_comp_room, l_comp_no, l_roi_sty;
+
+       // preconditions
+       assert(p_header_data != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       l_image = p_j2k->m_image;
+       l_nb_comp = l_image->numcomps;
+
+       if (l_nb_comp <= 256) {
+               l_comp_room = 1; }
+       else {
+               l_comp_room = 2; }
+
+       if (p_header_size != 2 + l_comp_room) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading RGN marker\n");
+               return OPJ_FALSE;
+       }
+
+       l_cp = &(p_j2k->m_cp);
+       l_tcp = (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH) ?
+                               &l_cp->tcps[p_j2k->m_current_tile_number] :
+                               p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       opj_read_bytes(p_header_data,&l_comp_no,l_comp_room);           /* Crgn */
+       p_header_data+=l_comp_room;
+       opj_read_bytes(p_header_data,&l_roi_sty,1);                                     /* Srgn */
+       ++p_header_data;
+
+#ifdef USE_JPWL
+       if (l_cp->correct) {
+               /* totlen is negative or larger than the bytes left!!! */
+               if (l_comp_room >= l_nb_comp) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR,
+                               "JPWL: bad component number in RGN (%d when there are only %d)\n",
+                               l_comp_room, l_nb_comp);
+                       if (!JPWL_ASSUME || JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+               }
+       };
+#endif /* USE_JPWL */
+
+       opj_read_bytes(p_header_data,(OPJ_UINT32 *) (&(l_tcp->tccps[l_comp_no].roishift)),1);   /* SPrgn */
+       ++p_header_data;
+
+       return OPJ_TRUE;
+
 }
 
 static void j2k_read_eoc(opj_j2k_t *j2k) {
        int i, tileno;
+       opj_bool success;
 
        /* if packets should be decoded */
        if (j2k->cp->limit_decoding != DECODE_ALL_BUT_PACKETS) {
                opj_tcd_t *tcd = tcd_create(j2k->cinfo);
                tcd_malloc_decode(tcd, j2k->image, j2k->cp);
                for (i = 0; i < j2k->cp->tileno_size; i++) {
+                       tcd_malloc_decode_tile(tcd, j2k->image, j2k->cp, i, j2k->cstr_info);
                        tileno = j2k->cp->tileno[i];
-                       tcd_decode_tile(tcd, j2k->tile_data[tileno], j2k->tile_len[tileno], tileno);
+                       success = tcd_decode_tile(tcd, j2k->tile_data[tileno], j2k->tile_len[tileno], tileno, j2k->cstr_info);
                        opj_free(j2k->tile_data[tileno]);
                        j2k->tile_data[tileno] = NULL;
+                       tcd_free_decode_tile(tcd, i);
+                       if (success == OPJ_FALSE) {
+                               j2k->state |= J2K_STATE_ERR;
+                               break;
+                       }
                }
                tcd_free_decode(tcd);
                tcd_destroy(tcd);
@@ -1499,9 +3780,68 @@ static void j2k_read_eoc(opj_j2k_t *j2k) {
                        opj_free(j2k->tile_data[tileno]);
                        j2k->tile_data[tileno] = NULL;
                }
+       }       
+       if (j2k->state & J2K_STATE_ERR)
+               j2k->state = J2K_STATE_MT + J2K_STATE_ERR;
+       else
+               j2k->state = J2K_STATE_MT; 
+}
+
+/**
+ * Reads a EOC marker (End Of Codestream)
+ *
+ * @param      p_header_data   the data contained in the SOD box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the SOD marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_eoc_v2 (     opj_j2k_v2_t *p_j2k,
+                                                       struct opj_stream_private *p_stream,
+                                                       struct opj_event_mgr * p_manager )
+{
+       OPJ_UINT32 i;
+       opj_tcd_v2_t * l_tcd = 00;
+       OPJ_UINT32 l_nb_tiles;
+       opj_tcp_v2_t * l_tcp = 00;
+       opj_bool l_success;
+
+       /* preconditions */
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       assert(p_stream != 00);
+
+       l_nb_tiles = p_j2k->m_cp.th * p_j2k->m_cp.tw;
+       l_tcp = p_j2k->m_cp.tcps;
+
+       l_tcd = tcd_create_v2(OPJ_TRUE);
+       if (l_tcd == 00) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Cannot decode tile, memory error\n");
+               return OPJ_FALSE;
        }
-       
-       j2k->state = J2K_STATE_MT;
+
+       for (i = 0; i < l_nb_tiles; ++i) {
+               if (l_tcp->m_data) {
+                       if (! tcd_init_decode_tile(l_tcd, i)) {
+                               tcd_destroy_v2(l_tcd);
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Cannot decode tile, memory error\n");
+                               return OPJ_FALSE;
+                       }
+
+                       l_success = tcd_decode_tile_v2(l_tcd, l_tcp->m_data, l_tcp->m_data_size, i, p_j2k->cstr_index);
+                       /* cleanup */
+
+                       if (! l_success) {
+                               p_j2k->m_specific_param.m_decoder.m_state |= J2K_STATE_ERR;
+                               break;
+                       }
+               }
+
+               j2k_tcp_destroy(l_tcp);
+               ++l_tcp;
+       }
+
+       tcd_destroy_v2(l_tcd);
+       return OPJ_TRUE;
 }
 
 typedef struct opj_dec_mstabent {
@@ -1540,6 +3880,10 @@ opj_dec_mstabent_t j2k_dec_mstab[] = {
   {J2K_MS_ESD, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_esd},
   {J2K_MS_RED, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_red},
 #endif /* USE_JPWL */
+#ifdef USE_JPSEC
+  {J2K_MS_SEC, J2K_STATE_MH, j2k_read_sec},
+  {J2K_MS_INSEC, 0, j2k_read_insec},
+#endif /* USE_JPSEC */
 
   {0, J2K_STATE_MH | J2K_STATE_TPH, j2k_read_unk}
 };
@@ -1606,6 +3950,71 @@ static void j2k_read_unk(opj_j2k_t *j2k) {
 
 }
 
+/**
+ * Reads an unknown marker
+ *
+ * @param      p_stream                                the stream object to read from.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_manager               the user event manager.
+ *
+ * @return     true                    if the marker could be deduced.
+*/
+opj_bool j2k_read_unk_v2 (     opj_j2k_v2_t *p_j2k,
+                                                       struct opj_stream_private *p_stream,
+                                                       OPJ_UINT32 *output_marker,
+                                                       struct opj_event_mgr * p_manager
+                                                       )
+{
+       OPJ_UINT32 l_unknown_marker;
+       const opj_dec_memory_marker_handler_t * l_marker_handler;
+       OPJ_UINT32 l_size_unk = 2;
+
+       // preconditions
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       assert(p_stream != 00);
+
+       opj_event_msg_v2(p_manager, EVT_WARNING, "Unknown marker\n");
+
+       while(1) {
+               // Try to read 2 bytes (the next marker ID) from stream and copy them into the buffer
+               if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                       return OPJ_FALSE;
+               }
+
+               // read 2 bytes as the new marker ID
+               opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_unknown_marker,2);
+
+               if (!(l_unknown_marker < 0xff00)) {
+
+                       // Get the marker handler from the marker ID
+                       l_marker_handler = j2k_get_marker_handler(l_unknown_marker);
+
+                       if (!(p_j2k->m_specific_param.m_decoder.m_state & l_marker_handler->states)) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Marker is not compliant with its position\n");
+                               return OPJ_FALSE;
+                       }
+                       else {
+                               if (l_marker_handler->id != J2K_MS_UNK) {
+                                       /* Add the marker to the codestream index*/
+                                       if (l_marker_handler->id != J2K_MS_SOT)
+                                               j2k_add_mhmarker_v2(p_j2k->cstr_index, J2K_MS_UNK,
+                                                                                       (OPJ_UINT32) opj_stream_tell(p_stream) - l_size_unk,
+                                                                                       l_size_unk);
+                                       break; /* next marker is known and well located */
+                               }
+                               else
+                                       l_size_unk += 2;
+                       }
+               }
+       }
+
+       *output_marker = l_marker_handler->id ;
+
+       return OPJ_TRUE;
+}
+
 /**
 Read the lookup table containing all the marker, status and action
 @param id Marker value
@@ -1625,15 +4034,19 @@ static opj_dec_mstabent_t *j2k_dec_mstab_lookup(int id) {
 /* ----------------------------------------------------------------------- */
 
 opj_j2k_t* j2k_create_decompress(opj_common_ptr cinfo) {
-       opj_j2k_t *j2k = (opj_j2k_t*)opj_malloc(sizeof(opj_j2k_t));
-       if(j2k) {
-               j2k->cinfo = cinfo;
-               j2k->default_tcp = (opj_tcp_t*)opj_malloc(sizeof(opj_tcp_t));
-               if(!j2k->default_tcp) {
-                       opj_free(j2k);
-                       return NULL;
-               }
+       opj_j2k_t *j2k = (opj_j2k_t*) opj_calloc(1, sizeof(opj_j2k_t));
+       if(!j2k)
+               return NULL;
+
+       j2k->default_tcp = (opj_tcp_t*) opj_calloc(1, sizeof(opj_tcp_t));
+       if(!j2k->default_tcp) {
+               opj_free(j2k);
+               return NULL;
        }
+
+       j2k->cinfo = cinfo;
+       j2k->tile_data = NULL;
+
        return j2k;
 }
 
@@ -1681,14 +4094,13 @@ void j2k_destroy_decompress(opj_j2k_t *j2k) {
 
                opj_free(cp);
        }
-
        opj_free(j2k);
 }
 
 void j2k_setup_decoder(opj_j2k_t *j2k, opj_dparameters_t *parameters) {
        if(j2k && parameters) {
                /* create and initialize the coding parameters structure */
-               opj_cp_t *cp = (opj_cp_t*)opj_malloc(sizeof(opj_cp_t));
+               opj_cp_t *cp = (opj_cp_t*) opj_calloc(1, sizeof(opj_cp_t));
                cp->reduce = parameters->cp_reduce;     
                cp->layer = parameters->cp_layer;
                cp->limit_decoding = parameters->cp_limit_decoding;
@@ -1705,12 +4117,29 @@ void j2k_setup_decoder(opj_j2k_t *j2k, opj_dparameters_t *parameters) {
        }
 }
 
-opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio) {
+void j2k_setup_decoder_v2(opj_j2k_v2_t *j2k, opj_dparameters_t *parameters)
+{
+       if(j2k && parameters) {
+               j2k->m_cp.m_specific_param.m_dec.m_layer = parameters->cp_layer;
+               j2k->m_cp.m_specific_param.m_dec.m_reduce = parameters->cp_reduce;
+
+#ifdef USE_JPWL
+               j2k->m_cp.correct = parameters->jpwl_correct;
+               j2k->m_cp.exp_comps = parameters->jpwl_exp_comps;
+               j2k->m_cp.max_tiles = parameters->jpwl_max_tiles;
+#endif /* USE_JPWL */
+       }
+}
+
+opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio, opj_codestream_info_t *cstr_info) {
        opj_image_t *image = NULL;
 
-       opj_common_ptr cinfo = j2k->cinfo;
+       opj_common_ptr cinfo = j2k->cinfo;      
 
        j2k->cio = cio;
+       j2k->cstr_info = cstr_info;
+       if (cstr_info)
+               memset(cstr_info, 0, sizeof(opj_codestream_info_t));
 
        /* create an empty image */
        image = opj_image_create0();
@@ -1722,13 +4151,12 @@ opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio) {
                opj_dec_mstabent_t *e;
                int id = cio_read(cio, 2);
 
-
 #ifdef USE_JPWL
                /* we try to honor JPWL correction power */
                if (j2k->cp->correct) {
 
                        int orig_pos = cio_tell(cio);
-                       bool status;
+                       opj_bool status;
 
                        /* call the corrector */
                        status = jpwl_correct(j2k);
@@ -1782,6 +4210,9 @@ opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio) {
                if (e->handler) {
                        (*e->handler)(j2k);
                }
+               if (j2k->state & J2K_STATE_ERR) 
+                       return NULL;    
+
                if (j2k->state == J2K_STATE_MT) {
                        break;
                }
@@ -1796,7 +4227,6 @@ opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio) {
        if (j2k->state != J2K_STATE_MT) {
                opj_event_msg(cinfo, EVT_WARNING, "Incomplete bitstream\n");
        }
-
        return image;
 }
 
@@ -1804,17 +4234,19 @@ opj_image_t* j2k_decode(opj_j2k_t *j2k, opj_cio_t *cio) {
 * Read a JPT-stream and decode file
 *
 */
-opj_image_t* j2k_decode_jpt_stream(opj_j2k_t *j2k, opj_cio_t *cio) {
+opj_image_t* j2k_decode_jpt_stream(opj_j2k_t *j2k, opj_cio_t *cio,  opj_codestream_info_t *cstr_info) {
        opj_image_t *image = NULL;
        opj_jpt_msg_header_t header;
        int position;
-
        opj_common_ptr cinfo = j2k->cinfo;
-       
+
+       OPJ_ARG_NOT_USED(cstr_info);
+
        j2k->cio = cio;
 
        /* create an empty image */
        image = opj_image_create0();
+       j2k->image = image;
 
        j2k->state = J2K_STATE_MHSOC;
        
@@ -1887,31 +4319,54 @@ opj_image_t* j2k_decode_jpt_stream(opj_j2k_t *j2k, opj_cio_t *cio) {
 /* ----------------------------------------------------------------------- */
 
 opj_j2k_t* j2k_create_compress(opj_common_ptr cinfo) {
-       opj_j2k_t *j2k = (opj_j2k_t*)opj_malloc(sizeof(opj_j2k_t));
+       opj_j2k_t *j2k = (opj_j2k_t*) opj_calloc(1, sizeof(opj_j2k_t));
        if(j2k) {
                j2k->cinfo = cinfo;
        }
        return j2k;
 }
 
+opj_j2k_v2_t* j2k_create_compress_v2()
+{
+       opj_j2k_v2_t *l_j2k = (opj_j2k_v2_t*) opj_malloc(sizeof(opj_j2k_v2_t));
+       if (!l_j2k) {
+               return NULL;
+       }
+
+       memset(l_j2k,0,sizeof(opj_j2k_v2_t));
+
+       l_j2k->m_is_decoder = 0;
+       l_j2k->m_cp.m_is_decoder = 0;
+
+       l_j2k->m_specific_param.m_encoder.m_header_tile_data = (OPJ_BYTE *) opj_malloc(J2K_DEFAULT_HEADER_SIZE);
+       if (! l_j2k->m_specific_param.m_encoder.m_header_tile_data) {
+               j2k_destroy(l_j2k);
+               return NULL;
+       }
+
+       l_j2k->m_specific_param.m_encoder.m_header_tile_data_size = J2K_DEFAULT_HEADER_SIZE;
+
+       // validation list creation
+       l_j2k->m_validation_list = opj_procedure_list_create();
+       if (! l_j2k->m_validation_list) {
+               j2k_destroy(l_j2k);
+               return NULL;
+       }
+
+       // execution list creation
+       l_j2k->m_procedure_list = opj_procedure_list_create();
+       if (! l_j2k->m_procedure_list) {
+               j2k_destroy(l_j2k);
+               return NULL;
+       }
+
+       return l_j2k;
+}
+
 void j2k_destroy_compress(opj_j2k_t *j2k) {
        int tileno;
 
        if(!j2k) return;
-
-       if(j2k->image_info != NULL) {
-               opj_image_info_t *image_info = j2k->image_info;
-               if (image_info->index_on && j2k->cp) {
-                       opj_cp_t *cp = j2k->cp;
-                       for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
-                               opj_tile_info_t *tile_info = &image_info->tile[tileno];
-                               opj_free(tile_info->thresh);
-                               opj_free(tile_info->packet);
-                       }
-                       opj_free(image_info->tile);
-               }
-               opj_free(image_info);
-       }
        if(j2k->cp != NULL) {
                opj_cp_t *cp = j2k->cp;
 
@@ -1940,7 +4395,7 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
        }
 
        /* create and initialize the coding parameters structure */
-       cp = (opj_cp_t*)opj_malloc(sizeof(opj_cp_t));
+       cp = (opj_cp_t*) opj_calloc(1, sizeof(opj_cp_t));
 
        /* keep a link to cp so that we can destroy it later in j2k_destroy_compress */
        j2k->cp = cp;
@@ -1966,12 +4421,6 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                memcpy(cp->matrice, parameters->cp_matrice, array_size);
        }
 
-       /* creation of an index file ? */
-       cp->index_on = parameters->index_on;
-       if(cp->index_on) {
-               j2k->image_info = (opj_image_info_t*)opj_malloc(sizeof(opj_image_info_t));
-       }
-
        /* tiles */
        cp->tdx = parameters->cp_tdx;
        cp->tdy = parameters->cp_tdy;
@@ -2020,12 +4469,12 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                int i;
 
                /* set JPWL on */
-               cp->epc_on = true;
-               cp->info_on = false; /* no informative technique */
+               cp->epc_on = OPJ_TRUE;
+               cp->info_on = OPJ_FALSE; /* no informative technique */
 
                /* set EPB on */
                if ((parameters->jpwl_hprot_MH > 0) || (parameters->jpwl_hprot_TPH[0] > 0)) {
-                       cp->epb_on = true;
+                       cp->epb_on = OPJ_TRUE;
                        
                        cp->hprot_MH = parameters->jpwl_hprot_MH;
                        for (i = 0; i < JPWL_MAX_NO_TILESPECS; i++) {
@@ -2046,7 +4495,7 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
 
                /* set ESD writing */
                if ((parameters->jpwl_sens_size == 1) || (parameters->jpwl_sens_size == 2)) {
-                       cp->esd_on = true;
+                       cp->esd_on = OPJ_TRUE;
 
                        cp->sens_size = parameters->jpwl_sens_size;
                        cp->sens_addr = parameters->jpwl_sens_addr;
@@ -2060,17 +4509,17 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                }
 
                /* always set RED writing to false: we are at the encoder */
-               cp->red_on = false;
+               cp->red_on = OPJ_FALSE;
 
        } else {
-               cp->epc_on = false;
+               cp->epc_on = OPJ_FALSE;
        }
 #endif /* USE_JPWL */
 
 
        /* initialize the mutiple tiles */
        /* ---------------------------- */
-       cp->tcps = (opj_tcp_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tcp_t));
+       cp->tcps = (opj_tcp_t*) opj_calloc(cp->tw * cp->th, sizeof(opj_tcp_t));
 
        for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
                opj_tcp_t *tcp = &cp->tcps[tileno];
@@ -2098,7 +4547,6 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                if (parameters->numpocs) {
                        /* initialisation of POC */
                        tcp->POC = 1;
-                       j2k_check_poc_val(parameters, image->numcomps, tcp->numlayers);
                        for (i = 0; i < parameters->numpocs; i++) {
                                if((tileno == parameters->POC[i].tile - 1) || (parameters->POC[i].tile == -1)) {
                                        opj_poc_t *tcp_poc = &tcp->pocs[numpocs_tile];
@@ -2117,8 +4565,8 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
                        tcp->numpocs = 0;
                }
 
-               tcp->tccps = (opj_tccp_t *) opj_malloc(image->numcomps * sizeof(opj_tccp_t));
-               
+               tcp->tccps = (opj_tccp_t*) opj_calloc(image->numcomps, sizeof(opj_tccp_t));
+
                for (i = 0; i < image->numcomps; i++) {
                        opj_tccp_t *tccp = &tcp->tccps[i];
                        tccp->csty = parameters->csty & 0x01;   /* 0 => one precinct || 1 => custom precinct  */
@@ -2196,241 +4644,8 @@ void j2k_setup_encoder(opj_j2k_t *j2k, opj_cparameters_t *parameters, opj_image_
        }
 }
 
-/**
-Create an index file
-@param j2k
-@param cio
-@param image_info
-@param index Index filename
-@return Returns 1 if successful, returns 0 otherwise
-*/
-static int j2k_create_index(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_info_t *image_info, char *index) {
-       int tileno, compno, layno, resno, precno, pack_nb, x, y;
-       FILE *stream = NULL;
-       double total_disto = 0;
-
-       image_info->codestream_size = cio_tell(cio) + j2k->pos_correction;      /* Correction 14/4/03 suite rmq de Patrick */
-
-
-#ifdef USE_JPWL
-       /* if JPWL is enabled and the name coincides with our own set
-          then discard the creation of the file: this was just done to
-          enable indexing, we do not want an index file
-       */
-       if (j2k->cp->epc_on && !strcmp(index, JPWL_PRIVATEINDEX_NAME))
-               return 1;
-#endif /* USE_JPWL */
-
-
-       stream = fopen(index, "w");
-       if (!stream) {
-               opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to open %s for writing\n", index);
-               return 0;
-       }
-       
-       fprintf(stream, "%d %d\n", image_info->image_w, image_info->image_h);
-       fprintf(stream, "%d\n", image_info->prog);
-       fprintf(stream, "%d %d\n", image_info->tile_x, image_info->tile_y);
-       fprintf(stream, "%d %d\n", image_info->tw, image_info->th);
-       fprintf(stream, "%d\n", image_info->comp);
-       fprintf(stream, "%d\n", image_info->layer);
-       fprintf(stream, "%d\n", image_info->decomposition);
-       
-       for (resno = image_info->decomposition; resno >= 0; resno--) {
-               fprintf(stream, "[%d,%d] ", 
-                       (1 << image_info->tile[0].pdx[resno]), (1 << image_info->tile[0].pdx[resno]));  /* based on tile 0 */
-       }
-       fprintf(stream, "\n");
-       fprintf(stream, "%d\n", image_info->main_head_end);
-       fprintf(stream, "%d\n", image_info->codestream_size);
-       
-       for (tileno = 0; tileno < image_info->tw * image_info->th; tileno++) {
-               fprintf(stream, "%4d %9d %9d %9d %9e %9d %9e\n",
-                       image_info->tile[tileno].num_tile,
-                       image_info->tile[tileno].start_pos,
-                       image_info->tile[tileno].end_header,
-                       image_info->tile[tileno].end_pos,
-                       image_info->tile[tileno].distotile, image_info->tile[tileno].nbpix,
-                       image_info->tile[tileno].distotile / image_info->tile[tileno].nbpix);
-       }
-       
-       for (tileno = 0; tileno < image_info->tw * image_info->th; tileno++) {
-               int start_pos, end_pos;
-               double disto = 0;
-               pack_nb = 0;
-                               
-               if (image_info->prog == LRCP) { /* LRCP */
-
-                       fprintf(stream, "pack_nb tileno layno resno compno precno start_pos  end_pos disto\n");
-
-                       for (layno = 0; layno < image_info->layer; layno++) {
-                               for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                                       for (compno = 0; compno < image_info->comp; compno++) {
-                                               int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                               for (precno = 0; precno < prec_max; precno++) {
-                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                       fprintf(stream, "%4d %6d %7d %5d %6d %6d %9d %9d %8e\n",
-                                                               pack_nb, tileno, layno, resno, compno, precno, start_pos, end_pos, disto);
-                                                       total_disto += disto;
-                                                       pack_nb++;
-                                               }
-                                       }
-                               }
-                       }
-               } /* LRCP */
-               else if (image_info->prog == RLCP) {    /* RLCP */
-
-                       fprintf(stream, "pack_nb tileno resno layno compno precno start_pos  end_pos disto\n");
-
-                       for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                               for (layno = 0; layno < image_info->layer; layno++) {
-                                       for (compno = 0; compno < image_info->comp; compno++) {
-                                               int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                               for (precno = 0; precno < prec_max; precno++) {
-                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                       fprintf(stream, "%4d %6d %5d %7d %6d %6d %9d %9d %8e\n",
-                                                               pack_nb, tileno, resno, layno, compno, precno, start_pos, end_pos, disto);
-                                                       total_disto += disto;
-                                                       pack_nb++;
-                                               }
-                                       }
-                               }
-                       }
-               } /* RLCP */
-               else if (image_info->prog == RPCL) {    /* RPCL */
-
-                       fprintf(stream, "\npack_nb tileno resno precno compno layno start_pos  end_pos disto\n"); 
-
-                       for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                               /* I suppose components have same XRsiz, YRsiz */
-                               int x0 = image_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tw * image_info->tile_x;
-                               int y0 = image_info->tile_Ox + (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tile_y;
-                               int x1 = x0 + image_info->tile_x;
-                               int y1 = y0 + image_info->tile_y;
-                               for (compno = 0; compno < image_info->comp; compno++) {
-                                       int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                       for (precno = 0; precno < prec_max; precno++) {
-                                               int pcnx = image_info->tile[tileno].pw[resno];
-                                               int pcx = (int) pow( 2, image_info->tile[tileno].pdx[resno] + image_info->decomposition - resno );
-                                               int pcy = (int) pow( 2, image_info->tile[tileno].pdy[resno] + image_info->decomposition - resno );
-                                               int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
-                                               int precno_y = (int) floor( (float)precno/(float)pcnx );
-                                               for(y = y0; y < y1; y++) {                                                      
-                                                       if (precno_y*pcy == y ) {
-                                                               for (x = x0; x < x1; x++) {                                                                     
-                                                                       if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < image_info->layer; layno++) {
-                                                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                                                       fprintf(stream, "%4d %6d %5d %6d %6d %7d %9d %9d %8e\n",
-                                                                                               pack_nb, tileno, resno, precno, compno, layno, start_pos, end_pos, disto); 
-                                                                                       total_disto += disto;
-                                                                                       pack_nb++; 
-                                                                               }
-                                                                       }
-                                                               }/* x = x0..x1 */
-                                                       } 
-                                               }  /* y = y0..y1 */
-                                       } /* precno */
-                               } /* compno */
-                       } /* resno */
-               } /* RPCL */
-               else if (image_info->prog == PCRL) {    /* PCRL */
-                       /* I suppose components have same XRsiz, YRsiz */
-                       int x0 = image_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tw * image_info->tile_x;
-                       int y0 = image_info->tile_Ox + (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tile_y;
-                       int x1 = x0 + image_info->tile_x;
-                       int y1 = y0 + image_info->tile_y;
-
-                       fprintf(stream, "\npack_nb tileno precno compno resno layno start_pos  end_pos disto\n"); 
-
-                       for (compno = 0; compno < image_info->comp; compno++) {
-                               for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                                       int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                       for (precno = 0; precno < prec_max; precno++) {
-                                               int pcnx = image_info->tile[tileno].pw[resno];
-                                               int pcx = (int) pow( 2, image_info->tile[tileno].pdx[resno] + image_info->decomposition - resno );
-                                               int pcy = (int) pow( 2, image_info->tile[tileno].pdy[resno] + image_info->decomposition - resno );
-                                               int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
-                                               int precno_y = (int) floor( (float)precno/(float)pcnx );
-                                               for(y = y0; y < y1; y++) {                                                      
-                                                       if (precno_y*pcy == y ) {
-                                                               for (x = x0; x < x1; x++) {                                                                     
-                                                                       if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < image_info->layer; layno++) {
-                                                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                                                       fprintf(stream, "%4d %6d %6d %6d %5d %7d %9d %9d %8e\n",
-                                                                                               pack_nb, tileno, precno, compno, resno, layno, start_pos, end_pos, disto); 
-                                                                                       total_disto += disto;
-                                                                                       pack_nb++; 
-                                                                               }
-                                                                       }
-                                                               }/* x = x0..x1 */
-                                                       } 
-                                               }  /* y = y0..y1 */
-                                       } /* precno */
-                               } /* resno */
-                       } /* compno */
-               } /* PCRL */
-               else {  /* CPRL */
-
-                       fprintf(stream, "\npack_nb tileno compno precno resno layno start_pos  end_pos disto\n"); 
-
-                       for (compno = 0; compno < image_info->comp; compno++) {
-                               /* I suppose components have same XRsiz, YRsiz */
-                               int x0 = image_info->tile_Ox + tileno - (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tw * image_info->tile_x;
-                               int y0 = image_info->tile_Ox + (int)floor( (float)tileno/(float)image_info->tw ) * image_info->tile_y;
-                               int x1 = x0 + image_info->tile_x;
-                               int y1 = y0 + image_info->tile_y;
-                               
-                               for (resno = 0; resno < image_info->decomposition + 1; resno++) {
-                                       int prec_max = image_info->tile[tileno].pw[resno] * image_info->tile[tileno].ph[resno];
-                                       for (precno = 0; precno < prec_max; precno++) {
-                                               int pcnx = image_info->tile[tileno].pw[resno];
-                                               int pcx = (int) pow( 2, image_info->tile[tileno].pdx[resno] + image_info->decomposition - resno );
-                                               int pcy = (int) pow( 2, image_info->tile[tileno].pdy[resno] + image_info->decomposition - resno );
-                                               int precno_x = precno - (int) floor( (float)precno/(float)pcnx ) * pcnx;
-                                               int precno_y = (int) floor( (float)precno/(float)pcnx );
-                                               for(y = y0; y < y1; y++) {
-                                                       if (precno_y*pcy == y ) {
-                                                               for (x = x0; x < x1; x++) {
-                                                                       if (precno_x*pcx == x ) {
-                                                                               for (layno = 0; layno < image_info->layer; layno++) {
-                                                                                       start_pos = image_info->tile[tileno].packet[pack_nb].start_pos;
-                                                                                       end_pos = image_info->tile[tileno].packet[pack_nb].end_pos;
-                                                                                       disto = image_info->tile[tileno].packet[pack_nb].disto;
-                                                                                       fprintf(stream, "%4d %6d %6d %6d %5d %7d %9d %9d %8e\n",
-                                                                                               pack_nb, tileno, compno, precno, resno, layno, start_pos, end_pos, disto); 
-                                                                                       total_disto += disto;
-                                                                                       pack_nb++; 
-                                                                               }
-                                                                       }
-                                                               }/* x = x0..x1 */
-                                                       }
-                                               } /* y = y0..y1 */
-                                       } /* precno */
-                               } /* resno */
-                       } /* compno */
-               } /* CPRL */   
-       } /* tileno */
-       
-       fprintf(stream, "%8e\n", image_info->D_max); /* SE max */
-       fprintf(stream, "%.8e\n", total_disto); /* SE totale */
-       fclose(stream);
-
-       return 1;
-}
-
-bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index) {
+opj_bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, opj_codestream_info_t *cstr_info) {
        int tileno, compno;
-       opj_image_info_t *image_info = NULL;
        opj_cp_t *cp = NULL;
 
        opj_tcd_t *tcd = NULL;  /* TCD component */
@@ -2440,34 +4655,39 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
 
        cp = j2k->cp;
 
-       /* j2k_dump_cp(stdout, image, cp); */
-
        /* INDEX >> */
-       image_info = j2k->image_info;
-       if (image_info && cp->index_on) {
-               image_info->index_on = cp->index_on;
-               image_info->tile = (opj_tile_info_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tile_info_t));
-               image_info->image_w = image->x1 - image->x0;
-               image_info->image_h = image->y1 - image->y0;
-               image_info->prog = (&cp->tcps[0])->prg;
-               image_info->tw = cp->tw;
-               image_info->th = cp->th;
-               image_info->tile_x = cp->tdx;   /* new version parser */
-               image_info->tile_y = cp->tdy;   /* new version parser */
-               image_info->tile_Ox = cp->tx0;  /* new version parser */
-               image_info->tile_Oy = cp->ty0;  /* new version parser */
-               image_info->comp = image->numcomps;
-               image_info->layer = (&cp->tcps[0])->numlayers;
-               image_info->decomposition = (&cp->tcps[0])->tccps->numresolutions - 1;
-               image_info->D_max = 0;          /* ADD Marcela */
+       j2k->cstr_info = cstr_info;
+       if (cstr_info) {
+               int compno;
+               cstr_info->tile = (opj_tile_info_t *) opj_malloc(cp->tw * cp->th * sizeof(opj_tile_info_t));
+               cstr_info->image_w = image->x1 - image->x0;
+               cstr_info->image_h = image->y1 - image->y0;
+               cstr_info->prog = (&cp->tcps[0])->prg;
+               cstr_info->tw = cp->tw;
+               cstr_info->th = cp->th;
+               cstr_info->tile_x = cp->tdx;    /* new version parser */
+               cstr_info->tile_y = cp->tdy;    /* new version parser */
+               cstr_info->tile_Ox = cp->tx0;   /* new version parser */
+               cstr_info->tile_Oy = cp->ty0;   /* new version parser */
+               cstr_info->numcomps = image->numcomps;
+               cstr_info->numlayers = (&cp->tcps[0])->numlayers;
+               cstr_info->numdecompos = (int*) opj_malloc(image->numcomps * sizeof(int));
+               for (compno=0; compno < image->numcomps; compno++) {
+                       cstr_info->numdecompos[compno] = (&cp->tcps[0])->tccps->numresolutions - 1;
+               }
+               cstr_info->D_max = 0.0;         /* ADD Marcela */
+               cstr_info->main_head_start = cio_tell(cio); /* position of SOC */
+               cstr_info->maxmarknum = 100;
+               cstr_info->marker = (opj_marker_info_t *) opj_malloc(cstr_info->maxmarknum * sizeof(opj_marker_info_t));
+               cstr_info->marknum = 0;
        }
        /* << INDEX */
-       
+
        j2k_write_soc(j2k);
        j2k_write_siz(j2k);
        j2k_write_cod(j2k);
        j2k_write_qcd(j2k);
-       
+
        if(cp->cinema){
                for (compno = 1; compno < image->numcomps; compno++) {
                        j2k_write_coc(j2k, compno);
@@ -2483,61 +4703,75 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
        if (cp->comment != NULL) {
                j2k_write_com(j2k);
        }
-       /* INDEX >> */
-       if(image_info && image_info->index_on) {
-               image_info->main_head_end = cio_tell(cio) - 1;
-       }
-       /* << INDEX */
 
        j2k->totnum_tp = j2k_calculate_tp(cp,image->numcomps,image,j2k);
-       /*       TLM Marker*/
+       /* TLM Marker*/
        if(cp->cinema){
                j2k_write_tlm(j2k);
                if (cp->cinema == CINEMA4K_24) {
                        j2k_write_poc(j2k);
                }
        }
+
+       /* uncomment only for testing JPSEC marker writing */
+       /* j2k_write_sec(j2k); */
+
+       /* INDEX >> */
+       if(cstr_info) {
+               cstr_info->main_head_end = cio_tell(cio) - 1;
+       }
+       /* << INDEX */
        /**** Main Header ENDS here ***/
-       
+
        /* create the tile encoder */
        tcd = tcd_create(j2k->cinfo);
 
        /* encode each tile */
-
        for (tileno = 0; tileno < cp->tw * cp->th; tileno++) {
                int pino;
                int tilepartno=0;
+               /* UniPG>> */
+               int acc_pack_num = 0;
+               /* <<UniPG */
+
 
                opj_tcp_t *tcp = &cp->tcps[tileno];
                opj_event_msg(j2k->cinfo, EVT_INFO, "tile number %d / %d\n", tileno + 1, cp->tw * cp->th);
-               
+
                j2k->curtileno = tileno;
                j2k->cur_tp_num = 0;
-
+               tcd->cur_totnum_tp = j2k->cur_totnum_tp[j2k->curtileno];
                /* initialisation before tile encoding  */
                if (tileno == 0) {
                        tcd_malloc_encode(tcd, image, cp, j2k->curtileno);
                } else {
                        tcd_init_encode(tcd, image, cp, j2k->curtileno);
                }
-               
+
                /* INDEX >> */
-               if(image_info && image_info->index_on) {
-                       image_info->tile[j2k->curtileno].num_tile = j2k->curtileno;
-                       image_info->tile[j2k->curtileno].start_pos = cio_tell(cio) + j2k->pos_correction;
+               if(cstr_info) {
+                       cstr_info->tile[j2k->curtileno].start_pos = cio_tell(cio) + j2k->pos_correction;
+                       cstr_info->tile[j2k->curtileno].maxmarknum = 10;
+                       cstr_info->tile[j2k->curtileno].marker = (opj_marker_info_t *) opj_malloc(cstr_info->tile[j2k->curtileno].maxmarknum * sizeof(opj_marker_info_t));
+                       cstr_info->tile[j2k->curtileno].marknum = 0;
                }
                /* << INDEX */
 
                for(pino = 0; pino <= tcp->numpocs; pino++) {
                        int tot_num_tp;
                        tcd->cur_pino=pino;
-                       
+
                        /*Get number of tile parts*/
                        tot_num_tp = j2k_get_num_tp(cp,pino,tileno);
                        tcd->tp_pos = cp->tp_pos;
 
                        for(tilepartno = 0; tilepartno < tot_num_tp ; tilepartno++){
                                j2k->tp_num = tilepartno;
+                               /* INDEX >> */
+                               if(cstr_info)
+                                       cstr_info->tile[j2k->curtileno].tp[j2k->cur_tp_num].tp_start_pos =
+                                       cio_tell(cio) + j2k->pos_correction;
+                               /* << INDEX */
                                j2k_write_sot(j2k);
 
                                if(j2k->cur_tp_num == 0 && cp->cinema == 0){
@@ -2550,83 +4784,2162 @@ bool j2k_encode(opj_j2k_t *j2k, opj_cio_t *cio, opj_image_t *image, char *index)
                                        }
                                }
 
+                               /* INDEX >> */
+                               if(cstr_info)
+                                       cstr_info->tile[j2k->curtileno].tp[j2k->cur_tp_num].tp_end_header =
+                                       cio_tell(cio) + j2k->pos_correction + 1;
+                               /* << INDEX */
+
                                j2k_write_sod(j2k, tcd);
-                               j2k->cur_tp_num ++;
-                       }
-                       
+
+                               /* INDEX >> */
+                               if(cstr_info) {
+                                       cstr_info->tile[j2k->curtileno].tp[j2k->cur_tp_num].tp_end_pos =
+                                               cio_tell(cio) + j2k->pos_correction - 1;
+                                       cstr_info->tile[j2k->curtileno].tp[j2k->cur_tp_num].tp_start_pack =
+                                               acc_pack_num;
+                                       cstr_info->tile[j2k->curtileno].tp[j2k->cur_tp_num].tp_numpacks =
+                                               cstr_info->packno - acc_pack_num;
+                                       acc_pack_num = cstr_info->packno;
+                               }
+                               /* << INDEX */
+
+                               j2k->cur_tp_num++;
+                       }                       
                }
-               /* INDEX >> */
-               if(image_info && image_info->index_on) {
-                       image_info->tile[j2k->curtileno].end_pos = cio_tell(cio) + j2k->pos_correction - 1;
+               if(cstr_info) {
+                       cstr_info->tile[j2k->curtileno].end_pos = cio_tell(cio) + j2k->pos_correction - 1;
                }
-               /* << INDEX */
-               
-               
+
+
                /*
                if (tile->PPT) { // BAD PPT !!! 
-                       FILE *PPT_file;
-                       int i;
-                       PPT_file=fopen("PPT","rb");
-                       fprintf(stderr,"%c%c%c%c",255,97,tile->len_ppt/256,tile->len_ppt%256);
-                       for (i=0;i<tile->len_ppt;i++) {
-                               unsigned char elmt;
-                               fread(&elmt, 1, 1, PPT_file);
-                               fwrite(&elmt,1,1,f);
-                       }
-                       fclose(PPT_file);
-                       unlink("PPT");
+               FILE *PPT_file;
+               int i;
+               PPT_file=fopen("PPT","rb");
+               fprintf(stderr,"%c%c%c%c",255,97,tile->len_ppt/256,tile->len_ppt%256);
+               for (i=0;i<tile->len_ppt;i++) {
+               unsigned char elmt;
+               fread(&elmt, 1, 1, PPT_file);
+               fwrite(&elmt,1,1,f);
+               }
+               fclose(PPT_file);
+               unlink("PPT");
                }
                */
-               
+
        }
-       
+
        /* destroy the tile encoder */
        tcd_free_encode(tcd);
        tcd_destroy(tcd);
 
+       opj_free(j2k->cur_totnum_tp);
+
        j2k_write_eoc(j2k);
-       
-       /* Creation of the index file */
-       if(image_info && image_info->index_on) {
-               if(!j2k_create_index(j2k, cio, image_info, index)) {
-                       opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to create index file %s\n", index);
-                       return false;
-               }
-       }
 
+       if(cstr_info) {
+               cstr_info->codestream_size = cio_tell(cio) + j2k->pos_correction;
+               /* UniPG>> */
+               /* The following adjustment is done to adjust the codestream size */
+               /* if SOD is not at 0 in the buffer. Useful in case of JP2, where */
+               /* the first bunch of bytes is not in the codestream              */
+               cstr_info->codestream_size -= cstr_info->main_head_start;
+               /* <<UniPG */
+       }
 
 #ifdef USE_JPWL
        /*
-       preparation of JPWL marker segments: can be finalized only when the whole
-       codestream is known
+       preparation of JPWL marker segments
        */
-       if(image_info && image_info->index_on && cp->epc_on) {
+       if(cp->epc_on) {
+
+               /* encode according to JPWL */
+               jpwl_encode(j2k, cio, image);
 
-               /* let's begin creating a marker list, according to user wishes */
-               jpwl_prepare_marks(j2k, cio, image);
+       }
+#endif /* USE_JPWL */
 
-               /* now we dump the JPWL markers on the codestream */
-               jpwl_dump_marks(j2k, cio, image);
+       return OPJ_TRUE;
+}
 
-               /* do not know exactly what is this for,
-               but it gets called during index creation */
-               j2k->pos_correction = 0;
+static void j2k_add_mhmarker(opj_codestream_info_t *cstr_info, unsigned short int type, int pos, int len) {
 
-               /* Re-creation of the index file, with updated info */
-               if(image_info && image_info->index_on) {
-                       if(!j2k_create_index(j2k, cio, image_info, index)) {
-                               opj_event_msg(j2k->cinfo, EVT_ERROR, "failed to re-create index file %s\n", index);
-                               return false;
-                       }
-               }
+       if (!cstr_info)
+               return;
+
+       /* expand the list? */
+       if ((cstr_info->marknum + 1) > cstr_info->maxmarknum) {
+               cstr_info->maxmarknum = 100 + (int) ((float) cstr_info->maxmarknum * 1.0F);
+               cstr_info->marker = (opj_marker_info_t*)opj_realloc(cstr_info->marker, cstr_info->maxmarknum);
+       }
 
-               /* now we finalize the marker contents */
-               /*jpwl_finalize_marks(j2k, cio, image);*/
+       /* add the marker */
+       cstr_info->marker[cstr_info->marknum].type = type;
+       cstr_info->marker[cstr_info->marknum].pos = pos;
+       cstr_info->marker[cstr_info->marknum].len = len;
+       cstr_info->marknum++;
 
+}
+
+static void j2k_add_mhmarker_v2(opj_codestream_index_t *cstr_index, OPJ_UINT32 type, OPJ_UINT32 pos, OPJ_UINT32 len) {
+
+       if (!cstr_index)
+               return;
+
+       /* expand the list? */
+       if ((cstr_index->marknum + 1) > cstr_index->maxmarknum) {
+               cstr_index->maxmarknum = 100 + (int) ((float) cstr_index->maxmarknum * 1.0F);
+               cstr_index->marker = (opj_marker_info_t*)opj_realloc(cstr_index->marker, cstr_index->maxmarknum *sizeof(opj_marker_info_t));
        }
-#endif /* USE_JPWL */
 
-       
-       return true;
+       /* add the marker */
+       cstr_index->marker[cstr_index->marknum].type = (OPJ_UINT16)type;
+       cstr_index->marker[cstr_index->marknum].pos = (OPJ_INT32)pos;
+       cstr_index->marker[cstr_index->marknum].len = (OPJ_INT32)len;
+       cstr_index->marknum++;
+
 }
 
+static void j2k_add_tlmarker( int tileno, opj_codestream_info_t *cstr_info, unsigned short int type, int pos, int len) {
+
+
+  opj_marker_info_t *marker;
+
+       if (!cstr_info)
+               return;
+
+       /* expand the list? */
+       if ((cstr_info->tile[tileno].marknum + 1) > cstr_info->tile[tileno].maxmarknum) {
+               cstr_info->tile[tileno].maxmarknum = 100 + (int) ((float) cstr_info->tile[tileno].maxmarknum * 1.0F);
+               cstr_info->tile[tileno].marker = (opj_marker_info_t*)opj_realloc(cstr_info->tile[tileno].marker, cstr_info->maxmarknum);
+       }
+
+       marker = &(cstr_info->tile[tileno].marker[cstr_info->tile[tileno].marknum]);
+
+       /* add the marker */
+       marker->type = type;
+       marker->pos = pos;
+       marker->len = len;
+       cstr_info->tile[tileno].marknum++;
+}
+
+
+
+
+/*
+ * -----------------------------------------------------------------------
+ * -----------------------------------------------------------------------
+ * -----------------------------------------------------------------------
+ */
+
+/**
+ * Ends the decompression procedures and possibiliy add data to be read after the
+ * codestream.
+ */
+opj_bool j2k_end_decompress(
+                                               opj_j2k_v2_t *p_j2k,
+                                               opj_stream_private_t *p_stream,
+                                               opj_event_mgr_t * p_manager)
+{
+       return OPJ_TRUE;
+}
+
+/**
+ * Reads a jpeg2000 codestream header structure.
+
+ *
+ * @param p_stream the stream to read data from.
+ * @param p_j2k the jpeg2000 codec.
+ * @param p_manager the user event manager.
+ *
+ * @return true if the box is valid.
+ */
+opj_bool j2k_read_header(      struct opj_stream_private *p_stream,
+                                                       opj_j2k_v2_t* p_j2k,
+                                                       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 */
+       p_j2k->m_image = opj_image_create0();
+       if (! p_j2k->m_image) {
+               return OPJ_FALSE;
+       }
+
+       /* customization of the validation */
+       j2k_setup_decoding_validation(p_j2k);
+
+       /* validation of the parameters codec */
+       if (! j2k_exec(p_j2k, p_j2k->m_validation_list, p_stream,p_manager)) {
+               opj_image_destroy(p_j2k->m_image);
+               p_j2k->m_image = NULL;
+               return OPJ_FALSE;
+       }
+
+       /* customization of the encoding */
+       j2k_setup_header_reading(p_j2k);
+
+       /* read header */
+       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;
+       }
+
+       if (! j2k_copy_img_header(p_j2k, p_image)){
+               opj_image_destroy(p_j2k->m_image);
+               p_j2k->m_image = NULL;
+               return OPJ_FALSE;
+       }
+
+
+       return OPJ_TRUE;
+}
+
+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->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]),
+                               &(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->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,
+                               p_j2k->m_image->icc_profile_buf,
+                               p_j2k->m_image->icc_profile_len);
+       }
+       else
+               p_image->icc_profile_buf = NULL;
+
+       return OPJ_TRUE;
+}
+
+
+
+/**
+ * Sets up the procedures to do on reading header. Developpers wanting to extend the library can add their own reading procedures.
+ */
+void j2k_setup_header_reading (opj_j2k_v2_t *p_j2k)
+{
+       // preconditions
+       assert(p_j2k != 00);
+
+       opj_procedure_list_add_procedure(p_j2k->m_procedure_list,(void*)j2k_read_header_procedure);
+
+       /* DEVELOPER CORNER, add your custom procedures */
+       opj_procedure_list_add_procedure(p_j2k->m_procedure_list,(void*)j2k_copy_default_tcp_and_create_tcd);
+
+}
+
+/**
+ * Sets up the validation ,i.e. adds the procedures to lauch to make sure the codec parameters
+ * are valid. Developpers wanting to extend the library can add their own validation procedures.
+ */
+void j2k_setup_decoding_validation (opj_j2k_v2_t *p_j2k)
+{
+       // preconditions
+       assert(p_j2k != 00);
+
+       opj_procedure_list_add_procedure(p_j2k->m_validation_list, (void*)j2k_build_decoder);
+       opj_procedure_list_add_procedure(p_j2k->m_validation_list, (void*)j2k_decoding_validation);
+       /* DEVELOPER CORNER, add your custom validation procedure */
+
+}
+
+/**
+ * Builds the cp decoder parameters to use to decode tile.
+ */
+opj_bool j2k_build_decoder (
+                                               opj_j2k_v2_t * p_j2k,
+                                               opj_stream_private_t *p_stream,
+                                               opj_event_mgr_t * p_manager
+                                               )
+{
+       // add here initialization of cp
+       // copy paste of setup_decoder
+       return OPJ_TRUE;
+}
+
+/**
+ * The default decoding validation procedure without any extension.
+ *
+ * @param      p_j2k                   the jpeg2000 codec to validate.
+ * @param      p_stream                                the input stream to validate.
+ * @param      p_manager               the user event manager.
+ *
+ * @return true if the parameters are correct.
+ */
+opj_bool j2k_decoding_validation (
+                                                               opj_j2k_v2_t *p_j2k,
+                                                               opj_stream_private_t *p_stream,
+                                                               opj_event_mgr_t * p_manager
+                                                         )
+{
+       opj_bool l_is_valid = OPJ_TRUE;
+
+       // preconditions
+       assert(p_j2k != 00);
+       assert(p_stream != 00);
+       assert(p_manager != 00);
+
+
+       /* STATE checking */
+       /* make sure the state is at 0 */
+#ifdef TODO_MSD
+       l_is_valid &= (p_j2k->m_specific_param.m_decoder.m_state == J2K_DEC_STATE_NONE);
+#endif
+       l_is_valid &= (p_j2k->m_specific_param.m_decoder.m_state == 0x0000);
+
+       /* POINTER validation */
+       /* make sure a p_j2k codec is present */
+       /* make sure a procedure list is present */
+       l_is_valid &= (p_j2k->m_procedure_list != 00);
+       /* make sure a validation list is present */
+       l_is_valid &= (p_j2k->m_validation_list != 00);
+
+       /* PARAMETER VALIDATION */
+       return l_is_valid;
+}
+
+opj_bool j2k_read_header_procedure(    opj_j2k_v2_t *p_j2k,
+                                                                       struct opj_stream_private *p_stream,
+                                                                       struct opj_event_mgr * p_manager)
+{
+       OPJ_UINT32 l_current_marker;
+       OPJ_UINT32 l_marker_size;
+       const opj_dec_memory_marker_handler_t * l_marker_handler = 00;
+
+       /* preconditions */
+       assert(p_stream != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       /*  We enter in the main header */
+       p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_MHSOC;
+
+       /* Try to read the SOC marker, the codestream must begin with SOC marker */
+       if (! j2k_read_soc_v2(p_j2k,p_stream,p_manager)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Expected a SOC marker \n");
+               return OPJ_FALSE;
+       }
+
+       /* Try to read 2 bytes (the next marker ID) from stream and copy them into the buffer */
+       if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+               return OPJ_FALSE;
+       }
+
+       /* Read 2 bytes as the new marker ID */
+       opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_current_marker,2);
+
+       /* Try to read until the SOT is detected */
+       while (l_current_marker != J2K_MS_SOT) {
+
+               /* Check if the current marker ID is valid */
+               if (l_current_marker < 0xff00) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "We expected read a marker ID (0xff--) instead of %.8x\n", l_current_marker);
+                       return OPJ_FALSE;
+               }
+
+               /* Get the marker handler from the marker ID */
+               l_marker_handler = j2k_get_marker_handler(l_current_marker);
+
+               /* Manage case where marker is unknown */
+               if (l_marker_handler->id == J2K_MS_UNK) {
+                       if (! j2k_read_unk_v2(p_j2k, p_stream, &l_current_marker, p_manager)){
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Unknow marker have been detected and generated error.\n");
+                               return OPJ_FALSE;
+                       }
+
+                       if (l_current_marker == J2K_MS_SOT)
+                               break; /* SOT marker is detected main header is completely read */
+                       else    /* Get the marker handler from the marker ID */
+                               l_marker_handler = j2k_get_marker_handler(l_current_marker);
+               }
+
+               /* Check if the marker is known and if it is the right place to find it */
+               if (! (p_j2k->m_specific_param.m_decoder.m_state & l_marker_handler->states) ) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Marker is not compliant with its position\n");
+                       return OPJ_FALSE;
+               }
+
+               /* Try to read 2 bytes (the marker size) from stream and copy them into the buffer */
+               if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                       return OPJ_FALSE;
+               }
+
+               /* read 2 bytes as the marker size */
+               opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_marker_size,2);
+               l_marker_size -= 2; // Subtract the size of the marker ID already read
+
+               /* Check if the marker size is compatible with the header data size */
+               if (l_marker_size > p_j2k->m_specific_param.m_decoder.m_header_data_size) {
+                       p_j2k->m_specific_param.m_decoder.m_header_data = (OPJ_BYTE*)
+                                       opj_realloc(p_j2k->m_specific_param.m_decoder.m_header_data,l_marker_size);
+                       if (p_j2k->m_specific_param.m_decoder.m_header_data == 00) {
+                               return OPJ_FALSE;
+                       }
+                       p_j2k->m_specific_param.m_decoder.m_header_data_size = l_marker_size;
+               }
+
+               /* Try to read the rest of the marker segment from stream and copy them into the buffer */
+               if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,l_marker_size,p_manager) != l_marker_size) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                       return OPJ_FALSE;
+               }
+
+               /* Read the marker segment with the correct marker handler */
+               if (! (*(l_marker_handler->handler))(p_j2k,p_j2k->m_specific_param.m_decoder.m_header_data,l_marker_size,p_manager)) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Marker handler function failed to read the marker segment\n");
+                       return OPJ_FALSE;
+               }
+
+               /* Add the marker to the codestream index*/
+               j2k_add_mhmarker_v2(p_j2k->cstr_index,
+                                                       l_marker_handler->id,
+                                                       (OPJ_UINT32) opj_stream_tell(p_stream) - l_marker_size - 4,
+                                                       l_marker_size + 4 );
+
+               /* Try to read 2 bytes (the next marker ID) from stream and copy them into the buffer */
+               if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                       return OPJ_FALSE;
+               }
+
+               /* read 2 bytes as the new marker ID */
+               opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_current_marker,2);
+       }
+
+       opj_event_msg_v2(p_manager, EVT_INFO, "Main header has been correctly decode.\n");
+
+       /* Position of the last element if the main header */
+       p_j2k->cstr_index->main_head_end = (OPJ_UINT32) opj_stream_tell(p_stream) - 2;
+
+       /* Next step: read a tile-part header */
+       p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_TPHSOT;
+
+       return OPJ_TRUE;
+}
+
+/**
+ * Excutes the given procedures on the given codec.
+ *
+ * @param      p_procedure_list        the list of procedures to execute
+ * @param      p_j2k                                   the jpeg2000 codec to execute the procedures on.
+ * @param      p_stream                                        the stream to execute the procedures on.
+ * @param      p_manager                       the user manager.
+ *
+ * @return     true                            if all the procedures were successfully executed.
+ */
+opj_bool j2k_exec (    opj_j2k_v2_t * p_j2k,
+                                       opj_procedure_list_t * p_procedure_list,
+                                       opj_stream_private_t *p_stream,
+                                       opj_event_mgr_t * p_manager )
+{
+       opj_bool (** l_procedure) (opj_j2k_v2_t * ,opj_stream_private_t *,opj_event_mgr_t *) = 00;
+       opj_bool l_result = OPJ_TRUE;
+       OPJ_UINT32 l_nb_proc, i;
+
+       // preconditions
+       assert(p_procedure_list != 00);
+       assert(p_j2k != 00);
+       assert(p_stream != 00);
+       assert(p_manager != 00);
+
+
+       l_nb_proc = opj_procedure_list_get_nb_procedures(p_procedure_list);
+       l_procedure = (opj_bool (**) (opj_j2k_v2_t * ,opj_stream_private_t *,opj_event_mgr_t *)) opj_procedure_list_get_first_procedure(p_procedure_list);
+
+       for     (i=0;i<l_nb_proc;++i) {
+               l_result = l_result && ((*l_procedure) (p_j2k,p_stream,p_manager));
+               ++l_procedure;
+       }
+
+       // and clear the procedure list at the end.
+       opj_procedure_list_clear(p_procedure_list);
+       return l_result;
+}
+
+/* FIXME DOC*/
+opj_bool j2k_copy_default_tcp_and_create_tcd
+                                               (
+                                               opj_j2k_v2_t * p_j2k,
+                                               opj_stream_private_t *p_stream,
+                                               opj_event_mgr_t * p_manager
+                                               )
+{
+       opj_tcp_v2_t * l_tcp = 00;
+       opj_tcp_v2_t * l_default_tcp = 00;
+       OPJ_UINT32 l_nb_tiles;
+       OPJ_UINT32 i,j;
+       opj_tccp_t *l_current_tccp = 00;
+       OPJ_UINT32 l_tccp_size;
+       OPJ_UINT32 l_mct_size;
+       opj_image_t * l_image;
+       OPJ_UINT32 l_mcc_records_size,l_mct_records_size;
+       opj_mct_data_t * l_src_mct_rec, *l_dest_mct_rec;
+       opj_simple_mcc_decorrelation_data_t * l_src_mcc_rec, *l_dest_mcc_rec;
+       OPJ_UINT32 l_offset;
+
+       /* preconditions */
+       assert(p_j2k != 00);
+       assert(p_stream != 00);
+       assert(p_manager != 00);
+
+       l_image = p_j2k->m_image;
+       l_nb_tiles = p_j2k->m_cp.th * p_j2k->m_cp.tw;
+       l_tcp = p_j2k->m_cp.tcps;
+       l_tccp_size = l_image->numcomps * sizeof(opj_tccp_t);
+       l_default_tcp = p_j2k->m_specific_param.m_decoder.m_default_tcp;
+       l_mct_size = l_image->numcomps * l_image->numcomps * sizeof(OPJ_FLOAT32);
+
+       /* For each tile */
+       for (i=0; i<l_nb_tiles; ++i) {
+               /* keep the tile-compo coding parameters pointer of the current tile coding parameters*/
+               l_current_tccp = l_tcp->tccps;
+               /*Copy default coding parameters into the current tile coding parameters*/
+               memcpy(l_tcp, l_default_tcp, sizeof(opj_tcp_v2_t));
+               /* Initialize some values of the current tile coding parameters*/
+               l_tcp->ppt = 0;
+               l_tcp->ppt_data = 00;
+               /* Reconnect the tile-compo coding parameters pointer to the current tile coding parameters*/
+               l_tcp->tccps = l_current_tccp;
+
+               /* Get the mct_decoding_matrix of the dflt_tile_cp and copy them into the current tile cp*/
+               if (l_default_tcp->m_mct_decoding_matrix) {
+                       l_tcp->m_mct_decoding_matrix = (OPJ_FLOAT32*)opj_malloc(l_mct_size);
+                       if (! l_tcp->m_mct_decoding_matrix ) {
+                               return OPJ_FALSE;
+                       }
+                       memcpy(l_tcp->m_mct_decoding_matrix,l_default_tcp->m_mct_decoding_matrix,l_mct_size);
+               }
+
+               /* Get the mct_record of the dflt_tile_cp and copy them into the current tile cp*/
+               l_mct_records_size = l_default_tcp->m_nb_max_mct_records * sizeof(opj_mct_data_t);
+               l_tcp->m_mct_records = (opj_mct_data_t*)opj_malloc(l_mct_records_size);
+               if (! l_tcp->m_mct_records) {
+                       return OPJ_FALSE;
+               }
+               memcpy(l_tcp->m_mct_records, l_default_tcp->m_mct_records,l_mct_records_size);
+
+               /* Copy the mct record data from dflt_tile_cp to the current tile*/
+               l_src_mct_rec = l_default_tcp->m_mct_records;
+               l_dest_mct_rec = l_tcp->m_mct_records;
+
+               for (j=0;j<l_default_tcp->m_nb_mct_records;++j) {
+
+                       if (l_src_mct_rec->m_data) {
+
+                               l_dest_mct_rec->m_data = (OPJ_BYTE*) opj_malloc(l_src_mct_rec->m_data_size);
+                               if(! l_dest_mct_rec->m_data) {
+                                       return OPJ_FALSE;
+                               }
+                               memcpy(l_dest_mct_rec->m_data,l_src_mct_rec->m_data,l_src_mct_rec->m_data_size);
+                       }
+
+                       ++l_src_mct_rec;
+                       ++l_dest_mct_rec;
+               }
+
+               /* Get the mcc_record of the dflt_tile_cp and copy them into the current tile cp*/
+               l_mcc_records_size = l_default_tcp->m_nb_max_mcc_records * sizeof(opj_simple_mcc_decorrelation_data_t);
+               l_tcp->m_mcc_records = (opj_simple_mcc_decorrelation_data_t*) opj_malloc(l_mcc_records_size);
+               if (! l_tcp->m_mcc_records) {
+                       return OPJ_FALSE;
+               }
+               memcpy(l_tcp->m_mcc_records,l_default_tcp->m_mcc_records,l_mcc_records_size);
+
+               /* Copy the mcc record data from dflt_tile_cp to the current tile*/
+               l_src_mcc_rec = l_default_tcp->m_mcc_records;
+               l_dest_mcc_rec = l_tcp->m_mcc_records;
+
+               for (j=0;j<l_default_tcp->m_nb_max_mcc_records;++j) {
+
+                       if (l_src_mcc_rec->m_decorrelation_array) {
+                               l_offset = l_src_mcc_rec->m_decorrelation_array - l_default_tcp->m_mct_records;
+                               l_dest_mcc_rec->m_decorrelation_array = l_tcp->m_mct_records + l_offset;
+                       }
+
+                       if (l_src_mcc_rec->m_offset_array) {
+                               l_offset = l_src_mcc_rec->m_offset_array - l_default_tcp->m_mct_records;
+                               l_dest_mcc_rec->m_offset_array = l_tcp->m_mct_records + l_offset;
+                       }
+
+                       ++l_src_mcc_rec;
+                       ++l_dest_mcc_rec;
+               }
+
+               /* Copy all the dflt_tile_compo_cp to the current tile cp */
+               memcpy(l_current_tccp,l_default_tcp->tccps,l_tccp_size);
+
+               /* Move to next tile cp*/
+               ++l_tcp;
+       }
+
+       /* Create the current tile decoder*/
+       p_j2k->m_tcd = (opj_tcd_v2_t*)tcd_create_v2(OPJ_TRUE); // FIXME why a cast ?
+       if (! p_j2k->m_tcd ) {
+               return OPJ_FALSE;
+       }
+
+       if ( !tcd_init_v2(p_j2k->m_tcd, l_image, &(p_j2k->m_cp)) ) {
+               tcd_destroy_v2(p_j2k->m_tcd);
+               p_j2k->m_tcd = 00;
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Cannot decode tile, memory error\n");
+               return OPJ_FALSE;
+       }
+
+       return OPJ_TRUE;
+}
+
+/**
+ * Reads the lookup table containing all the marker, status and action, and returns the handler associated
+ * with the marker value.
+ * @param      p_id            Marker value to look up
+ *
+ * @return     the handler associated with the id.
+*/
+const opj_dec_memory_marker_handler_t * j2k_get_marker_handler (const OPJ_UINT32 p_id)
+{
+       const opj_dec_memory_marker_handler_t *e;
+       for (e = j2k_memory_marker_handler_tab; e->id != 0; ++e) {
+               if (e->id == p_id) {
+                       break; // we find a handler corresponding to the marker ID
+               }
+       }
+       return e;
+}
+
+
+/**
+ * Destroys a jpeg2000 codec.
+ *
+ * @param      p_j2k   the jpeg20000 structure to destroy.
+ */
+void j2k_destroy (opj_j2k_v2_t *p_j2k)
+{
+       if (p_j2k == 00) {
+               return;
+       }
+
+       if (p_j2k->m_is_decoder) {
+
+               if (p_j2k->m_specific_param.m_decoder.m_default_tcp != 00) {
+                       j2k_tcp_destroy(p_j2k->m_specific_param.m_decoder.m_default_tcp);
+                       opj_free(p_j2k->m_specific_param.m_decoder.m_default_tcp);
+                       p_j2k->m_specific_param.m_decoder.m_default_tcp = 00;
+               }
+
+               if (p_j2k->m_specific_param.m_decoder.m_header_data != 00) {
+                       opj_free(p_j2k->m_specific_param.m_decoder.m_header_data);
+                       p_j2k->m_specific_param.m_decoder.m_header_data = 00;
+                       p_j2k->m_specific_param.m_decoder.m_header_data_size = 0;
+               }
+       }
+       else {
+
+               if (p_j2k->m_specific_param.m_encoder.m_encoded_tile_data) {
+                       opj_free(p_j2k->m_specific_param.m_encoder.m_encoded_tile_data);
+                       p_j2k->m_specific_param.m_encoder.m_encoded_tile_data = 00;
+               }
+
+               if (p_j2k->m_specific_param.m_encoder.m_tlm_sot_offsets_buffer) {
+                       opj_free(p_j2k->m_specific_param.m_encoder.m_tlm_sot_offsets_buffer);
+                       p_j2k->m_specific_param.m_encoder.m_tlm_sot_offsets_buffer = 00;
+                       p_j2k->m_specific_param.m_encoder.m_tlm_sot_offsets_current = 00;
+               }
+
+               if (p_j2k->m_specific_param.m_encoder.m_header_tile_data) {
+                       opj_free(p_j2k->m_specific_param.m_encoder.m_header_tile_data);
+                       p_j2k->m_specific_param.m_encoder.m_header_tile_data = 00;
+                       p_j2k->m_specific_param.m_encoder.m_header_tile_data_size = 0;
+               }
+       }
+
+       tcd_destroy_v2(p_j2k->m_tcd);
+
+       j2k_cp_destroy(&(p_j2k->m_cp));
+       memset(&(p_j2k->m_cp),0,sizeof(opj_cp_t));
+
+       opj_procedure_list_destroy(p_j2k->m_procedure_list);
+       p_j2k->m_procedure_list = 00;
+
+       opj_procedure_list_destroy(p_j2k->m_validation_list);
+       p_j2k->m_procedure_list = 00;
+
+       j2k_destroy_cstr_index(p_j2k->cstr_index);
+       p_j2k->cstr_index = 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->marker) {
+               opj_free(p_cstr_ind->marker);
+               p_cstr_ind->marker = NULL;
+       }
+
+       if (p_cstr_ind->tile_index) {
+               // FIXME
+       }
+}
+
+
+
+/**
+ * Destroys a tile coding parameter structure.
+ *
+ * @param      p_tcp           the tile coding parameter to destroy.
+ */
+void j2k_tcp_destroy (opj_tcp_v2_t *p_tcp)
+{
+       if (p_tcp == 00) {
+               return;
+       }
+
+       if (p_tcp->ppt_buffer != 00) {
+               opj_free(p_tcp->ppt_buffer);
+               p_tcp->ppt_buffer = 00;
+       }
+
+       if (p_tcp->tccps != 00) {
+               opj_free(p_tcp->tccps);
+               p_tcp->tccps = 00;
+       }
+
+       if (p_tcp->m_mct_coding_matrix != 00) {
+               opj_free(p_tcp->m_mct_coding_matrix);
+               p_tcp->m_mct_coding_matrix = 00;
+       }
+
+       if (p_tcp->m_mct_decoding_matrix != 00) {
+               opj_free(p_tcp->m_mct_decoding_matrix);
+               p_tcp->m_mct_decoding_matrix = 00;
+       }
+
+       if (p_tcp->m_mcc_records) {
+               opj_free(p_tcp->m_mcc_records);
+               p_tcp->m_mcc_records = 00;
+               p_tcp->m_nb_max_mcc_records = 0;
+               p_tcp->m_nb_mcc_records = 0;
+       }
+
+       if (p_tcp->m_mct_records) {
+               opj_mct_data_t * l_mct_data = p_tcp->m_mct_records;
+               OPJ_UINT32 i;
+
+               for (i=0;i<p_tcp->m_nb_mct_records;++i) {
+                       if (l_mct_data->m_data) {
+                               opj_free(l_mct_data->m_data);
+                               l_mct_data->m_data = 00;
+                       }
+
+                       ++l_mct_data;
+               }
+
+               opj_free(p_tcp->m_mct_records);
+               p_tcp->m_mct_records = 00;
+       }
+
+       if (p_tcp->mct_norms != 00) {
+               opj_free(p_tcp->mct_norms);
+               p_tcp->mct_norms = 00;
+       }
+
+       if (p_tcp->m_data) {
+               opj_free(p_tcp->m_data);
+               p_tcp->m_data = 00;
+       }
+}
+
+
+/**
+ * Destroys a coding parameter structure.
+ *
+ * @param      p_cp            the coding parameter to destroy.
+ */
+void j2k_cp_destroy (opj_cp_v2_t *p_cp)
+{
+       OPJ_UINT32 l_nb_tiles;
+       opj_tcp_v2_t * l_current_tile = 00;
+       OPJ_UINT32 i;
+
+       if
+               (p_cp == 00)
+       {
+               return;
+       }
+       if
+               (p_cp->tcps != 00)
+       {
+               l_current_tile = p_cp->tcps;
+               l_nb_tiles = p_cp->th * p_cp->tw;
+
+               for
+                       (i = 0; i < l_nb_tiles; ++i)
+               {
+                       j2k_tcp_destroy(l_current_tile);
+                       ++l_current_tile;
+               }
+               opj_free(p_cp->tcps);
+               p_cp->tcps = 00;
+       }
+       if
+               (p_cp->ppm_buffer != 00)
+       {
+               opj_free(p_cp->ppm_buffer);
+               p_cp->ppm_buffer = 00;
+       }
+       if
+               (p_cp->comment != 00)
+       {
+               opj_free(p_cp->comment);
+               p_cp->comment = 00;
+       }
+       if
+               (! p_cp->m_is_decoder)
+       {
+               if
+                       (p_cp->m_specific_param.m_enc.m_matrice)
+               {
+                       opj_free(p_cp->m_specific_param.m_enc.m_matrice);
+                       p_cp->m_specific_param.m_enc.m_matrice = 00;
+               }
+       }
+}
+
+
+
+/**
+ * Reads a tile header.
+ * @param      p_j2k           the jpeg2000 codec.
+ * @param      p_stream                        the stream to write data to.
+ * @param      p_manager       the user event manager.
+ */
+opj_bool j2k_read_tile_header( opj_j2k_v2_t * p_j2k,
+                                                               OPJ_UINT32 * p_tile_index,
+                                                               OPJ_UINT32 * p_data_size,
+                                                               OPJ_INT32 * p_tile_x0, OPJ_INT32 * p_tile_y0,
+                                                               OPJ_INT32 * p_tile_x1, OPJ_INT32 * p_tile_y1,
+                                                               OPJ_UINT32 * p_nb_comps,
+                                                               opj_bool * p_go_on,
+                                                               opj_stream_private_t *p_stream,
+                                                               opj_event_mgr_t * p_manager )
+{
+       OPJ_UINT32 l_current_marker = J2K_MS_SOT;
+       OPJ_UINT32 l_marker_size;
+       const opj_dec_memory_marker_handler_t * l_marker_handler = 00;
+       opj_tcp_v2_t * l_tcp = NULL;
+       OPJ_UINT32 l_nb_tiles;
+
+       /* preconditions */
+       assert(p_stream != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       /* Reach the End Of Codestream ?*/
+       if (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_EOC){
+               l_current_marker = J2K_MS_EOC;
+       }
+       /* We need to encounter a SOT marker (a new tile-part header) */
+       else if (p_j2k->m_specific_param.m_decoder.m_state != J2K_STATE_TPHSOT){
+               return OPJ_FALSE;
+       }
+
+       /* Read into the codestream until reach the EOC or ! can_decode ??? FIXME */
+       while ( (!p_j2k->m_specific_param.m_decoder.m_can_decode) && (l_current_marker != J2K_MS_EOC) ) {
+
+               /* Try to read until the Start Of Data is detected */
+               while (l_current_marker != J2K_MS_SOD) {
+
+                       /* Try to read 2 bytes (the marker size) from stream and copy them into the buffer */
+                       if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                               return OPJ_FALSE;
+                       }
+
+                       /* Read 2 bytes from the buffer as the marker size */
+                       opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_marker_size,2);
+
+                       /* Why this condition? FIXME */
+                       if (p_j2k->m_specific_param.m_decoder.m_state & J2K_STATE_TPH){
+                               p_j2k->m_specific_param.m_decoder.m_sot_length -= (l_marker_size + 2);
+                       }
+                       l_marker_size -= 2; /* Subtract the size of the marker ID already read */
+
+                       /* Get the marker handler from the marker ID */
+                       l_marker_handler = j2k_get_marker_handler(l_current_marker);
+
+                       /* Check if the marker is known and if it is the right place to find it */
+                       if (! (p_j2k->m_specific_param.m_decoder.m_state & l_marker_handler->states) ) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Marker is not compliant with its position\n");
+                               return OPJ_FALSE;
+                       }
+/* FIXME manage case of unknown marker as in the main header ? */
+
+                       /* Check if the marker size is compatible with the header data size */
+                       if (l_marker_size > p_j2k->m_specific_param.m_decoder.m_header_data_size) {
+                               p_j2k->m_specific_param.m_decoder.m_header_data = (OPJ_BYTE*)
+                                       opj_realloc(p_j2k->m_specific_param.m_decoder.m_header_data,l_marker_size);
+                               if (p_j2k->m_specific_param.m_decoder.m_header_data == 00) {
+                                       return OPJ_FALSE;
+                               }
+                               p_j2k->m_specific_param.m_decoder.m_header_data_size = l_marker_size;
+                       }
+
+                       /* Try to read the rest of the marker segment from stream and copy them into the buffer */
+                       if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,l_marker_size,p_manager) != l_marker_size) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                               return OPJ_FALSE;
+                       }
+
+                       /* Read the marker segment with the correct marker handler */
+                       if (! (*(l_marker_handler->handler))(p_j2k,p_j2k->m_specific_param.m_decoder.m_header_data,l_marker_size,p_manager)) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Fail to read the current marker segment (%#x)\n", l_current_marker);
+                               return OPJ_FALSE;
+                       }
+
+                       if (p_j2k->m_specific_param.m_decoder.m_skip_data) {
+                               /* Skip the rest of the tile part header*/
+                               if (opj_stream_skip(p_stream,p_j2k->m_specific_param.m_decoder.m_sot_length,p_manager) != p_j2k->m_specific_param.m_decoder.m_sot_length) {
+                                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                                       return OPJ_FALSE;
+                               }
+                               l_current_marker = J2K_MS_SOD; /* Normally we reached a SOD */
+                       }
+                       else {
+                               /* Try to read 2 bytes (the next marker ID) from stream and copy them into the buffer*/
+                               if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+                                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                                       return OPJ_FALSE;
+                               }
+                               /* Read 2 bytes from the buffer as the new marker ID */
+                               opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_current_marker,2);
+                       }
+               }
+
+               /* If we didn't skip data before, we need to read the SOD marker*/
+               if (! p_j2k->m_specific_param.m_decoder.m_skip_data) {
+                       /* Try to read the SOD marker and skip data ? FIXME */
+                       if (! j2k_read_sod_v2(p_j2k, p_stream, p_manager)) {
+                               return OPJ_FALSE;
+                       }
+
+                       if (! p_j2k->m_specific_param.m_decoder.m_can_decode){
+                               /* Try to read 2 bytes (the next marker ID) from stream and copy them into the buffer */
+                               if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+                                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                                       return OPJ_FALSE;
+                               }
+
+                               /* Read 2 bytes from buffer as the new marker ID */
+                               opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_current_marker,2);
+                       }
+               }
+               else {
+                       /* Indicate we will try to read a new tile-part header*/
+                       p_j2k->m_specific_param.m_decoder.m_skip_data = 0;
+                       p_j2k->m_specific_param.m_decoder.m_can_decode = 0;
+                       p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_TPHSOT;
+
+                       /* Try to read 2 bytes (the next marker ID) from stream and copy them into the buffer */
+                       if (opj_stream_read_data(p_stream,p_j2k->m_specific_param.m_decoder.m_header_data,2,p_manager) != 2) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                               return OPJ_FALSE;
+                       }
+
+                       /* Read 2 bytes from buffer as the new marker ID */
+                       opj_read_bytes(p_j2k->m_specific_param.m_decoder.m_header_data,&l_current_marker,2);
+               }
+       }
+
+       /* Current marker is the EOC marker ?*/
+       if (l_current_marker == J2K_MS_EOC) {
+               if (p_j2k->m_specific_param.m_decoder.m_state != J2K_STATE_EOC ){
+                       p_j2k->m_current_tile_number = 0;
+                       p_j2k->m_specific_param.m_decoder.m_state = J2K_STATE_EOC;
+               }
+       }
+
+       /* FIXME DOC ???*/
+       if ( ! p_j2k->m_specific_param.m_decoder.m_can_decode) {
+               l_tcp = p_j2k->m_cp.tcps + p_j2k->m_current_tile_number;
+               l_nb_tiles = p_j2k->m_cp.th * p_j2k->m_cp.tw;
+
+               while( (p_j2k->m_current_tile_number < l_nb_tiles) && (l_tcp->m_data == 00) ) {
+                       ++p_j2k->m_current_tile_number;
+                       ++l_tcp;
+               }
+
+               if (p_j2k->m_current_tile_number == l_nb_tiles) {
+                       *p_go_on = OPJ_FALSE;
+                       return OPJ_TRUE;
+               }
+       }
+
+       /*FIXME ???*/
+       if (! tcd_init_decode_tile(p_j2k->m_tcd, p_j2k->m_current_tile_number)) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Cannot decode tile, memory error\n");
+               return OPJ_FALSE;
+       }
+
+       opj_event_msg_v2(p_manager, EVT_INFO, "Header of tile %d / %d has been read.\n",
+                       p_j2k->m_current_tile_number +1, p_j2k->m_cp.th * p_j2k->m_cp.tw);
+
+       *p_tile_index = p_j2k->m_current_tile_number;
+       *p_go_on = OPJ_TRUE;
+       *p_data_size = tcd_get_decoded_tile_size(p_j2k->m_tcd);
+       *p_tile_x0 = p_j2k->m_tcd->tcd_image->tiles->x0;
+       *p_tile_y0 = p_j2k->m_tcd->tcd_image->tiles->y0;
+       *p_tile_x1 = p_j2k->m_tcd->tcd_image->tiles->x1;
+       *p_tile_y1 = p_j2k->m_tcd->tcd_image->tiles->y1;
+       *p_nb_comps = p_j2k->m_tcd->tcd_image->tiles->numcomps;
+
+        p_j2k->m_specific_param.m_decoder.m_state |= 0x0080;// FIXME J2K_DEC_STATE_DATA;
+
+       return OPJ_TRUE;
+}
+
+
+opj_bool j2k_decode_tile (     opj_j2k_v2_t * p_j2k,
+                                                       OPJ_UINT32 p_tile_index,
+                                                       OPJ_BYTE * p_data,
+                                                       OPJ_UINT32 p_data_size,
+                                                       opj_stream_private_t *p_stream,
+                                                       opj_event_mgr_t * p_manager )
+{
+       OPJ_UINT32 l_current_marker;
+       OPJ_BYTE l_data [2];
+       opj_tcp_v2_t * l_tcp;
+
+       /* preconditions */
+       assert(p_stream != 00);
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+
+       if ( !(p_j2k->m_specific_param.m_decoder.m_state & 0x0080/*FIXME J2K_DEC_STATE_DATA*/)
+               || (p_tile_index != p_j2k->m_current_tile_number) ) {
+               return OPJ_FALSE;
+       }
+
+       l_tcp = &(p_j2k->m_cp.tcps[p_tile_index]);
+       if (! l_tcp->m_data) {
+               j2k_tcp_destroy(&(p_j2k->m_cp.tcps[p_tile_index]));
+               return OPJ_FALSE;
+       }
+
+       if (! tcd_decode_tile_v2(       p_j2k->m_tcd,
+                                                               l_tcp->m_data,
+                                                               l_tcp->m_data_size,
+                                                               p_tile_index,
+                                                               p_j2k->cstr_index) ) {
+               j2k_tcp_destroy(l_tcp);
+               p_j2k->m_specific_param.m_decoder.m_state |= 0x8000;//FIXME J2K_DEC_STATE_ERR;
+               return OPJ_FALSE;
+       }
+
+       if (! tcd_update_tile_data(p_j2k->m_tcd,p_data,p_data_size)) {
+               return OPJ_FALSE;
+       }
+
+       j2k_tcp_destroy(l_tcp);
+       p_j2k->m_tcd->tcp = 0;
+
+       p_j2k->m_specific_param.m_decoder.m_can_decode = 0;
+       p_j2k->m_specific_param.m_decoder.m_state &= (~ (0x0080));// FIXME J2K_DEC_STATE_DATA);
+
+       if (p_j2k->m_specific_param.m_decoder.m_state != 0x0100){ //FIXME J2K_DEC_STATE_EOC)
+               if (opj_stream_read_data(p_stream,l_data,2,p_manager) != 2) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short\n");
+                       return OPJ_FALSE;
+               }
+
+               opj_read_bytes(l_data,&l_current_marker,2);
+
+               if (l_current_marker == J2K_MS_EOC) {
+                       p_j2k->m_current_tile_number = 0;
+                       p_j2k->m_specific_param.m_decoder.m_state =  0x0100;//FIXME J2K_DEC_STATE_EOC;
+               }
+               else if (l_current_marker != J2K_MS_SOT)
+               {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Stream too short, expected SOT\n");
+                       return OPJ_FALSE;
+               }
+       }
+
+       return OPJ_TRUE;
+}
+
+
+opj_bool j2k_update_image_data (opj_tcd_v2_t * p_tcd, OPJ_BYTE * p_data)
+{
+       OPJ_UINT32 i,j,k = 0;
+       OPJ_UINT32 l_width,l_height,l_offset_x,l_offset_y;
+       opj_image_comp_t * l_img_comp = 00;
+       opj_tcd_tilecomp_v2_t * l_tilec = 00;
+       opj_image_t * l_image = 00;
+       OPJ_UINT32 l_size_comp, l_remaining;
+       OPJ_UINT32 l_dest_stride;
+       OPJ_INT32 * l_dest_ptr;
+       opj_tcd_resolution_v2_t* l_res= 00;
+
+
+       l_tilec = p_tcd->tcd_image->tiles->comps;
+       l_image = p_tcd->image;
+       l_img_comp = l_image->comps;
+
+       for (i=0;i<p_tcd->image->numcomps;++i) {
+
+               if (!l_img_comp->data) {
+
+                       l_img_comp->data = (OPJ_INT32*) opj_malloc(l_img_comp->w * l_img_comp->h * sizeof(OPJ_INT32));
+                       if (! l_img_comp->data) {
+                               return OPJ_FALSE;
+                       }
+                       memset(l_img_comp->data,0,l_img_comp->w * l_img_comp->h * sizeof(OPJ_INT32));
+               }
+
+               l_size_comp = l_img_comp->prec >> 3; /*(/ 8)*/
+               l_remaining = l_img_comp->prec & 7;  /* (%8) */
+               l_res = l_tilec->resolutions + l_img_comp->resno_decoded;
+
+               if (l_remaining) {
+                       ++l_size_comp;
+               }
+
+               if (l_size_comp == 3) {
+                       l_size_comp = 4;
+               }
+
+               l_width = (l_res->x1 - l_res->x0);
+               l_height = (l_res->y1 - l_res->y0);
+
+               l_dest_stride = (l_img_comp->w) - l_width;
+
+               l_offset_x = int_ceildivpow2(l_img_comp->x0, l_img_comp->factor);
+               l_offset_y = int_ceildivpow2(l_img_comp->y0, l_img_comp->factor);
+               l_dest_ptr = l_img_comp->data + (l_res->x0 - l_offset_x) + (l_res->y0 - l_offset_y) * l_img_comp->w;
+
+               switch (l_size_comp) {
+                       case 1:
+                               {
+                                       OPJ_CHAR * l_src_ptr = (OPJ_CHAR*) p_data;
+                                       if (l_img_comp->sgnd) {
+                                               for (j=0;j<l_height;++j) {
+                                                       for (k=0;k<l_width;++k) {
+                                                               *(l_dest_ptr++) = (OPJ_INT32) (*(l_src_ptr++));
+                                                       }
+
+                                                       l_dest_ptr += l_dest_stride;
+                                               }
+                                       }
+                                       else {
+                                               for (j=0;j<l_height;++j) {
+                                                       for (k=0;k<l_width;++k) {
+                                                               *(l_dest_ptr++) = (OPJ_INT32) ((*(l_src_ptr++))&0xff);
+                                                       }
+
+                                                       l_dest_ptr += l_dest_stride;
+                                               }
+                                       }
+
+                                       p_data = (OPJ_BYTE*) l_src_ptr;
+                               }
+                               break;
+                       case 2:
+                               {
+                                       OPJ_INT16 * l_src_ptr = (OPJ_INT16 *) p_data;
+
+                                       if (l_img_comp->sgnd) {
+                                               for (j=0;j<l_height;++j) {
+                                                       for (k=0;k<l_width;++k) {
+                                                               *(l_dest_ptr++) = *(l_src_ptr++);
+                                                       }
+
+                                                       l_dest_ptr += l_dest_stride;
+                                               }
+                                       }
+                                       else {
+                                               for (j=0;j<l_height;++j) {
+                                                       for (k=0;k<l_width;++k) {
+                                                               *(l_dest_ptr++) = (*(l_src_ptr++))&0xffff;
+                                                       }
+
+                                                       l_dest_ptr += l_dest_stride;
+                                               }
+                                       }
+
+                                       p_data = (OPJ_BYTE*) l_src_ptr;
+                               }
+                               break;
+                       case 4:
+                               {
+                                       OPJ_INT32 * l_src_ptr = (OPJ_INT32 *) p_data;
+                                       for (j=0;j<l_height;++j) {
+                                               for (k=0;k<l_width;++k) {
+                                                       *(l_dest_ptr++) = (*(l_src_ptr++));
+                                               }
+
+                                               l_dest_ptr += l_dest_stride;
+                                       }
+
+                                       p_data = (OPJ_BYTE*) l_src_ptr;
+                               }
+                               break;
+               }
+
+               ++l_img_comp;
+               ++l_tilec;
+       }
+       return OPJ_TRUE;
+}
+
+/**
+ * Sets the given area to be decoded. This function should be called right after opj_read_header and before any tile header reading.
+ *
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_start_x               the left position of the rectangle to decode (in image coordinates).
+ * @param      p_end_x                 the right position of the rectangle to decode (in image coordinates).
+ * @param      p_start_y               the up position of the rectangle to decode (in image coordinates).
+ * @param      p_end_y                 the bottom position of the rectangle to decode (in image coordinates).
+ * @param      p_manager               the user event manager
+ *
+ * @return     true                    if the area could be set.
+ */
+opj_bool j2k_set_decode_area(  opj_j2k_v2_t *p_j2k,
+                                                               OPJ_INT32 p_start_x, OPJ_INT32 p_start_y,
+                                                               OPJ_INT32 p_end_x, OPJ_INT32 p_end_y,
+                                                               struct opj_event_mgr * p_manager )
+{
+       opj_cp_v2_t * l_cp = &(p_j2k->m_cp);
+
+       /* Check if we are read the main header */
+       if (p_j2k->m_specific_param.m_decoder.m_state != J2K_STATE_TPHSOT) { // FIXME J2K_DEC_STATE_TPHSOT)
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Need to decode the main header before begin to decode the remaining codestream");
+               return OPJ_FALSE;
+       }
+
+       /* ----- */
+       /* Check if the positions provided by the user are correct */
+
+       /* Left */
+       if (p_start_x > l_cp->tx0 + l_cp->tw * l_cp->tdx) {
+               opj_event_msg_v2(p_manager, EVT_ERROR,
+                       "Left position of the decoded area (ROI_x0=%d) is outside the tile area (XTOsiz + nb_tw x XTsiz=%d).\n",
+                       p_start_x, l_cp->tx0 + l_cp->tw * l_cp->tdx);
+               return OPJ_FALSE;
+       }
+       else if (p_start_x < l_cp->tx0){
+               opj_event_msg_v2(p_manager, EVT_WARNING,
+                               "Left position of the decoded area (ROI_x0=%d) is outside the tile area (XTOsiz=%d).\n",
+                               p_start_x, l_cp->tx0);
+               p_j2k->m_specific_param.m_decoder.m_start_tile_x = 0;
+       }
+       else
+               p_j2k->m_specific_param.m_decoder.m_start_tile_x = (p_start_x - l_cp->tx0) / l_cp->tdx;
+
+       /* Up */
+       if (p_start_y > l_cp->ty0 + l_cp->th * l_cp->tdy){
+               opj_event_msg_v2(p_manager, EVT_ERROR,
+                               "Up position of the decoded area (ROI_y0=%d) is outside the tile area (YTOsiz + nb_th x YTsiz=%d).\n",
+                               p_start_y, l_cp->ty0 + l_cp->th * l_cp->tdy);
+               return OPJ_FALSE;
+       }
+       else if (p_start_y < l_cp->ty0){
+               opj_event_msg_v2(p_manager, EVT_WARNING,
+                               "Up position of the decoded area (ROI_y0=%d) is outside the tile area (YTOsiz=%d).\n",
+                               p_start_y, l_cp->ty0);
+               p_j2k->m_specific_param.m_decoder.m_start_tile_y = 0;
+       }
+       else
+               p_j2k->m_specific_param.m_decoder.m_start_tile_y = (p_start_y - l_cp->ty0) / l_cp->tdy;
+
+       /* Right */
+       if (p_end_x < l_cp->tx0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR,
+                       "Right position of the decoded area (ROI_x1=%d) is outside the tile area (XTOsiz=%d).\n",
+                       p_end_x, l_cp->tx0);
+               return OPJ_FALSE;
+       }
+       else if (p_end_x > l_cp->tx0 + l_cp->tw * l_cp->tdx) {
+               opj_event_msg_v2(p_manager, EVT_WARNING,
+                       "Right position of the decoded area (ROI_x1=%d) is outside the tile area (XTOsiz + nb_tw x XTsiz=%d).\n",
+                       p_end_x, l_cp->tx0 + l_cp->tw * l_cp->tdx);
+               p_j2k->m_specific_param.m_decoder.m_end_tile_x = l_cp->tw; // FIXME (-1) ???
+       }
+       else
+               p_j2k->m_specific_param.m_decoder.m_end_tile_x = int_ceildiv((p_end_x - l_cp->tx0), l_cp->tdx);
+
+       /* Bottom */
+       if (p_end_y < l_cp->ty0) {
+               opj_event_msg_v2(p_manager, EVT_ERROR,
+                       "Right position of the decoded area (ROI_y1=%d) is outside the tile area (YTOsiz=%d).\n",
+                       p_end_y, l_cp->ty0);
+               return OPJ_FALSE;
+       }
+       if (p_end_y > l_cp->ty0 + l_cp->th * l_cp->tdy){
+               opj_event_msg_v2(p_manager, EVT_WARNING,
+                       "Bottom position of the decoded area (ROI_y1=%d) is outside the tile area (YTOsiz + nb_th x YTsiz=%d).\n",
+                       p_end_y, l_cp->ty0 + l_cp->th * l_cp->tdy);
+               p_j2k->m_specific_param.m_decoder.m_start_tile_y = l_cp->th; // FIXME (-1) ???
+       }
+       else
+               p_j2k->m_specific_param.m_decoder.m_end_tile_y = int_ceildiv((p_end_y - l_cp->ty0), l_cp->tdy);
+       /* ----- */
+
+       p_j2k->m_specific_param.m_decoder.m_discard_tiles = 1;
+
+       return OPJ_TRUE;
+}
+
+
+/* ----------------------------------------------------------------------- */
+/* J2K / JPT decoder interface                                             */
+/* ----------------------------------------------------------------------- */
+/**
+ * Creates a J2K decompression structure.
+ *
+ * @return a handle to a J2K decompressor if successful, NULL otherwise.
+*/
+opj_j2k_v2_t* j2k_create_decompress_v2()
+{
+       opj_j2k_v2_t *l_j2k = (opj_j2k_v2_t*) opj_malloc(sizeof(opj_j2k_v2_t));
+       if (!l_j2k) {
+               return 00;
+       }
+       memset(l_j2k,0,sizeof(opj_j2k_v2_t));
+
+       l_j2k->m_is_decoder = 1;
+       l_j2k->m_cp.m_is_decoder = 1;
+
+       l_j2k->m_specific_param.m_decoder.m_default_tcp = (opj_tcp_v2_t*) opj_malloc(sizeof(opj_tcp_v2_t));
+       if (!l_j2k->m_specific_param.m_decoder.m_default_tcp) {
+               j2k_destroy(l_j2k);
+               return 00;
+       }
+       memset(l_j2k->m_specific_param.m_decoder.m_default_tcp,0,sizeof(opj_tcp_v2_t));
+
+       l_j2k->m_specific_param.m_decoder.m_header_data = (OPJ_BYTE *) opj_malloc(J2K_DEFAULT_HEADER_SIZE);
+       if (! l_j2k->m_specific_param.m_decoder.m_header_data) {
+               j2k_destroy(l_j2k);
+               return 00;
+       }
+
+       l_j2k->m_specific_param.m_decoder.m_header_data_size = J2K_DEFAULT_HEADER_SIZE;
+
+       /* codestream index creation */
+       l_j2k->cstr_index = j2k_create_cstr_index();
+
+                       /*(opj_codestream_index_t*) opj_malloc(sizeof(opj_codestream_index_t));
+       if (!l_j2k->cstr_index){
+               j2k_destroy(l_j2k);
+               return NULL;
+       }
+
+       l_j2k->cstr_index->marker = (opj_marker_info_t*) opj_malloc(100 * sizeof(opj_marker_info_t));
+*/
+
+       /* validation list creation */
+       l_j2k->m_validation_list = opj_procedure_list_create();
+       if (! l_j2k->m_validation_list) {
+               j2k_destroy(l_j2k);
+               return 00;
+       }
+
+       /* execution list creation */
+       l_j2k->m_procedure_list = opj_procedure_list_create();
+       if (! l_j2k->m_procedure_list) {
+               j2k_destroy(l_j2k);
+               return 00;
+       }
+
+       return l_j2k;
+}
+
+
+opj_codestream_index_t* j2k_create_cstr_index(void)
+{
+       opj_codestream_index_t* cstr_index = (opj_codestream_index_t*)
+                       opj_calloc(1,sizeof(opj_codestream_index_t));
+       if (!cstr_index)
+               return NULL;
+
+       cstr_index->maxmarknum = 100;
+       cstr_index->marknum = 0;
+       cstr_index->marker = (opj_marker_info_t*)
+                       opj_calloc(cstr_index->maxmarknum, sizeof(opj_marker_info_t));
+       if (!cstr_index-> marker)
+               return NULL;
+
+       cstr_index->tile_index = NULL;
+
+       return cstr_index;
+}
+
+/**
+ * Reads a SPCod or SPCoc element, i.e. the coding style of a given component of a tile.
+ * @param      p_header_data   the data contained in the COM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the COM marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_SPCod_SPCoc(
+                                                   opj_j2k_v2_t *p_j2k,
+                                                       OPJ_UINT32 compno,
+                                                       OPJ_BYTE * p_header_data,
+                                                       OPJ_UINT32 * p_header_size,
+                                                       struct opj_event_mgr * p_manager
+                                                       )
+{
+       OPJ_UINT32 i, l_tmp;
+       opj_cp_v2_t *l_cp = NULL;
+       opj_tcp_v2_t *l_tcp = NULL;
+       opj_tccp_t *l_tccp = NULL;
+       OPJ_BYTE * l_current_ptr = NULL;
+
+       /* preconditions */
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       assert(p_header_data != 00);
+
+       l_cp = &(p_j2k->m_cp);
+       l_tcp = (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH) ?
+                               &l_cp->tcps[p_j2k->m_current_tile_number] :
+                               p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       /* precondition again */
+       assert(compno < p_j2k->m_image->numcomps);
+
+       l_tccp = &l_tcp->tccps[compno];
+       l_current_ptr = p_header_data;
+
+       /* make sure room is sufficient */
+       if (*p_header_size < 5) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading SPCod SPCoc element\n");
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(l_current_ptr, &l_tccp->numresolutions ,1);              /* SPcox (D) */
+       ++l_tccp->numresolutions;                                                                               /* tccp->numresolutions = read() + 1 */
+       ++l_current_ptr;
+
+       /* If user wants to remove more resolutions than the codestream contains, return error */
+       if (l_cp->m_specific_param.m_dec.m_reduce >= l_tccp->numresolutions) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error decoding component %d.\nThe number of resolutions to remove is higher than the number "
+                                       "of resolutions of this component\nModify the cp_reduce parameter.\n\n", compno);
+               p_j2k->m_specific_param.m_decoder.m_state |= 0x8000;// FIXME J2K_DEC_STATE_ERR;
+               return OPJ_FALSE;
+       }
+
+       opj_read_bytes(l_current_ptr,&l_tccp->cblkw ,1);                /* SPcoc (E) */
+       ++l_current_ptr;
+       l_tccp->cblkw += 2;
+
+       opj_read_bytes(l_current_ptr,&l_tccp->cblkh ,1);                /* SPcoc (F) */
+       ++l_current_ptr;
+       l_tccp->cblkh += 2;
+
+       opj_read_bytes(l_current_ptr,&l_tccp->cblksty ,1);              /* SPcoc (G) */
+       ++l_current_ptr;
+
+       opj_read_bytes(l_current_ptr,&l_tccp->qmfbid ,1);               /* SPcoc (H) */
+       ++l_current_ptr;
+
+       *p_header_size = *p_header_size - 5;
+
+       /* use custom precinct size ? */
+       if (l_tccp->csty & J2K_CCP_CSTY_PRT) {
+               if (*p_header_size < l_tccp->numresolutions) {
+                       opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading SPCod SPCoc element\n");
+                       return OPJ_FALSE;
+               }
+
+               for     (i = 0; i < l_tccp->numresolutions; ++i) {
+                       opj_read_bytes(l_current_ptr,&l_tmp ,1);                /* SPcoc (I_i) */
+                       ++l_current_ptr;
+                       l_tccp->prcw[i] = l_tmp & 0xf;
+                       l_tccp->prch[i] = l_tmp >> 4;
+               }
+
+               *p_header_size = *p_header_size - l_tccp->numresolutions;
+       }
+       else {
+               /* set default size for the precinct width and height */
+               for     (i = 0; i < l_tccp->numresolutions; ++i) {
+                       l_tccp->prcw[i] = 15;
+                       l_tccp->prch[i] = 15;
+               }
+       }
+
+#ifdef WIP_REMOVE_MSD
+       /* INDEX >> */
+       if (p_j2k->cstr_info && compno == 0) {
+               OPJ_UINT32 l_data_size = l_tccp->numresolutions * sizeof(OPJ_UINT32);
+
+               p_j2k->cstr_info->tile[p_j2k->m_current_tile_number].tccp_info[compno].cblkh = l_tccp->cblkh;
+               p_j2k->cstr_info->tile[p_j2k->m_current_tile_number].tccp_info[compno].cblkw = l_tccp->cblkw;
+               p_j2k->cstr_info->tile[p_j2k->m_current_tile_number].tccp_info[compno].numresolutions = l_tccp->numresolutions;
+               p_j2k->cstr_info->tile[p_j2k->m_current_tile_number].tccp_info[compno].cblksty = l_tccp->cblksty;
+               p_j2k->cstr_info->tile[p_j2k->m_current_tile_number].tccp_info[compno].qmfbid = l_tccp->qmfbid;
+
+
+               memcpy(p_j2k->cstr_info->tile[p_j2k->m_current_tile_number].pdx,l_tccp->prcw, l_data_size);
+               memcpy(p_j2k->cstr_info->tile[p_j2k->m_current_tile_number].pdy,l_tccp->prch, l_data_size);
+       }
+       /* << INDEX */
+#endif
+
+       return OPJ_TRUE;
+}
+
+/**
+ * Copies the tile component parameters of all the component from the first tile component.
+ *
+ * @param              p_j2k           the J2k codec.
+ */
+void j2k_copy_tile_component_parameters( opj_j2k_v2_t *p_j2k )
+{
+       // loop
+       OPJ_UINT32 i;
+       opj_cp_v2_t *l_cp = NULL;
+       opj_tcp_v2_t *l_tcp = NULL;
+       opj_tccp_t *l_ref_tccp = NULL, *l_copied_tccp = NULL;
+       OPJ_UINT32 l_prc_size;
+
+       /* preconditions */
+       assert(p_j2k != 00);
+
+       l_cp = &(p_j2k->m_cp);
+       l_tcp = (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH) ? /* FIXME J2K_DEC_STATE_TPH*/
+                               &l_cp->tcps[p_j2k->m_current_tile_number] :
+                               p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       l_ref_tccp = &l_tcp->tccps[0];
+       l_copied_tccp = l_ref_tccp + 1;
+       l_prc_size = l_ref_tccp->numresolutions * sizeof(OPJ_UINT32);
+
+       for     (i=1; i<p_j2k->m_image->numcomps; ++i) {
+               l_copied_tccp->numresolutions = l_ref_tccp->numresolutions;
+               l_copied_tccp->cblkw = l_ref_tccp->cblkw;
+               l_copied_tccp->cblkh = l_ref_tccp->cblkh;
+               l_copied_tccp->cblksty = l_ref_tccp->cblksty;
+               l_copied_tccp->qmfbid = l_ref_tccp->qmfbid;
+               memcpy(l_copied_tccp->prcw,l_ref_tccp->prcw,l_prc_size);
+               memcpy(l_copied_tccp->prch,l_ref_tccp->prch,l_prc_size);
+               ++l_copied_tccp;
+       }
+}
+
+/**
+ * Reads a SQcd or SQcc element, i.e. the quantization values of a band.
+ *
+ * @param      p_comp_no               the component being targeted.
+ * @param      p_header_data   the data contained in the COM box.
+ * @param      p_j2k                   the jpeg2000 codec.
+ * @param      p_header_size   the size of the data contained in the COM marker.
+ * @param      p_manager               the user event manager.
+*/
+opj_bool j2k_read_SQcd_SQcc(
+                                                       opj_j2k_v2_t *p_j2k,
+                                                       OPJ_UINT32 p_comp_no,
+                                                       OPJ_BYTE* p_header_data,
+                                                       OPJ_UINT32 * p_header_size,
+                                                       struct opj_event_mgr * p_manager
+                                                       )
+{
+       // loop
+       OPJ_UINT32 l_band_no;
+       opj_cp_v2_t *l_cp = 00;
+       opj_tcp_v2_t *l_tcp = 00;
+       opj_tccp_t *l_tccp = 00;
+       OPJ_BYTE * l_current_ptr = 00;
+       OPJ_UINT32 l_tmp, l_num_band;
+
+       // preconditions
+       assert(p_j2k != 00);
+       assert(p_manager != 00);
+       assert(p_header_data != 00);
+
+       l_cp = &(p_j2k->m_cp);
+       // come from tile part header or main header ?
+       l_tcp = (p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH) ? /*FIXME J2K_DEC_STATE_TPH*/
+                               &l_cp->tcps[p_j2k->m_current_tile_number] :
+                               p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       // precondition again
+       assert(p_comp_no <  p_j2k->m_image->numcomps);
+
+       l_tccp = &l_tcp->tccps[p_comp_no];
+       l_current_ptr = p_header_data;
+
+       if (*p_header_size < 1) {
+               opj_event_msg_v2(p_manager, EVT_ERROR, "Error reading SQcd or SQcc element\n");
+               return OPJ_FALSE;
+       }
+       *p_header_size -= 1;
+
+       opj_read_bytes(l_current_ptr, &l_tmp ,1);                       /* Sqcx */
+       ++l_current_ptr;
+
+       l_tccp->qntsty = l_tmp & 0x1f;
+       l_tccp->numgbits = l_tmp >> 5;
+       if (l_tccp->qntsty == J2K_CCP_QNTSTY_SIQNT) {
+        l_num_band = 1;
+       }
+       else {
+               l_num_band = (l_tccp->qntsty == J2K_CCP_QNTSTY_NOQNT) ?
+                       (*p_header_size) :
+                       (*p_header_size) / 2;
+
+               if( l_num_band > J2K_MAXBANDS ) {
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "While reading CCP_QNTSTY element inside QCD or QCC marker segment, "
+                               "number of subbands (%d) is greater to J2K_MAXBANDS (%d). So we limiting the number of elements stored to "
+                               "J2K_MAXBANDS (%d) and skip the other. \n", l_num_band, J2K_MAXBANDS, J2K_MAXBANDS);
+                       //return OPJ_FALSE;
+               }
+       }
+
+#ifdef USE_JPWL
+       if (l_cp->correct) {
+
+               /* if JPWL is on, we check whether there are too many subbands */
+               if ((l_num_band < 0) || (l_num_band >= J2K_MAXBANDS)) {
+                       opj_event_msg_v2(p_manager, JPWL_ASSUME ? EVT_WARNING : EVT_ERROR,
+                               "JPWL: bad number of subbands in Sqcx (%d)\n",
+                               l_num_band);
+                       if (!JPWL_ASSUME) {
+                               opj_event_msg_v2(p_manager, EVT_ERROR, "JPWL: giving up\n");
+                               return OPJ_FALSE;
+                       }
+                       /* we try to correct */
+                       l_num_band = 1;
+                       opj_event_msg_v2(p_manager, EVT_WARNING, "- trying to adjust them\n"
+                               "- setting number of bands to %d => HYPOTHESIS!!!\n",
+                               l_num_band);
+               };
+
+       };
+#endif /* USE_JPWL */
+
+       if (l_tccp->qntsty == J2K_CCP_QNTSTY_NOQNT) {
+               for     (l_band_no = 0; l_band_no < l_num_band; l_band_no++) {
+                       opj_read_bytes(l_current_ptr, &l_tmp ,1);                       /* SPqcx_i */
+                       ++l_current_ptr;
+                       if (l_band_no < J2K_MAXBANDS){
+                               l_tccp->stepsizes[l_band_no].expn = l_tmp>>3;
+                               l_tccp->stepsizes[l_band_no].mant = 0;
+                       }
+               }
+               *p_header_size = *p_header_size - l_num_band;
+       }
+       else {
+               for     (l_band_no = 0; l_band_no < l_num_band; l_band_no++) {
+                       opj_read_bytes(l_current_ptr, &l_tmp ,2);                       /* SPqcx_i */
+                       l_current_ptr+=2;
+                       if (l_band_no < J2K_MAXBANDS){
+                               l_tccp->stepsizes[l_band_no].expn = l_tmp >> 11;
+                               l_tccp->stepsizes[l_band_no].mant = l_tmp & 0x7ff;
+                       }
+               }
+               *p_header_size = *p_header_size - 2*l_num_band;
+       }
+
+       /* Add Antonin : if scalar_derived -> compute other stepsizes */
+       if (l_tccp->qntsty == J2K_CCP_QNTSTY_SIQNT) {
+               for (l_band_no = 1; l_band_no < J2K_MAXBANDS; l_band_no++) {
+                       l_tccp->stepsizes[l_band_no].expn =
+                               ((l_tccp->stepsizes[0].expn) - ((l_band_no - 1) / 3) > 0) ?
+                                       (l_tccp->stepsizes[0].expn) - ((l_band_no - 1) / 3) : 0;
+                       l_tccp->stepsizes[l_band_no].mant = l_tccp->stepsizes[0].mant;
+               }
+       }
+
+       return OPJ_TRUE;
+}
+
+/**
+ * Copies the tile component parameters of all the component from the first tile component.
+ *
+ * @param              p_j2k           the J2k codec.
+ */
+void j2k_copy_tile_quantization_parameters( opj_j2k_v2_t *p_j2k )
+{
+       OPJ_UINT32 i;
+       opj_cp_v2_t *l_cp = NULL;
+       opj_tcp_v2_t *l_tcp = NULL;
+       opj_tccp_t *l_ref_tccp = NULL;
+       opj_tccp_t *l_copied_tccp = NULL;
+       OPJ_UINT32 l_size;
+
+       /* preconditions */
+       assert(p_j2k != 00);
+
+       l_cp = &(p_j2k->m_cp);
+       l_tcp = p_j2k->m_specific_param.m_decoder.m_state == J2K_STATE_TPH ?
+                       &l_cp->tcps[p_j2k->m_current_tile_number] :
+                       p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       l_ref_tccp = &l_tcp->tccps[0];
+       l_copied_tccp = l_ref_tccp + 1;
+       l_size = J2K_MAXBANDS * sizeof(opj_stepsize_t);
+
+       for     (i=1;i<p_j2k->m_image->numcomps;++i) {
+               l_copied_tccp->qntsty = l_ref_tccp->qntsty;
+               l_copied_tccp->numgbits = l_ref_tccp->numgbits;
+               memcpy(l_copied_tccp->stepsizes,l_ref_tccp->stepsizes,l_size);
+               ++l_copied_tccp;
+       }
+}
+
+/**
+ * Dump some elements from the J2K decompression structure .
+ *
+ *@param p_j2k                         the jpeg2000 codec.
+ *@param flag                          flag to describe what elments are dump.
+ *@param out_stream                    output stream where dump the elements.
+ *
+*/
+void j2k_dump (opj_j2k_v2_t* p_j2k, OPJ_INT32 flag, FILE* out_stream)
+{
+       /* Check if the flag is compatible with j2k file*/
+       if ( (flag & OPJ_JP2_INFO) || (flag & OPJ_JP2_IND)){
+               fprintf(out_stream, "Wrong flag\n");
+               return;
+       }
+
+       /* Dump the image_header */
+       if (flag & OPJ_IMG_INFO){
+               if (p_j2k->m_image)
+                       j2k_dump_image_header(p_j2k->m_image, 0, out_stream);
+       }
+
+       /* Dump the codestream info from main header */
+       if (flag & OPJ_J2K_MH_INFO){
+               j2k_dump_MH_info(p_j2k, out_stream);
+       }
+
+
+       /* Dump the codestream info of the current tile */
+       if (flag & OPJ_J2K_TH_INFO){
+
+       }
+
+       /* Dump the codestream index from main header */
+       if (flag & OPJ_J2K_MH_IND){
+               j2k_dump_MH_index(p_j2k, out_stream);
+       }
+
+       /* Dump the codestream index of the current tile */
+       if (flag & OPJ_J2K_TH_IND){
+
+       }
+
+}
+
+/**
+ * Dump index elements of the codestream extract from the main header.
+ *
+ *@param p_j2k                         the jpeg2000 codec.
+ *@param out_stream                    output stream where dump the elements.
+ *
+*/
+void j2k_dump_MH_index(opj_j2k_v2_t* p_j2k, FILE* out_stream)
+{
+       opj_codestream_index_t* cstr_index = p_j2k->cstr_index;
+       OPJ_UINT32 it_marker;
+
+       fprintf(out_stream, "Codestream index from main header: {\n");
+
+       fprintf(out_stream, "\t Main header start position=%d\n\t Main header end position=%d\n",
+                       cstr_index->main_head_start, cstr_index->main_head_end);
+
+       fprintf(out_stream, "\t Marker list: {\n");
+
+       for (it_marker=0; it_marker < cstr_index->marknum ; it_marker++){
+               fprintf(out_stream, "\t\t type=%#x, pos=%d, len=%d\n",
+                               cstr_index->marker[it_marker].type,
+                               cstr_index->marker[it_marker].pos,
+                               cstr_index->marker[it_marker].len );
+       }
+
+       fprintf(out_stream, "\t }\n}\n");
+
+}
+
+/**
+ * Dump info elements of the codestream extract from the main header.
+ *
+ *@param p_j2k                         the jpeg2000 codec.
+ *@param out_stream                    output stream where dump the elements.
+ *
+*/
+void j2k_dump_MH_info(opj_j2k_v2_t* p_j2k, FILE* out_stream)
+{
+       opj_tcp_v2_t * l_default_tile=NULL;
+
+       fprintf(out_stream, "Codestream info from main header: {\n");
+
+       fprintf(out_stream, "\t tx0=%d, ty0=%d\n", p_j2k->m_cp.tx0, p_j2k->m_cp.ty0);
+       fprintf(out_stream, "\t tdx=%d, tdy=%d\n", p_j2k->m_cp.tdx, p_j2k->m_cp.tdy);
+       fprintf(out_stream, "\t tw=%d, th=%d\n", p_j2k->m_cp.tw, p_j2k->m_cp.th);
+
+       l_default_tile = p_j2k->m_specific_param.m_decoder.m_default_tcp;
+       if (l_default_tile)
+       {
+               OPJ_INT32 compno;
+               OPJ_INT32 numcomps = p_j2k->m_image->numcomps;
+
+               fprintf(out_stream, "\t default tile {\n");
+               fprintf(out_stream, "\t\t csty=%#x\n", l_default_tile->csty);
+               fprintf(out_stream, "\t\t prg=%#x\n", l_default_tile->prg);
+               fprintf(out_stream, "\t\t numlayers=%d\n", l_default_tile->numlayers);
+               fprintf(out_stream, "\t\t mct=%x\n", l_default_tile->mct);
+
+               for (compno = 0; compno < numcomps; compno++) {
+                       opj_tccp_t *l_tccp = &(l_default_tile->tccps[compno]);
+                       OPJ_INT32 resno, bandno, numbands;
+
+                       /* coding style*/
+                       fprintf(out_stream, "\t\t comp %d {\n", compno);
+                       fprintf(out_stream, "\t\t\t csty=%#x\n", l_tccp->csty);
+                       fprintf(out_stream, "\t\t\t numresolutions=%d\n", l_tccp->numresolutions);
+                       fprintf(out_stream, "\t\t\t cblkw=2^%d\n", l_tccp->cblkw);
+                       fprintf(out_stream, "\t\t\t cblkh=2^%d\n", l_tccp->cblkh);
+                       fprintf(out_stream, "\t\t\t cblksty=%#x\n", l_tccp->cblksty);
+                       fprintf(out_stream, "\t\t\t qmfbid=%d\n", l_tccp->qmfbid);
+
+                       fprintf(out_stream, "\t\t\t preccintsize (w,h)=");
+                       for (resno = 0; resno < l_tccp->numresolutions; resno++) {
+                               fprintf(out_stream, "(%d,%d) ", l_tccp->prcw[resno], l_tccp->prch[resno]);
+                       }
+                       fprintf(out_stream, "\n");
+
+                       /* quantization style*/
+                       fprintf(out_stream, "\t\t\t qntsty=%d\n", l_tccp->qntsty);
+                       fprintf(out_stream, "\t\t\t numgbits=%d\n", l_tccp->numgbits);
+                       fprintf(out_stream, "\t\t\t stepsizes (m,e)=");
+                       numbands = (l_tccp->qntsty == J2K_CCP_QNTSTY_SIQNT) ? 1 : l_tccp->numresolutions * 3 - 2;
+                       for (bandno = 0; bandno < numbands; bandno++) {
+                               fprintf(out_stream, "(%d,%d) ", l_tccp->stepsizes[bandno].mant,
+                                       l_tccp->stepsizes[bandno].expn);
+                       }
+                       fprintf(out_stream, "\n");
+
+                       /* RGN value*/
+                       fprintf(out_stream, "\t\t\t roishift=%d\n", l_tccp->roishift);
+
+                       fprintf(out_stream, "\t\t }\n");
+               } /*end of component of default tile*/
+               fprintf(out_stream, "\t }\n"); /*end of default tile*/
+
+       }
+
+       fprintf(out_stream, "}\n");
+
+}
+
+/**
+ * Dump an image header structure.
+ *
+ *@param img_header                    the image header to dump.
+ *@param dev_dump_flag         flag to describe if we are in the case of this function is use outside j2k_dump function
+ *@param out_stream                    output stream where dump the elements.
+ */
+void j2k_dump_image_header(opj_image_t* img_header, opj_bool dev_dump_flag, FILE* out_stream)
+{
+       char tab[2];
+
+       if (dev_dump_flag){
+               fprintf(stdout, "[DEV] Dump a image_header struct {\n");
+               tab[0] = '\0';
+       }
+       else {
+               fprintf(out_stream, "Image info {\n");
+               tab[0] = '\t';tab[1] = '\0';
+       }
+
+       fprintf(out_stream, "%s x0=%d, y0=%d\n", tab, img_header->x0, img_header->y0);
+       fprintf(out_stream,     "%s x1=%d, y1=%d\n", tab, img_header->x1, img_header->y1);
+       fprintf(out_stream, "%s numcomps=%d\n", tab, img_header->numcomps);
+
+       if (img_header->comps){
+               int compno;
+               for (compno = 0; compno < img_header->numcomps; compno++) {
+                       fprintf(out_stream, "%s\t component %d {\n", tab, compno);
+                       j2k_dump_image_comp_header(&(img_header->comps[compno]), dev_dump_flag, out_stream);
+                       fprintf(out_stream,"%s}\n",tab);
+               }
+       }
+
+       fprintf(out_stream, "}\n");
+}
+
+/**
+ * Dump a component image header structure.
+ *
+ *@param comp_header           the component image header to dump.
+ *@param dev_dump_flag         flag to describe if we are in the case of this function is use outside j2k_dump function
+ *@param out_stream                    output stream where dump the elements.
+ */
+void j2k_dump_image_comp_header(opj_image_comp_t* comp_header, opj_bool dev_dump_flag, FILE* out_stream)
+{
+       char tab[3];
+
+       if (dev_dump_flag){
+               fprintf(stdout, "[DEV] Dump a image_comp_header struct {\n");
+               tab[0] = '\0';
+       }       else {
+               tab[0] = '\t';tab[1] = '\t';tab[2] = '\0';
+       }
+
+       fprintf(out_stream, "%s dx=%d, dy=%d\n", tab, comp_header->dx, comp_header->dy);
+       fprintf(out_stream, "%s prec=%d\n", tab, comp_header->prec);
+       fprintf(out_stream, "%s sgnd=%d\n", tab, comp_header->sgnd);
+
+       if (dev_dump_flag)
+               fprintf(out_stream, "}\n");
+}
+
+
+/**
+ * Get the codestream info from a JPEG2000 codec.
+ *
+ *@param       p_j2k                           the component image header to dump.
+ *
+ *@return      the codestream information extract from the jpg2000 codec
+ */
+opj_codestream_info_v2_t* j2k_get_cstr_info(opj_j2k_v2_t* p_j2k)
+{
+       OPJ_UINT16 compno;
+       OPJ_UINT16 numcomps = p_j2k->m_image->numcomps;
+       opj_tcp_v2_t *l_default_tile;
+       opj_codestream_info_v2_t* cstr_info = (opj_codestream_info_v2_t*) opj_calloc(1,sizeof(opj_codestream_info_v2_t));
+
+       cstr_info->nbcomps = p_j2k->m_image->numcomps;
+
+       cstr_info->tx0 = p_j2k->m_cp.tx0;
+       cstr_info->ty0 = p_j2k->m_cp.ty0;
+       cstr_info->tdx = p_j2k->m_cp.tdx;
+       cstr_info->tdy = p_j2k->m_cp.tdy;
+       cstr_info->tw = p_j2k->m_cp.tw;
+       cstr_info->th = p_j2k->m_cp.th;
+
+       l_default_tile = p_j2k->m_specific_param.m_decoder.m_default_tcp;
+
+       cstr_info->m_default_tile_info.csty = l_default_tile->csty;
+       cstr_info->m_default_tile_info.prg = l_default_tile->prg;
+       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));
+
+       for (compno = 0; compno < numcomps; compno++) {
+               opj_tccp_t *l_tccp = &(l_default_tile->tccps[compno]);
+               opj_tccp_info_t *l_tccp_info = &(cstr_info->m_default_tile_info.tccp_info[compno]);
+               OPJ_INT32 bandno, numbands;
+
+               /* coding style*/
+               l_tccp_info->csty = l_tccp->csty;
+               l_tccp_info->numresolutions = l_tccp->numresolutions;
+               l_tccp_info->cblkw = l_tccp->cblkw;
+               l_tccp_info->cblkh = l_tccp->cblkh;
+               l_tccp_info->cblksty = l_tccp->cblksty;
+               l_tccp_info->qmfbid = l_tccp->qmfbid;
+               if (l_tccp->numresolutions < J2K_MAXRLVLS)
+               {
+                       memcpy(l_tccp_info->prch, l_tccp->prch, l_tccp->numresolutions);
+                       memcpy(l_tccp_info->prcw, l_tccp->prcw, l_tccp->numresolutions);
+               }
+
+               /* quantization style*/
+               l_tccp_info->qntsty = l_tccp->qntsty;
+               l_tccp_info->numgbits = l_tccp->numgbits;
+
+               numbands = (l_tccp->qntsty == J2K_CCP_QNTSTY_SIQNT) ? 1 : l_tccp->numresolutions * 3 - 2;
+               if (numbands < J2K_MAXBANDS) {
+                       for (bandno = 0; bandno < numbands; bandno++) {
+                               l_tccp_info->stepsizes_mant[bandno] = l_tccp->stepsizes[bandno].mant;
+                               l_tccp_info->stepsizes_expn[bandno] = l_tccp->stepsizes[bandno].expn;
+                       }
+               }
+
+               /* RGN value*/
+               l_tccp_info->roishift = l_tccp->roishift;
+       }
+
+
+       return cstr_info;
+}
+
+/**
+ * Get the codestream index from a JPEG2000 codec.
+ *
+ *@param       p_j2k                           the component image header to dump.
+ *
+ *@return      the codestream index extract from the jpg2000 codec
+ */
+opj_codestream_index_t* j2k_get_cstr_index(opj_j2k_v2_t* p_j2k)
+{
+       opj_codestream_index_t* l_cstr_index = (opj_codestream_index_t*)
+                       opj_calloc(1,sizeof(opj_codestream_index_t));
+       if (!l_cstr_index)
+               return NULL;
+
+       l_cstr_index->main_head_start = p_j2k->cstr_index->main_head_start;
+       l_cstr_index->main_head_end = p_j2k->cstr_index->main_head_end;
+       l_cstr_index->codestream_size = p_j2k->cstr_index->codestream_size;
+
+       l_cstr_index->maxmarknum = p_j2k->cstr_index->maxmarknum;
+       l_cstr_index->marknum = p_j2k->cstr_index->marknum;
+       l_cstr_index->marker = (opj_marker_info_t*)opj_malloc(l_cstr_index->marknum*sizeof(opj_marker_info_t));
+       if (!l_cstr_index->marker)
+               return NULL;
+
+       memcpy(l_cstr_index->marker, p_j2k->cstr_index->marker, l_cstr_index->marknum * sizeof(opj_marker_info_t) );
+
+
+       return l_cstr_index;
+}
+
+
+
+/**
+ * Reads the tiles.
+ */
+opj_bool j2k_decode_tiles (    opj_j2k_v2_t *p_j2k,
+                                                       opj_stream_private_t *p_stream,
+                                                       opj_event_mgr_t * p_manager)
+{
+       opj_bool l_go_on = OPJ_TRUE;
+       OPJ_UINT32 l_current_tile_no;
+       OPJ_UINT32 l_data_size,l_max_data_size;
+       OPJ_INT32 l_tile_x0,l_tile_y0,l_tile_x1,l_tile_y1;
+       OPJ_UINT32 l_nb_comps;
+       OPJ_BYTE * l_current_data;
+
+       l_current_data = (OPJ_BYTE*)opj_malloc(1000);
+       if (! l_current_data) {
+               return OPJ_FALSE;
+       }
+       l_max_data_size = 1000;
+
+       while (OPJ_TRUE) {
+               if (! j2k_read_tile_header(     p_j2k,
+                                                                       &l_current_tile_no,
+                                                                       &l_data_size,
+                                                                       &l_tile_x0, &l_tile_y0,
+                                                                       &l_tile_x1, &l_tile_y1,
+                                                                       &l_nb_comps,
+                                                                       &l_go_on,
+                                                                       p_stream,
+                                                                       p_manager)) {
+                       return OPJ_FALSE;
+               }
+
+               if (! l_go_on) {
+                       break;
+               }
+
+               if (l_data_size > l_max_data_size) {
+                       l_current_data = (OPJ_BYTE*)opj_realloc(l_current_data,l_data_size);
+                       if (! l_current_data) {
+                               return OPJ_FALSE;
+                       }
+
+                       l_max_data_size = l_data_size;
+               }
+
+               if (! j2k_decode_tile(p_j2k,l_current_tile_no,l_current_data,l_data_size,p_stream,p_manager)) {
+                       opj_free(l_current_data);
+                       return OPJ_FALSE;
+               }
+               opj_event_msg_v2(p_manager, EVT_INFO, "Tile %d/%d has been decode.\n", l_current_tile_no +1, p_j2k->m_cp.th * p_j2k->m_cp.tw);
+
+               if (! j2k_update_image_data(p_j2k->m_tcd,l_current_data)) {
+                       opj_free(l_current_data);
+                       return OPJ_FALSE;
+               }
+               opj_event_msg_v2(p_manager, EVT_INFO, "Image data has been updated with tile %d.\n\n", l_current_tile_no + 1);
+
+       }
+
+       opj_free(l_current_data);
+       return OPJ_TRUE;
+}
+
+
+/**
+ * Sets up the procedures to do on decoding data. Developpers wanting to extend the library can add their own reading procedures.
+ */
+void j2k_setup_decoding (opj_j2k_v2_t *p_j2k)
+{
+       // preconditions
+       assert(p_j2k != 00);
+
+       opj_procedure_list_add_procedure(p_j2k->m_procedure_list,(void*)j2k_decode_tiles);
+       /* DEVELOPER CORNER, add your custom procedures */
+
+}
+
+
+/**
+ * Decodes the tiles of the stream.
+ */
+opj_bool j2k_decode_v2(        opj_j2k_v2_t * p_j2k,
+                                               opj_stream_private_t * p_stream,
+                                               opj_image_t * p_image,
+                                               opj_event_mgr_t * p_manager)
+{
+       OPJ_UINT32 compno;
+
+       /* customization of the decoding */
+       j2k_setup_decoding(p_j2k);
+
+       /* write header */
+       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;
+       }
+
+       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*/;
+}