留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于光学滤波的水面弱纹理提取方法研究

张祥翔 傅雨田 陈永和

张祥翔, 傅雨田, 陈永和. 基于光学滤波的水面弱纹理提取方法研究[J]. 红外与激光工程, 2021, 50(3): 20200215. doi: 10.3788/IRLA20200215
引用本文: 张祥翔, 傅雨田, 陈永和. 基于光学滤波的水面弱纹理提取方法研究[J]. 红外与激光工程, 2021, 50(3): 20200215. doi: 10.3788/IRLA20200215
Zhang Xiangxiang, Fu Yutian, Chen Yonghe. Research on extracting water weak texture based on optical filtering[J]. Infrared and Laser Engineering, 2021, 50(3): 20200215. doi: 10.3788/IRLA20200215
Citation: Zhang Xiangxiang, Fu Yutian, Chen Yonghe. Research on extracting water weak texture based on optical filtering[J]. Infrared and Laser Engineering, 2021, 50(3): 20200215. doi: 10.3788/IRLA20200215

基于光学滤波的水面弱纹理提取方法研究

doi: 10.3788/IRLA20200215
详细信息
    作者简介:

    张祥翔(1981−),女,博士生,研究方向为红外系统成像技术,主要从事星载、机载海洋遥感探测红外技术,系统分析集成和数据仿真等方面的研究

    傅雨田(1967−),男,博士生导师,现任上海技术物理研究所工程九室主任,研究员,博士生导师,上海宇航学会会员。主要从事红外成像遥感技术的研究。曾负责完成国家863项目等相关课题。在国内外杂志上发表研究论文多篇。目前正负责国家重大项目的研究

    陈永和(1982−),男,副主任设计师,现任上海技术物理研究所工程九室主任助理。主要从事空间光学成像遥感技术、低温光学技术、遥感仪器杂散光抑制技术方面的研究

    通讯作者: 陈永和
  • 中图分类号: O43

