激光与光电子学进展, 2019, 56 (5): 051203, 网络出版: 2019-07-31   

基于非合作目标图像处理技术的大间距轴线一致性检测方法 下载: 1119次

Consistency Detection Method for Large Spacing Axis Based on Non-Cooperative Target Image Processing Technology
张勇 1,*吴浩 2牛刚 1
作者单位
1 中国人民解放军32181部队, 河北 石家庄 050000
2 陆军工程大学石家庄校区电子与光学工程系, 河北 石家庄 050003
摘要
在分析国内外大间距轴线一致性检测方法的基础上,提出一种基于非合作目标图像处理技术的轴线一致性检测方法。该方法任意选择远场中具有典型特征的景物作为非合作目标,利用不同光学传感器获取非合作目标图像,比较非合作目标在图像中的空间位置差异,进而得到轴线一致性检测结果。实验结果表明,该方法可满足大间距轴线一致性检测需求。将CCD成像系统夹持在被测对象机械轴线上,可实现机械轴线与光学传感器轴线间的一致性检测。该方法不受光学传感器工作波段限制,避免了其他检测方法体积重量大、对使用环境要求高等缺点,因此应用前景广阔。
Abstract
On the basis of analyzing the current methods for large spacing axes detection at home and abroad, a method of axial consistency detection based on non-cooperative target image processing technology is proposed, which selects scenes with typical features in far field as non-cooperative targets arbitrarily, utilizes different optical sensors to acquire non-cooperative target images, and obtains axial consistency detection results by comparing the spatial position differences of non-cooperative targets in images. Experimental results show that this method can meet requirement of large spacing axis consistency detection. While the CCD imaging system is reliably clamped on the mechanical axis of the target to be detected, the consistency detection between the mechanical axis and the optical sensor axis can also be achieved. This method is not limited by the optical sensor's working band and avoids the disadvantages of large volume, heavy carrying and high requirements on environment for other detection methods, so it has broad application prospects.

1 引言

基于火炮、雷达、导弹以及无人机等某一特定平台,由多个信息互联的任务单元相互协作,完成瞄准、侦察、跟踪、制导和标定等作战任务的平台光电装备,在各型**装备上得到了广泛应用。平台光电装备含可见光、微光、激光、红外和电视等光学传感器,而多光学传感器之间的轴线一致性是确保整个**系统作战效能的重要保证。**装备在安装、修理过程中需要进行轴线检测和调试,长时间野外环境影响也会造成轴线失调,需要适时进行校验。大跨度、宽光谱、多轴线是平台光电装备轴线一致性检测的主要特点。

目前,常用轴线一致性检测方法有投影靶板法[1-2]、平行光管法[3-6]、相交校靶法、惯性测量法[7]以及摄影测量方法[8]等。投影靶板法易受气候和场地条件限制,自动化程度低;平行光管法多用于实验室,要求被测对象中各光学传感器的物镜全部或部分包含在测试光路中,故需大通光口径,为适应宽波段且有效减轻重量,常采用折反式光学系统,并联合使用斜方棱镜、五角棱镜或平面反射镜组扩展平行光管口径;相交校靶法是一种解决大间距轴线一致性检测的有效手段,典型代表有以色列CI公司的先进**轴线检测系统O-AWBS[9-10]、德国Carl Zeiss公司的**装备视频校轴调整系统WASVB[11]和瑞典SCHILL公司的系列舰船轴线检测系统aligner-308等[12];惯性测量法基于陀螺稳定与光路自准直技术,无需架设平台或三脚架调校被测对象,解决了大口径平行光管无法完成的大间距机械轴线与光电轴线间一致性测量,美国AAI公司研制的先进校轴设备ABE-301A检测精度可达20″[13-14];摄影测量法具有非接触、测量范围大、测量快速和精度较高等特点[15-16],该方法应用于轴线检测的关键是大视场范围内基于多靶板的位姿跟踪测量技术,如挪威METRONOR公司研制的HarmoLign**校准系统测量精度优于20″,已推广应用于20个国家的空军[17-19]

