2.793

2018影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

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

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

OTH雷达图像的粗糙度指标及用于射频干扰自适应抑制

罗忠涛 郭人铭 郭杰 何子述 卢琨

罗忠涛, 郭人铭, 郭杰, 何子述, 卢琨. OTH雷达图像的粗糙度指标及用于射频干扰自适应抑制. 自动化学报, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286
引用本文: 罗忠涛, 郭人铭, 郭杰, 何子述, 卢琨. OTH雷达图像的粗糙度指标及用于射频干扰自适应抑制. 自动化学报, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286
Luo Zhong-Tao, Guo Ren-Ming, Guo Jie, He Zi-Shu, Lu Kun. Tamura coarseness for OTH radar image evaluation and its application in adaptive optimization of interference suppression. Acta Automatica Sinica, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286
Citation: Luo Zhong-Tao, Guo Ren-Ming, Guo Jie, He Zi-Shu, Lu Kun. Tamura coarseness for OTH radar image evaluation and its application in adaptive optimization of interference suppression. Acta Automatica Sinica, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286

OTH雷达图像的粗糙度指标及用于射频干扰自适应抑制


DOI: 10.16383/j.aas.c190286
详细信息
    作者简介:

    重庆邮电大学通信与信息工程学院讲师, 主要研究方向为信号与信息处理, 统计信号处理, 数字图像处理.E-mail: luozt@cqupt.edu.cn

    重庆邮电大学通信与信息工程学院硕士研究生, 主要研究方向为信号与信息处理, 数字图像处理.E-mail: s180131036@stu.cqupt.edu.cn

    重庆邮电大学通信与信息工程学院学士, 主要研究方向为数字图像处理.E-mail: guojie1072050774@163.com

    电子科技大学电子工程学院教授, 主要研究方向为阵列信号处理, 自适应信号处理, MIMO雷达与通信.E-mail: zshe@uestc.edu.cn

    中国电子科技集团有限公司第十四研究所研究员级高级工程师, 主要研究方向为超视距雷达系统设计和信息处理.E-mail: mimimomoba@gmail.com

  • 基金项目:  国家自然科学基金 (61701067, 61702065), 重庆市教育委员会科研基金 (KJ1600427, KJ1600429)资助

Tamura Coarseness for OTH Radar Image Evaluation and Its Application in Adaptive Optimization of Interference Suppression

