一种基于单目视觉的低频振动校准方法与流程

本发明属于振动计量与测试领域，尤其适用于宽频率范围内高精度、稳定、可靠的低频振动校准。

背景技术：

低频振动传感器及低频振动测量仪正在越来越多地应用于地震预测、桥梁与建筑监测、矿质勘探等领域。低频振动传感器及测量仪的灵敏度对于其准确的振动测量极其重要。低频振动校准用于确定低频振动传感器及测量仪的灵敏度，其是确保振动测量数据准确、可靠、有效的前提。因此，对于低频振动校准的研究具有十分重要的意义。

常用的低频振动校准方法分为比较法与绝对法两大类，其中，最为常用的绝对法包括激光干涉法与地球重力法。比较法利用灵敏度已知的标准低频振动传感器与被校低频振动传感器及测量仪背靠背安装，使它们具有相同的输入激励。然而，其校准精度受限于标准传感器的校准精度，通常不能满足高精度校准的需求；地球重力法通过转台向被校低频振动传感器及测量仪提供精确已知的当地重力加速度作为激励，以实现高精度的校准。但受到转台离心加速度的影响，其最大校准频率通常较低；激光干涉法通过激光干涉仪测量长冲程振动台向被校低频振动传感器及测量仪提供的输入激励，当激励信号频率大于0.2hz时校准精度高，当小于0.2hz时由于存在由导轨弯曲引入的校准误差，校准精度下降。传统的低频振动校准方法不能满足宽频率范围内高精度的低频振动传感器及测量仪校准需求。基于单目视觉的测量方法已广泛用于许多精密测量领域，通过单目视觉方法的高精度位移测量可改善低频振动校准精度。

因此，针对目前低频振动校准方法对于低频振动传感器及测量仪校准的校准精度有限、过程繁琐、系统成本高及校准频率范围受限等不足，本发明提出一种可适用于宽频率范围、过程简单、成本低、校准精度高的低频振动校准方法。

技术实现要素：

针对目前低频振动校准方法无法适用于宽频率范围、精度有限、过程复杂及成本高等不足，本发明实施实例提供一种高精度的低频振动校准方法，包括：

(1)基于径向畸变的非线性模型摄像机标定：通过棋盘格靶标图像的x角点亚像素检测，实现高精度的摄像机标定；

(2)平面运动序列图像的增强与特征边缘提取：用于序列图像特征边缘的准确提取；包括：基于高斯曲线拟合特征边缘邻域灰度梯度，基于梯度光流法检测模糊图像运动方向，实现不同运动方向模糊图像的增强，基于zernike矩方法的增强序列图像亚像素特征边缘提取；

(3)特征边缘的空间运动位移测量与输出电压信号测量：通过摄像机标定确定的矩阵h与提取的序列图像特征边缘计算特征边缘的空间运动位移，并计算低频振动传感器及低频振动测量仪的输出电压信号峰值；

(4)低频振动传感器及低频振动测量仪的校准：利用测量的序列图像的特征边缘的空间运动位移计算相应的输入激励加速度峰值，通过输出电压信号峰值与激励加速度峰值，确定其灵敏度与幅频特性。

一种基于单目视觉的低频振动校准方法，所述校准方法包括以下步骤，

s1：利用棋盘格靶标图像x角点的亚像素坐标检测，实现基于径向畸变的非线性模型摄像机标定；

s2：针对运动产生的图像模糊，基于特征边缘邻域灰度梯度的高斯曲线拟合与基于梯度光流法的相邻两帧图像运动方向检测，实现不同运动方向图像的增强，利用基于zernike矩的方法提取增强图像的亚像素特征边缘；

s3：利用标定的摄像机模型参数及提取的序列图像特征边缘计算特征边缘的空间运动位移，通过二次微分获得被校低频振动传感器及低频振动测量仪的输入激励加速度峰值，并获取低频振动传感器及测量仪的输出电压峰值；

s4：最后通过获取的输入激励加速度峰值与输出电压峰值，确定低频振动传感器及测量仪的灵敏度与幅频特性。

所述摄像机标定用于确定摄像机的模型参数，具体包括：

(1)棋盘格靶标图像的x角点检测

针对采集的棋盘格靶标图像，利用自动x角点检测方法检测棋盘格靶标图像的x角点亚像素坐标(xd,yd)；

(2)径向畸变的非线性模型摄像机标定

选用基于径向畸变的非线性摄像机模型，则理想像点(xu,yu)与实际像点(xd,yd)满足如下公式：

其中，k1与k2为径向畸变系数。利用(xd,yd)与相应世界坐标(xw,yw)确定基于线性模型的摄像机参数，(x'w,y'w)为(xd,yd)的重投影世界坐标。通过如下公式可求解k1与k2。

其中，r与c为x角点阵列的行与列。利用求解的k1与k2校正(xd,yd)可得到(xu,yu)，以确定无畸变像点与世界坐标对应的摄像机模型参数h。

所述平面运动序列图像的增强与特征边缘提取，包括：基于高斯曲线拟合canny算子检测的像素级边缘邻域灰度梯度，高斯函数梯度拟合的计算如下：

其中，g(p)与x(p)分别为所选择邻域像素的灰度梯度与像素的横坐标，p为选择的像素数，a、μ及σ分别为拟合峰值、均值及标准差。当σ大于非模糊图像梯度高斯函数拟合的σt，利用梯度光流法检测其运动方向。若f(x,y)的特征边缘位置像素的横坐标相对于前一帧沿垂直方向减少，则其增强为：

否则，其增强为：

其中，fe(x,y)为增强图像，fmax(x,y)与fmin(x,y)分别为最大与最小灰度值，为f(x,y)的归一化灰度，t1与t2为两不同的阈值。

对于所述的增强序列图像，利用三灰度边缘模型的zernike矩方法实现其亚像素级特征边缘的提取。消除kxk正方形zernike矩模板的放大效应，像素级边缘点(x0,y0)的亚像素坐标为：

其中，d1与d2及φ为计算的距离及旋转角度边缘参数。

所述特征边缘随时间的空间运动位移为正弦，通过确定的摄像机模型参数h与提取的序列图像特征边缘计算特征边缘的空间运动位移，再基于正弦逼近法拟合该空间运动位移，以获得相应的位移峰值。

同理，利用正弦逼近法拟合被校低频振动传感器及低频振动测量仪的输出电压，获得相应的电压峰值。

通过所述特征边缘的位移峰值计算被校低频振动传感器及低频振动测量仪的输入激励加速度峰值，被校低频振动传感器及测量仪的灵敏度为输出电压信号峰值与激励加速度峰值之比。被校低频振动传感器及测量仪在不同频率下的灵敏度为其幅频特性。

低频振动校准方法的校准装置，该装置主要包括：水平长冲程振动台1、特征标志2、低频振动传感器及低频振动测量仪3、光源4、摄像机固定装置5、摄像机6、图像传输设备7、电压信号采集与传输设备8、处理与显示设备9。

水平长冲程振动台1用于提供低频振动传感器及低频振动测量仪3的输入激励；特征标志2与低频振动传感器及低频振动测量仪3紧固于水平长冲程振动台1的工作台面；光源4为摄像机6提供照明；摄像机固定装置5用于固定摄像机6，使其光轴垂直于水平长冲程振动台1的工作台面；摄像机6用于采集运动工作台面的序列图像；图像传输设备7传输序列图像；电压信号采集与传输设备8用于低频振动传感器及低频振动测量仪3输出电压信号的采集与传输；处理与显示设备9处理序列图像与电压信号，保存与显示校准结果。

本发明低频振动校准方法具有如下优势：

⑴本发明方法稳定、可靠、实用，可同时适用于多个不同类型的低频振动传感器及低频振动测量仪的校准。

⑵本发明方法校准过程简单、系统成本低，对于低频率范围内振动传感器及测量仪的校准只需一个工业摄像机。

⑶本发明方法采用棋盘格靶标x角点的亚像素定位与径向畸变的非线性模型保证摄像机的标定精度。

⑷本发明方法通过对运动序列图像的模糊判断，采用不同的增强函数对不同运动方向模糊图像的增强，保证特征边缘的准确提取。基于zernike矩方法实现增强序列图像的特征边缘亚像素提取。

⑸本发明方法属于低频振动校准方法，能够实现宽频率范围的高精度低频振动传感器及测量仪的校准。

附图说明

图1为本发明方法具体实施实例安装装置示意图；

图2为一种基于单目视觉的低频振动校准方法流程图；

图3为基于径向畸变的非线性模型摄像机标定流程图；

图4为平面运动序列图像的增强与基于zernike矩方法的特征边缘提取流程图；

图5-7为本发明方法具体实施实例对低频三轴向加速度传感器的校准结果图。

具体实施方式

为了解决现有低频振动校准方法的校准精度有限、系统复杂且成本高、不适用于宽频率范围高精度校准的问题，本发明提供了一种基于单目视觉的低频振动校准方法，本发明方法对于宽频率范围的低频振动校准可获得较高的精度，下面结合附图和具体的实施实例对本发明做出详细描述。

参考图1为本发明方法的实施实例装置示意图，该装置主要包括：水平长冲程振动台(1)、特征标志(2)、低频振动传感器及低频振动测量仪(3)、光源(4)、摄像机固定装置(5)、摄像机(6)、图像传输设备(7)、电压信号采集与传输设备(8)、处理与显示设备(9)。水平长冲程振动台(1)用于提供低频振动传感器及低频振动测量仪(3)的输入激励；特征标志(2)与低频振动传感器及低频振动测量仪(3)紧固于水平长冲程振动台(1)的工作台面；光源(4)为摄像机(6)提供照明；摄像机固定装置(5)用于固定摄像机(6)，使其光轴垂直于水平长冲程振动台(1)的工作台面；摄像机(6)用于采集运动工作台面的序列图像；图像传输设备(7)传输序列图像；电压信号采集与传输设备(8)用于低频振动传感器及低频振动测量仪(3)输出信号的采集与传输；处理与显示设备(9)处理序列图像与电压信号，保存与显示校准结果。

