试验九 图像分割.pdf

试验九 图像分割 实验目的 1.掌握线检测 2.边缘检测 3.多种边缘检测器的比较 4.使用梯度的分水岭分割 5.控制标记符的分水岭分割 实验内容： 1. 线检测 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function B=pixeldup(A, m, n) if nargin2; error(At least two inputs are required); end if nargin==2; n=m; end u=1:size(A,1); m=round(m); u=u(ones(1,m), :); u=u(:); v = 1:size(A,1); n=round(n); v=v(ones(1,n), :); v=v(:); B=A(u,v); %%%%%%%%%%%%%%%%%%%%%%%%%%%% clear all clc f=imread(D:\ 图像库\DIP3E_CH10\Fig1005(a)(wirebond_mask).tif); figure,imshow(f); w=[2 -1 -1;-1 2 -1;-1 -1 2]; g=imfilter(double(f),w); figure,imshow(g,[])%fig 2 gtop=g(1:120,1:120); gtop=pixeldup(gtop, 4); figure,imshow(gtop,[])%fig 3 gbot=g(end-119:end,end-119:end); gbot=pixeldup(gbot,4); figure,imshow(gbot,[])%fig 4 g=abs(g); figure,imshow(g,[])%fig 5 T=max(g(:)); g=g=T; figure,imshow(g)%fig 6 %%%%%%%%%%%%%%%%%%%%%%%% 2. 边缘检测 %%%%%%%%%%%%%%%%%%%%%%%%%%% clear all clc f=imread(D:\ 图像库\DIP3E_CH10\Fig1016(a)(building_original).tif); figure,imshow(f); [gv,t]=edge(f,sobel,vertical); figure,imshow(gv); t gv=edge(f,sobel,0.15,vertical); figure,imshow(gv); gboth=edge(f,sobel,0.15); figure,imshow(gboth); w45=[-2 -1 0;-1 0 1;0 1 2] g45=imfilter(double(f),w45,replicate); T=0.3*max(abs(g45(:))); g45=g45=T; figure,imshow(g45); %%%%%%%%%%%%%%%%%%%%%%%% %Sobel，LoG 和Canny 边缘检测器的比较 clear all clc f=imread(D:\ 图像库\DIP3E_CH10\Fig1016(a)(building_original).tif); figure,imshow(f); [g_sobel_default,ts]=edge(f,sobel); figure,imshow(g_sobel_default); [g_log_default,tlog]=edge(f,log); figure,imshow(g_log_default); [g_canny_default,tc]=edge(f,canny); figure,imshow(g_canny_default); [g_sobel_best,ts]=edge(f,sobel,0.05); figure,imshow(g_sobel_best); [g_log_best,tlog]=edge(f,log,0.003,2.25); figure,imshow(g_log_best); [g_canny_best,tc]=edge(f,canny,[0.04 0.10],1.5); figure,imshow(g_canny_best); 4.使用梯度的分水岭分割 clear all clc f=imread(D:\ 图像库\DIP3E_CH10\Fig1057(a)(small_blobs-origin