%% 3. NCPR score and UACI values for image data
function [n1 n2 n3 u1 u2 u3]=npcr_uaci(A,B)
[M N C]=size(A);
[M N C]=size(B);
 A=uint8(A);
 B=uint8(B);
A255=double(A)/double(255);
B255=double(B)/double(255);
D_R=0;
D_G=0;
D_B=0;
%
C_R=0;
C_G=0;
C_B=0;
for i=1:M
    for j=1:N
      if A(i,j,1)~=B(i,j,1)
          D_R=D_R+1;
          C_R=C_R+abs(A255(i,j,1)-B255(i,j,1));
         % C_R=C_R+abs(double(A(i,j,1))-double(B(i,j,1)))/double(255);%or
      end
      if A(i,j,2)~=B(i,j,2) 
          D_G=D_G+1;
          C_G=C_G+abs(A255(i,j,2)-B255(i,j,2));
      end
      if A(i,j,3)~=B(i,j,3) 
          D_B=D_B+1;
          C_B=C_B+abs(A255(i,j,3)-B255(i,j,3));
      end
    end
end
n1=D_R*100/(M*N);
n2=D_G*100/(M*N);
n3=D_B*100/(M*N);
u1=C_R*100/(M*N);
u2=C_G*100/(M*N);
u3=C_B*100/(M*N);
end

