CN112465755B - 一种初始子区细分方法、装置、计算机设备及存储介质 - Google Patents

Abstract

本发明公开了一种初始子区细分方法、装置、存储介质及计算机设备，包括采集被测件的当前状态图像与下一状态图像后，分别识别当前状态图像与下一状态图像上的标记点；依据标记点将当前状态图像、下一状态图像初步划分为若干子区，匹配第一类子区与第二类子区以形成多个子区对，并确定每个子区对的误差估计特征值P_k；比较误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；若子区对的误差估计特征值P_k大于特征值阈值，则将相应的子区对分别再次细化后继续判断对应的子区划分是否达标，兼顾了精度需求和处理速度，为后续各像点的快速准确初值估计提供了基础。

Description

一种初始子区细分方法、装置、计算机设备及存储介质

技术领域

本发明涉及散斑应变测量技术领域，具体涉及一种初始子区细分方法、装置、计算机设备及存储介质。

背景技术

三维数字图像相关技术(Three Digital Image Correlation,3D-DIC)，三维散斑应变测量方法由于具有非接触、非干涉、高精度、可三维建模、可全场测量等诸多优点，从其问世伊始，就逐渐从传统的材料力学测量扩展到建筑工程、生命科学、航空航天、工业检测等诸多领域。

长久以来，如何进一步提高散斑应变测量的精度、计算速度、测量的实时性，并保证测量结果的鲁棒性一直都是研究的热点。针对被测对象变形前后对应像点的匹配问题，学者们先后提出了多种不同的匹配算法。这些算法基本都符合先进行粗略的匹配初值估计，然后再用各种迭代算法进行精确匹配的流程。要完成快速初值估计，就需要将相机拍摄的变形前后的图像划分为若干个子区，这些子区内的像点具有相似的形变趋势，从而利用连续变形理论，通过像点的对应关系完成子区内所有像点的快速匹配。但是子区划分太过稀疏，则子区内初值估计难度越大，初值估计误差越大；子区划分太过密集，则需要精确匹配的初始点太多，计算速度降低。

发明内容

本发明的目的是提供一种初始子区细分方法，保证计算速度的前提下，控制子区的划分精度，从而提高初值估计精度，以解决现有技术存在的问题。

为实现上述目的，本发明提供一种初始子区细分方法，包括如下步骤：

S1采集被测件的当前状态图像与下一状态图像，所述被测件表面设有带标记点的散斑图案，分别识别所述当前状态图像与所述下一状态图像上的标记点；

S2依据所述标记点将所述当前状态图像、下一状态图像初步划分为若干子区，其中，所述当前状态图像的子区为第一类子区，所述下一状态图像的子区为第二类子区；

S3匹配所述第一类子区与所述第二类子区以形成多个子区对，并确定每个所述子区对的误差估计特征值P_k；

S4比较所述误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；若所述子区对的误差估计特征值P_k小于所述特征值阈值，则判断相应的子区对的划分已达标；若所述子区对的误差估计特征值P_k大于所述特征值阈值，则将相应的子区对分别再次细化后转入步骤S3,直至判断全部的所述子区对划分达标。

优选的，所述步骤S1中，对被测工件表面制斑后，进行标记点喷涂以获取所述被测件。

优选的，所述步骤S3中，结合子区顶点形变向量差异评估函数

和兴趣场函数δ(x,y)以确定所述误差估计特征值P_k的大小。

进一步的，所述误差估计特征值P_k的计算公式为

其中：

为角度参数；η为向量大小参数；

d_k为第k个子区对的大小度量参数

ρ₁,ρ₂为平衡子区形变误差和兴趣场的权值。

优选的，所述兴趣场函数δ(x,y)为二维高斯分布函数：

或，所述兴趣场函数δ(x,y)为或二维分段函数；

或，当无法确定变形分布情况或无需设置特别关注区域时，所述兴趣场函数δ(x,y)的取值为恒定值。

其中：

κ为均一化参数；

子区顶点A发生形变后的对应顶点A′形成的向量

子区顶点B形变后的对应顶点B′形成的向量

子区顶点C形变后的对应顶点C’形成的向量

子区顶点D形变后的对应顶点D’形成的向量

—表示向量i和向量j之间所夹的角

向量i与向量j分别代指向量

