人脸识别matlab

1、1目录目录第一章 摘要.21.1 前言 .21.2 显微图像的特点 .2第二章边缘检测.32.1 边缘检测的定义 .32.2 细胞图像边缘检测算法的研究内容.4第三章边缘检测算子.63.1 SOBEL算子.63.2 CANNY算子.8第四章基于 MATLAB 的实验结果与分析.114.1 实验程序 .114.2 实验结果 .124.3 结果分析 .13参考文献.142第一章第一章 摘要摘要根据血液细胞边缘附近的灰度分布梯度较大的特性，运用 Sobel 算子、Canny 算子对血细胞进行图像的边缘提取。实验结果表明，Canny 算子进行边缘提取算法对于血细胞图像边缘提取有很好的效果。1.11.1

2、 前言前言长期以来，细胞和组织形态特征的分析一直是研究生物细胞学和病理学的主要手段。观察和计算血细胞的体积、表面积 、细胞及细胞核的形态，仍然是许多血液疾病临床和病理的常规检查。目前图像分析已经成为细胞学定量分析研究的有力工具，已用于细胞分类，细胞计数，肿瘤细胞的识别等各个方面 。不仅大大减轻了人的劳动， 而且为临床诊断提供了定量、客观的依据。然而要得到好的识别效果，有赖于细胞图像分割及边缘提取的准确性，即边缘提取是细胞图像识别的第一步，也是关键的一步。 本文根据细胞图像的特点，运用Sobel算子、Canny算子对血细胞进行图像的边缘提取，比较两种边缘提取方法相比较，验证了Canny算子对细胞

3、图像边缘提取的优越性。1.2 显微图像的特点显微图像的特点显微图像是指将组织、胞制成染色切片或涂片在显微镜下所看到的图像。本文以缺铁性贫血的血液涂片为研究对象。缺铁性贫血是贫血中最常见的一种，是机体对铁的需要增加、摄入不足或丢失过多等造成体内铁的缺乏，影响血红蛋白的合成而致的贫血。其典型血象为低色素小细胞性贫血，血红蛋白的减少比红细胞的减少更为明显。在染色的血涂片上可见红细胞的体积一般较小，大小和形态不一， 由于每个红细胞内所含的血红蛋白很少，红细胞表现为中心浅染区较大。一幅染色显微血液图像中有红细胞，各种白细胞，血小板，同时还可能含有杂质，细胞形状不固定，细胞之间可能会有交叠，也可能因光照不

4、均而使亮度分布不均匀。 因此在作血液细胞图像识别时，要能区分出图像中的红细胞、白细胞、血小板及杂质和交叠的细胞。对于白细胞还应当能够正确分割出它的细胞核，此外由于缺铁性贫血的红细胞的中心浅染区较大，因此进行细胞边缘检测时需要正确找到红细胞的中心浅染区的边缘。3第二章第二章边缘检测边缘检测2.1 边缘检测的定义边缘检测的定义图像边缘是图像最基本的特征，边缘在图像分析中起着重要的用。所谓边缘（edge）是指图像局部特征的不连续性。灰度或结构信息的突变称为边缘，例如：灰度级的突变、颜色的突变、纹理结的突变。边缘是一个区域的结束，也是另一个区域的开始，利用该特可以分割图像。当人们看一个有边缘的物体时，

5、首先感觉到的便是边缘，如 一条理想的边缘应该具有如图 2.1(a) 所示模型的特性。每个像素都处在灰度级跃变的一个垂直的台阶上（例如图形中所示的水平线通过图像的灰度剖面图） 。而实际上，诸如图像采集系统的性能、采样频率和获得图像的照明条件等因素的影响，得到的边缘往往是模糊的，边缘被模拟成具有“斜坡面”的剖面，如图2.1(b) 所示，在这个模型中不再有细线（宽为一个像素的线条） ，而是出现了边缘的点包含斜坡中任意点的情况。由此可以看到：模糊的边缘使边缘的“宽度”较大，面清晰的边缘使由此可以看到：模糊的边缘使边缘的“宽度”较大，面清晰的边缘使边缘的宽度较小。 图像的边缘有方向的幅度两个属性，沿边缘

