src/lib/colour_conversion.cc

   1 /*
   2     Copyright (C) 2013 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 "config.h"
  21 #include "colour_conversion.h"
  22 #include "util.h"
  23 #include "md5_digester.h"
  24 #include "raw_convert.h"
  25 #include <dcp/colour_matrix.h>
  26 #include <dcp/gamma_transfer_function.h>
  27 #include <dcp/modified_gamma_transfer_function.h>
  28 #include <libcxml/cxml.h>
  29 #include <libxml++/libxml++.h>
  30
  31 #include "i18n.h"
  32
  33 using std::list;
  34 using std::string;
  35 using std::cout;
  36 using std::vector;
  37 using boost::shared_ptr;
  38 using boost::optional;
  39 using boost::dynamic_pointer_cast;
  40
  41 ColourConversion::ColourConversion ()
  42         : dcp::ColourConversion (dcp::ColourConversion::srgb_to_xyz ())
  43 {
  44
  45 }
  46
  47 ColourConversion::ColourConversion (dcp::ColourConversion conversion_)
  48         : dcp::ColourConversion (conversion_)
  49 {
  50
  51 }
  52
  53 ColourConversion::ColourConversion (cxml::NodePtr node, int version)
  54 {
  55         shared_ptr<dcp::TransferFunction> in;
  56
  57         if (version >= 32) {
  58
  59                 /* Version 2.x */
  60
  61                 cxml::ConstNodePtr in_node = node->node_child ("InputTransferFunction");
  62                 string in_type = in_node->string_child ("Type");
  63                 if (in_type == "Gamma") {
  64                         _in.reset (new dcp::GammaTransferFunction (false, in_node->number_child<double> ("Gamma")));
  65                 } else if (in_type == "ModifiedGamma") {
  66                         _in.reset (new dcp::ModifiedGammaTransferFunction (
  67                                            false,
  68                                            in_node->number_child<double> ("Power"),
  69                                            in_node->number_child<double> ("Threshold"),
  70                                            in_node->number_child<double> ("A"),
  71                                            in_node->number_child<double> ("B")
  72                                            ));
  73                 }
  74
  75         } else {
  76
  77                 /* Version 1.x */
  78
  79                 if (node->bool_child ("InputGammaLinearised")) {
  80                         _in.reset (new dcp::ModifiedGammaTransferFunction (false, node->number_child<float> ("InputGamma"), 0.04045, 0.055, 12.92));
  81                 } else {
  82                         _in.reset (new dcp::GammaTransferFunction (false, node->number_child<float> ("InputGamma")));
  83                 }
  84         }
  85
  86         list<cxml::NodePtr> m = node->node_children ("Matrix");
  87         for (list<cxml::NodePtr>::iterator i = m.begin(); i != m.end(); ++i) {
  88                 int const ti = (*i)->number_attribute<int> ("i");
  89                 int const tj = (*i)->number_attribute<int> ("j");
  90                 _matrix(ti, tj) = raw_convert<double> ((*i)->content ());
  91         }
  92
  93         _out.reset (new dcp::GammaTransferFunction (true, node->number_child<double> ("OutputGamma")));
  94 }
  95
  96 boost::optional<ColourConversion>
  97 ColourConversion::from_xml (cxml::NodePtr node, int version)
  98 {
  99         if (!node->optional_node_child ("InputTransferFunction")) {
 100                 return boost::optional<ColourConversion> ();
 101         }
 102
 103         return ColourConversion (node, version);
 104 }
 105
 106 void
 107 ColourConversion::as_xml (xmlpp::Node* node) const
 108 {
 109         xmlpp::Node* in_node = node->add_child ("InputTransferFunction");
 110         if (dynamic_pointer_cast<const dcp::GammaTransferFunction> (_in)) {
 111                 shared_ptr<const dcp::GammaTransferFunction> tf = dynamic_pointer_cast<const dcp::GammaTransferFunction> (_in);
 112                 in_node->add_child("Type")->add_child_text ("Gamma");
 113                 in_node->add_child("Gamma")->add_child_text (raw_convert<string> (tf->gamma ()));
 114         } else if (dynamic_pointer_cast<const dcp::ModifiedGammaTransferFunction> (_in)) {
 115                 shared_ptr<const dcp::ModifiedGammaTransferFunction> tf = dynamic_pointer_cast<const dcp::ModifiedGammaTransferFunction> (_in);
 116                 in_node->add_child("Type")->add_child_text ("ModifiedGamma");
 117                 in_node->add_child("Power")->add_child_text (raw_convert<string> (tf->power ()));
 118                 in_node->add_child("Threshold")->add_child_text (raw_convert<string> (tf->threshold ()));
 119                 in_node->add_child("A")->add_child_text (raw_convert<string> (tf->A ()));
 120                 in_node->add_child("B")->add_child_text (raw_convert<string> (tf->B ()));
 121         }
 122
 123         boost::numeric::ublas::matrix<double> matrix = _matrix;
 124         for (int i = 0; i < 3; ++i) {
 125                 for (int j = 0; j < 3; ++j) {
 126                         xmlpp::Element* m = node->add_child("Matrix");
 127                         m->set_attribute ("i", raw_convert<string> (i));
 128                         m->set_attribute ("j", raw_convert<string> (j));
 129                         m->add_child_text (raw_convert<string> (matrix (i, j)));
 130                 }
 131         }
 132
 133         node->add_child("OutputGamma")->add_child_text (raw_convert<string> (dynamic_pointer_cast<const dcp::GammaTransferFunction> (_out)->gamma ()));
 134 }
 135
 136 optional<size_t>
 137 ColourConversion::preset () const
 138 {
 139         vector<PresetColourConversion> presets = Config::instance()->colour_conversions ();
 140         size_t i = 0;
 141         while (i < presets.size() && (presets[i].conversion != *this)) {
 142                 ++i;
 143         }
 144
 145         if (i >= presets.size ()) {
 146                 return optional<size_t> ();
 147         }
 148
 149         return i;
 150 }
 151
 152 string
 153 ColourConversion::identifier () const
 154 {
 155         MD5Digester digester;
 156
 157         if (dynamic_pointer_cast<const dcp::GammaTransferFunction> (_in)) {
 158                 shared_ptr<const dcp::GammaTransferFunction> tf = dynamic_pointer_cast<const dcp::GammaTransferFunction> (_in);
 159                 digester.add (tf->gamma ());
 160         } else if (dynamic_pointer_cast<const dcp::ModifiedGammaTransferFunction> (_in)) {
 161                 shared_ptr<const dcp::ModifiedGammaTransferFunction> tf = dynamic_pointer_cast<const dcp::ModifiedGammaTransferFunction> (_in);
 162                 digester.add (tf->power ());
 163                 digester.add (tf->threshold ());
 164                 digester.add (tf->A ());
 165                 digester.add (tf->B ());
 166         }
 167
 168         boost::numeric::ublas::matrix<double> matrix = _matrix;
 169         for (int i = 0; i < 3; ++i) {
 170                 for (int j = 0; j < 3; ++j) {
 171                         digester.add (matrix (i, j));
 172                 }
 173         }
 174
 175         digester.add (dynamic_pointer_cast<const dcp::GammaTransferFunction> (_out)->gamma ());
 176
 177         return digester.get ();
 178 }
 179
 180 PresetColourConversion::PresetColourConversion ()
 181         : name (_("Untitled"))
 182 {
 183
 184 }
 185
 186 PresetColourConversion::PresetColourConversion (string n, dcp::ColourConversion conversion_)
 187         : conversion (conversion_)
 188         , name (n)
 189 {
 190
 191 }
 192
 193 PresetColourConversion::PresetColourConversion (cxml::NodePtr node, int version)
 194         : conversion (node, version)
 195         , name (node->string_child ("Name"))
 196 {
 197
 198 }
 199
 200 void
 201 PresetColourConversion::as_xml (xmlpp::Node* node) const
 202 {
 203         conversion.as_xml (node);
 204         node->add_child("Name")->add_child_text (name);
 205 }
 206
 207 bool
 208 operator== (ColourConversion const & a, ColourConversion const & b)
 209 {
 210         return a.about_equal (b, 1e-6);
 211 }
 212
 213 bool
 214 operator!= (ColourConversion const & a, ColourConversion const & b)
 215 {
 216         return !(a == b);
 217 }
 218
 219 bool
 220 operator== (PresetColourConversion const & a, PresetColourConversion const & b)
 221 {
 222         return a.name == b.name && a.conversion == b.conversion;
 223 }