Fix incorrect colourspace conversion of XYZ content
[dcpomatic.git] / src / lib / j2k_image_proxy.cc
1 /*
2     Copyright (C) 2014-2015 Carl Hetherington <cth@carlh.net>
3
4     This program is free software; you can redistribute it and/or modify
5     it under the terms of the GNU General Public License as published by
6     the Free Software Foundation; either version 2 of the License, or
7     (at your option) any later version.
8
9     This program is distributed in the hope that it will be useful,
10     but WITHOUT ANY WARRANTY; without even the implied warranty of
11     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12     GNU General Public License for more details.
13
14     You should have received a copy of the GNU General Public License
15     along with this program; if not, write to the Free Software
16     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
17
18 */
19
20 #include "j2k_image_proxy.h"
21 #include "dcpomatic_socket.h"
22 #include "image.h"
23 #include "data.h"
24 #include "raw_convert.h"
25 #include <dcp/openjpeg_image.h>
26 #include <dcp/mono_picture_frame.h>
27 #include <dcp/stereo_picture_frame.h>
28 #include <dcp/colour_conversion.h>
29 #include <dcp/rgb_xyz.h>
30 #include <libcxml/cxml.h>
31 #include <libxml++/libxml++.h>
32 #include <iostream>
33
34 #include "i18n.h"
35
36 using std::string;
37 using std::cout;
38 using boost::shared_ptr;
39 using boost::optional;
40 using boost::dynamic_pointer_cast;
41
42 /** Construct a J2KImageProxy from a JPEG2000 file */
43 J2KImageProxy::J2KImageProxy (boost::filesystem::path path, dcp::Size size)
44         : _data (path)
45         , _size (size)
46 {
47
48 }
49
50 J2KImageProxy::J2KImageProxy (shared_ptr<const dcp::MonoPictureFrame> frame, dcp::Size size)
51         : _data (frame->j2k_size ())
52         , _size (size)
53 {
54         memcpy (_data.data().get(), frame->j2k_data(), _data.size ());
55 }
56
57 J2KImageProxy::J2KImageProxy (shared_ptr<const dcp::StereoPictureFrame> frame, dcp::Size size, dcp::Eye eye)
58         : _size (size)
59         , _eye (eye)
60 {
61         switch (eye) {
62         case dcp::EYE_LEFT:
63                 _data = Data (frame->left_j2k_size ());
64                 memcpy (_data.data().get(), frame->left_j2k_data(), _data.size ());
65                 break;
66         case dcp::EYE_RIGHT:
67                 _data = Data (frame->right_j2k_size ());
68                 memcpy (_data.data().get(), frame->right_j2k_data(), _data.size ());
69                 break;
70         }
71 }
72
73 J2KImageProxy::J2KImageProxy (shared_ptr<cxml::Node> xml, shared_ptr<Socket> socket)
74 {
75         _size = dcp::Size (xml->number_child<int> ("Width"), xml->number_child<int> ("Height"));
76         if (xml->optional_number_child<int> ("Eye")) {
77                 _eye = static_cast<dcp::Eye> (xml->number_child<int> ("Eye"));
78         }
79         _data = Data (xml->number_child<int> ("Size"));
80         socket->read (_data.data().get (), _data.size ());
81 }
82
83 void
84 J2KImageProxy::ensure_j2k () const
85 {
86         if (!_j2k) {
87                 _j2k = dcp::decompress_j2k (const_cast<uint8_t*> (_data.data().get()), _data.size (), 0);
88         }
89 }
90
91 shared_ptr<Image>
92 J2KImageProxy::image (optional<dcp::NoteHandler> note) const
93 {
94         ensure_j2k ();
95
96         if (_j2k->opj_image()->comps[0].prec < 12) {
97                 int const shift = 12 - _j2k->opj_image()->comps[0].prec;
98                 for (int c = 0; c < 3; ++c) {
99                         int* p = _j2k->data (c);
100                         for (int y = 0; y < _j2k->size().height; ++y) {
101                                 for (int x = 0; x < _j2k->size().width; ++x) {
102                                         *p++ <<= shift;
103                                 }
104                         }
105                 }
106         }
107
108         shared_ptr<Image> image (new Image (pixel_format(), _size, true));
109
110         /* Copy data in whatever format (sRGB or XYZ) into our Image; I'm assuming
111            the data is 12-bit either way.
112         */
113
114         int p = 0;
115         for (int y = 0; y < _j2k->size().height; ++y) {
116                 uint16_t* q = (uint16_t *) (image->data()[0] + y * image->stride()[0]);
117                 for (int x = 0; x < _j2k->size().width; ++x) {
118                         for (int c = 0; c < 3; ++c) {
119                                 *q++ = _j2k->data(c)[p] << 4;
120                         }
121                         ++p;
122                 }
123         }
124
125         return image;
126 }
127
128 void
129 J2KImageProxy::add_metadata (xmlpp::Node* node) const
130 {
131         node->add_child("Type")->add_child_text (N_("J2K"));
132         node->add_child("Width")->add_child_text (raw_convert<string> (_size.width));
133         node->add_child("Height")->add_child_text (raw_convert<string> (_size.height));
134         if (_eye) {
135                 node->add_child("Eye")->add_child_text (raw_convert<string> (_eye.get ()));
136         }
137         node->add_child("Size")->add_child_text (raw_convert<string> (_data.size ()));
138 }
139
140 void
141 J2KImageProxy::send_binary (shared_ptr<Socket> socket) const
142 {
143         socket->write (_data.data().get(), _data.size());
144 }
145
146 bool
147 J2KImageProxy::same (shared_ptr<const ImageProxy> other) const
148 {
149         shared_ptr<const J2KImageProxy> jp = dynamic_pointer_cast<const J2KImageProxy> (other);
150         if (!jp) {
151                 return false;
152         }
153
154         if (_data.size() != jp->_data.size()) {
155                 return false;
156         }
157
158         return memcmp (_data.data().get(), jp->_data.data().get(), _data.size()) == 0;
159 }
160
161 AVPixelFormat
162 J2KImageProxy::pixel_format () const
163 {
164         ensure_j2k ();
165
166         if (_j2k->opj_image()->color_space == CLRSPC_SRGB) {
167                 return AV_PIX_FMT_RGB48LE;
168         }
169
170         return AV_PIX_FMT_XYZ12LE;
171 }
172
173 J2KImageProxy::J2KImageProxy (Data data, dcp::Size size)
174         : _data (data)
175         , _size (size)
176 {
177
178 }