基于强度的分析方法精品课件

基于强度的分析方法

22021/4/1732021/4/1742021/4/17chapter1边缘检测52021/4/17

图像边缘图像最基本的特征

本质是图像局部的不连续性

包括灰度级的突变，颜色的突变，纹理结构的突变

组成由一些灰度函数的导数值超过预先设定的阈值的

像素组成62021/4/17

在实际图像中，由于图像传感器的性能、成像过程中的噪声等因素的影响，理想的阶梯状和屋顶状边缘是很少见的，而是在灰度变化的上升和下降沿都比较缓慢，表现为斜坡状。图像的边缘大致可以分为两类：阶梯状此类边缘位于图像灰度存在差异的两个区域之间屋顶状此类边缘是图像灰度突然从一个值变化到另一个值，保持一个比较小的行程之后又变回原来的值72021/4/17边缘检测概念查找图像的特征量急剧变化的位置的图像处理方法研究目的

对于光电检测来说，很多时候干涉条纹是我们检测结果的一个重要反映，那么对于干涉条纹的处理就显得尤为重要，直接影响着我们检测的准确性。方法变形曲线模型（轮廓法）导数法（一阶、二阶）

82021/4/17基于变形曲线模型（DeformableModel）提取方法：从外向内逐步逼近变形曲线模型定义在图像上的曲线，在与曲线相关的内部力和与图像相关的外部力的共同作用下向目标移动保持平滑性吸引轮廓曲线向目标或感兴趣区域运动。参数主动轮廓模型(snake)，使用参数来描述曲线几何主动轮廓模型，利用水平集的方法来实现先验模型，对目标形状和特征的先验性假设能量函数，反映这种先验知识以及曲线自身的描述92021/4/17

Snake模型

基本理念

在图像目标附近定义一条具有能量的曲线，在内部能量和外部能量的共同作用下，不断变形以寻求与图像中对应的能量极小的位置。

本质

求解满足能量最小的一条曲线，在能量最小化过程中产生了内部势力和外部势力。

能量最小化问题都能转化为求解泛函极值的问题。102021/4/179、人的价值，在招收诱惑的一瞬间被决定。2023/2/32023/2/3Friday,February3,202310、低头要有勇气，抬头要有低气。2023/2/32023/2/32023/2/32/3/20234:34:36PM11、人总是珍惜为得到。2023/2/32023/2/32023/2/3Feb-2303-Feb-2312、人乱于心，不宽余请。2023/2/32023/2/32023/2/3Friday,February3,202313、生气是拿别人做错的事来惩罚自己。2023/2/32023/2/32023/2/32023/2/32/3/202314、抱最大的希望，作最大的努力。03二月20232023/2/32023/2/32023/2/315、一个人炫耀什么，说明他内心缺少什么。。二月232023/2/32023/2/32023/2/32/3/202316、业余生活要有意义，不要越轨。2023/2/32023/2/303February202317、一个人即使已登上顶峰，也仍要自强不息。2023/2/32023/2/32023/2/32023/2/3

构造方法利用物理的概念

外部能量和内部能量

形变外部力和内部力的作用

可以用力的平衡原理来解释Snake模型的形变过程

实质上就是内部势力和外部势力趋于平衡灰度分布、梯度以及边缘形状构造

Snake模型的能量函数设置初始轮廓线逼近目标边缘122021/4/17Snake模型边缘检测过程图为Snake模型的边缘检测过程图(a)：实线为目标边缘，虚线为初始轮廓，图(b)：能量最小化过程，轮廓曲线在外部力和内部力的共同作用

下，不断下目标边缘靠近图(c)：停止在目标边缘上132021/4/17计算过程

在数学上：X(s)=(x(s),y(s))，s∈[0,1]

能量泛函：

α，β不为零时，曲线是连续光滑的当α(s)为零时，曲线不连续会出现断点当β(s)为零时，曲线会出现角点，即曲率断点，而β值很大时，E(X)的最小值对应的闭合曲线是一个圆，对应的非闭合曲线是一条直线。x(s)和

y(s)表示曲线在图像中的坐标位置。弹力系数强度系数代表曲线的弹性能量阻止轮廓曲线伸展代表曲线的刚性能量阻止轮廓曲线弯曲来自曲线的内部能量，确定了轮廓的连续性和光滑性，表示了轮廓拉伸和弯曲程度，其取值与图像无关142021/4/17

