數(shù)字圖像處理實(shí)驗(yàn)報(bào)告-直方圖規(guī)定化0528

1、數(shù)字圖像處理實(shí)驗(yàn)報(bào)告-直方 圖規(guī)定化數(shù)字圖像處理實(shí)驗(yàn)報(bào)告Sanz1直方圖規(guī)定化實(shí)驗(yàn)名稱: 姓名：班級(jí)學(xué)號(hào)專(zhuān)業(yè)電子信息工程（2+2）指導(dǎo)教師：陳華華實(shí)驗(yàn)日期：2012年5月24日直方圖匹配(規(guī)定化)直方圖均衡化能夠自動(dòng)增強(qiáng)整個(gè)圖像的對(duì)比度，但它的具體增強(qiáng)效果 不容易控制，處理的結(jié)果總是得到全局均勻化的直方圖。實(shí)際上有時(shí)需要 變換直方圖，使之成為某個(gè)特定的形狀，從而有選擇地增強(qiáng)某個(gè)灰度值范 圍內(nèi)的對(duì)比度。這時(shí)可以采用比較靈活的直方圖規(guī)定化。一般來(lái)說(shuō)正確地選擇規(guī)定化的函數(shù)可以獲得比直方圖均衡化更好的效果。所謂直方圖規(guī)定化，就是通過(guò)一個(gè)灰度映像函數(shù)，將原灰度直方圖改 造成所希望的直方圖。所以，直方圖修

