CN109544639A - 一种多镜面单相机三维振动测试装置及方法 - Google Patents
一种多镜面单相机三维振动测试装置及方法 Download PDFInfo
- Publication number
- CN109544639A CN109544639A CN201811274673.0A CN201811274673A CN109544639A CN 109544639 A CN109544639 A CN 109544639A CN 201811274673 A CN201811274673 A CN 201811274673A CN 109544639 A CN109544639 A CN 109544639A
- Authority
- CN
- China
- Prior art keywords
- camera
- image
- coordinate
- point
- mirror
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种多镜面单相机三维振动测试装置及方法,测试方法包括:采用平面标定板标定相机内部参数;在物体表面每个待测点位置粘贴标识点,并在每个标识点旁边固定反射镜;采用激励装置让物体发生振动,相机采集整个振动过程中的所有标识点实像与虚像;利用相机内部参数及每个标识点实像与虚像的对应图像坐标,求解加载过程中任意状态的相机坐标系与每个镜面坐标系的旋转矩阵及归一化平移向量,采用反射镜中任意两点之间的实际距离计算平移向量系数;计算物体表面待测点在每一时刻的三维空间坐标,即可测试物体振动中的频率及模态等振动特征。本发明利用单个相机及多个镜面实现多点式三维振动测试,具有非接触、测量精度高、方便实用等优点。
Description
技术领域
本发明涉及振动测试技术,具体涉及一种多镜面单相机三维振动测试装置及方法。
背景技术
接触式传感器由于需要与被测物接触且有附加质量限制了其在振动过程中的应用。常用的基于双相机的三维坐标测试装置则需要相机同步装置,额外增加了硬件成本;并且双相机的标定过程非常复杂,限制了其在三维坐标测试中的使用。基于单相机的单镜面光学测试方法由于相机景深的问题无法测量大型物体的三维坐标及位移。
发明内容
本发明的目的在于提供一种多镜面单相机三维振动测试装置及方法,克服了单镜面大视场测试中的相机景深问题,操作简单灵活,易于实现。
实现本发明目的的技术解决方案为:一种多镜面单相机三维振动测试装置,包括一台相机、高分辨率镜头、计算机、与待测点数目相同的标识点、与待测点数目相同的前镀膜反射镜,其中每个反射镜上贴有两个标识点,相机采样频率为物体振动频率2倍以上;
采用激励装置使得待测物发生振动,相机采集整个振动状态图像,保存至计算机,得到的任意一幅图像中均包含所有待测点的实像与镜面反射的虚像两部分;
通过标识点识别与定位算法确定任意一幅图像中每个实像与虚像的对应图像坐标,利用整个振动过程中对应图像坐标、镜面上标识点图像坐标以及相机内部参数计算相机坐标系与每个镜面坐标系的外部参数;根据相机的内部参数及振动过程中的对应图像坐标及外部参数计算待测物体上每个标识点的三维坐标,并进行模态及频率振动分析。
一种多镜面单相机三维振动测试方法,包括以下步骤:
步骤1、相机内部参数标定:将平面标定板在相机的视场范围内平移及旋转至少8个姿态,相机采集平面标定板图像,利用标定板图像计算相机内部参数矩阵及畸变系数;
步骤2、标识点粘贴与反射镜固定:在待测位置粘贴标识点,并在其旁边固定贴有两个标识点的前镀膜反射镜;
步骤3、相机图像采集:采用激励装置让待测物发生振动,相机采集整个振动状态图像,保存至计算机,得到的任意一幅图像中均包含所有待测点的实像与镜面反射的虚像两部分;
步骤4、相机与每个镜面坐标系外部参数确定:通过标识点识别与定位算法确定任意一幅图像中每个实像与虚像的对应图像坐标,利用整个振动过程中对应图像坐标、镜面上标识点图像坐标以及相机内部参数计算相机坐标系与每个镜面坐标系的外部参数;
步骤5、待测物体标识点的三维坐标计算:根据相机的内部参数及振动过程中的对应图像坐标及外部参数计算待测物体上每个标识点的三维坐标,并进行模态及频率振动分析。
本发明与现有技术相比,其显著优点为:(1)相机的内部参数固定,本发明可以实现一次标定,多次测量的目的;(2)本发明采用单个相机与多镜面结合的方式,解决了单镜面引起的相机景深问题;(3)本发明克服了常规双目系统设备复杂,且需要复杂标定的缺点;(4)相比于传统接触式传感器,本发明具有无损、非接触、高精度等优点。
附图说明
图1为本发明测量装置示意图。
图2为本发明测量方法流程图。
图3为本发明测量的阻尼衰减振动过程示意图。
具体实施方式
下面结合附图和具体实施例对本发明方案作进一步的说明。
如图1所示,一种多镜面单相机三维振动测试装置,由一台工业相机4、高分辨率镜头3、三脚架5、计算机6、与待测点数目相同的标识点1、与待测点数目相同的前镀膜反射镜2组成,反射镜放置与待测点附近且待测点及其反射镜虚像能同时被相机观测。其中每个反射镜上贴有两个标识点,选择采样频率超过物体振动频率2倍的快速相机或高速相机。高分辨率镜头像素在500万以上。
如图2所示,一种多镜面单相机三维振动测试方法,包括以下步骤:
步骤1、相机内部参数标定:将平面标定板在相机的视场范围内平移及旋转至少8个姿态,相机采集平面标定板图像,利用标定板图像计算相机内部参数矩阵及畸变系数;得到的相机内部参数矩阵为其中fx、fy为两个方向的等效焦距,cx、cy为光轴与相机靶面的交点图像坐标,一阶畸变系数用k1表示。
步骤2、标识点粘贴与反射镜固定:在待测位置粘贴标识点,并在其旁边固定贴有两个标识点的前镀膜反射镜;每个待测点位置均粘贴一个标识点,且每个标识点旁边均固定有一前镀膜反射镜,在反射镜上贴有两个标识点,且两个标识点之间的距离已知为l。
步骤3、相机图像采集:采用激励装置让待测物发生振动,相机采集整个振动状态图像,保存至计算机,得到的任意一幅图像中均包含所有待测点的实像与镜面反射的虚像两部分;
步骤4、相机与每个镜面坐标系外部参数确定:通过标识点识别与定位算法确定任意一幅图像中每个实像与虚像的对应图像坐标,利用整个振动过程中对应图像坐标、镜面上标识点图像坐标以及相机内部参数计算相机坐标系与每个镜面坐标系的外部参数;采用整个振动过程中的对应图像坐标ui、vi和u′i、v′i,其中i表示任意一个振动状态,相机坐标系与世界坐标系的外部参数中旋转矩阵的计算方法如下:
在反射镜平面上建立世界坐标系,设垂直于激光平面的方向为世界坐标系的X轴,相机坐标系与世界坐标系之间的旋转矩阵为则通过如下公式进行Levenberg-Marquardt算法迭代优化,精确计算旋转矩阵R:
(yiux′iu-xiuy′iu)(R3R7-R4R6)+(x′iu-xiu)(R3R8-R5R6)+(y′iu-yiu)(R4R8-R5R7)=0
其中,归一化畸变图像坐标:xid=(ui-cx)/fx,yid=(vi-cy)/fy,xid′=(ui′-cx)/fx,yid′=(vi′-cy)/fy;归一化去畸变图像坐标:xiu=xid(1+k1R2),yiu=yid(1+k1R2),x′iu=x′id(1+k1R′2),y′iu=y′id(1+k1R′2),旋转矩阵R的其他分量根据其正交性求解。
相机坐标系与每个镜面坐标系的外部参数中平移向量的计算方法如下:
根据反射镜位置的不同确定归一化平移向量t′=[t′x t′y t′z]T,若反射镜竖直放置,则归一化平移向量为t′=[1 0 0]T;若反射镜水平放置,则归一化平移向量为t′=[01 0]T;假设反射镜上两点的图像坐标分别为u1、v1和u2、v2,通过如下公式计算这两点的三维坐标[Xi Yi Zi],其中i=1,2:
其中,xiu、yiu为归一化去畸变图像坐标。采用归一化平移向量计算反射镜上两点的距离平移向量t=l/l1×t′。
步骤5、待测物体标识点的三维坐标计算:根据相机的内部参数及振动过程中的对应图像坐标及外部参数计算待测物体上每个标识点的三维坐标,并进行模态及频率振动分析。
计算待测物体表面某一标识点三维坐标的方法为:
其中,ti=[txi tyi tzi]T为计算得到的对应相机与镜面坐标系的平移向量,对应相机与镜面坐标系的旋转矩阵i代表某特定标识点。实像的对应图像坐标为u、v,虚像的对应图像坐标为u′、v′。归一化畸变图像坐标:xd=(u-cx)/fx,yd=(v-cy)/fy,xd′=(u′-cx)/fx,yd′=(v′-cy)/fy;归一化去畸变图像坐标:xu=xd(1+k1R2),yu=yd(1+k1R2),x′u=x′d(1+k1R′2),y′u=y′d(1+k1R′2),
通过相机采样帧频及所测得的三维数据利用傅里叶变换即可进行频率分析,对全部标识点的三维数据分析可以得到振动过程中的模态信息;具体为:通过相机采样帧频及所测得的三维数据即可得到待测点的位移时程曲线,利用傅里叶变换方法即可分析振动频率,振动过程中的模态信息可以基于分析所有待测点的位移数据得到。
如图3所示,本发明的测量精度高,可以达到亚微米的位移测试量级,采用单个相机可以同时获得多个待测点的振动信息,采样帧频高,利用通用的高速相机已经可以实现2000帧以上的高速动态测试,避免了常规双相机三维测试时的相机同步等技术难题。
Claims (9)
1.一种多镜面单相机三维振动测试装置,其特征在于,包括一台相机、高分辨率镜头、计算机、与待测点数目相同的标识点、与待测点数目相同的前镀膜反射镜,其中每个反射镜上贴有两个标识点,相机采样频率为物体振动频率2倍以上;
采用激励装置使得待测物发生振动,相机采集整个振动状态图像,保存至计算机,得到的任意一幅图像中均包含所有待测点的实像与镜面反射的虚像两部分;
通过标识点识别与定位算法确定任意一幅图像中每个实像与虚像的对应图像坐标,利用整个振动过程中对应图像坐标、镜面上标识点图像坐标以及相机内部参数计算相机坐标系与每个镜面坐标系的外部参数;根据相机的内部参数及振动过程中的对应图像坐标及外部参数计算待测物体上每个标识点的三维坐标,并进行模态及频率振动分析。
2.根据权利要求1所述的多镜面单相机三维振动测试装置,其特征在于,测试装置还包括用于固定相机的三脚架。
3.一种基于权利要求1所述多镜面单相机三维振动测试装置的测试方法,其特征在于,包括以下步骤:
步骤1、相机内部参数标定:将平面标定板在相机的视场范围内平移及旋转至少8个姿态,相机采集平面标定板图像,利用标定板图像计算相机内部参数矩阵及畸变系数;
步骤2、标识点粘贴与反射镜固定:在待测位置粘贴标识点,并在其旁边固定贴有两个标识点的前镀膜反射镜;
步骤3、相机图像采集:采用激励装置让待测物发生振动,相机采集整个振动状态图像,保存至计算机,得到的任意一幅图像中均包含所有待测点的实像与镜面反射的虚像两部分;
步骤4、相机与每个镜面坐标系外部参数确定:通过标识点识别与定位算法确定任意一幅图像中每个实像与虚像的对应图像坐标,利用整个振动过程中对应图像坐标、镜面上标识点图像坐标以及相机内部参数计算相机坐标系与每个镜面坐标系的外部参数;
步骤5、待测物体标识点的三维坐标计算:根据相机的内部参数及振动过程中的对应图像坐标及外部参数计算待测物体上每个标识点的三维坐标,并进行模态及频率振动分析。
4.根据权利要求3所述的多镜面单相机三维振动测试方法,其特征在于,步骤1中得到的相机内部参数矩阵为其中fx、fy为两个方向的等效焦距,cx、cy为光轴与相机靶面的交点图像坐标,一阶畸变系数用k1表示。
5.根据权利要求3所述的多镜面单相机三维振动测试方法,其特征在于,步骤2中每个待测点位置均粘贴一个标识点,且每个标识点旁边均固定有一前镀膜反射镜,在反射镜上贴有两个标识点,且两个标识点之间的距离已知为l。
6.根据权利要求3所述的多镜面单相机三维振动测试方法,其特征在于,所述步骤4中采用整个振动过程中的对应图像坐标ui、vi和u'i、v'i,其中i表示任意一个振动状态,相机坐标系与世界坐标系的外部参数中旋转矩阵的计算方法如下:
在反射镜平面上建立世界坐标系,设垂直于激光平面的方向为世界坐标系的X轴,相机坐标系与世界坐标系之间的旋转矩阵为则通过如下公式进行Levenberg-Marquardt算法迭代优化,精确计算旋转矩阵R:
(yiux'iu-xiuy'iu)(R3R7-R4R6)+(x'iu-xiu)(R3R8-R5R6)+(y'iu-yiu)(R4R8-R5R7)=0
其中,归一化畸变图像坐标:xid=(ui-cx)/fx,yid=(vi-cy)/fy,xid'=(ui'-cx)/fx,yid'=(vi'-cy)/fy;归一化去畸变图像坐标:xiu=xid(1+k1R2),yiu=yid(1+k1R2),x'iu=x'id(1+k1R'2),y'iu=y'id(1+k1R'2),旋转矩阵R的其他分量根据其正交性求解。
7.根据权利要求3所述的多镜面单相机三维振动测试方法,其特征在于,步骤4中相机坐标系与每个镜面坐标系的外部参数中平移向量的计算方法如下:
根据反射镜位置的不同确定归一化平移向量t'=[t'x t'y t'z]T,若反射镜竖直放置,则归一化平移向量为t'=[1 0 0]T;若反射镜水平放置,则归一化平移向量为t'=[0 1 0]T;假设反射镜上两点的图像坐标分别为u1、v1和u2、v2,通过如下公式计算这两点的三维坐标[Xi Yi Zi],其中i=1,2:
其中,xiu、yiu为归一化去畸变图像坐标;采用归一化平移向量计算反射镜上两点的距离平移向量t=l/l1×t'。
8.根据权利要求3所述的单相机三维视频引伸计测量方法,其特征在于,所述步骤5中计算待测物体表面某一标识点三维坐标的方法为:
其中,ti=[txi tyi tzi]T为计算得到的对应相机与镜面坐标系的平移向量,对应相机与镜面坐标系的旋转矩阵i代表某特定标识点。实像的对应图像坐标为u、v,虚像的对应图像坐标为u'、v'。归一化畸变图像坐标:xd=(u-cx)/fx,yd=(v-cy)/fy,xd'=(u'-cx)/fx,yd'=(v'-cy)/fy;归一化去畸变图像坐标:xu=xd(1+k1R2),yu=yd(1+k1R2),x'u=x'd(1+k1R'2),y'u=y'd(1+k1R'2),
9.根据权利要求3所述的单相机三维视频引伸计测量方法,其特征在于,所述步骤5中频率及模态振动分析方法为:
通过相机采样帧频及所测得的三维数据利用傅里叶变换即可进行频率分析,对全部标识点的三维数据分析可以得到振动过程中的模态信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811274673.0A CN109544639B (zh) | 2018-10-30 | 2018-10-30 | 一种多镜面单相机三维振动测试装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811274673.0A CN109544639B (zh) | 2018-10-30 | 2018-10-30 | 一种多镜面单相机三维振动测试装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109544639A true CN109544639A (zh) | 2019-03-29 |
CN109544639B CN109544639B (zh) | 2022-08-16 |
Family
ID=65845976
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811274673.0A Active CN109544639B (zh) | 2018-10-30 | 2018-10-30 | 一种多镜面单相机三维振动测试装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109544639B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113375584A (zh) * | 2021-06-03 | 2021-09-10 | 深圳市海塞姆科技有限公司 | 单目三维采集结构及单目三维多标距视觉系统、引伸计 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107014312A (zh) * | 2017-04-25 | 2017-08-04 | 西安交通大学 | 一种振镜式线激光扫描三维测量系统的整体标定方法 |
CN107687821A (zh) * | 2017-08-24 | 2018-02-13 | 南京理工大学 | 用于深孔洞旋转件形貌测量的多相机光刀系统标定方法 |
-
2018
- 2018-10-30 CN CN201811274673.0A patent/CN109544639B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107014312A (zh) * | 2017-04-25 | 2017-08-04 | 西安交通大学 | 一种振镜式线激光扫描三维测量系统的整体标定方法 |
CN107687821A (zh) * | 2017-08-24 | 2018-02-13 | 南京理工大学 | 用于深孔洞旋转件形貌测量的多相机光刀系统标定方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113375584A (zh) * | 2021-06-03 | 2021-09-10 | 深圳市海塞姆科技有限公司 | 单目三维采集结构及单目三维多标距视觉系统、引伸计 |
Also Published As
Publication number | Publication date |
---|---|
CN109544639B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8497901B2 (en) | Method and device for exact measurement of objects | |
CN104990515B (zh) | 大型物体三维形状测量系统及其测量方法 | |
CN106989812A (zh) | 基于摄影测量技术的大型风机叶片模态测试方法 | |
EP2538241A2 (en) | Advanced remote nondestructive inspection system and process | |
CN109099883A (zh) | 高精度大视场机器视觉测量与标定装置及方法 | |
Pan et al. | A novel mirror-assisted multi-view digital image correlation for dual-surface shape and deformation measurements of sheet samples | |
CN107941168B (zh) | 基于散斑位置标定的反射式条纹面形测量方法与装置 | |
CN112540089B (zh) | 数字成像系统在混凝土桥梁裂缝检测和分析中的应用方法 | |
US20140174149A1 (en) | Calibration method for trigonometric-based ranging systems in multiple media | |
CN108195314B (zh) | 基于多视场拼接的反射式条纹三维面形测量方法 | |
CN106885622A (zh) | 一种大视场多点三维振动测量方法 | |
US10401145B2 (en) | Method for calibrating an optical arrangement | |
CN110806182A (zh) | 基于远心镜头的高精度光学引伸计及测量方法 | |
Chen et al. | Mirror-assisted multi-view digital image correlation with improved spatial resolution | |
CN110139094A (zh) | 一种光心对齐方法、光心坐标计算方法及装置 | |
CN109556528B (zh) | 一种单相机三维视频引伸计及测量方法 | |
CN110044266B (zh) | 基于散斑投影的摄影测量系统 | |
WO2010013289A1 (ja) | カメラ校正画像作成装置およびカメラ校正画像作成プログラム | |
CN109544639A (zh) | 一种多镜面单相机三维振动测试装置及方法 | |
CN110763146A (zh) | 基于双相机的高精度光学引伸计及测量方法 | |
CN114062265A (zh) | 一种用于视觉系统的支撑结构稳定性的评估方法 | |
CN110146032A (zh) | 基于光场分布的合成孔径相机标定方法 | |
CN105783735A (zh) | 一种基于二维数字图像相关补偿算法的实时引伸计测量方法 | |
CN113899528A (zh) | 一种基于3d-dic原理的液面动态波高测量方法 | |
CN110702505B (zh) | 一种基于远心镜头和立方棱镜的双视场视频引伸计 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |