實(shí)驗(yàn)3形態(tài)學(xué)圖像處理

1、 實(shí)驗(yàn)3 圖像增強(qiáng)及形態(tài)學(xué)圖像處理實(shí)驗(yàn)?zāi)康模?. 掌握均值濾波和中值濾波的原理及實(shí)現(xiàn)；2. 掌握銳化模板prewitt，sobel和laplacian的使用方法；3. 了解形態(tài)學(xué)的基本理論和方法；4. 掌握對圖像進(jìn)行膨脹/腐蝕的方法；5. 掌握開閉運(yùn)算；實(shí)驗(yàn)內(nèi)容：1、 （1）給圖像headCT分別添加椒鹽噪聲和高斯噪聲，分別采用線性的均值濾波函數(shù)imfilter和非線性的中值濾波函數(shù)medfilt2濾波進(jìn)行處理，兩種濾波的掩模都分別嘗試定義為3*3和7*7。由得到的實(shí)驗(yàn)結(jié)果，分析哪種濾波對這種噪聲處理的效果比較好？掩模大小對噪聲處理效果有什么影響？f=imread(headCT.tif);g=

2、imnoise(f,salt & pepper);t=imnoise(f,gaussian);w1=ones(3);w2=ones(7);g1=imfilter(g,w1);g2=imfilter(g,w2);g3=medfilt2(g,3 3);g4=medfilt2(g,7 7);subplot(1,5,1);imshow(g);subplot(1,5,2);imshow(g1);subplot(1,5,3);imshow(g2);subplot(1,5,4);imshow(g3);subplot(1,5,5);imshow(g4); 椒鹽噪聲處理圖：t1=imfilter(t,w1);t

3、2=imfilter(t,w2);t3=medfilt2(g,3 3);t4=medfilt2(g,7 7);subplot(1,5,1);imshow(t);subplot(1,5,2);imshow(t1);subplot(1,5,3);imshow(t2);subplot(1,5,4);imshow(t3);subplot(1,5,5);imshow(t4);高斯噪聲處理圖：從對比可以看出非線性的中值濾波函數(shù)medfilt2濾波進(jìn)行處理的效果更好；掩模大小影響著圖片的平滑效果，掩模越大，平滑效果越好，但是容易造成邊緣信息缺失。2、 編寫m文件，實(shí)現(xiàn)同時(shí)用sobel算子對圖片的水平和垂直方

4、向進(jìn)行處理？function g=mysobel(a);%自己編寫m文件，實(shí)現(xiàn)同時(shí)用sobel算子對圖片的水平和垂直方向進(jìn)行處理M文件代碼：function g=mysobel(a);a1=im2double(a);f=rgb2gray(a1);gx=f;gy=f;m n=size(f);wx=fspecial(sobel);wy=wx;for i=2:m-1 for j=2:n-1 temp=double(f(i-1:i+1,j-1:j+1);gx(i,j)=abs(sum(sum(temp.*wx); gy(i,j)=abs(sum(sum(temp.*wy); endendg=gx+gy

5、;imshow(g);3、 對圖片text.tif進(jìn)行膨脹操作，選擇結(jié)構(gòu)元素為0 1 0;1 1 1;0 1 0和1 1 1;1 1 1;1 1 1觀察其處理效果的異同，如果一次膨脹的效果不明顯，可多次膨脹再比較；a=imread(text.tif);q=0 1 0;1 1 1;0 1 0;w=1 1 1;1 1 1;1 1 1;a1=imdilate(a,q);a2=imdilate(a1,q);a3=imdilate(a2,q);b1=imdilate(a,w);b2=imdilate(b1,w);b3=imdilate(b2,w);subplot(2,4,1);imshow(a);sub

6、plot(2,4,2);imshow(a1);subplot(2,4,3);imshow(a2);subplot(2,4,4);imshow(a3);subplot(2,4,6);imshow(b1);subplot(2,4,7);imshow(b2);subplot(2,4,8);imshow(b3);兩者膨化一次的效果較接近，在第二次膨化的時(shí)候前者效果更好，第三次膨化處理時(shí)前者變得比較模糊，后者無法完整識別。4、 選做題：對圖片shapes.tif進(jìn)行開、閉運(yùn)算，比較其處理的效果，采用方形結(jié)構(gòu)元素（可由strel（square,3）得到）。f=imread(shapes.tif);se=s

7、trel(square,3);f1=imopen(f,se);f2=imclose(f,se);f3=imclose(f1,se);subplot(2,2,1);imshow(f);title(原圖);subplot(2,2,2);imshow(f1);title(開運(yùn)算處理圖);subplot(2,2,3);imshow(f2);title(閉運(yùn)算處理圖);subplot(2,2,4);imshow(f3);title(先開再閉運(yùn)算處理圖);se=strel(square,3);se=strel(square,20);5、 選做題：仿照課本P272例9.7，計(jì)算和顯示圖像pattern.tif中最后一個(gè)字母”a”的質(zhì)心f=imread(pattern.tif);f=im2bw(f);imshow(f);L,num=bwlabel(f);figure,imshow(L,)g=zeros(size(f);g(L=1048)=1;figure,imsh