参考图2为一种基于单目视觉的低频振动校准方法流程图。本发明低频振动校准方法主要包括以下步骤：

步骤s20：基于径向畸变的非线性模型摄像机标定；

步骤s40：平面运动序列图像的增强与特征边缘提取，其包括：基于高斯函数拟合特征边缘邻域灰度梯度与基于梯度光流法的模糊图像运动方向检测，实现不同运动方向模糊图像的增强，基于zernike矩方法的增强序列图像亚像素特征边缘提取；

步骤s60：特征边缘的空间运动位移测量与输出电压信号测量，其包括：计算提取序列图像特征边缘的空间运动位移，以获取低频振动传感器及低频振动测量仪的输入激励加速度峰值，并计算对应的输出电压信号峰值；

步骤s80：低频振动传感器及低频振动测量仪的灵敏度与幅频特性：其包括：利用低频振动传感器及测量仪的输出电压信号峰值与输入激励加速度峰值计算其灵敏度，获取其幅频特性。

参考图3为基于径向畸变的非线性模型摄像机标定流程图。本发明非线性模型的摄像机标定包括如下步骤：

步骤s21：读入棋盘格靶标图像；

步骤s22：利用自动x角点检测方法检测棋盘格靶标图像x角点的亚像素坐标；

步骤s23：通过检测的x角点亚像素坐标与相应的世界坐标实现基于线性模型的摄像机标定；

步骤s24：基于最小二乘法优化x角点的世界坐标与其重投影世界坐标间距离最小的目标函数，求解径向畸变系数；

步骤s25：利用求解的畸变系数校正x角点的坐标，使其与相应的世界坐标满足线性关系；

步骤s26：通过校正的x角点坐标与世界坐标实现线性模型的摄像机标定，确定无畸变像点与世界坐标的摄像机模型参数h。

参考图4为平面运动序列图像的增强与基于zernike矩方法的特征边缘提取流程图。本发明序列图像的增强与基于zernike矩方法的特征边缘提取包括如下步骤：

步骤s41：读入运动工作台面的序列图像；

步骤s42：基于canny算子的特征边缘检测；

步骤s43：选取特征边缘邻域内对称分布的50个像素点的灰度梯度，通过高斯曲线拟合计算其高斯核；

步骤s44：判断灰度梯度拟合的高斯核σ是否大于阈值σt(σt＝2)，若满足该条件，则跳至步骤s45，否则跳至步骤s49；

步骤s45：基于梯度光流法检测当前帧与上一帧图像的梯度光流场，依据光流场判断特征边缘的运动方向；

步骤s46：特征边缘位置像素的横坐标是否增加，若满足条件，则跳至步骤s47，否则跳至步骤s48；

步骤s47：利用边缘位置像素增加方向的增强函数实现模糊图像增强；

步骤s48：利用边缘位置像素减小方向的增强函数实现模糊图像增强；

步骤s49：计算增强图像的不同阶次zernike矩；

步骤s50：通过不同阶次的zernike矩计算特征边缘像素点的距离及旋转角度边缘参数；

步骤s51：消除zernike矩模板的放大效应，利用计算的边缘参数获得像素级特征边缘的亚像素坐标。

本实施实例装置的具体参数为：频率范围为0.01-200hz、最大峰-峰位移为400mm的水平长冲程振动台，特征标志选用与振动台工作台面呈高对比度的矩形金属板，被校低频振动传感器选用msv3000-02三轴向加速度传感器，分辨率为1292x964像素、最大帧率为30fps的avtmantag-125b工业摄像机，镜头焦距为8mm，光源选用60w白炽灯，电压信号采集与传输设备选用采样频率范围为1hz-216khz的inv3062c。

为了验证本发明低频振动校准方法的校准精度，利用本发明校准方法实现了0.04-8hz范围内的加速度传感器校准。同时使用激光干涉法与地球重力法对该加速度传感器在相同频率范围内进行校准。在0.04-2hz范围内可选择以地球重力法的校准结果作为参考，在0.2-8hz范围内可选择以激光干涉法的校准结果作为参考。

参考图5-7分别为本发明方法具体实施实例对三轴向加速度传感器的x、y、z轴向的灵敏度校准结果图，本次校准的长冲程振动台提供的最大峰-峰位移为360mm。由结果图可知，本发明低频振动校准方法与地球重力法在0.04-2hz范围内的校准结果相比，其x、y、z轴向灵敏度的相对误差分别小于1.29％、1.75％、1.23％；本发明低频振动校准方法与激光干涉法在2-8hz范围内的最大相对误差分别为0.48％、0.53％、0.36，说明本发明针对宽频率范围内的低频振动校准具有较高的校准精度。

上述描述为本发明实施实例的详细介绍，其并非用于对本发明作任何形式上的限定。本领域相关技术人员可在本发明的基础上可做出一系列的优化、改进及修改等。因此，本发明的保护范围应由所附权利要求来限定。