6、方向像素变化平缓垂直于边缘方4向像素变化剧烈。边缘上的这种变化可以用微分算子检测出来，通常用一阶导数或二阶导数来检测边缘，认为最大值对应边缘位置，而二阶导数以过零点对应边缘位置。实际不同的是一阶导数上，对于图像中的任意方向上的边缘都可以进行类似的分析。图像边缘检测中对任意点的一阶导数可以利用该点梯度的幅度来获得，二阶导数可以用拉普拉斯算子得到。2.2 细胞图像边缘检测算法的研究内容细胞图像边缘检测算法的研究内容细胞图像边缘检测和分析可定义为应用一系列方法获取、校正、增强、变换、检测或压缩可视图像的技术。其目的是提高信息的相对质量，以便提取有用信息。图像边缘检测中的变换属于图像输入图像输出模式，

7、图像边缘检测是一种超越具体应用的过程，任何为解决某一特殊问题而开发的图像边缘检测新技术或新方法，几乎肯定都能找到其他完全不同的应用领域。细胞图像边缘检测的主要研究内容包括：（1）图像获得和抽样，其中通过人眼观察的视野获取图像的问题有：最常用的图像获取装置电视（TV）摄像机问题，对所获得信号进行独立的采样和数字化就可用数字形式表达景物中全部彩色内容；电荷耦合装置，用作图像传感器，对景物每次扫描一行，或通过平行扫描获得图像；选择正确的分辨力或采样密度，一幅图像实质上是二维空间中的信号，所以适用于信号处理的法则同样适用于图像边缘检测，在放射学中常常需要高分辨力，要求图像至少达到 2048 像素204

8、8 像素；灰度量化，图像强度也必须进行数字化，通常以 256级（按 1 字节编码）覆盖整个灰度，一般一幅灰度分辨力为 8 位，空间分辨力为 512 像素512 像素的图像需 0.25 兆字节的存贮容量。（2）图像分割，目的是把一个图像分解成它的构成成分，以便对每一目标进行测量。图像分割是一个十分困难的过程。但其测量结果的质量却极大地依赖于图像分割的质量。有两类不同的图像分割方法。一种方法是假设图像各成分的强度值是均匀的并利用这种均匀性；另一种方法寻找图像成分之间的边界，因而是利用图像的不均匀性。主要有直方图分割，区域生长，梯度法等。（3）边界查索，用于检测图像中线状局部结构，通常是作为图像分割

9、的一5个预处理步骤。大多数图像边缘检测技术应用某种形式的梯度算子，可应用对水平方向、垂直方向或对角线方向的梯度敏感的梯度算子，用它们的复合结果可检测任意方向的边界。（4）图像增强和复原，用于改进图像的质量。不同的增强技术可以用于不同的目的，这取决于应用的类型。如果打算直接观察图像，可以增强对比度。如果是为了进一步对图像作数字处理，可以选择分割（一种突出各图像成分之间的边界和线状结构的运算） 。该技术可以是整体的或局部的，也可以在某个频域或者空间域中进行。图像增强和复原的目的是为了提高图像的质量，如去除噪声，提高图像的清晰度等。图像增强不考虑图像降质的原因，突出图像中所感兴趣的部分。 （5）图像

10、分类（识别）图像分类（识别）属于模式识别的范畴，其主要内容是图像经过某些预处理（增强、复原、压缩）后，进行图像分割和特征提取，从而进行判决分类。图像分类常采用经典的模式识别方法，有统计模式分类和句法（结构）模式分类，近年来新发展起来的模糊模式识别和人工神经网络模式分类在图像识别中也越来越受到重视。 （6）图像变换：由于图像阵列很大，直接在空间域中进行处理，涉及计算量很大。因此，往往采用各种图像变换的方法，如傅立叶变换、沃尔什变换、离散余弦变换等间接处理技术，将空间域的处理转换为变换域处理，不仅可减少计算量，而且可获得更有效的处理（如傅立叶变换可在频域中进行数字滤波处理） 。目前新兴研究的小波变

