一种高动态范围激光散斑数字图像相关变形测量方法与流程

1.本发明涉及变形测量技术领域，特别是指一种高动态范围激光散斑数字图像相关变形测量方法。


背景技术：

2.近年来，基于数字图像相关方法的非接触式变形测量技术在光测力学、变形测量和各类工程测量领域得到了广泛的应用。数字图像相关法最早是在上世纪80年代初由日本的山口一郎
1.和美国南卡罗来纳大学的peter和ranson等人
2.相互独立提出，通过使用数字相机拍摄图像，然后使用算法比较变形前后的数字图像中的子区相关度，来确定变形前的子区在变形后的图像中的位置。
3.数字图像相关方法需要被测试件表面具有随机分布的特征散斑图样。除了常用的利用人工随机喷涂、机器制作或散斑模板等工艺手段来制作的“人工散斑”之外，也有学者研究使用了利用单波长高相关性的激光光源，照射在具有光学粗糙表面的被测试件上从而得到反射回来的特征图样，这种特征图样被称为“激光散斑
[3]”，同样具有位置随机分布、灰度随机分布的特点，适用于数字图像相关方法的相关函数运算。
[0004]
相比人工散斑，激光散斑
[4]
更适合于高温应用，因为它不需要寻找耐高温的特殊粘合剂或散斑涂料，并且能够应用与人工散斑材料无法达到的超高温度。但是由于激光光源的高相干特性，其在某些强反射表面，如低粗糙度表面、结构边缘凸起位置等，容易出现光的强反射现象，从而导致上述区域拍摄得到的激光散斑图像中出现了过曝光问题，具体为灰度值超出了工业相机的动态范围，从而采集到的激光散斑图像呈现“全白现象”，使得该区域不具有“位置随机分布、灰度随机分布的特征图样”。因此需要一种高动态范围激光散斑数字图像相关变形测量方法和系统，即hdr
‑
dic方法。其中，高动态范围(high dynamic range,hdr)是相对于低动态范围(low dynamic range,ldr)的一种概念，一幅hdr激光散斑图像需要使用多张ldr激光散斑图像经过算法合成得到。
[0005]
2020年，有学者提出使用dmd装置对强反射表面的“人工散斑”进行hdr合成，从而部分抑制人工散斑的强反射现象
[5]
。相比较于该方法，本发明与其的主要不同点在于：1、该方法需要借助dmd装置来实现不同照度，本发明不需要dmd装置，本发明使用的是单波长激光器，并且本发明使用了窄带宽滤波镜组。使用装置设备存在本质不同。2、该方法在图像融合算法使用了文献中提到的m1&m2两幅mask图像融合函数(公式1
‑
4)，本发明使用的图像融合函数是本发明提出的“自适应灰度级数映射hdr激光散斑图像融合算法”，并且在特征提取函数使用了本发明提出的“激光散斑hdr图像序列分割算法”。3、该方法使用的是人工散斑图样，人工散斑与激光散斑存在本质的不同，散斑的评价方法和影响散斑效果的因素也完全不同。4、该方法只使用了2幅ldr图像用于合成，本发明使用了4幅ldr图像合成hdr激光散斑图像。
[0006]
参考文献：
[0007]
[1]yamaguchi i.a laser
‑
speckle strain gauge[j].journal of physics e:
scientific instruments.1981,14:1270～1273.
[0008]
[2]w.h.peters,w.f.ranson.digital imaging techniques in experimental stress analysis[j].optical engineering.1981，21:427～431.
[0009]
[3]anwander m,zagar b g,weiss b,et al.noncontacting strain measurements at high temperatures by the digital laser speckle technique[j].experimental mechanics,2000,40(1):98
‑
105.
[0010]
[4]pradille c,bellet m,chastel y.a laser speckle method for measuring displacement field.application to resistance heating tensile test on steel[c]//applied mechanics and materials.trans tech publications,2010,24:135
‑
140.
[0011]
[5]suresh,vignesh,and beiwen li."digital image correlation for highly reflective objects using digital micro
‑
mirror device."emerging digital micromirror device based systems and applications xii.vol.11294.international society for optics and photonics,2020.


技术实现要素：