向量

向量

向量

且向量i与向量j相异设置；

η_ij—表示向量i和向量j的模大小之和；

M—为η_ij中的最大值。

优选的，所述步骤S2中，基于网格细分将所述当前状态图像、下一状态图像初步划分为若干子区。

优选的，所述步骤S4中，基于插值算法将相应的子区对分别再次细化。

本发明公开了一种初始子区细分装置，包括标记点识别模块，用于采集被测件的当前状态图像与下一状态图像后，分别识别所述当前状态图像与所述下一状态图像上的标记点，所述被测件表面设有带标记点的散斑图案；

初步子区划分模块，用于依据所述标记点将所述当前状态图像、下一状态图像初步划分为若干子区，其中，所述当前状态图像的子区为第一类子区，所述下一状态图像的子区为第二类子区；

误差估计特征值确定模块，用于匹配所述第一类子区与所述第二类子区以形成多个子区对，并确定每个所述子区对的误差估计特征值P_k；

子区细分模块，用于比较所述误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；若所述子区对的误差估计特征值P_k小于所述特征值阈值，则判断相应的子区对的划分已达标；若所述子区对的误差估计特征值P_k大于所述特征值阈值，则将相应的子区对分别再次细化后转入误差估计特征值确定模块,直至判断全部的所述子区对划分达标。

作为一优选方案，所述误差估计特征值确定模块中结合子区顶点形变向量差异评估函数和兴趣场函数以确认所述误差估计特征值P_k的大小。

作为一优选方案，所述误差估计特征值确定模块中，所述兴趣场函数δ(x,y)对被测件可能发生较大变形的区域或工件比较关注部分的变形设置较大的值。

作为一优选方案，所述误差估计特征值确定模块中，所述兴趣场函数δ(x,y)为二维高斯分布函数：

或，所述兴趣场函数为或二维分段函数；

或，在无法确定变形分布情况/无需设置特别关注区域时，所述兴趣场函数的取值为恒定值。

作为一优选方案，所述误差估计特征值确定模块中，所述误差估计特征值Pk的计算公式为

其中：

为角度参数，η为向量大小参数；

d_k为第k个子区对的大小度量参数；

ρ₁,ρ₂—平衡子区形变误差和兴趣场的权值。

进一步的，

计算方法如下所示：

其中：

κ为均一化参数；

子区顶点A与形变后的对应顶点A′形成的向量

子区顶点B与形变后的对应顶点B′形成的向量

子区顶点C与形变后的对应顶点C′形成的向量

子区顶点D与形变后的对应顶点D′形成的向量

—表示向量i和向量j之间所夹的角

向量i与向量j分别代指向量

向量

向量

向量

且向量i与向量j相异设置；

η_ij—表示向量i和向量j的模大小之和；

M—为η_ij中的最大值。

本发明还提供了一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任一项所述方法的步骤。

本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述方法的步骤。

本发明公开的初始子区细分方法、装置、计算机设备及存储介质，结合连续变形理论，充分考虑了子区各位置差异性可能造成的匹配误差，并提出采用子区顶点形变向量差异来衡量这种误差，指导子区细分过程。同时，兼顾考虑了实际工程测量过程中对于零件各部分关注度不同的状况，并结合模拟仿真等现有预测方法，提出了兴趣场函数对子区细分过程产生影响。此外，还添加了权衡因子，平衡子区顶点形变向量差异和兴趣场函数对子区细分的影响权重。

附图说明

图1为本发明所示的初始子区细分方法一实施例的工作流程示意图；

图2为图1所示实施例中子区细化图；

图3为本发明所示的初始子区细分方法另外一实施例的工作流程示意图；

图4为制斑流程示意图；

图5为图4所示实施例中标记点的位置示意图；

图6为发明所示的初始子区细分方法的硬件结构示意图

图7(a)为第一类子区的结构示意图；图7(b)为第二类子区的结构示意图；图7(c)为形变向量的示意图；

图8为兴趣场函数的形成示意图；

图9为图1为本发明所示的初始子区细分装置一实施例的结构示意图；

图10为是本发明实施例提供的一种计算机设备的结构示意图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例对本发明进行进一步的详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，而非以任何方式限制本发明的保护范围。

