两幅图像的互信息的matlab代码

% 计算互信息的函数
function mi = MI（ab）
%Caculate MI of a and b in the region of the overlap part

%计算重叠部分
[MaNa] = size（a）;
[MbNb] = size（b）;
M=min（MaMb）;
N=min（NaNb）;

%初始化直方图数组
hab = zeros（256256）;
ha = zeros（1256）;
hb = zeros（1256）;

%归一化
if max（max（a））~=min（min（a））
    a = （a-min（min（a）））/（max（max（a））-min（min（a）））;
else
    a = zeros（MN）;
end

if max（max（b））-min（min（b））
    b = （b-min（min（b）））/（max（max（b））-min（min（b）））;
else
    b = zeros（MN）;
end

a = double（int16（a*255））+1;
b = double（int16（b*255））+1;

%统计直方图
for i=1:M
    for j=1:N
        indexx =