aboutsummaryrefslogtreecommitdiffstats
path: root/Software/Visual_Studio/SideChains/Colourful/Conversion/VonKriesChromaticAdaptation.cs
blob: 5a6b4999cf1cf95c6ca3fc0a4489213ce3e88a4f (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
using System;
using Colourful.Implementation;
using Colourful.Implementation.Conversion;
using Matrix = System.Collections.Generic.IReadOnlyList<System.Collections.Generic.IReadOnlyList<double>>;

namespace Colourful.Conversion
{
    /// <summary>
    /// Basic implementation of the von Kries chromatic adaptation model
    /// </summary>
    /// <remarks>
    /// Transformation described here:
    /// http://www.brucelindbloom.com/index.html?Eqn_ChromAdapt.html
    /// </remarks>
    public sealed class VonKriesChromaticAdaptation : IChromaticAdaptation
    {
        private readonly IColorConversion<XYZColor, LMSColor> _conversionToLMS;
        private readonly IColorConversion<LMSColor, XYZColor> _conversionToXYZ;
        private Matrix _cachedDiagonalMatrix;

        private XYZColor _lastSourceWhitePoint;
        private XYZColor _lastTargetWhitePoint;

        /// <summary>
        /// Constructs von Kries chromatic adaptation using default <see cref="XYZAndLMSConverter" />
        /// </summary>
        public VonKriesChromaticAdaptation() : this(new XYZAndLMSConverter())
        {
        }

        /// <summary>
        /// Transformation matrix used for the conversion (definition of the cone response domain).
        /// <see cref="LMSTransformationMatrix" />
        /// </summary>
        public VonKriesChromaticAdaptation(Matrix transformationMatrix) : this(new XYZAndLMSConverter(transformationMatrix))
        {
        }

        private VonKriesChromaticAdaptation(XYZAndLMSConverter converter) : this(converter, converter)
        {
        }

        /// <summary>
        /// Constructs von Kries chromatic adaptation using given converters
        /// </summary>
        public VonKriesChromaticAdaptation(IColorConversion<XYZColor, LMSColor> conversionToLMS, IColorConversion<LMSColor, XYZColor> conversionToXYZ)
        {
            _conversionToLMS = conversionToLMS ?? throw new ArgumentNullException(nameof(conversionToLMS));
            _conversionToXYZ = conversionToXYZ ?? throw new ArgumentNullException(nameof(conversionToXYZ));
        }

        /// <summary>
        /// Transforms XYZ color to destination reference white.
        /// </summary>
        public XYZColor Transform(in XYZColor sourceColor, in XYZColor sourceWhitePoint, in XYZColor targetWhitePoint)
        {
            if (sourceWhitePoint.Equals(targetWhitePoint))
                return sourceColor;

            var sourceColorLMS = _conversionToLMS.Convert(sourceColor);

            if (sourceWhitePoint != _lastSourceWhitePoint || targetWhitePoint != _lastTargetWhitePoint)
            {
                var sourceWhitePointLMS = _conversionToLMS.Convert(sourceWhitePoint);
                var targetWhitePointLMS = _conversionToLMS.Convert(targetWhitePoint);

                _cachedDiagonalMatrix = MatrixFactory.CreateDiagonal(targetWhitePointLMS.L / sourceWhitePointLMS.L, targetWhitePointLMS.M / sourceWhitePointLMS.M, targetWhitePointLMS.S / sourceWhitePointLMS.S);
                _lastSourceWhitePoint = sourceWhitePoint;
                _lastTargetWhitePoint = targetWhitePoint;
            }

            var targetColorLMS = new LMSColor(_cachedDiagonalMatrix.MultiplyBy(sourceColorLMS.Vector));
            var targetColor = _conversionToXYZ.Convert(targetColorLMS);
            return targetColor;
        }
    }
}