一种电磁干扰源成像多分辨率分区算法.pdfVIP

第27卷第10期 强 激 光 与 粒 子 束 V01．27，No．10 BEAMS Oct．。2015 2015年10月 HIGHPOWERLASERANDPARTICLE 一种电磁干扰源成像多分辨率分区算法‘ 杜 鑫， 谢树果， 郝旭春， 王 超 (北京航空航天大学电子信息工程学院电磁兼容技术研究所，北京 100191) 摘要： 在利用抛物反射面对电磁干扰源成像过程中，由于系统衍射受限导致干扰源成像模糊，分辨率 低，难以分辨，由于不同频率不同区域干扰源所成图像分辨率不同，具有分区域多分辨率的特征，采用已有超分 辨算法难以提高分辨率。利用Meanshift算法，在原有算法基础上改进使其能够适应多分辨率的电磁干扰源 成像，在图像分割的基础上对多分辨率图像进行分块抽离，并采用基于L—R迭代的盲反卷积算法分别对各区 域进行分辨率的提高，仿真结果表明算法能够适应对干扰源的多分辨率电磁成像并提高分辨率。 关键词：多分辨率； 电磁成像；MeanShift算法；盲反卷积 中图分类号： TP391 文献标志码： A doi：10．11884／HPLPB201527．103223 目前对电磁干扰源的检测手段存在很多缺点，如检测速度慢、难以全面检测等，而利用抛物反射面进行电 磁成像的方法可以快速准确地检测干扰源的位置。电磁成像研究目前在高频段有太赫兹成像口屯j，在较低频段 主要集中在被动毫米波成像领域，受天线尺寸、衍射受限等因素的影响，所成图像分辨率不高，相关的提高分辨 率的超分辨算法成为了研究的热点[3。5]。而对于电磁干扰源的检测成像，与毫米波成像相比波长更长，衍射受 限更为严重，且所成图像中包含多种干扰源的点扩展函数，如何提高其分辨率使得干扰源能够清晰分辨成了新 的难点。抛物反射面对电磁干扰源的成像难以分辨是由衍射受限引起的，按照傅里叶光学的观点，抛物反射面 相当于低通滤波器，截止频率以上的高频信息“丢失”，对于此类问题需要利用先验知识对高频信息复原，称为 求成像系统需满足奈奎斯特采样条件，无法处理欠采样系统。近年来超分辨算法基于成像帧数可以分为单帧 超分辨算法与多帧超分辨算法，而基于运算方式可以分为频域算法及空域的算法，空域算法又可以细分为非均 匀差值法、统计论方法、凸集投影法、正则化法、盲重建法等。这些超分辨算法均只能对包含一个点扩展函数 (PSF)的图像进行处理，而对于衍射受限的抛物反射面系统，所成图像中干扰源的点扩展函数由于频率、位置 等因素导致多种多样，现有算法难以满足需求。本文利用MeanShift算法对不同分辨率干扰源进行提取，并 结合基于L—R迭代的盲反卷积算法，对多分辨率的电磁干扰源成像进行分辨率的提高。 1 算法理论 对于抛物反射面衍射受限系统所成干扰源图像，由于各干扰源分辨率不同，即各个源在衍射系统中所对应 的点扩展函数不同，简单采用现有超分辨算法难以完成分辨率的提高。但干扰源成像存在相近干扰源点扩展 函数相近的特点，即同一区域内干扰源的点扩展函数基本相同。如果我们可以将各区域内干扰源抽离出来，分 别进行提高分辨率的运算，再将运算后的各区域干扰源分别复原，则可以完成对多分辨率电磁成像的分辨率提 高。 区域抽离的基础是进行图像的分割。简单的图像分割方法包括均匀分割与非均匀分割均。采用均匀分割 的方法即对所成图像横向分割或纵向分割，规定所分区域的个数，将图像均匀的分割，这种方法显然是不适用 的，无法准确找到各干扰源位置，并且易将一个干扰源分到不同区域。非均匀分割法则根据各干扰源位置和大 小进行分割，可以准确找到分辨率相同的干扰源所在区域，但需要手动进行，耗时长且没有统一的分割方法。 作为图像处理领域的热点问题，近年来已有许多有效的图像分割方法，而其中基于MeanShift即均值漂移的 方法适用于对电磁辐射源成像的分割。 Mean Shift算法在一篇基于概率密度函数的梯度估计文献中被提出[8]，之后经过改进扩大了其使用范 围凹]，并被有效地利用到了图像分割、图像平滑等方面并取得了不错的效果[1”13]。其基本思想为假设一个d *收稿日期：2015一06—1l；修订日期：2015一06—30 基金项