Research on extracting water weak texture based on optical filtering

  • 摘要: 水下潜体运动产生的内波反映到水表会形成一定的红外弱纹理信号,这使得利用红外遥感手段探测成为可能。但是这种纹理信号的对比度很弱,而且往往和振幅很强的背景杂波混在一起,给信号的提取造成很大困难。比较了已有水面弱纹理探测方法的优缺点,提出光学滤波的方法增强水面弱纹理信号的信杂比,使这种信号的提取更容易。基于非相干光学滤波的原理和特定信杂比评估弱纹理信号提取性能的方法,对不同特定频率信号的光学滤波方法采用信杂比进行分析和评估,然后在此基础上,针对特定频率的信号优化了相应的非相干光学传递函数,通过仿真和实验验证了该方法的可行性,同时表明优化的双孔径非相干光系统对弱纹理模型和信号具有较好的提取能力,该优化的系统可以得到7%以上的信杂比增强,再结合后续的数字对比度增强技术可以进一步有效提高对弱纹理信号的提取能力。
  • 图  1  单孔径及其光学传递函数

    Figure  1.  Single-aperture and its MTF

    图  2  双孔径及其光学传递函数

    Figure  2.  Double-aperture and its MTF

    图  3  双孔径和单孔径光学传递函数比较

    Figure  3.  MTF comparison between double-aperture and single-aperture

    图  4  影响信杂比因素分析

    Figure  4.  Factors impacting SCR

    图  5  信号和杂波的单孔径和双孔径MTF比较

    Figure  5.  MTF comparison for signal and clutter in double aperture and single aperture

    图  6  理想模型双孔径和单孔径比较

    Figure  6.  Ideal model comparison of double aperture and single aperture

    图  7  Kelvin尾迹和弱纹理信号图像模拟

    Figure  7.  Kelvin wake and weak texture signal image simulation

    图  8  光学系统图

    Figure  8.  Optical layout

    图  9  信号和杂波在双孔径和单孔径MTF中的比较

    Figure  9.  Comparison of signal and clutter in MTF of double aperture and single aperture

    图  10  实验装置

    Figure  10.  Experiment setup

    图  11  双孔径光阑,单孔径光阑和分辨率靶标

    Figure  11.  Double aperture diaphragm, single aperture diaphragm and resolution target

    图  12  单孔径和双孔径实际成像效果比较

    Figure  12.  Imaging effect comparison between double aperture and single aperture

    表  1  理想模型信杂比增强

    Table  1.   SCR enhancement for ideal model

    ItemSingle apertureDouble apertureSCR enhancement
    SCR/dB−44.6−34.2+23.3%
    下载: 导出CSV

    表  2  双孔径对Kelvin尾迹的低频和高频部分信杂比增强比较

    Table  2.   SCR enhancement comparison between low frequency and high frequency of Kelvin wake through double aperture

    ItemSingle apertureDouble aperture

    Imaging

    Aperture function

    Non-interference MTF

    Digital contrast
    enhancement
    下载: 导出CSV

    表  3  光学系统设计参数

    Table  3.   Specification of optical system

    ItemParameter of optical systemSpecification
    1Numeric aperture0.1
    2Focal length/mm153
    3Object distance/mm6 000
    4Magnification0.03
    5Camera parameterSensor size 1/3 in
    Pixel size 1.67 μm
    下载: 导出CSV

    表  4  不同孔径信杂比增强比较

    Table  4.   SCR enhancement comparison of different apertures

    Type of apertureSingle apertureDouble aperture
    Imaging
    SCR enhancment-+11.2%
    下载: 导出CSV
  • [1] 杨卫平, 张志龙, 李吉成, 等. 基于缩比模型的水下目标热尾流可探测性研究[J]. 红外与激光工程, 2016, 45(3): 0302002. doi:  10.3788/IRLA201645.0302002

    Yang Weiping, Zhang Zhilong, Li Jicheng, et al. Thermal wakes detectability of submerged objects based on scale model [J]. Infrared and Laser Engineering, 2016, 45(3): 0302002. (in Chinese) doi:  10.3788/IRLA201645.0302002
    [2] 王平, 杜永成, 杨立. 基于重叠网格技术和VOF模型的潜艇热纹理浮升扩散规律的数值与实验研究[J]. 红外与激光工程, 2019, 48(4): 0404002. doi:  10.3788/IRLA201948.0404002

    Wang Ping, Du Yongcheng, Yang Li. Numerical and experimental study on the buoyancy and diffusion laws of submarine thermal wake based on overset grid technology and VOF model [J]. Infrared and Laser Engineering, 2019, 48(4): 0404002. (in Chinese) doi:  10.3788/IRLA201948.0404002
    [3] 来庆志. 海洋分层流中潜航体热尾流扩散特征分析[D]. 哈尔滨: 哈尔滨工业大学, 2019.

    Lai Qingzhi. Analysis on diffusion characteristics of submerged body thermal wake in ocean stratified flow[D]. Harbin: Harbin Engineering University, 2019. (in Chinese)
    [4] 张旭升, 郭亮, 胡日查, 等. 红外探测中潜艇冷热纹理的传热传质特性[J]. 光学 精密工程, 2017, 25(01): 107-114. doi:  10.3788/OPE.20172501.0107

    Zhang Xudong, Guo Liang, Hu Richa, et al. Heat and mass transfer characteristics of submarine cold-thermal wake in the infrared detection [J]. Optics and Precision Engineering, 2017, 25(01): 107-114. (in Chinese) doi:  10.3788/OPE.20172501.0107
    [5] 周哲, 白宗良, 史径丞, 等. 基于温度插值技术的潜艇热尾流浮升规律及海面热特征仿真方法研究[J]. 红外技术, 2019(41): 1025-1032.

    Zhou Zhe, Bai Zongliang, Shi Jingcheng. The study on the simulation method based on temperature interpolation technology for the buoyancy law and ocean-surface thermal characteristic of submarine wake [J]. Infrared Technology, 2019(41): 1025-1032. (in Chinese)
    [6] 陈雄. 潜艇水动力尾迹与热尾迹耦合作用下海面红外特性分析[D]. 南京: 南京理工大学, 2017.

    Cheng Xiong. Infrared chracteristics analysis of the sea surface under coupling effects of the submarine’s hydrodynamic wake and thermal wake[D]. Nanjing: Nanjing University of Science & Technology, 2017. (in Chinese)
    [7] 黄苗苗, 张楠, 朱爱军. 内波作用下水下航行体水动力载荷及运动特性研究[J]. 船舶力学, 2019, 23(5): 531-540. doi:  10.3969/j.issn.1007-7294.2019.05.004

    Huang Miaomiao, Zhang Nan, Zhu Aijun. Hydrodynamic loads and motion features of a submarine with interaction of internal solitary waves [J]. Journal of Ship Mechanics, 2019, 23(5): 531-540. (in Chinese) doi:  10.3969/j.issn.1007-7294.2019.05.004
    [8] 李丹梦, 金伟其, 李力, 等. 水下运动目标的水面波纹数值模拟及分析[J]. 红外与激光工程, 2018, 47(11): 1126004. doi:  10.3788/IRLA201847.1126004

    Li Danmeng, Jin Weiqi, Li Li, et al. Numerical simulation and analysis of free-surface wake generated by moving submerged target [J]. Infrared and Laser Engineering, 2018, 47(11): 1126004. (in Chinese) doi:  10.3788/IRLA201847.1126004
    [9] Zhang Xusheng, Guo Liang, Hu Richa, et al. Cold-thermal wake characteristics of submarine in temperature-density stratified seawater[C]//SPIE Proceedings, 2019, 11023: 1102356.
    [10] 张士成, 杨桢, 杨立, 等. 水下航行体自由表面波浪尾迹红外特征及探测[J]. 红外与激光工程, 2012, 41(10): 2615-2620. doi:  10.3969/j.issn.1007-2276.2012.10.011

    Zhang Shicheng, Yang Zhen, Yang Li, et al. Infrared characterization and detection of free surface wave wake of underwater vehicle [J]. Infrared and Laser Engineering, 2012, 41(10): 2615-2620. (in Chinese) doi:  10.3969/j.issn.1007-2276.2012.10.011
    [11] 徐曼, 裘溯, 金伟其, 等. 基于水面特征波纹的水下运动目标Radon变换探测方法[J]. 光学学报, 2019, 39(10): 1001003. doi:  1001003

    Xu Man, Qiu Su, Jin Weiqi, et al. Radon transform detection method for underwater moving target based on water surface characteristic wave [J]. Acta Optics Sinica, 2019, 39(10): 1001003. (in Chinese) doi:  1001003
    [12] Zhong Rui, Yang Li, Du Yongcheng. Deep transfer learning for underwater vehicle wake recognition in infrared imagery[C]//SPIE Proceedings, 2019, 11333: 11333U.
    [13] 胡建明, 乔凯, 智喜洋. 天基观测条件下复杂环境对空中弱小目标可探测性的影响[J]. 红外与毫米波学报, 2019, 38(3): 351-357. doi:  10.11972/j.issn.1001-9014.2019.03.016

    Hu Jianming, Qiao Kai, Zhi Xiyang. Influence of complex environment on the detectability of weak And small aerial target under space-based observation mode [J]. J Infrared Millim Waves, 2019, 38(3): 351-357. (in Chinese) doi:  10.11972/j.issn.1001-9014.2019.03.016
    [14] 苏显渝, 李继陶. 信息光学[M]. 第2版. 北京: 科学出版社, 2011: 235-268.

    Su Xianyu, Li Jitao. Information Optics[M]. 2nd ed. Beijing: Science Press, 2011: 235-268. (in Chinese)
    [15] 郝未倩, 梁忠诚, 刘肖尧, 等. 分形结构稀疏孔径阵列的成像性能[J]. 物理学报, 2019, 68(19): 199501. doi:  199501

    Hao Weiqian, Liang Zhongcheng, Liu Xiaoyao, et al. Imaging performance of fractal structure sparse aperture arrays [J]. Acta Phys Sin, 2019, 68(19): 199501. (in Chinese) doi:  199501
  • [1] 颜士玲, 朱学艺, 王洪飞, 解碧野, 张秀娟, 卢明辉, 詹鹏, 张利剑, 芦红.  基于非厄米和拓扑效应的光场调控机制与光学器件研究进展(特邀) . 红外与激光工程, 2021, 50(1): 20211001-1-20211001-7. doi: 10.3788/IRLA20211001
    [2] 武锦辉, 凌秀兰, 刘吉, 陈鑫.  缺陷诱导光学薄膜光场增强损伤分析 . 红外与激光工程, 2021, 50(8): 20210357-1-20210357-6. doi: 10.3788/IRLA20210357
    [3] 李东康, 高丽媛, 王涛, 田杏霞, 付长宝.  三束相干光场驱动的Ru原子蒸汽中可调谐的光吸收 . 红外与激光工程, 2020, 49(9): 20190528-1-20190528-5. doi: 10.3788/IRLA20190528
    [4] 顾忠政, 殷达, 聂守平, 冯少彤, 邢芳俭, 马骏, 袁操今.  基于涡旋滤波的图像边缘增强研究进展 . 红外与激光工程, 2019, 48(6): 603015-0603015(14). doi: 10.3788/IRLA201948.0603015
    [5] 张敏敏, 田珍耘, 熊元康, 马凤英, 弓巧侠, 潘志峰, 杜艳丽.  非相干光自干涉数字全息成像技术研究 . 红外与激光工程, 2019, 48(12): 1224001-1224001(8). doi: 10.3788/IRLA201948.1224001
    [6] 钱润达, 赵东, 周慧鑫, 于君娜, 王士成, 荣生辉.  基于加权引导滤波与时域高通滤波的非均匀性校正算法 . 红外与激光工程, 2018, 47(12): 1204001-1204001(6). doi: 10.3788/IRLA201847.1204001
    [7] 赵太飞, 王秀峰, 王花, 余叙叙, 李永明.  弱湍流中紫外光非直视分集接收技术的研究 . 红外与激光工程, 2018, 47(12): 1222002-1222002(6). doi: 10.3788/IRLA201847.1222002
    [8] 赵新亮, 王海霞, 李同海, 崔莉.  部分空间相干光经准均匀介质的散射 . 红外与激光工程, 2017, 46(8): 818003-0818003(6). doi: 10.3788/IRLA201746.0818003
    [9] 李东康, 刘玉娟, 田杏霞, 王涛, 付长宝.  多色相干光学烧孔的动态生成与调控 . 红外与激光工程, 2017, 46(11): 1106007-1106007(4). doi: 10.3788/IRLA201746.1106007
    [10] 潘继环, 苏安, 吕琳诗, 农慧霞, 韦爱珍, 钟发辉, 韦宇森, 高英俊.  四元异质结构光子晶体的可调通道滤波功能 . 红外与激光工程, 2017, 46(S1): 110-114. doi: 10.3788/IRLA201746.S121001
    [11] 唐秀福, 苏安, 何小超, 吕琳诗, 黎玉婷, 高英俊.  四元异质结构光子晶体的光开关功能 . 红外与激光工程, 2017, 46(S1): 120-124. doi: 10.3788/IRLA201746.S121003
    [12] 米仁杰, 万助军, 汪涵.  MEMS可调谐平顶窄带光学滤波器 . 红外与激光工程, 2016, 45(7): 720001-0720001(5). doi: 10.3788/IRLA201645.0720001
    [13] 孔英秀, 柯熙政, 杨媛.  空间相干光通信中本振光功率对信噪比的影响 . 红外与激光工程, 2016, 45(2): 222002-0222002(6). doi: 10.3788/IRLA201645.0222002
    [14] 许轰烈, 陈钱, 顾国华, 张玉珍.  利用导向滤波的宽动态范围图像增强技术 . 红外与激光工程, 2015, 44(12): 3843-3849.
    [15] 张珂卫, 汪伟, 赵卫, 谢小平.  基于FFT和多重信号分类算法的高精度相干光相移键控信号频率偏移估计算法研究 . 红外与激光工程, 2015, 44(5): 1593-1597.
    [16] 王利国, 吴振森, 王明军, 王万君, 张耿.  湍流大气中部分相干光二阶统计特性的三参数模型及其应用 . 红外与激光工程, 2015, 44(1): 317-320.
    [17] 赵知诚, 沈为民, 唐敏学, 贺虎成.  非相干叠加激光发射系统光学设计 . 红外与激光工程, 2014, 43(6): 1728-1734.
    [18] 贺明, 王亚弟, 贺菁, 陈向军, 路建方, 寇广, 沈鹏.  双边滤波直方图均衡的非均匀性校正算法 . 红外与激光工程, 2013, 42(10): 2620-2624.
    [19] 方翔, 江新光, 吴逢铁, 程治明, 范丹丹.  非相干LED白光产生无衍射光的光源设计 . 红外与激光工程, 2013, 42(1): 159-162.
    [20] 刘华松, 傅翾, 王利栓, 姜玉刚, 冷健, 庄克文, 季一勤.  弱吸收单面薄膜光学特性的表征方法 . 红外与激光工程, 2013, 42(8): 2108-2114.
  • 加载中
