用于数字图像相关位移测量的一种改进爬山算法 下载: 944次
1 引言
数字图像相关(DIC)位移测量方法是20世纪80年代由日本Yamaguchi、美国的Peters和Ranson等共同提出的[1]。DIC方法因具有光路简单、非接触、全场测量、对环境要求低等优点,在光测力学相关研究中得到越来越广泛的应用。提高相关搜索速度是DIC方法研究的重要方向。在早期的研究中,大都采用逐点搜索法,搜索时间长速度慢。后来经过不断改进,提出了许多基于经典数学理论的相关搜索方法,如粗-细搜索法、Newton-Raphson迭代算法、双参数法等,以及后来基于智能的搜索算法,如爬山法、人工鱼群算、粒子群算法、遗传算法等[2-5]。
1983年,Peters等[6]提出粗细搜索法,将整个搜索过程分两部分进行:首先用粗搜索定位搜索范围,然后利用细搜索进行精确搜索,计算物体实际位移,提高算法的精度。1988年,Sutton等[7]提出Newton-Raphson算法,将数字散斑相关方法的测量精度提高了一个数量级。1989年,Sutton等[8]将Newton-Raphson迭代算法加以改进,对散斑图进行二元三次样条插值,提高了DIC法的测量速度和精度。1991年,陈绍春等[9]提出双参数法,通过参数拟协调元的直接分析来计算位移,后来其被用于DIC分析中。1996年,Zhao等[10]提出了爬山搜索算法,提高了DIC的测量精度。2012年,张小莲等[11]提出了一种引入停止机制的爬山法,并应用在了风力发电领域。2018年,姚万业等[12]提出了基于风速预测的改进爬山法。
粗-细搜索法是对位移及位移导数同时搜索,计算量大,搜索时间长; Newton-Raphson迭代算法相对于粗-细搜索法比较省时,但当计算目标函数的Hessian矩阵以及其逆矩阵的维数较高时计算量也很大[13];双参数法因公式推导复杂,计算量大,已经很少用于DIC测量中;爬山法采用局部寻优方法,只需少量的搜索步骤,计算速度相对较快,但易陷入局部极值。智能算法在亚像素分辨率位移测量时优势明显,2018年,葛朋祥等[14]将遗传算法、粒子群算法、人工鱼群算法等智能算法在匹配精度和搜索速度等方面进行了对比分析,认为微小位移测量时,遗传算法的综合性能较优。但并没有解决测量大位移时时间较长且误差较大的问题。
为了改善传统DIC方法易陷入局部最优及大位移计算速度过慢等缺点,本文对爬山法进行了改进,通过使用自适应的改变步长方法,避免陷入局部寻优,并通过模拟散斑图对基于改进爬山法的DIC方法在大位移测量中的适用性进行验证。
2 DIC方法介绍
2.1 DIC基本原理
选择变形前后的两幅图像作为研究对象,变形前图像为参考图像,变形后图像为目标图像。在参考图像中以待测点P(x,y)为中心选择一个大小为(2N+1)×(2N+1)的矩形参考子区,在目标图像中选择相同尺寸的目标子区。对参考子区与目标子区进行相关匹配,用相关系数的大小来衡量匹配程度,相关系数最大值所对应的点P'(x',y')即为目标子区的中心点,根据点P(x,y)和P'(x',y')的坐标差值可求出两幅图像的相对位移量d=
目前,衡量匹配程度的相关系数有数十种,根据Tong等[18-19]对相关系数的研究,零均值归一化互相关函数对图像灰度值的线性变化不敏感,相关计算的精度和稳健性均较高,因此本文以其作为评价参考子区和变形子区匹配程度的函数,表达式为
式中:M为子区尺寸;f(x,y)为参考图中点(x,y)处的灰度值;g(x',y')为目标图中点(x',y')处的灰度值;
2.2 爬山法原理
爬山算法是一种局部择优的方法,采用启发式方法对深度优先搜索进行改进,它利用反馈信息帮助生成解的决策, 属于人工智能算法的一种。但也因为爬山算法局部寻优以及不断缩小寻优区域的特性,当测量范围较大及爬山距离较远时,容易陷入局部极值。
3 改进爬山法
3.1 内部自适应循环
1)在参考图像中选择以待测点P(x,y)为中心的参考子区,并在目标图像中以相同坐标选择起算点Q1(x,y)为中心的目标子区。应用(1)式计算两个子区的相关系数,以作为目前最大相关系数,用C'max(k)表示,其中k为内循环计数,其对应的Qk(x,y)作为第k次内循环的起算点。
2)每次内循环以内步长Lin为间隔,以起算点Qk(x,y)为中心选择八邻点,分别计算以每个点为中心的子区与目标子区的相关系数Ci(i=1,…,8)。八邻点的像素坐标位置如
3)比较C'max(k)与 Ci(i=1,…,8)的大小,并改变Lin(k)与C'max(k)。当前起算点相关系数小于八邻点相关系数中的最大值时,以这个最大值点作为新的起算点,以这个最大值作为目前最大相关系数。否则,起算点与目前最大相关系数不变,改变内步长找到新的八邻点,新的最大相关系数为
当前起算点相关系数大于八邻点相关系数中的最大值时,内步长改变过程可表示为
为了减少多余计算,内步长根据理想相关系数的最大值1与当前起算点相关系数C'max(k)的比值自适应改变,重复循环直到Lin(k)=0,得到局部极值Cfmax=C'max(k)。
将内部循环的最后结果Cfmax作为重要依据输入外循环中。
3.2 基于阈值的自适应外循环
改进爬山法使用3.1节中得到的Cfmax,以设定好的阈值Cth为判断依据,这个阈值在使用模拟散斑图时一般设为0.85,即当相关系数大于0.85时认为搜索成功。并提出了一个外循环步长Lout为内循环步长Lin赋值。Lout以整像素最小单位1 pixel为初始值,根据理想相关系数最大值1与当前局部极值Cfmax的差值进行自适应改变,并且为了加快搜索速度设置了一个步长改变系数S。
当Cfmax<Cth,改变Lout的值:
为Lin重新赋值,进入内循环:
式中:K为外循环计数;S为外步长改变系数。
当Cfmax>Cth时,将Cfmax作为正确结果输出。
改进爬山法流程图如
4 实验验证
4.1 传统爬山法与改进爬山法结果对比
将一幅模拟散斑图向固定方向平移1~50 pixel位移,得到50幅位移图像,如
对这一系列散斑图分别用传统爬山法与改进爬山法进行对比测试,测试时使用的子区尺寸为41 pixel×41 pixel,部分测试结果如
可以看出,传统爬山法搜索时间波动较大,并且不能保证搜索精度,而改进爬山法的搜索速度相比传统爬山法有明显提升,并且能够保持稳定的搜索精度。从中提取位移为50 pixel的一组结果进行对比,在参考图中随机选择坐标为(94,94)的感兴趣点,目标图像中的正确搜索结果坐标应为(144,144)。
图 5. 系列位移图像。(a)参考图;(b) 1~50 pixel位移
Fig. 5. Serial displacement images. (a) Reference map; (b) 1--50 pixel displacement
传统爬山法相关系数搜索分布图如
改进爬山法搜索结果如
表 1. 改进爬山法与传统爬山法对比实验结果
Table 1. Comparison of experimental results between improved climbing method and traditional climbing method
|
4.2 步长改变系数S对计算时间的影响
利用上面已有的一系列位移图,在目标图像中以坐标点(103,103)确定目标子区,在所有位移图像中对这一目标子区进行搜索。同一组参数进行多次计算,对30次正确结果所用时间取均值。所用电脑配置如下:处理器Intel(R) Core(TM) i5_7200U,2.5 GHz,运行软件环境为MATLAB r2014a。
图 8. 步长改变系数S对搜索时间的影响。(a) 位移为25 pixel;(b) 子区尺寸51 pixel× 51 pixel
Fig. 8. Influence of step-size change coefficient S on search time. (a) Displacement is 25 pixel; (b) subarea size is 51 pixel×51 pixel
由
4.3 测量位移大小对计算时间的影响
图 9. 测量位移对搜索时间的影响。(a) 子区尺寸51 pixel×51 pixel;(b) S=9
Fig. 9. Influence of measured displacement on search time. (a) Sub area size is 51 pixel×51 pixel; (b) S=9
可以看出搜索时间随搜索位移的增大呈增大趋势,在保证搜索结果正确的基础上,当搜索位移在30 pixel以下时,能够保证搜索时间在0.1 s以内,当搜索位移大于30 pixel时,搜索时间偶尔会大于0.2 s。以上实验体现了改进爬山法能够满足1~50 pixel范围内的位移测量,并且能够保证较快的搜索速度。
4.4 更大位移的测量
为证明改进爬山法对较大位移测量的适用性,使用更大位移模拟散斑图对改进爬山法进行测试。对于尺寸为256 pixel×256 pixel的图片来说,100 pixel已经是相对较大的位移。如果位移大于100 pixel,就会出现搜索子区移动到图片范围之外的问题,导致无法搜索到正确结果。
表 2. 位移为90~100 pixel时改进爬山法的测量结果
Table 2. Measurement results of improved climbing method for displacement of 90--100 pixel
|
5 结论
基于位移量为100 pixel以内的测量实验,对改进爬山法的性能进行了实验测试。实验结果表明,改进爬山法解决了传统爬山法易陷入局部极值的缺点。在测量50 pixel以内位移时,当步长改变系数S小于12时,改进爬山法搜索时间随着S的增大而减小。改进爬山法能够满足较大位移的测量(小于100 pixel),并能够保证一定的搜索速度。
[1] 芮嘉白, 金观昌, 徐秉业. 一种新的数字散斑相关方法及其应用[J]. 力学学报, 1994, 26(5): 599-607.
Rui J B, Jin G C, Xu B Y. A new digital speckle correlation method and its application[J]. Acta Mechanica Sinica, 1994, 26(5): 599-607.
[2] Chu T C, Ranson W F, Sutton M A. Applications of digital-image-correlation techniques to experimental mechanics[J]. Experimental Mechanics, 1985, 25(3): 232-244.
[3] Bruck H A. McNeill S R, Sutton M A, et al. Digital image correlation using Newton-Raphson method of partial differential correction[J]. Experimental Mechanics, 1989, 29(3): 261-267.
[4] Yamaguchi I. A laser-speckle strain gauge[J]. Journal of Physics E: Scientific Instruments, 1981, 14(11): 1270-1273.
[5] Peters W H, Ranson W F. Digital imaging techniques in experimental stress analysis[J]. Optical Engineering, 1982, 21(3): 213427.
[6] Peters W H, Ranson W F, Sutton M A, et al. Application of digital correlation methods to rigid body mechanics[J]. Optical Engineering, 1983, 22(6): 226738.
[7] Sutton M A. McNeill S R, Jang J, et al. Effects of subpixel image restoration on digital correlation error estimates[J]. Optical Engineering, 1988, 27(10): 271070.
[8] Sutton M A, Turner J L, Bruck H A, et al. Full-field representation of discretely sampled surface deformation for displacement and strain analysis[J]. Experimental Mechanics, 1991, 31(2): 168-177.
[9] 陈绍春, 石钟慈. 构造单元刚度矩阵的双参数法[J]. 计算数学, 1991, 13(3): 286-296.
Chen S C, Shi Z C. Two-parameter method for constructing element stiffness matrix[J]. Computational mathematics, 1991, 13(3): 286-296.
[10] Zhao W Z, Jin G C. An experimental study on measurement of Poisson's ratio with digital correlation method[J]. Journal of Applied Polymer Science, 1996, 60(8): 1083-1088.
[11] 张小莲, 李群, 殷明慧, 等. 一种引入停止机制的改进爬山算法[J]. 中国电机工程学报, 2012, 32(14): 128-134.
Zhang X L, Li Q, Yin M H, et al. An improved hill-climbing searching method based on halt mechanism[J]. Proceedings of the CSEE, 2012, 32(14): 128-134.
[12] 姚万业, 贾昭鑫, 黄璞. 基于风速预测的改进爬山法最大功率追踪策略[J]. 电力科学与工程, 2018, 34(2): 44-49.
Yao W Y, Jia Z X, Huang P. Improved MPPT strategy for climbing method based on wind speed prediction[J]. Electric Power Science and Engineering, 2018, 34(2): 44-49.
[13] 王怀文, 亢一澜, 谢和平. 数字散斑相关方法与应用研究进展[J]. 力学进展, 2005, 35(2): 195-203.
Wang H W, Kang Y L, Xie H P. Advance in digital speckle correlation method and its application[J]. Advances in Mechanics, 2005, 35(2): 195-203.
[14] 葛朋祥, 叶沛, 李桂华. 基于遗传算法的数字图像相关法在微位移测量中的应用[J]. 光学学报, 2018, 38(6): 0612006.
[15] 张华俊, 李桂华, 刘程, 等. 基于SURF特征匹配的数字图像相关变形初值可靠估计[J]. 光学学报, 2013, 33(11): 1112005.
[16] 葛朋祥, 李桂华. 一种组合式DIC算法在亚像素级位移测量中的应用[J]. 激光与光电子学进展, 2018, 55(11): 111202.
[17] 葛继科, 邱玉辉, 吴春明, 等. 遗传算法研究综述[J]. 计算机应用研究, 2008, 25(10): 2911-2916.
Ge J K, Qiu Y H, Wu C M, et al. Summary of genetic algorithms research[J]. Application Research of Computers, 2008, 25(10): 2911-2916.
[18] Tong W. An evaluation of digital image correlation criteria for strain mapping applications[J]. Strain, 2005, 41(4): 167-175.
[19] Pan B, Xie H M, Wang Z Y. Equivalence of digital image correlation criteria for pattern matching[J]. Applied Optics, 2010, 49(28): 5501-5509.
Article Outline
叶沛, 张梅, 马万龙, 朱天天, 李桂华. 用于数字图像相关位移测量的一种改进爬山算法[J]. 激光与光电子学进展, 2020, 57(2): 021506. Ye Pei, Zhang Mei, Ma Wanlong, Zhu Tiantian, Li Guihua. Improved Climbing Algorithm for Digital Image Correlation Displacement Measurements[J]. Laser & Optoelectronics Progress, 2020, 57(2): 021506.