CN106225668B - 基于多折射模型的风洞投放物高速位姿测量方法 - Google Patents
基于多折射模型的风洞投放物高速位姿测量方法 Download PDFInfo
- Publication number
- CN106225668B CN106225668B CN201610594960.4A CN201610594960A CN106225668B CN 106225668 B CN106225668 B CN 106225668B CN 201610594960 A CN201610594960 A CN 201610594960A CN 106225668 B CN106225668 B CN 106225668B
- Authority
- CN
- China
- Prior art keywords
- refraction
- wind
- tunnel
- missile
- models
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 238000002474 experimental method Methods 0.000 claims abstract description 9
- 238000000691 measurement method Methods 0.000 claims abstract description 4
- 238000011426 transformation method Methods 0.000 claims abstract 2
- 230000000007 visual effect Effects 0.000 claims description 11
- 238000013459 approach Methods 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 230000003287 optical effect Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 2
- 238000005096 rolling process Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 239000002828 fuel tank Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
- G01B11/005—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates coordinate measuring machines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明基于多折射模型的风洞投放物高速位姿测量方法属于计算机视觉测量技术领域，涉及一种适用于风洞多折射环境下投放物高速位姿测量方法。测量方法首先建立风洞实验环境下的多折射模型；然后利用多折射模型对视觉测量图像数据进行处理优化，得到左、右图像特征信息的正确匹配；最终利用坐标转换方法求的投放物位姿测量信息。该方法对风洞环境高速气流以及观察窗引起的光折射现象进行分析，引入了风洞实验环境下的光路折射模型，有效消除了观察窗、高速气流对光线传播的影响，建立更为精确的多折射模型对位姿测量过程进行校正优化，提高了存在多折射现象的风洞投放物高速位姿测量精度与稳定性。
Description
技术领域
本发明属于计算机视觉测量技术领域，涉及一种适用于风洞多折射环境下投放物高速位姿测量方法。
背景技术
风洞实验是航天航空领域中应用最为广泛的实验方法，投放物位姿测量是风洞实验中极其重要的实验内容。随着现代科技的不断进步，对投放物位姿测量精度、速度提出了全新的要求。针对这种测量需求，视觉测量由于其高频响、高精度成为非接触式测量中的首选。但由于风洞高速气流会使视觉设备震动，因此视觉设备多架设于风洞试验段外通过观察窗进行测量，这样势必引入观察窗与风洞高速气流的光折射现象，这将大大降低位姿测量的精度。因此，在存在光路多折射的情况下，获得高精度的位姿测量数据是现阶段要解决的主要难题。
目前，风洞位姿测量试验中考虑光折射问题研究相对较少，目前，大连理工大学刘巍等人2015年在航空学报第5期36卷发表的《基于彩色编码的副油箱风洞模型位姿测量方法》中提出了基于彩色编码图像的副油箱位姿测量方法，利用彩色编码的自发光标记点进行图像匹配，解决了低照度条件下目标标记亮度不足和由于物体滚转造成的标记点消隐现象。但此方法没有考虑观察视窗给光路带来的折射问题，系统精度会有所降低。大连理工大学刘巍等人申请的发明专利公开号为ZL 201310139656.7，“一种高速滚转体位姿测量方法”中提出了风洞内位置姿态视觉测量方法，可以在风洞内测量运动物体的六维的位置姿态信息，但是也没有考虑观察视窗等导致的光线偏折现象。
发明内容
本发明要解决的技术难题是克服现有技术的缺陷，发明一种基于多折射模型的风洞投放物高速位姿测量方法，解决由于被测物体速度高、自身滚转等问题所产生的测量难题。测量方法采用基于双目视觉对高速运动的滚转体位姿进行测量，在风洞环境下采用多折射模型代替传统的小孔成像模型，将传统的针对非滚转体位姿测量方法扩展为对于滚转体位姿测量方法，并提高了测量的精确度；解决了位姿测量过程中由于被测物体滚转所产生的标记点消隐问题，并且考虑多折射问题，进行新的图像匹配方法，解决了风洞环境下视觉测量系统特征匹配正确率低、测量精度较差等问题。
本发明所采用的技术方案是一种基于多折射模型的风洞投放物高速位姿测量方法，其特征是，测量方法首先建立风洞实验环境下的多折射模型；然后利用多折射模型对视觉测量图像数据进行处理优化，得到左、右图像特征信息的正确匹配；最终利用坐标转换方法求的投放物位姿测量信息；该方法的具体提取步骤如下：
步骤1建立风洞多折射模型
投放物特征标记所发出的光需要先后穿过风洞高速气流4、风洞观察窗3、镜头前非流动空气2才能进入左、右视觉测量设备1、5进行感光成像；高速气流4与风洞观察窗3以及镜头前分流动空气2将会引起光路的折射现象，导致图象畸变；设各种物体的折射率为μi，i＝0,1,2,其中，空气折射率为μ0、玻璃折射率为μ1、流动气体折射率为μ2，折射后光线传播方向组成的折射矢量vi，i＝0,1,2,时分别为：镜头前非流动空气2折射矢量为风洞观察窗3内折射矢量为高速气流4折射矢量为镜头前非流动气流2厚度为d0、风洞观察窗3厚度为d1、高速气流4厚度为d2；由斯涅尔定律可知：
μi sinθi＝μi+1sinθi+1，i＝0,1,2, (1)
其中，θi为各折射矢量vi与风洞视察窗(3)法矢量n的夹角；整理得：
vi+1＝ai+1vi+bi+1n (2)
其中，
根据多折射几何原理，得出P点的三维坐标为：
步骤2基于多折射模型的图像匹配方法
由于视窗所产生的多折射问题的时候，相机的左右视图的匹配点将不符合传统的极限约束，这样将会导致左右图像标记点特征匹配的准确性降低，影响最终测量精度。由于高速气流相对于左、右视觉测量设备的折射矢量v2 l，v2 r交于点P，所以风洞观察窗相对于左、右视觉测量设备的折射矢量v2 l，v2 r与P在同一个平面内，得到以下公式：
(q2 l-q2 r)·(v2 l×v2 r)＝0 (4)
其中，q2 l是左视觉测量设备的折射矢量v2 l与观察视窗的交点，q2 r是右视觉测量设备的折射矢量v2 r与观察视窗的交点，该方法首先根据光流方法找到合作标记的运动轨迹，然后取左图像运动轨迹上某些时刻的特征点与右图像的特征点利用公式(4)进行判定，满足公式(4)即表示左右图像标记点是匹配点，当左图像一条轨迹上多个标记点与右图像上的标记点分别匹配时，说明两条轨迹是匹配的，即两条轨迹上多有标记点是匹配的；
步骤3基于坐标转换的位姿求解方法
首先建立投放物局部坐标系OtXtYtZt为：世界坐标系中点pow＝(xp yp zp)T与局部坐标系下点的转换关系如下：
其中：Twb＝(x0 y0 z0)为投放物局部坐标系原点即投放物质心在世界坐标系中的坐标；Rwb为局部坐标系相对于世界坐标系的旋转矩阵；
其中，-θZ，-θX，-θY为投放物相对于世界坐标系的偏航角，俯仰角和滚转角，求得Rwb，Twb即单个投放物对于世界坐标系的位姿信息，以此求得所有投放物位姿信息。
本发明的有益效果对风洞实验环境中的多层光折射问题进行分析，利用光折射模型对测量过程进行优化从而提高风洞环境下投放物位姿测量精度。不仅矫正了图像由于光折射畸变，提高了测量的精度，而且提高了图像特征匹配的正确率。
附图说明
图1为基于多折射模型的风洞投放物高速位姿测量方法原理图。其中，1-视觉测量设备、2-镜头前非流动空气、3-风洞观察窗、4-高速气流。
图2为位姿测量方法流程图。
具体实施方式
以下结合技术方案和附图详细叙述本发明的具体实施方式。
附图1为基于多折射模型的风洞投放物高速位姿测量方法原理图。根据几何折射原理对实验环境进行分析，建立多折射模型，利用多折射模型对测量过程进行校对优化，最终求解位姿参数。该方法的流程如附图2所示，具体测量步骤如下：
步骤1建立风洞多折射模型
投放物特征标记所发出的光需要先后穿过风洞高速气流4、风洞观察窗3、镜头前非流动空气2才能进入视觉测量设备1进行感光成像；高速气流4与风洞观察窗3以及镜头前分流动空气2将会引起光路的折射现象，以镜头前非流动空气折射率为1、风洞观察窗折射率1.634以及高速气流1.315；各层的厚度为镜头前非流动气流100mm、风洞观察窗150mm、高速气流300mm；代入公式(1)、(2)、(3)求得风洞环境多折射模型。
步骤2基于多折射模型的图像匹配方法
考虑由于视窗所产生的多折射问题的时候，相机的左右视图的匹配点将不符合传统的极限约束，这样将会导致左右图像标记点特征匹配的准确性降低，影响最终测量精度。所以发明了基于多折射共面约束的图像特征匹配方法，由于高速气流相对于左、右视觉测量设备的折射矢量v2 l，v2 r交于点P，所以风洞观察窗相对于左、右视觉测量设备的折射矢量q2 l、q2 r与P在同一个平面内，取左图像运动轨迹上某些时刻的特征点与右图像的特征点利用公式(4)进行判定，满足公式(4)即表示左右图像标记点是匹配点，当左图像一条轨迹上多个标记点与右图像上的标记点分别匹配时，说明两条轨迹是匹配的，即两条轨迹上多有标记点是匹配的。
步骤3基于坐标转换的位姿求解方法，利用公式(5)、(6)求得投放物位姿参数为：
本发明对风洞环境高速气流以及观察窗引起的光折射现象进行分析，建立更为精确的多折射模型对位姿测量过程进行校正优化，并且提出新的图像匹配方法，以提高存在多折射现象的风洞投放物高速位姿测量精度与稳定性。
Claims (1)
1.一种基于多折射模型的风洞投放物高速位姿测量方法，其特征是，测量方法首先建立风洞实验环境下的多折射模型；然后利用多折射模型对视觉测量图像数据进行处理优化，得到左、右图像特征信息的正确匹配；最终利用坐标转换方法求得投放物位姿测量信息；该方法的具体步骤如下：
步骤1建立风洞多折射模型
高速气流(4)与风洞观察窗(3)以及镜头前非流动空气(2)会引起光路折射现象，设各种物体的折射率为μi，i＝0,1,2,其中，镜头前非流动空气(2)折射率为μ0、风洞观察窗(3)折射率为μ1、高速气流(4)折射率为μ2，折射后光线传播方向组成的各折射矢量vi，i＝0,1,2,其中，镜头前非流动空气(2)折射矢量为v0、风洞观察窗(3)内折射矢量为v1、高速气流(4)折射矢量为v2；镜头前非流动空气(2)厚度为d0、风洞观察窗(3)厚度为d1、高速气流(4)厚度为d2；由斯涅尔定律：
μisinθi＝μi+1sinθi+1，i＝0,1,2, (1)
其中，θi为各折射矢量vi与风洞视察窗(3)法矢量n的夹角；整理得：
vi+1＝ai+1vi+bi+1n (2)
其中，ai+1＝μi/μi+1,
根据多折射几何原理，得出P点的三维坐标为：
步骤2基于多折射模型的图像匹配方法
设高速气流(4)相对于左、右视觉测量设备(1、5)的折射矢量v2 l，v2 r交于点P，所以风洞观察窗(3)相对于左、右视觉测量设备(1、5)的折射矢量v2 l，v2 r与(P)在同一个平面内，得到以下公式：
(q2 l-q2 r)·(v2 l×v2 r)＝0 (4)
其中，q2 l是左视觉测量设备(1)的折射矢量v2 l与风洞观察窗(3)的交点，q2 r是右视觉测量设备(5)的折射矢量v2 r与风洞观察窗(3)的交点；
该方法首先根据光流方法找到合作标记的运动轨迹，然后取左图像运动轨迹上某些时刻的特征点与右图像的特征点利用公式(4)进行判定，满足公式(4)即表示左右图像标记点是匹配点，当左图像一条轨迹上多个标记点与右图像上的标记点分别匹配时，说明两条轨迹是匹配的，即两条轨迹上多个标记点是匹配的；
步骤3基于坐标转换的位姿求解方法
首先建立投放物局部坐标系OtXtYtZt为：世界坐标系中点pow＝(xp yp zp)T与局部坐标系下点的转换关系如下：
其中：Twb＝(x0 y0 z0)为投放物局部坐标系原点即投放物质心在世界坐标系中的坐标；Rwb为局部坐标系相对于世界坐标系的旋转矩阵；
其中,-θZ，-θX，-θY为投放物相对于世界坐标系的偏航角，俯仰角和滚转角，求得Rwb，Twb即单个投放物对于世界坐标系的位姿信息，以此求得所有投放物位姿信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610594960.4A CN106225668B (zh) | 2016-07-27 | 2016-07-27 | 基于多折射模型的风洞投放物高速位姿测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610594960.4A CN106225668B (zh) | 2016-07-27 | 2016-07-27 | 基于多折射模型的风洞投放物高速位姿测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106225668A CN106225668A (zh) | 2016-12-14 |
| CN106225668B true CN106225668B (zh) | 2018-11-09 |
Family
ID=57533888
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610594960.4A Active CN106225668B (zh) | 2016-07-27 | 2016-07-27 | 基于多折射模型的风洞投放物高速位姿测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106225668B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107560554A (zh) * | 2017-09-04 | 2018-01-09 | 大连理工大学 | 一种基于旋转透镜的三维信息视觉测量方法 |
| CN112001991B (zh) * | 2020-10-27 | 2021-01-26 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种高速风洞动态油流图谱图像处理方法 |
| CN112419410B (zh) * | 2020-11-20 | 2021-10-19 | 北京航空航天大学 | 一种基于水下斯涅尔窗口边缘辨识的水平姿态确定方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101699237A (zh) * | 2009-11-20 | 2010-04-28 | 中国航空工业空气动力研究院 | 用于风洞模型试验的三维模型姿态角视频测量系统 |
| CN103256896A (zh) * | 2013-04-19 | 2013-08-21 | 大连理工大学 | 一种高速滚转体位姿测量方法 |
| CN103364167A (zh) * | 2013-07-15 | 2013-10-23 | 中国航天空气动力技术研究院 | 一种观察窗折射偏移校正方法 |
| CN105571518A (zh) * | 2016-01-19 | 2016-05-11 | 大连理工大学 | 基于折射图像偏差的三维信息视觉测量方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11160021A (ja) * | 1997-11-27 | 1999-06-18 | Nippon Telegr & Teleph Corp <Ntt> | 広域３次元位置計測方法及び装置 |
| JP3491041B2 (ja) * | 2001-09-26 | 2004-01-26 | 独立行政法人航空宇宙技術研究所 | 磁力支持用風洞模型 |
-
2016
- 2016-07-27 CN CN201610594960.4A patent/CN106225668B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101699237A (zh) * | 2009-11-20 | 2010-04-28 | 中国航空工业空气动力研究院 | 用于风洞模型试验的三维模型姿态角视频测量系统 |
| CN103256896A (zh) * | 2013-04-19 | 2013-08-21 | 大连理工大学 | 一种高速滚转体位姿测量方法 |
| CN103364167A (zh) * | 2013-07-15 | 2013-10-23 | 中国航天空气动力技术研究院 | 一种观察窗折射偏移校正方法 |
| CN105571518A (zh) * | 2016-01-19 | 2016-05-11 | 大连理工大学 | 基于折射图像偏差的三维信息视觉测量方法 |
Non-Patent Citations (2)
| Title |
|---|
| A photogrammetry pose measurement method for moving targets in a wind tunnel;Xin Ma等;《Instrumentation and Measurement Technology Conference Proceedings,2016 IEEE International》;20160725;全文 * |
| 基于彩色编码的副油箱风洞模型位姿测量方法;刘巍等;《航空学报》;20150525;第36卷(第5期);第1556-1563页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106225668A (zh) | 2016-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103575227B (zh) | 一种基于数字散斑的视觉引伸计实现方法 | |
| CN104165750B (zh) | 立体视觉结合陀螺仪风洞模型位姿测量方法 | |
| CN105300316B (zh) | 基于灰度重心法的光条中心快速提取方法 | |
| CN109269466B (zh) | 基于特征点的靶面相对位姿测量方法及系统 | |
| CN105674880B (zh) | 基于双目原理的接触网几何参数测量方法及系统 | |
| CN106290256B (zh) | 基于视频测量的定量背景纹影方法 | |
| CN106225668B (zh) | 基于多折射模型的风洞投放物高速位姿测量方法 | |
| CN106596037B (zh) | 风洞试验模型绕流密度投影场的视频测量方法 | |
| CN103364171A (zh) | 一种高速风洞模型姿态视频测量系统及测量方法 | |
| CN110470226A (zh) | 一种基于无人机系统的桥梁结构位移测量方法 | |
| CN109766758A (zh) | 一种基于orb特征的视觉slam方法 | |
| CN104034269A (zh) | 一种单目视觉测量方法与装置 | |
| CN108692656B (zh) | 一种激光扫描数据获取方法及装置 | |
| CN111427451A (zh) | 采用扫描仪与眼动仪确定注视点在三维场景中位置的方法 | |
| Zhang et al. | Robust method for measuring the position and orientation of drogue based on stereo vision | |
| CN103247032B (zh) | 一种基于姿态补偿的微弱扩展目标定位方法 | |
| CN113724326B (zh) | 一种自主空中加油场景下锥套目标的单目视觉位姿解算方法 | |
| CN103954220A (zh) | 撞桥试验中船体运动状态数字图像测量方法 | |
| CN108132058A (zh) | 摄影测量系统动态位移测量误差校准装置及方法 | |
| CN109724586A (zh) | 一种融合深度图和点云的航天器相对位姿测量方法 | |
| CN112966571A (zh) | 一种基于机器视觉的立定跳远腾空高度量测方法 | |
| CN112712566B (zh) | 基于结构参数在线校正的双目立体视觉传感器测量方法 | |
| Du et al. | Autonomous measurement and semantic segmentation of non-cooperative targets with deep convolutional neural networks | |
| CN111912606B (zh) | 一种基于光学折射的相机物距标定方法 | |
| Liu et al. | Fusing dense features and pose consistency: A regression method for attitude measurement of aircraft landing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |