ColorMapCreator.m

Go to the documentation of this file.

/* (Autoinserted by mtoc++)
 * This source code has been filtered by the mtoc++ executable,
 * which generates code that can be processed by the doxygen documentation tool.
 *
 * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
 * Except for the comments, the function bodies of your M-file functions are untouched.
 * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
 * attached source files that are highly readable by humans.
 *
 * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
 * the correct locations in the source code browser.
 * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
 */

class ColorMapCreator
  :public handle {
    public: /* ( Constant ) */

        static const Resolution = 1000;

        
    public:

        LogPlot = false;

        
    private:

        vals = "[]";

        cols = {""};

        
    public:


        function  useJet(vals) {
            this.addColor(vals(1),[0 0 1]);
            this.addColor(vals(2),[0 1 1]);
            this.addColor(vals(3),[1 1 0]);
            this.addColor(vals(4),[1 0 0]);
        }


        function  addColor(value,color,colorbefore) {
            if nargin < 4
                colorbefore = color;
            end
            this.vals(end+1) = value;
            this.cols[end+1] = color;
            if ~isequal(color, colorbefore)
                this.vals(end+1) = value-1/this.Resolution;
                this.cols[end+1] = colorbefore;
            end
        }


        function cmap = create(Zdata) {
            dvals = this.vals;
            dcols = this.cols;

            if this.LogPlot
                dvals = log10(dvals);
            end

            if nargin > 1
                if this.LogPlot
                    Zdata = log10(Zdata);
                end
                minZ = min(Zdata(:));
                dvals(end+1) = minZ;
                dcols[end+1] = [0 0 .5]; /* zeros(3,1) */


                maxZ = max(Zdata(:));
                dvals(end+1) = maxZ;
                dcols[end+1] = [.5 0 0]; /* zeros(3,1); */


                rem = dvals < minZ;
                dvals(rem)  = [];
                dcols(rem) = [];
                rem = dvals > maxZ;
                dvals(rem)  = [];
                dcols(rem) = [];
            end

            /*  Colormap */
            cmap = [];
            if length(dvals) > 1
                [val, idx] = sort(dvals," ascend ");
                valrange = linspace(val(1),val(end),this.Resolution);
                ncol = length(valrange);
                r = zeros(ncol,1);
                g = zeros(ncol,1);
                b = zeros(ncol,1);

                curpos = 1;
                curcol = dcols[idx(1)];
                for k = 2:length(val)
                    nextcol = dcols[idx(k)];
                    [~,nextpos] = min(abs(valrange-val(k)));

                    pos = curpos:nextpos;
                    r(pos) = linspace(curcol(1),nextcol(1),length(pos));
                    g(pos) = linspace(curcol(2),nextcol(2),length(pos));
                    b(pos) = linspace(curcol(3),nextcol(3),length(pos));

                    curpos = nextpos;
                    curcol = nextcol;
                end
                cmap = [r g b];
            end
        }

    

};