基于时序变形预测的数字图像相关加速方法.pdfVIP

第 29卷 第 6期 实 验 力 学 V01．29 No．6 2014年 12月 joURNALOFEXPERIMENTAIMECHANICS Dec．2O14 文章编号 ：i001—4888(2014)06—711-08 基于时序变形预测的数字图像相关加速方法 席涛 ，熊宸 ，孔繁羽 ，许峰 ，缪泓h (L中国科学技术大学 中科院材料力学行为和设计重点实验室，安徽合肥 230027；2．四川农业大学 理学院，四川I雅安 625014) 摘要 ：利用数字图像相关求解连续变形物体的位移场和应变场时，会遇到处理速度非常慢的问 题 ，原 因是相关算法计算量大且忽略了物体在时间轴上的变形规律。本文提 出了一种基于变形 预测的数字图像相关方法，该方法利用物体在 时间轴上的变形规律，通过 已经得到 的变形值来 预测后面时刻的变形初值。首先分析 了物体在时间轴上的变形规律 ，然后结合 实际应用对其进 行修正，得到物体下一步变形的初值 ，最后通过 NR方法 (Newton-RaphsonMethod)得到物体的 位移场和应变场。通过计算机模拟和金属试件拉伸 实验并结合 GPU(GraphicProcessingUnit) 编程验证 了该方法的有效性 。计算机模拟和金属试件拉伸 实验 的计算结果表明，该方法能够在 计算精度保持不变的情况下，使计算速度提高4～7倍 。 关键词：数字图像相关；变形预测；初值搜 索 ‘ 中图分类号 ：O34 文献标识码 ：A DOI：10．7520／1001-4888—13—146 O 引言 数字图像相关方法 【2_(又称数字散斑相关方法)是在上世纪 80年代初提出的，利用被测试件变形 前后灰度的变化来测量其位移场和应变场 的光学测量方法。 在数字图像相关的发展过程 中，高精度和快速度一直是追求的目标 。现在 已经有许多亚像素算法 ， 比如 ：灰度差值法[3]、曲面拟合法 ]、梯度法[5]、梯度迭代法 ]、牛顿拉普逊法 (Newton—RaphsonMeth— od)L7]、拟牛顿法 (Quasi—NewtonMethod)_8]、列文伯格一马夸尔特法 (Lenvenberg—MarquardtMethod)[。] 等，这些算法的计算精度从 0．005 0．1像素不等口 。通常情况下，这 已经能够满足实际情况对精度的 要求 。但是由于这些算法都是局部优化算法 ，因此必须提供适 当的初值才能得到精确的解。而全局优 化算法 比如 ：遗传算法n 、粒子群算法口胡等都需要巨大的计算量 ，使其很难达到快速度的要求。 在局部优化算法 中，可以通过整像素搜索法口引、光流法f“等提供初值 ，也可以通过现在常用的种子 点法n阳来提供初值 。但是这些方法要么计算量太大，要么由于是串行计算，很难达到快速的要求 。特 别是 当要快速得到被测试件在一个时间段 内的连续位移场和变形场时 ，对计算速度的要求则更高 。 针对上述情况 ，本文提出了一种基于变形预测的方法——变形预测 图像相关方法。该方法利用物 体在连续变形过程中，相邻时刻的位移场和变形场的改变量相近的特征 ，通过上个时刻 已经得到的位移 场来预测试件当前位移场的初值 ，并通过测试点进行修正 ，从而得到可靠的变形初值，最后利用 NR方 法计算变形物体的位移场。该方法不仅能简单快速地找到变形初值而且适合并行运算，很大程度上提 高了计算速度 ，为数字图像相关方法的实时性处理打下了基础。 本文首先介绍了数字图像相关方法的基本原理 ，接着介绍了基于变形预测的图像相关方法 ，并通过 * 收稿 日期：2013—09—28；修订 日期：2014—0124 基金项 目：国家重点基础研究发展计划资助项 目(2012CB937500)，国家 自然科学基金项 目11472267) 通讯作者 ：缪泓(1968一)，男 ，教授，主要从事实验力学测试技术、数字图像分析等领域 的研究 E—mail：miaohong@ustc．edu．Crl 712 实 验 力 学 (2014年)第 29卷 模拟散斑 图旋转实验和