[0012]
本发明要解决的技术问题是提供一种高动态范围激光散斑数字图像相关变形测量方法，该方法使用大功率单波长激光器作为主动光源，使用空间滤波扩束镜组用来在一定距离下投射激光从而生成激光散斑，使用窄带宽滤波镜组用来过滤杂光等噪声干扰，使用灰度工业相机采集多幅ldr图像后续融合使用，使用“激光散斑hdr图像序列分割算法”和“自适应灰度级数映射hdr激光散斑图像融合算法”两种方法来获得高动态范围的激光散斑图像，即hdr激光散斑图像。该方法和系统经验证可以有效改善激光散斑在高反射表面的过曝光问题，实现全场dic变形计算。本发明经过验证可以应用与高温环境下的棒状高温合金试样单轴拉伸实验。
[0013]
该方法包括步骤如下：
[0014]
s1：搭建采集系统，调试激光投射面大小，调试镜头光圈到激光散斑清晰成像的位置；
[0015]
s2：选取4组不同的ldr图像曝光时间t1～t4，由低到高输入到采集系统中，等待后续自动采集图像；
[0016]
s3：根据设定好的拍摄参数，变形前拍摄4幅ldr激光散斑图像，保存为i1～4，变形后拍摄4幅ldr激光散斑图像，保存为j1～4；
[0017]
s4：使用激光散斑hdr图像序列分割算法，把图像每一行设置为一个组，从每行中间像素位置，向左右逐渐累加像素的灰度值，作为局部灰度总和值，基于局部灰度总和阈值不超过50(8bit图像，0
‑
255灰度值)控制分割像素区域，把上述累加灰度值总和不超过50的像素区域作为过曝光的无效像素区域，同理图像每一行进行上述分割，得到整张图像的过曝光无效像素区域。最终得到变形前图像i1～4和变形后图像j1～4两组激光散斑图中过曝光的无效像素区域；
[0018]
s5：使用自适应灰度级数映射hdr激光散斑图像融合算法，以ldr激光散斑图像的正常曝光的激光散斑区域为参考，统计该区域的灰度分布情况，绘制灰度分布直方图，获得该区域的灰度级数分布，将分布作为映射参量，对变形前图像i1～4和变形后图像j1～4进
行映射融合(即将ldr激光散斑的正常曝光区域作为输入图像，以前面获取的映射参量为灰度级数参考值，把正常曝光区域做直方图均衡化，然后i1～4和j1～4各组均衡化后的图像对应像素位置直接累加，可以得到2张i图和j图)得到2幅分别代表变形前和变形后图像组的hdr激光散斑图像，然后使用多因子融合散斑评价指标(multi
‑
factor fusion index,mffi)评价融合区域和正常曝光区域的激光散斑质量，如果融合区域的数值低于正常曝光区域数值的30％，则返回重新升高或降低映射参量再次融合hdr激光散斑图像；
[0019]
s6：使用融合得到的两幅变形前后hdr激光散斑图像进行数字图像相关运算，最终得到全场变形信息。
[0020]
做数字图像相关运算时，可使用包括互相关算法，如零均值归一化互相关算法(zncc)或距离和算法，如零均值最小平方距离和算法(znssd)。
[0021]
其中，s1中采集系统包括激光器、空间滤波扩束镜组、窄带宽滤波镜组和工业相机镜头，空间滤波扩束镜组位于激光器和待测试件之间，前端设置窄带宽滤波镜组的工业相机镜头置于空间滤波扩束镜组的上方，且倾斜放置，激光器光束直线和相机拍摄光轴之间的倾角不超过30度。
[0022]
空间滤波扩束镜组通过透镜的组合实现对激光器发射的单线激光的均匀扩束，射出的是一个圆形面光源，面光源内部各处的激光功率近似均匀。
[0023]
本发明的上述技术方案的有益效果如下：
[0024]
上述方案中，使用了大功率单波长蓝光激光器、空间滤波扩束镜组、窄带宽滤波镜组和高精度工业相机等激光散斑生成和图像采集设备，能够得到质量更好的激光散斑图案，并且更适应于高温环境下的激光散斑生成和图像采集任务。
[0025]
本发明设计了一种“激光散斑hdr图像序列分割算法”和“自适应灰度级数映射hdr激光散斑图像融合算法”，能够基于局部区域灰度值自动拾取过曝光区域，能够根据正常曝光区域自动获取灰度技术的映射参量，能够使用mffi评价因子动态评价融合的hdr激光散斑图像质量。
[0026]
综上，本发明为基于激光散斑的数字图像相关方法提供了有效的高动态范围成像改进方法和测量系统，弥补了传统方法在高反射表面存在的激光散斑过曝光从而导致测量失效的现象。并且能够满足数字图像相关测量所需要的测量精度，将会很好的推动基于数字图像相关方法的非接触式材料/构件变形检测技术的发展。
附图说明
[0027]
图1为本发明的高动态范围激光散斑数字图像相关变形测量方法流程图；
[0028]
图2为本发明所使用硬件系统布置简图；
[0029]
图3为本发明所使用的“激光散斑hdr图像序列分割算法”和“自适应灰度级数映射hdr激光散斑图像融合算法”示意图。
具体实施方式
[0030]
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
[0031]
本发明提供一种高动态范围激光散斑数字图像相关变形测量方法。
[0032]
如图1所示，该方法包括步骤如下：
[0033]
s1：搭建采集系统，调试激光投射面大小，调试镜头光圈到激光散斑合适位置，使能够清晰成像；
[0034]
s2：选取4组不同的ldr图像曝光时间t1～t4，由低到高输入到采集系统中，等待后续自动采集图像；
[0035]
s3：根据设定好的拍摄参数，变形前拍摄4幅ldr激光散斑图像，保存为i1～4，变形后拍摄4幅ldr激光散斑图像，保存为j1～4；
[0036]
s4：使用激光散斑hdr图像序列分割算法，基于局部灰度总和阈值控制对变形前图像i1～4和变形后图像j1～4两组激光散斑图中过曝光的无效区域进行分割；
[0037]
s5：使用自适应灰度级数映射hdr激光散斑图像融合算法，以ldr激光散斑图像的正常曝光的激光散斑区域为参考，求解其散斑灰度分布直方图，获得该区域的灰度级数分布，将分布作为映射参量，对变形前图像i1～4和变形后图像j1～4进行映射融合得到2幅分别代表变形前和变形后图像组的hdr激光散斑图像，然后使用多因子融合指标评价融合区域和正常曝光区域的激光散斑质量，如果融合区域的数值低于正常曝光区域数值的30％，则返回重新调整映射参量再次融合hdr激光散斑图像；
[0038]
s6：使用融合得到的两幅变形前后hdr激光散斑图像进行数字图像相关运算，最终得到全场变形信息。
[0039]
s1中采集系统包括激光器、空间滤波扩束镜组、窄带宽滤波镜组和工业相机镜头，如图2所示，空间滤波扩束镜组位于激光器和待测试件之间，前端设置窄带宽滤波镜组的工业相机镜头置于空间滤波扩束镜组的上方，且倾斜放置，保证激光器光束直线和相机拍摄光轴之间的倾角不超过30度。采集系统中，大功率单波长激光器输出的光线需要经过空间滤波扩束镜组进行扩束，然后照射到待测试件的表面，最后待测试件的激光散斑通过窄带宽滤波镜组，然后进入工业相机，采集成为数字图像。
[0040]
其中，空间滤波扩束镜组通过透镜的组合实现对激光器发射的单线激光的均匀扩束，射出的是一个圆形面光源，面光源内部各处的激光功率近似均匀。
[0041]
如图3所示，为本发明所使用从s4到s5的融合过程。
[0042]
如图3，采集到的ldr图像序列由于动态范围较低，会存在部分过曝光区域，由黑色矩形标识出。实际上每一幅ldr激光散斑图像的过曝光区域是不同的，需要使用自动提取算法得到图中黑色实线矩形框的过曝光区域。通过提取黑色虚线中的正常曝光区域图像，得到自动提取直方图从而获得灰度映射所需要的灰度级数，然后再把ldr序列的过曝光区域映射合成。最终合成的hdr激光散斑图还需要对于“合成区域”和“原始区域”使用mffi评估方法进行评估判断。
[0043]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。