一种光学单站姿态处理二义性正确解确定方法

胡小丽, 唐明刚, 蔡文泽, 吴海英, 张思琪, 张三喜

胡小丽, 唐明刚, 蔡文泽, 吴海英, 张思琪, 张三喜. 一种光学单站姿态处理二义性正确解确定方法[J]. 应用光学, 2023, 44(1): 55-60. DOI: 10.5768/JAO202344.0101008
引用本文: 胡小丽, 唐明刚, 蔡文泽, 吴海英, 张思琪, 张三喜. 一种光学单站姿态处理二义性正确解确定方法[J]. 应用光学, 2023, 44(1): 55-60. DOI: 10.5768/JAO202344.0101008
HU Xiaoli, TANG Minggang, CAI Wenze, WU Haiying, ZHANG Siqi, ZHANG Sanxi. Determination method for ambiguous correct solution of optical single-station pose processing[J]. Journal of Applied Optics, 2023, 44(1): 55-60. DOI: 10.5768/JAO202344.0101008
Citation: HU Xiaoli, TANG Minggang, CAI Wenze, WU Haiying, ZHANG Siqi, ZHANG Sanxi. Determination method for ambiguous correct solution of optical single-station pose processing[J]. Journal of Applied Optics, 2023, 44(1): 55-60. DOI: 10.5768/JAO202344.0101008

一种光学单站姿态处理二义性正确解确定方法

基金项目: 一般科研项目(012016018200A22304)
详细信息
    作者简介:

    胡小丽(1979—),女,工程师,主要从事光学姿态测量及数据处理技术研究。E-mail:huxiaoli32@163.com

  • 中图分类号: TN29;TP394.1