图(12) / 表(4)
计量
  • 文章访问数:  432
  • HTML全文浏览量:  96
  • PDF下载量:  31
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-11-15
  • 修回日期:  2020-12-20
  • 网络出版日期:  2021-05-12
  • 刊出日期:  2021-03-15

基于光学滤波的水面弱纹理提取方法研究

doi: 10.3788/IRLA20200215
    作者简介:

    张祥翔(1981−),女,博士生,研究方向为红外系统成像技术,主要从事星载、机载海洋遥感探测红外技术,系统分析集成和数据仿真等方面的研究

    傅雨田(1967−),男,博士生导师,现任上海技术物理研究所工程九室主任,研究员,博士生导师,上海宇航学会会员。主要从事红外成像遥感技术的研究。曾负责完成国家863项目等相关课题。在国内外杂志上发表研究论文多篇。目前正负责国家重大项目的研究

    陈永和(1982−),男,副主任设计师,现任上海技术物理研究所工程九室主任助理。主要从事空间光学成像遥感技术、低温光学技术、遥感仪器杂散光抑制技术方面的研究

    通讯作者: 陈永和
  • 中图分类号: O43

摘要: 水下潜体运动产生的内波反映到水表会形成一定的红外弱纹理信号,这使得利用红外遥感手段探测成为可能。但是这种纹理信号的对比度很弱,而且往往和振幅很强的背景杂波混在一起,给信号的提取造成很大困难。比较了已有水面弱纹理探测方法的优缺点,提出光学滤波的方法增强水面弱纹理信号的信杂比,使这种信号的提取更容易。基于非相干光学滤波的原理和特定信杂比评估弱纹理信号提取性能的方法,对不同特定频率信号的光学滤波方法采用信杂比进行分析和评估,然后在此基础上,针对特定频率的信号优化了相应的非相干光学传递函数,通过仿真和实验验证了该方法的可行性,同时表明优化的双孔径非相干光系统对弱纹理模型和信号具有较好的提取能力,该优化的系统可以得到7%以上的信杂比增强,再结合后续的数字对比度增强技术可以进一步有效提高对弱纹理信号的提取能力。