More Information
  • Fund Project:  Supported by National Natural Science Foundation of China (61701067, 61702065), Project supported by Scientific Research Foundation of the Chongqing Education Committee (KJ1600427, KJ1600429)
  • 摘要: 针对OTH雷达距离-多普勒(Range-Doppler, RD)图, 本文首次提出采用纹理粗糙度作为RD图质量的评价指标, 即计算RD图所转化灰度图的Tamura纹理粗糙度. 分析表明, 粗糙度指标能准确反映RD图受干扰情况, 对于不同灰度转换函数具有稳健性. 作为应用举例, 本文将图像粗糙度用于改进射频干扰抑制算法, 使干扰抑制达到自适应优化. 实验结果表明, Tamura粗糙度能够正确反映RD图干扰抑制情况, 优化粗糙度指标能够使干扰抑制自适应达到最优.
  • 图  1  OTH雷达RD图

    Fig.  1  RD maps of OTH radar

    图  2  不同干扰情形下的RD灰度图: (a)-(c)分别表示无动态范围转换的无干扰、有窄带射频干扰、有宽带射频干扰的RD灰度图; (d)-(f)分别表示对应的动态范围为80 dB转换的RD灰度图

    Fig.  2  RD gray-scale images: (a)-(c) are images without dynamic range conversion, for no RFI, narrowband RFI, and wideband RFI, respectively; (d)-(f) are corresponding images with dynamic range conversion for 80 dB

    图  3  不同干扰抑制结果的RD灰度图: (a)-(e)分别表示无动态范围转换的无干扰、有宽带干扰、干扰未完全抑制、干扰完全抑制、干扰过度抑制的RD灰度图; (f)-(j)分别表示对应的动态范围为80 dB转换的RD灰度图

    Fig.  3  RD gray-scale image for various interference suppression results: (a)-(e) are images without interference, broadband RFI, incomplete suppression, complete suppression, and excessive suppression, without dynamic range conversion; (f)-(j) are the corresponding images with dynamic range conversion for 80 dB

    图  4  在无干扰的海杂波数据中加入目标回波后, RD图中20个目标的距离-多普勒单元位置

    Fig.  4  The position of 20 targets in the RD map without any RFI

    图  5  干扰抑制自适应优化流程框图

    Fig.  5  Flow chart of adaptive RFI suppression

    图  6  基于不同$ \varepsilon$参数下的RD灰度图与粗糙度

    Fig.  6  RD gray-scale image and coarseness for varying $ \varepsilon$

    图  7  粗糙度$ F$$ \varepsilon$的变化趋势

    Fig.  7  The coarseness $ F$ versus $ \varepsilon$

    图  8  不同$ K$值下$ \varepsilon$变化时的$ k_{best}$比例

    Fig.  8  The ratio of $ k_{best}$ versus $ \varepsilon$ for varying $ K$

    表  1  不同情况下RD图的粗糙度

    Table  1  Coarseness of RD images in various cases

    无动态转换时灰度级有动态转换时灰度级
    25664322566432
    无干扰6.926.926.916.846.846.83
    窄带RFI7.377.377.367.317.317.30
    宽带RFI8.588.588.558.588.588.56
    下载: 导出CSV

    表  2  不同情况下RD图的粗糙度

    Table  2  Coarseness of RD images in various cases

    无动态转换时灰度级有动态转换时灰度级
    25664322566432
    无干扰6.926.926.916.846.846.83
    有干扰8.588.588.558.588.588.56
    未完全抑制7.687.687.667.687.687.67
    完全抑制7.127.127.106.896.896.89
    过度抑制7.687.687.677.677.677.66
    下载: 导出CSV

    表  3  不同目标情况下RD图的粗糙度

    Table  3  Coarseness of RD images of various targets

    目标个数015101520
    无干扰6.846.846.856.866.876.89
    有干扰8.588.588.588.588.588.58
    未完全抑制7.687.687.697.697.697.70
    完全抑制6.896.896.906.916.926.94
    过度抑制7.677.677.697.697.707.71
    下载: 导出CSV

    表  4  滤波器自适应优化算法

    Table  4  Adaptive algorithm of filter optimization

    步骤操作流程
    步1设置滤波器参数$\varepsilon$的优化步长为$\triangle\varepsilon,$阈值比率$\eta\!=\!0.01.$初始化迭代次数$j\!=\!0,$对应滤波器参数为$\varepsilon_0\!=\!1,$计算干扰噪声区域的粗糙度值$F_0.$
    步2$j=j+1,\varepsilon_j=\varepsilon_{j-1}-\triangle\varepsilon,$根据$\varepsilon_j$设计滤波器并抑制干扰, 计算其RD灰度图的粗糙度$F_j.$
    步3计算粗糙度差值$\triangle F=F_{j-1}-F_j.$判断: 若$\triangle F/{F_{j-1}}<\eta,$取得最优抑制, 转至步4; 否则, 转至步2.
    步4$\triangle F<0,\varepsilon_{opt}=\varepsilon_{j-1};$否则, $\varepsilon_{opt}=\varepsilon_{j}.$输出$\varepsilon_{opt}$对应最优滤波器和干扰抑制结果.
    下载: 导出CSV
  • [1] 周文瑜, 焦培南. 超视距雷达技术. 北京:电子工业出版社, 2008. 132−134

    Zhou Wen-Yu, Jiao Pei-Nan. Over-the-horizon Radar Technology. Beijing: Publishing House of Electronics Industry, 2008. 132−134
    [2] Wang W, Wyatt L R. Radio frequency interference cancellation for sea-state remote sensing by high-frequency radar. IET Radar, Sonar and Navigation, 2011, 5(4): 405−415 doi:  10.1049/iet-rsn.2010.0041
    [3] 孔舒亚, 王亮. 压缩感知合成孔径雷达自适应射频干扰抑制方法. 电子测量技术, 2019, 42(13): 107−111

    Kong Shu-Ya, Wang Liang. Adaptive suppression method for SAR narrowband interference based on compressed sensing. Electronic Measurement Technology, 2019, 42(13): 107−111
    [4] 文必洋, 韩金柱, 周企豪, 李艳. 高频地波雷达空域射频干扰抑制算法. 西安电子科技大学学报, 2017, 44(4): 132−137

    Wen Bi-Yang, Han Jin-Zhu, Zhou Qi-Hao, Li Yan. Radio frequency interference suppression algorithm in the spatial domain for high-frequency radars. Journal of Xidian University, 2017, 44(4): 132−137
    [5] 罗忠涛. 新体制天波超视距雷达信号处理研究[Ph. D. dissertation], 电子科技大学, 中国, 2015

    Luo Zhong-Tao. Study on signal processing for future sky-wave over-the-horizon radar system[Ph. D. dissertation], University of Electronic Science and Technology of China, China, 2015
    [6] 朱永恒, 黎明, 牛炯, 张玲, 纪永刚. 基于方向分析的高频地波雷达射频干扰抑制. 中国海洋大学学报, 2018, 48(7): 142−148

    Zhu Yong-Heng, Li Ming, Niu Jiong, Zhang Ling, Ji Yong-Gang. Reasearch on radio frequency interference suppression for HFSWR. Periodical Of Ocean University Of China, 2018, 48(7): 142−148
    [7] 刘建成, 全厚德, 李召瑞, 刘东林, 赵宏志. 多路延迟正交合成的多径信道射频干扰对消. 电子与信息学报, 2017, 39(3): 654−661

    Liu Jian-Cheng, Quan Hou-De, Li Zhao-Rui, Liu Dong-Lin, Zhao Hong-Zhi. RF interference cancellation based on multi-tap delay and orthogonal combination in multipath channel. Journal of Electronics and Information Technology, 2017, 39(3): 654−661
    [8] Luo Z T, He Z S, Li J. AAn effective scheme for radio frequency interference suppression in high-frequency radar. In: Proceedings of IEEE Radar Conference. Arlington, VA, USA: IEEE, 2015. 0539−0544
    [9] Luo Z T, Song T C, He Z S, Hu J F. Approach for transient interference detection based on straight line extraction for high-frequency sky-wave radar. Electronics Letters, 2017, 53(9): 618−620 doi:  10.1049/el.2016.4125
    [10] 罗忠涛, 吴太锋, 何子述, 陈绪元. 基于图像分割的高频雷达射频干扰提取算法. 系统工程与电子技术, 2018, 40(4): 776−781 doi:  10.3969/j.issn.1001-506X.2018.04.10

    Luo Zhong-Tao, Wu Tai-Feng, He Zi-Shu, Chen Xu-Yuan. Extration of radio frequency interference based on image segmentation for high-frequency radar. Systems Engineering And Electronics, 2018, 40(4): 776−781 doi:  10.3969/j.issn.1001-506X.2018.04.10
    [11] 11 Tamura H, Mori S, Yamawaki T. Textural features corresponding to visual perception. IEEE Transactions On Systems, Man, and Cybernetics, 1978, 8(6): 460−473
    [12] Liu Y Q, Huang T Q. Exposing video inter-frame forgery by Zernike opponent chromaticity moments and coarseness analysis. Multimedia Systems, 2017, 23(2): 223−238 doi:  10.1007/s00530-015-0478-1
    [13] 陈苏婷, 胡海锋, 张闯. 基于激光散斑成像的零件表面粗糙度建模. 物理学报, 2015, 64(23): 105−113

    Chen Su-Ting, Hu Hai-Feng, Zhang Chuang. Surface roughness modeling based on laser speckle imaging. Acta Physica Sinica, 2015, 64(23): 105−113
    [14] 王顺杰, 齐春, 程玉胜. Tamura纹理特征在水下目标分类中的应用. 应用声学, 2012, 31(2): 135−139 doi:  10.11684/j.issn.1000-310X.2012.02.013

    Wang Shun-Jie, Qi Chun, Cheng Yu-Sheng. Application of Tamura texture feature to classify underwater targets. Applied Acoustics, 2012, 31(2): 135−139 doi:  10.11684/j.issn.1000-310X.2012.02.013
    [15] Zhang X D, Shen P Y, Gao J R, X X, Qi D, Zhang L, et al.. A license plate recognition system based on tamura texture in complex conditions. In: Proceedings of IEEE International Conference on Information and Automation. Harbin, China: IEEE, 2010. 1947−1952
    [16] Gupta A, Garg M, Mittal A. A comparative performance evaluation of segmented image with obstacle for textural coarseness. In: Proceedings of IEEE International Conference on Computational Intelligence and Computing Research. Coimbatore, India: IEEE, 2010. 1−6
    [17] Majtner T, Svoboda D. Extension of tamura texture features for 3D fluorescence microscopy. In: Proceedings of 2012 Second International Conference on 3D Imaging, Modeling, Processing, Visualization and Transmission. Zurich, Switzerland: IEEE, 2012. 301−307
    [18] Luo Z T, He Z S, Lu K, Chen X Y. Optimal receive filter design under similarity constraint in coloured noise. ET Radar, Sonar and Navigation, 2015, 9(7): 888−899 doi:  10.1049/iet-rsn.2014.0421
  • [1] 聂祥丽, 黄夏渊, 张波, 乔红. 极化合成孔径雷达图像相干斑抑制和分类方法综述[J]. 自动化学报, doi: 10.16383/j.aas.c180097
    [2] 张学峰, 王鹏辉, 冯博, 杜兰, 刘宏伟. 基于多分类器融合的雷达高分辨距离像目标识别与拒判新方法[J]. 自动化学报, doi: 10.3724/SP.J.1004.2014.00348
    [3] 康春玉. 盲源分离与波束形成融合抑制方向性强干扰研究[J]. 自动化学报, doi: 10.3724/SP.J.1004.2014.00983
    [4] 赵涛, 肖建. 基于包含度的区间二型模糊粗糙集[J]. 自动化学报, doi: 10.3724/SP.J.1004.2013.01714
    [5] 蔡秀珊, 高虹, 刘洋. 一类多输入非线性系统的同时H镇定[J]. 自动化学报, doi: 10.3724/SP.J.1004.2012.00473
    [6] 徐怡, 李龙澍. 基于(α, λ)联系度容差关系的变精度粗糙集模型[J]. 自动化学报, doi: 10.3724/SP.J.1004.2011.00303
    [7] 黄超, 何衍, 叶旭东. 基于经典理论的协作跟踪与干扰抑制技术[J]. 自动化学报, doi: 10.3724/SP.J.1004.2011.00766
    [8] 祝依龙, 范红旗, 卢再奇, 付强. 基于高分辨一维多普勒像的雷达目标机动检测算法[J]. 自动化学报, doi: 10.3724/SP.J.1004.2011.00901
    [9] 蔡秀珊. 鲁棒镇定受到干扰的非线性微分包含系统[J]. 自动化学报, doi: 10.3724/SP.J.1004.2010.01327
    [10] 金小峥, 杨光红. 考虑有执行器故障和有界扰动的鲁棒自适应容错补偿控制[J]. 自动化学报, doi: 10.3724/SP.J.1004.2009.00305
    [11] 陈刚, 王树青. 约束环境中多指手操作的鲁棒控制[J]. 自动化学报
    [12] 刘翔, 王文海, 熊斌, 孙优贤. 典型大时变时滞工业过程的鲁棒数字PⅡ调节控制[J]. 自动化学报
    [13] 张卫东, 许晓鸣, 孙优贤. 一类不稳定时滞过程的最优抗干扰控制[J]. 自动化学报
    [14] 吴方向, 史忠科, 戴冠中. 区间系统的H∞鲁棒控制[J]. 自动化学报
    [15] 王晓军, 袁梅, 吴立德. 基于多视角距离图象的三维物体建模及其在识别中的应用[J]. 自动化学报
    [16] 吴俊, 方华京. 差分有界干扰的最优抑制[J]. 自动化学报
    [17] 袁璞, 左信, 郑海涛. 状态反馈预估控制[J]. 自动化学报
    [18] 丁原志, 宋晓梅. 雷达跟踪系统的随机干扰自适应[J]. 自动化学报
    [19] 刘其真. 用锁相式采样法抑制周期性干扰噪声[J]. 自动化学报
    [20] 王世林, 许可康, 韩京清. 干扰解耦问题[J]. 自动化学报
  • 加载中
