CN100342209C - 一种激光自准直测角***的标定方法 - Google Patents
一种激光自准直测角***的标定方法 Download PDFInfo
- Publication number
- CN100342209C CN100342209C CNB2005100694314A CN200510069431A CN100342209C CN 100342209 C CN100342209 C CN 100342209C CN B2005100694314 A CNB2005100694314 A CN B2005100694314A CN 200510069431 A CN200510069431 A CN 200510069431A CN 100342209 C CN100342209 C CN 100342209C
- Authority
- CN
- China
- Prior art keywords
- angle
- model
- beta
- collimation
- measuring system
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于测量技术,涉及对激光自准直测角***标定方法的改进。本发明基于激光自准直三自由度动态角度测量***,其步骤为:建立激光自准直测角***的简化标定模型,获取标定数据对,标定模型的参数,得到***的测量模型。本发明标定过程简单,标定精度高、标定数据量较少,***的测量时间较短,适合现场标定。
Description
技术领域
本发明属于测量技术,涉及对激光自准直测角***标定方法的改进。
背景技术
激光自准直测角技术自从60年代出现以来,已成为一门日益重要的小角度精密测试技术,在实现小角度的多维、非接触测量中具有独特的优点,基于该技术研制的测角***具有灵敏度高、结构简单、使用方便、可以实现二维角度的静态和动态测量等特点。任何一个激光自准直测角***在构建完成之后,都需要经过一个标定过程来确定***测量模型中未知参数的精确值,然后***才能正确行使其功能。激光自准直测角***模型参数的有效标定方法一直是一个重要的研究内容。目前主要有以下方法:
多项式法。
采用多项式法标定***的目的是通过获得大量标定数据对,采用多项式拟合的方法得到***的测量模型。尚鸿雁、张广军在文章“三自由度动态角度实时测量***的研究”(光电工程,已录用)中描述了这种标定方法。该方法主要是通过将被测物体安装在精密三轴旋转台上,被测物体表面贴上反射镜。以精密转台转动的角度为真值,获得大量的标定数据对,然后采用多项式拟和的方法标定出激光自准直测量数学模型中的未知参数。采用多项式法标定测量***,标定精度较高,但是标定数据量较大,并且在拟合之前需要对标定点数据进行处理或者在拟合之后对标定误差进行补偿,计算过程复杂。如果减少标定数据,则会大大的降低标定精度,难以获得令人满意的结果。
查表法。
杨勇等在文章“一种新型用于陀螺定向的小角度测量***的研究”(航空精密制造技术,2000,36(6):31-34)中描述了这种方法:***标定时,通过高精度的测量仪器获得多组的标定数据对,即被测量和***的输出值。将标定数据对输入到计算机中,制成表格,存储在ROM中。实际测量时,根据测量***的输出结果进行查表从而得到被测量的值。这种标定方法需要的标定点数据量很大,标定精度难以保证,并且在查表过程中耗费时间较长,不适用于实时的动态测量。
发明内容
本发明的目的是:针对现有的激光自准直测角***标定方法中存在的不足,提出一种精度高、标定过程简单的标定方法即参数优化法,进一步提高测量***的标定精度和减少测量时间,改善其工程化应用的便捷性。
本发明的技术方案是:一种激光自准直测角***的标定方法,基于由以下物理设备组成的激光自准直三自由度动态角度测量***:一个精度在角秒级的三自由度旋转台,在旋转台的旋转中心有一个平面反射镜,激光测角传感器即测角***的光学部分,光电接收器件,数据采集卡和计算机,其特征在于,
1、建立激光自准直测角***的简化标定模型,其步骤如下:
1.1、在三自由度旋转台和光学***焦平面上建立三自由度旋转台坐标系和图像坐标系,均为右手系;
1.2、设测量距离为L,光学***的焦距为f,平面反射镜的法线向量坐标为{nx,ny,nz},平面反射镜初始空间位置为(o,θo),光学***与三自由度旋转台之间的安装位置为:俯仰方向偏了θv角,偏航方向偏了v角,光学***中心在三自由度旋转台坐标系中的坐标为(ΔHx,ΔHy,L)。光电接收器件的安装位置为:俯仰方向倾斜了u角,偏航方向倾斜了v角,光敏面中心点在图像坐标系中的坐标为(xuo,yuo,zuo);
1.3、获得激光自准直测角***的标定模型,将平面反射镜空间旋转,θ后,根据基于透视投影变换建立激光自准直测量数学模型的方法,得到激光自准直测角***的标定模型;
1.4、简化标定模型,
①标定模型中假设ΔHx=ΔHy=0,
②把光电接收器件的位置旋转关系近似成一个旋转角度β,即像平面绕光轴转动的角度,
③设定未知参数J标定像点坐标的横纵比,
基于上述模型参数的简化,得到激光自准直测角***的简化标定模型为:
其中:
r1~r9九个系数的表达式为
2、获取标定数据对,将被测物体安装在三自由度旋转台上,被测物体表面贴上平面反射镜,被测物体分别做偏航/俯仰转动,角度转动范围≤1.0000°,以t为测量点间隔,0<t<0.1000°,获取P个角度,P>10,每个角度下,分别进行左偏航、右偏航、上俯仰和下俯仰转动,读取光电接收器件的两个输出电压值,共获取4×P个数据对;
3、标定模型的参数,
根据获得的标定数据对,采用Levenberg-Marquardt优化算法把上述简化标定模型中9个未知参数(θo,o,θv,v,β,xuo,yuo,zuo,J)标定出来;
4、***的测量模型如下:当光电接收器件探测到一个像点坐标(xIm,yIm)时,根据下面给出的***测量模型,就可以计算出平面反射镜空间转动的角度(,θ):
式中
本发明的优点是:与现有的激光自准直测角***的标定方法相比,该方法的标定过程简单,标定精度高、标定数据量较少,***的测量时间较短。本方法适合现场标定,通过该方法得到的测量模型适合***动态角度的实时测量。
附图说明
图1是本发明激光自准直三自由度动态角度测量***示意图。图中1被测物体;2平面反射镜;3准直目镜;4激光器;5光电接收器件;6半透半反镜;7三轴转台;8支架;9数据采集与显示部分。
具体实施方式
下面对本发明做进一步详细说明。参见图1,本发明激光自准直测角***的简化标定模型的建立过程如下:
在三自由度旋转台和光学***焦平面上建立转台坐标系和图像坐标系,均为右手系。
设测量距离为L,光学***的焦距为f,平面镜的法线向量坐标为{nx,ny,nz},平面镜初始空间位置为(o,θo),光学***与转台之间的安装位置为:俯仰方向偏了θv角,偏航方向偏了v角,光学***中心在转台坐标系中的坐标为(ΔHx,ΔHy,L)。光电接收器件的安装位置为:俯仰方向倾斜了u角,偏航方向倾斜了v角,光敏面中心点在图像坐标系中的坐标为(xuo,yuo,zuo)。
当平面镜空间旋转,θ后,根据专利“一种激光自准直二维角度测量的建模方法”(申请号:200410102586.9)中描述的基于透视投影变换建立激光自准直测量数学模型的方法,得到激光自准直测角***的标定模型。
直接采用上述建立的标定模型进行***标定,存在两个问题:①通过多组标定数据对,无法直接求解得到标定模型中未知参数的值。如果采用非线性优化法求解得到标定参数的值,经过反算,得到的标定精度特别差。②无法从标定模型中推导出***测量模型的显式表达式。
本发明通过化简标定模型达到获得***测量模型的目的。
标定模型的化简。
①标定模型中假设ΔHx=ΔHy=0。
②把光电接收器件的旋转位置关系近似成一个旋转角度β,即像平面绕光轴转动的角度。
③设定未知参数J标定像点坐标的横纵比。
基于上述模型参数的简化,得到激光自准直测角***的简化标定模型为:
其中:
获取标定数据对,将被测物体安装在精密三轴旋转台上,被测物体表面贴上反射镜,被测物体分别做偏航/俯仰转动,角度转动范围≤1.0000°,以t为测量点间隔,0<t<0.1000°,获取P个角度,P>10。每个角度下,分别进行左偏航、右偏航、上俯仰和下俯仰转动,读取PSD的两个输出电压值,共获取4×P个数据对。
模型参数的标定。
通过获得的标定数据对,采用Levenberg-Marquardt优化算法把上述简化标定模型中9个未知参数(θo,o,θv,v,β,xuo,yuo,zuo,J)标定出来。(陈宝林,最优化理论与算法,清华大学出版社,1989)
***的测量模型。
当光电接收器件探测到一个像点坐标(xIm,yIm)时,根据下面给出的***测量模型,就可以计算出平面镜空间转动的角度(,θ):
式中
实施例
图1中的激光自准直三自由度动态角度测量***主要由以下物理设备组成:①一个精度在角秒级的三自由度旋转台,并在旋转台的旋转中心放置一平面反射镜。②激光测角传感部分,即测角***的光学部分,测量范围±1°,焦距f=168.44mm。光学***到被测物体的距离为L=1.2m。③光电接收器件为日本滨松公司生产的二维枕型位置敏感探测器PSD,型号为S1880。④数据采集与显示部分:主要由数据采集卡、计算机等组成。
标定数据对的获取。
将被测物体安装在精密三轴旋转台上,被测物体表面贴上反射镜。被测物体分别做偏航/俯仰转动,角度转动范围≤0.8500°。以0.1500°测量点开始,并以测量点间隔0.0500°分别进行左偏航、右偏航、上俯仰和下俯仰,同时读取PSD的两个输出电压值,直到转动角度为0.8500°。共获取15个角度,60个数据对,如表1所示。
表1 60个数据对
序号 | 物体转角 | 左偏航(V) | 右偏航(V) | 上俯仰(V) | 下俯仰(V) |
1 | 0.1500 | -0.9570,-0.1758 | 0.7129,0.1611 | 0.0488,-0.8447 | -0.2881,0.8252 |
2 | 0.2000 | -1.2305,-0.2393 | 0.9814,0.2246 | 0.0977,-1.1133 | -0.3418,1.1035 |
3 | 0.2500 | -1.5039,-0.2930 | 1.2598,0.2832 | 0.1514,-1.3965 | -3.3955,1.3818 |
4 | 0.3000 | -1.7773,-0.3613 | 1.5381,0.3320 | 0.2100,-1.6797 | -0.4443,1.6357 |
5 | 0.3500 | -2.0605,-0.4053 | 1.8115,0.3906 | 0.2637,-1.9434 | -0.4883,1.8848 |
6 | 0.4000 | -2.3291,-0.4736 | 2.0850,0.4395 | 0.3223,-2.2266 | -0.5420,2.1338 |
7 | 0.4500 | -2.6074,-0.5420 | 2.3682,0.4932 | 0.3809,-2.4951 | -0.5908,2.3828 |
8 | 0.5000 | -2.8809,-0.5957 | 2.6416,0.5469 | 0.4443,-2.7637 | -0.6348,2.6416 |
9 | 0.5500 | -3.1396,-0.6641 | 2.9150,0.6055 | 0.4932,-3.0420 | -0.6934,2.9053 |
10 | 0.6000 | -3.4131,-0.7178 | 3.1885,0.6543 | 0.5518,-3.3057 | -0.7373,3.1494 |
11 | 0.6500 | -3.6865,-0.7666 | 3.4570,0.7129 | 0.6006,-3.5596 | -0.7910,3.4229 |
12 | 0.7000 | -3.9502,-0.8350 | 3.7256,0.7666 | 0.6592,-3.8184 | -0.8398,3.6816 |
13 | 0.7500 | -4.2090,-0.8789 | 3.9795,0.8154 | 0.7129,-4.0674 | -0.9033,3.9795 |
14 | 0.8000 | -4.4580,-0.9424 | 4.2383,0.8691 | 0.7568,-4.3213 | -0.9616,4.2676 |
15 | 0.8500 | -4.7266,-0.9912 | 4.5020,0.9326 | 0.8203,-4.5703 | -1.0205,4.5654 |
标定结果。
基于简化的标定模型即式1,采用Levenberg-Marquardt优化算法优化的结果如下:
o=0.1036°,θo=0.0648°,v=0.0994°,θv=0.1831°,xuo=0.188mm,yuo=0.660mm,zuo=4.465mm,β=-11.9523°,J=1.0164。
将上述获得的标定参数的值带入到***的测量模型即式2中,把标定数据对中像点坐标(xIm,yIm)作为测量模型的输入值,计算模型的输出值。利用下式:
计算出平面镜转动的俯仰角和偏航角的RMS误差为:俯仰角的RMS等于0.23,偏航角的RMS等于0.45,测量时间约3μs。
Claims (1)
1、一种激光自准直测角***的标定方法,基于由以下物理设备组成的激光自准直三自由度动态角度测量***:一个精度在角秒级的三自由度旋转台,在旋转台的旋转中心有一个平面反射镜,激光测角传感器即测角***的光学部分,光电接收器件,数据采集卡和计算机,其特征在于,
1.1、建立激光自准直测角***的简化标定模型,其步骤如下:
1.1.1、在三自由度旋转台和光学***焦平面上建立三自由度旋转台坐标系和图像坐标系,均为右手系;
1.1.2、设测量距离为L,光学***的焦距为f,平面反射镜的法线向量坐标为{nx,ny,nz},平面反射镜初始空间位置为(o,θo),光学***与三自由度旋转台之间的安装位置为:俯仰方向偏了θv角,偏航方向偏了v角,光学***中心在三自由度旋转台坐标系中的坐标为(ΔHx,ΔHy,L)。光电接收器件的安装位置为:俯仰方向倾斜了u角,偏航方向倾斜了v角,光敏面中心点在图像坐标系中的坐标为(xuo,yuo,zuo);
1.1.3、获得激光自准直测角***的标定模型,将平面反射镜空间旋转,θ后,根据基于透视投影变换建立激光自准直测量数学模型的方法,得到激光自准直测角***的标定模型;
1.1.4、简化标定模型,
①标定模型中假设ΔHx=ΔHy=0,
②把光电接收器件的位置旋转关系近似成一个旋转角度β,即像平面绕光轴转动的角度,
③设定未知参数J标定像点坐标的横纵比,
基于上述模型参数的简化,得到激光自准直测角***的简化标定模型为:
其中:
1.2、获取标定数据对,将被测物体安装在三自由度旋转台上,被测物体表面贴上平面反射镜,被测物体分别做偏航/俯仰转动,角度转动范围≤1.0000°,以t为测量点间隔,0<t<0.1000°,获取P个角度,P>10,每个角度下,分别进行左偏航、右偏航、上俯仰和下俯仰转动,读取光电接收器件的两个输出电压值,共获取4×P个数据对;
1.3、标定模型的参数,
根据获得的标定数据对,采用Levenberg-Marquardt优化算法把上述简化标定模型中9个未知参数(θo,o,θv,v,β,xuo,yuo,zuo,J)标定出来;
1.4、***的测量模型如下:当光电接收器件探测到一个像点坐标(xIm,yIm)时,根据下面给出的***测量模型,就可以计算出平面反射镜空间转动的角度(,θ):
式中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100694314A CN100342209C (zh) | 2005-05-10 | 2005-05-10 | 一种激光自准直测角***的标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100694314A CN100342209C (zh) | 2005-05-10 | 2005-05-10 | 一种激光自准直测角***的标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1862221A CN1862221A (zh) | 2006-11-15 |
CN100342209C true CN100342209C (zh) | 2007-10-10 |
Family
ID=37389700
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100694314A Expired - Fee Related CN100342209C (zh) | 2005-05-10 | 2005-05-10 | 一种激光自准直测角***的标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100342209C (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101169320B (zh) * | 2007-10-29 | 2010-10-13 | 北京交通大学 | 运动物体的瞬时姿态测量装置 |
CN100559114C (zh) * | 2008-08-07 | 2009-11-11 | 北京空间机电研究所 | 一种高精度测角方法 |
CN102141386B (zh) * | 2010-12-29 | 2012-11-28 | 哈尔滨工业大学 | 卫星光通信终端光轴与终端基准面间夹角的测量方法 |
DE102011089837A1 (de) * | 2011-12-23 | 2013-06-27 | Hilti Aktiengesellschaft | Optisches System |
CN104296655B (zh) * | 2014-09-26 | 2017-03-29 | 中国科学院光电研究院 | 一种激光跟踪仪像旋公式初始角的标定方法 |
CN104501715B (zh) * | 2014-12-05 | 2017-06-27 | 中国科学院西安光学精密机械研究所 | 一种激光对中仪接收***及方法 |
CN104647144B (zh) * | 2014-12-12 | 2017-05-17 | 西安航天计量测试研究所 | 基于坐标转换的光学测角装置及方法 |
CN105423917B (zh) * | 2015-12-01 | 2018-05-29 | 中国科学院西安光学精密机械研究所 | 位置敏感探测器定位误差的标定方法 |
CN106840043B (zh) * | 2016-12-30 | 2019-09-10 | 青岛海润特长新电器科技有限公司 | 一种基于到达时间的角度校准装置及校准方法 |
CN107728565B (zh) * | 2017-09-29 | 2020-09-04 | 中国神华能源股份有限公司 | 臂架俯仰角度计算方法、装置以及机器可读存储介质 |
CN111750773B (zh) * | 2019-03-29 | 2022-02-18 | 南京理工大学 | 测量不同光斑点在位置敏感探测器上响应的方法 |
CN111504344B (zh) * | 2020-05-15 | 2022-03-11 | 天津时空经纬测控技术有限公司 | 用于对非接触姿态测量设备进行标定的标定***及方法 |
CN111623775B (zh) * | 2020-05-15 | 2022-10-04 | 天津时空经纬测控技术有限公司 | 运载体姿态测量***、方法以及存储介质 |
RU2757856C1 (ru) * | 2021-03-01 | 2021-10-21 | Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкопорация "Росатом") | Устройство программно-аппаратной имитации лазерного датчика угловых скоростей |
CN115854835B (zh) * | 2023-02-15 | 2023-05-19 | 中国空气动力研究与发展中心空天技术研究所 | 一种用于无人机活动喷管角度标定的工装 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2105701U (zh) * | 1991-07-22 | 1992-05-27 | 北京建筑工程研究所 | 准直激光束角误差矫正器 |
CN1068418A (zh) * | 1991-07-09 | 1993-01-27 | 北京市建筑工程研究所 | 一种激光准直装置 |
JP2000299521A (ja) * | 1999-02-15 | 2000-10-24 | Photonetics Sa | 波長のフィルタリング用の単一方向の自動整合した後方反射光学系 |
WO2003025499A1 (en) * | 2001-09-17 | 2003-03-27 | Sharp Kabushiki Kaisha | Tilt sensing apparatus |
-
2005
- 2005-05-10 CN CNB2005100694314A patent/CN100342209C/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1068418A (zh) * | 1991-07-09 | 1993-01-27 | 北京市建筑工程研究所 | 一种激光准直装置 |
CN2105701U (zh) * | 1991-07-22 | 1992-05-27 | 北京建筑工程研究所 | 准直激光束角误差矫正器 |
JP2000299521A (ja) * | 1999-02-15 | 2000-10-24 | Photonetics Sa | 波長のフィルタリング用の単一方向の自動整合した後方反射光学系 |
WO2003025499A1 (en) * | 2001-09-17 | 2003-03-27 | Sharp Kabushiki Kaisha | Tilt sensing apparatus |
Also Published As
Publication number | Publication date |
---|---|
CN1862221A (zh) | 2006-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100342209C (zh) | 一种激光自准直测角***的标定方法 | |
CN100342210C (zh) | 一种激光自准直测角零位基准误差测量方法 | |
CN101334267B (zh) | 数字影像测头矢量坐标变换标定与误差修正方法及其装置 | |
CN101078627A (zh) | 一种基于陀螺全站仪-激光标靶的盾构机自动导向***的在线标定方法 | |
CN1308656C (zh) | 测量激光光束平行性的装置 | |
CN1687707A (zh) | 光纤陀螺惯测装置快速启动和精度保证的工程实现方法 | |
CN106468544B (zh) | 基于光电自准直仪的卫星高精度测角方法 | |
CN104677280B (zh) | 摆臂式轮廓仪旋转轴空间状态标定方法 | |
CN103292738A (zh) | 一种球面面形误差绝对检测方法 | |
CN101046386A (zh) | 太阳敏感器测量基准转换方法与装置 | |
CN102207380B (zh) | 一种高精度的横轴倾斜误差补偿方法 | |
CN110375694A (zh) | 基于便携关节式坐标测量机圆光栅测角误差的自校准方法 | |
CN109506555A (zh) | 两轴u型架轴承孔同轴度的测量装置和测量方法 | |
CN113790682A (zh) | 一种全量程角度误差检测方法、***和存储介质 | |
CN1793811A (zh) | 数码相机可量测化检测方法 | |
CN112683443B (zh) | 一种气浮型动态扭矩校准装置及校准方法 | |
CN111044077A (zh) | 星敏感器测量坐标系与基准立方镜坐标系间的标定方法 | |
CN102175267B (zh) | 一种光电经纬仪水平角高精度补偿方法 | |
CN207456381U (zh) | 提高激光跟踪仪测量精度的装置 | |
CN104647144B (zh) | 基于坐标转换的光学测角装置及方法 | |
CN102155956A (zh) | 一种高精度垂直角横轴倾斜误差补偿方法 | |
CN102707092A (zh) | 一种基于角速率台的单波束激光测速仪标定方法 | |
CN113899324B (zh) | 基于单轴激光陀螺测角仪的多轴转台垂直度误差检测方法 | |
CN115326107A (zh) | 一种扫描反射镜镜面法线误差数字化标校方法 | |
CN102221371A (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20071010 Termination date: 20110510 |