11、换在时域和频域中都具有良好的局部化特性，它在图像边缘检测中也有着广泛而有效的应用。 6第三章边缘检测算子第三章边缘检测算子3.1 Sobel 算子算子索贝尔算子（Sobel operator）是图像处理中的算子之一，主要用作边缘检测。在技术上，它是一离散性差分算子，用来运算图像亮度函数的梯度之近似值。在图像的任何一点使用此算子，将会产生对应的梯度矢量或是其法矢量。该算子包含两组 3x3 的矩阵，分别为横向及纵向，将之与图像作平面卷积，即可分别得出横向及纵向的亮度差分近似值。如果以 A 代表原始图像，Gx 及Gy 分别代表经横向及纵向边缘检测的图像，其公式如下图像的每一个像素的横向及纵向梯度近似

12、值可用以下的公式结合，来计算梯度的大小。在以上例子中，如果以上的角度 等于零，即代表图像该处拥有纵向边缘，左方较右方暗。Sobel 算子有两个，一个是检测水平边沿的；另一个是检测垂直平边沿的。7Sobel 算子另一种形式是各向同性 Sobel(IsotropicSobel)算子，也有两个，一个是检测水平边沿的，另一个是检测垂直平边沿的。各向同性 Sobel 算子比普通Sobel 算子的位置加权系数更为准确，在检测不同方向的边沿时梯度的幅度一致。由于 Sobel 算子是滤波算子的形式，用于提取边缘，可以利用快速卷积函数，简单有效，因此应用广泛。美中不足的是，Sobel 算子并没有将图像的主体与背

13、景严格地区分开来，换言之就是 Sobel 算子没有基于图像灰度进行处理，由于Sobel 算子没有严格地模拟人的视觉生理特征，所以提取的图像轮廓有时并不能令人满意。在观测一幅图像的时候，我们往往首先注意的是图像与背景不同的部分，正是这个部分将主体突出显示，基于该理论，我们给出了下面阈值化轮廓提取算法，该算法已在数学上证明当像素点满足正态分布时所求解是最优的。Sobel边缘算子的卷积和图 3.2 所示图像中的每个像素都用这两个核做卷积。这两个核分别对垂直边缘和水平边缘响应最大，两个卷积的最大值作为该点的输出位。运算结果是一幅边缘幅度图像。Sobel 算子认为邻域的像素对当前像素产生的影响不是等价的

14、，所以距离不同的像素具有不同的权值，对算子结果产生的影响也不同。一般来说，距离越大，产生的影响越小。第一章第一章 摘要摘要 11.1 前言 11.2 显微图像的特点 1第二章边缘检测第二章边缘检测 282.1 边缘检测的定义 22.2 细胞图像边缘检测算法的研究内容 3第三章边缘检测算子第三章边缘检测算子 53.1 SOBEL算子 53.2 CANNY算子 6第四章基于第四章基于 MATLAB 的实验结果与分析的实验结果与分析 104.1 实验程序 104.2 实验结果 114.3 结果分析 12参考文献参考文献 123.2 Canny 算子算子检测阶跃边缘的基本思想是在图像中找出具有局部最大

15、梯度幅值的像素点。图像边缘检测必须满足两个条件：一是必须能有效地抑制噪声；二是必须尽量精确确定边缘的位置。既要提高边缘检测算子对边缘的敏感性，同时也提高了对噪声的敏感。1.Canny 边缘检测基本原理： (1)具有既能滤去噪声又保持边缘特性的边缘检测最优滤波器，其采用一阶微分滤波器。采用二维高斯函数的任意方向上的一阶方向导数为噪声滤波器，通过与图像卷积进行滤波；然后对滤波后的图像寻找图像梯度的局部最大值，以此来确定图像边缘。根据对信噪比与定位乘积进行测度，得到最优化逼近算子。这就是 Canny 边缘检测算子。(2)类似与 Marr（LOG）边缘检测方法，也属于先平滑后求导数的方法。2.Cann