综上所述,大间距轴线一致性检测方法在测量范围、精度、便携性及价格方面各有优缺点,相比之下,相交校靶法结构简单、操作便携,但需较长的通视距离,若能在保证测量精度前提下,缩短通视距离,将会大幅提高设备野外在线检测的适用性。为此,本文提出一种基于非合作目标图像处理技术的轴线一致性检测方法,通过选择远场中具有典型特征的景物作为非合作目标,采集不同光学传感器获取非合作目标图像,比较非合作目标在图像中的空间位置差异,进而得到轴线一致性检测结果。

2 有限远图像测量法

有限远图像测量方法利用被测光学传感器采集远场中某一固定景物作为非合作目标,通过计算不同光学传感器对图像中目标特征点的成像位置差异进行轴线一致性检测。图1为非合作目标位于有限远D处时的示意图。

图 1. 有限远图像测量法示意图(D为有限远)。(a)整体;(b)截面

Fig. 1. Diagram of finite distance image measurement method (D is finite). (a) Overall schematic; (b) section

下载图片 查看所有图片

被测对象两根光轴分别为O1O2,轴间距为L,在轴线方向相距W,与目标距离D。通过其中一个光学传感器瞄准远场中明显特征目标T(如建筑物屋顶、窗台棱角等),在与光轴垂直的图像平面P1P2内,O1O2的分划坐标分别为(x0,y0)和(x1,y1)。理想状态下,若O1O2平行且T位于无穷远时,TP1P2中的图像均应位于分划中心;若O1O2平行而T位于有限远时,TP1P2中的图像存在位置差异,该差异值应和投影靶板法在规定的测量距离下被测对象镜像轴线间距相同。如图2所示,匹配计算P1P2,求得O1O2分划中心在P2中的坐标分别为O1(x0,y0)和O2(x1,y2)。结合轴间距L,在某一特定观测距离(D)上,若O1O2不平行,则O1分划中心位于O'2(xa,ya)处,O1O'2的欧氏距离为两根光轴偏差,即

a=Rx0-xa2+y0-ya2,(1)

式中R为图像像素角分辨率。

图 2. 图像测量法原理图

Fig. 2. Schematic of image measurement

下载图片 查看所有图片

2.1 图像匹配算法

实际应用中,两根光轴的光学系统参数并不相同,所成像的放大倍率也不同,若要对两幅图像进行匹配计算,图像处理算法应具备仿射不变性。采用尺度不变特征变换(SIFT)算法匹配图像P1P2。首先,建立图像的尺度空间表示;然后,在尺度空间中搜索图像极值点,建立特征描述向量;最后,寻找两幅图像间的特征点对应关系。以P2为基准图像,从P1中选取含有光轴O1十字分划的图像块(包含尽可能多的特征或细节信息)作为待匹配模板,与P2进行模板匹配计算,选用标准平方差作为相似性度量函数,得到

σ=x,yI(x,y)-I'(x,y)2/x,yI(x,y)2·I'(x,y)2,(2)

式中I(x,y)和I'(x,y)分别为选取的图像块和待匹配图像块灰度值。图3给出了计算结果,当两个图形块达到最大相似度时,认为两个图像块达到了最佳匹配关系。

P2中求得最佳匹配图像块后,利用随机抽样一致性算法(RANSAC)求得变换矩阵,点对点配准两幅图像获得图像点的对应关系,所得结果如图4所示,A'点为O1分划中心在图像P2上的坐标,B点为O2分划中心坐标。

图 3. SIFT计算结果

Fig. 3. Calculation result by SIFT

下载图片 查看所有图片

图 4. 光学传感器采集图像及RANSAC配准后十字分划中心位置示意图。(a) O1采集图像;(b) O2采集图像;(c) O2采集图像放大图

