一种用于浅埋目标高分辨率透视成像雷达杂波抑制方法
【专利摘要】本发明公开了一种用于浅埋目标高分辨率透视成像雷达的杂波抑制方法，将回波信号进行SVD分解，求得奇异值矩阵S，将第一阶奇异值置为0，可以去掉强杂波，保留目标回波信息；进一步的，对去掉强杂波后的回波信号进行PCA分解，得到特征值矩阵，通过累次方差贡献率，剔除噪声杂散，保留目标信息，对回波数据进行抑制，可去除介质表面强的反射杂波，提高目标信号的信噪比；可滤除墙体不均匀和杂质引起较小杂波，使成像效果更清晰。
【专利说明】
一种用于浅埋目标高分辨率透视成像雷达杂波抑制方法
技术领域
[0001] 本发明涉及一种用于浅埋目标高分辨率透视成像雷达的杂波抑制方法，该方法适 用于浅埋目标高分辨率透视成像雷达中的介质表面杂波抑制处理，属于雷达信号处理技术 领域。
【背景技术】
[0002] 浅埋目标高分辨率透视成像雷达工作过程中向介质内部发射电磁波，回波中除了 目标回波外，还包括介质分界表面的杂波，并且由于介质表面与探测天线距离近、波束接触 面积大，介质表面反射杂波远远强于隐蔽在介质内的目标回波，介质表面杂波抑制性能是 影响透视成像雷达成像效果的重要因素，传统深层目标探测雷达，如探地雷达，由于采用的 波长长，对介质表面的粗糙程度不敏感，可近似认为介质表面反射杂波没有空变性和时变 性，对于浅埋目标高分辨率透视成像雷达，由于需要完成埋藏目标的高分辨成像，需采用较 高频段工作频率，因而不能忽略介质表面反射杂波的空变性和时变性，因此，介质表面杂波 抑制是浅埋目标高分辨率透视成像雷达的一个难点问题。
[0003] 传统的杂波抑制方法主要有抑制法、滤波法和子空间分解法三大类方法，其中抑 制法包括空场景抑制法和均值抑制法;滤波法主要包括距离窗选通、FK滤波法;子空间分解 法分为主成分分析法(PCA)、奇异值分解法(SVD)和因子分解法(FA)。
【发明内容】
[0004] 本发明涉及一种用于浅埋目标高分辨率透视成像雷达杂波抑制方法，克服了浅埋 目标高分辨率透视成像雷达体制下介质表面杂波难以抑制的问题，抑制了介质表面反射波 干扰，提高了目标回波信号的信杂比，使成像效果更加清晰。
[0005] -种用于浅埋目标探测雷达的杂波抑制方法，包括如下步骤：
[0006] S00:对回波信号矩阵进行奇异值分解，具体为：
[0007] 对接收的mXn阶回波信号矩阵A进行奇异值分解，分解后得到m阶酉阵U，n阶酉阵V 和奇异值矩阵S，即：
[0008] A = UXSXVh (1)
[0009] 其中H表示共辄转置；
[0010] ATA的特征值为入i彡入2彡…彡入n;
[0011] 记k=min{m，n}，min{ ? }表示取最小值运算;在ATA的特征值中保留前k个特征值， 则巧=#为回波信号矩阵A的奇异值，1 = 1，2，…，k;将奇异值〇1组成奇异值矩阵S;
[0012] S01:将奇异值矩阵S中的第一阶奇异值置为0,得新奇异值矩阵S' ；
[0013] S02:根据置零后的奇异值矩阵S'对回波信号矩阵进行重构，得到去掉强杂波的回 波信号矩阵sig:
[0014] sig = UXS，XVH (3)
[0015] S03:对步骤S02得到的回波信号矩阵sig进行PCA分解，得到：
[0016] :s.ig = TmxilPl0<n =(必 + t2p2 + …+ tfPi H "nP" (:4)
[0017]其中，i =卜_1!山为第i个主元向量，各主元向量两两之间相互正交;Pl为第i个投 影向量，且为单位向量，各投影向量两两之间相互正交；
[0018] 将式(4)的等号两侧同时右乘Pi，有ti = sig ? Pi;
[0019]对于回波信号sig的协方差矩阵R，根据以上结论得：
[0021 ]对协方差矩阵R进行特征值分解，有： 乂 0 0 ? -
[0022] [病為??? k] q 言 〇 g v【](6) 0 0 0忠
[0023] 比较(5)和(6)两式，得到协方差矩阵R的特征值：
(7) 0 Q 0：" 〇 乂： 〇 0
[0025] 其中，i = l^_n;[巾1巾2…<K]表示前向矩阵，n. " A表示特征值矩 0 0 V,. 0 _0:00( 阵，vf… 表示后向矩阵；
[0026] S04:将协方差矩阵R的特征值由大到小排列，根据式(8)计算前r个特征值对应的 累次方差贡献率CVCR:
⑶
[0028]同时，使得协方差矩阵R的前r个特征值对应的累次方差贡献率CVCR小于或等于设 定阈值，保留该前r个特征值，将其它特征值剔除，形成协方差矩阵R的新的特征值矩阵;其 中，所述设定阈值根据实际需求确定；
[0029] S05:按照S02的方法，根据步骤S04中获得的协方差矩阵R的新的特征值矩阵，对回 波信号矩阵进行重构，即为杂波抑制之后的回波数据。
[0030]本发明与现有的杂波抑制方法相比，具有以下有益效果：
[0031]由于接收信号中强杂波本身具有很强的相关性，而目标回波信息与强杂波相关性 很差，所以本发明将回波信号进行SVD分解，求得奇异值矩阵S，将第一阶奇异值置为0,可以 去掉强杂波，保留目标回波信息;进一步的，对去掉强杂波后的回波信号进行PCA分解，得到 特征值矩阵，通过累次方差贡献率，剔除噪声杂散，保留目标信息，对回波数据进行抑制，可 去除介质表面强的反射杂波，提高目标信号的信噪比；可滤除墙体不均匀和杂质引起较小 杂波，使成像效果更清晰。
【附图说明】
[0032]图1为本发明的SVD-PCA级联杂波抑制处理流程图；
[0033]图2为使用现有技术中SVD方法杂波抑制处理成像结果图；
[0034]图3为使用现有技术中PCA方法杂波抑制处理成像结果图；
[0035]图4为使用本发明的SVD-PCA级联杂波抑制处理结果图。
【具体实施方式】
[0036]下面结合附图和实施例，对本发明技术方案做进一步解释。
[0037]本发明所述一种用于浅埋目标探测雷达的杂波抑制方法，如图1所示，包括如下步 骤：
[0038] S00:回波信号矩阵奇异值分解
[0039] 奇异值分解抗干扰能力较强，数据矩阵被施加较小的扰动时，数据的奇异值不会 有明显改变。针对强杂波特点，首先使用SVD分解，初步去除直达波。
[0040] 由奇异值分解定理可得:信号矩阵A进行奇异值分解，已知A为mXn阶复矩阵，分解 后得到，m阶酉阵U，n阶酉阵V和奇异值矩阵S，即：
[0041] A = UXSXVh (1)
[0042] 其中h表示共辄转置；S = diag(〇i，〇2，."，〇r )，〇i>0(i = l，???，!)，r = rank(A) ;r = rank( ?)表示求维数;ATA的特征值为心彡\2彡…彡An。酉阵U和V是复矩阵，并且满足： \UBU -1
[0043] ⑵ {VVH = I
[0044] 其中I为单位矩阵。记k = min{m，n}，min{ ? }表示取最小值运算；由奇异值分解引 理可知，在特征值Xi彡h彡…彡中保留前k个特征值，即心彡\2彡…彡Xk彡Xk+1 =…=入n = 0，则^为回波信号矩阵A的奇异值，1 = 1，2，…，k;将奇异值〇1组成奇异值矩阵S;
[0045] S01:奇异值矩阵置零处理
[0046] 在数据处理时，S的奇异值是按从大到小排列，奇异值的大小反映对应信息的大 小，最大奇异值W表示其第一阶特征值的相关成分信息。
[0047] 接收信号中，强杂波本身具有很强的相关性，而目标回波信息与强杂波相关性很 差，所以可以按此特性去掉强杂波，保留目标回波信息。求得奇异值矩阵S，将第一阶奇异值 置为0，即去掉最大杂波项，得新奇异值矩阵S '。
[0048] S02:奇异值矩阵重构回波信息
[0049]根据置零后的奇异值矩阵S'对强杂波滤除后的回波信号矩阵进行重构，得对应的 回波信号矩阵Sig:
[0050] sig = UXS，XVH (3)
[0051] S03:对步骤S02得到的回波信号矩阵sig进行PCA分解，具体为：
[0052]主元分析法(PCA)利用信号的二阶距，将信号分解到不同的子空间中，寻求不相关 的投影信息，不同子空间投影信息反映出目标的回波信息，而其权重关系按照能量大小清 晰地区分。利用该性质，对B-scan矩阵（sig矩阵)进行分解，再按照能量权重关系进行进一 步筛选、滤除，从而实现进一步去除杂波的目的，具体包括如下步骤：
[0053] 将SVD处理后的每个B-scan信息按以下形式分解：
[0054] sigmxn = TmMPln = txp[ -VUpl ++ tnpl (4)
[0055] 其中，m为酉阵U的阶数，也表示每个变量采集序列的长度，n为酉阵V，也表示数据 集中的变量个数。称ti为第i个主元向量（或称得分向量），pi为第i个投影向量（或负荷向 量)。
[0056]并做如下规定：
[0057] 1.各主元向量两两之间相互正交；
[0058] 2.各投影向量为单位向量；
[0059] 3.各投影向量两两之间相互正交。
[0060]将礙_ =^<+^^+...+4^两侧同时右乘Pl，结合三条约束条件，有 t1 = sig ? Pl，即第i个得分向量即为数据矩阵在第i个负荷向量Pl方向上的投影。对于回波 信号sig的协方差矩阵R，根据以上结论可得：
[0062]另，对R进行特征值分解，有： ~a； 0 0 0"
[0063] R :电織T …毳]:〇 ^ [vf g …vf] (6) 0 0 0 2:_ -1； 0 0 0, 0. 0. .0
[0064] 其中，i = l"_n;[巾1巾2…<K]表示前向矩阵，n A rt表示特征值矩 0 0 ?，_.0. 0 0 0 1: 阵，[vf <… 表示后向矩阵；
[0065]比较(5)和(6)两式，于是有： (7)
[0067]综上，对回波信息进行主元分析，可通过对回波数据的协方差矩阵进行分解实现。 其协方差矩阵的特征值即为数据的主元长度的期望，协方差矩阵的特征向量对应为相应主 元的投影方向。
[0068] S04:计算累次方差贡献率
[0069] 反映了 sig在Pi方向上的投影长度，即反映了 sig在Pi方向上变化的能量大小。 将特征值由大到小排列人/多A/多…多，则可去掉反映强杂波和噪声的主元投影，实现 降维的目的。若选取前r个特征值，则降维后的信息保留度可通过累次方差贡献率 (Cumulative Variance Contribution Ratio，CVCR)石角定：
(8)
[0071] 保留贡献率CVCR小于或等于设定阈值的特征值，其中，设定阈值根据实际需求确 定，本实施例中，可选取设定阈值在80 %至90 %的范围。
[0072] S05:B_Scan 数据重构
[0073] 按照S02的方法，根据步骤S04中剔除噪声杂散干扰后的特征值矩阵，对B-scan矩 阵即回波信号矩阵进行重构，即为杂波抑制之后的回波数据。
[0074] 实施例：
[0075]测试场景中，穿透的介质为5cm木板，目标为两根平行的2mm直径，5cm长度铁棍。 [0076]设备采集数据后对回波数据进行杂波抑制处理步骤如下：
[0077] S00.将接收信号矩阵奇异值分解，形式为W=UXSXVH。
[0078] SOI.将第一阶奇异值置为0,即去掉最大杂波项，得新奇异值矩阵S'。
[0079] S02.将置零后的奇异值矩阵还原，得S'对应的新信息矩阵
[0080] sig = UXS，XVH
[0081 ] S03.将重构矩阵PCA分解处理： 0 0 0 ~ ^ r -I 0 ^ri 0 0 「!T, ~I
[0082] 及= (&[=[>, 4 - k] 〇 〇 ^ 〇 v27 …v；：] _〇 0 .0 .14」
[0083] S04计算累次方差贡献率：
[0084] 可通过CVCR数值的大小，看出该特征所包含能量的大小。根据每个B-scan数据的 CVCR，由第一主元回波能量强弱，判断是否包含较多强杂波信息，若回波能量过强则剔除， 能量适中予以保留。同时，对于主元投影中代表噪声信息的极小能量，也可进一步去除。
[0085] S05 步骤 S05:B_Scan 数据重构
[0086] 对进行B-scan重构，可得：
[0087] Sig recon Sig[伞 ii ??? <K][伞 ii ??? <K]t
[0088]重构后即为杂波抑制之后数据。
[0089] 单独使用SVD进行杂波抑制处理，其成像后的结果如附图2所示，能够分辨目标，但 其仍留有较大杂波旁瓣，视觉上可以明确区分。单独使用PCA进行杂波抑制处理，其成像后 的结果如附图3所示，采用SVD-PCA级联处理后的成像结果如附图4所示。
[0090] 对于三种方法的成像效果，采用信杂比（SCR)指标进行定量评价，其结果如下表所 示，从表中可以看出采用SVD-PCA级联的杂波抑制方法后的成像结果SCR值明显高于单独使 用SVD或者PCA的杂波抑制方法。
[0091] 表1 SVD、PCA、SVD-PCA三种方法抑制性能评估
【主权项】
1. 一种用于浅埋目标探测雷达的杂波抑制方法，其特征在于，包括如下步骤： S00:对回波信号矩阵进行奇异值分解，具体为： 对接收的mXn阶回波信号矩阵A进行奇异值分解，分解后得到m阶酉阵U，n阶酉阵V和奇 异值矩阵S，即： A = UXSXVh (1) 其中Η表示共辄转置； ΑΤΑ的特征值为λι多λ2多…多λη; 记k = min{m，n}，min{ · }表示取最小值运算；在ΑΤΑ的特征值中保留前k个特征值，贝lj为回波信号矩阵A的奇异值，1 = 1，2，…，k;将奇异值σι组成奇异值矩阵S; SOI:将奇异值矩阵S中的第一阶奇异值置为0，得新奇异值矩阵S' ； S02:根据置零后的奇异值矩阵S'对回波信号矩阵进行重构，得到去掉强杂波的回波信 号矩阵sig: sig = UXS，XVH (3) S03:对步骤S02得到的回波信号矩阵sig进行PCA分解，得到： 啦= 7Lx"€L …+U)/ +…+ 匕尺 (4) 其中，i = 1H_n，tl为第i个主元向量，各主元向量两两之间相互正交;Pl为第i个投影向 量，且为单位向量，各投影向量两两之间相互正交； 将式(4)的等号两侧同时右乘pi，有ti = sig · pi; 对于回波信号sig的协方差矩阵R，根据以上结论得：对协方差矩阵R进行特征值分解，有：比较(5)和(6)两式，得到协方差矩阵R的特征值：其中，i = 1···η ; [ Φ 1 Φ 2 ··· Φ η]表示前向矩阵表示特征值矩阵， [X <…表示后向矩阵； S04:将协方差矩阵R的特征值由大到小排列，根据式(8)计算前r个特征值对应的累次 方差贡献率CVCR:同时，使得协方差矩阵R的前r个特征值对应的累次方差贡献率CVCR小于或等于设定阈 值，保留该前r个特征值，将其它特征值剔除，形成协方差矩阵R的新的特征值矩阵;其中，所 述设定阈值根据实际需求确定； S05:按照S02的方法，根据步骤S04中获得的协方差矩阵R的新的特征值矩阵，对回波信 号矩阵进行重构，即为杂波抑制之后的回波数据。
【文档编号】G01S13/89GK105929373SQ201610239480
【公开日】2016年9月7日
【申请日】2016年4月18日
【发明人】刘海波, 丁泽刚, 龙腾, 王珣
【申请人】北京理工大学