16、y 边缘检测算法： step1:用高斯滤波器平滑图象；step2:用一阶偏导的有限差分来计算梯度的幅值和方向；step3:对梯度幅值进行非极大值抑制；step4:用双阈值算法检测和连接边缘。 其数学描述如下：step1: 二维为高斯函数为：9在某一方向 n 上是 G（x,y）的一阶方向导数为：式中：n 式方向矢量，G 是梯度矢量。将图像 f(x,y)与 Gn 作卷积，同时改变 n 的方向，Gn* f(x,y)取得最大值时的n 就是正交于检测边缘的方向。. 式中：n 式方向矢量，G 是梯度矢量。将图像 f(x,y)与 Gn 作卷积，同时改变 n 的方向，Gn* f(x,y)取得最大值时的n 就是

17、正交于检测边缘的方向。step2:A(x,y) 反映了图像(x,y)点处的边缘强度，是图像(x,y)点处的法向矢量。10step3:仅仅得到全局的梯度并不足以确定边缘，因此为确定边缘，必须保留局部梯度最大的点而抑制非极大值。 （non-MaxiMa suppression,NMS）解决方法：利用梯度的方向。四个扇区的标号为 0 到 3，对应 3*3 邻域的四种可能组合。在每一点上，邻域的中心像素 M 与沿着梯度线的两个像素相比。如果 M 的梯度值不比沿梯度线的两个相邻像素梯度值大，则令 M=0。step4:减少假边缘段数量的典型方法是对 G(x,y)使用一个阈值。将低于阈值的所有值赋零值。但问

18、题是如何选取阈值？解决方法：双阈值算法进行边缘判别和连接边缘。首先是边缘判别：凡是边缘强度大于高阈值的一定是边缘点；凡是边缘强度小于低阈值的一定不是边缘点；如果边缘强度大于低阈：值又小于高阈值，则看这个像素的邻接像素中有没有超过高阈值的边缘点，如果有，它就是边缘点，如果没有，它就不是边缘点。其次是连接边缘：双阈值算法对非极大值抑制图像作用两个阈值 1 和2，且 212，从而可以得到两个阈值边缘图像 G1(x,y)和 G2(x,y)由于G2(x,y)使用高阈值得到，因而含有很少的假边缘，但有间断(不闭合)。双阈值法要在G2(x,y)中把边缘连接成轮廓，当到达轮廓的端点时，该算法就在 G1(x,y

19、)的 8 邻点位置寻找可以连接到轮廓上的边缘，这样，算法不断地在 G1(x,y)中收集边缘，直到将 G1(x,y)连接起来为止。实际上，还有多种边缘点判别方法，如：将边缘的梯度分为四种：水平、11竖直、45 度方向、135 度方向。各个方向用不同的邻接像素进行比较，以决定局部极大值。若某个像素的灰度值与其梯度方向上前后两个像素的灰度值相比并不是最大的，那么将该像素置为零，即不是边缘。此外，在实际应用中，检测效果还与滤波模板大小有关，当 N=b2 +1 时有较好的检测效果。Canny 算子检测方法的优点：低误码率，很少把边缘点误认为非边缘点；高定位精度，即精确地把边缘点定位在灰度变化最大的像素上

20、；抑制虚假边缘。在这几种算法中除 Roberts 算子外都使用了图像模板，模板运算是图像的一种处理手段邻域处理，有许多图像增强效果都可以采用模板运算实现，如平滑效果，中值滤波，油画效果，图像的凹凸效果等等。在模板运算中，首先定义一个模板，模板的大小以 3*3 的较常见，也有 2*2, 5*5 或更大尺寸的。运算时，把模板中心对应到图像的每一个像素位置，然后按照模板对应的公式对中心像素和它周围的像素进行数学运算，算出的结果作为输出图像对应像素点的值。这些经典的边缘提取算子，虽然各自不同，有不同的长处，但是它们也有共同的特点每种算子对应的预定义的边缘是它们最适合使用的情形，也就是说它们有针对性。这一点在应用中是有优越性的，它们的针对性可以帮助我们完成特定的任务。同时这也是算子的局限性，对于一般的问题或者情况未知的问题，预定义边缘的方法可能不会达到最佳效果。第四章基于第四章基于 Matlab 的实验结果与分析的实验结果与分析4.1 实验程序实验程序qq=imread(xibao.jpg); figure(1) 12imshow(qq);title( 原图像);qq=r