2003-11-27 11:10:17 +01:00
|
|
|
/*
|
2017-05-15 12:21:30 +02:00
|
|
|
* The copyright in this software is being made available under the 2-clauses
|
|
|
|
* BSD License, included below. This software may be subject to other third
|
2014-04-03 17:30:57 +02:00
|
|
|
* party and contributor rights, including patent rights, and no such rights
|
|
|
|
* are granted under this license.
|
|
|
|
*
|
|
|
|
* Copyright (c) 2002-2014, Universite catholique de Louvain (UCL), Belgium
|
|
|
|
* Copyright (c) 2002-2014, Professor Benoit Macq
|
2005-12-02 14:34:15 +01:00
|
|
|
* Copyright (c) 2001-2003, David Janssens
|
|
|
|
* Copyright (c) 2002-2003, Yannick Verschueren
|
2017-05-15 12:21:30 +02:00
|
|
|
* Copyright (c) 2003-2007, Francois-Olivier Devaux
|
2014-04-03 17:30:57 +02:00
|
|
|
* Copyright (c) 2003-2014, Antonin Descampe
|
2007-01-15 10:55:40 +01:00
|
|
|
* Copyright (c) 2005, Herve Drolon, FreeImage Team
|
2003-11-27 11:10:17 +01:00
|
|
|
* All rights reserved.
|
|
|
|
*
|
|
|
|
* Redistribution and use in source and binary forms, with or without
|
|
|
|
* modification, are permitted provided that the following conditions
|
|
|
|
* are met:
|
|
|
|
* 1. Redistributions of source code must retain the above copyright
|
|
|
|
* notice, this list of conditions and the following disclaimer.
|
|
|
|
* 2. Redistributions in binary form must reproduce the above copyright
|
|
|
|
* notice, this list of conditions and the following disclaimer in the
|
|
|
|
* documentation and/or other materials provided with the distribution.
|
|
|
|
*
|
|
|
|
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS `AS IS'
|
|
|
|
* AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
|
|
|
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
|
|
|
* ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
|
|
|
|
* LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
|
|
|
|
* CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
|
|
|
|
* SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
|
|
|
|
* INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
|
|
|
|
* CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
|
|
|
|
* ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
|
|
|
|
* POSSIBILITY OF SUCH DAMAGE.
|
|
|
|
*/
|
|
|
|
|
2017-07-30 19:46:52 +02:00
|
|
|
#ifndef OPJ_PI_H
|
|
|
|
#define OPJ_PI_H
|
2005-12-02 14:34:15 +01:00
|
|
|
/**
|
|
|
|
@file pi.h
|
|
|
|
@brief Implementation of a packet iterator (PI)
|
2003-11-27 11:10:17 +01:00
|
|
|
|
2005-12-02 14:34:15 +01:00
|
|
|
The functions in PI.C have for goal to realize a packet iterator that permits to get the next
|
|
|
|
packet following the progression order and change of it. The functions in PI.C are used
|
|
|
|
by some function in T2.C.
|
|
|
|
*/
|
2003-11-27 11:10:17 +01:00
|
|
|
|
2005-12-02 14:34:15 +01:00
|
|
|
/** @defgroup PI PI - Implementation of a packet iterator */
|
|
|
|
/*@{*/
|
|
|
|
|
|
|
|
/**
|
2012-10-03 14:39:17 +02:00
|
|
|
FIXME DOC
|
2005-12-02 14:34:15 +01:00
|
|
|
*/
|
|
|
|
typedef struct opj_pi_resolution {
|
2017-05-15 12:21:30 +02:00
|
|
|
OPJ_UINT32 pdx, pdy;
|
|
|
|
OPJ_UINT32 pw, ph;
|
2005-12-02 14:34:15 +01:00
|
|
|
} opj_pi_resolution_t;
|
2003-11-27 11:10:17 +01:00
|
|
|
|
2005-12-02 14:34:15 +01:00
|
|
|
/**
|
2012-10-03 14:39:17 +02:00
|
|
|
FIXME DOC
|
2005-12-02 14:34:15 +01:00
|
|
|
*/
|
|
|
|
typedef struct opj_pi_comp {
|
2017-05-15 12:21:30 +02:00
|
|
|
OPJ_UINT32 dx, dy;
|
|
|
|
/** number of resolution levels */
|
|
|
|
OPJ_UINT32 numresolutions;
|
|
|
|
opj_pi_resolution_t *resolutions;
|
2005-12-02 14:34:15 +01:00
|
|
|
} opj_pi_comp_t;
|
2003-11-27 11:10:17 +01:00
|
|
|
|
2012-08-30 18:56:31 +02:00
|
|
|
/**
|
|
|
|
Packet iterator
|
2005-12-02 14:34:15 +01:00
|
|
|
*/
|
|
|
|
typedef struct opj_pi_iterator {
|
2017-05-15 12:21:30 +02:00
|
|
|
/** Enabling Tile part generation*/
|
|
|
|
OPJ_BYTE tp_on;
|
|
|
|
/** precise if the packet has been already used (useful for progression order change) */
|
|
|
|
OPJ_INT16 *include;
|
2017-07-26 12:24:26 +02:00
|
|
|
/** Number of elements in include array */
|
|
|
|
OPJ_UINT32 include_size;
|
2017-05-15 12:21:30 +02:00
|
|
|
/** layer step used to localize the packet in the include vector */
|
|
|
|
OPJ_UINT32 step_l;
|
|
|
|
/** resolution step used to localize the packet in the include vector */
|
|
|
|
OPJ_UINT32 step_r;
|
|
|
|
/** component step used to localize the packet in the include vector */
|
|
|
|
OPJ_UINT32 step_c;
|
|
|
|
/** precinct step used to localize the packet in the include vector */
|
|
|
|
OPJ_UINT32 step_p;
|
|
|
|
/** component that identify the packet */
|
|
|
|
OPJ_UINT32 compno;
|
|
|
|
/** resolution that identify the packet */
|
|
|
|
OPJ_UINT32 resno;
|
|
|
|
/** precinct that identify the packet */
|
|
|
|
OPJ_UINT32 precno;
|
|
|
|
/** layer that identify the packet */
|
|
|
|
OPJ_UINT32 layno;
|
|
|
|
/** 0 if the first packet */
|
|
|
|
OPJ_BOOL first;
|
|
|
|
/** progression order change information */
|
|
|
|
opj_poc_t poc;
|
|
|
|
/** number of components in the image */
|
|
|
|
OPJ_UINT32 numcomps;
|
|
|
|
/** Components*/
|
|
|
|
opj_pi_comp_t *comps;
|
|
|
|
/** FIXME DOC*/
|
2019-11-17 01:18:26 +01:00
|
|
|
OPJ_UINT32 tx0, ty0, tx1, ty1;
|
2017-05-15 12:21:30 +02:00
|
|
|
/** FIXME DOC*/
|
2019-11-17 01:18:26 +01:00
|
|
|
OPJ_UINT32 x, y;
|
2017-05-15 12:21:30 +02:00
|
|
|
/** FIXME DOC*/
|
|
|
|
OPJ_UINT32 dx, dy;
|
2020-12-02 14:02:17 +01:00
|
|
|
/** event manager */
|
|
|
|
opj_event_mgr_t* manager;
|
2005-12-02 14:34:15 +01:00
|
|
|
} opj_pi_iterator_t;
|
2003-11-27 11:10:17 +01:00
|
|
|
|
2005-12-02 14:34:15 +01:00
|
|
|
/** @name Exported functions */
|
|
|
|
/*@{*/
|
|
|
|
/* ----------------------------------------------------------------------- */
|
2012-03-15 11:23:20 +01:00
|
|
|
/**
|
|
|
|
* Creates a packet iterator for encoding.
|
|
|
|
*
|
2017-05-15 12:21:30 +02:00
|
|
|
* @param image the image being encoded.
|
|
|
|
* @param cp the coding parameters.
|
|
|
|
* @param tileno index of the tile being encoded.
|
|
|
|
* @param t2_mode the type of pass for generating the packet iterator
|
2020-12-02 14:02:17 +01:00
|
|
|
* @param manager Event manager
|
2012-08-30 18:56:31 +02:00
|
|
|
*
|
2017-05-15 12:21:30 +02:00
|
|
|
* @return a list of packet iterator that points to the first packet of the tile (not true).
|
2012-03-15 11:23:20 +01:00
|
|
|
*/
|
2012-10-03 16:29:26 +02:00
|
|
|
opj_pi_iterator_t *opj_pi_initialise_encode(const opj_image_t *image,
|
2017-05-15 12:21:30 +02:00
|
|
|
opj_cp_t *cp,
|
|
|
|
OPJ_UINT32 tileno,
|
2020-12-02 14:02:17 +01:00
|
|
|
J2K_T2_MODE t2_mode,
|
|
|
|
opj_event_mgr_t* manager);
|
2012-03-15 11:23:20 +01:00
|
|
|
|
|
|
|
/**
|
|
|
|
* Updates the encoding parameters of the codec.
|
|
|
|
*
|
2017-05-15 12:21:30 +02:00
|
|
|
* @param p_image the image being encoded.
|
|
|
|
* @param p_cp the coding parameters.
|
|
|
|
* @param p_tile_no index of the tile being encoded.
|
2012-03-15 11:23:20 +01:00
|
|
|
*/
|
2017-05-15 12:21:30 +02:00
|
|
|
void opj_pi_update_encoding_parameters(const opj_image_t *p_image,
|
|
|
|
opj_cp_t *p_cp,
|
|
|
|
OPJ_UINT32 p_tile_no);
|
2012-03-15 11:23:20 +01:00
|
|
|
|
|
|
|
/**
|
|
|
|
Modify the packet iterator for enabling tile part generation
|
|
|
|
@param pi Handle to the packet iterator generated in pi_initialise_encode
|
|
|
|
@param cp Coding parameters
|
|
|
|
@param tileno Number that identifies the tile for which to list the packets
|
2012-08-30 18:56:31 +02:00
|
|
|
@param pino FIXME DOC
|
2012-03-15 11:23:20 +01:00
|
|
|
@param tpnum Tile part number of the current tile
|
|
|
|
@param tppos The position of the tile part flag in the progression order
|
2012-08-30 18:56:31 +02:00
|
|
|
@param t2_mode FIXME DOC
|
2012-03-15 11:23:20 +01:00
|
|
|
*/
|
2017-05-15 12:21:30 +02:00
|
|
|
void opj_pi_create_encode(opj_pi_iterator_t *pi,
|
|
|
|
opj_cp_t *cp,
|
|
|
|
OPJ_UINT32 tileno,
|
|
|
|
OPJ_UINT32 pino,
|
|
|
|
OPJ_UINT32 tpnum,
|
|
|
|
OPJ_INT32 tppos,
|
|
|
|
J2K_T2_MODE t2_mode);
|
2005-12-02 14:34:15 +01:00
|
|
|
|
2011-09-19 15:04:04 +02:00
|
|
|
/**
|
|
|
|
Create a packet iterator for Decoder
|
|
|
|
@param image Raw image for which the packets will be listed
|
|
|
|
@param cp Coding parameters
|
|
|
|
@param tileno Number that identifies the tile for which to list the packets
|
2020-12-02 14:02:17 +01:00
|
|
|
@param manager Event manager
|
2011-09-19 15:04:04 +02:00
|
|
|
@return Returns a packet iterator that points to the first packet of the tile
|
2012-10-03 16:29:26 +02:00
|
|
|
@see opj_pi_destroy
|
2011-09-19 15:04:04 +02:00
|
|
|
*/
|
2017-05-15 12:21:30 +02:00
|
|
|
opj_pi_iterator_t *opj_pi_create_decode(opj_image_t * image,
|
2012-10-25 17:09:31 +02:00
|
|
|
opj_cp_t * cp,
|
2020-12-02 14:02:17 +01:00
|
|
|
OPJ_UINT32 tileno,
|
|
|
|
opj_event_mgr_t* manager);
|
2011-09-19 15:04:04 +02:00
|
|
|
/**
|
|
|
|
* Destroys a packet iterator array.
|
|
|
|
*
|
2017-05-15 12:21:30 +02:00
|
|
|
* @param p_pi the packet iterator array to destroy.
|
|
|
|
* @param p_nb_elements the number of elements in the array.
|
2011-09-19 15:04:04 +02:00
|
|
|
*/
|
2012-10-03 16:29:26 +02:00
|
|
|
void opj_pi_destroy(opj_pi_iterator_t *p_pi,
|
|
|
|
OPJ_UINT32 p_nb_elements);
|
2011-09-19 15:04:04 +02:00
|
|
|
|
2005-12-02 14:34:15 +01:00
|
|
|
/**
|
|
|
|
Modify the packet iterator to point to the next packet
|
|
|
|
@param pi Packet iterator to modify
|
2012-08-30 18:56:31 +02:00
|
|
|
@return Returns false if pi pointed to the last packet or else returns true
|
2005-12-02 14:34:15 +01:00
|
|
|
*/
|
2012-11-15 14:13:36 +01:00
|
|
|
OPJ_BOOL opj_pi_next(opj_pi_iterator_t * pi);
|
Add support for generation of PLT markers in encoder
* -PLT switch added to opj_compress
* Add a opj_encoder_set_extra_options() function that
accepts a PLT=YES option, and could be expanded later
for other uses.
-------
Testing with a Sentinel2 10m band, T36JTT_20160914T074612_B02.jp2,
coming from S2A_MSIL1C_20160914T074612_N0204_R135_T36JTT_20160914T081456.SAFE
Decompress it to TIFF:
```
opj_uncompress -i T36JTT_20160914T074612_B02.jp2 -o T36JTT_20160914T074612_B02.tif
```
Recompress it with similar parameters as original:
```
opj_compress -n 5 -c [256,256],[256,256],[256,256],[256,256],[256,256] -t 1024,1024 -PLT -i T36JTT_20160914T074612_B02.tif -o T36JTT_20160914T074612_B02_PLT.jp2
```
Dump codestream detail with GDAL dump_jp2.py utility (https://github.com/OSGeo/gdal/blob/master/gdal/swig/python/samples/dump_jp2.py)
```
python dump_jp2.py T36JTT_20160914T074612_B02.jp2 > /tmp/dump_sentinel2_ori.txt
python dump_jp2.py T36JTT_20160914T074612_B02_PLT.jp2 > /tmp/dump_sentinel2_openjpeg_plt.txt
```
The diff between both show very similar structure, and identical number of packets in PLT markers
Now testing with Kakadu (KDU803_Demo_Apps_for_Linux-x86-64_200210)
Full file decompression:
```
kdu_expand -i T36JTT_20160914T074612_B02_PLT.jp2 -o tmp.tif
Consumed 121 tile-part(s) from a total of 121 tile(s).
Consumed 80,318,806 codestream bytes (excluding any file format) = 5.329697
bits/pel.
Processed using the multi-threaded environment, with
8 parallel threads of execution
```
Partial decompresson (presumably using PLT markers):
```
kdu_expand -i T36JTT_20160914T074612_B02.jp2 -o tmp.pgm -region "{0.5,0.5},{0.01,0.01}"
kdu_expand -i T36JTT_20160914T074612_B02_PLT.jp2 -o tmp2.pgm -region "{0.5,0.5},{0.01,0.01}"
diff tmp.pgm tmp2.pgm && echo "same !"
```
-------
Funded by ESA for S2-MPC project
2020-04-21 15:55:44 +02:00
|
|
|
|
|
|
|
/**
|
|
|
|
* Return the number of packets in the tile.
|
|
|
|
* @param image the image being encoded.
|
|
|
|
* @param cp Coding parameters
|
|
|
|
* @param tileno Number that identifies the tile.
|
|
|
|
*/
|
|
|
|
OPJ_UINT32 opj_get_encoding_packet_count(const opj_image_t *p_image,
|
|
|
|
const opj_cp_t *p_cp,
|
|
|
|
OPJ_UINT32 p_tile_no);
|
|
|
|
|
2005-12-02 14:34:15 +01:00
|
|
|
/* ----------------------------------------------------------------------- */
|
|
|
|
/*@}*/
|
|
|
|
|
|
|
|
/*@}*/
|
2003-11-27 11:10:17 +01:00
|
|
|
|
2017-07-30 19:46:52 +02:00
|
|
|
#endif /* OPJ_PI_H */
|