%?計算互信息的函數
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?=?