在说明书的全文中，相同的附图标号指代相同的元件。表述“和/或”包括相关联的所列相目中的一个或多个的任何和全部组合。在附图中，为了便于说明，已稍微夸大了物体的厚度、尺寸和形状。附图仅为示例而非严格按比例绘制。

还应理解的是，用语“包括”、“包括有”、“具有”、“包含”和/或“包含有”，当在本说明书中使用时表示存在所陈述的特征、步骤、整体、操作、元件和/或部件，但不排除存在或附加有一个或多个其它特征、步骤、整体、操作、元件、部件和/或它们的组合。

如在说明书中使用的用语“基本上”、“大约”以及类似的用于用作表示近似的用语，而不用作表示程度的用语，并且旨在说明将由本领域普通技术人员认识到的、测量值或计算值中的固有偏差。

除非另有限定，否则本文中使用的所有用语(包括技术用语和科学用语)均具有与本申请所属领域普通技术人员的通常理解相同的含义。还应理解的是，用语(例如在常用词典中定义的用语)应被解释为具有与它们在相关技术的上下文中的含义一致的含义，并且将不被以理想化或过度正式意义解释，除非本文中明确如此限定。

实施例一

如图1所示，本发明公开了一种初始子区细分方法，包括如下步骤：

S1采集被测件的当前状态图像与下一状态图像，所述被测件表面设有带标记点的散斑图案，分别识别所述当前状态图像与所述下一状态图像上的标记点；

S2依据所述标记点将所述当前状态图像、下一状态图像初步划分为若干子区，其中，所述当前状态图像的子区为第一类子区，所述下一状态图像的子区为第二类子区；

S3匹配所述第一类子区与所述第二类子区以形成多个子区对，并确定每个所述子区对的误差估计特征值P_k；

S4比较所述误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；若所述子区对的误差估计特征值P_k小于所述特征值阈值，则判断相应的子区对的划分已达标；如图2所示，若所述子区对的误差估计特征值P_k大于所述特征值阈值，则将相应的子区对分别再次细化后转入步骤S3,直至判断全部的所述子区对划分达标。

本发明所示的一种初始子区细分方法，以标记点构成的子区为基础，利用形变误差估计特征值控制子区细分过程，兼顾了精度需求和处理速度，为后续各像点的快速、准确初值估计提供了基础。

实施例二

如图3所示，本发明公开了一种初始子区细分方法，包括如下步骤：

S1采集被测件的当前状态图像与下一状态图像，所述被测件表面设有带标记点的散斑图案，分别识别所述当前状态图像与所述下一状态图像上的标记点；

本发明创新性的将标记点引入到了散斑图像应变测量方法中，由于标记点可以被系统轻易识别，因此相比现有的种子点扩散算法等技术中的初始子区划分，不需经过复杂的计算过程，且由于标记点位置的确定性，这些点具有绝对的匹配正确率，精度较其他方法高。

作为一优选方案，如图4所示，为保证标记点不被遮盖或破坏，可先对被测件表面制斑，然后进行标记点喷涂，喷涂方法可选用水转印、喷墨打印等。

作为一优选方案，标记点的样式和密度、排布方式等可根据实际测量需求决定。标记点不宜过密，否则将给标记点在被测件的喷涂过程造成一定的困扰。如图5所示，标记点可由任意具备尺度不变性、旋转不变性、受光照强度影响较小中一个或多个性质的图案组成，也可通过颜色区别等手段达到识别的目的。选定所需的标记点图案后，通过计算机控制带标记点的散斑图案生成，并将其喷涂在已经制斑的被测件表面即可。

完成被测件的制斑流程后，采用图6所示相机系统对变形前后的被测件进行图像采集。相机系统的硬件组成如图6所示，其包括两个或多个相机[1-1、1-2]、光照灯[2]、计算机及控制逻辑单元[3]和带标记点散斑图案的被测件[4]。相机可在被测件发生形变前后动态采集被测件表面的图像信息，光照灯用以提高光照质量。计算机及控制逻辑通过控制光照灯、相机来完成整个图像采集过程，以及对采集的散斑图像进行处理。

