added LCH2RGB

setups/vestibular chair/vPrime/trial/experiment/pb_vRunVC.m

@@ -8,8 +8,8 @@ function Dat = pb_vRunVC(signal)
 % PBToolbox (2018): JJH: j.heckman@donders.ru.nl
 
    %% CREATE SIGNALS
-   vSignal = pb_vCreateSignal(1, signal(1).duration, 10, 1, signal(1).type);
-   hSignal = pb_vCreateSignal(1, signal(2).duration, 10, 1, signal(2).type);
+   vSignal = pb_vCreateSignal(1, signal(1).duration, 10, .1, signal(1).type);
+   hSignal = pb_vCreateSignal(1, signal(2).duration, 10, .1, signal(2).type);
  
    Dat.v.x = vSignal.x .* signal(1).amplitude;
    Dat.v.t = (0:1:length(Dat.v.x)-1)/10;

utilities/graphics/pb_LCH2RGB.m

+function RGB = pb_LCH2RGB(LCH)
+% PB_LCH2RGB()
+%
+% PB_LCH2RGB()  ...
+%
+% See also ...
+
+% PBToolbox (2018): JJH: j.heckman@donders.ru.nl
+
+   if nargin<1
+      ncol	= 64;
+      l		= linspace(0,1,ncol);
+      H		= repmat(300,1,ncol);
+      C		= zeros(1,ncol);
+      L		= 90-l*30;
+      LCH		= [L;C;H]';
+      LCH = LCH(1:15,:);
+      close all
+   end
+   
+   %% Convert CIE L*CH to CIE L*ab
+   L = LCH(:,1);
+   C = LCH(:,2);
+   H = LCH(:,3);
+   a = cos(H*pi/180).*C;
+   b = sin(H*pi/180).*C;
+
+   %% CIE-L*ab -> XYZ
+   % https://en.wikipedia.org/wiki/Lab_color_space
+   Y	= (L+16)/116;
+   X	= Y+a/500;
+   Z	= Y-b/200;
+   XYZ = [X Y Z];
+
+   sel			= XYZ>6/29;
+   XYZ(sel)	= XYZ(sel).^3;
+   XYZ(~sel)	= (XYZ(~sel)-4/29)*3*(6/29)^2;
+
+   % Reference white point
+   ref	=  [95.05, 100.000, 108.90]/100;
+   XYZ =  bsxfun(@times,XYZ,ref);
+
+   %% XYZ -> RGB
+   % Conversion/transformation matrix
+   % Matrix multiplication: C(i,j) = A(i,p)*B(p,j);
+   % https://en.wikipedia.org/wiki/CIE_1931_color_space
+   % http://en.wikipedia.org/wiki/SRGB
+   T		=	[3.2406, -1.5372, -0.4986;...
+            -0.9689, 1.8758, 0.0415;...
+            0.0557, -0.2040, 1.0570];
+   RGB		= XYZ*T';
+
+   %% Gamma correction to convert RGB to sRGB
+   sel			= RGB>0.0031308; % colorspace: 0.0031306684425005883
+   RGB(sel)	= 1.055*(RGB(sel).^(1/2.4))-0.055;
+   RGB(~sel)	= 12.92*RGB(~sel);
+
+   %% Final check
+   % [I,J]			= find(RGB>1); %#ok<*NASGU>
+   % RGB(I,:)		= repmat([1 1 1],numel(I),1);
+   % [I,J]			= find(RGB<0);
+   % RGB(I,:)		= repmat([1 1 1],numel(I),1);
+
+   %% Clip values
+   sel			= RGB>1;
+   RGB(sel)	= 1;
+   sel			= RGB<0;
+   RGB(sel)	= 0;
+
+
+
+
+ 
+end
+ 
+% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 
+%                                                           %
+%       Part of Programmeer Beer Toolbox (PBToolbox)        %
+%       Written by: Jesse J. Heckman (2018)                 %
+%                                                           %
+% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 
+

utilities/graphics/pb_statcolor.m

@@ -192,7 +192,7 @@ function Col = pb_statcolor(ncol,statmap,palette,Par,varargin)
          LCH		= [LCH1; LCH2];
 
    end
-   Col		= LCH2RGB(LCH);
+   Col		= pb_LCH2RGB(LCH);
 
    if dispFlag
       figure();