function?mi=CalculateMI（ab）
%{
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%








%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%}
a=rgb2gray（a）;
b=rgb2gray（b）;
a=double（a）;
b=double（b）;
[MN]=size（a）;
hab=zeros（256256）;
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=