matlab计算图像间的相似度程序

%%%%  读入图片   %%%%%%%%%%%%%%%%%%%%%%%%
% clc;clear;close;
F1=imread（‘SE3_t613.bmp‘）;
G1=imread（‘SE3_t613.bmp‘）;
%%%%%  计算图像F1和G1的互信息   %%%%%%%%%%%%%%%%%%%%
F1_1=rgb2gray（F1）;
G1_1=rgb2gray（G1）;
Pf_1=imhist（F1_1）./numel（F1_1）;
Pg_1=imhist（G1_1）./numel（G1_1）;
Pfg_1=zeros（256256）;
for i=1:512
    for j=1:512
        Pfg_1（F1_1（ij）+1G1_1（ij）+1）=Pfg_1（F1_1（ij）+1G1_1（ij）+1）+1;
    end
end
num=sum（Pfg_1（:））;
Pfg_1=Pfg_1/num;
W1_1=zeros（1256）;
W2_1=zeros（1256）;
W12_1=zeros（1256）;
for m=1:256
    if Pf_1（m）~=0
       W1_1（m）=Pf_1（m）.*log10（Pf_1（m））;
    end
end
for m=1:256
    if Pg_1（m）~=0
       W2_1（m）=Pg_1（m）.*log10（Pg_1（m））;
    end
end
for m=1:256
    if Pfg_1（m）~=0
       W12_1（m）=Pfg_1（m）.*log10（Pfg_1（m））;
    end
end
N1_1=-sum（W1_1（:））;
N2_1=-sum（W2_1（:））;
N12_1=-sum（W12_1（:））;