被测件通过表面制斑，从而被测件的不同区域位置可以被轻易识别。所不同的是，涂敷在被测件表面的散斑被人为的添加了可轻易被系统识别的标记图案或标记点(以下将标记图案或标记点统称为标记点)，从而通过简单的图像处理就可定位识别这些标记位置。根据本步骤中所添加的标记点类型的不同，可以依据颜色识别、模式识别等方法进行标记点的快速识别变形前后图像中的各个标记点。

S2依据所述标记点将所述当前状态图像、下一状态图像初步划分为若干子区，其中，所述当前状态图像的子区为第一类子区，所述下一状态图像的子区为第二类子区；

要完成快速初值估计，就需要将相机拍摄图像划分为若干个子区，这些子区内的像点具有相似的形变趋势，从而利用连续变形理论，通过像点的对应关系完成子区内所有像点的快速匹配。标记点作为可快速识别的单元，在形变前后的初值估计及后续的立体匹配过程中可作为具有绝对匹配精度的位置点。

本实施例中，变形前后图像的标记点识别出来之后，可经由网格划分技术将识别到的标记点连接成特定形状的网络(即将对应的标记点连接)，如划分为多个三角形网格、四边形网格或六边形网格等，从而以将进行变形前后的图像初步划分为若干子区，其中当前状态图像的各个子区记为第一类子区，下一状态图像中的各个子区记为第二类子区；有多种方法可实现网格划分，如可采用有限元进行网格划分。

S3匹配所述第一类子区与所述第二类子区以形成多个子区对，并确定每个所述子区对的误差估计特征值P_k。其中，k为为各子区对的编号。

依据连续变形理论，在同一方向上，任何两个标记点在变形前后都不可能交叠错位，因此，虽然变形前后标记点位置发生了改变，但可以根据标记点的相邻分布关系，确定变形后的标记点和变形前的哪一个标记点对应。此外，也可以通过深度学习方法将变形前后图像中的子区一一匹配对应。

子区匹配完毕之后，继续确定每个子区对的误差估计特征值P_k，以度量子区是否需要再细分。确定每个子区形变误差估计特征值P_k时，本实施例结合实际工程测量经验方法，引入兴趣场函数，将子区顶点形变向量差异评估函数

和兴趣场函数δ(x,y)结合来控制子区细分过程。

作为一优选方案，步骤S3中，误差估计特征值P_k的计算公式为

其中：

为角度参数，η为向量大小参数；

d_k为第k个子区对的大小度量参数；

ρ₁,ρ₂为平衡子区形变误差和兴趣场的权值。

本实施例中，d_k为第k个子区对中对应的第一类子区或第二类子区的边长、半径或者面积，其中，若选择第一类子区的边长、半径或者面积作为d_k的取值，后续计算时，则一直沿用第一类子区的参数；同理，若选择第二类子区的边长、半径或者面积作为d_k的取值，后续计算时，则一直沿用第二类子区中的参数。

ρ₁,ρ₂权值大小取决于对子区形变误差和对于某一区域的感兴趣程度。更进一步的，可使两者之和为1，且两者均不小于0。

函数

对子区在变形前后的各个顶点形变向量的差异性进行计算，当各个向量之间的角度差异和向量模大小差异越大时，函数

得到的值越大，被细分的可能性也就越大。

如图7所示，本实施例中以步骤S2中采用四边形网格进行子区划分、当前图像为未变形的原始图像，下一状态图像为变形后图像为例对

的大小确定进行说明。

对应的，依据如下公式确定

的大小：

其中：

κ为均一化参数；

子区顶点A与形变后的对应顶点A′形成的向量

子区顶点B与形变后的对应顶点B′形成的向量

子区顶点C与形变后的对应顶点C′形成的向量

子区顶点D与形变后的对应顶点D′形成的向量

—表示向量i和向量j之间所夹的角

向量i与向量j分别代指向量

向量

向量

向量

且向量i与向量j相异设置；

η_ij—表示向量i和向量j的模大小之和；

M—为η_ij中的最大值。

如图7所示，顶点ABCD为第一类子区的顶点，顶点A’B’C’D’为第二类子区对应的顶点，上述公式中

表示向量

与向量

之间所夹的角；

表示向量

与向量

之间所夹的角；

表示向量

与向量

之间所夹的角；

表示向量

与向量

之间所夹的角；

表示向量

与向量

之间所夹的角；

表示向量

与向量

之间所夹的角；

