openjpeg/src/lib/openjp2/j2k.c, branch floating_point_irreversible_encoding Fork of openjpeg library for JPEG2000 encoding and decoding https://git.carlh.net/cgit/openjpeg/atom?h=floating_point_irreversible_encoding 2020-05-20T18:31:28Z opj_j2k_setup_encoder(): add validation of tile width and height to avoid potential division by zero 2020-05-20T18:31:28Z Even Rouault even.rouault@spatialys.com 2020-05-19T20:26:22Z urn:sha1:4ab2ed090747ad2f22e78d45a1db5640d40e6a34 Add support for generation of PLT markers in encoder 2020-04-21T13:55:44Z Even Rouault even.rouault@spatialys.com 2020-04-21T13:55:44Z urn:sha1:4edb8c83374f52cd6a8f2c7c875e8ffacccb5fa5 * -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 Rename mis-named function opj_tcd_get_encoded_tile_size() to opj_tcd_get_encoder_input_buffer_size() 2020-04-16T18:33:22Z Even Rouault even.rouault@spatialys.com 2020-04-16T18:33:22Z urn:sha1:221a801a97a3ea968a311f7905c18a1eb7f034c4 Implement writing of IMF profiles 2020-02-12T14:55:25Z Even Rouault even.rouault@spatialys.com 2020-02-12T14:55:16Z urn:sha1:84f3bebbff515f2b00ccf0c817930ebb10b91760 Add -IMF switch to opj_compress as well opj_j2k_update_image_dimensions(): reject images whose coordinates are beyond INT_MAX (fixes #1228) 2020-01-11T00:51:58Z Even Rouault even.rouault@spatialys.com 2020-01-11T00:51:19Z urn:sha1:024b8407392cb0b82b04b58ed256094ed5799e04 Merge pull request #1164 from sebras/master 2019-09-03T15:03:54Z Even Rouault even.rouault@spatialys.com 2019-09-03T15:03:54Z urn:sha1:e66125fe260deee49fdf6e9978d9bd29871dd5bb openjp2/j2k: Report error if all wanted components are not decoded. Change opj_j2k_check_poc_val() to take into account tile number 2019-04-25T13:06:45Z Even Rouault even.rouault@spatialys.com 2019-04-25T12:07:46Z urn:sha1:a94cfbd5334922ca5b63cfac9d2e5e0ec98155be Fix POC in multi-tile scenarios: avoid almost endless loop when a tile has no POC settings 2019-04-25T12:40:56Z Even Rouault even.rouault@spatialys.com 2019-04-24T23:27:02Z urn:sha1:6423163141412cb93364de4e33d90bcffefa0885 opj_j2k_check_poc_val(): prevent potential write outside of allocated array 2019-04-25T12:40:56Z Even Rouault even.rouault@spatialys.com 2019-04-24T22:34:44Z urn:sha1:23883458b9de2c57fc1890b42efbd0832c8fbe3b opj_j2k_check_poc_val(): fix starting index for checking layer dimension 2019-04-25T12:40:55Z Even Rouault even.rouault@spatialys.com 2019-04-24T22:28:05Z urn:sha1:6589c609f6d6b3743715fceefbdac6e4ecb76aee The standard mandates that the layer index always starts at zero for every progression.