File: LjgToXYZ.m

package info (click to toggle)
psychtoolbox-3 3.0.19.14.dfsg1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 86,796 kB
  • sloc: ansic: 176,245; cpp: 20,103; objc: 5,393; sh: 2,753; python: 1,397; php: 384; makefile: 193; java: 113
file content (103 lines) | stat: -rw-r--r-- 3,198 bytes parent folder | download | duplicates (2)
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
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
function XYZ = LjgToXYZ(Ljg)
% XYZ = LjgToXYZ(Ljg)
%
% Convert OSA Ljg to XYZ (10 degree).  Works by using numerical
% search to invert XYZToLjg.  See XYZToLjg for details on
% formulae used.  See also OSAUCSTest.
%
% This can return imaginary values if you pass XYZ values
% that are outside reasonable physical gamut limits.
%
% 3/27/01  dhb      Wrote it.
% 3/4/05   dhb      Handle new version of optimization toolbox, too.
% 9/23/12  dhb, ms  Update options for current Matlab versions.

% On Octave we need the 'optim' package v1.6.0 or later for fmincon() support,
% and Octave 6 or later for support for handles to nested functions like @LjgToXYZFun below:
if IsOctave
    v = version;
    if str2num(v(1)) < 6
        error('For use with Octave, you need at least Octave version 6.');
    end

    try
        % Try loading the optim package with the optimization functions:
        pkg load optim;
    catch
        error('For use with Octave, you must install the ''optim'' package from Octave Forge. See ''help pkg''.');
    end

    % Got optim package loaded. Does it support fmincon()?
    if ~exist('fmincon')
        error('For use with Octave, you need at least version 1.6.0 of the ''optim'' package from Octave Forge.');
    end
end



% Check for needed optimization toolbox, and version.
if (exist('fmincon') == 2)
    % Search options
    options = optimset('fmincon');
    options = optimset(options, 'Display', 'off', 'Algorithm', 'active-set');

    % Search bounds -- XYZ must be positive.
    vlb = zeros(3,1);
    vub = [200 100 200]';

    % Do the search for each passed value.
    XYZ = zeros(size(Ljg));
    n = size(XYZ,2);
    for i = 1:n
        x0 = xyYToXYZ([.28 .28 30]');
        LjgCurrent = Ljg(:,i);
        x = fmincon(@LjgToXYZFun,x0,[],[],[],[],vlb,vub,[],options);
        XYZ(:,i) = x;
    end
elseif (exist('constr') == 2)
    % Search options
    options = foptions;
    options(1) = 0;

    % Search bounds -- XYZ must be positive.
    vlb = zeros(3,1);
    vub = [200 100 200]';

    % Do the search for each passed value.
    XYZ = zeros(size(Ljg));
    n = size(XYZ,2);
    for i = 1:n
        x0 = xyYToXYZ([.28 .28 30]');
        LjgCurrent = Ljg(:,i);
        x = constr(@LjgToXYZFun,x0,options,vlb,vub);
        XYZ(:,i) = x;
    end
else
    error('LjgToXYZ requires the optional Matlab Optimization Toolbox from Mathworks');
end

        function [f,g] = LjgToXYZFun(XYZ)
            % [f,g] = LjgToXYZFun(XYZ)
            %
            % Nested optimization function for LjgToXYZ numerical
            % search. Not vectorized. Uses LjgCurrent from parent function.
            %
            % 3/27/01  dhb  Wrote it.
            % 12/4/20  mk   Turned it into nested function.

            % Convert back.
            Ljg1 = XYZToLjg(XYZ);

            % Handle case where XYZ is so weird that
            % an imaginary value is returned.
            if (any(~isreal(Ljg1)))
                sdiff = (abs(Ljg1)-LjgCurrent).^2;
                f = sum(sdiff);
                g = 10;
            else
                sdiff = (Ljg1-LjgCurrent).^2;
                f = sum(sdiff);
                g = -XYZ;
            end
        end
end