查 月,曹玉音
(1.中国人民解放军92941部队44分队,辽宁 葫芦岛 125001;2.宜昌测试技术研究所,湖北 宜昌 443003)
20世纪以来,红外成像技术有了很大的发展,并迅速在精确制导武器上得到了广泛应用[1-2]。红外成像制导武器是利用探测目标和背景间微小的温差或辐射频率差引起的热辐射分布情况来实现对目标的跟踪,能够实现在复杂战场中对目标的全向攻击[3]。然而,如何准确评估红外成像导引头在干扰条件下的工作效能,为红外成像导引头抗干扰效能的提高及其验收鉴定提供重要的依据,是一项亟待完成的工作[4-6]。
国内外目前评估系统作战效能的主观评估方法主要有专家调查法、层次分析法、菲尔德法;客观评估法主要有因子分析法、回归分析法、主成分分析法;综合评估法主要有神经网络评估法、灰色评估法、模糊综合评估法[7-14]。目前,对红外成像导引头进行抗干扰效能评估多是利用层次分析法和模糊综合评估法,得到综合抗干扰概率指标。李凡等利用多层模糊分析法对红外制导的抗干扰效能进行了评估,构建了抗干扰效能评估指标体系,提高了评估结果的可信度和准确性[8];钱欢等使用多层次灰色评估模型评估了电子对抗装备不同工作状态下的作战效能[15];葛辰杰等利用线性支持向量机对红外导引头的抗干扰效能进行评估,适用于小数据样本情况[16]。
由于作战场景下的红外干扰因素复杂多样性,对于红外成像制导系统抗干扰效能的评估方法也尚未规范统一[17-18]。本文提出了基于目标识别、图像特征、跟踪精度的多维度反舰导弹抗红外干扰效能评估体系,并通过仿真实验,验证了两种红外成像识别跟踪算法抗干扰表现与体系评估指标的一致性。
红外导引头的工作窗口波段通常是指3~5 μm的中波红外和8~12 μm的长波红外,红外感知系统的工作原理是通过特殊的材料如熔融石英辅以冷光栅和冷屏来捕获红外波长内的热能量。根据已知的材料,可以得到某红外成像探测器的黑体温度与相应灰度值之间的关系,根据目标温度,可以采用黑体的辐射公式计算出黑体的光谱辐射M如下:
(1)
式中,c1为第一辐射常数,c2为第二辐射常数。
红外成像探测器的工作波长为[λ1,λ2],可以根据其工作波长得到黑体的光谱辐射W为:
(2)
辐射照度L的计算公式为:
(3)
即单位面积上的辐射能量。从此处可以得到黑体的辐射照度L和红外成像探测器红外灰度图灰度值的对应关系。最终可以通过普朗克定律和兰伯特定律计算得出远距离目标的探测辐射亮度如下:
(4)
式中,h为普朗克常数,c为光速,k为波尔兹曼常数,T为温度。
红外成像制导系统主要由三部分构成,即红外成像系统、图像处理系统及随动系统。目标的红外辐射信号首先进入到光学系统,通过扫描机构在红外探测器的光敏面上聚焦,经过光电转换将光信号转换为电信号,再对电信号进行滤波和放大处理,得到目标的位置误差信息;或由成像系统解析图像,经过非均匀校正和预处理后得到原始的目标红外图像,再传输至随动系统联动光学系统,引导导引头向目标方位运动,形成角跟踪回路,实现导引系统对目标的连续跟踪。同时生成自动驾驶仪所需的导引信号,指挥控制导弹飞向目标[11]。红外成像制导系统组成如图1所示。
图1 红外成像制导系统构成
红外成像导引头在在识别跟踪过程中具有三种工作状态,分别为检测状态、跟踪状态以及判断状态。
1)检测状态:在检测状态下,红外成像导引头对目标生成的红外图像进行识别检测,寻找所要追踪的目标红外信号。如果导引头发现符合目标特性的红外信号,会通过识别算法确定目标信息并锁定其位置,在确定该信号有效后保存目标信息与位置并转入跟踪状态,否则保持检测状态。
2)跟踪状态:在跟踪状态下,红外成像导引头锁定目标并调整自身方向使得导弹朝向目标所在方位偏移,从而实现改变飞行轨迹与飞行目的地,同时导引头也将记录目标物体的红外成像特征,用于出现抗干扰判断。根据自身采取的跟踪算法对目标置信度进行评估,出现置信度过低或者丢失目标信息的情况,则说明出现干扰,导引头即转入判断状态。
3)判断状态:红外成像导引头判断状态主要出现于多个真假目标物体同时出现造成跟踪算法受到干扰无法正常工作,或者目标信息丢失导致跟踪算法无法继续跟踪目标的情况。此时跟踪算法难以确定跟踪目标,需要通过红外成像导引头抗干扰算法重新选择要追踪的目标物体或者返回检测状态。
红外成像导引头抗干扰效能评估是在复杂战场环境下对红外成像制导系统的捕获、识别跟踪目标的水平进行评估的准则和方法。它对于在复杂海战场环境试验中如何构设红蓝双方光电对抗的博弈态势有着重要指导作用。
评估红外成像导引头效能首先构建合适的评估指标,其次建立数学模型并设置典型场景,通过仿真实验对系统进行抗干扰效果评估,最后根据实验结果对抗干扰效能进行分析评价。通过对红外成像制导工作原理的分析,可以从目标识别、图像特征、跟踪精度三个维度构建评估指标体系[19-20]。
红外成像导引头处于搜索检测工作状态时,对其视场内的目标进行识别,通过目标红外辐射特征判断其类型,并给出置信度指标。目标抗干扰识别能力是指红外成像导引头运用相应的目标识别算法,在打击目标的过程中,从其所采集到的每一帧图像中识别出正确目标的能力[4]。目标识别概率是反映导引头由探测转入截获过程的核心参数,易受自然环境、人为干扰影响,主要表现为规定时间内导引头识别概率下降将导致跟踪目标错误。
平均目标识别概率Pt是指在干扰环境下对检测目标的置信度,定义为:
(5)
式中,pi为单次目标识别的置信度;n为导引头检测识别总次数。
红外成像导引头的核心功能是目标探测、识别和跟踪,而其探测、识别和跟踪的效能依赖于目标图像特征的强弱[22]。当导引头受到干扰时,会导致目标在海天背景下的红外成像质量下降,目标与背景的对比度降低。因此可以根据目标红外成像特征的对比度变化来评估导引头的抗干扰能力[20]。目标成像对比度是指目标成像平均灰度值与背景平均灰度值之比,计算公式如下:
(6)
式中,σ1为目标成像平均灰度值,σ2为背景平均灰度值。
红外成像导引头在截获目标后转入对目标的自动跟踪状态,并连续测量目标的运动方位,控制导弹指向目标。在目标跟踪阶段,目标实施光电干扰措施或进行紧急机动都会对红外成像制导系统的跟踪精度造成影响,具体体现在对目标水平或俯仰方位上的跟踪角度误差。红外导引头方位/俯仰跟踪角度误差定义如下:
(7)
(8)
式中,θ1为目标与导弹连线方位角,θ2为导引头光轴指向方位角,α1为目标与导弹连线俯仰角,α2为导引头光轴指向俯仰角,n为跟踪次数。
红外成像制导系统在跟踪角误差逐渐累积,或目标置信度严重下降的情况下,导引头可能会丢失目标,打断连续跟踪状态而转入判断状态。因此,平均连续跟踪时间比是评价导引头持续抗干扰和稳定跟踪能力的重要指标之一。平均连续跟踪时间比定义为:
(9)
式中,Ti为导引头单次连续跟踪时间,n为目标丢失次数,n为跟踪次数,T0为导引头开机时刻,Ts为导引头销毁时刻。
本文使用的光电对抗仿真平台基于场景想定数据创建作战任务,在整个任务仿真过程中,根据红外环境信息以及目标、目标的红外辐射特征进行建模与红外可视成像仿真,实时进行目标识别和跟踪,并且模拟测试红外成像导引头对抗面源红外干扰信号的情况下打击目标的能力。光电对抗仿真系统通常由场景想定模块、红外环境模块、平台运动模块、红外成像导引头模块和干扰弹模块组成,其中红外成像导引头模块为核心组件,主要具有红外成像、目标识别和目标跟踪功能。
红外成像导引头模块通过光栅化渲染的方式,获取场景中的所有物体位置以及物体对应的模型文件数据,读入其中的网格体数据,然后通过坐标变换到对应三维空间中,通过将红外成像导引头看作虚拟相机进行渲染成像,最终成像结果由物体的温度信息以及位置信息,结合红外辐射成像计算结果进行赋值。
整个仿真场景主要由天空和背景、多个潜在目标物体以及可能的干扰目标组成。红外可视渲染视图与三维渲染视图的视角一致,因此图像中各个物体位置信息以及轮廓等特征均保持一致,区别在于图像颜色。红外可视渲染场景中不考虑光照信息、物体纹理信息等关于颜色信息的资源与数据,而是通过场景中的物体与环境的温度信息,通过真实辐射亮度的计算公式,代入场景当前温度信息,渲染得出贴合真实物理规则的红外灰度图像。最终红外成像模块所得到的信息即为红外场景渲染图,红外图像无RGB信息,以灰度图的形式保存。
由于仿真生成的红外图像信息与真实红外成像采到的数据之间的差异,以及使用原始的红外图像直接传入识别、追踪算法中效果无法达到最佳,因此在信号处理仿真模块中添加高斯噪声、白噪声以及滤波处理操作,将处理后的红外图像数据提供给目标识别跟踪模块,可以获得更好的目标识别、跟踪效果。
红外成像导引头识别目标的过程是先找出目标和背景的差异,提取目标特征,将目标特征进行比较,选则最匹配结果特征,再进行决策分类处理。从兼顾效果与效率上进行选择,在目标识别算法上主要采用YOLO和SSD算法。
1)YOLO算法:YOLO算法(you only look once)是一个one stage的实时目标检测算法,采用单个神经网络直接预测类别概率和物体边界,融合了目标类别和目标区域预测,能够实现端到端的物体检测。YOLO相对于基于Region Proposal如R-CNN的算法准确度较低一些,但它的检测速度更快。YOLO算法训练过程的流程如图2所示。
图2 YOLO算法训练流程
YOLO算法首先使用卷积神经网络提取图像特征,其次预测候选区域的类别和位置。设定每个预测框为一个样本,标注真实框的相对类别和位置并得到标签值,通过网络模型预测其类别和位置并比较网络预测值和标签值,建立损失函数。在图像上依据一定规则产生一系列的候选区域,再根据图像上物体的真实框和候选区域间的位置关系标注候选区域。标注与真实框接近的候选区域为正样本,并设定真实框的位置为正样本的位置目标。标注偏离真实框较大的候选区域为负样本,负样本不用预测类别或位置。
2)SSD算法:SSD(single shot multibox detector)与YOLO一样,也是一个one stage的实时目标检测算法,译为“单激发多框探测器”,思路是在图像多层特征图上密集均匀地抽样,可以采用不同尺度和长宽比的抽样方式,通过CNN提取特征后直接分类和回归。
SSD沿用了YOLOv1的基本思路,糅合了Faster R-CNN的anchor-box,又开创性地提出了多分辨率预测,整个过程只需要一步,所以其优势是速度快,精度高于YOLO,速度胜于faster R-CNN。实时SSD300模型的运行速度,比最初版本的实时YOLO更快,同时产生了显著优越的检测精度。其训练流程如图3所示,最后通过极值抑制得到置信度最高的目标。训练还涉及到选择一些默认框检测尺度,负难例挖掘,数据增强策略。
图3 SSD算法训练流程
SSD预测过程和训练过程相比较为简单,首先将待预测的图片输入网络,将初始框根据输出的目标类别和偏移量进行对应调整,得到预测框,并确定预测框的类别,然后过滤掉属于背景的预测框和阈值较低的预测框。对于留下的预测框执行NMS算法,将重叠度较大的预测框再次过滤掉,留下的预测框为最终的检测结果。预测过程如图4所示。使用中先输入图像信息,将得到置信度最高的目标作为识别目标,若置信度不足则认为识别失败。
图4 SSD算法预测流程
目标跟踪算法主要运用到了Byte Track以及基于Mean Shift的带有模板匹配思路的目标追踪算法。其中Mean Shift算法仅支持单目标追踪,但具有较强的抗干扰能力;Byte Track算法鲁棒性强,能完成多目标追踪。
1)Byte Track算法:Byte Track是基于 tracking-by-detection范式的跟踪方法,是一种简单高效的数据关联方法Byte。它是利用跟踪轨迹和检测框之间的相似性,保留住高分检测结果,并除去低分检测结果中的背景,识别出需要的目标,Byte Track能够降低目标漏检并提高目标轨迹的连贯性。由于Byte Track没有采用外表特征进行匹配,所以跟踪的效果非常依赖检测的效果,若检测器的效果很好,跟踪也会取得不错的效果,若检测的效果不好,就会严重影响跟踪的效果。为降低该算法缺陷,采用的检测算法都是基于神经网络准确率更高的算法,在跟踪出现置信度不足时会进入判断状态,重新使用准确度更高的检测算法重新确定跟踪目标信息。
2)Mean Shift算法:Mean Shift算法是一种基于密度的非参数聚类算法,思路是假设不同簇类的数据集合不同的概率密度分布,寻找样本数据点密度增大最快方向,样本密度高的局部区域对应于该分布局部极值,样本点会在局部密度极值处收敛,收敛到相同局部极值的点即是同一簇类的成员。
状态初始化阶段,算法的主要工作是对导引头探测到的图像进行初步处理,以提高下一步中导引头追踪目标的成功率。在这一阶段,除去一般的像素展开操作外,算法还会使用包括二值化、众数滤波等一系列的传统图像处理手段。在初始化步骤中完成红外图像的初步处理之后,导引头便根据预设参数更新自身的跟踪制导行为。通过均值漂移操作和四向延拓之后,导引头将标记出目标的追踪框图,进而调整导弹的运动方向等行动参数。
通过光电对抗平台仿真软件构设两种实验场景,一种为无干扰场景,一种为红外诱骗干扰场景。在两种场景下分别对红外成像导引头的“SSD+Mean Shift”(以下称算法1)及“YOLO+Byte Track”(以下称算法2)识别跟踪算法进行数字仿真验证,获取效能评估指标样本数据。仿真设置红外成像制导导弹速度为0.7 Ma,距离目标5 km,目标直航速度为20 kn,并通过执行红外面源诱骗干扰战术对红外成像导引头实施干扰,干扰设置参数如表1所示。
表1 仿真参数设置
仿真实验过程中实时记录仿真时间、目标识别概率、目标成像灰度值、目标跟踪方位/俯仰角误差数据,仿真结果如图5~8所示。
图5 目标识别概率仿真曲线
图6 平均目标成像对比度仿真曲线
图7 目标跟踪方位角误差仿真曲线
图8 目标跟踪俯仰角误差仿真曲线
4.2.1 目标识别概率分析
从目标识别概率仿真曲线可以看出,两种算法在无干扰场景下都能够对目标进行准确识别,平均识别概率可以保持在80%以上。面对干扰对抗场景,算法1对红外诱饵的识别反应有一定延迟,并且在仿真时间15 s左右被干扰物引偏,致使目标脱离导引头视场范围,仿真后段持续跟踪干扰物,识别概率维持在57%;算法2在仿真时间第5 s时对红外诱饵完成抗干扰识别,第8 s后继续转向跟踪目标,识别概率稳定在70%以上。通过以上分析可知,算法2的抗干扰识别能力强于算法1。
4.2.2 平均目标成像对比度
在红外成像灰度图中,红外诱饵和目标相对于海天背景具有较高的灰度值,并且诱饵灰度值大于目标。从平均目标成像对比度仿真曲线可以看出,在无干扰场景下,导引头稳定跟踪目标,随着弹目距离的减小,目标成像在识别跟踪框中的图像占比提高,因此平均成像对比度稳定增大。面对干扰对抗场景,算法1导引头被引偏后跟踪诱饵,由于诱饵红外辐射强度更大,短时间能具有较高的成像对比度,随着距离的接近和诱饵的消散,海天背景凸显,成像对比度急剧下降;算法2在干扰对抗期间,识别跟踪框内被面源诱饵占据,成像对比度增大,摆脱干扰后成像对比度曲线与无干扰场景几乎吻合。
4.2.3 目标跟踪方位角误差
导引头跟踪目标的过程中,图像处理系统实时计算弹目方位和俯仰误差,并输入导弹随动系统,导弹根据自身的制导方法输出机动指令,常见的有比例引导和前置角引导等方法。仿真曲线中体现了在无干扰场景下,面对目标,导弹跟踪目标的方位角误差按照导引律稳定变化。面对干扰对抗场景,算法1导引头被引偏后跟踪诱饵,由于诱饵红外在水平方位几乎没有移动速度,因此算法1的跟踪方位角误差在0°附近振荡;算法2在干扰对抗期间,跟踪方位急剧摆动,干扰对抗成功后稳定跟踪目标,此时方位角误差按照导引律稳定变化。通过以上分析可知,跟踪方位角误差并不是越小越好,而是要符合相应的制导规律。
4.2.4 目标跟踪俯仰角误差
从目标跟踪俯仰角误差曲线中可以看出,导引头稳定跟踪目标时几乎没有俯仰角误差,那是因为目标在高度俯仰方向没有运动趋势,而被诱饵诱骗后的算法1导引头,受诱饵自由落体运动的影响,在俯仰指向上出现较大的误差摆动。通过以上分析可知,跟踪俯仰角误差越小,越能体现导引头识别跟踪算法的优越性。
仿真数据表明,算法2的平均目标识别概率大于算法1,因此算法2的目标识别能力更强;由于跟踪诱饵时目标成像对比度会先增大后减小,产生较大变化,因此成像对比度需要根据实时数据进行分析,其平均值难以表明导引头的抗干扰效能;两种算法的平均跟踪角度误差结果也如分析结论一致,容易被引偏并跟踪水平方位几乎静止的诱饵的算法1具有更小的平均跟踪方位角误差,而跟踪俯仰高度稳定的目标的算法2具有更小的平均跟踪俯仰角误差,说明算法2的目标跟踪能力更强;综上所述,算法2具有更好的抗干扰效能。
本文从目标识别、图像特征和跟踪精度三个维度构建评估了效能评估指标体系。通过构设典型场景,对两种红外成像识别跟踪算法进行了红外诱骗干扰对抗仿真实验,通过对抗仿真结果与效能评估指标的数据表现,对比分析可知算法2具有更好的抗干扰效能,验证了效能评估体系的有效性。