基于强度的分析方法

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,#,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,Page,#,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,1/30/2013,#,基于强度的分析方法,图象的数字化处理技术,基于条纹亮度,强度分析法,基于相位,相位分析法,(,主流,),前景(目标),图像中特定的、具有独特性质的区域,背景,图像的其他部分,目标特征,表征图像最基本的属性或特征，是目标识别和分析的基础,目标特征,边缘、灰度、形状、纹理、色彩等多个方面,chapter 1,边缘检测,图像边缘,图像最基本的特征,本质,是图像局部的不连续性,包括,灰度级的突变，颜色的突变，纹理结构的突变,组成,由一些灰度函数的导数值超过预先设定的阈值的,像素组成,在实际图像中，由于图像传感器的性能、成像过程中的噪声等因素的影响，理想的阶梯状和屋顶状边缘是很少见的，而是在灰度变化的上升和下降沿都比较缓慢，表现为斜坡状。,图像的边缘大致可以分为两类：,阶梯状,此类边缘位于图像,灰度存在差异的两个,区域之间,屋顶状,此类边缘是图像灰度突 然从一个值变化到另一个值，保持一个比较小的行程之后又变回原来的值,边缘检测,概念,查找图像的特征量急剧变化的位置的图像处理方法,研究目的,对于光电检测来说，很多时候干涉条纹是我们检测结果的 一个重要反映，那么对于干涉条纹的处理就显得尤为重要，直接影响着我们检测的准确性。,方法,变形曲线模型（轮廓法）,导数法（一阶、二阶）,基于变形曲线模型（,Deformable Model,）,提取方法：从外向内逐步逼近,变形曲线模型,定义在图像上的曲线，在与曲线相关的内部力和与图像相关的外部力的共同作用下向目标移动,保持平滑性,吸引轮廓曲线向目标或感兴趣区域运动。,参数主动轮廓模型,(snake),，使用参数来描述曲线,几何主动轮廓模型，利用水平集的方法来实现,先验模型，对目标形状和特征的先验性假设,能量函数，反映这种先验知识以及曲线自身的描述,Snake,模型,基本理念,在图像目标附近定义一条具有能量的曲线，在,内部,能量和,外部,能量的共同作用下，不断变形以寻求与图像中对应的能量极小的位置。,本质,求解满足能量最小的一条曲线，在能量最小化过程中产生了内部势力和外部势力。,能量最小化问题都能转化为求解泛函极值的问题。,构造方法,利用物理的概念,外部能量和内部能量,形变,外部力和内部力的作用,可以用力的平衡原理来解释,Snake,模型的形变过程,实质上就是内部势力和外部势力趋于平衡,灰度分布、梯度以及边缘形状,构造,Snake,模型的能量函数,设置初始轮廓线,逼近目标边缘,Snake,模型边缘检测过程,图,为,Snake,模型的边缘检测过程,图,(a),：实线为目标边缘，虚线为初始轮廓，,图,(b),：能量最小化过程，轮廓曲线在外部力和内部力的共同作用,下，不断下目标边缘靠近,图,(c),：停止在目标边缘上,计算过程,在数学上,：,X(s)=(x(s),y(s),，,s,0,1,能量泛函：,，,不为零时，曲线是连续光滑的,当,(s),为零时，曲线不连续会出现断点,当,(s),为零时，曲线会出现角点，即曲率断点，而,值很大时，,E(X),的最小值对应的闭合曲线是一个圆，对应的非闭合曲线是一条直线。,x(s),和,y(s),表示曲线在图像中的坐标位置。,弹力系数,强度系数,代表曲线的弹性能量,阻止轮廓曲线伸展,代表曲线的刚性能量,阻止轮廓曲线弯曲,来自曲线的内部能量，确定了轮廓的连续性和光滑性，表示了轮廓拉伸和弯曲程度，其取值与图像无关,虽然较大的,会使图像的边缘模糊，但为了扩大初始轮廓线的捕捉区域，适当的增加,是有必要的。,综上，能量泛函可以表示为：,权重系数，为具有标准差是,的二维高斯函数,外部能量，基于图像数据决定轮廓的区域范围。它吸引曲线向目标运动,计算过程就是求解上式极小值的过程，初始轮廓线在内部能量和外部能量的作用下不断逼近目标轮廓，两种能量均衡的结果就是所求的目标边缘。,为使总能量最小，曲线,X(s),应满足欧拉方程：,即：,上式可以看作是一个力的平衡方程：,Fint+Fext=0,其中，内部力为：,外部力为：,将,X(s),看作时间,t,的函数,X(s,t),，欧拉方程式将变为：,能量最小化的过程就是将初始轮廓放在图像空间，按式欧拉方程进行变形，当上式的解趋于平稳时，轮廓线将收敛到目标边缘。,读取图像,图像预处理,设置初始轮廓线及参数,结束,计算能量函数,控制点调整,能量变化是否小于,5%,？,是,否,主动轮廓模型收敛过程,框架示意图,导数法,基本理念,边缘,：,由灰度级和邻域点不同的像素构成,是灰度不连续,的反映,若想检测边缘就应该突出相邻的灰度级的变化,微分运算就成为图像边缘清晰的重要工具,基本思想,1.,利用边缘增强算子,突出条纹图像中的局部边缘,2.,定义像素的边缘强度,通过设置阈值的方法提取边缘点集,边缘处,一阶导数,存在一个阶跃,可以用一阶导数的幅度值来检测边缘的存在,幅度的峰值一般对应边缘的位置,二阶导数,有一个向上的脉冲和一个向下的脉冲，两个脉冲之间有一个过零点，对应边缘位置,可以用二阶导数的过零点检测边缘的位置,二阶导数在过零点附近的符号确定图像边缘两侧的明区和暗区,对于干涉图样来说，如要对其进行强度分析，可以利用差分近似微分得到，也就是需要一些空域微分算子。,sobel,算子,改进的,Laplacian,算子,Sobel,算子,导数算子具有突出灰度变化的作用,对条纹图像运用导数算子,灰度变化较大的点处算得的值较高。,基于一阶微分,估计条纹图像灰度变化梯度方向,增强条纹的这些变化区域,对其设,定阈值,不是边缘点,判为边缘点,判断梯度模值是否大于阈值,小于,大于,二元图像函数,f(x,y),的,梯度函数,是矢量：,梯度值,大小：,梯度方向,：,梯度方向为变化率最大方向,替代依据,对于干涉条纹图像来说，是由,CCD,采集的数字,图像,是离散的,可以用一阶差分直接代替条纹,图像的偏导数,f(i,j),梯度算子,f,i,(i,j)=f(i,j)-f(i-1,j),f,i,(i,j)=f(i+1,j)-f(i,j),f,i,(i,j)=,f,j,(i,j)=f(i,j)-f(i,j-1),f,j,(i,j)=f(i,j+1)-f(i,j),f,j,(i,j)=,f,i,(i,j),f,j,(i,j),一般来说，任意正交方向都可以定义出对应的梯度：,比如,=45,时，有,f,45,(i,j)=,f,135,(i,j)=,构造出一个模板（单位化）,i,方向,j,方向,这就是,roberts,算子，也叫交叉算子,-1,0,0,1,0,-1,1,0,检测结果,Roberts,算子通过对角线方向上相邻的两个像素之差近似梯度幅值。计算,出来的梯度近似值位置相同，点位于内插点,i+1/2,j+1/2,处，即在,2,2,邻域的四个像素之间。,结论,该算子仅对噪声干扰小且边缘较为陡峭的图像有着较为理想的检测效果。,优势,：,提取,水平边缘,和,垂直边缘,时,效果好，定位精度高,劣势,：,对,斜边缘,的提取效果不理想，存在漏检。,不能有效的抑制,噪声,Sobel,算子,改善目的,为改善该算子中会增强边缘和噪音的特性,改善方法,上述差分式分别求出了灰度在,x,和,y,方向上的变化率,但是,要对每一个像素进行以上的运算,运算量较大,所以在实际,中采用小型模板利用卷积来近似计算,对,x,方向和,y,方向分,别使用一个模板。,Sobel,算子是在,33,邻域内计算,x,方向和,y,方向的偏导数,该方法赋予了上下左右四个像素点更大的权重。,f,x,=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),f,y,=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),Sobel,算子模板,i,方向,j,方向,条纹中的每个点都与图中的两个模板作卷积,第一个模板对水平边缘影响最大,;,第二模板对垂直边缘影响最大。,两个卷积的最大值作为该点的输出,运算结果是一幅边缘幅度图像。,-1,0,1,-2,0,2,-1,0,1,-1,-2,-1,0,0,0,1,2,1,总结,Sobel,边缘检测算子是一种非线性边缘算子,本质,上是通过计算一阶导数来检测边缘的,同时也可以给出边缘点的梯度方向。该算子在求梯度值前，先进行邻域的加权平均，再进行微分，它是边缘检测中最常用的算子。,优势,：,适合于直干涉条纹的处理,受噪声影响小,使用大的邻域时,抗干扰特性会更好,劣势,：,会使边缘较粗,既有大小又有方向,因此,数据存储量较大,Laplacian,算子,基本思想,图像边缘点除了对应于一阶微分幅度大的特点外,也对应于二阶微分的零交叉点,即就是在图像拐点位置处的二阶导数为零。,检测方法,通过寻找二阶导数的零交叉点来寻找边缘,基于二阶微分,二元图像函数,f(x,y),的,Laplacian,变换为：,用差分代替偏微分，可以在,i,方向和,j,方向上得到二次偏微分：,f,i,2,(i,j)=f,i,(i+,j)-f,i,(i,j),=f(i+1,j)-2f(i,j)+f(i-1,j),f,i,2,(i,j),f,j,2,(i,j)=f,j,(i,j+)-f,j,(i,j),=f(i,j+1)-2f(i,j)+f(i,j-1),f,j,2,(i,j),Laplacian,算子是一个标量,具有各向同向性，同时具有线性和位移不变性，其离散形式为,:,Laplacian,算子也是借助各种模板卷积实现的,(,邻域中心值具有较大权重,),缺陷,一阶导数对噪声敏感,因而不稳定,由此,二阶导数对噪声会更加敏感,因而会更不稳定。,0,-1,0,-1,4,-1,0,-1,0,-1,-1,-1,-1,8,-1,-1,-1,-1,Laplacian,算子改进,改进目的,在实际的干涉图中,不仅有有规律的灰度值分布,而且有噪声存在。由于求导运算起到了,噪声放大,的作用,因而这类方法效果并不好。,改进方法,改进的方法是先对图像进行适当的,平滑,以抑制,噪声,然后求导数。,平滑方法,对图像进行线性平滑,在数学上是进行一次卷积,运算，一般来说，我们用,Gauss,函数,来进行平滑。,这种将高斯滤波和拉普拉斯边缘检测结合在一起的方法就称为,LoG(Laplacian of Gauss),算子法。,对图像的卷积为：,p,为原来像素的灰度值,为权函数,P,为平滑后的灰度值,即求某个邻域中的灰度值的加权平均。,高斯脉冲函数和相应的高斯拉普拉斯算子,LoG,如下：,LOG,算子的函数,图形如墨西哥草帽，,也叫墨西哥草帽算子。,最后通过零交叉点的位置确定边缘点，边缘点的集合,P(x,y),可表示为：,P(x,y)=(x,y,)|f(x,y)*G(x,y,)=0,典型的模板为,：,5,5LOG,算子模板,0,0,-1,0,0,0,-1,-2,-1,0,-1,-2,16,-2,-1,0,-1,-2,-1,0,0,0,-1,0,0,-2,-4,-4,-4,-2,-4,0,8,0,-4,-4,8,24,8,-4,-4,0,8,0,-4,-2,-4,-4,-4,-2,两种等效计算方法,：,1.,图像与高斯函数卷积，再求卷积的拉普拉斯微分,2.,求高斯函数的拉普拉斯微分，再与图像卷积,实例对比,原始,sobel,算子,Laplacian,算子,LOG,算子,结果分析,Sobel,算子,对直干涉条纹有较好的效果,但使直条纹边缘变得较粗,使屋顶状条纹出现响应飘移不定的问题,Laplacian,算子,使两种干涉条纹变模糊,LoG,算子,克服了,Sobel,算子和,Laplaci