/* Copyright (C) 2012-2015 Carl Hetherington 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 . 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 #include using std::min; using std::pow; using std::shared_ptr; using std::string; using namespace dcp; shared_ptr dcp::decompress_j2k(Data const& data, int reduce) { return dcp::decompress_j2k (data.data(), data.size(), reduce); } shared_ptr dcp::decompress_j2k(shared_ptr 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::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(¶meters); parameters.cp_reduce = reduce; opj_setup_decoder(decoder, ¶meters); auto cio = opj_cio_open((opj_common_ptr) decoder, const_cast(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(image); } ArrayData dcp::compress_j2k(shared_ptr 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(¶meters); 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, ¶meters, 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; }