資源描述:
《直方圖均衡報(bào)告》由會(huì)員上傳分享,免費(fèi)在線(xiàn)閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、這幾次上機(jī)編程,由于我基礎(chǔ)太差,都沒(méi)有做出自己的程序,因此我改變策略,先讀懂老師給的程序,這也是算是一點(diǎn)收獲吧。下面就老師的程序,我解讀一下。%直方圖均衡clcpicture=input('請(qǐng)輸入圖片文件名:','s');輸入一個(gè)黑白圖像文件fid=fopen(picture);打開(kāi)該文件dat=fread(fid);將該文件中的像素點(diǎn)讀出來(lái)生成一個(gè)單列矩陣。d=dat';轉(zhuǎn)置成一個(gè)一行的行矩陣。h=length(d);newname=input('輸入新的文件名:','s');為后面生成的新圖像文件命名N=ze
2、ros(1,256);該循環(huán)將h這個(gè)行矩陣中各灰度級(jí)的像素點(diǎn)加權(quán)統(tǒng)計(jì)存入N這個(gè)零矩陣中。N矩陣存放的就是各灰度級(jí)的像素?cái)?shù)。fori=0:255forj=1:hifi==d(j)N(i+1)=N(i+1)+1;endendendfigure(1)通過(guò)對(duì)N矩陣的循環(huán)訪(fǎng)問(wèn),各灰度級(jí)像素?cái)?shù)除以總像素?cái)?shù)h,就得到各灰度級(jí)的概率。也就是原始直方圖。i=0:255;p=N/h;stem(i,p,'.')title('原始直方圖')m=zeros(1,256);fork=1:256該循環(huán)是求Sk的過(guò)程,很簡(jiǎn)單。ifk==1m(k
3、)=p(k);elsem(k)=m(k-1)+p(k);endendfork=1:256r=round(m*255);量化Sk得到Sk‘ends=zeros(1,256);fora=1:256forb=1:256ifr(a)==r(b)s(a)=N(a)+N(b);該循環(huán)統(tǒng)計(jì)量化后個(gè)灰度級(jí)像素?cái)?shù)。endendendq=s/h;量化后各灰度級(jí)概率figure(2)forj=1:256k(j)=r(j);endstem(k,q,'.');均衡直方圖title('均衡直方圖');M=zeros(1,h);將均衡后的像素
4、存入M中fori=1:hk=d(i);M(i)=r(k+1);這一句不理解,希望得到指導(dǎo)。endjht=fopen(newname,'wb');創(chuàng)建新圖片文件g=fwrite(jht,M);fclose('all');我對(duì)直方圖均衡的理解:在數(shù)學(xué)上通過(guò)對(duì)圖像中各灰度級(jí)像素依次往后的累加概率進(jìn)行量化,使得原先靠近的灰度級(jí)合并,從而使得合并后該灰度級(jí)像素的增多,該灰度級(jí)得到加強(qiáng)。實(shí)質(zhì)是犧牲了灰度級(jí),拉大了圖像對(duì)比度。電信科0803班呂鵬輝0807080324我的郵箱P616942659@163.com