基于matlab边缘提取的几种方法的比较,基于matlab边缘检测与提取的几种方法的比较.doc...

基于matlab边缘检测与提取的几种方法的比较.doc

基于MATLAB边缘检测与提取的几种方法的比较数字图像边缘检测DIGITALIMAGEPROCESSING又称为计算机图像边缘检测，它是指将图像信号转换成数字信号并利用计算机对其进行处理的过程。由于图像边缘是图像最基本的特征之一，往往携带着一幅图像的大部分信息。而边缘存在于图像的不规则结构和不平稳现象中，也即存在于信号的突变点处，这些点给出了图像轮廓的位置，这些轮廓常常是我们在图像边缘检测时所需要的非常重要的一些特征条件，这就需要我们对一幅图像检测并提取出它的边缘。在通常情况下，我们可以将信号中的奇异点和突变点认为是图像中的边缘点，其附近灰度的变化情况可从它相邻像素灰度分布的梯度来反映。根据这一特点，提出了多种边缘检测算子如ROBERT算子、SOBEL算子、PREWITT算子、LAPLACIAN算子，CANNY算子等。这些方法多是以待处理像素为中心的邻域作为进行灰度分析的基础，实现对图像边缘的提取并已经取得了较好的处理效果。经典的边界提取技术大都基于微分运算。首先通过平滑来滤除图像中的噪声，然后进行一阶微分或二阶微分运算，求得梯度最大值或二阶导数的过零点，最后选取适当的阈值来提取边界。本文主要介绍几种经典的边缘提取算法，选取两种用MATLAB语言编程实现，对提取结果进行比较和分析。图像边缘检测的基本步骤(1)滤波。边缘检测主要基于导数计算，但受噪声影响。但滤波器在降低噪声的同时也导致边缘强度的损失。(2)增强。增强算法将邻域中灰度有显著变化的点突出显示。一般通过计算梯度幅值完成。(3)检测。但在有些图像中梯度幅值较大的并不是边缘点。最简单的边缘检测是梯度幅值阈值判定。(4)定位。精确确定边缘的位置。几种边缘算子的比较以柚子的图片为例1、ROBERTS算子是一种利用局部差分算子寻找边缘的算子，ROBERTS算子边缘定位准，但是对噪声敏感。适用于边缘明显而且噪声较少的图像分割，在应用中经常用ROBERTS算子来提取道路。程序如下IIMREAD D\研一\阳建宏信号\作业\柚子JPG 读取图像JRGB2GRAYI转化为灰度图像KIMADJUSTJ,40/2551调整灰度值BWEDGEK, ROBERTS 边缘检测FIGURE1IMSHOWBW显示图像2、SOBEL边缘算子图像中的每个像素都用这核做卷积。这两个核分别对垂直边缘和水平边缘响应最大，两个卷积的最大值作为该点的输出位。运算结果是一幅边缘幅度图像。SOBEL算子认为邻域的像素对当前像素产生的影响不是等价的，所以距离不同的像素具有不同的权值，对算子结果产生的影响也不同。一般来说，距离越大，产生的影响越小。程序如下IIMREAD D\研一\阳建宏信号\作业\柚子JPG 读取图像JRGB2GRAYI转化为灰度图像KIMADJUSTJ,40/2551调整灰度值BWEDGEK, SOBEL 边缘检测FIGURE1IMSHOWBW显示图像3、PREWITT算子在一个方向求微分，而在另一个方向求平均，因而对噪声相对不敏感，有抑制噪声作用。但是像素平均相当于对图像的低通滤波，所以PREWITT算子对边缘的定位不如ROBERTS算子。程序如下IIMREAD D\研一\阳建宏信号\作业\柚子JPG 读取图像JRGB2GRAYI转化为灰度图像KIMADJUSTJ,40/2551调整灰度值BWEDGEK, PREWITT 边缘检测FIGURE2IMSHOWBW显示图像这是未经滤波的效果图。4、LAPLACIAN算子利用二阶导数信息，具有各向同性，即与坐标轴方向无关，坐标轴旋转后梯度结果不变。使得图像经过二阶微分后，在边缘处产生一个陡峭的零交叉点，根据这个对零交叉点判断边缘。LAPLACIAN算子对噪声比较敏感，LAPLACIAN算子有一个缺点是它对图像中的某些边缘产生双重响应。所以图像一般先经过平滑处理，通常把LAPLACIAN算子和平滑算子结合起来生成一个新的模板。5、LOG(LAPLACIANOFGASSIAN)算法将高斯滤波和拉普拉斯检测算子结合在一起进行边缘检测的方法。也称之为拉普拉斯高斯算法。该算法的主要思路和步骤是滤波(滤波函数根据人类视觉特性选为高斯函数)、增强(对平滑图像进行拉普拉斯运算)、检测(边缘检测判据是二阶导数的零交叉点并对应一阶导数的较大峰值)。这种方法的特点是图像首先与高斯滤波器进行卷积，这样既平滑了图像又降低了噪声，孤立的噪声点和较小的结构组织将被滤除。但是由于平滑会造成图像边缘的延伸，因此边缘检测器只考虑那些具有局部梯度最大值的点为边缘点。高斯拉普拉斯算子把高斯平滑滤波器和拉普拉斯锐化滤波器结合起来，先平化掉噪声，再进行边缘检测，所以效果更好。程序如下IIMREAD D\研一\阳建宏信号\作业\柚子JPG JRGB2GRAYI转化为灰度图像KIMADJUSTJ,40/2551调整灰度值BW1EDGEK, LOG ,001FIGURE1SUBPLOT3,2,1IMSHOWBW1TITLE 阈值为001的LOG算子边缘检测图像 BW11EDGEK, LOG ,0001SUBPLOT3,2,5IMSHOWBW11TITLE 阈值为0001的LOG算子边缘检测图像 BW2EDGEK, LOG ,0005SUBPLOT3,2,4IMSHOWBW2TITLE 阈值为0005的LOG算子边缘检测图像 BW20EDGEK, LOG ,0007SUBPLOT3,2,3IMSHOWBW20TITLE 阈值为0007的LOG算子边缘检测图像 BW22EDGEK, LOG ,0009SUBPLOT3,2,2IMSHOWBW22TITLE 阈值为0009的LOG算子边缘检测图像 BW23EDGEK, LOG ,00001SUBPLOT3,2,6IMSHOWBW23TITLE 阈值为00001的LOG算子边缘检测图像 结果分析如图所示采用的阈值越小则图像的边缘处理效果越清晰，且边缘点条理显著。但是当阈值超过00001时图像边缘的处理效果越模糊如图520所示，但是仍能进行边缘的判断6、CANNY算子类似与MARR(LOG)边缘检测方法，也属于先平滑后求导数的方法。采用二维高斯函数的任意方向上的一阶方向导数为噪声滤波器，通过与图像卷积进行滤波；然后对滤波后的图像寻找图像梯度的局部最大值，以此来确定图像边缘。根据对信噪比与定位乘积进行测度，得到最优化逼近算子。这就是CANNY边缘检测算子。用双阈值算法检测和连接边缘。检测方法的优点①低误码率，很少把边缘点误认为非边缘点；②高定位精度，即精确地把边缘点定位在灰度变化最大的像素上；③抑制虚假边缘。程序如下IIMREAD D\研一\阳建宏信号\作业\柚子JPG JRGB2GRAYI转化为灰度图像KIMADJUSTJ,40/2551调整灰度值BW1EDGEK, CANNY ,001FIGURE2SUBPLOT3,2,1IMSHOWBW1TITLE 阈值为001的CANNY算子边缘检测图像 BW11EDGEK, CANNY ,0009SUBPLOT3,2,2IMSHOWBW11TITLE 阈值为0009的CANNY算子边缘检测图像 BW2EDGEK, CANNY ,0007SUBPLOT3,2,3IMSHOWBW2TITL