Fig. 4. Images acquired by optical sensors and center position of cross-section after registration of RANSAC. (a) Image acquired by O1; (b) image acquired by O2; (c) enlarged image by O2

下载图片 查看所有图片

2.2 目镜图像采集装置

光学传感器成像有目镜观察和视频输出两种方式。人眼通过目镜观察光学系统不能存储成像结果,无法实现事后回放、图像处理和分析等功能,因此,设计目镜图像采集装置[20],利用磁性表座将其固定在光学系统目镜出瞳位置代替人眼,CCD相机将视频图像无线传输至数据处理系统处理并显示,利用接管组件及万向节组件调整CCD相机位置,确保CCD相机对准目镜出瞳。图5给出了目镜图像采集装置外观及示意图。

图 5. 目镜图像采集装置外观及示意图。(a)装置外观;(b)装置结构

Fig. 5. Appearance and diagram of eyepiece image acquiring device. (a) Appearance of the device; (b) device structure

下载图片 查看所有图片

2.3 机械轴线可视化

实际应用中,不但要求平台光电装备内的各光学传感器共轴,有时还需调校火炮身管与平台光电装备之间的镜炮一致性。为测量机械轴线和多光学传感器轴线一致性,还需将机械轴线可视化。因此,设计带有准直光学系统的机械连接插轴,将机械轴线转换为光轴[21],保证准直系统分划中心出射光线与火炮身管轴线平行。为实现与机械轴线内壁的紧密配合,可根据机械轴线形式和尺寸设计专用的膨胀式插轴,利用插轴膨胀顶起实现刚性连接,机械连接轴外观和内部结构如图6所示。其中顶块和弹簧通过连接轴固定在连接板上,顶块与连接轴压下后在弹簧弹力作用下自行弹起,连接板与插杆采用螺纹连接固定,插杆通过连接筒以螺纹连接方式与CCD相机连接固定。使用时,将插杆插入机械轴线内腔,顶块受压后在弹簧弹力作用下自动顶紧内壁,从而实现紧定功能;拔出插杆时,顶块在弹簧弹力作用下自动弹起。

图 6. 机械连接插轴外观及内部结构。(a)外观;(b)装置结构

Fig. 6. Appearance and internal structure of mechanical connection plug. (a) Appearance; (b) device structure

下载图片 查看所有图片

配合误差是影响轴线检测的主要误差源。设备出厂前必须进行自校准,即将机械连接插轴插入标准通孔中,前面放置平行光管,使其十字分划成像于CCD相机十字分划中心,然后顺时针转动360°,若平行光管十字分划像相对机械连接插轴CCD相机十字分划中心摆动量小于阈值,即认为符合出厂精度要求。

3 实验结果与分析

将激光测距机和望远镜安装在光学平板上,望远镜可做高低、方位及俯仰调整,利用大口径平行光管调整两台设备共轴。然后将实验装置对准远场中的典型目标,利用激光测距机测量目标和实验装置的距离,利用目镜图像采集装置分别采集激光测距机光学瞄准系统目镜出瞳图像和望远镜右目镜出瞳图像。激光测距机和望远镜间距L=300 mm。图7为实验装置图。在图像坐标系下进行计算,得到不同距离下的轴线偏差计算结果如表1所示。

图 7. 实验装置

Fig. 7. Experimental equipment

下载图片 查看所有图片

表1中坐标A'为O1分划中心在O2图像中所对应的坐标位置,坐标B为光轴O2图像分划中心位置,参考坐标为光轴O2分划密位读数为2的位置。分别计算LA'BLBC的像素长度,由图8可知, lBC对应10 mil,则在某一距离上,lA'B的空间长度即为像高H,H可表示为

H=D·lA'B/lBC·100,(3)

式中D为实验装置与目标的测距值。由于两根光轴平行,则测量误差角为

aarctanH-L/D,(4)

表 1. 两根光轴平行时轴线偏差计算结果

Table 1. Results of axis deviation when two optical axes are parallel

