close?all;
subplot（222）;
%一，圖像的預處理，讀入彩色圖像將其灰度化
PS=imread（‘F:\lenna.bmp‘）;?????????????????%讀入彩色圖像文件?????????????????????????????
imshow（PS）
title（‘輸入的彩色圖像‘）
imwrite（rgb2gray（PS）‘PicSampleGray.bmp‘）;?%將彩色圖片灰度化并保存
PS=rgb2gray（PS）;????????????????????????????%灰度化后的數據存入數組

%二，繪制直方圖
[mn]=size（PS）;?????????????????????????????%測量圖像尺寸參數
GP=zeros（1256）;????????????????????????????%預創建存放灰度出現概率的向量
for?k=0:255
?????GP（k+1）=length（find（PS==k））/（m*n）;??????%計算每級灰度出現的概率，將其存入GP中相應位置
end
subplot（224）;
bar（0:255GP‘g‘）????????????????????%繪制直方圖
title（‘原圖像直方圖‘）
xlabel（‘灰度值‘）
ylabel（‘出現概率‘）
%三，直方圖均衡化
S1=zeros（1256）;
for?i=1:256
?????for?j=1:i
??????????S1（i）=GP（j）+S1（i）;?????????????????%計算Sk
?????end
end
S2=round（（S1*256）+0.5）;??????????????????????????%將Sk歸到相近級的灰度
for?i=1:256
?????GPeq（i）=sum（GP（find（S2==i）））;???????????%計算現有每個灰度級出現的概率
end
subplot（223）;
bar（0:255GPeq‘b‘）??????????????????%顯示均衡化后的直方圖
title（‘均衡化后的直方圖‘）
xlabel（‘灰度值‘）
ylabel（‘出現概率‘）

%四，圖像均衡化
PA=PS;
for?i=0:255
?????PA（find（PS==i））=S2（i+1）;????????????????%將各個像素歸一化后的灰度值賦給這個像素
end
subplot（221）;
imshow（PA）???????????????????????????%顯示均衡化后的圖像?
title（‘均衡化后圖像‘）
imwrite（PA‘PicEqual.bmp‘）;

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

?????文件???????1432??2009-05-31?01:33??3206003189丘譯明\juhenghua.m

?????文件????????168??2009-05-31?23:26??3206003189丘譯明\juhenghuaqym.m

?????文件?????196662??2009-04-10?22:39??3206003189丘譯明\3206003189丘譯明\lenna.bmp

?????文件???????2646??2009-06-01?14:57??3206003189丘譯明\3206003189丘譯明\test\test.clw

?????文件???????4263??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\ReadMe.txt

?????文件???????1334??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\test.h

?????文件???????4173??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\test.cpp

?????文件???????1054??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\StdAfx.h

?????文件????????206??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\StdAfx.cpp

?????文件???????1581??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\MainFrm.h

?????文件???????2505??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\MainFrm.cpp

?????文件???????1453??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\testDoc.h

?????文件???????1702??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\testDoc.cpp

?????文件????????531??2009-04-10?15:31??3206003189丘譯明\3206003189丘譯明\test\test.dsw

?????文件?????173056??2009-06-01?14:57??3206003189丘譯明\3206003189丘譯明\test\test.ncb

?????文件????????753??2009-06-01?14:55??3206003189丘譯明\3206003189丘譯明\test\test.plg

?????文件??????12358??2009-05-29?01:30??3206003189丘譯明\3206003189丘譯明\test\test.rc

?????文件????????824??2009-05-29?01:30??3206003189丘譯明\3206003189丘譯明\test\Resource.h

?????文件??????29592??2009-06-01?00:20??3206003189丘譯明\3206003189丘譯明\test\test.aps

?????文件???????1867??2009-05-29?01:45??3206003189丘譯明\3206003189丘譯明\test\testView.h

?????文件???????4758??2009-05-29?01:46??3206003189丘譯明\3206003189丘譯明\test\test.dsp

?????文件???????1325??2009-05-29?01:56??3206003189丘譯明\3206003189丘譯明\test\StatChart.h

?????文件????????971??2009-05-29?18:17??3206003189丘譯明\3206003189丘譯明\test\dib.h

?????文件???????8400??2009-05-29?18:17??3206003189丘譯明\3206003189丘譯明\test\dib.cpp

?????文件???????6908??2009-05-29?18:19??3206003189丘譯明\3206003189丘譯明\test\testView.cpp

?????文件?????230400??2009-06-01?14:55??3206003189丘譯明\3206003189丘譯明\test\Debug\vc60.idb

?????文件????5507044??2009-05-29?01:19??3206003189丘譯明\3206003189丘譯明\test\Debug\test.pch

?????文件?????372736??2009-06-01?00:28??3206003189丘譯明\3206003189丘譯明\test\Debug\vc60.pdb

?????文件?????105407??2009-05-29?01:19??3206003189丘譯明\3206003189丘譯明\test\Debug\StdAfx.obj

?????文件??????22607??2009-05-29?01:45??3206003189丘譯明\3206003189丘譯明\test\Debug\test.obj

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