表示向量

和向量

的模大小之和；

表示向量

和向量

的模大小之和；

表示向量

和向量

的模大小之和；

表示向量

和向量

的模大小之和；

表示向量

和向量d的模大小之和；

表示向量

和向量

的模大小之和。

M为每两个向量模大小之和的最大值，即M＝Max(η_ab,η_ac,η_bd,η_cd,η_ad,η_bc)

当初步划分的子区判断不达标被再次细化后，根据选择的子区划分方法不同，子区的顶点可以是标记点，也可以是细分位置点。

此外，若当步骤S2中采用三角形、六边形网格进行子区划分时，

计算公式对应变化，为等同的方案，此处不在赘述。

如图8所示，兴趣场函数标志了对测量件不同区域的关注程度，可由实际测量环境和先验知识决定。比如，在一些场景下，测量件的大致变形趋势可由以往经验或软件模拟获得。此时，可通过兴趣场函数对测量件可能发生较大变形的区域设置较大的值以提高关注度。再如，有时对工件测量过程中会比较关注工件的某一部分的变形情况。此时同样可以通过兴趣场函数对这些部分设置较大的值。兴趣场函数的值越大，窗口细化等级也就越高，初值估计的结果就越精确。

其中，兴趣场函数可以是任意合理的函数，如二维高斯分布函数或二维分段函数，此外，在不清楚变形分布情况或无需设置特别关注区域时，可将兴趣函数规定为恒定值。

作为一优选方案，当兴趣场函数为二维高斯分布函数时，计算公式如下所示：

x，y是指图像卡迪尔坐标的两个轴向。

作为一优选方案，当兴趣场函数为二维分段函数时，计算公式如下所示：

S4比较所述误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；

S41若所述子区对的误差估计特征值P_k小于所述特征值阈值，则判断相应的子区对的划分已达标；

S42若所述子区对的误差估计特征值P_k大于所述特征值阈值，则将相应的子区对分别再次细化后转入步骤S3,直至判断全部的所述子区对划分达标。

子区窗口划分太过稀疏，则子区内初值估计难度越大，初值估计误差越大；子区窗口划分太过密集，则需要精确匹配的初始点太多，计算速度降低。本步骤中，对子区初步划分结果进行评估，比较每个所述子区对的误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对是否需要细化，其中特征值阈值为一个预设的数值，对不同的待测件，根据实际经验设置不同的数值。若所述误差估计特征值P_k大于所述特征值阈值，则依据所述标记点将对应的子区对别划分为若干子区后转入步骤S3,直至判断全部的所述子区对误差估计特征值P_k小于所述特征值阈值。

本步骤中，作为一优选方案，经由精确匹配算法(插值算法)对误差估计特征值Pk大于特征值阈值的子区进行划分(包括对变形前后的子区均进行划分)，即利用插值算法计算子区内某点变形前后位置的方法，将相应的子区对分别细化。

本发明所公开的初始子区细分方法，结合连续变形理论，充分考虑了子区各位置差异性可能造成的匹配误差，并提出采用子区顶点形变向量差异来衡量这种误差，指导子区细分过程。同时，兼顾考虑了实际工程测量过程中对于零件各部分关注度不同的状况，并结合模拟仿真等现有预测方法，提出了兴趣场函数对子区细分过程产生影响。此外还添加了权衡因子，平衡子区顶点形变向量差异和兴趣场函数对子区细分的影响权重。上述子区细化完成后，对任意的子区可利用各种快速初值匹配方法进行初值估计，如种子点扩散算法。

实施例三

如图9所示，本发明公开了一种初始子区细分装置10，包括：

标记点识别模块11，用于采集被测件的当前状态图像与下一状态图像后，分别识别所述当前状态图像与所述下一状态图像上的标记点，所述被测件表面设有带标记点的散斑图案；

初步子区划分模块12，用于依据所述标记点将所述当前状态图像、下一状态图像初步划分为若干子区，其中，所述当前状态图像的子区为第一类子区，所述下一状态图像的子区为第二类子区；

误差估计特征值确定模块13，用于匹配所述第一类子区与所述第二类子区以形成多个子区对，并确定每个所述子区对的误差估计特征值P_k；

