Page 132 - Computational Colour Science Using MATLAB
P. 132

IMPLEMENTATIONS AND EXAMPLES                      119

                  plot(DRGB(:,3),NRGB(:,3),’b*’)
                  y1 = compgog(gogvals(3,:),x1);
                  hold on
                  plot(x1,y1,’b-’)

                  disp(’gog values’)
                  disp(gogvals)


                  RGB = zeros(3+length(N(:,1)),3);
                  RDACS = [R(1); G(1); B(1); N(:,1)]/255;
                  RGB(:,1) = compgog(gogvals(1,:), RDACS);
                  RDACS = [R(2); G(2); B(2); N(:,2)]/255;
                  RGB(:,2) = compgog(gogvals(2,:), RDACS);
                  RDACS = [R(3); G(3); B(3); N(:,3)]/255;
                  RGB(:,3) = compgog(gogvals(3,:), RDACS);

                  XYZ = (A*RGB’)’;


                  AXYZ = [R(4:6); G(4:6); B(4:6); N(:,4:6)];

                  for i = 1:8
                    [lab1] = xyz2lab(XYZ(i,:),’d65___64’);
                    [lab2] = xyz2lab(AXYZ(i,:),’d65___64’);
                    [thisDE] = cielabde(lab1,lab2);
                    de(i) = thisDE;
                  end

                  disp (’known values’)
                  disp(de)

                  RGB = zeros(length(T(:,1)),3);
                  RDACS = [T(:,1)]/255;
                  RGB(:,1) = compgog(gogvals(1,:), RDACS);
                  RDACS = [T(:,2)]/255;
                  RGB(:,2) = compgog(gogvals(2,:), RDACS);
                  RDACS = [T(:,3)]/255;
                  RGB(:,3) = compgog(gogvals(3,:), RDACS);

                  XYZ = (A*RGB’)’;

                  AXYZ = [T(:,4:6)];
                  clear de
   127   128   129   130   131   132   133   134   135   136   137