虽然较大的σ会使图像的边缘模糊，但为了扩大初始轮廓线的捕捉区域，适当的增加σ是有必要的。综上，能量泛函可以表示为：权重系数，为具有标准差是σ的二维高斯函数外部能量，基于图像数据决定轮廓的区域范围。它吸引曲线向目标运动152021/4/17

计算过程就是求解上式极小值的过程，初始轮廓线在内部能量和外部能量的作用下不断逼近目标轮廓，两种能量均衡的结果就是所求的目标边缘。

为使总能量最小，曲线X(s)应满足欧拉方程：即：上式可以看作是一个力的平衡方程：Fint+Fext=0162021/4/17其中，内部力为：外部力为：将X(s)看作时间t的函数X(s,t)，欧拉方程式将变为：

能量最小化的过程就是将初始轮廓放在图像空间，按式欧拉方程进行变形，当上式的解趋于平稳时，轮廓线将收敛到目标边缘。172021/4/17

读取图像图像预处理设置初始轮廓线及参数结束计算能量函数控制点调整能量变化是否小于5%？是否主动轮廓模型收敛过程框架示意图182021/4/17导数法

基本理念

边缘：由灰度级和邻域点不同的像素构成,是灰度不连续

的反映

若想检测边缘就应该突出相邻的灰度级的变化

微分运算就成为图像边缘清晰的重要工具

基本思想1.利用边缘增强算子,突出条纹图像中的局部边缘2.定义像素的边缘强度,通过设置阈值的方法提取边缘点集192021/4/17

边缘处一阶导数存在一个阶跃可以用一阶导数的幅度值来检测边缘的存在幅度的峰值一般对应边缘的位置二阶导数有一个向上的脉冲和一个向下的脉冲，两个脉冲之间有一个过零点，对应边缘位置可以用二阶导数的过零点检测边缘的位置二阶导数在过零点附近的符号确定图像边缘两侧的明区和暗区202021/4/17212021/4/17

对于干涉图样来说，如要对其进行强度分析，可以利用差分近似微分得到，也就是需要一些空域微分算子。sobel算子改进的Laplacian算子222021/4/17Sobel算子

导数算子具有突出灰度变化的作用,对条纹图像运用导数算子,灰度变化较大的点处算得的值较高。基于一阶微分估计条纹图像灰度变化梯度方向增强条纹的这些变化区域对其设定阈值不是边缘点判为边缘点判断梯度模值是否大于阈值小于大于232021/4/17二元图像函数f(x,y)的梯度函数是矢量：梯度值大小：梯度方向：梯度方向为变化率最大方向242021/4/17

替代依据对于干涉条纹图像来说，是由CCD采集的数字

图像,是离散的,可以用一阶差分直接代替条纹

图像的偏导数f(i,j)梯度算子fi(i,j)=f(i,j)-f(i-1,j)fi(i,j)=f(i+1,j)-f(i,j)fi(i,j)=fj(i,j)=f(i,j)-f(i,j-1)fj(i,j)=f(i,j+1)-f(i,j)fj(i,j)=fi(i,j)fj(i,j)252021/4/17一般来说，任意正交方向都可以定义出对应的梯度：比如α=45°时，有

f45°(i,j)=

f135°(i,j)=262021/4/17构造出一个模板（单位化）i方向j方向这就是roberts算子，也叫交叉算子-10010-110272021/4/17

检测结果

Roberts算子通过对角线方向上相邻的两个像素之差近似梯度幅值。计算出来的梯度近似值位置相同，点位于内插点[i+1/2,j+1/2]处，即在2×2邻域的四个像素之间。结论

该算子仅对噪声干扰小且边缘较为陡峭的图像有着较为理想的检测效果。优势：提取水平边缘和垂直边缘时效果好，定位精度高劣势：对斜边缘的提取效果不理想，存在漏检。不能有效的抑制噪声282021/4/17Sobel算子

改善目的为改善该算子中会增强边缘和噪音的特性

改善方法上述差分式分别求出了灰度在x和y方向上的变化率,但是

要对每一个像素进行以上的运算,运算量较大,所以在实际

中采用小型模板利用卷积来近似计算,对x方向和y方向分