2、正的關(guān)鍵就是灰度映像函數(shù)。直方圖匹配方法主要有 3個(gè)步驟(這里設(shè) M和N分別為原始圖和規(guī)定 圖中的灰度級(jí)數(shù)，且只考慮NKM的情況)：(1)如同均衡化方法中，對(duì)原始圖的直方圖進(jìn)行灰度均衡化：(2)規(guī)定需要的直方圖，并計(jì)算能使規(guī)定的直方圖均衡化的變換：(3)將第1個(gè)步驟得到的變換反轉(zhuǎn)過(guò)來(lái)，即將原始直方圖對(duì)應(yīng)映射到 規(guī)定的直方圖，也就是將所有pf (fi )對(duì)應(yīng)到pu( uj)去。A圖直方圖規(guī)定B圖Fil Cd Vic Inj( Toe Desk, Wind fto *膈人的屈像尸F(xiàn)J FigureF住髭立工壬R舊 Edit事處r-ttll Tcok4 W ild口璃中口聞昌ll瞰國(guó)力 及 口國(guó)產(chǎn):

3、-舊消Matlab 程序:%直方圖規(guī)定化clear allA=imread(C:UsershpDesktopA.tif);%讀入A圖像imshow(A)%顯示出來(lái)title(輸入的A圖像)%繪制直方圖m,n=size(A);%測(cè)量圖像尺寸B=zeros(1,256);%預(yù)創(chuàng)建存放灰度出現(xiàn)概率的向量for k=0:255B(k+1)=length(find(A=k)/(m*n);% 計(jì)算每級(jí)灰度出現(xiàn)的概率，將其存入B中相應(yīng)位置 endfigure,bar(0:255,B,g);%繪制直方圖title(A圖像直方圖)xlabel(灰度值)ylabel(出現(xiàn)概率)axis(0,260,0,0.015

4、)C=imread(C:UsershpDesktopB.tif);% 讀入 B圖像imshow(C)%顯示出來(lái)title(輸入的B圖像)%繪制直方圖m,n=size(C);%測(cè)量圖像尺寸D=zeros(1,256);%預(yù)創(chuàng)建存放灰度出現(xiàn)概率的向量for k=0:255D(k+1)=length(find(C=k)/(m*n);% 計(jì)算每級(jí)灰度出現(xiàn)的概率，將其存入 D中相應(yīng)位 置 endfigure,bar(0:255,D,g);%繪制直方圖 title(B圖像直方圖) xlabel(灰度值) ylabel(出現(xiàn)概率) axis(0,260,0,0.015)S1=zeros(1,256);fo門(mén)

5、=1:256for j=1:iS1(i)=D(j)+S1(i);% 計(jì)算 B 灰度圖累計(jì)直方圖 end endcounts=Bfigure,bar(0:255,counts,r)title(A圖像直方圖) S2=zeros(1,256);fo門(mén)=1:256for j=1:iS2(i)=counts(j)+S2(i); endend;%累計(jì)”規(guī)定化直方圖%對(duì)比直方圖，找到相差最小的灰度級(jí) fo門(mén)=1:256for j=1:255if S1(j)=S2(i) ifabs(S1(j)-S2(i)0for k=(T(i-1)+1):T(i)H(i)=H(i)+D(k);endelse H(i)=0;e

6、ndendfigure,bar(0:255,H,g)%顯示規(guī)定化后的直方圖title(A規(guī)定B后的直方圖)xlabel(灰度值)ylabel(出現(xiàn)概率)axis(0,260,0,0.6)%顯示規(guī)定圖PA=C;%將各個(gè)像素歸一化后的灰度值賦給這個(gè)像素for i=1:mfor j=1:nfor k=2:255if T(k-1)=C(i,j)&C(i,j)=T(k)PA(i,j尸k； break;endendendendfigure,imshow(PA)%顯示規(guī)定化后的圖像title(A規(guī)定B后圖像)imwrite(PA,guidinghua.bmp);二、用已知直方圖規(guī)定A圖規(guī)定灰度為,4zero

7、s(1,49),03zeros(1,49),0.1,zeros(1,49),0.2,zeros(1,49),0.1,zeros(19),0.2,zeros(1,49),0.1. Figure- 2ini. Edrf ，；f ：i%e 。& *4ta& ehDIFie Ed I 5m wen TcqI& DeUitoc Window. Hlp，Ur. 盛*星羽E擊 yw biCiri rode DhllepHlp d 口白k做戲心g 口目口嶷范化良方爸 。等Me idit 中電 kxrf Tgk DmL r&* 出履H臣Ie段反9啟K二I我定七后井通方三D1Q5Matlab 程序：clear a

8、llA=imread(C:UsershpDesktopB.tif);%讀入A圖像imshow(A)%顯示出來(lái)title(輸入的A圖像)%繪制直方圖m,n=size(A);%測(cè)量圖像尺寸B=zeros(1,256);%預(yù)創(chuàng)建存放灰度出現(xiàn)概率的向量%計(jì)B中相應(yīng)位置for k=0:255B(k+1)=length(find(A=k)/(m*n);算每級(jí)灰度出現(xiàn)的概率，將其存入endfigure,bar(0:255,B,g);繪制直方圖title(A圖像直方圖,)xlabel(灰度值)ylabel(出現(xiàn)概率)axis(0,260,0,0.015)S1=zeros(1,256);fo門(mén)=1:256for

9、 j=1:iS1(i尸B(j)+S1(i);%計(jì)算原灰度圖累計(jì)直方圖endendcounts=zeros(1,49),0.3,zeros(1,49),0.1,zeros(1,49 ),02zeros(1,49),0.1,zeros(1,49),02zeros(1,49),0. 1;%規(guī)定化直方圖figure,bar(1:300,counts,r)title(規(guī)定化直方圖)S2=zeros(1,256);for i=1:256for j=1:iS2(i)=counts(j)+S2(i);endend;%累計(jì)規(guī)定化直方圖%對(duì)比直方圖，找到相差最小的灰度級(jí)fo門(mén)=1:256for j=1:256if S1(j)=S2ifabs(S1(j)-S2(i)0for k=(T(i-1)+1):T(i)H(i尸H(i)+B(k);endelse H(i)=0;endendfigure,bar(0:255,H,g)%顯示規(guī)定化后的直方圖title(規(guī)定化后的直方圖)xlabel(灰度值)ylabel(出現(xiàn)概率)axis(0,260,0,0.6)%顯示規(guī)定圖PA=A；%將各個(gè)像素歸一化后的灰度值賦給這個(gè)像素for