English Abstract

    • 水下潜体的运动产生的内波反映到水表会形成一定的红外弱纹理信号,使得红外探测成为可能[1],研究表明,水下潜体运动产生的水面弱纹理信号主要有V型纹理、涡型纹理等,它们的特点是背景杂波强、振幅弱,呈现一定空间频率分布[2-8]。目前所采用的光学探测手段主要是使用相机获取海面纹理图像,再利用数字滤波算法识别目标纹理[9-12],所用相机一般为单孔径光瞳。由于这种纹理信号的振幅弱,对比度很低,往往和海浪等振幅很强的背景杂波混在一起,给信号的提取造成很大困难,尤其是随着潜体降噪、深潜技术的发展,这些纹理信号的探测难度更大。

      根据这些研究基础,文中提出了一种基于非相干光学滤波提升水面弱纹理信号提取能力的方法。该方法针对水面弱纹理目标特性,采用特定的信杂比评价方式,在非相干光学滤波原理的基础上优化了适合特定频率信杂比增强的光学传递函数,提高了弱纹理信号相关分量在整个图像频谱中的占比,从理论和实验上验证了基于非相干光学滤波增强弱纹理信号信杂比的可行性。

    • 弱纹理探测可以采用信杂比进行评价[13]。目标信号是具有特定频率分布的弱纹理,其他频率分量则属于杂波,表示目标纹理的物理量ζsignal除以表示杂波纹理的物理量ζclutter,用分贝(dB)表示,称之为信杂比SCRtexture,如公式(1)所示:

      $$ {\rm{SCR}}_{\rm texture}=10\lg\left(\dfrac{{\zeta}_{\rm signal}}{{\zeta }_{\rm clutter}}\right) $$ (公式(1))

      对图像直接进行分析时采用信号方差σs2和杂波方差σc2分别表示信号和杂波的物理量,并用分贝(dB)表示信杂比SCRtexture,如公式(2)所示:

      $$ {\rm{SCR}}_{\rm texture}=10\lg\left(\frac{{{\sigma }_{\tilde {\rm s}}}^{2}}{{{\sigma }_{\rm c}}^{2}}\right) $$ (公式(2))

      信杂比越高,说明信号越强,背景越弱,这样信号提取就更容易。提高信杂比,可以通过提高信号分量,也可以通过抑制杂波分量。

    • 用光学遥感手段进行海面纹理探测是基于海面辐射的被动探测,属于非相干光探测。

      非相干光探测中的光学滤波通过对非相干光学传递函数(MTF)进行调制,可以在一定程度上提取特定频率的信号。对衍射受限系统,非相干光学传递函数H(u,v)就是光瞳函数P(ξ,η)的自相关,所以衍射受限的非相干光学传递函数可以表示为[14]

      $$ H\left(u,v\right)=\dfrac{{ \displaystyle\iint }_{-\infty }^{+\infty }P(\xi ,\eta )P(\xi +\alpha ')\lambda l'(v+\beta '){\rm d}\alpha '{\rm d}\beta '}{{ \displaystyle\iint }_{-\infty }^{+\infty }{\left[P\right(\xi ,\eta \left)\right]}^{2}{\rm d}\xi {\rm d}\eta } $$ (公式(3))

      普通单孔径系统,如图1(a)所示,其MTF用归一化频率表达如公式(4)所示,MTF随归一化频率γ变大而减少,如图1(b)所示。

      图  1  单孔径及其光学传递函数

      Figure 1.  Single-aperture and its MTF

      $$ {\rm MTF}_{\rm single}\left(\gamma \right)=\frac{2}{\pi }\left[{\rm arccos}\gamma -\gamma {\left(1-{\gamma }^{2}\right)}^{\frac{1}{2}}\right]\;\;\;\;\;0\leqslant \gamma \leqslant 1 $$ (公式(4))

      根据公式(1),具有特定频率γtypical的目标纹理的物理量用MTFsingledγ表示,其他频率分量为杂波,其物理量用Ssingle表示,所以单孔径系统的信杂比SCRsingle可以表示为:

      $$ {\rm SCR}_{\rm single}=10 {\rm lg}\frac{{\rm MTF}_{\rm single}{\rm d}\gamma }{{S}_{\rm single}} $$ (公式(5))

      根据非相干光学滤波原理,适当改变光瞳的形状,可以提高目标频率的相对占比,减少非目标频率分量,也就是杂波分量占比,从而提高信杂比。减少杂波分量,最理想的是使杂波分量对应的传递函数为0,由于在非相干光系统中零频分量无法减少,所以主要考虑使其他非零频杂波分量尽量减少。

      多孔径光瞳可以抑制部分非零频杂波分量,其MTF函数MTFmulti-aperture[15]

      $$ \begin{split} &{\rm MTF}_{\rm multi-aperture}={\rm MTF}_{\rm d}+\\ & \frac{1}{N}{\rm MTF}_{\rm d}*\sum _{j=1}^{N-1}\sum _{k=j+1}^{N}\delta \left(\xi \pm \frac{{x}_{j}-{x}_{k}}{\lambda f},\eta \pm \frac{{y}_{j}-{y}_{k}}{\lambda f}\right) \end{split} $$ (公式(6))

      式中:MTFd为多孔径中每一个子孔径的MTF;$ {x}_{j}- $$ {x}_{k} $$ {y}_{j}-{y}_{k} $为子孔径的相对位置;N为子孔径的数目,其中N≥2。

      由公式(6)可知,当N=2时,即采用双孔径系统时,如图2(a)所示,选取合适的子孔径直径及其相对位置,在目标频率处可以得到最大50%的MTFm,同时可使部分非零频杂波分量的传递函数为0,如图2(b)所示。对应的双孔径系统非相干光学传递函数MTFdouble如公式(7)所示:

      图  2  双孔径及其光学传递函数

      Figure 2.  Double-aperture and its MTF

      $$ \begin{split} {\rm MTF}_{\rm double}=&{\rm MTF}_{\rm d}+\\ & \frac{1}{2}{\rm MTF}_{\rm d}\left(\xi +\dfrac{{x}_{j}-{x}_{k}}{\lambda f}\right){\rm MTF}_{\rm d}\left(\xi -\dfrac{{x}_{j}-{x}_{k}}{\lambda f}\right) \end{split} $$ (公式(7))

      根据公式(1),双孔径系统中目标频率的物理量可以用MTFmdγ表示,杂波可以用Sdouble表示,因此双孔径信杂比可以表示为:

      $$ {\rm SCR}_{\rm double}=10{\rm lg}\dfrac{0.5×{\rm d}\gamma }{{S}_{\rm double}}=10{\rm lg}\dfrac{0.5×{\rm d}\gamma }{2{\left(\dfrac{{d}_{\rm double}}{{D}_{\rm single}}\right)}^{2}{S}_{\rm single}} $$ (公式(8))

      式中:Dsingle为单孔径直径;ddouble为双孔径中子孔径直径,0.5为双孔径的最大MTF,即MTFm

      因此在目标频率${\rm{\gamma }}_{\rm typical}$下双孔径系统相对于单孔径系统信杂比的增强比例△S为:

      $$ \begin{split} \vartriangle S=&\dfrac{{\rm SCR}_{\rm double}-{\rm SCR}_{\rm single}}{\left|{\rm SCR}_{\rm single}\right|}=\\ &10{\rm lg}\dfrac{0.5}{{2{\left(\dfrac{{d}_{\rm double}}{{D}_{\rm single}}\right)}^{2}{\rm MTF}_{\rm single}}}/{|\rm SCR}_{\rm single}| \end{split} $$ (公式(9))

      由于ddouble/Dsingle≤0.5,MTFsingle≤1,因此ΔS>0,也就是双孔径相对单孔径总有信杂比增强。图3(a)是双孔径对中低频信号的信杂比增强,图3(b)是对中高频信号的信杂比增强,阴影部分是双孔径可以抑制的杂波频率分量。γ0.5是MTF为MTFm即0.5时的频率。

      图  3  双孔径和单孔径光学传递函数比较

      Figure 3.  MTF comparison between double-aperture and single-aperture

    • 根据以上分析,通过选取双孔径中子孔径的相对位置和大小,可以对特定频率的信号提高信杂比。基于该原理对理想模型进行仿真,然后对弱纹理模型进行验证。

    • 影响双孔径系统信杂比增强的因素有双孔径中子孔径直径和归一化频率。

      图4(a)所示,双孔径中归一化子孔径直径(ddouble/Dsingle)越小,信杂比增强越大;如图4(b)所示,相对于截止频率的归一化频率(γtypical/γcutoff)越高,信杂比增强越大。

      图  4  影响信杂比因素分析

      Figure 4.  Factors impacting SCR

      遥感相机为了获得清晰的图像,通常希望系统的MTF能够达到0.2以上,这要求相机光学系统的MTF不低于0.4。MTF越高,所需要的口径越大。文中以目标MTF达到0.4所需要的口径作为单孔径系统的参数。

      双孔径系统中双孔径的间距L由目标频率γtypical决定,如公式(10)所示:

      $$ {\gamma }_{\rm typical}=\dfrac{{d}_{\rm double}}{L+{d}_{\rm double}}×\dfrac{L}{\lambda F} $$ (公式(10))

      式中:F为镜头焦距;λ为系统波长;ddouble为双孔径中子孔径直径。

      子孔径直径的选取决定杂波抑制程度和积分时间。子孔径直径越小对杂波的抑制越有效,但所需的积分时间越长,图像的噪声越大;子孔径直径越大,杂波抑制程度越低,信杂比增强越少。综合考虑两方面因素,选取子孔径直径为单孔径直径的0.25。

    • 以周期50 lp/mm, 灰度20阶的信号作为目标信号,以20、 40、80、94 lp/mm为杂波信号。先对信号和杂波进行MTF分析。信号50 lp/mm在双孔径中略有增强,杂波20 、80、94 lp/mm在双孔径中基本全部抑制,40 lp/mm的杂波也有很大程度的抑制,如图5所示。

      图  5  信号和杂波的单孔径和双孔径MTF比较

      Figure 5.  MTF comparison for signal and clutter in double aperture and single aperture

      然后对此信号和杂波进行Matlab仿真。图6(a)是原始信号和杂波,图6(b)是经过双孔径的成像图,图6(c)是经过单孔径的成像图。和MTF分析相对应,在双孔径中20、80、94 lp/mm的杂波基本全部抑制,40 lp/mm的杂波也大部分抑制,信号50 lp/mm得以保留。

      图  6  理想模型双孔径和单孔径比较

      Figure 6.  Ideal model comparison of double aperture and single aperture

      最后采用公式(2)对图像的信杂比进行评价,双孔径相对于单孔径信杂比增加了23%,如表1所示。

      表 1  理想模型信杂比增强

      Table 1.  SCR enhancement for ideal model

      ItemSingle apertureDouble apertureSCR enhancement
      SCR/dB−44.6−34.2+23.3%
    • 文中基于参考文献[10]模拟潜艇航行时产生的水面纹理(周期大约为20 m),图7(a)示出考虑海浪和风速对纹理对比度的削弱作用,叠加上更多复杂的低频和高频背景杂波,图7(b)更真实模拟潜艇水面弱纹理信号,其特点是背景杂波信号强,纹理信号振幅弱,纹理有一定的频率分布。

      图  7  Kelvin尾迹和弱纹理信号图像模拟

      Figure 7.  Kelvin wake and weak texture signal image simulation

      根据公式(2)对图7(b)进行信杂比计算,得到双孔径相对于单孔径信杂比增强7%。从表2可以看出,第一行是图7(b)分别经过单孔径和双孔径的仿真图,第二行是单孔径和双孔径函数,第三行是孔径函数的自相关,即非相干光学传递函数,第四行是经过数字对比度增强的效果。可见双孔径相比单孔径可以更有效地提取微弱的纹理信号,再结合数字对比度增强使得弱纹理更加清晰。

      表 2  双孔径对Kelvin尾迹的低频和高频部分信杂比增强比较

      Table 2.  SCR enhancement comparison between low frequency and high frequency of Kelvin wake through double aperture

      ItemSingle apertureDouble aperture

      Imaging

      Aperture function

      Non-interference MTF

      Digital contrast
      enhancement
    • 为了验证双孔径光瞳滤波的效果,结合第3节的仿真进行了实验。

    • 探测器上的目标频率50 lp/mm,对应实际纹理的频率1.5 lp/mm,根据该条件设计孔径光阑,见图8

      图  8  光学系统图

      Figure 8.  Optical layout

      光学系统关键参数如表3所示(表中,1 in=2.54 cm)。

      表 3  光学系统设计参数

      Table 3.  Specification of optical system

      ItemParameter of optical systemSpecification
      1Numeric aperture0.1
      2Focal length/mm153
      3Object distance/mm6 000
      4Magnification0.03
      5Camera parameterSensor size 1/3 in
      Pixel size 1.67 μm

      将信号和杂波在该光学系统的MTF图中进行比较。物面信号频率1.5 lp/mm在双孔径中略有增强,0.5~1.0 lp/mm的杂波分量基本抑制,其余杂波分量也都有一定程度的抑制,如图9所示。

      图  9  信号和杂波在双孔径和单孔径MTF中的比较

      Figure 9.  Comparison of signal and clutter in MTF of double aperture and single aperture

      实验装置如图10所示,将光学系统(Optical system)对准分辨率靶标(reticle)成像,分辨率靶标用光源照明(LED light source)照亮,光学系统接收靶标(reticle)漫反射成像,用黑白相机(BW camera)拍摄。

      图  10  实验装置

      Figure 10.  Experiment setup

      在光学系统孔径光阑的位置,插入单孔径光阑和双孔径光阑进行对比成像,如图11(a)图11(b)所示,旋转双孔径光阑得到目标分辨率线对横方向或竖方向最清晰像。

      图  11  双孔径光阑,单孔径光阑和分辨率靶标

      Figure 11.  Double aperture diaphragm, single aperture diaphragm and resolution target

      实验靶标采用Edmund的53714灰度渐变靶标。选取第11个分辨率图形(OD为0.49),以第0组第5个线对(1.5 lp/mm)为目标,对应的像面频率就是50 lp/mm,其余频率都是杂波,如图11(c)所示。

    • 根据公式(2)对实验图像进行信杂比评价,双孔径相比于单孔径的信杂比增强为11.2%,见表4

      表 4  不同孔径信杂比增强比较

      Table 4.  SCR enhancement comparison of different apertures

      Type of apertureSingle apertureDouble aperture
      Imaging
      SCR enhancment-+11.2%

      对3.2节理想模型、3.3节弱纹理模型和本节的实验靶标,其信杂比增强不同,主要是由于杂波的频率相对于双孔径抑制杂波的频率区间不同导致的,如果杂波正好落在0.5~1.0 lp/mm或者2.0~3.0 lp/mm的频率区间,则会有比较好的信杂比增强效果,如果杂波落在0~0.5 lp/mm或者1.0~2.0 lp/mm,则信杂比增强效果减弱。

      3.2节中的理想模型,杂波分量在2.0~3.0 lp/mm的频率区间比较多,所以信杂比增强比较高。3.3节的弱纹理模型,杂波分量在0~0.5 lp/mm的频率区间比较多,所以信杂比增强较弱。此节中的杂波分布比较均匀,因此信杂比增强在两者之间。

      为了验证非相干光学滤波在实际波纹中的信杂比增强和杂波抑制效果,用螺旋桨搭建水波模型,如图12(a)所示 。用图8的光学系统拍摄单孔径和双孔径的成像图,分别是图12(b)图12(c)所示。目标频率为红色框内的纹理,杂波为绿色框内的纹理。双孔径相对于单孔径红框内的目标纹理得以保留并略有增强,绿框内的杂波纹理基本抑制。

      图  12  单孔径和双孔径实际成像效果比较

      Figure 12.  Imaging effect comparison between double aperture and single aperture

    • 文中根据水面弱纹理的特征,采用信杂比评价方式,提出了一种基于非相干光学滤波提升信杂比的方法,然后通过仿真和实验,验证了该方法的可行性。根据仿真和实验结果,弱纹理模型的信杂比有7%以上的提升。再结合后续的数字对比度增强技术可以有效提高对弱纹理信号的提取能力。

参考文献 (15)

目录

    /

    返回文章
    返回