一种穿墙雷达成像后多目标侧墙一次多径抑制方法
【专利摘要】本发明公开一种适用于穿墙雷达成像的多目标侧墙一次多径抑制方法，首先利用4连通域检测方法从二元原始雷达图像中提取出所有连通域；然后对所有连通域编号并依照它们的平均坐标离雷达阵列中心的欧式距离升序排序得到连通域向量；接着把连通域向量元素重组成两个子向量，分别按照到雷达的多径传输路径和直接回波传输路径对这两个子向量进行扩展；然后利用相似度函数计算两个子向量的相似度矩阵；最后使用相似度门限及判决规则的方法判定相似度矩阵的所有元素，得到真实目标所在的连通域并将多径假目标抑制掉。本发明能够对墙后隐蔽多目标进行检测，在保留目标的同时能够有效地将多径剔除。
【专利说明】
一种穿墙雷达成像后多目标侧墙一次多径抑制方法
技术领域
[0001 ]本发明涉及穿墙雷达成像技术领域，特别是涉及部分目标与其他目标的侧墙一次 多径重合时多径假目标的抑制方法。
【背景技术】
[0002] 穿墙雷达成像技术主要是利用电磁波穿透建筑墙体等障碍物，对密闭建筑物内隐 藏目标进行检测、定位和识别，同时对建筑物布局进行成像，在反恐、公安执法、灾难救援等 领域具有重大的应用价值。电磁波在密闭的建筑环境内传播时，存在许多不同的传播路径， 其中从发射雷达打到目标然后直接反射回接收雷达的路径称为直接路径，由于建筑物内光 滑墙体的存在，电磁波从发射雷达发射出去之后会在墙体表面进行一次或者多次反射后再 经过目标反射回接收雷达，或者先打到目标，然后在反射到内墙面，最后反射回接收雷达， 这都就会造成多径假目标的产生，因为电磁波在每一次墙面反射时强度都会衰减，所以二 次及高次多径强度都很弱，通常只考虑一次多径。多径假目标会导致虚警，严重影响穿墙雷 达的成像质量。因此，密闭建筑物中穿墙成像多径幻象抑制十分重要。
[0003] 国内外许多研究机构开展了密闭空间一次多径干扰抑制方法的研究。电子科技大 学提出了一种基于子孔径成像的多径抑制方法(Z.X.Li，Y.Jia,et al，"A novel approach of multi-path suppression based on sub-aperture imaging in through-wall-radar imaging"，IEEE Radar Conference,2013.)，它利用多径的位置随着雷达位置的改变而改 变这个特性，选取不同位置的子孔径进行成像，最后将这些子孔径图像融合在一起达到多 径抑制的目的。美国维拉诺瓦大学提出了一种基于点扩散函数的多径抑制方法(Setlur P， Alii G,Nuzzo L.Multipath exploitation in through-wall radar imaging via point spread functions[J]. IEEE Transaction on Image Processing. 2013 ·)，它把多径假目 标的值叠加到与它关联的真实目标上，在抑制多径的同时能够提高目标位置的信杂比。上 述方法只能够应对单目标场景，或者多径与目标之间不存在重叠的场景，当处理存在目标 与多径假目标重叠的多目标场景时，上述方法会造成多径假目标区域内的真实目标丢失， 出现漏警。因此，研究密闭空间中存在重叠现象的多目标多径干扰抑制方法在穿墙雷达成 像领域具有重要的价值。
【发明内容】
[0004] 本发明为解决上述技术问题，提出了一种穿墙雷达成像后多目标侧墙一次多径抑 制方法，首先通过门限检测将原始雷达图像进行二值化，使用连通域检测方法对该二值化 图像进行连通域检测得到连通域向量，基于多径的产生原理使用相似度方法分析该连通域 向量中两两元素间的关联性，判断出目标区域和多径区域，最后在保留目标的同时将多径 假目标剔除，从而实现了多径抑制。
[0005] 本发明采用的技术方案是：一种穿墙雷达成像后多目标侧墙一次多径抑制方法， 包括：
[0006] S1、连通区域检测，根据设置的第一门限对原始雷达图像进行门限检测得到二元 图像，对二元图像进行4-连通检测，并标记二元图像中的所有的N个连通域，按照平均坐标 值到阵列中心的欧式距离由小到大对所有的N个连通域进行排序，得到一个连通域向量； [0007] S2、构造相似度矩阵，将步骤S1得到的连通域向量分成两组，得到第一向量与第二 向量;根据侧墙一次多径的形成原理对第一向量进行扩展，得到包含侧墙一次多径传输距 离信息的第三向量;根据雷达的直接回波传输路径对第二向量进行扩展，得到包含直接回 波传输距离信息的第四向量;根据相似度函数计算得到第三向量与第四向量的相似度矩 阵；
[0008] S3、元素判决及多径抑制，根据第二门限判定相似度矩阵中的所有元素，得到目标 区域。
[0009] 进一步地，步骤S2中所述将步骤S1得到的所有连通域向量分成两组，具体为:所述 第一向量为按照各个连通域的平均坐标到阵列中心的欧式距离由小到大的顺序连续地在 连通域向量中取连通域，所述第二向量为按照平均坐标值到阵列中心的欧式距离由大到小 的顺序连续地在连通域向量中取连通域，所述第一向量与第二向量所取连通域个数相等， 记为Υ大于连通域向量中总的连通域个数的一半，并且小于总的连通域个数。
[0010] 进一步地，步骤S2中所述的根据侧墙一次多径的形成原理对第一向量进行扩展， 得到包含侧墙一次多径传输距离信息的第三向量;具体为:计算第一向量中每个连通域的 平均坐标到K个雷达位置的侧墙一次多径传输距离，每个连通域得到的K个距离值组成一个 向量ri; i = 1，2，3···，N、根据Μ个向量η得到包含侧墙一次多径传输距离信息的第三向量。
[0011] 步骤S2中所述的根据雷达的直接回波传输路径对第二向量进行扩展，得到包含直 接回波传输距离信息的第四向量;具体为:计算第二向量中每个连通域的平均坐标到K个雷 达位置的直接回波距离，每个连通域得到的K个距离值组成一个向量打'；V = 1，2，3···，Y， 根据Μ个向量Γι〃得到包含Y个直接回波传输距离信息的第四向量。
[0012] 进一步地，步骤S2中所述的根据相似度函数计算得到第三向量与第四向量的相似 度矩阵，具体为:根据相似度函数定义，第三向量与第四向量的相似度矩阵表达式为：
[0014] 其中，sinKri，;^ 2)表示计算向量ri和V 2之间相似度的相似度计算函数，其取值范 围为[0，1]。
[0015] 进一步地，步骤S3所述根据第二门限判定相似度矩阵中的所有元素，得到目标区 域，具体为：
[0016] A1、对相似度矩阵S进行预处理，根据多径回波路径始终长于目标回波路径的特 点，将相似度矩阵S中的上三角元素全部置零，得到一个下三角矩阵f
[0018] A2、将矩阵S'中的所有非零元素按照取值从大到小依次放入行向量Sv:
[0020] A3、假如只考虑侧墙一次多径的影响，那么对于连通域向量D中的N个元素，最多存 在N-1个目标，因此选取向量Sv中的前N-1个元素用来计算用于相似度检测的第二门限；
[0021] A4、将矩阵S'中大于第二门限的元素提取出来，并对提取出的元素按照以下步骤 进行判决：
[0022] B1、如果Sim(rk，rV)高于门限值，且矩阵S'中不存在第一个自变量为〇的元素， 则将rk对应的连通域判定为目标区域，将rV对应的连通域判定为多径区域；
[0023] 132、如果同时存在8；[111(^，1'/[/)和8；[111(1' 1)，1'%0高于门限值，且1'/1/与1'1)对应相同的 连通域;根据s im(rk，r7 p0将rk对应的连通域判定为目标区域，将r7 [/对应的连通域判定为 多径区域;然后根据将&对应的连通域判定为目标区域，将对应的连通域 判定为多径区域;最终得到rk对应的连通域与rV对应的连通域均为目标区域，V 对应的 连通为多径区域；
[0024] B3、当矩阵S'中所有高于第二门限的元素判定完成之后，将判定为目标区域的连 通域保留，同时将判定为多径区域的连通域置零，得到抑制了多径的后的二值化图像If (·);最后把原始图像1。（ ·）与IK ·)进行乘法图像融合得到一幅保留所有目标并且消除 了多径的雷达图像Ifirml( ·)。
[0025] 更进一步地，所述第二门限的计算式如下：
[0027] 其中，j'=1，2,3···Ν_1。
[0028] 本发明的有益效果：本发明的一种穿墙雷达成像后多目标侧墙一次多径抑制方 法，首先利用4连通域检测方法从二元原始雷达图像中提取出所有连通域;然后对所有连通 域编号并依照它们的平均坐标离雷达阵列中心的欧式距离升序排序得到连通域向量;接着 把连通域向量元素重组成两个子向量，分别按照到雷达的多径传输路径和直接回波传输路 径对这两个子向量进行扩展;然后利用相似度函数计算两个子向量的相似度矩阵；最后使 用相似度门限及判决规则的方法判定相似度矩阵的所有元素，得到真实目标所在的连通域 并将多径假目标抑制掉。本发明的方法能够有效识别雷达图像中的真实目标和它们的一阶 侧墙多径假目标，在抑制多径假目标的同时不会造成真实目标的丢失。相比于现有技术中 的多径抑制方法，本发明能够防止正常信噪比条件下目标的丢失，形成高质量的雷达图像， 保证了穿墙雷达在密闭建筑物探测环境中的实用效果，为操作人员做出正确的决策提供了 有力的保障。
【附图说明】
[0029] 图1为侧墙一次多径抑制方法处理流程图。
[0030] 图2为侧墙一次多径穿墙传播示意图。
[0031] 图3为【具体实施方式】中仿真场景示意图。
[0032] 图4为【具体实施方式】中原始雷达仿真图像。
[0033] 图5为【具体实施方式】中4连通域检测及标记结果图。
[0034] 图6为【具体实施方式】中多径抑制后二值化雷达图像。
[0035] 图7为【具体实施方式】中最终结果图像。
【具体实施方式】
[0036] 为便于本领域技术人员理解本发明的技术内容，下面结合附图对本
【发明内容】
进一 步阐释。
[0037] 如图1所示为本发明的处理流程图，本发明的技术方案为:一种穿墙雷达成像后多 目标侧墙一次多径抑制方法，包括：
[0038] S1、连通区域检测，根据设置的第一门限对原始雷达图像进行门限检测得到二元 图像，对二元图像进行4-连通检测，并标记二元图像中的所有连通域，按照平均坐标值到阵 列中心的欧式距离由小到大对所有连通域进行排序，得到一个连通域向量。
[0039]设置第一门限THb，用THb对原始雷达图像1。（ ·）进行门限检测得到二元图像Ib (·)，对Ib( ·)进行4-连通（同一像素在上、下、左、右四个方向有连通)检测，并标记图像Ib (·)中的所有连通域。假设有N个连通域，计算每个连通域的平均坐标(W,)，= …，， 并按照平均坐标值到阵列中心的欧式距离由小到大对所有连通域进行排序，将排好序的连 通域以符号〇1-〇〃进行编号，得到一个连通域向量0=[01，02，"_，0"。
[0040] S2、构造相似度矩阵，将步骤S1得到的连通域向量分成两组，得到第一向量与第二 向量;根据侧墙一次多径的形成原理对第一向量进行扩展，得到包含侧墙一次多径传输距 离信息的第三向量;根据雷达的直接回波传输路径对第二向量进行扩展，得到包含直接回 波传输距离信息的第四向量;根据相似度函数计算得到第三向量与第四向量的相似度矩 阵。多径穿墙传播如图2所示，第三向量即由虚线表示的一阶多径，第四向量即由实现表示 的直接路径。
[0041] 将连通域向量D分成两组，具体为:所述第一向量为按照各个连通域的平均坐标值 到阵列中心的欧式距离由小到大的顺序连续地在连通域向量中取连通域，所述第二向量为 按照平均坐标值到阵列中心的欧式距离由大到小的顺序连续地在连通域向量中取连通域， 所述第一向量与第二向量所取连通域个数相等，且第一向量或第二向量所取连通域的个数 Y大于连通域向量中总的连通域个数的一半，并且小于总的连通域个数。假设第一组中包 含元素 Di-Dn-i，用向量T = [Di，D2，…，Dn-i ]表示，即用T表示第一向量;第二组包含元素 D2_Dn， 用向量M= [D2，D3，…，Dn]表示，其中Dj代表第j个连通域，即用Μ表示第二向量;且第一向量T 与第二向量11所取连通域个数相等，记为Μ个，即在本实施例中Μ = Ν-1。
[0042]接下来将向量T和Μ进行扩展;首先，根据侧墙一次多径的形成原理，计算向量T中 每个元素所代表的连通域的平均坐标到各个雷达之间的侧墙一次多径传输距离。如:对于 向量τ中第i个连通域来说，它的平均坐标为(1，又）。假设雷达总数为Κ，第k个雷达的位置坐 标为( XRk，0)，则这个雷达到该连通域的多径传输距离为：
[0044] 其中，d和ε分别为前墙厚度以及相对介电常数，XwR表的是左侧墙内表面的横坐 标值，Axe为电磁波在从侧墙反射回接收雷达时在前墙中发生的水平偏移值，△ XT为电磁波 从目标反射回接收雷达时在前墙中发生的水平偏移值，天:表示向量T中第i个连通域的平均 横坐标，義.表示向量T中第i个连通域的平均纵坐标，XRk表示第k个雷达的横坐标，k = 1，2， 3， · · · ，K〇
[0045] 这样Τ中每个连通域对应Κ个雷达就有Κ个距离值，将第i个连通域计算得到的Κ个 距离值放进向量^中。对Τ中的每一个连通域都进行这样的距离计算后，得到一个包含侧墙 一次多径传输距离信息的向量f，1^的维度为1 X (N-l)K，将f表示如下：
[0047]其次，计算Μ中的每个连通域的平均坐标到K个雷达的直接回波距离，其中直接回 波距离的计算公式为：
[0049] 其中，5表示向量Μ中第i/个连通域的平均横坐标，JV表示向量Μ中第i/个连通域 的平均纵坐标。
[0050] 同样地，对Μ中的每一个元素都进行这样的距离计算后得到一个包含直接回波传 输距离信息的向量W的维度为1 X (N-l)K，将，表示如下：
[0052]定义任意两个向量X和y的相似度函数为：
[0054]其中，J代表向量中元素的个数。
[0055]基于这个相似度函数，f和W之间的相似度矩阵可以表示为：
[0057]由于多径回波路径始终长于目标回波路径，那么可将相似度矩阵S中的上三角元 素全部置零，得到一个下三角矩阵： ~乂一/一 J\ A
A _
[0059] S3、元素判决及多径抑制，根据第二门限判定相似度矩阵中的所有元素，得到目标 区域。
[0060] 为了便于后续分析，将矩阵S'的所有非零元素按照取值由大到小依次放入行向量 Sv：
[0062]假如只考虑侧墙一次多径的影响，那么对于连通域向量D中的N个元素，最多存在 N-1个目标，因此选取向量Sv中的前N-1个元素用来计算用于相似度检测的第二门限THsim， 第二门限的计算公式为：
[0064] 其中，太=1，2,3…N-1。
[0065] 将相似度矩阵f中大于第二门限的元素提取出来，对这些提取出来的元素按照下 面的判决规则进行判决。
[0066] B1、如果Sim(rk，rV)高于门限值，且矩阵S'中不存在第一个自变量为〇的元素， 则将r k对应的连通域判定为目标区域，将rV对应的连通域判定为多径区域；
[0067] 132、如果同时存在8；1111(^，1'/[/)和8；[111(1'1)，1'%0高于门限值，且1' /1/与1'1)对应相同的 连通域，即对应的是连通域向量D中的同一个连通域;根据^^(^，^。将^对应的 连通域判定为目标区域，将rV对应的连通域判定为多径区域;然后根据 对应的连通域判定为目标区域，将V 对应的连通域判定为多径区域;最终得到rk、rV各自 对应的连通域均为目标区域，V 对应的连通为多径区域；
[0068] B3、当矩阵S'中所有高于第二门限的元素判定完成之后，将判定为目标区域的连 通域保留，同时将判定为多径区域的连通域置零，得到抑制了多径后的二值化图像IK ·); 最后把原始图像1。（ ·）与If( ·)进行乘法图像融合得到一幅保留所有目标并且消除了多径 的雷达图像Ifinai( ·)。
[0069] 下面根据一个MATLAB仿真例子给出本发明的【具体实施方式】。
[0070] 对存在多个目标的待探测区域，其仿真场景如图3所示，坐标系零点位于前墙前表 面的中心位置，31个收发共置雷达紧贴前墙放置，相邻两个雷达之间的水平间隔为0.1m，雷 达的发射信号为中心频率1.5GHz、带宽1GHz的步进频信号，频率步进为2MHz。房间内表面的 长度和宽度均为6m，四面墙的厚度与相对介电常数相等且为0.24m和7.6。三个目标分别位 于(-0.5，3.0)m、（-1.5，4.0 )m和(-2.4，3.9 )m处，值得注意的是目标3与目标1的侧墙一次多 径重合，原始雷达仿真图像如图4所示。
[0071] S1、使用门限THb对原始雷达仿真图像进行二值化检测，门限值取为雷达图像最大 值像素点值的0.4倍，在这里它的取值为1.9134。门限检测后得到二值化的图像I b( ·)如图 5所示。通过4-连通域检测的方法对图像Ib( ·)进行检测并将所有连通域标记，可以得到五 个连通域，其中，连通域1包含104个像素点，这些像素点的坐标分别为(-3.1，4.065)111，（-3.085，4.065 )nr·· (-3.055，4.245 )m，然后根据这些像素点坐标，计算求得连通域1的均值坐 标为(_3.085,4.14)111，根据仿真场景示意图所建立的坐标系，雷达阵列中心位于(0,0)111，那 么可以求得连通域1到雷达阵列中心点的距离为5.1630m。同理，可以得到其他连通域的坐 标均值与到雷达阵列中心的距离，如表1所示：
[0072] 表1其他连通域的坐标均值与到雷达阵列中心的距离
[0074]将所有连通域按照其平均坐标到雷达阵列中心的欧式距离由小到大进行排序，得 到五个连通域，依次编号为Di-Ds，如图5所示。然后将这些连通域标号放进一个连通域向量D 中。则有下面的对应关系:Di对应连通域5，D2对应连通域4，D3对应连通域3，D4对应连通域1， D5对应连通域2。
[0075] S2、将向量D分成向量T和向量Μ，其中向量T包含连通域DX-D4，向量Μ包含连通域D2- D5。对向量T进行扩展，计算T中的每一个连通域的中心坐标到31个雷达的侧墙一次多径传 输距离。这样T的每一个元素就扩展成了一个1 X 31维的向量ri，i = 1，2，3，4，扩展之后的向 量记作f = [ Γ1Γ2Γ3Γ4]。同理，对向量Μ进行扩展，计算Μ中的每一个连通域的中心坐标到31 个雷达的直接路径传输距离，则Μ的每一个元素也扩展成了一个1 X 31维的向量V i，i = 2， 3，4，5，M向量扩展之后记作，=[V Y 4]。根据前面定义的相似度函数对向量f和， 进行相似度计算，并且得到一个4 X 4维的相似度矩阵S，S的值如表2所示：
[0076]表2相似度矩阵S的值
[0078]因为目标的多径传播距离总是长于它的直接路径传播距离，则相似度矩阵S中对 角线以上的元素全部近似为0后得到一个下三角矩阵的值如表3所示：
[0079] 表3矩阵S'的值
[0081] S3、元素判决及多径抑制
[0082] 将下三角矩阵f中的元素按值的大小进行降序排序得到一个IX 10维的行向量 Sv,
[0084] 然后对该向量的前4个元素进行求均值计算后得到相似度检测门限THsim的值为 0.8932，将矩阵S'中大于门限^心的元素提取出来，它们依次为0.9840、0.9150和0.8990。 所对应的连通域相似度函数为sinKn，;^ 5)和sinKn，:^ 4)，由判决准则可得连 通域D^DdPDs为真实目标，D4PD5分别为连通域D#PD2的侧墙一次多径假目标，连通域〇4勺 侧墙一次多径假目标与连通域D3重合。经过多径抑制后的二值化雷达图像如图6所示，原始 雷达仿真图像与多径抑制后的二值化雷达图像的非相干融合结果如图7所示。
[0085] 由仿真结果可知，本发明提供的适用于穿墙雷达成像后多目标侧墙一次多径抑制 方法不仅能够有效剔除多径，而且不会造成目标的丢失，验证了本发明的正确性和有效性。 [0086]本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发 明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领 域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的 任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。
【主权项】
1. 一种穿墙雷达成像后多目标侧墙一次多径抑制方法，其特征在于，包括： 51、 连通区域检测，根据设置的第一口限对原始雷达图像进行口限检测得到二元图像， 对二元图像进行4-连通检测，并标记二元图像中的所有的N个连通域，按照平均坐标值到阵 列中屯、的欧式距离由小到大对所有的N个连通域进行排序，得到一个连通域向量； 52、 构造相似度矩阵，将步骤S1得到的连通域向量分成两组，得到第一向量与第二向 量;根据侧墙一次多径的形成原理对第一向量进行扩展，得到包含侧墙一次多径传输距离 信息的第Ξ向量;根据雷达的直接回波传输路径对第二向量进行扩展，得到包含直接回波 传输距离信息的第四向量;根据相似度函数计算得到第Ξ向量与第四向量的相似度矩阵； 53、 元素判决及多径抑制，根据第二口限判定相似度矩阵中的所有元素，得到目标区 域。2. 根据权利要求1所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法，其特征 在于，步骤S2中所述将步骤S1得到的连通域向量分成两组，具体为:所述第一向量为按照平 均坐标值到阵列中屯、的欧式距离由小到大的顺序连续地在连通域向量中取连通域，所述第 二向量为按照平均坐标值到阵列中屯、的欧式距离由大到小的顺序连续地在连通域向量中 取连通域，所述第一向量与第二向量所取连通域个数相等，记为沪，沪大于连通域向量中总 的连通域个数的一半，并且小于总的连通域个数。3. 根据权利要求2所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法，其特征 在于，步骤S2中所述的根据侧墙一次多径的形成原理对第一向量进行扩展，得到包含侧墙 一次多径传输距离信息的第Ξ向量;具体为:计算第一向量中每个连通域的平均坐标到K个 雷达位置的侧墙一次多径传输距离，每个连通域得到的K个距离值组成一个向量ri;i = l， 2,3···，Ν/，根据妒个向量ri得到包含侧墙一次多径传输距离信息的第Ξ向量。4. 根据权利要求3所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法，其特征 在于，步骤S2中所述的根据雷达的直接回波传输路径对第二向量进行扩展，得到包含直接 回波传输距离信息的第四向量;具体为:计算第二向量中每个连通域的平均坐标到K个雷达 位置的直接回波距离，每个连通域得到的K个距离值组成一个向量ri'/ =1，2,3···，Ν/，根 据妒个向量η' /得到包含沪个直接回波传输距离信息的第四向量。5. 根据权利要求4所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法，其特征 在于，步骤S2中所述的根据相似度函数计算得到第Ξ向量与第四向量的相似度矩阵，具体 为:根据相似度函数定义，第Ξ向量与第四向量的相似度矩阵表达式为：其中，sim(ri，r/2)表示计算向量ri和之间相似度的相似度计算函数，其取值范围为 [〇，1]。6. 根据权利要求1所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法，其特征 在于，步骤S3所述根据第二口限判定相似度矩阵中的所有元素，得到目标区域，具体为： A1、对相似度矩阵S进行预处理，根据多径回波路径始终长于目标回波路径的特点，将 相似度矩阵s中的上Ξ角元素全部置零，得到一个下Ξ角矩阵S/A2、将矩阵S/中的所有非零元素按照取值从大到小依次放入行向量Sv:A3、假如只考虑侧墙一次多径的影响，那么对于连通域向量D中的N个元素，最多存在N- 1个目标，因此选取向量Sv中的前N-1个元素用来计算用于相似度检测的第二口限； A4、将矩阵S/中大于第二口限的元素提取出来，并对提取出的元素按照W下步骤进行 判决： B1、如果sim(rk，rV )高于口限值，且矩阵少中不存在第一个自变量为。的元素，则将η 对应的连通域判定为目标区域，将rV对应的连通域判定为多径区域； B2、如果同时存在sim(rk，rV )和sim(rp,;r%')高于口限值，且rV与。对应相同的连通 域;根据sim(rk，rV)将η对应的连通域判定为目标区域，将rV对应的连通域判定为多径 区域;然后根据sim(rp，rV)将。对应的连通域判定为目标区域，将rV对应的连通域判定 为多径区域;最终得到η对应的连通域与rV对应的连通域均为目标区域，rV对应的连通 为多径区域； B3、当矩阵S/中所有高于第二口限的元素判定完成之后，将判定为目标区域的连通域 保留，同时将判定为多径区域的连通域置零，得到抑制了多径后的二值化图像If( ·);最后 把原始图像1。（ ·）与If( ·)进行乘法图像融合得到一幅保留所有目标并且消除了多径的雷 达图像Ifinal( ·)。7.根据权利要求1或6所述的一种穿墙雷达成像后多目标侧墙一次多径抑制方法，其特 征在于，所述第二口限的计算式如下：其中，j'=l，2,3…N-1。
【文档编号】G01S13/88GK106093933SQ201610698090
【公开日】2016年11月9日
【申请日】2016年8月22日 公开号201610698090.5, CN 106093933 A, CN 106093933A, CN 201610698090, CN-A-106093933, CN106093933 A, CN106093933A, CN201610698090, CN201610698090.5
【发明人】崔国龙, 郭世盛, 刘健强, 卢金伟, 王明阳, 孔令讲, 杨晓波
【申请人】电子科技大学