计量
  • 文章访问数:  4
  • HTML全文浏览量:  2
  • 被引次数: 0
出版历程

OTH雷达图像的粗糙度指标及用于射频干扰自适应抑制

doi: 10.16383/j.aas.c190286
    基金项目:  国家自然科学基金 (61701067, 61702065), 重庆市教育委员会科研基金 (KJ1600427, KJ1600429)资助
    作者简介:

    重庆邮电大学通信与信息工程学院讲师, 主要研究方向为信号与信息处理, 统计信号处理, 数字图像处理.E-mail: luozt@cqupt.edu.cn

    重庆邮电大学通信与信息工程学院硕士研究生, 主要研究方向为信号与信息处理, 数字图像处理.E-mail: s180131036@stu.cqupt.edu.cn

    重庆邮电大学通信与信息工程学院学士, 主要研究方向为数字图像处理.E-mail: guojie1072050774@163.com

    电子科技大学电子工程学院教授, 主要研究方向为阵列信号处理, 自适应信号处理, MIMO雷达与通信.E-mail: zshe@uestc.edu.cn

    中国电子科技集团有限公司第十四研究所研究员级高级工程师, 主要研究方向为超视距雷达系统设计和信息处理.E-mail: mimimomoba@gmail.com

摘要: 针对OTH雷达距离-多普勒(Range-Doppler, RD)图, 本文首次提出采用纹理粗糙度作为RD图质量的评价指标, 即计算RD图所转化灰度图的Tamura纹理粗糙度. 分析表明, 粗糙度指标能准确反映RD图受干扰情况, 对于不同灰度转换函数具有稳健性. 作为应用举例, 本文将图像粗糙度用于改进射频干扰抑制算法, 使干扰抑制达到自适应优化. 实验结果表明, Tamura粗糙度能够正确反映RD图干扰抑制情况, 优化粗糙度指标能够使干扰抑制自适应达到最优.

