面阵探测器点目标探测随机误差度量方法研究 下载: 929次
1 引言
激光链路的构建是空间光通信的核心技术之一。空间光通信系统的信号光束散角较小,需要利用捕获(Acquisition)、跟踪(Tracking)、瞄准(Pointing)系统(ATP系统)[1]来实现激光链路的构建。ATP系统采用面阵探测器检测目标位置,目标信标光斑在面阵探测器上的位置变化与跟踪机构形成闭环控制,实现了外界扰动下的目标精确跟踪指向。因此,信标光斑的定位误差在终端光学系统焦距确定的情况下会直接影响目标位置的确定,进而影响整个ATP系统的探测精度[2]。
由于距离遥远,空间光通信中面阵探测器上接收的信标光斑目标一般信噪比较低,此外由于背景辐射、大气湍流,探测器材料和制作工艺等因素的影响[3-4],理论计算出的星点目标质心位置和实际质心位置之间存在误差,且这些误差按照来源可简单分为系统误差和随机误差两部分[5-6]。系统误差本质上是星点定位算法中采用像元几何中心坐标代替像元真实能量中心坐标进行计算而引入的误差。系统误差的分布呈现以像元尺寸大小为周期的规律性变化,因此找到该变化规律即可对其进行补偿[5]。随机误差主要来源于成像过程中产生的各种噪声,如响应非均匀性噪声、读出噪声、暗电流噪声和散粒噪声等。非均匀性响应噪声可以在后期通过两点校正、多项式校正、逐像素点校正等方法[7-11]进行补偿, 但这些方法得出的响应非均匀性受到光子散粒噪声、暗噪声、读出噪声等大量随机因素的影响,无法最大程度地进行平场校正[12]。其他随机噪声一般在设计图像采集电路时采取针对性的措施或采用暗场图像校正、多幅图像单向素时域均值处理等方法予以尽量去除[5,12],也可通过一些滤波和阈值算法[13-14]将目标信号与背景噪声分离以进一步提高信噪比。但是,基于目前的硬件发展水平以及信标光斑星点小目标的特征,这些方法在增加计算开销、提高系统复杂度的同时,容易使光斑质心产生偏移从而引入更大的定位误差[15]。
目前,已经实现了1000 km级的星地量子通信[16],但未来需要发展到10000 km及更远距离以提高实用化水平,链路指向精度将由目前的μrad提高到100 nrad以下,目标探测的随机误差将会成为ATP系统不可忽略的关键问题。为了对面阵探测器的目标定位随机误差进行高精度定量化的衡量,本文面向星地量子通信对链路建立的需求,考虑多种导致信标光斑位置探测随机误差的因素,包括探测器本底噪声、光照强度和光斑尺寸等,建立了质心定位的随机噪声等效角模型,典型工况仿真证明了此数学模型的合理性。
2 目标位置探测原理与信号模型
2.1 质心定位算法
质心定位算法具有计算简单、稳定性好的优点,是星点目标探测领域最常见的定位算法。在该方法中,探测器平面内信标光成像的实际位置等效于像平面的能量重心。由于面阵探测器像元大小的有限性,所得到的信标光斑图像中的每个像元的响应值实际是面阵探测器经模/数转换得到的采样值[17],因此质心定位算法可表示为
式中:D(xi,yj)是面阵探测器上横、纵坐标分别为xi、yj的像元响应值;s为参与星点信标光斑质心计算的空间区域;i、j分别为横、纵坐标轴方向的像元编号;N1,N2分别为横、纵坐标轴方向的像元个数。
2.2 噪声信号的分布模型
引起探测随机误差的噪声按照能量分布形式可进一步描述为均匀分布噪声、高斯分布噪声和泊松分布噪声。均匀分布噪声主要指量化噪声,是在图像由模拟输出电压转化为数字电压时引入,一旦探测器和模拟数字转换器( Analog-to-Digital Converter,ADC)模块电路确定后,量化噪声就为一定值。高斯分布噪声主要包括暗电流噪声和读出噪声[18]。泊松分布噪声主要包括光子散粒噪声和暗电流散粒噪声。
由于泊松噪声是一种信号依赖噪声,随着信号的增强而增大,且平均光子噪声等于入射光子数的平方根;高斯噪声是由光敏器件中电子热随机运动产生的一种加性白噪声,无论光照强度是多少,噪声的平均水平不变。因此为了便于分析,本文忽略在探测器和ADC模块电路确定时固定不变的量化噪声和可以在后期补偿的响应非均匀性噪声,仅考虑高斯分布噪声和泊松分布散粒噪声。
2.3 信标光斑点扩散函数模型
如果不考虑弥散斑的影响,根据目标面积是否可以充满观测相机的分辨率视场,将空间目标在探测器上的成像分为点目标或面目标。如果目标的面积不能充满相机瞬时视场的覆盖范围,空间目标的成像为点目标,否则成像为面目标[19]。根据相机实际所采集到的信标光斑图像,如
图 1. 原始星点图和星点灰度直方图。(a)原始星点图;(b)星点灰度直方图
Fig. 1. Original star spot images and star spot gray histogram. (a) Original star spot images; (b) star spot gray histogram
由于远距离下的光信号传输受光学系统中孔径光阑衍射效应的影响,将
式中:(x0,y0)为光斑能量的中心位置;rx、ry分别为横、纵坐标轴方向上的点扩散函数半径;D0为光斑的中心能量值。
3 随机噪声等效角模型的建立
3.1 传统噪声等效角模型
噪声等效角(NEA)模型是将目标位置探测的随机误差等效为目标角度信息的度量。工业界中已有的NEA模型[20]为
式中:NE为噪声等效角;RSN为信标光斑在焦平面上成像的功率信噪比;FS为探测器的斜率系数,与光斑的空间强度分布和空间尺寸有关。若光斑直径为θbeam,对于艾里光斑,FS=4.14/θbeam;对于高斯光斑,FS=1.27/θbeam。
3.2 基于质心定位的NEA模型
因为面阵探测器得到的信标光斑图像是D(xi,yj)在像元位置的采样值,所以各像元之间互不相关或相关系数很小,各像元的随机误差可近似是相互独立的[17]。因此,面阵探测器上信标光斑位置探测的随机误差可认为是由各像元探测的随机误差叠加而成的,进而可以按照函数随机误差计算方法[21]建立另一种形式的NEA模型。以x坐标方向上的NEA为例,令
式中:f(·)为噪声等效角和探测器像元响应值之间的函数关系。基于质心定位算法[17]原理,
式中:var(·)是光斑质心计算的随机误差;
式中:R1(xi,yj)是坐标为(xi,yj)的像元的散粒噪声;RF(xi,yj)是坐标为(xi,yj)的像元的读出噪声;RT(xi,yj)是坐标为(xi,yj)的像元的暗电流噪声;Δt为曝光时间。根据散粒噪声的信号依赖性,坐标为(xi,yj)的像元的散粒噪声所带来的随机误差近似为像元响应值大小[18],因此坐标为(xi,yj)的像元的随机误差可进一步表示为
此外,对于面阵探测器接收到的信标光斑图像,无论是艾里分布还是高斯分布,其像元响应值均满足以下条件:当中心位置像元的响应值最大时,与中心位置距离相同的像元的响应值相同;与中心像元所在的行或列距离相同的两行或两列像元的响应值相同。基于像元响应值分布的对称性,可对所建立的NEA模型作进一步化简。假设横、纵坐标轴方向上的光斑直径d均为 2N+1个像元,每个像元的坐标等于其坐标的下标,即xi=i,yj=j,并令
则有
由此可见,所推导的NEA计算模型仅与信标光斑半径N、目标信号强度S以及噪声
在特殊条件下,当曝光时间极短,暗电流导致的随机误差可以忽略不计,每个像元的高斯白噪声都为均值且目标光斑能量分布为均匀性分布时,所建立的NEA模型又可精简为
此种情况下容易得出,N和var[RF(xi,yj)]越大NEA越大,S越大NEA越小。但此种情况多适用于利用地球热图像作为探测目标以实现深空光通信终端的捕获、定位、跟踪功能的方案中[22-24]。
3.3 NEA模型对比
综上可知,基于质心定位的NEA模型在通用性和随机误差求解精度方面优于传统NEA模型。因为在传统NEA模型中,求解目标探测随机误差时不仅需要信噪比(RSN)、光斑直径(θbeam)等参数,还需要根据不同的目标光斑点扩散函数模型实时调整FS。当抖动或大气湍流等因素导致光斑图像模糊或形状不规则时,此种模型下还没有合适的方法计算FS参数,如果仍然按照某种典型分布方式计算FS,将会为NEA的求解带来误差。基于质心定位的NEA模型不需要根据信标光斑点扩散函数模型的变化实时调整某个或某些参数,其仅与光斑直径d、信号总强度S以及噪声
4 仿真计算与分析
为了验证所建立的信标光斑质心定位NEA模型的合理性,本文对不同光斑点扩散函数模型、不同光照强度、不同噪声水平和不同光斑尺寸条件下的NEA进行仿真。同时为了覆盖实际面阵探测器接收到的信标光斑的大小和信号强度,仿真参数的设置如
图 4. 三种信号能量分布模型。(a)均匀分布;(b)高斯分布;(c)艾里分布
Fig. 4. Three signal energy distribution models. (a) Uniform distribution; (b) Gaussian distribution; (c) Airy distribution
表 1. 求解NEA的仿真参数设置
Table 1. Simulation parameter setting for NEA calculation
|
1) 不同模型下的NEA都由于高斯噪声的加入而增大,噪声越大NEA越大。同一噪声下,不同模型的NEA数值不同。艾里分布的NEA小于高斯分布和均匀分布,噪声越小不同信号分布下的NEA差值越大,最大差值为29.8%。
2) S值变化时,三种模型的NEA变化趋势一致但数值不同,艾里分布的NEA数值小于其他两种分布。不同分布的NEA最大差异为18.2%,S越大差异越小。
3) d值变化时,三种模型的NEA变化趋势仍然一致。当d值小于5 pixel时不同模型的NEA数值几乎相同;当d值大于5 pixel时,艾里分布的NEA数值小于其他两种分布,且d值越大不同信号分布下的NEA差值越大,最大差值为25.7%。
图 5. 单个随机影响因素下的NEA。(a)噪声水平;(b)信号强度;(c)光斑大小
Fig. 5. NEA under single random influencing factor. (a) Noise level; (b) signal intensity; (c) spot diameter
根据上述分析,得到三种目标信号点扩散函数模型下的NEA随光斑直径d和光照信号强度S的变化趋势,如
1) 噪声条件不变时,不同点扩散函数模型下的NEA都随着d的增加而增加,随着S的增加而减少,因此减小d或增加S时NEA的数值会减小。
2) S值较小时,NEA随d值变化的幅度明显;S值较大时,NEA随d值变化的幅度不明显。因此,在受到光学衍射极限的限制,光斑尺寸不可改变的情况下,在S值较小范围内通过增大信号强度以降低NEA是可取的。
3) 不同模型下的NEA都由于噪声的加入而增大,故减少噪声也可以降低NEA。
4) 不同目标信号分布的NEA数值不同,最大差异小于30 %。当需要满足高帧频高计算速率时,可按照均匀分布模型下的噪声等效角求解公式计算出初始的NEA,加上相应模型间的差异值即为实际所需求解的NEA数值。
仿真结果与工业界NEA模型随信噪比和光斑尺寸变化的结果一致,符合探测器信噪比高或光斑直径小则信标光斑探测精度高的特性,间接证明了所建立的NEA模型的合理性。
5 结论
综合影响信标光斑位置探测的多个随机因素,建立了衡量光斑质心定位随机误差大小的数学模型NEA。仿真分析了不同因素对NEA的影响以及NEA随外界条件的变化趋势。结果表明,当光斑尺寸小于5 pixel,信号强度大于25000 e-1,高斯噪声低于10 e-1时,三种模型的NEA均小于0.01 pixel;信号强度越大,噪声越小,信标光斑半径越小,随机误差则越小,反之亦然,这符合面阵探测器信噪比高或光斑直径小则目标探测精度高的特性。
此外,所建立的NEA模型仅需信号总强度、光斑尺寸和高斯噪声三个参数,不需要根据光斑信号点扩散函数的改变实时调整参数,不仅为随机误差的抑制和消除提供了一种更加便捷的理论评定方法,而且对于模糊或形状不规则的光斑图像,仍可利用此种方法计算随机噪声等效角。所提随机噪声等效角模型可为更高精度ATP系统的发展提供一定参考。
[2] 曹阳, 李保权, 李海涛, 等. 基于星点像重采样的星敏感器高精度质心算法[J]. 光学学报, 2019, 39(7): 0712003.
[3] 陈少杰, 张亮, 王建宇. 探测器条状噪声对精跟踪系统光斑定位的影响[J]. 中国激光, 2017, 44(9): 0906008.
[4] 曹明华, 武鑫, 王惠琴, 等. Gamma-Gamma大气湍流下超奈奎斯特光通信系统性能[J]. 中国激光, 2020, 47(9): 0906003.
[5] 姜亮, 张立国, 张星祥, 等. 星敏感器星点定位系统误差补偿[J]. 中国激光, 2015, 42(3): 0314001.
[6] 杨君, 张涛, 宋靖雁, 等. 星点质心亚像元定位的高精度误差补偿法[J]. 光学精密工程, 2010, 18(4): 1002-1010.
[7] 董建婷, 杨小乐, 董杰, 等. 面阵CMOS图像传感器响应非均匀性校正方法研究[J]. 半导体光电, 2012, 33(4): 570-574, 578.
[11] 杨敏, 徐文斌, 田禹泽, 等. 面向运动目标检测的分时型红外偏振成像系统[J]. 光学学报, 2020, 40(15): 1511001.
[12] 孙慧, 徐抒岩, 孙守红, 等. 光电成像传感器光子响应非均匀性噪声评价方法研究[J]. 激光与光电子学进展, 2015, 52(4): 042302.
[13] 许威. 星点快速提取与高精度定位技术研究[D]. 浙江: 浙江大学, 2013.
XuW. Star fast extraction and high accuracy centroid estimation of star camera[D]. Zhejiang: Zhejiang University, 2013.
[16] Yin J, Li Y H, Liao S K, et al. Entanglement-based secure quantum cryptography over 1, 120 kilometres[J]. Nature, 2020, 582(7813): 501-505.
[17] 陈少杰, 张亮, 吴金才, 等. 空间激光通信中精跟踪系统的实现与优化[J]. 红外与毫米波学报, 2018, 37(1): 35-41.
[18] 刘南南. 高精度导星测量系统细分定位技术研究[D]. 北京: 中国科学院大学, 2014.
Liu NN. Research on hyper accuracy subdivision position technology of fine guidance sensor[D]. Beijing: University of Chinese Academy of Sciences, 2014.
[19] 翟永立. 星空背景下空间运动目标自主识别技术研究[D]. 北京: 中国科学院大学, 2018.
Zhai YL. Research on automatic identification of space motion targets under starry background[D]. Beijing: University of Chinese Academy of Sciences, 2018.
[20] 姜会林, 佟首峰, 张立中. 空间激光通信技术与系统[M]. 北京: 国防工业出版社, 2010: 243- 244.
Jiang HL, Tong SF, Zhang LZ. The technologies and systems of space laser communication technologies and systems of space laser communication[M]. Beijing: National Defense Industry Press, 2010: 243- 244.
[21] 费业泰, 陈晓怀, 秦岚, 等. 误差理论与数据处理[M]. 北京: 机械工业出版社, 1999: 61- 63.
Fei YT, Chen XH, QinL, et al.Error theory and data processing[M]. Beijing: China Machine Press, 1999: 61- 63.
[22] HamkinsJ, Simon MK. Deep-space communications and navigation series[M]. Hoboken: John Wiley & Sons, Inc., 2006.
[23] S.Lee, G. G.Ortiz, W. T.Roberts, et al., tracking, pointing using earth thermal images for deep-space ka-band and opticalcommunications[EB/OL]. ( 2003-11-15)[2020-06-11]. https:∥ipnpr.jpl.nasa.gov/progress_report/42-155/155E.pdf.
Article Outline
刘负负, 王旭, 张亮. 面阵探测器点目标探测随机误差度量方法研究[J]. 光学学报, 2021, 41(4): 0404001. Fufu Liu, Xu Wang, Liang Zhang. Measurement Method of Random Errors in Spot Target Detection by Flat-Panel Detector[J]. Acta Optica Sinica, 2021, 41(4): 0404001.