Coordinate A'Coordinate BCoordinate CW /mLA'B /pixellA'B /mmMeasurement /(″)
2247,15602249,15402277,201390020.10381.7818.74
2131,14642133,14462159,187190018.11382.8118.98
2135,14742136,14582159,183390016.03384.0319.26
1727,15621766,15351799,217639047.43288.226.23
1913,16311938,16121958,198439031.40328.7315.19
2151,15932166,15732186,193939025.00265.9917.98
2336,13852424,13472466,208022095.85287.2211.98
2299,16632369,16192400,220322082.68311.0310.34
2360,16532415,16142445,213122067.42286.4312.72
1925,15041984,14812002,177414563.32312.7918.20
1928,15261965,15121977,169714539.56309.4213.39

查看所有表

式中L为实验装置两根光轴的间距。由表1结果可知,测量误差角均值为14.82″,标准偏差为4.27″。

为进一步验证本文方法的测量精度,调整望远镜与激光测距机不共轴,具体方法如图8所示,将激光测距机十字分划中心对准远处某一典型目标,望远镜右目镜十字分划中心与目标错开一定距离,采集两幅图像并计算偏差量,利用望远镜分划密位线计算所得两根光轴的不平行偏差为1219″。实验装置至目标距离为2930 m,对同一目标计算5次并求均值。两根光轴不平行时,轴线偏差计算结果如表2所示,与理论值(1219″)之间的偏差均值为-15.37″,标准偏差为26.84″。对比表1表2结果可知,本文方法测量精度满足绝大部分被测对象检测需求。

图 8. 两根光轴不平行时的采集图像。(a) O1采集图像;(b) O2采集图像;(c) O1与O2分划中心偏差示意图

Fig. 8. Captured images when two optical axes are not parallel. (a) Images captured by O1; (b) images captured by O2; (c) center deviation of O1 and O2

下载图片 查看所有图片

表 2. 两根光轴不平行时轴线偏差

Table 2. Axis deviation when two optical axes are not parallel

Coordinate A'Coordinate BCoordinate CW /mLA'B /pixellA'B /mmMeasurement /(″)Deviation /(″)
1680,13181900,13151915,16779002205165.41183.8-35.2
1913,14322024,14292031,16139001115127.41175.1-43.9
1991,12582269,12492288,17049002785197.01191.1-27.9
2008,14202131,14162138,16003901232306.61219.90.9
1967,12652269,12382289,16913903032307.81220.61.6
1902,11022365,10532394,17443904652325.51229.910.9
1899,14252030,14192038,16042201311258.01179.5-39.5
1774,13232029,13092043,16702202551255.21176.8-42.2
1464,10931955,10691983,17622204921259.31180.7-38.3
1878,14312024,14192032,1602145146859.61222.83.8
1713,13301993,12972007,1661145282822.31169.7-49.3
1338,11261895,10571918,1749145561875.41245.326.3

查看所有表

4 误差分析

影响测量精度的误差主要包括机械轴线装配误差、图像匹配误差、目标测距误差以及被测对象各轴间距误差等。机械连接插轴装配误差由加工精度和出厂前自校准精度确定,一般情况下可控制σ1<30;图像匹配误差主要取决于CCD相机获取的图像显示精度,CCD镜头光学系统焦距f=80 mm,像元宽度w=10 μm,则CCD图像显示精度为

σ2=arctan(w/f)=25.78,(5)

观测目标距离D由激光测距机测距获得,误差为σ3<±1 m;轴间距L由机械加工图纸确定,误差为σ4<±10 mm。除σ1σ2外,轴间距误差σ3和测距误差σ4是影响精度的重要因素。将σ3σ4分别代入(4)式中进行蒙特卡罗误差分析,样本数取104,目标测距误差为正态分布,轴间距测量误差为均匀分布,所得结果如图9所示。

图 9. 轴间距测量误差为10 mm,测距误差为1 m时的误差分析结果