English Abstract

罗忠涛, 郭人铭, 郭杰, 何子述, 卢琨. OTH雷达图像的粗糙度指标及用于射频干扰自适应抑制. 自动化学报, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286
引用本文: 罗忠涛, 郭人铭, 郭杰, 何子述, 卢琨. OTH雷达图像的粗糙度指标及用于射频干扰自适应抑制. 自动化学报, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286
Luo Zhong-Tao, Guo Ren-Ming, Guo Jie, He Zi-Shu, Lu Kun. Tamura coarseness for OTH radar image evaluation and its application in adaptive optimization of interference suppression. Acta Automatica Sinica, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286
Citation: Luo Zhong-Tao, Guo Ren-Ming, Guo Jie, He Zi-Shu, Lu Kun. Tamura coarseness for OTH radar image evaluation and its application in adaptive optimization of interference suppression. Acta Automatica Sinica, 2020, 46(x): 1−9. doi: 10.16383/j.aas.c190286
  • 现代雷达普遍采用动目标检测技术, 信号检测通常在距离-多普勒(range-Doppler, RD)图中进行. 当雷达信号带宽内出现干扰时, 干扰会叠加到RD单元幅度上, 从而影响目标检测. 该现象在超视距(Over-the-Horizon, OTH)雷达中经常遇到[1].

    OTH雷达工作在高频段, 易受其它设备信号干扰的影响. 射频干扰被雷达接收并经信号处理后, 会在RD图上表现出不同形态, 影响目标检测甚至是海杂波谱探测[1, 2]. 干扰抑制一直受到高频雷达研究的关注, 学者们也提出了基于不同原理的干扰抑制方法[3-8]. 不过, 由于天波雷达中干扰和杂波特性均为未知, 传统干扰抑制方法都不能从理论上证明某方法最优, 也难以提出明确的优化方向.

    事实上, 基于分析RD图的经验发现, 尽管从数据特征上难以评判RD图的优劣, 但用人眼观察RD图却很容易有直观的评判. 如图1所示, 噪声区域像背景底色; 海杂波区域是不规则边缘的带状区域; 窄带干扰表现为平行于距离维的高亮直线; 宽带干扰表现为遍布全图的条纹状[5]. 当干扰抑制不充分或过度抑制时, RD图形态也会变化. 尽管此类变化从传统的信号处理角度来说难以表征和量化, 从人眼和图像的角度却可以直观发现.

    图  1  OTH雷达RD图

    Figure 1.  RD maps of OTH radar

    本文拟从图像分析的角度, 研究OTH雷达的RD图像质量及干扰问题. 之前, 作者研究直线检测技术在OTH雷达中的应用, 如在文献[9]中提出快-慢时间图中检测瞬态干扰直线, 在文献[10]中提出了RD灰度图中检测窄带干扰直线. 本文考虑从图像分析角度, 分析RD图像的特征与形态, 通过提取图像的纹理特征来判断RD图的质量.

    纹理是图像的一个重要特征, 可用于分析许多类型的图像, 如自然图像、遥感图像、医疗图像等. 作为一种视觉特征, 纹理分析的方法有多种, 如空间自相关法、共生矩阵法、Tamura方法等. 其中, Tamura纹理由Hideyuki Tamura于1978年提出[11], 它基于人类对纹理的视觉感知, 侧重对图像全局特征的提取, 其特征量中对应于心理学角度上纹理特征的6种属性, 即粗糙度、对比度、方向度、线性度、规整度和粗略度.

    在Tamura纹理中, 粗糙度(Coarseness)是最基本的纹理特征, 可以直观反映人眼对图像粗糙程度的感知, 常用于基于内容的图像检索与评估、异常纹理检测等. 例如, 文献[12]通过提取视频帧Tamura纹理特征来检测数字视频是否被复制移动或篡改; 文献[13]提出基于金属零件表面图像Tamura纹理的粗糙度、对比度和方向度来识别不同工艺加工的零件; 文献[14]提出基于水下目标功率谱图的Tamura纹理的粗糙度和对比度, 检测与识别水下目标速度的高低. 文献[15]提出使用Tamura纹理特征实现对嘈杂环境和不良照明下车牌的准确定位. 文献[16]提出基于纹理的障碍物检测算法, Tamura方法比空间自相关法和共生矩阵法更加稳健. 在生物医学图像分析中, Tamura纹理可用于表征研究对象的形状和内部结构, 并判断不同形态器官组织代表的信息[17]. 总之, Tamura粗糙度在图像分析中有着广泛应用.

    据作者所知, 目前尚无文献对Tamura纹理和RD图像结合展开研究. 本文将开展创新性工作, 首次讨论将Tamura纹理粗糙度用于OTH雷达RD图像, 并分析其合理性. 其次, 考虑RD图像可能受干扰的情况, 分析粗糙度如何反映受干扰情况. 最后给出粗糙度的使用举例, 以粗糙度作为指标实现射频干扰的自适应抑制.

    • 本小节介绍OTH雷达的RD图处理方法, 利用灰度变换函数将其转化为灰度图, 再给出Tamura粗糙度的计算方法.

    • 设OTH雷达接收阵列为$ K $个阵元组成的均匀线阵, 一个相干处理间隔内有$ N $个采样点, 雷达接收数据排列为矩阵$ {{X}}_{R}\in C^{K\times N}. $采用$ K $维权向量$ {{a}} $$ {{X}}_{R} $进行波束形成, 输出为时域数据$ {{x}}_{R} = {{a}}\cdot{{{X}}_{R}} .$对时域数据$ {{x}}_{R} $进行脉冲压缩和多普勒处理, 即可得到RD矩阵$ {{Z}}\in{C^{P\times M}}. $其中, $ P $为一个相干处理间隔内的积累周期数, $ M $为单周期内的采样点数. 将RD矩阵$ {{Z}} $的值按照幅度dB值显示, 即可得到RD图[5].

      考虑获得RD图后进入图像分析阶段. 首先将RD图转化为灰度图. 人们在观察时, 一般采用dB值来呈现RD图, 以便显示绝对幅度比海杂波小得多的噪声和目标.

      本文采用基于对数函数的灰度变换方法, 将常规RD图转化为RD灰度图. 对数变换中可设置合理的最小值, 以控制原数据的动态范围. 例如, 设置动态范围为$ b $ dB, 则RD矩阵元素不会小于最小值

      $$ {\left| {{{Z}}} \right|_{\min }} = {\left| {{{Z}}} \right|_{\max }} \times {10^{ - b/20}} $$ (1)

      其中, $ \left| {{{Z}}} \right|_{\max } $表示矩阵$ {{Z}} $中的最大模值. 总的灰度变换函数可以表示为

      $$ {{Y}}(p,m) = c \cdot \log\{\max[|{{Z}}(p,m)|,|{{Z}}|_{\min}]/{|{{Z}}|_{\min}}\} $$ (2)

      其中$ c $控制灰度级, 例如256或64.

      经过灰度变换, 原RD矩阵$ {{Z}} $或RD图转化为RD灰度图, 原距离-多普勒单元$ (p,m) $的幅度转换成灰度图中像素点$ (p,m) $的灰度值.

    • 粗糙度是反映纹理中粒度的一个量, 是最基本的纹理特征. 粗糙度描述了一幅图像平均的粗糙程度. 针对灰度图$ {{Y}}(p,m), $粗糙度的计算可分为3步.

      (i) 计算图像中大小为$ 2^k\times{2^k} $个像素的活动窗口中像素的平均强度值, 即:

      $$ {{{{A}}}_k}(x,y) = \sum\limits_{p = x - {2^{k - 1}}}^{x + {2^{k - 1}} - 1} {\sum\limits_{m = y - {2^{k - 1}}}^{y + {2^{k - 1}} - 1} {{{{\bf{Y}}(p,m)} / {{2^{2k}}}}} } $$ (3)

      其中: $k = 1,\cdots,K$控制活动窗口大小, 一般$ K $值在2至6之间; $ {{Y}}(p,m) $为像素点$ (p,m) $的灰度值.

      (ii) 计算每个像素点在水平和垂直方向上互不重叠的窗口之间的平均强度差

      $$ {{{E}}_{k,h}}(x,y) = |{{{{A}}}_k}(x+2^{k-1},y)-{{{{A}}}_k}(x-2^{k-1},y)| $$ (4)
      $$ {{{E}}_{k,v}}(x,y) = |{{{{A}}}_k}(x,y+2^{k-1})-{{{{A}}}_k}(x,y-2^{k-1})| $$ (5)

      其中$ {{{E}}_{k,h}}(x,y) $代表水平方向差值, $ {{{E}}_{k,v}}(x,y) $代表垂直方向差值. 对于每个像素, 假设$ k_{best} $是使$ {{{E}}_{k,h}}(x,y) $$ {{{E}}_{k,v}}(x,y) $中较大者达到最大的$ k $值, 对应最佳尺寸

      $$ {{S}_{best}}(p,m) = 2^{k_{best}} $$ (6)

      (iii) 粗糙度计算为整幅图最佳尺寸的平均值

      $$ F = \frac{1}{{P \times M}}\sum\limits_{p = 1}^P {\sum\limits_{m = 1}^M {{S_{best}}(p,m)} } $$ (7)

      基于OTH雷达RD图转化的灰度图, 经过(i)-(iii)三步运算, 即可得到该图的粗糙度. 考虑到一般情况下, OTH雷达重点关注目标检测, 而RD图中非海杂波即噪声区域才是潜在的目标检测区域. 因此, 可采用RD灰度图的海杂波区域以外(例如±20倍Bragg频率外[5])区域来计算粗糙度.

    • 本文提出将粗糙度作为一个指标, 来衡量RD灰度图质量. 那么, 需要讨论粗糙度的影响因素和应具有的性质. 首先, 粗糙度应当稳健, 因为各相干处理间隔的数据功率有所不同, 生成灰度图时灰度转换函数含有人为参数. 其次, 粗糙度需要能够表征RD图质量, 即准确反映RD图是否受干扰影响、完全抑制或过度抑制等情况.

    • 首先分析灰度变换对粗糙度指标的影响. 将RD图经灰度变换函数转换为灰度图像的过程中, 动态范围转换和灰度级都会影响RD灰度图, 从而影响灰度图的粗糙度.

      仿真一组相干处理间隔为5.12秒的海杂波数据, 再结合无干扰、宽带干扰和窄带干扰数据, 分别处理得到RD图, 再由灰度变换生成RD灰度图, 设置灰度级为256. 当未设置动态范围时, 不同干扰情况的RD灰度图为图2(a)(b)(c); 设置动态范围$ b = 80 $ dB时, 对应的灰度图分别是图2(d)(e)(f). 可以看到, 设置动态范围后, 大部分像素灰度值下降, 海杂波、目标与背景噪声的对比更为明显, 更利于视觉观察.

      图  2  不同干扰情形下的RD灰度图: (a)-(c)分别表示无动态范围转换的无干扰、有窄带射频干扰、有宽带射频干扰的RD灰度图; (d)-(f)分别表示对应的动态范围为80 dB转换的RD灰度图

      Figure 2.  RD gray-scale images: (a)-(c) are images without dynamic range conversion, for no RFI, narrowband RFI, and wideband RFI, respectively; (d)-(f) are corresponding images with dynamic range conversion for 80 dB

      横向比较图2中的三幅图. 在图1(a)(d)中, RD图无干扰, 噪声区域的幅度分布均匀, 是理想的噪声形态. 目标清晰显示, 可认为是RD图像质量最佳. 在图2(b)(e)中, RD图存在窄带射频干扰直线, 目标检测也会受影响, 图像质量较差. 在图2(c)(f)中, 宽带射频干扰遍布RD图, 全图区域布满条纹, 目标难以观察, 可认为图像质量很差.

      针对图2中各图数据, 采用$ K = 4 $和灰度级256、64和32, 粗糙度计算结果如表1.

      表 1  不同情况下RD图的粗糙度

      Table 1.  Coarseness of RD images in various cases

      无动态转换时灰度级有动态转换时灰度级
      25664322566432
      无干扰6.926.926.916.846.846.83
      窄带RFI7.377.377.367.317.317.30
      宽带RFI8.588.588.558.588.588.56

      表1进行横向分析: 对于每种干扰状态的粗糙度, 在不同灰度级256或64或32下, 有动态范围比无动态范围时稍小一点, 可认为是基本相同. 说明在不同灰度图转换方式下, 每种干扰状态的粗糙度是稳健的.

      表1进行纵向分析: 不含干扰时粗糙度最低, 含窄带射频干扰时次之, 宽带射频干扰的粗糙度最高. 从是否有利于目标检测的图像质量观点来看, 可发现一个规律: 粗糙度越低, 图像质量越高, 噪声区域越干净, 越利于目标检测; 反之, 粗糙度越高, 图像质量越差, 越不利于目标检测.

    • 下面分析干扰情况对粗糙度的影响, 观察粗糙度是否能够准确表现干扰严重程度和抑制程度. 当OTH雷达RD图中含干扰时, 干扰严重程度及是否充分抑制, 会在RD灰度图上表现出不同的形态, 而不同的形态所对应的粗糙度值也不同. 图3仿真不同干扰情况和抑制情况, 分别为无干扰、有宽带射频干扰、及干扰抑制程度不同的RD灰度图. 图3中, $ (a)-(e) $表示灰度转换时无动态范围, $ (f)-(j) $表示灰度转换有动态范围80 dB.干扰抑制算法采用了文献[18]中的滤波器设计方法, 未完全抑制或完全抑制效果是由过大或过小的相似度值产生.

      图  3  不同干扰抑制结果的RD灰度图: (a)-(e)分别表示无动态范围转换的无干扰、有宽带干扰、干扰未完全抑制、干扰完全抑制、干扰过度抑制的RD灰度图; (f)-(j)分别表示对应的动态范围为80 dB转换的RD灰度图

      Figure 3.  RD gray-scale image for various interference suppression results: (a)-(e) are images without interference, broadband RFI, incomplete suppression, complete suppression, and excessive suppression, without dynamic range conversion; (f)-(j) are the corresponding images with dynamic range conversion for 80 dB

      可以观察到, 当图3(a)(f)中RD图无干扰时, 噪声区域的幅度分布较均匀, 目标可观察. 当图3(b)(g)中RD图含宽带射频干扰时, 目标难以观察. 图3(a)(f)中经滤波器设计抑制干扰后, 图3(c)(h)是干扰未完全抑制的RD灰度图, 可见噪声区域的纹理粗糙. 图3(d)(i)是干扰完全抑制的RD灰度图, 可见噪声区域幅度分布较均匀, 目标可观察. 图3(e)(j)是干扰过度抑制的RD灰度图, 可见噪声区域纹理较粗糙, 目标被掩盖.

      图3分别计算当灰度级为256、64和32时的粗糙度, 计算结果如表2所示. 从表2可以看到, 无干扰时, 粗糙度数值处于最低; 有干扰但被完全抑制时, 粗糙度次之; 有干扰且未抑制时, 粗糙度最高; 当干扰未完全抑制或者过度抑制时, 噪声区域的粗糙度介于最优抑制和未抑制之间. 此外, 有无动态范围转换或灰度级大小, 对粗糙度基本无影响.

      表 2  不同情况下RD图的粗糙度

      Table 2.  Coarseness of RD images in various cases

      无动态转换时灰度级有动态转换时灰度级
      25664322566432
      无干扰6.926.926.916.846.846.83
      有干扰8.588.588.558.588.588.56
      未完全抑制7.687.687.667.687.687.67
      完全抑制7.127.127.106.896.896.89
      过度抑制7.687.687.677.677.677.66
    • 接下来分析目标数量对粗糙度的影响. 模拟不同径向速度和距离的目标回波, 与海杂波一起组成雷达回波. 在无干扰的情况下, 采用传统脉冲压缩和多普勒处理, 输出RD图如图4. 总体上目标位置呈现出随机分布, 每个目标形成了一个光斑, 占据了多个距离-多普勒单元, 可视为反映了目标回波在电离层等影响下距离和多普勒维的旁瓣展宽. 我们想知道, 目标像素点会在多大程度上影响粗糙度.

      图  4  在无干扰的海杂波数据中加入目标回波后, RD图中20个目标的距离-多普勒单元位置

      Figure 4.  The position of 20 targets in the RD map without any RFI

      测量不同目标情况对应的粗糙度, 观察目标数量对粗糙度的影响. 考虑图3中的干扰有无和抑制程度的不同情况, 再加上目标数量变化的因素, 组成不同干扰、目标和抑制的多种组合, 处理得到各自的RD灰度图, 并计算其粗糙度.

      分别以目标数量为0, 1, 5, 10, 15, 20为例, 考虑图3(f)-(j)中干扰和抑制处理情况, 得到256灰度级动态范围为$ 80 $ dB的RD灰度图, 计算粗糙度的结果如表3.

      表 3  不同目标情况下RD图的粗糙度

      Table 3.  Coarseness of RD images of various targets

      目标个数015101520
      无干扰6.846.846.856.866.876.89
      有干扰8.588.588.588.588.588.58
      未完全抑制7.687.687.697.697.697.70
      完全抑制6.896.896.906.916.926.94
      过度抑制7.677.677.697.697.707.71

      表3中, 当存在干扰但无抑制或未充分抑制时, 目标对粗糙度基本无影响. 这并不难理解, 目标幅度往往较干扰更低故而被干扰掩盖掉, 因此对整图的粗糙度未形成明显影响.

      当干扰不存在, 以及干扰存在且完全或过度抑制时, 目标数量对RD图的粗糙度有了一定的影响. 由表3可以看到, 当目标数量由0增至20, 无干扰下的粗糙度由6.84增长到6.89; 干扰完全抑制下的粗糙度由6.89增长至6.94. 总体而言, 随目标数据增加, 粗糙度的增长幅度很小, 不到粗糙度值的1%.

      最后从纵向看表3. 在相同目标个数下, 无干扰情况具有最小的粗糙度; 有干扰但被完全抑制时, 粗糙度次之; 有干扰且未抑制时, 粗糙度最高. 这些仿真结果与2.2节中无目标时是一致的. 此外, 改变动态范围或灰度级大小, 对于目标存在下的粗糙度结果基本无影响. 因此, 可以说目标对于粗糙度的影响很小.

    • 除以上仿真外, 作者基于实测的OTH雷达数据, 进行了RD灰度图转换和纹理粗糙度分析, 其结果与上面的仿真结果基本吻合. 结合已有的研究, OTH雷达图像粗糙度分析的结论主要有以下三点.

      (1)一组OTH雷达相干处理间隔内数据, 生成RD灰度图并计算Tamura纹理粗糙度, 该粗糙度数值在灰度转换函数参数时基本不变. 这说明, 针对不同类型和状态下的OTH雷达数据, 粗糙度可以作为一个稳健的衡量指标.

      (2)以目标检测环境是否理想的视角来看, 纹理粗糙度能够正确反映RD灰度图的理想程度, 不受干扰时的理想RD图粗糙度小, 受干扰时的RD图粗糙度数值大, 最优抑制的RD图粗糙度与无干扰时粗糙度很接近.

      (3)纹理粗糙度作为衡量指标, 基本不受目标存在与否及其数量的影响. 目标的存在既不会明显改变粗糙度, 也没有改变干扰抑制情况与粗糙度大小的对应次序. 因此, 粗糙度可以正确地衡量有目标时的RD图理想程度.

      综上, 纹理粗糙度表现稳健, 能够正确反映RD灰度图是否理想, 可作为RD图质量的指标. 以粗糙度为指标, OTH雷达不仅可以评估雷达工作状态, 还可帮助提升信号处理算法, 实现自适应优化处理.

    • 本文以基于滤波器设计的射频干扰抑制为例, 展示粗糙度指标如何使干扰抑制变得自适应. 以此例说明粗糙度在OTH雷达信号处理中具有很大的应用潜力.

      针对射频干扰, 传统抑制算法虽多, 但由于干扰、噪声及杂波参数未知, 没有理论指标作参考, 不能确认是否取得自适应最优效果. 在引入粗糙度指标后, 以实现噪声区域(非杂波区域)的指标最小化为目标, 可对射频干扰进行自适应抑制处理, 避免干扰抑制不充分或过度抑制的情况.

    • 粗糙度指标能够自适应地判断和评价干扰情况. 以文献[18]中固定相似度约束下设计滤波器的干扰抑制方法为例, 基于粗糙度指标优化的自适应干扰抑制方案如图5. 其中, 对RD图粗糙度的计算与优化, 构成了一个循环迭代过程. 如图5中阴影部分所示, 首先将传统的RD图转化为灰度图像, 然后计算Tamura纹理粗糙度. 当粗糙度指标未达到最小时, 该方案会继续优化滤波器参数, 循环迭代直至粗糙度达到最小.

      图  5  干扰抑制自适应优化流程框图

      Figure 5.  Flow chart of adaptive RFI suppression

      借助于图像粗糙度, 自适应干扰抑制算法具体流程如表4. 其中, 干扰抑制算法的参数是相似度$ \varepsilon, $取值范围为$ (0,1]. $本文采用一个简单的迭代方法: 设置初始值$ \varepsilon_0 = 1 ,$步长$ \triangle\varepsilon $为一个较小值, 如0.01; 令$ \varepsilon $步进减小, 表示抑制程度在增加, 直到取得最小粗糙度值. 相比传统干扰抑制方法[8]设定固定相似度0.8, 本文方法以取得最优相似度$ \varepsilon_{opt}, $对应于最小粗糙度值, 实现干扰抑制最优化.

      表 4  滤波器自适应优化算法

      Table 4.  Adaptive algorithm of filter optimization

      步骤操作流程
      步1设置滤波器参数$\varepsilon$的优化步长为$\triangle\varepsilon,$阈值比率$\eta\!=\!0.01.$初始化迭代次数$j\!=\!0,$对应滤波器参数为$\varepsilon_0\!=\!1,$计算干扰噪声区域的粗糙度值$F_0.$
      步2$j=j+1,\varepsilon_j=\varepsilon_{j-1}-\triangle\varepsilon,$根据$\varepsilon_j$设计滤波器并抑制干扰, 计算其RD灰度图的粗糙度$F_j.$
      步3计算粗糙度差值$\triangle F=F_{j-1}-F_j.$判断: 若$\triangle F/{F_{j-1}}<\eta,$取得最优抑制, 转至步4; 否则, 转至步2.
      步4若$\triangle F<0,\varepsilon_{opt}=\varepsilon_{j-1};$否则, $\varepsilon_{opt}=\varepsilon_{j}.$输出$\varepsilon_{opt}$对应最优滤波器和干扰抑制结果.
    • 采用MATLAB软件仿真和处理天波雷达数据, 验证本文提出的灰度纹理特征判断图像质量和自适应干扰抑制方法. 设雷达发射线性调频连续波, 波形参数为带宽$ 20k $ Hz, 采样频率$ 50k $ Hz, 波形周期$ 10\;ms $, 相干积累周期数512. 模拟射频干扰, 与目标回波及海杂波相加为雷达接收数据, 采用脉冲压缩和快速傅里叶变换做多普勒处理, 设置动态范围为$ 80 $ dB, 所得RD图如图6(a)所示. 图中每个像素点即距离-多普勒单元的多普勒维长度约为$ 0.2 $ Hz, 距离维长度约为$ 3\;km. $注意, 图1图2图3中RD图均基于此组参数产生.

      图  6  基于不同$ \varepsilon$参数下的RD灰度图与粗糙度

      Figure 6.  RD gray-scale image and coarseness for varying $ \varepsilon$

      图6(a)可见, 宽带射频干扰存在时图像纹理粗糙, 噪声区域的纹理特性不均匀. 采用本文方法计算噪声区域的粗糙度. 首先, 以动态范围$ 80 $ dB和灰度级256, 将RD图转换为灰度图. 然后, 取20倍Bragg频率以外的区域为噪声区域, 计算噪声区域粗糙度为$ F_0 = 8.58 $.

      采用表4中自适应干扰抑制算法, 设置粗糙度计算参数$ K = 4 $, 相似度初始值$ \varepsilon_0 = 1 $, 优化步长$ \triangle\varepsilon = 0.01 $, 对干扰抑制进行自适应优化. 图7中用正方形标记的虚线给出了自适应算法的迭代过程. 当迭代次数$ j $由0增加到5时, $ \varepsilon $由1减小到0.95, 粗糙度在逐步减小, 对应干扰抑制效果也在优化; 当$ j = 6 $$ \varepsilon = 0.94 $时, 粗糙度增大, 说明进入过度抑制阶段. 因此迭代终止, 最优抑制停在$ \varepsilon = 0.95 $, 此时粗糙度值为6.89, 基本等于无干扰时粗糙度6.84. 比较图6(c)图2(d), 可以看到, 本文最优抑制结果与无干扰时RD图像非常接近.

      图  7  粗糙度$ F$$ \varepsilon$的变化趋势

      Figure 7.  The coarseness $ F$ versus $ \varepsilon$

      仿真其它相似度$ \varepsilon $值对应的粗糙度, 如图7中实线所示. 可以看到, 粗糙度的变化趋势是先减后增, 表明粗糙度确实存在最小值且容易寻找. 四个典型相似度参数$ \varepsilon $值所对应的RD灰度图如图6, 其中图6(d)为白化滤波器即$ \varepsilon = 0.76 $的结果. 可以看到, 粗糙度确实能够正确反映人们对RD灰度图质量的视觉效果. 此外, 改变$ K $值为3和5, 粗糙度与相似度的对应关系保持不变, 如图7所示.

      为观察粗糙度更多细节, 图8画出了$ K $分别为3,4,5时最佳尺寸参数$ k_{best} $的比例随$ \varepsilon $的变化. 三幅图具有共同特点: $ k_{best} = 1 $比例总是随干扰抑制优化($ \varepsilon $从1到0.95)而增大, 随抑制恶化($ \varepsilon $从0.95到0.94)而减小. 当$ K = 5 $时, $ k_{best} = 4,5 $比例的增减趋势与$ k_{best} = 1,2 $恰好相反; 当$ K = 4 $时, $ k_{best} = 4 $比例与$ k_{best} = 1 $恰好相反. 式$(6)(7)$可见, $ k_{best} $越小则粗糙度越小, 故小$ k_{best} $值的比例与粗糙度增减趋势正好相反. 这与仿真图相吻合: 针对同一$ K $值, 图7$ k_{best} = 1 $比例越大, 对应图7中粗糙度越小. 结合图7图8可见, $ K = 4 $和5对应的粗糙度变化趋势比$ K = 3 $更加明显, 故更适合用来计算粗糙度.

      图  8  不同$ K$值下$ \varepsilon$变化时的$ k_{best}$比例

      Figure 8.  The ratio of $ k_{best}$ versus $ \varepsilon$ for varying $ K$

      此外, 针对多目标场景的仿真结果表明, 目标对粗糙度影响很小, 不改变粗糙度与相似度的变化趋势, 本文干扰抑制方法依然可行. 另外, 本文算法也用于实测OTH雷达数据的干扰抑制, 处理结果与本文仿真相似. 这说明, 本文提出的粗糙度指标和自适应滤波器设计算法, 可以有效抑制射频干扰, 达到最优抑制效果. 对于本文抑制方案, 未来工作可考虑运算量上的优化, 例如缩小粗糙度计算的图像区域, 或采用线搜索算法来迭代优化参数.

    • 针对OTH雷达图像, 本文首次提出计算RD图中非杂波区域的Tamura纹理粗糙度, 作为RD图像质量的评估指标. 粗糙度指标能够反映干扰抑制的不同程度, 为干扰抑制提供了合理的优化方向. 结合RD灰度图的Tamura纹理粗糙度, 能够设计出具有自适应优化能力的干扰抑制算法. 实验结果表明, 借助于Tamura纹理粗糙度优化, 自适应干扰抑制算法能够取得最佳干扰抑制效果. 作为OTH雷达RD图像质量评估指标, Tamura粗糙度有着很好的应用潜力, 将在未来研究中发展出更多应用.

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回