%計算MutualEntropy函數計算兩幅圖像的互信息熵
%輸入：已讀入matlab的兩幅圖像，的變量名
function?y=MutualEntropy（x1x2）
a=x1;
b=x2;
[MaNa]=size（a）;
[MbNb]=size（b）;
M=min（MaMb）;
N=min（NaNb）;
%初始化直方圖數組
hab=zeros（256256）;
ha=zeros（1256）;
hb=zeros（1256）;
%歸一化
maxA=max（max（a））;
maxA=maxA（1）;
minA=min（min（a））;
minA=minA（1）;
if?maxA~=minA;
????a=（a-minA）/（maxA-minA）;
else
????a=zeros（MN）;
end
maxB=max（max（b））;
maxB=maxB（1）;
minB=min（min（b））;
minB=minB（1）;
if?maxB~=minB;
????b=（b-minB）/（maxB-minB）;
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=a（ij）;
????????indexy=b（ij）;
????????hab（indexxindexy）=hab（indexxindexy）+1;%聯合直方圖
????????ha（indexx）=ha（indexx）+1;%a圖直方圖
????????hb（indexy）=hb（indexy）+1;%b圖直方圖
????end
end
%計算a的信息熵
hsum=sum（sum（ha））;
index=find（ha~=0）;
p=ha/hsum;
Ha=sum（sum（-p（index）.*log（p（index））））;
%計算b圖信息熵
hsum=sum（sum（hb））;
index=find（hb~=0）;
p=hb/hsum;
Hb=sum（sum（-p（index）.*log（p（index））））;

%計算聯合信息熵
hsum=sum（sum（hab））;
index=find（hab~=0）;
p=hab/hsum;
Hab=sum（sum（-p（index）.*log（p（index））））;
%計算互信息熵
y=Ha+Hb-Hab;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????766??2012-06-19?03:17??car?Identification\1.jpg

?????文件???????2143??2012-06-19?03:17??car?Identification\1.車牌灰度圖像.jpg

?????文件????????791??2012-06-19?03:17??car?Identification\2.jpg

?????文件???????3226??2012-06-19?03:17??car?Identification\2.車牌邊緣檢測.jpg

?????文件????????607??2012-06-19?03:17??car?Identification\3.jpg

?????文件???????3149??2012-06-19?03:17??car?Identification\3.車牌二值處理.jpg

?????文件????????773??2012-06-19?03:17??car?Identification\4.jpg

?????文件???????3043??2012-06-19?03:17??car?Identification\4.車牌均值濾波.jpg

?????文件????????789??2012-06-19?03:17??car?Identification\5.jpg

?????文件???????2974??2012-06-19?03:17??car?Identification\5.車牌腐蝕處理.jpg

?????文件????????777??2012-06-19?03:17??car?Identification\6.jpg

?????文件????????653??2012-06-19?03:17??car?Identification\7.jpg

?????文件?????108279??2004-01-15?05:14??car?Identification\Car.jpg

?????文件???????1239??2011-10-26?23:18??car?Identification\MutualEntropy.m

?????文件??????11751??2012-06-19?00:13??car?Identification\untitled1.fig

?????文件??????19111??2012-06-19?02:38??car?Identification\untitled1.m

?????文件???????4890??2011-10-27?17:22??car?Identification\字母模板\0.bmp

?????文件???????4578??2011-10-27?17:22??car?Identification\字母模板\1.bmp

?????文件???????4578??2011-10-27?17:26??car?Identification\字母模板\10.bmp

?????文件???????3954??2011-10-27?17:26??car?Identification\字母模板\11.bmp

?????文件???????6918??2011-10-27?17:27??car?Identification\字母模板\12.bmp

?????文件???????5358??2011-10-27?17:27??car?Identification\字母模板\13.bmp

?????文件???????5670??2011-10-27?17:28??car?Identification\字母模板\14.bmp

?????文件???????5202??2011-10-27?17:28??car?Identification\字母模板\15.bmp

?????文件???????5670??2011-10-27?17:28??car?Identification\字母模板\16.bmp

?????文件???????5514??2011-10-27?17:29??car?Identification\字母模板\17.bmp

?????文件???????4422??2011-10-27?17:29??car?Identification\字母模板\18.bmp

?????文件???????5202??2011-10-27?17:30??car?Identification\字母模板\19.bmp

?????文件???????4578??2011-10-27?17:23??car?Identification\字母模板\2.bmp

?????文件???????4422??2011-10-27?18:41??car?Identification\字母模板\20.bmp

............此處省略64個文件信息