光信数字图像处理图像的边缘检测.pdf

实验三 图像的边缘检测 一、实验目的 1．进一步理解边沿检测的基本原理。 2 ．掌握对图像边沿检测的基本方法。 3．学习利用 Matlab 图像工具箱对图像进行边沿检测。 二、实验原理 边缘检测在图像处理中可分为空域下的操作和频域下的操作。本实验以空域下 的操作为主要实验内容，其他方法同学们可以在课后自行实验。空域下的图像 边缘检测主要应用的是领域运算。 三、实验要求 对边缘检测的要求：使用 Matlab 图像处理工具箱中的不同方法对下图的边缘 进行提取；注意观察不同操作对垂直方向、水平方向、斜方向的提取效果有何 区别；注意观察提取后的边界是否连续，若不连续可采用什么方法使其连续。 四、实验步骤 1．打开计算机，启动 MATLAB 程序； 2．调入实验指导书所示的数字图像，利用 MATLAB 图像工具箱中已有函数 进行编程以实现图像的边缘检测； 3．显示原图和处理过的图像。 4．记录和整理实验报告 五、实验仪器 1．计算机； 2．MATLAB 程序； 3．记录用的笔、纸。 六、实验报告要求 1．对所使用的图像边沿检测方法的原理进行说明。 2 ．给出程序清单和注释。 3．对不同处理方法进行说明、比较。 ******************** 自己的方法 *********************** QQ = imread('Image1.bmp'>。 blood=rgb2gray(QQ>。 [x,y]=size(blood> 。 % 求出图象大小 b=double(blood>。 N =sqrt(100> * randn(x,y> 。 % 生成方差为 10 的白噪声 I=b+N 。 % 噪声干扰图象 for i=1:x % 实际图象的灰度为 0～255 for j=1:y if (I(i,j>>255> I(i,j>=255 。 end if (I(i,j><0> I(i,j>=0 。 end end end z0=max(max(I>> 。 % 求出图象中最大的灰度 z1=min(min(I>> 。 % 最小的灰度 T=(z0+z1>/2 。 TT=0 。 S0=0。 n0=0。 S1=0。 n1=0。 allow=0.5 。 % 新旧阈值的允许接近程度 d=abs(T-TT> 。 count=0。 %记录几次循环 while(d>=allow> % 迭代最佳阈值分割算法 count=count+1。 for i=1:x for j=1:y if (I(i,j>>=T> S0=S0+I(i,j> 。 2 / 7 n0=n0+1。 end if (I(i,j><T> S1=S1+I(i,j> 。 n1=n1+1。 end end end T0=S0/n0。 T1=S1/n1。 TT=(T0+T1>/2 。 d=abs(T-TT> 。 T=TT 。 end Seg=zeros(x,y>。 for i=1:x for j=1:y if(I(i,j>>=T> Seg(i,j>=1。 % 阈值分割的图象 end end end SI=1-Seg。 %阈值分割后的图象求反，便于用腐蚀算法求边缘 se1=str