CN103364167A - 一种观察窗折射偏移校正方法 - Google Patents
一种观察窗折射偏移校正方法 Download PDFInfo
- Publication number
- CN103364167A CN103364167A CN2013103002964A CN201310300296A CN103364167A CN 103364167 A CN103364167 A CN 103364167A CN 2013103002964 A CN2013103002964 A CN 2013103002964A CN 201310300296 A CN201310300296 A CN 201310300296A CN 103364167 A CN103364167 A CN 103364167A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- point
- coordinate
- refraction
- target
- 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
Images
Landscapes
- Image Analysis (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
一种观察窗折射偏移校正方法,是对高速风洞中利用摄影测量***测量模型姿态的测量值校正。先标定出摄影测量***的摄像机坐标系与观察窗玻璃平面所在坐标系之间的空间位置关系,确定光线经玻璃折射后的折射偏移量,精确获得风洞内部模型上标识点的三维坐标;利用本发明的方法可有效地提高摄影测量***在高速风洞模型姿态测量应用中的测量精度。
Description
技术领域
本发明涉及一种用于摄影测量***中玻璃窗的折射偏移校正方法,特别是对高速风洞模型姿态测量中玻璃观察窗折射偏移校正方法。
背景技术
模型在风洞吹风过程中的姿态角是风洞试验的基础数据,测量结果的一个微小变化,常常会对实际应用价值造成巨大的影响。诸如在攻角测量过程中,0.01°的测量误差能够给阻力系数带来0.0001的不确定度,这将对飞行器实际巡航过程中的有效载荷带来1%的变化。在铰链力矩试验和测压试验中,无法用弹性角去修正名义攻角,支杆和天平的变形直接影响测量结果的准确程度。动态试验中,模型姿态角的变化和模型的振动对常规加速度计响应产生扰动,给攻角测量带来比较明显的误差。同时,风洞试验模拟要求的日益提高,以前将模型视为刚性体的假设也难以满足高精度的要求。
基于摄影测量技术的***作为一种“非接触式”的测量***,可以实现动态瞬时测量和空间立体测量,在低速风洞中,可以将定位***置于风洞内部,利用测量的标识点坐标直接计算模型的姿态。而在高速风洞中,由于风洞结构尺寸的限制,相机置于风洞内部通常难以取得满意的视场,且高速气流会导致定位基准的不断变化。因此,需要将摄像机放置于风洞外部,摄像机通过光学玻璃观察窗追踪模型上标识点的空间位置。光线通过观察窗玻璃时存在折射现象,而定位***的成像模型是以光线沿直线传播为前提的,如果直接利用定位***的成像模型计算标识点的坐标,测量结果根本无法满足高精度的风洞试验要求。
发明内容
本发明的目的在于:克服现有技术的不足,提供一种高速风洞模型姿态测量中观察窗折射偏移校正方法,消除光线折射引起的标识点坐标偏差,实现风洞内部模型上标准点三维坐标的准确测量。
本发明的技术解决方案:
一种观察窗折射偏移校正方法,用于高速风洞模型姿态测量时摄影测量***成像模型的校正,包括如下步骤:
(1)将标定靶固定在所述高速风洞的观察窗玻璃平面上,将摄影测量***的摄像机置于观察窗外部对所述标定靶成像,得到标定靶上标识点在摄影测量***的摄像机坐标系下的三维坐标;
(2)根据步骤(1)中得到的所述三维坐标计算出摄像机坐标系与所述高速风洞的观察窗玻璃平面所在坐标系之间的空间位置关系,得到坐标系间的转换矩阵;
(3)确定光线经玻璃折射后的折射偏移量;
(4)令摄像机直接对所述高速风洞内部的模型进行成像,得到该模型上标识点的三维坐标;
(5)将步骤(3)中得到的折射偏移量叠加到步骤(4)中得到的模型上标识点的三维坐标上,得到校正后的模型上标识点的三维坐标值。
所述步骤(2)中计算出摄像机坐标系与所述高速风洞的观察窗玻璃平面所在坐标系之间的空间位置关系,得到坐标系间的转换矩阵M=R+T,其中,R和T通过公式 计算得到,[X Y Z]T为标识点在标定靶坐标系下的坐标,[XC YC ZC]T为标识点在摄像机坐标系下的坐标。
所述步骤(3)中确定光线经玻璃折射后的折射偏移量具体为:
(3.1)取摄像机坐标系原点OC及入射光线与图像坐标系交点S,原点OC和交点S在摄像机坐标系下的坐标分别为(x,y,l)和(0,0,0),则交点S与原点OC在标定靶坐标系下的坐标为:
(3.3)计算出所述摄像机测量值相对真实值在标定靶坐标系下的偏移量如下:ΔX=-b·tgθDi,ΔY=d(tgθMi-tgθMr),ΔY=d(tgθMi-tgθMr),式中b为平行于Z轴且过摄像机直接获取位置点的直线和折射光线的交点到真实位置点沿Z方向的距离,ΔX=-b·tgθDi点到真实位置点沿Z方向的距离,点到真实位置点沿Z方向的距离;θDi为下部折射光线与Z轴的夹角;d为观察窗玻璃的厚度;θMi和θMr分别为XOY平面上入射光线的入射角和折射角;θUi为下部折射光线与Z轴的夹角;
(3.4)在摄像机坐标系下的偏移量ΔXc,ΔYc和ΔZc可表示为: 则标识点在摄像机坐标系下的真实坐标为(Xc+ΔXc,Yc+ΔYc,Zc+ΔZc)。
本发明与现有技术相比的优点为:
1、本发明一种观察窗折射偏移校正方法,将摄像机光线经过观察窗时的折射偏移叠加到成像模型,提高了模型姿态测量的精度;
2、本发明一种观察窗折射偏移校正方法,使得吹风过程中摄像机可以置于风洞外部测量模型姿态,避免了摄像机置于内部时存在测量基准抖动的问题,同时将摄影测量技术的应用拓展到高速风洞,尤其是小尺寸的高速风洞。
附图说明
图1为本发明校正方法流程图;
图2为标定靶形状、布点方式及坐标系设置示意图;
图3为图像坐标系、摄像机坐标系、标定靶坐标系与目标坐标系间相互位置关系示意图;
图4为折射偏移校正原理图;
图5、图6和图7分别为图4中局部区域I、II和III的放大图。
具体实施方式
图1为本发明一种观察窗折射偏移校正方法流程图,如图1所示,本发明一种观察窗折射偏移校正方法包括如下步骤:
(1)标定靶结构示意图如图2所示,建立标定靶坐标系O-XYZ,将标定靶固定在所述高速风洞的观察窗玻璃平面上,其中XOY面与玻璃平面重合,将摄影测量***的摄像机置于观察窗外部对所述标定靶成像,得到标定靶上标识点在摄影测量***的摄像机坐标系下的三维坐标;
(2)图2标定靶上各标识点在标定靶坐标系中的三维坐标已知,再根据步骤(1)中得到的所述三维坐标计算出摄像机坐标系与所述高速风洞的观察窗玻璃平面所在坐标系之间的空间位置关系,得到坐标系间的转换矩阵。
摄像机坐标系到标定靶坐标系的变换矩阵M通过一个正交变换矩阵R和一个平移变换矩阵T实现,即公式1
其中,[X Y Z]T为标识点在标定靶坐标系下的坐标,[XC YC ZC]T为标识点在摄像机坐标系下的坐标,R是正交旋转矩阵,T=[tx ty tz]T是摄像机坐标系原点在标定靶坐标系中的坐标。
(3)确定光线经玻璃折射后的折射偏移量
如图3所示,入射光线与图像坐标平面的交点S的坐标为(x,y),焦距OCO长度为l。取入射光线上点S与OC,其在摄像机坐标系下的坐标分别为(x,y,l)和(0,0,0),则入射光线上点S与OC在观察窗玻璃平面坐标系下的坐标为:
入射光线与标定靶坐标系XOY平面,即玻璃平面的夹角θi满足
标识点坐标偏移如图4所示,图中玻璃窗右侧上下两黑色粗实线为入射光线,左侧两黑色粗实线为定位***进行三角求交计算重心位置时所用光线,细虚线为入射光线经过玻璃窗折射后的实际光线。如局部放大图5所示定位***计算出的重心位置和标识点的实际位置分别为Pv和Pr,在标定靶坐标系下,计算结果相对实际位置沿Z轴的偏移量为沿X轴的偏移量为则
则有:
计算结果相对实际位置沿Y方向的偏移为定位***的中间摄像机投射光线经过玻璃窗时的偏移量,有:
ΔY=d(tgθMi-tgθMr) (6)
如图6和图7所示,当投射角、玻璃厚度及折射率已知时,单个光线折射角和折射偏移量利用公式(7)计算:
O=d(tgθi-tgθr)
由上述公式(1)~(7)可分别计算出三个镜头的折射角和折射偏移量θUr,OU、θMr,OM与θDr,OD。
(4)标识点坐标折射偏移校正
公式(5)和(6)计算出的ΔX,ΔY和ΔZ是标识点在玻璃平面坐标系下的坐标偏移量,最后需要将它们转换为标识点在摄像机坐标系下的偏移量ΔXc,ΔYc和ΔZc
令摄像机直接对所述高速风洞内部的模型进行成像,得到该模型上标识点在摄像机坐标系下的三维坐标;将公式(8)中得到的转换后的折射偏移量叠加到摄像机直接测量得到的标识点的三维坐标上,得到校正后的标识点三维坐标值为(Xc+ΔXc,Yc+ΔYc,Zc+ΔZc)。
Claims (4)
1.一种观察窗折射偏移校正方法,用于高速风洞模型姿态测量时摄影测量***成像模型的校正,其特征在于包括如下步骤:
(1)将标定靶固定在所述高速风洞的观察窗玻璃平面上,将摄影测量***的摄像机置于观察窗外部对所述标定靶成像,得到标定靶上标识点在摄影测量***的摄像机坐标系下的三维坐标;
(2)根据步骤(1)中得到的所述三维坐标计算出摄像机坐标系与所述高速风洞的观察窗玻璃平面所在坐标系之间的空间位置关系,得到坐标系间的转换矩阵;
(3)确定光线经玻璃折射后的折射偏移量;
(4)令摄像机直接对所述高速风洞内部的模型进行成像,得到该模型上标识点的三维坐标;
(5)将步骤(3)中得到的折射偏移量叠加到步骤(4)中得到的模型上标识点的三维坐标上,得到校正后的模型上标识点的三维坐标值。
2.根据权利要求1所述的一种观察窗折射偏移校正方法,其特征在于:所述步骤(2)中计算出摄像机坐标系与所述高速风洞的观察窗玻璃平面所在坐标系之间的空间位置关系,得到坐标系间的转换矩阵M=R+T,其中,R和T通过公式 计算得到,[X Y Z]T为标识点在标定靶坐标系下的坐标,[XC YC ZC]T为标识点在摄像机坐标系下的坐标。
3.根据权利要求1所述的一种观察窗折射偏移校正方法,其特征在于:所述步骤(3)中确定光线经玻璃折射后的折射偏移量具体为:
(3.1)取摄像机坐标系原点OC及入射光线与图像坐标系交点S,原点OC和交点S在摄像机坐标系下的坐标分别为(x,y,l)和(0,0,0),则交点S与原点OC在标定靶坐标系下的坐标为:
(3.2)通过公式计算出入射光线与观察窗玻璃平面的夹角θi;
(3.3)计算出所述摄像机测量值相对真实值在标定靶坐标系下的偏移量如下:ΔX=-b·tgθDi,ΔY=d(tgθMi-tgθMr),ΔZ=-OU·ctgθUi-b,式中b为平行于Z轴且过摄像机直接获取位置点的直线和折射光线的交点到真实位置点沿Z方向的距离;θDi为下部折射光线与Z轴的夹角;d为观察窗玻璃的厚度;θMi和θMr分别为XOY平面上入射光线的入射角和折射角;θUi为下部折射光线与Z轴的夹角;
(3.4)在摄像机坐标系下的偏移量ΔXc,ΔYc和ΔZc可表示为: 则标识点在摄像机坐标系下的真实坐标为(Xc+ΔXc,Yc+ΔYc,Zc+ΔZc)。
4.根据权利要求1所述的一种观察窗折射偏移校正方法,其特征在于:将标定靶固定在观察窗玻璃面上,其中标定靶坐标系的XOY面与玻璃平面重合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310300296.4A CN103364167B (zh) | 2013-07-15 | 2013-07-15 | 一种观察窗折射偏移校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310300296.4A CN103364167B (zh) | 2013-07-15 | 2013-07-15 | 一种观察窗折射偏移校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103364167A true CN103364167A (zh) | 2013-10-23 |
CN103364167B CN103364167B (zh) | 2015-09-09 |
Family
ID=49366023
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310300296.4A Active CN103364167B (zh) | 2013-07-15 | 2013-07-15 | 一种观察窗折射偏移校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103364167B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698100A (zh) * | 2013-12-27 | 2014-04-02 | 中国人民解放军国防科学技术大学 | 超声速风洞及其确定方法 |
CN106121823A (zh) * | 2016-06-22 | 2016-11-16 | 中国人民解放军63820部队吸气式高超声速技术研究中心 | 飞行器异型曲面内流道流场可视化玻璃观察窗及设计方法 |
CN106225668A (zh) * | 2016-07-27 | 2016-12-14 | 大连理工大学 | 基于多折射模型的风洞投放物高速位姿测量方法 |
CN107963563A (zh) * | 2017-11-17 | 2018-04-27 | 中国石油集团渤海钻探工程有限公司 | 一种使测井仪器车的绞车滚筒与井口对正的方法 |
CN108198219A (zh) * | 2017-11-21 | 2018-06-22 | 合肥工业大学 | 用于摄影测量的相机标定参数的误差补偿方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1566906A (zh) * | 2003-06-11 | 2005-01-19 | 北京航空航天大学 | 一种基于平面靶标的结构光视觉传感器标定方法 |
US20060145703A1 (en) * | 2004-12-20 | 2006-07-06 | Marcus Steinbichler | Automatic component testing |
CN101270982A (zh) * | 2007-03-20 | 2008-09-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于立体视觉的风洞模型外形监测方法 |
CN101334267A (zh) * | 2008-07-25 | 2008-12-31 | 西安交通大学 | 数字影像测头矢量坐标变换标定与误差修正方法及其装置 |
CN102288131A (zh) * | 2011-05-12 | 2011-12-21 | 上海大学 | 物体360°轮廓误差的自适应条纹测量装置和方法 |
-
2013
- 2013-07-15 CN CN201310300296.4A patent/CN103364167B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1566906A (zh) * | 2003-06-11 | 2005-01-19 | 北京航空航天大学 | 一种基于平面靶标的结构光视觉传感器标定方法 |
US20060145703A1 (en) * | 2004-12-20 | 2006-07-06 | Marcus Steinbichler | Automatic component testing |
CN101270982A (zh) * | 2007-03-20 | 2008-09-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于立体视觉的风洞模型外形监测方法 |
CN101334267A (zh) * | 2008-07-25 | 2008-12-31 | 西安交通大学 | 数字影像测头矢量坐标变换标定与误差修正方法及其装置 |
CN102288131A (zh) * | 2011-05-12 | 2011-12-21 | 上海大学 | 物体360°轮廓误差的自适应条纹测量装置和方法 |
Non-Patent Citations (2)
Title |
---|
张广军: "《机器视觉》", 30 June 2005, 科学出版社, article "摄像机标定", pages: 84-85 * |
陈元杰等: "基于双目视觉的水下定位***", 《机电工程》, vol. 28, no. 5, 31 May 2011 (2011-05-31), pages 568 - 569 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698100A (zh) * | 2013-12-27 | 2014-04-02 | 中国人民解放军国防科学技术大学 | 超声速风洞及其确定方法 |
CN103698100B (zh) * | 2013-12-27 | 2016-05-04 | 中国人民解放军国防科学技术大学 | 超声速风洞及其确定方法 |
CN106121823A (zh) * | 2016-06-22 | 2016-11-16 | 中国人民解放军63820部队吸气式高超声速技术研究中心 | 飞行器异型曲面内流道流场可视化玻璃观察窗及设计方法 |
CN106121823B (zh) * | 2016-06-22 | 2018-04-24 | 中国人民解放军63820部队吸气式高超声速技术研究中心 | 飞行器异型曲面内流道流场可视化玻璃观察窗及设计方法 |
CN106225668A (zh) * | 2016-07-27 | 2016-12-14 | 大连理工大学 | 基于多折射模型的风洞投放物高速位姿测量方法 |
CN106225668B (zh) * | 2016-07-27 | 2018-11-09 | 大连理工大学 | 基于多折射模型的风洞投放物高速位姿测量方法 |
CN107963563A (zh) * | 2017-11-17 | 2018-04-27 | 中国石油集团渤海钻探工程有限公司 | 一种使测井仪器车的绞车滚筒与井口对正的方法 |
CN108198219A (zh) * | 2017-11-21 | 2018-06-22 | 合肥工业大学 | 用于摄影测量的相机标定参数的误差补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103364167B (zh) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103364171B (zh) | 一种高速风洞模型姿态视频测量***及测量方法 | |
CN104406607B (zh) | 一种多视场复合光学敏感器的标定装置及方法 | |
CN101858755B (zh) | 一种星敏感器的标定方法 | |
CN111076880B (zh) | 一种考虑相机姿态变化的大跨桥梁多点挠度测量方法 | |
CN101334267B (zh) | 数字影像测头矢量坐标变换标定与误差修正方法及其装置 | |
CN101813465B (zh) | 非接触精密测转角的单目视觉测量方法 | |
CN102749047B (zh) | 一种轧辊轴承座尺寸测量方法 | |
CN103759669A (zh) | 一种大型零件的单目视觉测量方法 | |
CN103364167B (zh) | 一种观察窗折射偏移校正方法 | |
CN107256568A (zh) | 一种高精度机械臂手眼相机标定方法及标定*** | |
CN103559707B (zh) | 基于运动方靶标定物的工业定焦相机参数标定方法 | |
CN107356244B (zh) | 一种路侧单元天线的标定方法及装置 | |
CN101539397B (zh) | 物体三维姿态的精密光学测量方法 | |
CN102162738A (zh) | 摄像头与惯性传感器组合定位定姿***的标定方法 | |
CN107339935B (zh) | 用于全视角扫描测量***的靶标空间交会测量方法 | |
CN106290256A (zh) | 基于视频测量的定量背景纹影方法 | |
CN102721409B (zh) | 一种基于车身控制点的移动车辆三维运动轨迹的测定方法 | |
CN102136140B (zh) | 一种基于矩形图样的视频图像距离检测方法 | |
CN106483330A (zh) | 一种基于反光丝线姿态角视觉识别二维风速风向测试方法 | |
CN105486289A (zh) | 一种激光摄影测量***及相机标定方法 | |
CN110398208A (zh) | 基于摄影全站仪***的大数据变形监测方法 | |
CN101526338A (zh) | 线结构光视觉传感***结构参数的现场标定方法 | |
CN102012236A (zh) | 一种月球车双目视觉避障***标定方法 | |
CN103424124A (zh) | 基于图像测量技术的无磁惯导单元标定方法 | |
CN107478195A (zh) | 一种基于光学的空间物体姿态测量装置及其测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |