一种基于可变窗口匹配技术的数字图像相关方法技术

本发明专利技术涉及图像处理领域，尤其涉及一种基于可变窗口匹配技术的数字图像相关方法，其采用可变搜索窗口，实现了可作用于整像素搜索的快速搜索方法，实现快速准确的整像素测量，大大提高了像素点匹配效率，缩短了计算时间。所述基于可变窗口匹配技术的数字图像相关方法，其不需要进行全场整像素搜索，只针对第一个点和亚像素搜索不成功的点进行整像素搜索，提高了计算速度，并有矫正搜索误差的能力。并有矫正搜索误差的能力。并有矫正搜索误差的能力。

【技术实现步骤摘要】
一种基于可变窗口匹配技术的数字图像相关方法


[0001]本专利技术涉及图像处理领域，尤其涉及一种基于可变窗口匹配技术的数字图像相关方法。

技术介绍

[0002]数字图像相关方法，即DIC技术是一种非接触，非干涉式的全场光学测量方法，被广泛应用与全场位移和应变测量中。在物体表面自然散斑或者人工喷涂散斑，当物体发生变形或位移时候，散斑也发生变化。通过机拍摄变形前后散斑图像，并比较相关性，计算出出物体全场位移以及应变。
[0003]数字图像相关方法通常长分为整像素搜索和亚像素搜索两部分，整像素搜索的常用有逐点搜索，人工标定，爬山搜索，人工鱼群搜索等算法。其中逐点搜索时一种全局搜索算法，准确度高，但耗时大；爬山搜索，人工鱼群等智能搜索容易陷入局部最优的误差；人工标定的方法操作不便。
[0004]数字图像相关方法的测量精度主要取决于亚像素测量精度，现在较为常用的亚像素测量方法是高斯牛顿反向组合算法(ICGN)。它是一种迭代算法，需要提供较为准确的初始变形猜测，通常是将整像素的搜索结果作为初始猜测进行亚像素搜素。

技术实现思路

