uv0fun.m

自己写的计算色坐标的小程序

function [u,v,Y]=UV0fun(beta,measurelight,standardlight)
%%%%%%%%%%%%%%%%%%%%%        color science 1 %%%   ----caculate the  X Y Z K  W U V
%%%%%%%%%%%%%%%%%%%%%      processing  and   estimating  the  lighe bulb           
%%%%%%%%    measurelight matrix is  77X1
%%%%%%%%    standardlight matrix is 77X3
%%%%%%%     beta  matrix is 77X8
%%%         k123 is (x,y,z)  %%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%standardlight=loadstlight./10000        %%%% load 1931  CIE standard color watcher x()  y()  z()  
%measurelight=loadmeasurelight   
%beta=loadsample
k123=double(standardlight);
p=double(measurelight);

%%%   p is the power of the light source   %%%%%%%

[a,b]=size(k123);      %%%%%%%%%%  a,b is the 77X3
[a1,b1]=size(beta);
k2=zeros(1,8);         %%%% k2 is the middle variable 1X8

if beta==0 
    f=1:a;
    beta(f,1)=1;
else  
end

 for j=1:b1
   for i=1:a
       k2(1,j)=k2(1,j)+k123(i,2)*p(i)*beta(i,j);
   end
end
        

ky=zeros(1,b1);                   %%%%%%%%% ky is the cofficient  which is used to caculated Y X Z
 for i=1:b1
     if  k2(i)~=0
         ky(i)=100./k2(i);
     else
         ky(i)=0;
     end
 end

%%%%%%%%%%%%%%%%%%     X Y Z       %%%
X=zeros(b1,1);  Y=zeros(b1,1);  Z=zeros(b1,1);
for j=1:b1
for i=1:a1
X(j)=X(j)+ky(j)*k123(i,1)*p(i);
Y(j)=Y(j)+ky(j)*k123(i,2)*p(i);
Z(j)=Z(j)+ky(j)*k123(i,3)*p(i);
end
end

%%%%%%%%%%%%%%  x,y,z,  U ,  V  %%%%%%%%
x=zeros(b1,1);
y=zeros(b1,1);
z=zeros(b1,1);
U=zeros(b1,1);
V=zeros(b1,1);
for i=1:b1;
      if ky(i)~=0
        x(i,1)=X(i,1)/(X(i,1)+Y(i,1)+Z(i,1));
        y(i,1)=Y(i,1)/(X(i,1)+Y(i,1)+Z(i,1));
        z(i,1)=Z(i,1)/(X(i,1)+Y(i,1)+Z(i,1));

        U(i,1)=4*X(i,1)/(X(i,1)+15*Y(i,1)+3*Z(i,1));
        V(i,1)=6*Y(i,1)/(X(i,1)+15*Y(i,1)+3*Z(i,1));
    end
end
x
y
z
X
Y
Z

if U(2)==0
   u=U(1)
    v=V(1)
    Y=Y(1)
    
elseif U(2)~=0
    u=U
    v=V
    Y=Y
  end