Fig. 9. Error analysis result for axis distance measurement error of 10 mm and range error of 1 m

下载图片 查看所有图片

图9可知,当轴间距测量误差为10 mm,测距误差为1 m,轴间距为100~1000 mm,实验装置与远场目标距离100 m时,系统测量误差范围从20″增大至29″;实验装置与远场目标距离200 m时,系统测量误差范围从10″增大至11″;而在50 m时,系统测量误差范围从42″增大至90″,可见为保证测量精度,当轴间距增大时,要求系统具有更大的测量距离。另外,轴间距越小,两个光学传感器共有成像区域越大,此时被测对象与远场目标之间的检测距离较近;轴间距越大,两个光学传感器共有成像区域越小,检测距离则较远,该结果还与光学系统焦距、口径等参数有关。选择100 m左右的通视距离作为被测对象与目标的测量距离,即可满足绝大部分装备检测精度要求,同时,在较短的测量距离下,利用半导体激光器可将测距误差控制在更小范围内,进一步提高测量精度,而轴间距的获得主要依靠被测对象设计图纸,±10 mm的测量误差有较大余量。

5 结论

归纳分析多轴线一致性检测技术,提出了一种基于非合作目标图像处理的多光轴一致性检测技术,为平台光电装备光轴一致性在线检测提供了有效技术途径,该方法实现了大空间跨度机械轴线与光学传感器轴线的自动检测,能够满足复杂光电观瞄系统的野外在线检测需求,同时,结合后续的图像处理技术,可进一步完成复杂轴线检测任务,适用于多种光电观瞄系统轴线检测。

下一步将研究含激光光轴的平台光电装备轴线一致性检测;优化火力轴线可视化过程中的共轴问题,进一步减少配合误差;同时,寻找含红外光轴的平台光电系统中更优的共有特征,优化图像匹配算法,提高测量精度。

参考文献

[1] 凌军, 刘秉琦, 赵熙林. 几种光轴平行性测试方法的比较与探讨[J]. 应用光学, 2003, 24(1): 43-45.

    Ling J, Liu B Q, Zhao X L. The comparison and discussion of several testing methods about optical-axis parallelism[J]. Journal of Applied Optics, 2003, 24(1): 43-45.

[2] 姜宏滨. 用投影靶板调整光轴平行性[J]. 舰船科学技术, 1995(4): 61-65.

    Jiang H B. Adjustment of optical axis parallelism with projected target plate[J]. Ship Science and Technology, 1995(4): 61-65.

[3] 刘爱敏, 高立民, 吴易明, 等. ZEMAX辅助分析斜方棱镜面形误差对出射光平行度的影响[J]. 应用光学, 2009, 30(3): 491-495.

    Liu A M, Gao L M, Wu Y M, et al. ZEMAX auxiliary analysis for effect of rhombic prism surface shape error on parallelism of emergent beam[J]. Journal of Applied Optics, 2009, 30(3): 491-495.

[4] 常山, 曹益平, 陈永权. 五角棱镜的光束转向误差对波前测量的影响[J]. 应用光学, 2006, 27(3): 186-191.

    Chang S, Cao Y P, Chen Y Q. Influence of beam turning error of pentagonal prism on wave-front measurement[J]. Journal of Applied Optics, 2006, 27(3): 186-191.

[5] 马世帮, 杨红, 杨照金, 等. 光电系统多光轴平行性校准方法的研究[J]. 应用光学, 2011, 32(5): 917-921.

    Ma S B, Yang H, Yang Z J, et al. Multi-spectral axes parallelism calibration of electro-optical system[J]. Journal of Applied Optics, 2011, 32(5): 917-921.

[6] 王涛, 田留德, 赵建科, 等. 基于平行光管的空间目标姿态模拟研究[J]. 光学学报, 2018, 38(6): 0612003.

    Wang T, Tian L D, Zhao J K, et al. Research on space target attitude simulation based on collimator[J]. Acta Optica Sinica, 2018, 38(6): 0612003.

