灰度圖像直方圖均衡化變換函數(shù)算法

《灰度圖像直方圖均衡化變換函數(shù)算法》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、灰度圖像直方圖均衡化變換函數(shù)算法方案一：1)灰度圖像直方圖均衡化算法步驟：1、列出原始圖像和變換后圖像的灰度級I：j=0,1,-L-1,其中L是灰度級的個(gè)數(shù)；2、統(tǒng)計(jì)原始圖像個(gè)灰度級的像素個(gè)數(shù)Ni;3、計(jì)算原始圖像直方圖：p(i)=Ni/N,N為原始圖像像素總數(shù)；4、計(jì)算累計(jì)直方圖Pj=》p(k),k二0,l???j；5、利用灰度變換函數(shù)計(jì)算變換后的灰度值，并四舍五入：j=INT[(L-l)Pj+0.5];6、確定灰度變換關(guān)系i->j,據(jù)此將原圖像的灰度值f(m,n)=i修正為g(m,n)=j;7、統(tǒng)計(jì)變換后各

2、灰度級的像素個(gè)數(shù)Nj；8、計(jì)算變換后圖像的直方圖p(j)二Nj/N;2)m文件I=imread('football.jpg1);1=1(:,:,1);[mzn]=size(I);num=m*n;pre_mat=zeros(1,256);aft_mat=zeros(1,256);fori=l:mforj=l:npre_mat(I(izj)+1)=pre_mat(I(izj)+1)+1;endendaft_mat(I)=aft_mat(I)/num;fork=2:256s(k)=pre_mat(k)/num;aft

3、_mat(k)=aft_mat(k-l)+s(k);endM=zeros(叫n);aft_mat=aft_mat*255;fori=l:mforj=l:nM(izj)=aft_mat(I(i,j)+1);endendJ=uint8(M)subplot(3,2fl)imshow(I);圖1均衡化處理對比1.原始圖像T;subplot(3,2,2)imhist(I);title(*2.原始圖像直方圖');subplot(3,2^3)imshow(J);title(!3.均衡化圖像')；subplot(3,2,4)i

4、mhist(J);title「4?均衡化圖像直方圖*);A=histeq(I);subplot(3,2^5)imshow(A);title(*5.系統(tǒng)均衡化圖像T;subplot(3,2,6)imhist(A);title「6?系統(tǒng)均衡化圖像直方圖1)經(jīng)實(shí)際驗(yàn)證：該方案對黑白圖片均衡化處理有較好的效果。不能正確的處理彩色圖像。如上例中所示：該方案先將彩色圖像轉(zhuǎn)換成黑白圖像，再對其均衡化處理。為了良好的處理彩色圖像，設(shè)計(jì)方案二如下所示：方案二PS=imread(*football.jpg');%讀入JPG彩色圖像

5、文件subplot(3,2,1);imshow(PS);%顯不圖像title(*原始彩色JPG圖像1)imwrite(rgb2gray(PS),1PicSampleGray.bmp');%將彩色圖片灰度化并保存PS=rgb2gray(PS);%灰度花后的數(shù)據(jù)存入數(shù)組$繪制直方圖[m,n]=size(PS);%測量圖像尺寸參數(shù)GP=zeros(1,256);%預(yù)創(chuàng)建存放灰度出現(xiàn)概率的向量fork=0:255GP(k+l)=length(find(PS==k))/(m*n)；%計(jì)算每級灰度出現(xiàn)的概率，并將其存入GP

6、中相應(yīng)位置endsubplot(3,2,2),bar(0:255,GPZ)%繪制直方圖title(*原圖像直方圖Jylabel(f出現(xiàn)概率1)%直方圖均衡化Sl=zeros(1,256);fori=l:256forj=l:isi(i)=GP(j)+Sl(i);endendS2=round((Sl*256)+0.5);fori=l:256GPeq(i)=sum(GP(find(S2==i)));endsubplot(3,2,3)fbar(0:255,GPeq,fb1)title(*均衡化后的直方圖1)ylabel

7、(*出現(xiàn)概率')axis([025500?025])%圖像均衡化PA=PS;fori=0:255PA(find(PS==i))=S2(i+1);endsubplot(3’2,4)fimshow(PA)title(*均衡化后圖像Jimwrite(PAZ1PicEqual.bmp1);%計(jì)算Sk%將$1＜?xì)w到相近級的灰度%計(jì)算現(xiàn)有每個(gè)灰度級岀現(xiàn)的概率%顯示均衡化后的直方圖%將各個(gè)像素歸一化后的灰度值賦給這個(gè)像素%顯示均衡化后的圖像K=histeq(PS);subplot(3,2,5),imhist(K);咎顯示系統(tǒng)

8、均衡化處理圖像的直方圖title(係統(tǒng)均衡化圖像直方圖，)ylabel(1tB現(xiàn)個(gè)數(shù)')subplot(3,2/6)imshow(K);%顯示系統(tǒng)自帶均衡化函數(shù)才處理的圖像title(!系統(tǒng)均衡化圖像TfflBtt■方H!300002fi-vs系據(jù)均荷化用a該方案針對彩色圖像,圖2彩色圖像均衡化對比不能處理黑白圖像。