别使用一个模板。Sobel算子是在3×3邻域内计算x方向和y方向的偏导数该方法赋予了上下左右四个像素点更大的权重。fx=[f(x+1,y-1)+2f(x+1,y)+f(x+1,y+1)]-[f(x-1,y-1)+2f(x-1,y)+f(x-1,y+1)]fy=[f(x-1,y+1)+2f(x,y+1)+f(x+1,y+1)]-[f(x-1,y-1)+2f(x,y)+f(x+1,y-1)]292021/4/17Sobel算子模板i方向j方向

条纹中的每个点都与图中的两个模板作卷积

第一个模板对水平边缘影响最大;第二模板对垂直边缘影响最大。两个卷积的最大值作为该点的输出,运算结果是一幅边缘幅度图像。-101-202-101-1-2-1000121302021/4/17

总结

Sobel边缘检测算子是一种非线性边缘算子

本质上是通过计算一阶导数来检测边缘的,同时也可以给出边缘点的梯度方向。该算子在求梯度值前，先进行邻域的加权平均，再进行微分，它是边缘检测中最常用的算子。优势：适合于直干涉条纹的处理受噪声影响小使用大的邻域时,抗干扰特性会更好劣势：会使边缘较粗既有大小又有方向,因此,数据存储量较大312021/4/17Laplacian算子

基本思想

图像边缘点除了对应于一阶微分幅度大的特点外,也对应于二阶微分的零交叉点,即就是在图像拐点位置处的二阶导数为零。

检测方法

通过寻找二阶导数的零交叉点来寻找边缘基于二阶微分322021/4/17二元图像函数f(x,y)的Laplacian变换为：用差分代替偏微分，可以在i方向和j方向上得到二次偏微分：fi2(i,j)=fi(i+,j)-fi(i,j)

=f(i+1,j)-2f(i,j)+f(i-1,j)fi2(i,j)fj2(i,j)=fj(i,j+)-fj(i,j)

=f(i,j+1)-2f(i,j)+f(i,j-1)fj2(i,j)332021/4/17Laplacian算子是一个标量,具有各向同向性，同时具有线性和位移不变性，其离散形式为:Laplacian算子也是借助各种模板卷积实现的(邻域中心值具有较大权重)缺陷

一阶导数对噪声敏感,因而不稳定,由此,二阶导数对噪声会更加敏感,因而会更不稳定。0-10-14-10-10-1-1-1-18-1-1-1-1342021/4/17Laplacian算子改进

改进目的在实际的干涉图中,不仅有有规律的灰度值分布,

而且有噪声存在。由于求导运算起到了噪声放大

的作用,因而这类方法效果并不好。

改进方法改进的方法是先对图像进行适当的平滑,以抑制

噪声,然后求导数。

平滑方法对图像进行线性平滑,在数学上是进行一次卷积

运算，一般来说，我们用Gauss函数来进行平滑。

这种将高斯滤波和拉普拉斯边缘检测结合在一起的方法就称为LoG(LaplacianofGauss)算子法。352021/4/17对图像的卷积为：p为原来像素的灰度值

ω为权函数,P为平滑后的灰度值,即求某个邻域中的灰度值的加权平均。362021/4/17高斯脉冲函数和相应的高斯拉普拉斯算子LoG如下：LOG算子的函数图形如墨西哥草帽，也叫墨西哥草帽算子。372021/4/17

最后通过零交叉点的位置确定边缘点，边缘点的集合P(x,y)可表示为：P(x,y)={(x,y,δ)|[f(x,y)*G(x,y,δ)]=0}典型的模板为：5×5LOG算子模板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-2382021/4/17两种等效计算方法：1.图像与高斯函数卷积，再求卷积的拉普拉斯微分2.求高斯函数的拉普拉斯微分，再与图像卷积392021/4/17

实例对比原始sobel算子Laplacian算子LOG算子402021/4/17结果分析Sobel算子对直干涉条纹有较好的效果,但使直条纹边缘变得较粗,使屋顶状条纹出现响应飘移不定的问题Laplacian算子使两种干涉条纹变模糊LoG算子克服了Sobel算子和Laplacian算子的不足,使直条纹和屋顶形条纹都达到了清晰的边缘,可满足干涉条纹后续处理要求412021/4/17参考文献：[1]李文峰.图形图像处理与应用[M].北京:中国标准出版社.2006.[2]唐立群等.数字图像模式识别方法分析[M].哈尔滨:哈尔滨工程

大学出版社.2008.6.[3]彭真明等.光电图像处理及应用[M].成都:电子科技大学出版社.2008.3.[4]谢珊珊.基于