Matlab-DIP(瓦)ch10图像分割练习题.docVIP

1 %% 图像分割 2  3 %% 点检测 4 clc 5 clear 6 f=imread(.\images\dipum_images_ch10\Fig1002(a)(test_pattern_with_single_pixel).tif); 7 subplot(121),imshow(f),title(点检测原图); 8  9 w=[-1,-1,-1; 10 -1,8,-1; 11 -1,-1,-1] 12 g=abs(imfilter(double(f),w)); 13 subplot(122),imshow(g),title(点检测结果); 14%点检测图结果如下：  15  16 %% 非线性滤波点检测 17 clc 18 clear 19 f=imread(.\images\dipum_images_ch10\Fig1002(a)(test_pattern_with_single_pixel).tif); 20 subplot(121),imshow(f),title(非线性滤波点检测原图); 21  22 m=3; 23 n=3; 24 %ordfilt2(f,order,domin)表示的是用矩阵domin中为1的地方进行排序，然后选择地order个位置的值代替f中的值，属于非线性滤波 25 g=imsubtract(ordfilt2(f,m*n,ones(m,n)),ordfilt2(f,1,ones(m,n)));%滤波邻域的最大值减去最小值 26 T=max(g(:)); 27 g2=g=T; 28 subplot(122),imshow(g2); 29 title(非线性滤波点检测后图); 30%运行结果如下：  31  32 %% 检测指定方向的线 33 clc 34 clear 35 f=imread(.\images\dipum_images_ch10\Fig1004(a)(wirebond_mask).tif); 36 subplot(321),imshow(f); 37 title(检测指定方向线的原始图像); 38  39 w=[2 -1 -1; 40 -1 2 -1; 41 -1 -1 2];%矩阵用逗号或者空格隔开的效果是一样的 42 g=imfilter(double(f),w); 43 subplot(322),imshow(g,[]); 44 title(使用-45度检测器处理后的图像); 45  46 gtop=g(1:120,1:120);%取g的左上角图 47 gtop=pixeldup(gtop,4);%扩大4*4倍的图 48 subplot(323),imshow(gtop,[]); 49 title(-45度检测后左上角放大图); 50  51 gbot=g(end-119:end,end-119:end);%取右下角图 52 gbot=pixeldup(gbot,4);%扩大16倍的图 53 subplot(324),imshow(gbot,[]); 54 title(-45度检测后右下角后放大图); 55  56 g=abs(g); 57 subplot(325),imshow(g,[]); 58 title(-45度检测后的绝对值图); 59  60 T=max(g(:)); 61 g=g=T; 62 subplot(326),imshow(g); 63 title(-45度检测后取绝对值最大的图) 64%检测指定方向的线过程如下：  65  66 %% sobel检测器检测边缘 67 clc 68 clear 69 f=imread(.\images\dipum_images_ch10\Fig1006(a)(building).tif); 70 subplot(321),imshow(f); 71 title(sobel检测的原始图像); 72  73 [gv,t]=edge(f,sobel,vertical);%斜线因为具有垂直分量，所以也能够被检测出来 74 subplot(322),imshow(gv); 75 title(sobel垂直方向检测后图像); 76  77 gv=edge(f,sobel,0.15,vertical); 78 subplot(323),imshow(gv); 79 title(sobel垂直检测0.15阈值后图像); 80  81 gboth=edge(f,sobel,0.15); 82 subplot(324),imshow(