Determination method for ambiguous correct solution of optical single-station pose processing

  • 摘要:

    因单站成像的信息欠缺以及成像的对称性,光学单站姿态处理的二义性是普遍性难点问题,针对光学单站姿态处理二义性解问题,以目标主体直线矢量与姿态角对应关系为突破点,以主体直线矢量的物像映射关系为单站姿态处理基础,从处理结果界定局部二义性错误解,然后根据单站姿态求解关键步骤修正二义性解,修正前后处理结果与常规交会处理结果比较效果显著,在文中所述场景下,误差量级至少缩小20%,该方法为光学单站姿态处理获取正确解提供了可靠理论支持。

    Abstract:

    Due to the lack of information and the symmetrical imaging of single-station imaging, the ambiguity of the optical single-station pose processing is a common and difficult problem. Aiming at the ambiguous solution of optical single-station pose processing, the corresponding relationship between the target subject straight-line vector and the pose angle was taken as the breakthrough point, and the object-image mapping relationship of the target subject straight-line vector was taken as the basis of single-station pose processing. Defined the partial ambiguous error solution from the single-station pose processing results, and then the ambiguous solution was modified according to the key steps of single-station pose solution. The effects of processing results before and after ambiguous correction were significant compared with that of the conventional intersection processing results. In the described scenario, the error magnitude is reduced by at least 20%. This method provides a reliable theoretical support for obtaining the correct solution of optical single-station pose processing.

  • 为保证在连续变焦过程中获得目标区域的图像清晰一致,使变焦系统在长短焦极限位置均具有较高的调制传递函数[1-4],对光学系统的装调提出较高要求。目前,常用的装调方法有两种:一种是以光学定心加工为基准,根据平行光管星点像的状态,通过调整光学元件的不同空间姿态完成装调;另一种是以调焦组的入射光线光轴为基准,依次对变倍组、补偿组及后固定组完成共轴调整[5]。由于上述过程中光学元件数量较多,无法对光学系统中各光学元件的空间姿态进行量化,且装调周期长,装调精度低,已不再满足现阶段的装调需求。如何在变焦系统装调过程中定位影响像差的光学元件空间姿态,以及确保光学元件具有较高的共轴精度,是一项非常有价值的研究课题。

    本文所研究的变焦光学系统,通过机械补偿法来降低像差的影响[6-8]。该系统由4个不同组元相互作用完成成像,包括前固定镜组、变倍镜组、补偿镜组和后固定镜组。在成像过程中前固定镜组完成对光线的汇聚,为后面的变倍镜组提供参考目标;变倍镜组做线性运动,完成系统焦距的变化;补偿镜组通过光杆直线导轨实现与变倍组的相对运动,用以补偿运动过程中的像面位移。后固定镜组则是将补偿镜组所成物像汇聚在探测器靶面上,如图1所示。

    光学性能指标要求该系统焦距为13.8 mm~280 mm,像元大小为3.45 μm×3.45 μm,相对孔径为1/2.4~1/6.23。装调完成后要求空间频率140 lp/mm时,长短焦中心视场传递函数值分别大于0.45和0.55,长短焦轴外0.7视场传递函数值分别大于0.25和0.35,变焦全程光学系统相对畸变优于3.2%。为达到上述要求,需要严格保证系统的光轴一致性。

    图  1  变焦光学系统
    Figure  1.  Zoom optical system

    首先需要对光学元件完成微应力装调,通常使用Zygo干涉仪将其面形控制在RMS$ \leqslant $$0.025\lambda $,减小面型对像散等像差的影响。

    光学系统中各光学元件的空间姿态是影响该系统成像质量的主要因素,采用光学定心加工技术[9],对微应力光学元件共轴精度进行有效控制,保证各光学元件间同轴精度在0.01 mm以内,上下配合端面平行差$ \leqslant $0.005 mm。为了降低变焦光学系统的球差,要求调焦组和后固定组中各光学元件空气间隔公差$ \leqslant $±0.01 mm,变倍组与补偿组中各光学元件空气间隔公差$ \leqslant $±0.02 mm,变倍组中销钉中心与胶合镜2的间隔公差$ \leqslant $±0.015 mm,补偿组中销钉中心与透镜7的间隔公差$ \leqslant $±0.015 mm。

    使用1 m焦距平行光管对变焦系统在长短焦位置时的3#鉴别率板及$\varphi $0.10 mm星点像进行目视检测,所得鉴别率板图像锐利度较高,星点像较圆整。为量化光学系统指标,使用传递函数仪进一步检测,结果如表1所示。从表1可以看出,系统长短焦焦距满足设计要求,但长短焦传递函数MTF低于技术要求,表明初装完成的光学系统内部存在一定的像差,需要做进一步分析与调整。

    表  1  光学系统检测结果
    Table  1.  Detection results of optical system
    测量位置焦距/mm MTF(中心视场) MTF(0.7视场)
    理论值实际值 理论值实际值 理论值实际值
    长焦 13.8±0.2 13.85 0.45 0.36 0.25 0.133
    短焦 280±0.2 279.93 0.55 0.42 0.35 0.256
    下载: 导出CSV 
    | 显示表格

    图1可知,该变焦系统由多个组元构成,装调过程不可避免会引入人为误差,使光学系统不再具有对称性,从而产生各类像差。每一光学元件对于像差的影响各不相同,在装调过程中需要控制的变量也相对较多,增加了判断像差敏感元件的难度,降低了装调效率。本节在像差检测的基础上,使用CAA技术,根据敏感元件的分布,完成对影响像差敏感元件的定位工作。

    CAA(computer aided alignment)技术是指光学系统的计算机辅助装调[10],是一种针对传统光学系统装调过程中出现的光学元件空间姿态不可控、无量化、随机因子影响大等缺点提出的一种新方法,通常应用在离轴三反光学系统的装调中。CAA技术的核心是对光学系统中光学元件空间姿态进行量化及控制,其主要方法首先是使用干涉仪对光学系统进行波像差检测,然后根据检测结果及光学系统中心偏测量结果,进行Zemax或CodeV软件仿真,从而得出各光学元件的空间姿态位置及其调整量。最后通过使用各种调整工装,使各光学元件的空间姿态趋近于理论分析结果,从而完成光学系统装调。

    由定义可知,波像差是实际波面与理想波面之间的光程差。当光学元件面形或空间姿态发生微量变化时,便会引入波像差,导致成像质量发生相应变化,进一步表现为MTF指标不达标。对于精度要求较高的光学系统,通过定心加工校正各元件同轴度后,在系统装调中由于人为因素的存在,仍会引入装调误差,即光学元件偏心与倾斜,均会对该光学系统彗差、像散、场曲和畸变产生相应的变化[11-12],但是不会对系统的球差产生影响。因此在精密光学镜头装调过程中必须保证光学系统的共轴精度,消除由光学元件偏心或倾斜所造成的像差。针对上文初装后检测结果,采用Zygo干涉仪对其光学系统波像差进行检测,原理如图2所示。

    对初装完成的变焦光学系统进行波像差检测可知,在短焦极限位置的像散为0.013$\lambda $,彗差为1.282$\lambda $,球差为0.005$\lambda $;在长焦极限位置的像散为0.014$\lambda $,彗差为1.882$\lambda $,球差为0.005$\lambda $。由此可见,彗差是该系统主要的像差,是造成变焦光学系统长短焦MTF不达标的主要因素。

    图  2  干涉仪检测镜头波像差示意图
    Figure  2.  Schematic diagram of wave aberration testing by interferometer

    为了进一步有效定位产生彗差的光学元件,对变焦光学系统在长焦极限位置进行中心偏检测,结果如表2所示。

    表  2  偏心检测结果
    Table  2.  Results of eccentricity detection
    光学元件偏心/μm倾斜/(″)
    胶合镜1 12.324 54.8
    透镜2 8.268 102.3
    胶合镜2 3.764 71.4
    透镜4 5.955 32.7
    透镜5 2.015 66.2
    胶合镜3 3.906 35.9
    透镜7 1.554 80.3
    透镜8 5.676 66.8
    透镜9 4.534 42.3
    透镜10 10.357 25.3
    下载: 导出CSV 
    | 显示表格

    图3为变焦光学系统中光学元件在不同情况下偏心及倾斜对彗差的影响变化。结合现有装配精度,使用CodeV软件仿真建模,得到光学元件敏感等级由高至低分布(1~10)及相应MTF变化量,结果如表3所示。由此可知,胶合镜1、透镜2、透镜5、透镜7以及透镜10的不同空间姿态是变焦光学系统产生彗差的主要因素,即为该系统产生像差的敏感光学元件,需要做针对性调整。

    根据逆灵敏度公差分析方法[13-14],将上述变焦光学系统中心偏检测结果带入CodeV软件中构建正态分布模型并进行蒙特卡洛模拟仿真[15],确定误差理想分布区间,从而精确求出补偿量,即当胶合镜1、透镜5和透镜7的单边偏心小于0.01 mm以及透镜10单边偏心小于0.008 mm,胶合镜1倾斜小于$45''$,透镜2倾斜小于$60''$、透镜5倾斜小于$50''$,透镜7倾斜小于$60''$时,可满足理论设计对彗差的要求,使光学系统MTF达到指标要求。

    图  3  偏心与倾斜彗差仿真结果
    Figure  3.  Simulation results of eccentricity and tilt coma
    表  3  光学元件敏感等级
    Table  3.  Optical element sensitivity grades
    敏感等级偏心>±0.01 mm倾斜>0.003 rad($60''$)
    光学元件MTF变换量光学元件MTF变换量
    1 透镜10 ±0.241 透镜5 ±0.312
    2 透镜5 ±0.226 透镜7 ±0.283
    3 胶合镜1 ±0.218 胶合镜1 ±0.254
    4 透镜7 ±0.202 透镜2 ±0.226
    5 透镜4 ±0.093 胶合镜2 ±0.098
    6 胶合镜2 ±0.043 胶合镜3 ±0.063
    7 透镜2 ±0.036 透镜4 ±0.074
    8 胶合镜3 ±0.031 透镜8 ±0.031
    9 透镜9 ±0.025 透镜10 ±0.028
    10 透镜8 ±0.022 透镜9 ±0.025
    下载: 导出CSV 
    | 显示表格

    根据中心偏测量仪原理及应用方法[16],设计了一种共轴调试装置,通过芯轴与主镜筒内孔紧配合,完成对变焦相机主镜筒机械内孔轴与直线光杆导轨的平行性进行精确测量,保证动组组元光轴的共轴精度,从而降低彗差对光学系统成像质量的影响,测量原理如图4所示。

    图  4  测量原理
    Figure  4.  Principle of measurement

    表4为调整后的变焦光学系统光学元件偏心及倾斜调试后的测量结果。各光学元件偏心及倾斜均满足上文中对变焦系统敏感光学元件的分析结果。表5为变焦光学系统在长短焦极限位置时的MTF测量值。通过对变焦光学系统中敏感光学元件空间姿态的调整,测得长焦极限位置时轴上MTF值为0.453,短焦轴上MTF达到0.552;长焦轴外0.7视场MTF为0.256,短焦轴外0.7视场MTF为0.358,结果满足设计要求。调整后的变焦光学系统长短焦极限位置的彗差分别减小到0.011$\lambda $和0.008$\lambda $

    表  4  偏心调试结果
    Table  4.  Eccentric adjustment results
    光学元件偏心/μm倾斜/(″)
    胶合镜13.0599.1
    透镜22.69040.4
    胶合镜23.608202
    透镜44.36217.1
    透镜55.12717.7
    胶合镜35.24423.4
    透镜74.5076.7
    透镜83.72549.2
    透镜95.4152.3
    透镜101.88520.1
    下载: 导出CSV 
    | 显示表格
    表  5  长短焦MTF实测值
    Table  5.  Measured value of long and short focal length MTF
    测量位置MTF(中心视场) MTF(0.7视场)
    理论值实际值 理论值实际值
    长焦 0.45 0.453 0.25 0.256
    短焦 0.55 0.552 0.35 0.358
    下载: 导出CSV 
    | 显示表格

    本文采用CodeV软件分析了光学元件不同空间姿态变化时的像差特点,在此基础上确定了影响光学系统成像质量的敏感光学元件,并针对性地完成了调整。除此之外,采用了一种同轴调试装置,对变焦相机主镜筒机械内孔轴与直线光杆导轨的平行性进行了精确测量,保证了动组组元光轴的同轴精度,实现了变焦光学系统各组元偏心的快速调试。最终检测结果表明,在变焦光学系统装调过程中采用该方法,可降低装调过程中彗差的影响,使传递函数MTF满足设计要求。

  • 图  1   单站直线矢量姿态映射图示

    Figure  1.   Diagram of pose mapping for single-station straight line vector

    图  2   弹类、飞机类目标相对像面投影二义性说明

    Figure  2.   Ambiguous description of projection to image plane of missile and aircraft targets

    图  3   单站姿态处理二义性正确解修正步骤

    Figure  3.   Modified steps of ambiguous correct solution of single-station pose processing

    图  4   飞机姿态数据处理-平飞

    Figure  4.   Data processing of aircraft in level flight attitude

    图  5   飞机平飞段姿态数据经二义性修正后处理结果对比

    Figure  5.   Comparison of processing results of aircraft attitude data in level flight after ambiguous modification

    图  6   飞机姿态数据处理

    Figure  6.   Data processing of aircraft attitude

    图  7   飞机姿态数据经二义性修正后处理结果对比

    Figure  7.   Comparison of processing results of aircraft attitude data after ambiguous modification

    表  1   二义性修正前后结果指标对比

    Table  1   Comparison of outcome indexes before and after ambiguous modification

    相对误差/(°)修正前修正后
    ΔΨ RMSE23.832.11
    Δφ RMSE8.861.27
    Δγ RMSE5.631.20
    下载: 导出CSV

    表  2   二义性修正前后结果指标对比

    Table  2   Comparison of outcome indexes before and after ambiguous modification

    相对误差/(°)修正前修正后
    ΔΨ RMSE 11.18 2.65
    Δφ RMSE 18.55 3.15
    Δγ RMSE 27.03 2.19
    下载: 导出CSV
  • [1] 苏国中. 基于光电经纬仪影像的飞机姿态测量方法研究 [D]. 武汉: 武汉大学, 2005: 1-184.

    SU Guozhong. Research on the method of aircraft pose surveying from image sequences obtained by photoelectric theodolite [D]. Wuhan: Wuhan University, 2005: 1-184

    [2] 蔡伟柯;基于模型的单相机位姿参数估计方法研究[D];长沙: 国防科学技术大学, 2006: 1-37

    CAI Weike, Study on pose estimation using single camera basedon object model[D]. Changsha: National University of Defense Technology, 2006: 1-37.

    [3] 王晓明, 乔彦峰, 宋立维, 等. 用单站光电雷达一体化系统测量目标三维姿态[J]. 光学精密工程, 2011, 19(8): 1779-1786.

    WANG Xiaoming, QIAO Yanfeng, SONG Liwei, et al. Measurement of three-dimensional attitudeby single opto-electrical and radar integration instrument[J]. Optics and Precision Engineering, 2011, 19(8): 1779-1786.

    [4] 赵汝进, 刘恩海, 张文明, 等. 一种基于弱透视成像模型的目标姿态测量方法[J]. 光子学报,2014,43(5):0512002. doi: 10.3788/gzxb20144305.0512002

    ZHAO Rujin, LIU Enhai, ZHANG Wenming, et al. Measurements of 3D pose of objects based on weak perspective model[J]. Acta Photonica Sinica,2014,43(5):0512002. doi: 10.3788/gzxb20144305.0512002

    [5] 王鲲鹏, 张小虎, 于起峰. 基于目标特征点比例的单站图像定姿方法[J].应用光学, 2009, 30(6): 974-978.

    WANG Kunpeng, ZHANG Xiaohu, YU Qifeng, Measuring cooperative targets'pose based on proportions of feature points [J]. Journal of Applied Optics, 2009, 30(6): 974-978.

    [6] 唐自力, 马彩文, 刘波, 等. 单站光测图像确定空间目标三维姿态[J]. 光子学报,2004,33(12):1480-1485.

    TANG Zili, MA Caiwen, LIU Bo, et al. An approach of determination of objects 3D pose from mono-view[J]. Acta Photonica Sinica,2004,33(12):1480-1485.

    [7] 张世杰, 曹喜滨, 张凡, 等. 基于特征点的空间目标三维位姿单目视觉确定算法[J], 中国科学(信息科学), 2010, 40(4): 591-604.

    ZHANG Shijie, CAO Xibin, ZHANG Fan, et al. 3D position and pose determination method of spatial targets based on feature points from mono- view[J]. Scientia Sinica Informations, 2010, 40(4): 591-604.

    [8] 张三喜, 姚敏, 孙卫平. 高速摄像及其应用技术[M]. 北京: 国防工业出版社, 2006: 27–34.

    ZHANG Sanxi, YAO Min, SUN Weiping. High speed photogrammetry and application[M]. Beijing: National Defence Industry Press, 2006: 27–34.

    [9] 张三喜, 张伟光, 李宏凯, 等. 弹道特征参数摄像测量[M]. 北京: 国防工业出版社, 2014: 1–15.

    ZHANG Sanxi, ZHANG Weiguang, LI Hongkai, et al. Ballistic characteristic parameters of visionmeasurement[M]. Beijing: National Defence Industry Press, 2014: 1–15.

    [10] 赵学颜, 李迎春. 靶场光学测量[M]. 北京: 国防大学出版社, 2001: 168–230.

    ZHAO Xueyan, LI Yingchun. Optical measurement in shooting range[M]. Beijing: National Defense Press, 2001: 168–230.

    [11] 王之卓. 摄影测量原理[M]. 武汉: 武汉大学出版社, 2007.

    WANG Zhizhuo. The principles of photogrammetry[M]. Wuhan: Wuhan University Press, 2007.

    [12] 张祖勋, 张剑清. 数字摄影测量学[M]. 武汉: 武汉大学出版社, 1997.

    ZHANG Zuxun, ZHANG Jianqing. Digital photogrammetry[M]. Wuhan: Wuhan University Press, 2002.

    [13] LU Chenping , MEMBER I, GREGORY D H, et al. Fast and globally convergent pose estimation from video images[J]. IEEE Transactions ON Patten Analysisand Machine Intelligence, 2000, 22(6): 610-622.

    LU Chenping, MEMBER I, GREGORY D H, et al. Fast and globally convergent pose estimation from video images[J]. IEEE Transactions on Patten Analysisand Machine Intelligence,2000,22(6):610-622.

    [14] 张政, 张小虎, 傅丹. 一种高精度鲁棒的基于直线对应的位姿估计迭代算法[J]. 计算机应用,2008,28(2):326-329. doi: 10.3724/SP.J.1087.2008.00326

    ZHANG Zheng, ZHANG Xiaohu, FU Dan. Accurate and robust iterative pose estimation from line correspondences[J]. Journal of Computer Applications,2008,28(2):326-329. doi: 10.3724/SP.J.1087.2008.00326

    [15] 许允喜, 蒋云良, 陈方. 基于点和直线段对应的扩展正交迭代位姿估计算法[J]. 光学学报,2009,29(11):3129-3135. doi: 10.3788/AOS20092911.3129

    XU Yunxi, JIANG Yunliang, CHEN Fang. An extended orthogonal iteration algorithm for pose estimationfrom point and line segment correspondences[J]. Acta Optica Sinica,2009,29(11):3129-3135. doi: 10.3788/AOS20092911.3129

  • 期刊类型引用(5)

    1. 韦晓孝,戴振秋,何晓冰,周辰,胡蝶. 用于病理切片检测的大视场变焦显微镜光学系统设计. 光学技术. 2025(01): 7-12 . 百度学术
    2. 赵志明,王乃余,李治国. 基于霍夫变换的透镜中心偏差检测方法. 光学技术. 2024(04): 441-447 . 百度学术
    3. 韩颖,项建胜. 基于OptiCentric中心偏测量仪的红外光学镜组装调方法. 应用光学. 2022(04): 726-731 . 本站查看
    4. 徐彦,朱元,单良. 应用OPD-Scan像差仪对白内障超声乳化患者像差分析的意义. 中国医疗器械信息. 2021(21): 120-122 . 百度学术
    5. 赵希婷,张超,冀翼,刘辉,焦文春,黄阳,李重阳,张志飞. 超宽视场离轴光学系统畸变一致性校正技术. 应用光学. 2020(05): 1032-1036 . 本站查看

    其他类型引用(1)

图(7)  /  表(2)
计量
  • 文章访问数:  244
  • HTML全文浏览量:  135
  • PDF下载量:  30
  • 被引次数: 6
出版历程
  • 收稿日期:  2022-03-03
  • 修回日期:  2022-04-27
  • 网络出版日期:  2022-09-16
  • 刊出日期:  2023-01-14

目录

/

返回文章
返回