红外点目标双门限检测方法研究

刘辉, 刘太阳, 王成良, 许松

刘辉, 刘太阳, 王成良, 许松. 红外点目标双门限检测方法研究[J]. 红外技术, 2023, 45(6): 622-629.
引用本文: 刘辉, 刘太阳, 王成良, 许松. 红外点目标双门限检测方法研究[J]. 红外技术, 2023, 45(6): 622-629.
LIU Hui, LIU Taiyang, WANG Chengliang, XU Song. Double Threshold Detection Algorithm for Infrared Point Targets[J]. Infrared Technology , 2023, 45(6): 622-629.
Citation: LIU Hui, LIU Taiyang, WANG Chengliang, XU Song. Double Threshold Detection Algorithm for Infrared Point Targets[J]. Infrared Technology , 2023, 45(6): 622-629.

红外点目标双门限检测方法研究

详细信息
    作者简介:

    刘辉(1982-),男,四川省仁寿县人,讲师,博士,长期从事红外光电和航天遥感领域的教学科研工作。E-mail: 120524618@qq.com

  • 中图分类号: TP391

Double Threshold Detection Algorithm for Infrared Point Targets

  • 摘要: 检测概率、虚警概率和信噪比是红外点目标检测系统的3个核心指标,本文在恒虚警检测准则下,分析了红外点目标单帧检测(单一门限检测)的检测概率、虚警概率和信噪比三者之间的关系。重点分析计算了采用多帧累积检测的双门限检测方法时,第二门限的选取策略,实现了在不提高对系统信噪比要求的前提下,提高检测概率、降低虚警概率的目的,解决了威胁度大、误报和漏报代价高的高价值目标的检测问题。
    Abstract: The probability of detection, probability of false alarm, and signal-to-noise ratio (SNR) are three core indicators of the detection system for infrared point targets. Under the criterion of constant false alarm detection, this study analyzed the relationship between probability of detection, probability of false alarm, and signal-to-noise ratio(SNR) using single frame detection method, also called single threshold detection for infrared point target systems. This study focused on the analysis and calculation of the selection strategy of the second threshold under a multi-frame accumulation detection method (called the double threshold detection algorithm) and achieved a higher probability of detection and lower probability of false alarm with a lower SNR. A good solution was also found for the problem of detecting high-value targets with high threat and cost for false and missed alarms.
  • 非球面按照面形的不同可以分为两类,即旋转对称与非旋转对称非球面,后者通常称为自由曲面,其中轴对称非球面按照子午截面曲线类型,又分为二次非球面和高次非球面[1]。光学非球面元件与光学球面、传统平面元件相比形状多样且具有更大的自由度和灵活性[2]。然而在检测难度上非球面光学元件要远远难于平面及球面光学元件,所以其应用依然无法和球面及平面光学元件的应用规模相比。早期的非球面光学元件大多是简单的小矢高、小口径元件,加工和检测精度相对较低。20世纪八九十年代,在工业发达国家中,中、小型红外透镜的加工精度是表面粗糙度Ra 0.01~0.005 μm,形状精度0.5~0.1 μm[3]。随着基于计算机数控技术的单点金刚石车削、磁流变抛光以及离子束抛光等现代精密加工技术的出现,使加工纳米级面形精度的非球面成为可能,对非球面的检测技术提出了更大的需求和更高的挑战。近二十年来,受红外热成像技术、空间光学等领域的发展和推动,高精度、大矢高、大口径非球面光学元件的设计与加工需求迅速增大[4-8],面形检测成为制约非球面光学元件更大规模应用的关键因素之一,光学非球面的检测技术已成为超精密加工领域的一个重要研究方向。

    通过对光学非球面面形检测技术相关期刊的收集、整理和研究,对目前光学非球面常用的面形检测技术进行了概述。其中,直接面形轮廓法通用性强,应用范围广,世界上主要面形检测仪器生产商均推出了相应检测设备,进一步提高了该技术在光学非球面加工领域应用推广的程度,本文结合典型的仪器设备,对直接面形轮廓法进行较全面的介绍。另外,自由曲面及其面形检测技术的研究在世界范围内方兴未艾,直接面形轮廓法在自由曲面面形检测的应用,使自由曲面的研究重点从设计更多的扩展到实际应用。最后,总结了光学非球面检测技术的现状和发展趋势。

    非球面光学元件是表面形状偏离球面的光学元件[9-10],如图 1所示。一般地,非球面由以下通用方程表示:

    $$ \begin{gathered} Z = \frac{{C{r^2}}}{{\{ 1 + {{[1 - (1 + k)c{r^2}]}^{0.5}}\} }} + \hfill \\ A{r^4} + B{r^6} + C{r^8} + D{r^{10}} + \cdots \hfill \\ \end{gathered} $$
    图  1  非球面光学元件示意图
    Figure  1.  Diagram of aspherical optical elements

    式中:Z表示平行于对称轴的表面轮廓;r表示距非球面对称轴的水平距离;k为非球面系数;AB、…表示非球面修正系数;C为非球面的顶点曲率。

    按照测量原理的不同,非球面光学元件面形检测方法可以分为以下三大类:

    1)几何光线法。几何光线法是用几何学研究光的传播和成像规律的方法,运用光的直线传播原理,采用非接触式测量,包括哈特曼法、刀口阴影法、朗奇法等[11]。哈特曼法无需辅助元件,动态范围大[12]。刀口阴影法设备简单,速度快,结果直观,灵敏度高,适于现场检测[12]。朗奇法制作简单,使用方便,检测精度在微米或亚微米量级[13]

    2)干涉法。干涉法检测是基于光波叠加原理[14],光波在干涉场中产生亮暗交替的干涉条纹,通过对干涉条纹的分析处理来获取被测量零件的相关信息[15]。干涉法是精密抛光时期面形检测的重要方法。干涉法主要包括零位检测与非零位检测,零位检测主要有无像差点法、补偿镜法、全息法等;非零位检测主要有部分补偿法、子孔径拼接法、剪切干涉法等[11]。干涉法对震动和气流扰动敏感。

    3)直接面形轮廓法。直接面形轮廓法使用探针对光学元件的表面进行扫描,以此来获得光学非球面面形相关数据。目前,直接面形轮廓法是非球面面形检测的主流方法,虽然轮廓法相对效率较低,但应用范围广,可用作非球面成型、研磨、抛光、镀膜后各阶段的面形检测[11]

    综上所述,几何光线法检测仪器简单、检测速度较快,但检测精度不如干涉法。干涉法精度高和灵敏度较高、测试范围较广,但对检测设备、环境和操作要求也较高。直接面形轮廓法由于通用性强,操作方便,成为国内外各研究机构和企业广泛使用的光学非球面面形检测方法。然而,面形轮廓检测仪的设计和制造技术基本上都掌握在国外发达国家手中,对我国光学非球面和尖端产品的研制和发展构成很大的隐患。

    光学非球面面形检测方法多种多样,从古老的刀口阴影法,到当今以Luphoscan和NMF轮廓仪为代表的非接触三维轮廓法,无论是传统方法还是现代方法,都是采用离线检测的方式来完成,这样不仅降低效率,而且还由于二次装夹带来的误差影响加工精度。近年来有学者提出了在线检测的方法,在线检测最大的优势是不再需要对工件进行反复装夹,这样不仅可以消除装夹误差、对刀误差和跳动误差等,而且可以节省时间,提高效率。目前,在线检测还未被广泛应用的最主要原因就是测量精度较低,无法达到离线检测的精度水平。通过不断研究,人们总结出影响机床在线检测系统精度的误差来源主要为机床的热误差、几何误差以及测头误差[16],这将是在线检测技术需要解决的问题。

    根据探针测量方式的不同可以分为接触式和非接触式两种,前者使用机械探针,后者使用光学探针。

    采用轮廓法的测量仪器一般称轮廓仪,接触式轮廓仪的测量原理如图 2所示,非球面光学元件置于工作台上,当机械探针接触到光学元件表面后,在轮廓仪系统中设置滑行速度,机械探针匀速滑过元件表面时,传感器感知到测量表面的位置变化,分别在XZ方向上进行采样,探针滑行结束后,采样结果转化为电信号,最后经过处理,电信号转化为数据。这些数据与事先标定的数据进行比对,从而得出被测表面与理论表面的误差。轮廓仪得到的主要参数是面形峰谷值PV和均方根值RMS,以及表面粗糙度等。这种检测方法主要应用于中小口径的非球面光学元件。该方法在测量硫系玻璃和有色金属等较软材料时,会对已加工表面造成划伤,影响工件表面质量。接触式轮廓仪无法对非旋转对称的非球面光学元件面形进行测量。

    图  2  接触式轮廓仪测量原理
    Figure  2.  Measuring principle of contact profilometer

    当今世界上生产轮廓仪的企业主要有英国的泰勒霍普森、德国蔡司、马尔和日本东京精密等公司。泰勒霍普森公司无论在国外还是国内都占据了相当大的市场份额,是此类轮廓仪最著名的生产商,典型产品有PGI 1240和PGI Dimension轮廓仪(如图 3所示)。PGI 1240是一台适用于中小型尺寸的非球面测量仪器,测量最大直径Φ200 mm,可提供12.4 mm的测量范围,垂直分辨率高达0.8 nm。PGI Dimension测量直径超过Φ200 mm,测量几乎不受非球面陡度的限制,边缘测量角度最大达到85°。由于是接触式测量设备,即使测量探头与工件表面的接触应力再小,也有损伤表面的可能性,探头还存在损耗,而且针对不同的面形须配备不同的探头。

    图  3  接触式轮廓仪
    Figure  3.  Contact profilometer

    非接触轮廓仪和接触式轮廓仪的主要区别是将机械探针换成光学探针,对被测件材料的软硬程度和材料表面粗糙度没有选择,可实现对非球面轮廓的无损检测,将非接触式轮廓法成功应用于商业领域的典型产品有荷兰的NMF(如图 4(a)所示)和英国泰勒霍普森公司的Luphoscan(如图 5(a)所示),下面结合这两款产品对非接触式轮廓法进行介绍。

    图  4  NMF轮廓仪及测量原理
    Figure  4.  NMF profilometer and measuring principle
    图  5  Luphoscan轮廓仪及测量原理
    Figure  5.  Luphoscan profilometer and measurement principle

    NMF是基于新一代圆柱测量框架的超高精度三维光学轮廓绝对测量装备,是一款可用于非接触式高精度测量非球面、离轴面、自由曲面、阵列镜等各类面形的革命性测量设备。NMF目前可测量口径350/600/800/1000 mm的光学元件,集通用性、集成化、非接触、快速、高精度等特性于一体,集成差动共焦和双频激光干涉仪的光谱共焦测头具备极强的测试能力,可测量各种研磨、抛光、镀膜元件,高精度运动控制系统和校准系统与非接触光学测头相结合,实现了快速精准的三维光学轮廓形貌绝对测量。

    NMF三维轮廓仪测量原理如图 4(b)所示,被测元件安放在气浮主轴上,主轴以最大1转/s的速度旋转。非接触式共焦测头安装在旋转轴平轴上,测量时测头法向实时垂直于被测表面,同时,测头还可分别沿R轴(径向)和Z轴(垂直方向)移动。

    Luphoscan轮廓仪有260、420等系列,是一款基于多波长干涉传感技术的高速非接触式3D非球面光学面形测量系统,LuphoScan轮廓仪通过激光测距的方法来测量,测量原理如图 5(b)所示,传感器的位置由RZT轴控制,传感器在被测工件表面按照预设路径进行扫描,同时被测工件在C轴上自身进行旋转,传感器对整个面形进行螺旋式扫描。扫描过程中传感器始终垂直于被测面,从而保证了传感器与被测点之间的距离保持恒定。利用多波长干涉传感技术可以灵活地测量绝对距离,即使被测光学非球面元件表面存在缺口也可以实现检测,因为该技术不受光束间断的影响。

    Luphoscan轮廓仪采用光学探针,对镀膜后的光学表面,如果膜层的透过波段和探针的工作波段有重叠,探针发出的光束会在膜层之间发生透射和反射,干扰检测结果,使仪器得出错误的数据,这一定程度限制了此类轮廓仪对镀膜后表面的面形检测。

    光学自由曲面因面形复杂多变,矫正像差能力强,逐渐成为现代光学工程领域的研究热点[18-19]。在非成像领域,将自由曲面应用于汽车后视镜,在镜面大小不变的情况下获得更大的视野,为行车安全提供更大的保障[20]。在显示领域,自由曲面可用于头戴式显示器[21-23]、微型投影仪[24-25]等系统中。在成像领域,自由曲面主要应用于天文观测和空间光学系统中[26-27]

    尽管光学自由曲面有着突出的优点,但依旧没有广泛发展到现代光学系统中,最主要的原因就是光学自由曲面不具备旋转对称性,面形检测成为了最大的难点。前文所述的NMF和Luphoscan不但提供了对自由曲面的检测方法而且是在非球面面形检测技术基础上发展而来的,主要包括三坐标测量法、轮廓仪法、摆臂式轮廓扫描法,以下结合典型仪器设备介绍轮廓仪法在光学自由曲面面形检测的应用。

    自由曲面是常规非球面一个新的分支,从测量原理上来说,自由曲面的轮廓检测和常规非球面完全一致,所不同的是,自由曲面没有旋转对称性,为了实现对自由曲面的轮廓测量,就需要更多的运动维度、更复杂的运动机构和控制精度来支撑其测量精度。前文所述的LuphoScan和NMF两款仪器率先实现了自由曲面面形轮廓检测,推动了自由曲面从理论研究到实际应用的进程。

    LuphoScan是对自由曲面进行快速、非接触、全三维测量的理想设备,LuphoScan260检测自由曲面,需要在扩展软件里输入自由曲面参数,生成4轴(3个机械轴,1个空气轴)轨迹,然后在LuphoScan主界面调用轨迹,按照螺旋线、双波长的方法进行扫描测量。LuphoScan在径向轮廓扫面的基础上,增加了环形轮廓的测量,通过径向与环形轮廓数据的拼接得到整个面形。

    图 6为LuphoScan测量硅单晶自由曲面的实例,工件直径为164 mm。根据自由曲面的数学表达式生成相应的矢高表,测量过程中,信号传感器根据矢高表自动优化运动轨迹,减小测量误差。该轮廓仪只能测量面形与球面基或非球面基偏离度较小的自由曲面。

    图  6  LuphoScan测量硅自由曲面
    Figure  6.  LuphoScan measures silicon free-form surfaces

    NMF是全球第一台实现高精度自由曲面 & 非球面测试能力的设备。NMF系列的OS软件可实现自由曲面、离轴元件居中测量和分析,目前,支持双锥型(Biconic)和多项式型(Polynomial)两种自由曲面输入方式。NMF轮廓仪对自由曲面扫描时有螺旋、往复、步进共焦3种扫描方式,自由曲面测量不确定度<15 nm RMS。NMF轮廓仪对最大口径为500 mm×100 mm的自由曲面的测量不确定度达到30 nm。图 7为NMF测量凸自由曲面的实例[17]。自由曲面偏离度±2 mm,局部斜率11°,测试口径R29.7 mm。

    图  7  NMF测量凸自由曲面[17]
    Figure  7.  NMF measures convex free-form surfaces[17]

    NMF是全采用轮廓法测量自由曲面通用性较强,既可以测量凸面,也可以测量凹面,测量精度和效率较高,然而目前面形轮廓法只能测量梯度变化较小的自由曲面[18, 27]

    随着科技的进步,尖端产品和先进光电系统对光学非球面元件的需求越来越大。一方面,高精度非球面光学元件的制造离不开高精度的非球面检测技术;另一方面,非球面检测技术的发展又反过来促进了非球面制造技术的发展及其应用的普及,可以说非球面检测技术是未来超精密加工领域一个极其重要的发展方向。总结起来,光学非球面面形检测技术呈现出以下特点和发展趋势。

    1)面形轮廓法仍是现在和未来相当一段时期内的主流光学非球面测量方法,其中,非接触无损测量及全口径三维轮廓测量将是该技术的主要发展方向;

    2)检测数据和加工设备共享、输出结果和相应标准兼容已成为标配;

    3)检测设备应具备测量速度快、稳定性好、精度高、检测流程自动化等特点;

    4)在线检测技术将是光学非球面检测技术研究的一个重要方向;

    5)有必要加快国产非球面检测仪器的研制和发展,实现光学非球面应用和发展的自主可控;

    6)降低非球面检测设备的成本,将更好地推动非球面加工技术的发展和非球面的应用普及。

  • 图  1   过门限检测系统的检测概率和虚警概率的定义

    Figure  1.   The definition of detection probability and false alarm probability of crossing threshold detection system

    图  2   恒虚警检测准则下不同虚警概率对应的检测概率和信噪比的关系曲线

    Figure  2.   Relationship between detection probability and SNR corresponding to different false alarm probability under constant false alarm detection criteria

    图  3   双门限检测流程示意图

    Figure  3.   Double threshold detection flow diagram

    图  4   n=5时不同第二检测门限(k)对应的PDPFApdpfa的关系曲线

    Figure  4.   Curves of PD, PFA and pd, pfa while n= 5 at different second threshold

    图  5   双门限检测中图像信噪比SNR与双门限检测准则k/n之间的对应关系

    Figure  5.   Relationship between SNR and double threshold detection criteria k/n in double threshold detection

    图  6   双门限检测的各项技术指标与第二门限之间的对应关系

    Figure  6.   Relationship between the index and second threshold for double threshold detection

    图  7   表 1中几种典型信噪比SNR下的红外仿真图像

    Figure  7.   Infrared simulation image at several typical SNR in table 1

    图  8   虚警概率pfa=0.042的恒虚警单帧检测结果(连续10帧)

    Figure  8.   The detection results of constant false alarm detection while pfa=0.042 at sequential 10 frames

    图  9   6/10双门限检测准则下的检测结果

    Figure  9.   The detection result of double threshold with using 6/10 detection criteria

    表  1   双门限检测的各项技术指标与第二门限之间对应关系的典型计算值

    Table  1   Typical calculated value of relationship between the index and second threshold for double threshold detection

    Consecutive detection
    frames n
    Overall (tactics) indicator PD = 0.98,PFA =1×10-6
    Optimal second
    threshold k
    SNR Single detection
    probability pd
    Single alarm
    probability pfa
    1 1 7.3113 0.9800 1.00×10-6
    2 2 6.0442 0.9899 1.00×10-3
    3 2 5.2419 0.9160 5.77×10-4
    4 3 4.7384 0.9399 6.31×10-3
    5 3 4.3800 0.8647 4.65×10-3
    6 4 4.1050 0.8912 1.62×10-2
    7 4 3.8882 0.8273 1.31×10-2
    8 5 3.7152 0.8528 2.86×10-2
    9 5 3.5644 0.7990 2.44×10-2
    10 6 3.4447 0.8228 4.20×10-2
    11 6 3.3306 0.7768 3.69×10-2
    12 6 3.2391 0.7344 3.30×10-2
    13 7 3.1527 0.7588 4.98×10-2
    14 7 3.0774 0.7219 4.52×10-2
    15 8 3.0109 0.7439 6.25×10-2
    下载: 导出CSV
  • [1] 吴军辉, 朱竟成. 红外成像系统对点目标动态探测概率分析[J]. 红外技术, 2000, 22(1): 40-44. DOI: 10.3969/j.issn.1001-8891.2000.01.010

    WU Junhui, ZHU Jingcheng. An analysis on dynamic detecting probability of IR imaging system for point targets[J]. Infrared Technology, 2000, 22(1): 40-44. DOI: 10.3969/j.issn.1001-8891.2000.01.010

    [2] 郭杰峰. 针对点目标探测信号处理研究[D]. 上海: 中科院上海技术物理研究所, 2009.

    GUO Jiefeng. Research on Signal Processing for Point Target Detection[D]. Shanghai: Shanghai Institute of Technical Physics, Chinese Academy of Sciences, 2009.

    [3] 祁蒙. 红外搜索跟踪系统的探测概率研究[J]. 激光与红外, 2004, 34(4): . 269-271. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW200404009.htm

    QI Meng. Detection probability of IR search and track system[J]. Laser & Infrared, 2004, 34(4): 269-271. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW200404009.htm

    [4] 刘辉. 红外光电探测原理[M]. 北京: 国防工业出版社, 2016.

    LIU Hui. Principle of Infrared Photoelectric Detection[M]. Beijing: National Defense Industry Press, 2016.

    [5] 杨宜禾. 红外系统信号的积累检测问题[J]. 激光与红外, 1980(3): 61-71. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW198003013.htm

    YANG Yihe. The problem of signal accumulation and detection in infrared system[J]. Laser & Infrared, 1980(3): 61-71. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW198003013.htm

    [6] 杨宜禾. 红外系统[M]. 2版: 北京: 国防工业出版社, 1995.

    YANG Yihe. Infrared System[M]. Second edition: Beijing: National Defense Industry Press, 1995.

    [7] 魏宏建, 任浩, 沈同圣. 红外图像探测器噪声的数值模拟方法[J]. 航天电子对抗, 2002(6): 33-35. https://www.cnki.com.cn/Article/CJFDTOTAL-HTDZ200206011.htm

    WEI Hongjian, REN Hao, SHENG Tongshen. Numerical simulation of infrared image detector noise[J]. Space Electronic Countermeasures, 2002(6): 33-35. https://www.cnki.com.cn/Article/CJFDTOTAL-HTDZ200206011.htm

    [8] 钟伟, 刘辉, 许松, 等. 红外点目标探测系统中标准正态分布函数的工程计算[J]. 空军预警学院学报, 2016, 30(2): 102-105. https://www.cnki.com.cn/Article/CJFDTOTAL-KLDX201602008.htm

    ZHONG Wei, LIU Hui, XU Song, et al. Engineering computation of standard normal distribution function in infrared target detection[J]. Journal of Air Force Early Warning Academy, 2016, 30(2): 102-105. https://www.cnki.com.cn/Article/CJFDTOTAL-KLDX201602008.htm

    [9] 小哈得逊. 红外系统原理[M]. 红外系统原理翻译组, 译北京: 国防工业出版社, 1975.

    Hudson Jr P R. Principles of Infrared Systems[M]. Translated by Group on Principles of Infrared Systems, Beijing: National Defense Industry Press, 1975.

    [10] 张帆, 李智. 计算DSP预警系统对弹道导弹探测概率的一种方法[J]. 装备指挥技术学院学报, 2007, 18(3): 67-71. https://www.cnki.com.cn/Article/CJFDTOTAL-XYZH200703014.htm

    ZHANG Fan, LI Zhi. A method of calculating DSP early warning system's detect probability of ballistic missile[J]. Journal of the Academy of Equipment Command & Technology, 2007, 18(3): 67-71. https://www.cnki.com.cn/Article/CJFDTOTAL-XYZH200703014.htm

    [11] 张弘, 赵保军, 史彩成. 对低信噪比下的红外点目标高检测率的研究[J]. 系统工程与电子技术, 2001, 23(3): 58-60. https://www.cnki.com.cn/Article/CJFDTOTAL-XTYD200103018.htm

    ZHANG Hong, ZHAO Baojun, SHI Caicheng. The study of high detecting probability of infrared point target under low SNR[J]. Systems Engineering and Electronics, 2001, 23(3): 58-60. https://www.cnki.com.cn/Article/CJFDTOTAL-XTYD200103018.htm

    [12] 门雷. 红外系统中积累检测技术的理论推算与实施[J]. 红外研究(A辑), 1987(4): 297-302. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH198704009.htm

    MEN Lei. Accumulating detection technique in infrared system-calculation and practice[J]. Chin. J. Infrared Res. , 1987(4): 297-302. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH198704009.htm

  • 期刊类型引用(1)

    1. 肖永平,马世娟,倪晓东. 频率源中关键信号电磁兼容研究. 电子设计工程. 2024(19): 115-118+123 . 百度学术

    其他类型引用(1)

图(9)  /  表(1)
计量
  • 文章访问数:  127
  • HTML全文浏览量:  37
  • PDF下载量:  36
  • 被引次数: 2
出版历程
  • 收稿日期:  2022-04-10
  • 修回日期:  2022-05-03
  • 刊出日期:  2023-06-19

目录

/

返回文章
返回