一种适用于穿墙雷达多径幻象抑制的成像方法
【专利摘要】本发明提供一种适用于穿墙雷达多径幻象抑制的成像方法，包括子孔径成像步骤：将全孔径分为N个没有交集的子孔径；对每个子孔径，分别通过后向投影成像算法，得到N个独立的子图像；子图像融合步骤：对N个独立的子图像分别进行归一化处理，得到N个归一化子图像I′i；对于N幅归一化的子图像进行联乘得到融合后的图像If；最后对融合后的图像If进行门限检测与二值化处理得到最终的雷达图像I。本发明有效抑制了多径幻象干扰，提高了图像的信杂比和墙后隐蔽目标的探测性能。
【专利说明】
-种适用于穿墙雷达多径幻象抑制的成像方法
技术领域
[0001] 本发明属于穿墙雷达成像技术，设及多径幻象抑制技术。
【背景技术】
[0002] 穿墙雷达成像技术主要利用天线阵元发射特定电磁波穿透建筑物墙体，实现建筑 物布局成像和墙后隐蔽目标的检测、成像与定位，近年来在反恐、执法、灾难救援等领域具 有广泛应用和重大的社会研究价值。然而，在封闭建筑物探测环境中，电磁波在建筑物墙体 内表面产生多次反射，引入严重的多径杂波干扰，导致成像结果中产生大量多径幻象，影响 真实目标的检测成像性能。因此，多径幻象的定位与抑制是实现建筑物内隐蔽目标检测成 像的关键技术。
[0003] 对于穿墙雷达中的多径幻象干扰，国内外科研机构开展了很多研究。美国维拉诺 瓦大学提出两种多径幻象的抑制方法。第一种方法假设电磁波在墙体表面发生镜面反射， 计算成像结果中多径幻象出现的区域，然后利用2-D高斯权函数，将多径幻象区域的像素点 搬移并叠加到相应目标位置区域，从而实现多径幻象干扰的抑制，但该方法会造成多径幻 象区域真实目标的丢失，出现漏警;第二种是基于压缩感知理论的多径幻象抑制方法，利用 降采样后的回波信号数据，对原始场景进行稀疏重构得到量测矩阵，再通过混合范数最优 化问题的求解，还原雷达成像场景，但是该方法设及到庞大的矩阵运算，导致计算复杂度的 大大增加，难W满足实时性探测成像的要求。
[0004] 穿墙雷达成像中的多径幻象严重影响建筑物内隐蔽目标的探测成像，导致最终雷 达图像质量的恶化，因此，研究穿墙雷达中多径幻象抑制具有重要的理论和实践价值。
【发明内容】
[0005] 本发明所要解决的技术问题是，提供一种穿墙雷达探测成像中能抑制多径幻象的 方法。
[0006] 本发明为解决上述技术问题所采用的技术方案是，一种适用于穿墙雷达多径幻象 抑制的成像方法，包括W下步骤：
[0007] 子孔径成像步骤:将全孔径分为N个没有交集的子孔径;对每个子孔径，分别通过 后向投影成像算法，得到N个独立的子图像，分别表示为Ii, In, i GU，2,…，N}，其中Ii 代表第i个子孔径单独作用时形成的子图像；
[000引子图像融合步骤:对N个独立的子图像分别进行归一化处理，得到N个归一化子图 像1/，1￡{1，2，-，，的；对于训畐归一化的子图像进行联乘得到融合后的图像。，。^，7)二 I/l(x，y)oI/2(x，y)…oI/N(x，y);最后对融合后的图像If进行口限检测与二值化处理得 到最终的雷达图像I。
[0009]
【申请人】发现在不同的子孔径下的子图像中，目标总是位于相同的正确的位置，而 多径幻象则位于不同的位置。基于运个特性，本发明提出一种基于子孔径的多径幻象抑制 成像方法，通过对没有交集的子孔径下的子图像进行融合，能够抑制多径幻象同时保持目 标仍在正确的位置。
[0010] 本发明的有益效果是，有效抑制了多径幻象干扰，提高了图像的信杂比和墙后隐 蔽目标的探测性能。
【附图说明】
[0011] 图1为多径幻象定位与抑制方法的处理流程图。
[0012] 图2为穿墙雷达成像场景与多径模型示意图。
[0013] 图3为多径幻象定位示意图。
[0014] 图4为不同子孔径下多径幻象位置结果图。
[0015] 图5为穿墙雷达仿真场景示意图。
[0016] 图6为全孔径多径幻象定位结果图。
[0017] 图7为全孔径多径幻象成像结果图。
[001引图8为子孔径1成像结果图。
[0019] 图9为子孔径2成像结果图。
[0020] 图10为子孔径3成像结果图。
[0021] 图11为子图像融合后结果图。
【具体实施方式】
[0022] 穿墙雷达成像场景与多径模型如图2所示。房间的长度和宽度分别为Dx和Dy,墙体 厚度为d，相对介电常数是为e，静止点目标Pt=Ut, yT)位于房间内部，合成孔径天线阵列包 含N个收发共置天线。为简单起见，天线阵列贴前墙放置。合成孔径雷达SAR是利用一个天线 阵列，沿着直线轨迹移动，并发射相参信号，把在不同位置接收的回波进行相干处理，从而 获得高分辨率成像。本文中的全孔径就是天线阵列从起始位置到终止位置的长度，即天线 移动的长度。子孔径是在运个天线合成孔径长度中选取的一段。
[0023] 假定墙壁是光滑的，电磁波在墙壁内侧发生镜面反射。考虑一阶多径回波，即目标 直接回波路径(路径-1)和=条一次反射的多径路径(路径-2、路径-3和路径-4)的组合。高 阶多径回波经过多次反射，衰减比较强，因此在多径模型中将它剔除。
[0024] 多径幻象的定位如图3所示。假设墙体参数E = I，即先不考虑前墙的影响，运里只 分析侧墙-1对应的一阶多径回波，即图3中路径-1和路径-2的组合，该分析同样适用于其它 墙体的多径回波。首先分析两个天线独立工作时的情况，天线位置坐标分别为Rm= (Xm，0)， 3。=枯，0)，1?"'瓜'和?1'分别表示天线和目标关于侦鹏-1的对称位置。
[0025] 对于天线Rn,路径-2在侧墙-1的反射点表示为Sn,由几何关系可得：
[0026] RnPT+P 巧 n+SnRn = P 讯n+P 讯n'
[0027] 则单程多径距离rn可表示为：
[002引
[0029]在图3中，在考虑对称天线的情况下，侧墙-1对应的一阶多径轨迹形成一个楠圆， 表示为双基多径轨迹W，其焦点为Rn和Rn'，长轴长度为化n;等效到单个天线作用下，多径轨 迹形成一个圆，表示为单基多径轨迹W，其圆屯、和半径分别为Rn和'。。山是双基多径轨迹W 与侧墙-I的交点，则由几何关系可得：
[0030] Tn = RnBn
[0031] 类似的分析应用于当天线位于Rm时的情况。显然，单基多径轨迹W和单基多径轨 迹<"> 的交点就是两个天线独立工作时形成多径幻象的位置，用Pm表示，排除另外一个交点 因为它不在雷达成像范围之内。
[0032] 因此，当天线位于为Rm和Rn时，侧墙-1对应一阶多径的单基多径轨迹W和单基多径 轨迹与其半径rn和rm分别表示为：
[0037] 假设Xn〉Xm，由几何关系易得rn〉rm。求解上述方程组，得到多径幻象Pm的坐标为：
[0033]
[0034]
[0035]
[0036]
[00；3 引
[0039] 因此，当收发共置天线在N个位置进行合成孔径成像时，在每个天线位置，得到一 个多径圆形轨迹。从理论上讲，N个多径圆形轨迹两两相交，得到个交点，运些交点所形 成的区域即为全孔径下多径幻象的位置。
[0040] 当e声1，即考虑前墙的情况下，多径幻象的定位方法仍然适用。每个天线位置形成 的多径圆形轨迹和半径的公式可能略有不同。然而，如果给定前墙的厚度和墙体的介电常 数，使用上述定位方法，同样可W得到多径幻象的坐标表达式。
[0041] 通过分析多径幻象坐标的解析表达式，得到多径幻象位置依赖于天线阵列位置的 特性，即多径幻象的位置与子孔径选择有关。
[0042] 一个初步仿真如图4所示，显示了不同子孔径下多径幻象的位置。=个长度相同、 没有交集的子孔径，其位置分别用不同，用上述多径定位方法得到的不同子孔径下多径幻 象区域在图中标注。从图4中我们得出，在不同的子孔径下，目标总是位于相同的正确的位 置，而多径幻象则位于不同的位置。基于运个特性，本发明提出一种基于子孔径的多径幻象 抑制方法，处理流程如图1:
[0043] -、子孔径成像步骤：
[0044] 1.将全孔径分为N个独立没有交集的子孔径；
[0045] 2.对每个子孔径，分别通过后向投影成像算法，得到N个独立的子图像，分别表示 为Ii, 12…In, iG{l，2,…，N}，其中Ii代表第i个子孔径单独作用时形成的子图像。
[0046] 二、基于子图像融合的多径幻象抑制步骤：
[0047] 由于在不同子图像中，目标总是位于相同的位置，而多径幻象的位置不同，通过子 图像的融合，能够抑制多径幻象同时保持目标仍在正确的位置，具体步骤如下：
[0048] 1.为了数据处理方便，对于步骤2所得到的N个独立的子图像，分别进行归一化处 理，得到N个归一化子图像，对于第i个子图像I/ :
[0049]
[0050] Vky)表示任意(x，y);
[0051] 2.对于N幅归一化的子图像，通过联乘的图像融合算法，即对应像素点幅度值相 乘，得到融合后的图像If，其融合表达式为：
[0化2]
[0053] 3.对融合后的图像If进行口限检测与二值化处理，W便目标识别，处理过程表示 为：
[0化4]
[0055] 其中，A是预设的二值化口限，取A = O.5是有效的经验口限值，则得到最终的雷达 图像I。
[0056] 当选择的子孔径有交集时，不同的子孔径单独成像时形成多径幻象区域会有交 集，从而导致最终成像结果中多径幻象残留，严重影响多径幻象的抑制效果和目标的检测 成像性能。因此，在本发明中选取的是独立没有交集的子孔径。同时，为了保证获得高分辨 率成像，子孔径在保证不重合的同时长度尽量大。
[0057] 下面根据一个MATLAB仿真的例子给出本发明的实施列。
[0化引仿真场景如图5所示，房间的长度和宽度Dx=Dy=IOm,墙体厚度d = 0.2m，其相对介 电常数为e = 7.5,单个静止点目标坐标为(6,4)m。合成孔径天线阵列贴前墙放置，包含51个 收发共置的天线，相邻天线间距为0.2m，天线依次进行信号收发。仿真采用超宽带步进频信 号，信号载频1.5G化、带宽IG化。
[0059] 根据本发明的处理步骤：
[0060] 第一步:根据上述多径幻象定位方法，获得全孔径下多径幻象的定位结果，如图6 所示，通过后向投影成像算法获得全孔径多径幻象的成像结果，如图7所示。
[0061] 第二步:将全孔径分为3个独立没有交集的子孔径，子孔径1对应第1个到第17个天 线位置;子孔径2对应第18个到第34个天线位置;子孔径3对应第35个到第51个天线位置。对 每个子孔径分别采用后向投影成像算法，得到S幅子图像，如图8、图9和图10所示。
[0062] 第=步:对于得到的=幅子图像进行归一化处理，再通过联乘的图像融合算法、n 限检测与二值化处理后，得到最终的成像结果如图11所示。
[0063] 仿真结果表明，本发明能够准确得到了多径幻象的位置，有效的抑制了多径幻象， 提高了图像的信杂比和墙后隐蔽目标的探测性能。
【主权项】
1. 一种适用于穿墙雷达多径幻象抑制的成像方法，其特征在于，包括以下步骤： 子孔径成像步骤:将全孔径分为N个没有交集的子孔径;对每个子孔径，分别通过后向 投影成像算法，得到N个独立的子图像，分别表示为···，《，其中I1代表 第i个子孔径单独作用时形成的子图像； 子图像融合步骤:对N个独立的子图像分别进行归一化处理，得到N个归一化子图像疒η 1^{1，2，-_少};对于_畐归一化的子图像进行联乘得到融合后的图像1^14^7) = 1/1(^ y) 0I'2(x，y)…0I' N(x，y);最后对融合后的图像1￡进行门限检测与二值化处理得到最终 的雷达图像I。2. 如权利要求1所述一种适用于穿墙雷达多径幻象抑制的成像方法，其特征在于，对融 合后的图像If进行门限检测时，二值化门限取值为0.5。
【文档编号】G01S13/90GK105954745SQ201610284105
【公开日】2016年9月21日
【申请日】2016年4月29日
【发明人】崔国龙, 阎冬, 胡露, 卢金伟, 郭世盛, 王明阳, 刘健强, 孔令讲, 杨晓波
【申请人】电子科技大学