summaryrefslogtreecommitdiff
path: root/src/j2k_transcode_oj1.cc
diff options
context:
space:
mode:
authorCarl Hetherington <cth@carlh.net>2024-04-06 22:34:10 +0200
committerCarl Hetherington <cth@carlh.net>2024-04-07 19:06:24 +0200
commitd8504c2d684128af07589ebe99f69321daeef37a (patch)
tree46088f0722aba99e81f49f1cbca02d01bdb94c30 /src/j2k_transcode_oj1.cc
parent0eff1f460e187c5d0e3f3829a266d8dc9ee000db (diff)
Re-add option to build with openjpeg 1.5.x.
Diffstat (limited to 'src/j2k_transcode_oj1.cc')
-rw-r--r--src/j2k_transcode_oj1.cc222
1 files changed, 222 insertions, 0 deletions
diff --git a/src/j2k_transcode_oj1.cc b/src/j2k_transcode_oj1.cc
new file mode 100644
index 00000000..6cd6e16d
--- /dev/null
+++ b/src/j2k_transcode_oj1.cc
@@ -0,0 +1,222 @@
+/*
+ Copyright (C) 2012-2015 Carl Hetherington <cth@carlh.net>
+
+ This file is part of libdcp.
+
+ libdcp is free software; you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation; either version 2 of the License, or
+ (at your option) any later version.
+
+ libdcp is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with libdcp. If not, see <http://www.gnu.org/licenses/>.
+
+ In addition, as a special exception, the copyright holders give
+ permission to link the code of portions of this program with the
+ OpenSSL library under certain conditions as described in each
+ individual source file, and distribute linked combinations
+ including the two.
+
+ You must obey the GNU General Public License in all respects
+ for all of the code used other than OpenSSL. If you modify
+ file(s) with this exception, you may extend this exception to your
+ version of the file(s), but you are not obligated to do so. If you
+ do not wish to do so, delete this exception statement from your
+ version. If you delete this exception statement from all source
+ files in the program, then also delete it here.
+*/
+
+
+#include "compose.hpp"
+#include "data.h"
+#include "dcp_assert.h"
+#include "exceptions.h"
+#include "j2k_transcode.h"
+#include "openjpeg_image.h"
+#include <openjpeg.h>
+#include <cmath>
+
+
+using std::min;
+using std::pow;
+using std::shared_ptr;
+using std::string;
+using namespace dcp;
+
+
+shared_ptr<dcp::OpenJPEGImage>
+dcp::decompress_j2k(Data const& data, int reduce)
+{
+ return dcp::decompress_j2k (data.data(), data.size(), reduce);
+}
+
+
+shared_ptr<dcp::OpenJPEGImage>
+dcp::decompress_j2k(shared_ptr<const Data> data, int reduce)
+{
+ return dcp::decompress_j2k(data->data(), data->size(), reduce);
+}
+
+
+
+/** Decompress a JPEG2000 image to a bitmap.
+ * @param data JPEG2000 data.
+ * @param size Size of data in bytes.
+ * @param reduce A power of 2 by which to reduce the size of the decoded image;
+ * e.g. 0 reduces by (2^0 == 1), ie keeping the same size.
+ * 1 reduces by (2^1 == 2), ie halving the size of the image.
+ * This is useful for scaling 4K DCP images down to 2K.
+ * @return XYZ image.
+ */
+shared_ptr<dcp::OpenJPEGImage>
+dcp::decompress_j2k(uint8_t const* data, int64_t size, int reduce)
+{
+ auto decoder = opj_create_decompress(CODEC_J2K);
+ opj_dparameters_t parameters;
+ opj_set_default_decoder_parameters(&parameters);
+ parameters.cp_reduce = reduce;
+ opj_setup_decoder(decoder, &parameters);
+ auto cio = opj_cio_open((opj_common_ptr) decoder, const_cast<uint8_t*>(data), size);
+ auto image = opj_decode(decoder, cio);
+ if (!image) {
+ opj_destroy_decompress(decoder);
+ opj_cio_close (cio);
+ boost::throw_exception(ReadError(String::compose ("could not decode JPEG2000 codestream of %1 bytes.", size)));
+ }
+
+ opj_destroy_decompress(decoder);
+ opj_cio_close(cio);
+
+ image->x1 = rint(float(image->x1) / pow (2, reduce));
+ image->y1 = rint(float(image->y1) / pow (2, reduce));
+ return std::make_shared<OpenJPEGImage>(image);
+}
+
+
+ArrayData
+dcp::compress_j2k(shared_ptr<const OpenJPEGImage> xyz, int bandwidth, int frames_per_second, bool threed, bool fourk, string comment)
+{
+ /* Set the max image and component sizes based on frame_rate */
+ int max_cs_len = ((float) bandwidth) / 8 / frames_per_second;
+ if (threed) {
+ /* In 3D we have only half the normal bandwidth per eye */
+ max_cs_len /= 2;
+ }
+ int const max_comp_size = max_cs_len / 1.25;
+
+ /* get a J2K compressor handle */
+ auto cinfo = opj_create_compress(CODEC_J2K);
+ if (cinfo == 0) {
+ throw MiscError ("could not create JPEG2000 encoder");
+ }
+
+ if (comment.empty()) {
+ /* asdcplib complains with "Illegal data size" when reading frames encoded with an empty comment */
+ throw MiscError("compress_j2k comment can not be an empty string");
+ }
+
+ /* Set encoding parameters to default values */
+ opj_cparameters_t parameters;
+ opj_set_default_encoder_parameters(&parameters);
+ if (fourk) {
+ parameters.numresolution = 7;
+ }
+
+ /* Set default cinema parameters */
+ parameters.tile_size_on = false;
+ parameters.cp_tdx = 1;
+ parameters.cp_tdy = 1;
+
+ /* Tile part */
+ parameters.tp_flag = 'C';
+ parameters.tp_on = 1;
+
+ /* Tile and Image shall be at (0,0) */
+ parameters.cp_tx0 = 0;
+ parameters.cp_ty0 = 0;
+ parameters.image_offset_x0 = 0;
+ parameters.image_offset_y0 = 0;
+
+ /* Codeblock size = 32x32 */
+ parameters.cblockw_init = 32;
+ parameters.cblockh_init = 32;
+ parameters.csty |= 0x01;
+
+ /* The progression order shall be CPRL */
+ parameters.prog_order = CPRL;
+
+ /* No ROI */
+ parameters.roi_compno = -1;
+
+ parameters.subsampling_dx = 1;
+ parameters.subsampling_dy = 1;
+
+ /* 9-7 transform */
+ parameters.irreversible = 1;
+
+ parameters.tcp_rates[0] = 0;
+ parameters.tcp_numlayers++;
+ parameters.cp_disto_alloc = 1;
+ parameters.cp_rsiz = fourk ? CINEMA4K : CINEMA2K;
+ parameters.cp_comment = strdup(comment.c_str());
+ if (fourk) {
+ parameters.numpocs = 2;
+ parameters.POC[0].tile = 1;
+ parameters.POC[0].resno0 = 0;
+ parameters.POC[0].compno0 = 0;
+ parameters.POC[0].layno1 = 1;
+ parameters.POC[0].resno1 = parameters.numresolution - 1;
+ parameters.POC[0].compno1 = 3;
+ parameters.POC[0].prg1 = CPRL;
+ parameters.POC[1].tile = 1;
+ parameters.POC[1].resno0 = parameters.numresolution - 1;
+ parameters.POC[1].compno0 = 0;
+ parameters.POC[1].layno1 = 1;
+ parameters.POC[1].resno1 = parameters.numresolution;
+ parameters.POC[1].compno1 = 3;
+ parameters.POC[1].prg1 = CPRL;
+ }
+
+ parameters.cp_comment = strdup ("libdcp");
+ parameters.cp_cinema = fourk ? CINEMA4K_24 : CINEMA2K_24;
+
+ /* 3 components, so use MCT */
+ parameters.tcp_mct = 1;
+
+ /* set max image */
+ parameters.max_comp_size = max_comp_size;
+ parameters.tcp_rates[0] = ((float) (3 * xyz->size().width * xyz->size().height * 12)) / (max_cs_len * 8);
+
+ /* Set event manager to null (openjpeg 1.3 bug) */
+ cinfo->event_mgr = 0;
+
+ /* Setup the encoder parameters using the current image and user parameters */
+ opj_setup_encoder(cinfo, &parameters, xyz->opj_image ());
+
+ opj_cio_t* cio = opj_cio_open((opj_common_ptr) cinfo, 0, 0);
+ if (cio == 0) {
+ opj_destroy_compress (cinfo);
+ throw MiscError("could not open JPEG2000 stream");
+ }
+
+ int const r = opj_encode(cinfo, cio, xyz->opj_image(), 0);
+ if (r == 0) {
+ opj_cio_close (cio);
+ opj_destroy_compress (cinfo);
+ throw MiscError("JPEG2000 encoding failed");
+ }
+
+ ArrayData enc(cio->buffer, cio_tell(cio));
+
+ opj_cio_close (cio);
+ free (parameters.cp_comment);
+ opj_destroy_compress (cinfo);
+
+ return enc;
+}
+