clear;
close all;
A0=imread('1.bmp');
A0=im2double(A0);
A7=100*A0;
[a1,a2]=size(A7);
A2=sum(sum(A7));
a=A2/(a1*a2);
A3=ones(a1,a2);
A5=A3.*a;
%A6=uint8(A5);
A6=A5;
A=A7-A6;
Ah=imhist(A7);
non_zeros_index=find(Ah(:)~=0);
maxgrayA=max(non_zeros_index);
disp(maxgrayA);
B0=imread('2.bmp');
B0=im2double(B0);
B7=100*B0;
[b1,b2]=size(B7);
B2=sum(sum(B7));
b=B2/(b1*b2);
B3=ones(b1,b2);
B5=B3.*b;
%B6=uint8(B5);
B6=B5;
B=B7-B6;
Bh=imhist(B);
non_zeros_index1=find(Bh(:)~=0);
maxgrayB=max(non_zeros_index1);
disp(maxgrayB);
figure;
subplot(221),imshow(A7);
subplot(222),imshow(B7);
subplot(223),imshow(A);
subplot(224),imshow(B);
E=A.*B;
F=sum(sum(E));
K=sum(dot(A,A));
L=sum(dot(B,B));
M=K*L;
P=sqrt(M);
D=F/P;
%D=F/sqrt(sum(dot(A,A))*sum(dot(B,B)));
disp(D);
G=corr2(A0,B0);
disp(G);
%disp(F);
%disp(K);
%disp(L);
%disp(M);
%disp(P);
%whos;