數(shù)字圖像處理作業(yè)(第三次)直方圖處理

1、精選優(yōu)質(zhì)文檔-傾情為你奉上數(shù)字圖像處理作業(yè)（第三次）直方圖處理第一題第二題(例3.6)1. 問(wèn)題：1）編寫函數(shù)Y=MyHisteq(X), 對(duì)灰度圖象X進(jìn)行直方圖均衡化，得到均衡化后的圖像Y，假設(shè)X的數(shù)據(jù)類型為uint8。 2）利用自己編寫的直方圖均衡化函數(shù)，重現(xiàn)例3.6中的實(shí)驗(yàn)，即對(duì)圖3.16中的圖像進(jìn)行直方圖均衡化。2程序：(1) 函數(shù)程序：function Y = MyHisteq(X)%對(duì)灰度圖象X進(jìn)行直方圖均衡化，得到均衡化后的圖像Y m,n=size(X); A=zeros(1,256); for k=0:255 A(k+1)=length(find(X=k)/(m*n); end

2、 S1=zeros(1,256); for i=1:256 for j=1:i S1(i)=A(j)+S1(i); %算Sk end end S2=round(S1*256); for i=1:256 B(i)=sum(A(find(S2=i); end figure,bar(0:255,B,'b'); title('均衡化后的直方圖'); xlabel('灰度'),ylabel('出現(xiàn)概率') Y=X; for i=0:255 Y(find(X=i)=S2(i+1); endend(2) 主程序：X=imread('Fi

3、g0316(1)(top_left).tif');figure,imshow(X);Y=MyHisteq(X);figure,imshow(Y);3. 實(shí)驗(yàn)結(jié)果：第三題(例3.12)1. 問(wèn)題：編寫程序，重現(xiàn)例3.12中的實(shí)驗(yàn)，產(chǎn)生圖3.27。2程序：X=imread('Fig0327(a)(tungsten_original).tif');figure,imshow(X);Y=histeq(X);figure,imshow(Y);E=4.0;k0=0.4;k1=0.02;k2=0.4;m,n=size(X);X1=X;Mg=sum(X1(:)/(m*n);Vg=sum

4、(sum(X1-Mg).2)/(m*n);for i=2:m-1for j=2:n-1X2=X1(i-1:i+1,j-1:j+1);Ms=sum(X2(:)/9;Vs=sum(sum(X2-Ms).2)/9;if(Ms<=k0*Mg&&k1*Vg<=Vs&&Vs<=k2*Vg)X1(i,j)=E*X1(i,j);endendendfigure,imshow(X1);3. 實(shí)驗(yàn)結(jié)果：放大了約130倍的鎢絲的SME圖像：全局直方圖均衡的結(jié)果：使用局部直方圖統(tǒng)計(jì)增強(qiáng)后的圖像：1. 問(wèn)題：用自己編寫編程實(shí)現(xiàn)中值濾波，用96頁(yè)例3.14的圖測(cè)試。2 程序

5、：(1) 函數(shù)程序：function Y=mid_filter(X,n)%自編的中值濾波函數(shù)。X是需要濾波的圖像,n是模板大小(即n×n) height, width=size(X); %輸入圖像是p×q的,且p>n,q>n X1=double(X); X2=X1; for i=1:height-n+1 for j=1:height-n+1 C=X1(i:i+(n-1),j:j+(n-1); %取出X1中從(i,j)開(kāi)始的n行n列元素,即模板(n×n的) E=C(1,:); %是C矩陣的第一行 for k=2:n E=E,C(k,:); %將C矩陣變?yōu)橐粋€(gè)行矩陣 end mm=median(E); %mm是中值 X2(i+(n-1)/2,j+(n-1)/2)=mm;%將模板各元素的中值賦給模板中心位置的元素 end end Y=uint8(X2);end(2) 主程序：X=imread('Fig0335(a)(ckt_board_saltpep_prob_pt05).tif');figure,imshow(X);n=3;