Various small packages for Matlab.zip

function [ RGBOUT, LINEPROPOUT ] = colourscale( varargin ) %COLOURSCALE Make colourful, nice-looking plots % This function takes the current figure and applies a series of colours % to each "line". These colours are a spectrum of saturations and % intensities for a given colour hue. % % COLOURSCALE(...,'ch',CH) % Use CH to specify the mode of operation: % CH is empty - apply colours to any lines present in the active axis % except lines for which UserData = 'colourscale:ignore' % CH is array of line handles - apply colours to specified lines % CH is an integer - generate CH number of colours (RBGOUT) and line % properties (LINEPROPOUT) to be used later % % COLOURSCALE(...,'hue',H) % Use hue H for colour scheme (default 20). % H is a standard "HSV" hue, from 0 to 100, where approx.: % H = 0 - red % H = 10 - orange % H = 15 - yellow % H = 35 - green % H = 55 - light blue % H = 60 - dark blue % H = 75 - purple % H = 85 - magenta % H = 100 - red again % % COLOURSCALE(...,'chroma',C) % Use chroma C for colour scheme (default 70). % Chroma appears to be a nonlinear parameter with sensible maximum; values % around 40 (dull) to 100 (bright) appear to be best, although higher than % this produces brighter colours they also start clipping what is possible % represent in RGB. % % COLOURSCALE(...,'lumin',[l_N L_N]) % Use [l_N L_N] as the range for lumin values to vary over. Lumin values % can range from 0 (dark/black) to 100 (light). % % COLOURSCALE(...,'lumin',{[l_1 L_1] [l_2 L_2] ... [l_M L_M]}) % For N colours, use [l_N L_N] as the range for lumin values to vary over. % This approach isn't the most convenient for the user but allows the most % flexibility, as different lumin ranges appears best for different values % of N, and for different chroma/hue combinations. % If N>M then [l_M L_M] is used as the range. % % COLOURSCALE(...,'linewidth',[LW1 LW2]) % If not specified, the plots take on their "natural" linewidth as default % or as specified by the user. If set to a two-element vector, the % linewidths of the lines will be set to vary linearly from LW1 to LW2 as % the plots change colour from dark to light. (This is useful as lighter % lines often need to be slightly thicker to remain visible compared to % darker lines.) % % COLOURSCALE(...,'repeat',N) % An optional argument specifies the number of times to use the % colour space: e.g., colourscale(2) will turn, in a graph with 6 data % series, the first and fourth plot blue, the second and fifth % green, and the third and six red. The divisor of the number of % plots and the number of colour space repetitions must be an % integer. % % COLOURSCALE(...,'permute',P) % By default the lines are coloured in the order in which they were plot. % This order can be changed by specifying a permutation of the order using % indexing, such as in a four-plot graph: % colourscale(...,'permute',[1 3 2 4]) % % RGBOUT = COLOURSCALE(...) % As above, and also returns the colours in a CH x 3 array. % % [RGBOUT,LINEPROPSOUT] = COLOURSCALE(...) % As above, and also returns the line properties in a CH x 2 or CH x 4 % cell array. Example: % [RGBOUT,LINEPROPSOUT] = COLOURSCALE(...,'ch',2); % plot(x1,y1,LINEPROPSOUT{1,:}); % plot(x2,y2,LINEPROPSOUT{2,:}); %% COLOURSCALE v0.1 % % Copyright (c) 2017-2018 Will Robertson % All rights reserved. % Licence (BSD) appended. % % Please report bugs and feature requests for % this package at the development repository: % <http://github.com/wspr/matlab-plot-tools/> %% Option parsing p = inputParser; p.addOptional('ch',findobj(gca,'Type','line','-not','UserData','colourscale:ignore'),@(x) all(isgraphics(x,'line')) || (isscalar(x)&&isnumeric(x))); p.addOptional('hue',20,@isnumeric); p.addOptional('chroma',70,@isnumeric); p.addOptional('repeat',1,@isnumeric); p.addOptional('permute',[],@isnumeric); p.addOptional('lumin',{[65 65] [50 80] [40 80] [30 90]},@(x) iscell(x) || isnumeric(x)); p.addOptional('linewidth',[],@isnumeric); p.parse(varargin{:}); ch = p.Results.ch; hue = p.Results.hue; chroma = p.Results.chroma; Nseries = p.Results.repeat; permute = p.Results.permute; lumin = p.Results.lumin; lw_range = p.Results.linewidth; %% Option massaging if ~isempty(lw_range) if numel(lw_range) == 1 lw_range = lw_range([1 1]); end end if isnumeric(lumin) lumin = {lumin}; end if isscalar(ch) && isnumeric(ch) do_apply_colours = false; Nch = ch; else do_apply_colours = true; Nch = numel(ch); end Ncol = Nch/Nseries; Nlum = numel(lumin); if round(Ncol) ~= Ncol % Each set of data series must be the same length to avoid rounding problems!! disp(['There are ',num2str(Nch),' data series']) error('There must be an integer multiple of specified data series in the figure.') end if isempty(permute) permute = 1:Nch; else if ~isequal(sort(permute),1:Nch) error('2nd argument must be a permutation of 1:N where N is the number of colours.'); end end %% Calculate colours % indexing into lumin values needs a trick. % let's say we have lumin values of [65 50 40 30]; % for n=1, lumin=65; n=3, lumin=40; etc. % for n=6, the index is too high, so we want n=4, which is min([n,Nlum]): lumin_index = min([Ncol,Nlum]); lmin = lumin{lumin_index}(1); lmax = lumin{lumin_index}(2); % for linewidths we just do linear interpolation, no need for the indexing % as in the above: if ~isempty(lw_range) lw = linspace(lw_range(1),lw_range(2),Ncol); end hcl = nan(Ncol,3); hcl(:,1) = hue*3.6; hcl(:,2) = chroma; hcl(:,3) = linspace(lmin,lmax,Ncol)'; rgb = nan(size(hcl)); for ii = 1:Ncol rgb(ii,:) = hcl2rgb(hcl(ii,1),hcl(ii,2),hcl(ii,3))'; end rgb = rgb/255; %% Assign colours for ii = 1:Nch ind = mod(ii-1,Ncol)+1; if isempty(lw_range) lineprop(permute(ii),:) = { 'Color',rgb(ind,:),... }; else lineprop(permute(ii),:) = { 'Color',rgb(ind,:),... 'LineWidth',lw(ind),... }; end if do_apply_colours set(ch(permute(ii)), lineprop{permute(ii),:}); end end %% Fin if nargout > 0 RGBOUT = rgb; end if nargout > 1 LINEPROPOUT = lineprop; end return %% Licence % % Distributed under the BSD licence in accordance with the wishes of the % Matlab File Exchange. % % Redistribution and use in source and binary forms, with or without % modification, are permitted provided that the following conditions are met: % * Redistributions of source code must retain the above copyright % notice, this list of conditions and the following disclaimer. % * Redistributions in binary form must reproduce the above copyright % notice, this list of conditions and the following disclaimer in the % documentation and/or other materials provided with the distribution. % % THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDER ''AS IS'' AND ANY % EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED % WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE % DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER BE LIABLE FOR ANY % DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES % (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; % LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND % ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT % (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF % THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. function rgb = hcl2rgb(h, c, l) %HCL2RGB Convert a HCL (i.e., CIELUV) color space value to one % in sRGB space. % RGB = HCL2RGB(H, C, L) will convert the color (H, C, L) in % HCL color space to RGB = [R, G, B] in sRGB color space. % Values that lie outside sRGB space will be silently corrected. % Code