[7] 陈志斌, 肖文健, 马东玺, 等. 大间距多光轴一致性野外在线检测方法[J]. 光学学报, 2017, 37(1): 0112006.

    Chen Z B, Xiao W J, Ma D X, et al. A method for large distance multi-optical axis parallelism online detection[J]. Acta Optica Sinica, 2017, 37(1): 0112006.

[8] 李敏, 张昂, 张俊波, 等. 大口径望远镜主镜径向偏心误差检测技术[J]. 激光与光电子学进展, 2017, 54(9): 091103.

    Li M, Zhang A, Zhang J B, et al. Radial eccentricity error detection of primary mirror of large aperture telescopes[J]. Laser & Optoelectronics Progress, 2017, 54(9): 091103.

[9] CI systems. Advanced weapon boresight system (AWBS)[EB/OL].[2018-08-01]. https:∥www.ci-systems.com/CI∥userdata/SendFile.asp?DBID=1&LNGID=1&GID=608.

[10] Cabib D, Rahav A, Barak T. Broad-band optical test bench (OPTISHOP) to measure MTF and transmittance of visible and IR optical components[J]. Proceedings of SPIE, 2007, 6543: 654311.

[11] CarlZeiss. Calibration and alignment[EB/OL]. [2018-08-01].www.zeiss.com/optronics.

[12] SchillReglerteknik. Aligner 308 ship alignment system[EB/OL]. [2018-08-01].www.schill.se.

[13] Jaklitsch JJ, Ehart AF, Jones DA, et al. Gyroscopic system for boresighting equipment: US7065888[P]. 2005-07-14.

[14] Jaklitsch JJ, Paturzo VA. Non line of sight boresight based on inertial measurement technology[C]∥IEEE Systems Readiness Technology Conference, 2003: 527- 533.

[15] Ahmadabadian A H, Yazdan R, Karami A, et al. Clustering and selecting vantage images in a low-cost system for 3D reconstruction of texture-less objects[J]. Measurement, 2017, 99: 185-191.

[16] Cuesta F, Lopez-Rodriguez F, Esteban A. A new blondin system for surveying and photogrammetry[J]. Sensors, 2013, 13(12): 16894-16914.

[17] Metronor. Next generation boresight system[EB/OL].[2018-08-01]. www.metronor.com/military/.

[18] 杨博文. 大型装备装配位姿视觉检测的关键技术研究[D]. 南京: 南京航空航天大学, 2013: 22- 28.

    Yang BW. Research on the key technologies of vision-based assembly pose measurement for large-scale equipments[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2013: 22- 28.

[19] 翁璇, 叶南, 张丽艳. 单目视觉测量系统中目标位姿对图像点灵敏度影响的研究[J]. 机械科学与技术, 2015, 34(6): 969-973.

    Weng X, Ye N, Zhang L Y. Effects of the target position on the sensitivity of image point in monocular visual measurement system[J]. Mechanical Science and Technology for Aerospace Engineering, 2015, 34(6): 969-973.

[20] 张勇, 马东玺, 闫鹏程. 火炮身管轴线通用检测装置及其检测方法: ZL201310191768.7[P].2013-08-14.

    ZhangY, Ma DX, Yan P C. Common detection device and detection method of gun barrel axis: ZL201310191768.7[P].2013-08-14.

[21] Blevins WM, Fischer S G.Boresighting and pointing accuracy determination of gun systems: US8006427[P].2011-08-30.

张勇, 吴浩, 牛刚. 基于非合作目标图像处理技术的大间距轴线一致性检测方法[J]. 激光与光电子学进展, 2019, 56(5): 051203. Yong Zhang, Hao Wu, Gang Niu. Consistency Detection Method for Large Spacing Axis Based on Non-Cooperative Target Image Processing Technology[J]. Laser & Optoelectronics Progress, 2019, 56(5): 051203.

本文已被 3 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!