-
在实际测量中,目标由于受到大气湍流、光学系统像差衍射以及系统跟踪稳定性等因素的影响,常常会对点目标辐射源产生弥散效应,造成目标与背景的混搅,形成弥散模糊图像[5]。虽然此时成像像元数目增多,但由于点目标的能量分散,造成目标在红外靶面的图像出现模糊退化现象,且每个像元的灰度值偏小,目标与背景叠混在一起,辐射特性测量计算时如果将其当作面目标处理,会造成测量精度的降低[6]。因此在靶场测量中,一般将远距离飞行的导弹,飞机与轮船等在红外焦平上成像像元数量小于8×8的目标都考虑作为点目标处理[7]。
在成像观测过程中,理想点目标经过稳定大气传输后,在经过线性不变的光学成像系统,最后到达红外靶面阵列的红外图像模型可以表示为[8]:
$$ G\left( {x,y} \right) = P * S\left( {x,y} \right) + t\left( {x,y} \right) $$ (1) 式中:G(x,y)为目标落在靶面的像元灰度;P为点扩散函数(Point spread function, PSF),是光学系统传递函数的空间域表达;S(x, y)为目标理想图像;t(x, y)为光学成像系统的系统噪声;*表示卷积。
其模糊退化过程如图1所示,其中点扩散函数P模型的建立对目标的模糊退化恢复具有重要的意义,选取不同的模型会导致出现靶面图像的不同呈现状态,特别是会影响的目标与背景灰度位置的界定,进而影响目标的红外辐射测量精度。假设理想条件下红外靶面的能量分布函数为S(x0,y0),则目标辐射能经过点扩散函数P模型作用后,其弥散点目标像素灰度值分布G可以表示为:
$$ G({x}_{0},\text{ }{y}_{0})=K\times \left[P\ast S({x}_{0},\text{ }{y}_{0})\right]+B $$ (2) 式中:K为红外光学系统的响应增益;B为红外光学系统的响应偏置;(x0, y0)为弥散目标成像中心点位置[9]。
对于被测量的理想点目标在红外靶面成像时,所有像素点均可以分为目标像素和背景像素两部分组成。只有准确地识别出点目标成像的边缘位置,区分出目标和背景,才能有效的恢复出目标的真实灰度。通过长时间持续观察发现,目标能量分布规律呈现出连续的偏态样分布形状,联想到工程中常用的瑞利分布可以通过设置参数建立弥散斑中心位置、边界带与灰度值之间的关系,因此确定弥散中心以及边界问题可以归结于概率论中数据集离散度问题。文中在总结之前研究人员经验成果上,提出了基于瑞利分布的判别方法来分辨点目标和背景。
-
通过持续分析观察发现,X(Y)轴方向上不同亮度的点目标在红外靶面灰度值分布曲线图如图2所示的几种情况,从统计学角度来看,曲线整体分布与趋势线符合一个连续的瑞利分布概率密度函数,其定义可表示为:
图 2 不同亮度的点目标X(Y)轴灰度直方图
Figure 2. Gray level histogram of X(Y) axis of point target with different brightness
$$\begin{split} \\ P(x)= \left\{ \begin{array}{l} \dfrac{x}{{s}^{2}}\mathrm{exp}\left(-\dfrac{{x}^{2}}{2{s}^{2}}\right)\text{,}x \gt 0\\ \begin{array}{ccc}0& ,& x \leqslant 0\end{array} \end{array} \right. \end{split}$$ (3) 式中:x为随机变量;σ为参数变量,决定曲线形状的平坦程度,可对右偏斜数据、对称数据等建立数学模型,其期望、方差和得到累计分布函数分别如公式(4)~(6)所示[10]:
$$ E(X) = \int_{ - \infty }^{ + \infty } {xP(x)} {\rm{d}}x = \sqrt {\frac{\pi }{2}} \sigma \approx 1.253\sigma $$ (4) $$ D(X) = E({X^2}) - \left[ {E({X^2})} \right] = \frac{{4 - \pi }}{2}{\sigma ^2} \approx 0.429\sigma $$ (5) $$ L(X) = 1 - \exp \left( {\frac{{ - {x^2}}}{{2{\sigma ^2}}}} \right) $$ (6) 当分别设参数变量σ=0.5、σ=1、σ=2,σ=4时形成的瑞利分布概率密度曲线图如图3所示,其与图2中的目标灰度值分布曲线极为相似,曲线形状与现实状况相匹配。变量x可用来描述点目标像素点能量聚集程度中心位置,距离x越近,目标灰度值越大,成为目标灰度的概率越高;反之,则成为背景灰度的概率会越高;变量σ来描述数据变量集的离散程度,σ值越小曲线形状越窄,数据集越聚集;反之,σ值越大曲线形状越平坦,数据集越分散,符合典型的瑞利分布规律[11]。因此,实际中不同亮度的点目标在红外靶面能量分布情况可以使用瑞利分布来反映,确定能量中心位置与灰度值之间的关系。
对瑞利分布概率密度函数进行积分,得到其累计分布函数,从工程上可以理解为空间占比率,即曲线与数轴所包含的面积,当变量x是固定数值时,σ数值越大,覆盖面积越大,利用上述公式令σ=1,x=2,得到L(X)=87.3%;σ=1,x=3,得到L(X)=98.7%,同理可得到σ
$ \geqslant 2 $ 的情况,如图4所示。那么依据瑞利分布概率的计算,点目标灰度值落在区间(−3σ,3σ)概率为大于98.7%,因此对某像素点的灰度值P(x)可以定义为:如果$ P(x) - 3\sigma \geqslant {G_B} $ 则认为该像素灰度值为目标灰度值;如果$ P(x) - 3\sigma \lt {G_B} $ 则认为该像素灰度值为背景灰度值[12]。因此,可利用瑞利分布来描述实际中不同亮度的点目标弥散情况,确定弥散斑中心位置、目标与背景边界与灰度值之间相互关系。 -
为考核该方法是否具有良好工程适用性和较高的测量精度,通过辐射特性测量精度比对进行实验验证,实验选取对海作用距离试验中400 mm口径中波红外成像系统作为测试载体,其红外辐射测量系统参数设置如表1所示。另采用600 mm×600 mm的标准方形面源黑体作为真值目标,黑体架高固定在近海可移动的靶船上,黑体温度设置从60 ℃开始逐渐升温至130 ℃,每10 ℃设置为一个温度间隔,为减少外部环境影响,试验整体测量时间控制在1 h以内完成。
表 1 红外辐射测量系统参数设置
Table 1. Parameter setting of infrared radiation measurement system
Parameter name Value Aperture/mm 400 Band/μm 3.7-4.8 Pixel number 640×512 Pixel size/μm 15 Bit depth/bit 14 Integration time/μs 1500/2500 图5为实验中波红外实时拍摄目标截取图像,黑色方框内为传统图像二值化处理方法识别出的目标区域,圆形框内为瑞利分布法识别出的目标区域。并给出该基于瑞利分布方法处理该目标在X(Y)方向的灰度分布图,从图中可以看出两者很好的拟合在一起。
图 5 实验中红外目标图像及瑞利分布拟合情况
Figure 5. Infrared target image and Rayleigh distribution fitting in the experiment
考虑目标成像清晰且便于分析与判断的原则,观测距离设置为3.5 km的高山上,经理论计算目标成像约为4×4,满足大于10倍的辐射源尺寸的要求,因此可以将此标准辐射源黑体看成点目标处理。图6为该黑体不同温度时刻呈现在MATLAB上的灰度变化,从图像上可以看出,选取出的点目标在1500 μs积分情况下的70 ℃、90 ℃、110 ℃和130 ℃成像能量分布均呈现出典型瑞利分布形式[13]。而且目标成像灰度值区间在4000~14000之间,满足中波红外成像系统最低灰度可探测值约1500,最高灰度可探测值约15000要求,且成像像元数来说满足点目标测量计算的要求。
为检验基于瑞利分布对弱小目标能量判断方法是否具备较好的测量精度,将应用传统图像二值化处理[14]得到的辐射强度与应用基于瑞利分布图像处理得到的辐射强度进行比对实验,在相同的外部环境条件下,采用同一的红外辐射特性测量设备,选取相同的积分时间和温度,以面源黑体辐射源作为测量真值,利用两种不同方法下辐射强度测量值与黑体标定真值之差作为测量偏差,即测量精度,得到数据结果如表2所示。
表 2 不同计算方法测量结果数据对比
Table 2. Comparison of measurement results of different calculation methods
t/μs Temperature/℃ Blackbody radiation
intensity/×
10−2 W·sr−1Radiation intensity with
binaryzation/×
10−2 W·sr−1Error Radiation intensity with
Rayleigh distribution/×
10−2 W·sr−1Error 1500 60 3.0846 3.5106 13.81% 3.1572 2.35% 70 4.1615 4.5597 9.57% 4.1441 −0.41% 80 5.5209 5.9500 7.77% 5.5857 1.17% 90 7.2124 7.5829 5.14% 7.3901 2.46% 100 9.2898 9.7716 5.19% 9.3119 0.23% 110 11.8104 12.3958 4.96% 11.3885 −3.57% 120 14.8351 15.5974 5.14% 14.8336 −0.01% 130 19.1106 20.5432 4.42% 19.5578 0.70% 2500 60 3.0846 3.4156 10.73% 3.2901 6.66% 70 4.1615 4.4879 7.84% 4.2319 1.69% 80 5.5209 5.8564 6.08% 5.5331 0.22% 90 7.2124 7.5955 5.31% 7.0158 −2.72% 100 9.2898 9.8316 5.83% 9.3011 0.12% 110 11.8104 10.5275 −10.86% 11.9808 1.44% 120 14.8351 15.2391 2.72% 14.1766 −4.43% 130 18.4276 19.2415 4.42% 18.5578 0.70% 通过对比和分析不难发现,应用基于瑞利分布方法的辐射强度测量偏差绝对值可控制在8%以下,而应用传统图像二值化处理方法得到的辐射强度值偏差绝对值高达20%以上;而且从图7也可以看出,应用基于瑞利分布的计算方法得到的测量数据离散程度更小,集中度更高[15]。因此在探测距离较远、大气情况复杂、目标弥散严重、整体图像对比度低等目标与背景边界区别不明显情况下,应用基于瑞利分布目标能量计算方法可有效提高特性测量精度。
Research on energy calculation method of weak and small targets based on Rayleigh distribution
-
摘要: 随着科学技术的不断进步,军事目标也向着小型化、超高速、低可探测度的趋势发展,因而对目标的探测与识别能力也提出了更高的要求。文中在深入分析点目标成像过程及其能量分布的基础上,提出基于偏态瑞利分布点目标能量计算方法,并通过红外辐射特性测量精度对比实验进行验证,得出应用该方法的辐射强度测量偏差可控制在8%以下,而且测量离散程度更小,可有效区分弥散点目标和背景灰度值,证明该方法具有良好的工程适用性,较高的测量精度且广泛的应用前景。Abstract: With the continuous progress of science and technology, military targets are developing towards miniaturization, ultra-high speed and low detectivity, so the detection and identification ability of targets is also put forward higher requirements. Based on the deep analysis of point target imaging process and on the basis of energy distribution, based on abnormal Rayleigh distribution of point target energy calculation method, the infrared radiation characteristics and has a better precision than the experimental verification, it is concluded that the measurement deviation of radiation intensity can be controlled under 8%, and measuring the degree of discrete numerical smaller, can effectively distinguish between the point target and background. It is proved that the calculation method has good engineering applicability, high measurement accuracy and wide application prospect.
-
表 1 红外辐射测量系统参数设置
Table 1. Parameter setting of infrared radiation measurement system
Parameter name Value Aperture/mm 400 Band/μm 3.7-4.8 Pixel number 640×512 Pixel size/μm 15 Bit depth/bit 14 Integration time/μs 1500/2500 表 2 不同计算方法测量结果数据对比
Table 2. Comparison of measurement results of different calculation methods
t/μs Temperature/℃ Blackbody radiation
intensity/×
10−2 W·sr−1Radiation intensity with
binaryzation/×
10−2 W·sr−1Error Radiation intensity with
Rayleigh distribution/×
10−2 W·sr−1Error 1500 60 3.0846 3.5106 13.81% 3.1572 2.35% 70 4.1615 4.5597 9.57% 4.1441 −0.41% 80 5.5209 5.9500 7.77% 5.5857 1.17% 90 7.2124 7.5829 5.14% 7.3901 2.46% 100 9.2898 9.7716 5.19% 9.3119 0.23% 110 11.8104 12.3958 4.96% 11.3885 −3.57% 120 14.8351 15.5974 5.14% 14.8336 −0.01% 130 19.1106 20.5432 4.42% 19.5578 0.70% 2500 60 3.0846 3.4156 10.73% 3.2901 6.66% 70 4.1615 4.4879 7.84% 4.2319 1.69% 80 5.5209 5.8564 6.08% 5.5331 0.22% 90 7.2124 7.5955 5.31% 7.0158 −2.72% 100 9.2898 9.8316 5.83% 9.3011 0.12% 110 11.8104 10.5275 −10.86% 11.9808 1.44% 120 14.8351 15.2391 2.72% 14.1766 −4.43% 130 18.4276 19.2415 4.42% 18.5578 0.70% -
[1] 黄智国. 空间目标地基红外探测技术研究[D]. 长春: 中国科学院大学, 2018. Huang Zhiguo. Research on space target of ground-based infrared detection technology[D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, 2008. (in Chinese) [2] 张琪. 结合边缘检测的图像二值化算法[D]. 长春: 吉林大学, 2011. Zhang Qi. Image binarization algorithm with edge detection [D]. Changchun: Jilin University, 2011. (in Chinese) [3] Liu Zheng, Yang Dezhen, Li Jiangyong, et al. A review of infrared single frame dim small target detection algorithms [J]. Laser & Infrared, 2022, 52(2): 154-162. (in Chinese) doi: 10.3969/j.issn.1001-5078.2022.02.002 [4] 王涛. 复杂空间环境下红外弱小目标高可靠检测方法研究[D]. 中国科学院大学(中国科学院上海技术物理研究所), 2017. Wang Tao. Research on high reliability detection method of infrared dim target in complex space environment [D]. Shanghai: University of Chinese Academy of Sciences (Shanghai Institute of Technical Physics, Chinese Academy of Sciences), 2017. (in Chinese) [5] Ming Yong, Gan Xiaomin, Yang Fan. Infrared weak target detection algorithm based on time-space domain and high-order moment [J]. Foreign Electronic Measurement Technology, 2021(12): 1-6. (in Chinese) [6] 高波. 空间红外弱小目标检测算法研究及实现[D]. 成都: 电子科技大学, 2014. Gao Bo. Research and implementation of space infrared dim target detection algorithm [D]. Chengdu: University of Electronic Science and Technology of China, 2014. (in Chinese) [7] Cao Lihua, Wan Chunming, Zhang Yunfeng, et al. Infrared radiation characteristic measure method of point target [J]. J Infrared Millim Waves, 2015, 34(4): 460-464. (in Chinese) [8] Chang Songtao, Zhang Xiaoyu, Sun Zhiyuan, et al. Radiation measurement method for infrared small target [J]. Acta Optica Sinica, 2014, 34(5): 0504001. (in Chinese) [9] Wang Bing, Gao Qiong, Zhang Yanxiu. Characterization and inversion method of infrared radiation properties of helicopter in point source detection [J]. Infrared, 2018(5): 37-41. (in Chinese) [10] 李周. 地基靶场红外辐射特性测量系统宽动态辐射测量研究[D]. 长春: 中国科学院大学2018. Li Zhou. Research on ground-based infrared characteristics measurement systems in wide dynamic range radiometry [D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, 2008. (in Chinese) [11] Lu Xiaofei. The engineering arithmetic of the reflected infrared radiation characteristic for air-target [J]. Infrared Technology, 2013(5): 289-294. (in Chinese) [12] Han Jinhui, Wei Yantao, Peng Zhenming, et al. Infrared dim and small target detection: A review [J]. Infrared and Laser Engineering, 2022, 51(4): 20210393. (in Chinese) [13] 钱惟贤. 复杂背景下红外小目标探测与跟踪若干关键技术研究[D]. 南京: 南京理工大学, 2010. Qian Weixian. Research on some key technologies of infrared small target detection and tracking under complex background[D]. Nanjing: Nanjing University of Science & Technology, 2010. (in Chinese) [14] 曾雅琼. 复杂背景下红外弱小目标的检测[D]. 南京: 南京理工大学, 2012. Zeng Yaqiong. Detection of infrared dim target in complex background[D]. Nanjing: Nanjing University of Science & Technology, 2012. (in Chinese) [15] Wu Jie, Lin Suzhen, Lu Xiaofei, et al. Infrared dim target detection based on multi-feature fusion [J]. Electronic Measurement Technology, 2022(5): 67-71. (in Chinese)