一种标准的虚拟正弦直线振动测量方法技术

本发明专利技术公开了一种标准的虚拟正弦直线振动测量方法，包括计算图像的灰度梯度幅值，利用灰度梯度幅值的分布寻找区间内的幅值极大值，根据极大值的个数判断图像中的感兴趣区域是否存在不同运动位置目标的叠加；针对感兴趣区域的不同叠加情形，以提高后续图像处理的精度；使用一种亚像素边缘检测方法实现运动目标特征边缘的提取；引入目标振动模型，利用最小二乘法进行模型拟合；最后对目标振动模型进行位移误差补偿，获得虚拟正弦直线振动的时间‑位移测量曲线。本发明专利技术可避免因机械制造、运动控制等因素导致的振动台非理想正弦振动情形，通过Matlab编程输出标准的虚拟正弦直线振动，为基于机器视觉的平面运动测量提供了一种虚拟溯源途径。

A Standard Virtual Sinusoidal Linear Vibration Measurement Method

【技术实现步骤摘要】
一种标准的虚拟正弦直线振动测量方法
本专利技术属于图像处理与机器视觉检测领域，尤其涉及一种标准的虚拟正弦直线振动测量方法。
技术介绍
振动传感器及测量仪被广泛用于多个领域的实时监控与振动参数测量，如桥梁建筑、地震、汽车、航空航天等。为了保证振动传感器及测量仪所得测量结果的可靠性与准确性，需要定期对振动传感器及测量仪进行校准。典型的振动传感器及测量仪校准方法有激光干涉绝对法校准、比较法校准等。这些校准法使用不同的测量方法获得振动台的振动曲线图，通过比较这些测量结果来实现振动传感器的校准。而因振动台存在机械加工、运动控制等误差影响，振动台的振动形式并非标准的正弦直线振动，导致现有校准方法存在较大误差。本专利技术方法通过Matlab编程产生标准的虚拟正弦直线振动，通过机器视觉实现对标准的虚拟正弦直线振动的测量，避免了振动台引入的误差，进一步提高了振动测量精度。除此之外，本专利技术方法为平面振动测量方法提供了一种虚拟溯源途径。
技术实现思路
为了避免因机械加工、运动控制等外界因素引入的误差，使得振动台的振动形式为非标准正弦直线振动，导致振动测量精度不高，本专利技术提出一种标准的虚拟正弦直线振动测量方法，包括：图像有无目标叠加区域的判断：用于判断读取的序列图像中是否存在目标叠加区域，包括：图像中运动目标的感兴趣区域选择，感兴趣区域有无目标叠加区域的判断；图像复原方法的选择：用于提高后续图像处理的精度，包括：目标运动方向的判断，基于运动方向的图像复原方法的选择；运动目标特征边缘的高精度检测：用于提取运动目标特征边缘的亚像素检测，以提高特征边缘精度，保证运动测量结果的可靠性；目标振动模型的拟合：引入目标振动模型，使用最小二乘法进行模型拟合；目标振动模型的修正：对拟合得到的振动模型进行位移误差补偿，得到虚拟正弦直线振动测量关于位移-时间的最佳曲线。一种标准的虚拟正弦直线振动测量方法，包括以下步骤：S1：计算图像的灰度梯度幅值，利用灰度梯度幅值的分布寻找区域内的峰值，根据峰值个数判断图像感兴趣区域内是否存在运动目标叠加；S2：针对区域存在运动目标叠加的图像，利用光流梯度法判断运动目标的运动方向，采用不同的图像复原方法复原图像，以提高后续图像处理的精度；S3：利用亚像素边缘检测方法完成运动目标特征边缘的提取；S4：引入目标振动模型，使用最小二乘法进行模型拟合；S5：对目标振动模型进行位移误差补偿，得到虚拟正弦直线振动的时间-位移测量曲线。所述步骤S1具体包括：(1)运动目标感兴趣区域的确定使用帧间差分法，确定运动目标的运动区域；为了降低噪声的影响，对获得的运动区域进行连通域分析；计算每个连通域的面积，将面积最大的两个区域作为运动目标的感兴趣区域；(2)运动目标叠加区域的判断在运动目标的感兴趣区域，计算灰度梯度幅值：其中，f(x,y)为像素(x,y)处的灰度值；Ix和Iy分别为x方向和y方向上的灰度梯度；▽f(x,y)为像素(x,y)处的梯度矢量；mag(▽f)为梯度矢量对应的幅值；当区域无运动目标叠加时，遍历区域的梯度幅值，只存在两个局部极大值点；当区域存在运动目标叠加时，局部极大值点数量大于两个。所述步骤S2具体包括：(1)目标运动方向的确定利用光流梯度法计算运动目标的感兴趣区域光流，该算法假定相邻两帧图像对应像素的亮度恒定，计算公式如下：Ixu+Iyv+It＝0(2)其中，It为图像在t时刻的灰度梯度；利用x方向上的光流u和y方向上的光流v判断目标运动方向；(2)图像复原方法针对背景与运动目标灰度对比强烈的图像，其中运动目标为实心矩形；假定运动目标特征边缘始终为其右边缘，且当相邻帧图像运动目标特征边缘的x坐标呈增大趋势时认为目标运动方向为正方向；对目标特征边缘领域进行灰度区域划分；针对运动方向为正方向的感兴趣区域实行最右区域灰度增强，其余区域灰度减弱；针对运动方向为负方向的感兴趣区域实行最左区域灰度增强，其余区域灰度减弱。所述步骤S3具体包括：亚像素边缘检测通过Zernike矩法实现，图像f(x,y)的n阶m次Zernike矩定义为：其中，是在极坐标系下的单位圆内的正交n阶m次Zernike多项式；*表示复共轭；A'nm为图像旋转φ角度后的n阶m次Zernike矩；通过3个不同阶次的Zernike矩实现亚像素边缘检测，分别为A00、A11、A20，其对应的积分核函数为：根据公式(3)、(4)得：图像运动目标特征边缘的亚像素位置为：所述步骤S4具体包括：引入目标振动模型：利用最小二乘法进行模型拟合。所述步骤S5具体包括：图像尺寸以像素为单位，故相机获取的虚拟正弦直线振动序列图像存在量化误差；显示器输出位置sd(ti)与真实位置s(ti)满足如下公式：sd(ti)＝s(ti)+Δs(ti)(8)为了提高测量精度，需要对目标振动模型进行位移误差补偿，公式如下：Δsc(tj)＝int[sp-psin(ωv(NTtj)-π/2)]-sp-psin(ωv(NTtj)-π/2)(9)其中，sp-psin(ωv(NTtj)-π/2)表示为：振动模型在tj时刻对应的位移值；根据公式(8)、(9)可得任意位置sd(ti)的修正值：s'd(ti)＝sd(ti)-Δsc(tj)(10)对目标振动模型修正后，获得虚拟正弦直线振动关于时间-位移的测量曲线。本专利技术虚拟正弦直线振动测量方法具有如下优势：(1)本专利技术通过图像复原和模型修正方法提高了振动测量精度。(2)本专利技术采用标准的虚拟正弦直线振动，可避免振动台引入的误差。(3)本专利技术方法为平面振动测量方法提供了一种虚拟溯源途径。附图说明图1为本专利技术方法装置示意图；图2为一种标准的虚拟正弦直线振动测量方法流程图；图3为判断图像有无目标叠加区域的流程图；图4为基于目标运动方向的图像复原方法选取流程图；图5为不同运动方向、不同灰度分布下的运动目标特征边缘位置示意图；图6为摄像机采集的标准虚拟正弦直线振动序列图像；图7为本专利技术方法测量标准虚拟正弦直线振动得到的时间-位移结果图(频率为1Hz)；图8-11为本专利技术方法测量不同频率下的标准虚拟正弦直线振动得到的位移、相位的均值绝对值柱状图和均方差柱状图。具体实施方式为了避免振动台引入的误差，本专利技术提供了一种标准的虚拟正弦直线振动测量方法，本专利技术方法对于平面运动的测量具有较高的精度，下面结合附图和具体的实施实例对本专利技术做出详细描述。参考图1为本专利技术方法的装置示意图，该装置主要包括：光源1、摄像机2、处理及显示设备3。光源1为摄像机2提供照明；摄像机2用于采集标准的虚拟正弦直线振动的序列图像；处理及显示设备3用于呈现Matlab生成的标准虚拟正弦直线振动以及实现序列图像的处理。参考图2为一种标准的虚拟正弦直线振动测量方法流程图。本专利技术虚拟正弦直线振动测量方法主要包括以下步骤：步骤S160：通过Matlab编程生成标准的虚拟正弦直线振动；步骤S180：通过摄像机获取标准的虚拟正弦直线振动序列图像；步骤S200：读入虚拟正弦直线振动序列图像；步骤S220：图像有无目标叠加区域的判断，其包括：图像中运动目标的感兴趣区域检测，感兴趣区域有无目标叠加区域的判断；步骤S240：图像复原方法的选择，其包括：运动目标运动方向的判断，基于运动方向的图像复原方法的选择；步骤S260：序列图本文档来自技高网...

