基于强度的分析方法.ppt

Page 1基于强度的分析方法.Page 2.Page 3.Page 4.Page 5.Page 6 图像边缘 图像最基本的特征 本质 是图像局部的不连续性 包括 灰度级的突变，颜色的突变，纹理结构的突变 组成 由一些灰度函数的导数值超过预先设定的阈值的 像素组成.Page 7 在实际图像中，由于图像传感器的性能、成像过程中的噪声等因素的影响，理想的阶梯状和屋顶状边缘是很少见的，而是在灰度变化的上升和下降沿都比较缓慢，表现为斜坡状。图像的边缘大致可以分为两类：.Page 8边缘检测边缘检测概念 查找图像的特征量急剧变化的位置的图像处理方法研究目的 对于光电检测来说，很多时候干涉条纹是我们检测结果的 一个重要反映，那么对于干涉条纹的处理就显得尤为重要，直接影响着我们检测的准确性。方法 变形曲线模型（轮廓法）导数法（一阶、二阶）.Page 9基于变形曲线模型（Deformable Model）.Page 10 Snake 模型 基本理念 在图像目标附近定义一条具有能量的曲线，在内部能量和外部能量的共同作用下，不断变形以寻求与图像中对应的能量极小的位置。本质 求解满足能量最小的一条曲线，在能量最小化过程中产生了内部势力和外部势力。能量最小化问题都能转化为求解泛函极值的问题。.Page 11 构造方法利用物理的概念 外部能量和内部能量 形变 外部力和内部力的作用 可以用力的平衡原理来解释 Snake 模型的形变过程 实质上就是内部势力和外部势力趋于平衡.Page 12Snake模型边缘检测过程图为 Snake 模型的边缘检测过程图(a)：实线为目标边缘，虚线为初始轮廓，图(b)：能量最小化过程，轮廓曲线在外部力和内部力的共同作用 下，不断下目标边缘靠近图(c)：停止在目标边缘上.Page 13计算过程 在数学上：X(s)=(x(s),y(s)，s0,1 能量泛函：，不为零时，曲线是连续光滑的 当(s)为零时，曲线不连续会出现断点 当(s)为零时，曲线会出现角点，即曲率断点，而 值很大时，E(X)的最小值对应的闭合曲线是一个圆，对应的非闭合曲线是一条直线。.Page 14 虽然较大的 会使图像的边缘模糊，但为了扩大初始轮廓线的捕捉区域，适当的增加 是有必要的。综上，能量泛函可以表示为：.Page 15 计算过程就是求解上式极小值的过程，初始轮廓线在内部能量和外部能量的作用下不断逼近目标轮廓，两种能量均衡的结果就是所求的目标边缘。为使总能量最小，曲线 X(s)应满足欧拉方程：即：上式可以看作是一个力的平衡方程：Fint+Fext=0.Page 16其中，内部力为：外部力为：将 X(s)看作时间 t 的函数 X(s,t)，欧拉方程式将变为：能量最小化的过程就是将初始轮廓放在图像空间，按式欧拉方程进行变形，当上式的解趋于平稳时，轮廓线将收敛到目标边缘。.Page 17 主动轮廓模型收敛过程框架示意图.Page 18导数法 基本理念 边缘：由灰度级和邻域点不同的像素构成,是灰度不连续 的反映 若想检测边缘就应该突出相邻的灰度级的变化 微分运算就成为图像边缘清晰的重要工具 基本思想 1.利用边缘增强算子,突出条纹图像中的局部边缘 2.定义像素的边缘强度,通过设置阈值的方法提取边缘点集.Page 19 边缘处.Page 20.Page 21 对于干涉图样来说，如要对其进行强度分析，可以利用差分近似微分得到，也就是需要一些空域微分算子。.Page 22Sobel算子 导数算子具有突出灰度变化的作用,对条纹图像运用导数算子,灰度变化较大的点处算得的值较高。.Page 23二元图像函数f(x,y)的梯度函数是矢量：梯度值大小：梯度方向：梯度方向为变化率最大方向.Page 24 替代依据对于干涉条纹图像来说，是由CCD采集的数字 图像,是离散的,可以用一阶差分直接代替条纹 图像的偏导数f(i,j)梯度算子.Page 25一般来说，任意正交方向都可以定义出对应的梯度：比如=45时，有 .Page 26构造出一个模板（单位化）i方向 j方向这就是roberts算子，也叫交叉算子-10010-110.Page 27 检测结果 Roberts 算子通过对角线方向上相邻的两个像素之差近似梯度幅值。计算出来的梯度近似值位置相同，点位于内插点i+1/2,j+1/2处，即在 22 邻域的四个像素之间。结论 该算子仅对噪声干扰小且边缘较为陡峭的图像有着较为理想的检测效果。.Page 28Sobel算子 改善目的 为改善该算子中会增强边缘和噪音的特性 改善方法 上述差分式分别求出了灰度在x 和y 方向上的变化率,但是 要对每一个像素进行以上的运算,运算量较大,所以在实际 中采用小型模板利用卷积来近似计算,对x 方向和y 方向分 别使用一个模板。Sobel 算子是在33邻域内计算x 方向和y 方向的偏导数该方法赋予了上下左右四个像素点更大的权重。.Page 29 Sobel 算子模板 i方向 j方向 条纹中的每个点都与图中的两个模板作卷积 第一个模板对水平边缘影响最大;第二模板对垂直边缘影响最大。两个卷积的最大值作为该点的输出,运算结果是一幅边缘幅度图像。-101-202-101-1-2-1000121.Page 30 总结 Sobel 边缘检测算子是一种非线性边缘算子 本质上是通过计算一阶导数来检测边缘的,同时也可以给出边缘点的梯度方向。该算子在求梯度值前，先进行邻域的加权平均，再进行微分，它是边缘检测中最常用的算子。.Page 31Laplacian 算子 基本思想 图像边缘点除了对应于一阶微分幅度大的特点外,也对应于二阶微分的零交叉点,即就是在图像拐点位置处的二阶导数为零。检测方法 通过寻找二阶导数的零交叉点来寻找边缘.Page 32二元图像函数f(x,y)的Laplacian 变换为：用差分代替偏微分，可以在i方向和j方向上得到二次偏微分：.Page 33 Laplacian 算子是一个标量,具有各向同向性，同时具有线性和位移不变性，其离散形式为:Laplacian 算子也是借助各种模板卷积实现的(邻域中心值具有较大权重)缺陷 一阶导数对噪声敏感,因而不稳定,由此,二阶导数对噪声会更加敏感,因而会更不稳定。0-10-14-10-10-1-1-1-18-1-1-1-1.Page 34Laplacian 算子改进 改进目的 在实际的干涉图中,不仅有有规律的灰度值分布,而且有噪声存在。由于求导运算起到了噪声放大噪声放大 的作用,因而这类方法效果并不好。改进方法 改进的方法是先对图像进行适当的平滑平滑,以抑制 噪声,然后求导数。平滑方法 对图像进行线性平滑,在数学上是进行一次卷积 运算，一般来说，我们用Gauss函数函数来进行平滑。这种将高斯滤波和拉普拉斯边缘检测结合在一起的方法就称为LoG(Laplacian of Gauss)算子法。.Page 35对图像的卷积为：p 为原来像素的灰度值 为权函数,P为平滑后的灰度值,即求某个邻域中的灰度值的加权平均。.Page 36高斯脉冲函数和相应的高斯拉普拉斯算子LoG如下：LOG算子的函数图形如墨西哥草帽，也叫墨西哥草帽算子。.Page 37 最后通过零交叉点的位置确定边缘点，边缘点的集合P(x,y)可表示为：P(x,y)=(x,y,)|f(x,y)*G(x,y,)=0典型的模板为：55LOG算子模板00-1000-1-2-10-1-216-2-10-1-2-1000-100-2-4-4-4-2-4080-4-48248-4-4080-4-2-4-4-4-2.Page 38两种等效计算方法：1.图像与高斯函数卷积，再求卷积的拉普拉斯微分 2.求高斯函数的拉普拉斯微分，再与图像卷积.Page 39 实例对比 原始 sobel算子 Laplacian算子 LOG算子.Page 40 结果分析.Page 41参考文献：1 李文峰.图形图像处理与应用M.北京:中国标准出版社.2006.2 唐立群等.数字图像模式识别方法分析M.哈尔滨:哈尔滨工程 大学出版社.2008.6.3彭真明等.光电图像处理及应用M.成都:电子科技大学出版社.2008.3.4谢珊珊.基于主动轮廓模型的图像边缘检测研究D.天津理工大学学位 论文.2011.1.5戴磊.轮廓检测的图像检索方法J.计算机仿真.2009.10.6周战荣等.干涉条纹边缘检测方法的比较与改进J.光学技术.2007.11.7吴东楼等.干涉条纹的处理方法研究J.光学学报.1999.1.8范志刚等.干涉条纹图的处理方法研究J.光学技术.2000.5.Page 42Thanks!.