子区细分模块14，用于比较所述误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；若所述子区对的误差估计特征值P_k小于所述特征值阈值，则判断相应的子区对的划分已达标；若所述子区对的误差估计特征值P_k大于所述特征值阈值，则将相应的子区对分别再次细化后转入误差估计特征值P_k确定模块,直至判断全部的所述子区对划分达标。

作为一优选方案，所述误差估计特征值确定模块13中，结合子区顶点形变向量差异评估函数和兴趣场函数以确认所述误差估计特征值P_k的大小。

作为一优选方案，所述误差估计特征值确定模块13中，所述兴趣场函数δ(x,y)对被测件可能发生较大变形的区域或工件比较关注部分的变形设置较大的值。

作为一优选方案，所述误差估计特征值确定模块13中，所述兴趣场函数δ(x,y)为二维高斯分布函数：

或，所述兴趣场函数为或二维分段函数；

或，在无法确定变形分布情况/无需设置特别关注区域时，所述兴趣场函数的取值为恒定值。

作为一优选方案，所述误差估计特征值确定模块13，所述误差估计特征值P_k的计算公式为

其中：

为角度参数，η为向量大小参数；

d_k为第k个子区对的大小度量参数；

ρ₁,ρ₂为平衡子区形变误差和兴趣场的权值。

进一步的，

计算方法如下所示：

其中：

κ为均一化参数；

子区顶点A与形变后的对应顶点A′形成的向量

子区顶点B与形变后的对应顶点B′形成的向量

子区顶点C与形变后的对应顶点C′形成的向量

子区顶点D与形变后的对应顶点D′形成的向量

—表示向量i和向量j之间所夹的角

向量i与向量j分别代指向量

向量

向量

向量

且向量i与向量j相异设置；

η_ij—表示向量i和向量j的模大小之和；

M—为η_ij中的最大值。

作为一优选方案，所述标记点识别模块11中，先对被测件表面制斑，然后进行标记点喷涂，喷涂方法可选用水转印、喷墨打印等。

作为一优选方案，标记点识别模块11包括两个或多个相机[1-1、1-2]、光照灯[2]、计算机及控制逻辑单元[3]和带标记点散斑图案的被测件[4]。相机可在被测件发生形变前后动态采集被测件表面的图像信息，光照灯用以提高光照质量。计算机及控制逻辑通过控制光照灯、相机来完成整个图像采集过程，以及对采集的散斑图像进行处理。

实施例四

图10所示是本发明实施例提供的一种计算机设备的结构示意图，如可以执行程序的智能手机、平板电脑、笔记本电脑、台式计算机、机架式服务器、刀片式服务器、塔式服务器或机柜式服务器(包括独立的服务器，或者多个服务器所组成的服务器集群)等。本实施例的计算机设备20至少包括但不限于：可通过系统总线相互通信连接的存储器21、处理器22，如图10所示。需要指出的是，图10仅示出了具有组件21-22的计算机设备20，但是应理解的是，并不要求实施所有示出的组件，可以替代的实施更多或者更少的组件。

本实施例中，存储器21(即可读存储介质)包括闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)，存储器21也可以是计算机设备20的外部存储设备，例如该计算机设备20上配备的插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)等。当然，存储器21还可以既包括计算机设备20的内部存储单元也包括其外部存储设备。本实施例中，存储器21通常用于存储安装于计算机设备20的操作系统和各类应用软件，例如方法实施例中的初始子区细分装置的程序代码等。此外，存储器21还可以用于暂时地存储已经输出或者将要输出的各类数据。

处理器22在一些实施例中可以是中央处理器(Central Processing Unit，CPU)、控制器、微控制器、微处理器、或其他数据处理芯片。该处理器22通常用于控制计算机设备20的总体操作。本实施例中，处理器22用于运行存储器21中存储的程序代码或者处理数据，例如运行初始子区细分装置，以实现方法实施例中的初始子区细分方法。

实施例五

