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 <libxml++/libxml++.h>
  21 #include <dcp/colour_matrix.h>
  22 #include <dcp/raw_convert.h>
  23 #include <libcxml/cxml.h>
  24 #include "config.h"
  25 #include "colour_conversion.h"
  26 #include "util.h"
  27 #include "md5_digester.h"
  28
  29 #include "i18n.h"
  30
  31 using std::list;
  32 using std::string;
  33 using std::stringstream;
  34 using std::cout;
  35 using std::vector;
  36 using boost::shared_ptr;
  37 using boost::optional;
  38 using dcp::raw_convert;
  39
  40 ColourConversion::ColourConversion ()
  41         : input_gamma (2.4)
  42         , input_gamma_linearised (true)
  43         , matrix (3, 3)
  44         , output_gamma (2.6)
  45 {
  46         for (int i = 0; i < 3; ++i) {
  47                 for (int j = 0; j < 3; ++j) {
  48                         matrix (i, j) = dcp::colour_matrix::srgb_to_xyz[i][j];
  49                 }
  50         }
  51 }
  52
  53 ColourConversion::ColourConversion (double i, bool il, double const m[3][3], double o)
  54         : input_gamma (i)
  55         , input_gamma_linearised (il)
  56         , matrix (3, 3)
  57         , output_gamma (o)
  58 {
  59         for (int i = 0; i < 3; ++i) {
  60                 for (int j = 0; j < 3; ++j) {
  61                         matrix (i, j) = m[i][j];
  62                 }
  63         }
  64 }
  65
  66 ColourConversion::ColourConversion (cxml::NodePtr node)
  67         : matrix (3, 3)
  68 {
  69         input_gamma = node->number_child<double> ("InputGamma");
  70         input_gamma_linearised = node->bool_child ("InputGammaLinearised");
  71
  72         for (int i = 0; i < 3; ++i) {
  73                 for (int j = 0; j < 3; ++j) {
  74                         matrix (i, j) = 0;
  75                 }
  76         }
  77
  78         list<cxml::NodePtr> m = node->node_children ("Matrix");
  79         for (list<cxml::NodePtr>::iterator i = m.begin(); i != m.end(); ++i) {
  80                 int const ti = (*i)->number_attribute<int> ("i");
  81                 int const tj = (*i)->number_attribute<int> ("j");
  82                 matrix(ti, tj) = raw_convert<double> ((*i)->content ());
  83         }
  84
  85         output_gamma = node->number_child<double> ("OutputGamma");
  86 }
  87
  88 void
  89 ColourConversion::as_xml (xmlpp::Node* node) const
  90 {
  91         node->add_child("InputGamma")->add_child_text (raw_convert<string> (input_gamma));
  92         node->add_child("InputGammaLinearised")->add_child_text (input_gamma_linearised ? "1" : "0");
  93
  94         for (int i = 0; i < 3; ++i) {
  95                 for (int j = 0; j < 3; ++j) {
  96                         xmlpp::Element* m = node->add_child("Matrix");
  97                         m->set_attribute ("i", raw_convert<string> (i));
  98                         m->set_attribute ("j", raw_convert<string> (j));
  99                         m->add_child_text (raw_convert<string> (matrix (i, j)));
 100                 }
 101         }
 102
 103         node->add_child("OutputGamma")->add_child_text (raw_convert<string> (output_gamma));
 104 }
 105
 106 optional<size_t>
 107 ColourConversion::preset () const
 108 {
 109         vector<PresetColourConversion> presets = Config::instance()->colour_conversions ();
 110         size_t i = 0;
 111         while (i < presets.size() && (presets[i].conversion != *this)) {
 112                 ++i;
 113         }
 114
 115         if (i >= presets.size ()) {
 116                 return optional<size_t> ();
 117         }
 118
 119         return i;
 120 }
 121
 122 string
 123 ColourConversion::identifier () const
 124 {
 125         MD5Digester digester;
 126
 127         digester.add (input_gamma);
 128         digester.add (input_gamma_linearised);
 129         for (int i = 0; i < 3; ++i) {
 130                 for (int j = 0; j < 3; ++j) {
 131                         digester.add (matrix (i, j));
 132                 }
 133         }
 134         digester.add (output_gamma);
 135
 136         return digester.get ();
 137 }
 138
 139 PresetColourConversion::PresetColourConversion ()
 140         : name (_("Untitled"))
 141 {
 142
 143 }
 144
 145 PresetColourConversion::PresetColourConversion (string n, double i, bool il, double const m[3][3], double o)
 146         : name (n)
 147         , conversion (i, il, m, o)
 148 {
 149
 150 }
 151
 152 PresetColourConversion::PresetColourConversion (cxml::NodePtr node)
 153         : conversion (node)
 154 {
 155         name = node->string_child ("Name");
 156 }
 157
 158 void
 159 PresetColourConversion::as_xml (xmlpp::Node* node) const
 160 {
 161         conversion.as_xml (node);
 162         node->add_child("Name")->add_child_text (name);
 163 }
 164
 165 static bool
 166 about_equal (double a, double b)
 167 {
 168         static const double eps = 1e-6;
 169         return fabs (a - b) < eps;
 170 }
 171
 172 bool
 173 operator== (ColourConversion const & a, ColourConversion const & b)
 174 {
 175         if (
 176                 !about_equal (a.input_gamma, b.input_gamma) ||
 177                 a.input_gamma_linearised != b.input_gamma_linearised ||
 178                 !about_equal (a.output_gamma, b.output_gamma)) {
 179                 return false;
 180         }
 181
 182         for (int i = 0; i < 3; ++i) {
 183                 for (int j = 0; j < 3; ++j) {
 184                         if (!about_equal (a.matrix (i, j), b.matrix (i, j))) {
 185                                 return false;
 186                         }
 187                 }
 188         }
 189
 190         return true;
 191 }
 192
 193 bool
 194 operator!= (ColourConversion const & a, ColourConversion const & b)
 195 {
 196         return !(a == b);
 197 }