%Caculate?MI?of?a?and?b?in?the?region?of?the?size（a）

function?[mi]?=?EMI（ab）
%
[MN]?=?size（a）;
%?change?a?b?to?range?1:n
if?min（a（:））<1
????a?=?a-min（a（:））+1;
end
if?min（b（:））<1
????b?=?b-min（b（:））+1;
end
????

%?ncolora?=?max（a（:））-min（a（:））+1;
%?ncolorb?=?max（b（:））-min（b（:））+1;
ncolora?=?max（a（:））;
ncolorb?=?max（b（:））;

hab?=?zeros（ncolorancolorb）;
ha?=?zeros（1ncolora）;
hb?=?zeros（1ncolorb）;

%count?begins
[MN]=size（a）;
for?i=1:M
????for?j=1:N
???????indexx?=??a（ij）;
???????indexy?=??b（ij）;

???????hab（indexxindexy）?=?hab（indexxindexy）+1;
???????ha（indexx）?=?ha（indexx）+1;
???????hb（indexy）?=?hb（indexy）+1;

???end
end

%cal?hahbhab
hsum?=?sum（sum（hab））;
index?=?find（hab~=0）;
p?=?hab/hsum;
Hab?=?sum（sum（-p（index）.*log（p（index））））;

hsum?=?sum（sum（ha））;
index?=?find（ha~=0）;
p?=?ha/hsum;
Ha?=?sum（sum（-p（index）.*log（p（index））））;

hsum?=?sum（sum（hb））;
index?=?find（hb~=0）;
p?=?hb/hsum;
Hb?=?sum（sum（-p（index）.*log（p（index））））;

mi?=?Ha+Hb-Hab;
%?cal?EMI
%?feature?efficiency?coefficient
n?=2/5;
if?Hab?~=?0
????e?=?（mi^n）?/?（?Hab^（1-n））;
else%?not?usual?happen?just?for?test
????e?=?-1;%（mi^n）?/?（?Hab^（1-n）+1e-10）;
end
mi?=?e;
%?%?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2005-11-27?20:08??pso圖像配準(zhǔn)\
?????文件????????1218??2005-10-19?20:09??pso圖像配準(zhǔn)\EMI.m
?????文件?????????697??2005-10-19?20:09??pso圖像配準(zhǔn)\FPMI.m
?????文件????????3338??2005-10-19?20:10??pso圖像配準(zhǔn)\GMI.m
?????文件?????????228??2005-10-19?20:10??pso圖像配準(zhǔn)\gobad.m
?????文件????????1404??2005-10-19?20:10??pso圖像配準(zhǔn)\HiMI.m
?????文件????????1034??2005-10-19?20:10??pso圖像配準(zhǔn)\MI.m
?????文件?????????653??2005-10-19?20:10??pso圖像配準(zhǔn)\myMI.m
?????文件?????????439??2005-10-19?20:10??pso圖像配準(zhǔn)\oneDimSearch.m
?????文件????????2528??2005-10-19?20:10??pso圖像配準(zhǔn)\POWELL.m
?????文件???????35888??2005-10-19?20:10??pso圖像配準(zhǔn)\proj.fig
?????文件???????12544??2005-10-19?20:10??pso圖像配準(zhǔn)\proj.m
?????文件????????5902??2005-10-19?20:10??pso圖像配準(zhǔn)\PSO.m
?????文件?????????390??2005-10-19?20:10??pso圖像配準(zhǔn)\restore.m
?????文件????????3719??2005-11-27?20:04??pso圖像配準(zhǔn)\RMI.m
?????文件??????????67??2005-10-19?20:10??pso圖像配準(zhǔn)\test.m
?????文件?????????238??2005-10-19?20:10??pso圖像配準(zhǔn)\test2.m
?????文件????????7648??2005-10-19?20:10??pso圖像配準(zhǔn)\testListBoxlj.fig
?????文件????????7936??2005-11-27?20:02??pso圖像配準(zhǔn)\testListBoxlj.m
?????文件?????????326??2005-10-19?20:10??pso圖像配準(zhǔn)\說明.txt