【技术保护点】
1.一种标准的虚拟正弦直线振动测量方法，其特征在于，包括以下步骤：S1：计算图像的灰度梯度幅值，利用灰度梯度幅值的分布寻找区域内的峰值，根据峰值个数判断图像感兴趣区域内是否存在运动目标叠加；S2：针对区域存在运动目标叠加的图像，利用光流梯度法判断运动目标的运动方向，采用不同的图像复原方法复原图像，以提高后续图像处理的精度；S3：利用亚像素边缘检测方法完成运动目标特征边缘的提取；S4：引入目标振动模型，使用最小二乘法进行模型拟合；S5：对目标振动模型进行位移误差补偿，得到虚拟正弦直线振动的时间‑位移测量曲线。

【技术特征摘要】
1.一种标准的虚拟正弦直线振动测量方法，其特征在于，包括以下步骤：S1：计算图像的灰度梯度幅值，利用灰度梯度幅值的分布寻找区域内的峰值，根据峰值个数判断图像感兴趣区域内是否存在运动目标叠加；S2：针对区域存在运动目标叠加的图像，利用光流梯度法判断运动目标的运动方向，采用不同的图像复原方法复原图像，以提高后续图像处理的精度；S3：利用亚像素边缘检测方法完成运动目标特征边缘的提取；S4：引入目标振动模型，使用最小二乘法进行模型拟合；S5：对目标振动模型进行位移误差补偿，得到虚拟正弦直线振动的时间-位移测量曲线。2.根据权利要求1所述的一种标准的虚拟正弦直线振动测量方法，其特征在于：所述步骤S1具体包括：(1)运动目标感兴趣区域的确定使用帧间差分法，确定运动目标的运动区域；为了降低噪声的影响，对获得的运动区域进行连通域分析；计算每个连通域的面积，将面积最大的两个区域作为运动目标的感兴趣区域；(2)运动目标叠加区域的判断在运动目标的感兴趣区域，计算灰度梯度幅值：其中，f(x,y)为像素(x,y)处的灰度值；Ix和Iy分别为x方向和y方向上的灰度梯度；为像素(x,y)处的梯度矢量；为梯度矢量对应的幅值；当区域无运动目标叠加时，遍历区域的梯度幅值，只存在两个局部极大值点；当区域存在运动目标叠加时，局部极大值点数量大于两个。3.根据权利要求1所述的一种标准的虚拟正弦直线振动测量方法，其特征在于：所述步骤S2具体包括：(1)目标运动方向的确定利用光流梯度法计算运动目标的感兴趣区域光流，该算法假定相邻两帧图像对应像素的亮度恒定，计算公式如下：Ixu+Iyv+It＝0(2)其中，It为图像在t时刻的灰度梯度；利用x方向上的光流u和y方向上的光流v判断目标运动方向；(2)图像复原方法针对背景与运动目标灰度对比强烈的图像，其中运动目标为实心矩形；假定运动目标特征边缘始终为其右边缘，且当相邻...

【专利技术属性】
技术研发人员：蔡晨光，张颖，刘志华，
申请(专利权)人：中国计量科学研究院，
类型：发明
国别省市：北京,11