本申请还提供一种计算机可读存储介质，如闪存、硬盘、多媒体卡、卡型存储器(例如，SD或DX存储器等)、随机访问存储器(RAM)、静态随机访问存储器(SRAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、可编程只读存储器(PROM)、磁性存储器、磁盘、光盘、服务器、App应用商城等等，其上存储有计算机程序，程序被处理器执行时实现相应功能。本实施例的计算机可读存储介质用于存储初始子区细分装置的程序代码，被处理器执行时实现方法实施例中的初始子区细分方法。

需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。

上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。

以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (9)

1.一种初始子区细分方法，其特征在于：包括如下步骤：

S1采集被测件的当前状态图像与下一状态图像，所述被测件表面设有带标记点的散斑图案，分别识别所述当前状态图像与所述下一状态图像上的标记点；

S2依据所述标记点将所述当前状态图像、下一状态图像初步划分为若干子区，其中，所述当前状态图像的子区为第一类子区，所述下一状态图像的子区为第二类子区；

S3匹配所述第一类子区与所述第二类子区以形成多个子区对，并确定每个所述子区对的误差估计特征值P_k；

S4比较所述误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；若所述子区对的误差估计特征值P_k小于所述特征值阈值，则判断相应的子区对的划分已达标；若所述子区对的误差估计特征值P_k大于所述特征值阈值，则将相应的子区对分别再次细化后转入步骤S3,直至判断全部的所述子区对划分达标。

2.根据权利要求1所述的一种初始子区细分方法，其特征在于，所述步骤S1中，对被测工件表面制斑后，进行标记点喷涂以获取所述被测件；

和/或，所述当前状态图像为变形前图像，所述下一状态图像为终态变形图像。

3.根据权利要求1所述的一种初始子区细分方法，其特征在于：所述步骤S3中，结合相似度误差估算函数

和兴趣场函数δ(x,y)以确定所述误差估计特征值P_k的大小 ；所述误差估计特征值P_k的计算公式为：

其中：

为角度参数；η为向量大小参数；

d_k为第k个子区对的大小度量参数；

ρ₁,ρ₂为平衡子区形变误差和兴趣场的权值。

4.根据权利要求3所述的一种初始子区细分方法，其特征在于：所述兴趣场函数δ(x,y)为二维高斯分布函数：

或，所述兴趣场函数δ(x,y)为二维分段函数；

或，当无法确定变形分布情况或无需设置特别关注区域时，所述兴趣场函数δ(x,y)的取值为恒定值。

5.根据权利要求3所述的一种初始子区细分方法，其特征在于，所述相似度误差估算函数

计算公式为：

其中：

κ为均一化参数；

子区顶点A与形变后的对应顶点A′形成的向量

子区顶点B与形变后的对应顶点B′形成的向量

子区顶点C与形变后的对应顶点C′形成的向量

子区顶点D与形变后的对应顶点D′形成的向量

—表示向量i和向量j之间所夹的角

向量i与向量j分别代指向量

向量

向量

向量

且向量i与向量j相异设置；

η_ij—表示向量i和向量j的模大小之和；

M—为η_ij中的最大值。

6.根据权利要求1所述的一种初始子区细分方法，其特征在于：

所述步骤S2中，基于网格细分将所述当前状态图像、下一状态图像初步划分为若干子区；

和/或，所述步骤S4中，基于插值算法将相应的子区对分别再次细化。

7.一种初始子区细分装置，其特征在于：包括：

标记点识别模块，用于采集被测件的当前状态图像与下一状态图像后，分别识别所述当前状态图像与所述下一状态图像上的标记点，所述被测件表面设有带标记点的散斑图案；

初步子区划分模块，用于依据所述标记点将所述当前状态图像、下一状态图像初步划分为若干子区，其中，所述当前状态图像的子区为第一类子区，所述下一状态图像的子区为第二类子区；

误差估计特征值确定模块，用于匹配所述第一类子区与所述第二类子区以形成多个子区对，并确定每个所述子区对的误差估计特征值P_k；

子区细分模块，用于比较所述误差估计特征值P_k与预设的特征值阈值大小，判断对应的子区对的划分是否达标；若所述子区对的误差估计特征值P_k小于所述特征值阈值，则判断相应的子区对的划分已达标；若所述子区对的误差估计特征值P_k大于所述特征值阈值，则将相应的子区对分别再次细化后转入误差估计特征值确定模块,直至判断全部的所述子区对划分达标。

8.一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述处理器执行所述计算机程序时实现权利要求1至6任一项所述方法的步骤。

9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于：所述计算机程序被处理器执行时实现权利要求1至6任一项所述方法的步骤。

Images