[0005]为了解决上述现有技术中数字图像相关方法的缺陷，本专利技术提出了一种基于可变窗口匹配技术的数字图像相关方法，是一种可作用于整像素搜索的快速搜索方法，其采用可变搜索窗口，实现快速准确的整像素测量，大大提高了像素点匹配效率，缩短了计算时间。
[0006]本专利技术提出了一种基于可变窗口匹配技术的数字图像相关方法，首先获取参考图像上至少一个任意计算点在目标图像上的最佳匹配像素点；然后基于最佳匹配像素点计算目标图像与参考图像之间的变形参数取值；参考图像为试件变形前的拍摄图像，目标图像为试件变形后的拍摄图像；最佳匹配像素点的获取包括以下步骤：
[0007]SA1、令计算点在参考图像上的像素坐标记作(x0，y0)，计算点记作P(x0，y0)；获取参考图像上以P(x0，y0)为中心的第一设定区域作为参考子区R；获取目标图像上以像素坐标(x0，y0)为中心的第二设定区域作为搜索区域；
[0008]SA2、获取搜索区域内各个像素点与计算点P(x0，y0)的相关系数；令搜索区域以像素点Q(x1,y1)为中心的第一设定区域记作像素点Q(x1,y1)对应的目标子区；搜索区域内像素点Q(x1,y1)与计算点P(x0，y0)的相关系数，为Q(x1,y1)对应的目标子区与参考子区R的相关系数；
[0009]SA3、判断搜索区域内是否存在对应的相关系数大于设定的相关阈值的像素点；
[0010]SA4、否，则以像素坐标(x0，y0)为中心扩大搜索区域范围，以形成新的搜索区域，获取新的搜索区域内各个像素点与计算点P(x0，y0)的相关系数，然后返回步骤SA3；
[0011]SA5、是，则获取搜索区域内对应的相关系数最大的像素点作为计算点的最佳匹配像素点。
[0012]优选的，步骤SA2中目标子区与参考子区的相关系数的计算公式如下：
[0013][0014]其中，R(x,y)为参考子区R中像素坐标为(x,y)的像素点的灰度值，T(x',y')为目标子区T中像素坐标为(x
’
,y
’
)的像素点的灰度值，为参考子区R的平均灰度，为目标子区T的平均灰度；表示目标子区T中像素坐标为(x
’
,y
’
)的像素点在目标子区T上位置等同于参考子区R中像素坐标为(x,y)像素点在参考子区R上的位置；即目标子区T上像素点(x
’
,y
’
)与目标子区T中心点构成的矢量坐标等于参考子区R上像素点(x,y)与参考子区R中心点构成的矢量坐标。
[0015]优选的，步骤SA2还包括构建节点与搜索区域内像素点一一对应的Mx矩阵，SA2中将搜索区域内各像素点对应的相关系数填充到Mx矩阵中对应的节点上；步骤SA3具体为判断Mx矩阵中是否存在赋值大于设定的相关阈值的节点；否，则执行步骤SA4，SA4中还包括根据搜索区域扩充结果扩充Mx矩阵；是，则执行步骤SA5。
[0016]优选的，步骤SA4具体为：以像素坐标(x0，y0)为中心扩大搜索区域范围，以形成新的搜索区域；根据新的搜索区域扩充Mx矩阵，使得新的搜索区域中的像素点与扩充后的Mx矩阵的节点一一对应；遍历Mx矩阵中的节点，当节点赋值为初始值，则计算搜索区域中对应该节点的像素点与计算点P(x0，y0)的相关系数并填充到Mx矩阵中该节点上。
[0017]优选的，SA2中，Mx矩阵中各节点的初始值为
‑
1；SA4中，Mx矩阵中各新增节点的初始值为
‑
1。
[0018]优选的，包括以下步骤：
[0019]S1、拍摄试件的变形前图像和变形后图像，将变形前图像作为参考图像，将变形后图像作为目标图像；
[0020]S2、在参考图像上设置感兴趣区域，在感兴趣区域内筛选计算点，获取参考图像上以计算点为中心的参考子区的变形参数组其中u
x
表示参考子区的水平位移，v
y
表示参考子区的竖直位移，表示参考子区的水平拉伸，表示参考子区的水平剪切，表示参考子区的竖直剪切，表示参考子区的竖直拉伸；
[0021]对各参考子区的变形参数组进行迭代，获取迭代收敛时的参数组取值作为参考子区的变形参数取值；
[0022]部分计算点对应的变形参数组的初始赋值为{u0，v0，0，0，0，0}，u0为计算点和对应的最佳匹配像素点之间的水平位移，v0为计算点和对应的最佳匹配像素点之间的竖直位移；
[0023]另一部分计算点对应的变形参数组的初始赋值为
已计算的相邻计算点对应的变形参数取值。
[0024]优选的，步骤S2具体包括以下步骤：
[0025]S21、在参考图像上设置感兴趣区域，在感兴趣区域内筛选计算点；设置变形参数组令首个计算点作为检测点；
[0026]S23、获取所述检测点在目标图像上的最佳匹配像素点；结合所述检测点和对应的最佳匹配像素点计算水平位移u0和竖直位移v0；
[0027]S24、令检测点的变形参数组的初始赋值为{u0，v0，0，0，0，0}；然后采用ICGN算法对变形参数组进行迭代，获取迭代收敛时的参数组取值作为以检测点为中心点的参考子区的变形参数取值；
[0028]S25、如果感兴趣区域内还存在未计算变形参数取值的计算点，则将与检测点相邻的计算点作为新的检测点，令检测点的变形参数组的计算点作为新的检测点，令检测点的变形参数组的初始赋值为上一个检测点的变形参数取值；
[0029]S26、采用ICGN算法对检测点的变形参数组进行迭代；
[0030]如果变形参数组在设定的迭代次数内收敛，则获取迭代收敛时的参数组取值作为以检测点为中心点的参考子区的变形参数取值，然后返回步骤S25；
[0031]如果变形参数组在设定的迭代次数内未收敛，则令参考图像上的像素点(x2+u,y2+v)作本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于可变窗口匹配技术的数字图像相关方法，其特征在于，首先获取参考图像上至少一个任意计算点在目标图像上的最佳匹配像素点；然后基于最佳匹配像素点计算目标图像与参考图像之间的变形参数取值；参考图像为试件变形前的拍摄图像，目标图像为试件变形后的拍摄图像；最佳匹配像素点的获取包括以下步骤：SA1、令计算点在参考图像上的像素坐标记作(x0，y0)，计算点记作P(x0，y0)；获取参考图像上以P(x0，y0)为中心的第一设定区域作为参考子区R；获取目标图像上以像素坐标(x0，y0)为中心的第二设定区域作为搜索区域；SA2、获取搜索区域内各个像素点与计算点P(x0，y0)的相关系数；令搜索区域以像素点Q(x1,y1)为中心的第一设定区域记作像素点Q(x1,y1)对应的目标子区；搜索区域内像素点Q(x1,y1)与计算点P(x0，y0)的相关系数，为Q(x1,y1)对应的目标子区与参考子区R的相关系数；SA3、判断搜索区域内是否存在对应的相关系数大于设定的相关阈值的像素点；SA4、否，则以像素坐标(x0，y0)为中心扩大搜索区域范围，以形成新的搜索区域，获取新的搜索区域内各个像素点与计算点P(x0，y0)的相关系数，然后返回步骤SA3；SA5、是，则获取搜索区域内对应的相关系数最大的像素点作为计算点的最佳匹配像素点。2.如权利要求1所述的基于可变窗口匹配技术的数字图像相关方法，其特征在于，步骤SA2中目标子区与参考子区的相关系数的计算公式如下：其中，R(x,y)为参考子区R中像素坐标为(x,y)的像素点的灰度值，T(x',y')为目标子区T中像素坐标为(x
’
,y
’
)的像素点的灰度值，为参考子区R的平均灰度，T为目标子区T的平均灰度；表示目标子区T中像素坐标为(x
’
,y
’
)的像素点在目标子区T上位置等同于参考子区R中像素坐标为(x,y)像素点在参考子区R上的位置；即目标子区T上像素点(x
’
,y
’
)与目标子区T中心点构成的矢量坐标等于参考子区R上像素点(x,y)与参考子区R中心点构成的矢量坐标。3.如权利要求1所述的基于可变窗口匹配技术的数字图像相关方法，其特征在于，步骤SA2还包括构建节点与搜索区域内像素点一一对应的Mx矩阵，SA2中将搜索区域内各像素点对应的相关系数填充到Mx矩阵中对应的节点上；步骤SA3具体为判断Mx矩阵中是否存在赋值大于设定的相关阈值的节点；否，则执行步骤SA4，SA4中还包括根据搜索区域扩充结果扩充Mx矩阵；是，则执行步骤SA5。4.如权利要求3所述的基于可变窗口匹配技术的数字图像相关方法，其特征在于，步骤SA4具体为：以像素坐标(x0，y0)为中心扩大搜索区域范围，以形成新的搜索区域；根据新的搜索区域扩充Mx矩阵，使得新的搜索区域中的像素点与扩充后的Mx矩阵的节点一一对应；遍历Mx矩阵中的节点，当节点赋值为初始值，则计算搜索区域中对应该节点的像素点与计
算点P(x0，y0)的相关系数并填充到Mx矩阵中该节点上。5.如权利要求4所述的基于可变窗口匹配技术的数字图像相关方法，其特征在于，SA2中，Mx矩阵中各节点的初始值为
‑
...

【专利技术属性】
技术研发人员：张育中，杨乔飞，李娟，董敬涛，舒双宝，张腾达，
申请(专利权)人：合肥工业大学，
类型：发明
国别省市：