一种用于合成孔径雷达的目标特征提取方法
【专利摘要】本发明属于雷达技术领域，涉及一种用于合成孔径雷达的目标特征提取方法。本发明联合图像域的峰值区域分割与信号域的稀疏贝叶斯学习进行特征提取，首先通过PRS算法分割强散射区域，以避免强散射点对SBL算法的影响，然后利用SBL算法自主学习提取弱散射特征点，弥补了基于图像域的PRS算法对弱目标分割的不足。该方法充分利用了PRS算法和SBL算法各自的优势，对弱目标特征提取表现出了很好的效果，具有可调参数少，鲁棒性强等优点；并且在PRS算法分割强散射区域后，测量矩阵的维数也相应地下降，减小了运算量，提高了运算效率，从而可以广泛应用于SAR目标特征提取领域，为后续自动目标识别(ATR)提供基础。
【专利说明】
一种用于合成孔径雷达的目标特征提取方法
技术领域
[0001]本发明属于雷达技术领域，涉及一种用于合成孔径雷达(SAR)的目标特征提取方 法。
【背景技术】
[0002] 合成孔径雷达(Synthetic Aperture Radar，SAR)目标特征提取是SAR信号处理的 一个非常重要的研究领域，其准确性直接影响目标识别的准确性，因此精确提取目标特征 具有重要意义。特征提取的本质是将采样数据从高维空间映射到低维特征空间，现有的SAR 特征提取方法可以分为两大类:第一类基于图像域，通过奈奎斯特采样定理和匹配滤波理 论进行特征提取，如分水岭算法等。这类算法仅使用了幅度信息而没有用相位信息，因此很 容易受到杂波和噪声的影响，从而使特征提取准确率下降。第二类方法基于信号域，通过某 种映射直接从信号中获得特征，如基于压缩感知(Compressive Sensing，CS)理论的稀疏重 构算法等。目前研究较多的CS算法主要分为三大类:凸松弛类算法、贪婪类算法与贝叶斯类 算法。对于凸松弛类算法，其主要思想是将4范数非凸问题等价为4范数凸问题来求解，如 基追踪算法(B a s i s P u r s u i t，B P)、内点法等，这类算法计算复杂，用于大数据时计算效率 低，不能满足实时性需求;对于贪婪类算法，其主要思想是通过每次迭代得到一个局部最优 解来逐步逼近原始信号，如匹配追踪算法(Matching Pursuit，MP)、正交匹配追踪算法 (Orthogonal Matching Pursuit ,ΟΜΡ)，这类算法计算速度快，实现容易，但往往需要知道 目标场景的稀疏度，特征提取结果很大依赖于稀疏度的设置，且其不考虑场景的空间结构 信息，易受噪声影响；而贝叶斯类算法，如稀疏贝叶斯学习算法（Sparse Bayesian Learning，SBL)，其主要思想是基于高斯先验假设，可以自主的学习并确定算法中的所有超 参数，获得更稀疏的解，该算法引入了目标的空间结构信息，有一定的抗噪能力，且与凸优 化算法相比具有更小的计算复杂度。然而，实际情况下，CS算法很容易受到强点的干扰，使 得对弱散射特征的提取性能不太理想，同时，当目标场景过大时，构造的测量矩阵的维数过 于庞大，使得实际处理时计算量较大，效率较低，不便于有效应用。。
【发明内容】
[0003] 本发明所要解决的，就是针对上述问题，提出一种用于合成孔径雷达的目标特征 提取方法。
[0004] 本发明的技术方案是:一种用于合成孔径雷达的目标特征提取方法，其特征在于， 包括以下步骤：
[0005] a.构建聚束SAR的成像模型如下公式1所示：
[0006] g=?y+n(公式 1)
[0007]公式1中，g为回波向量，μ是场景散射系数向量，Θ为观测矩阵，η表示噪声向量；
[0008] b.由极坐标成像算法得到目标场景的成像图像；
[0009] c.根据获得的成像图像，基于图像域的PRS(Peak Region Segmentation,联合图 像域的峰值区域分割)算法对强散射点进行特征提取，具体方法为：
[0010] Cl.对获得的成像图像进行分割，具体为：
[0011] Cl 1.设置强散射点门限τι，利用分水岭算法对强散射区域进行初始分割，得到心个 初始分割区域；
[0012] C12.设置合并门限τ2,对步骤ell得到的初始分割区域进行合并，当任意两个相邻 区域边界上像素点的幅值低于合并门限时，认为这两个相邻区域是由噪声和纹理杂波引起 的过分割导致的，因此合并这两个相邻区域，最后得到R 2个合并后的分割区域，R2彡R1;
[0013] c2.获取分割区域所对应的强散射点位置和幅度估计，并将位置索引放入集合Ω 中；
[0014] d.获取成像图像残余散射区域回波，具体为：
[0015] dl.设原始回波信号的矩阵形式如下公式2:
[0016] g = Θ (Ai+A2)y+n = gi+g2+n(公式2)
[0017] 公式2中，如为区域选择矩阵，其形式为将单位矩阵I中与集合Ω中的原子所对应 的列置为零，Ai+ApI，由此得到维数较小的测量矩阵， gl= ΘΑιμ为强散射区域回波，g2= Θ Α2μ为残余散射区域回波；
[0018] d2.将残余散射区域回波方程等价描述为如下公式3:
[0019] g2=?〇w(公式 3)
[0020] 公式3中，Θ 〇表示残余散射区域的测量矩阵，W为相应的散射函数；
[0021 ] e.利用基于信号域的SBL(Sparse Bayesian Learning,稀疏贝叶斯学习）算法对 残余信号的弱散射点进行特征提取，具体为：
[0022] el.基于高斯先验假设可得，所述残余散射区域的散射函数w的分布如下公式4:
[0023] />:〇; V)~.V.)(公式 4 >
[0024] 公式4中，V = 尽(ν,?，M为目标场景中残余散射区域等效散射中心个数，Vl为一 个非负的参数，控制着w的稀疏度，当vi = 0时，表示Wi为零；
[0025] e2.假设噪声也满足高斯分布，结合公式4可得，w的后验分布可以表示成均值为 q =g2，方差为G = [V -1 J-1的高斯分布，其中，β是噪声的方差；当超参数v 和β被估计出来后，w的点估计由其后验分布均值确定;超参数ν和β由第二类最大似然估计 获得，即等效于计算公式5所示函数的极值：
[0026]
[0027] 公式5中，I为单位矩阵；
[0028] f.采用期望最大化算法，得到参数ν的学习规则如下公式6和β的学习规则如下公
式7:
[0029]
[0030]
[0031] 其中，Ν是残余散射区域回波向量g2的维数；
[0032] g.根本步骤c中获得的强散射点和步骤e-f中获得的弱散射点，得到SAR目标特征 提取结果。
[0033]进一步的，所述步骤a的具体方法为：
[0034] al.根据电磁理论，当目标散射体的电尺寸远大于波长时，目标的高频电磁散射特 性可以由局部散射特性的合成来表示，这些局部散射通常被称为等效多散射中心，由此得 到高频区的回波表达式如下公式8所示：
[0035]
[0036] 公式8中，K为目标场景中等效散射中心个数，y(xi，yi)表示场景中点（xi，yi)的散 射函数，f为雷达发射频率，Θ为雷达脉冲与场景y轴形成的夹角，c为光速，j为虚数单位；
[0037] a2.将单个散射中心的回波进行频率和方位的离散化采样，频率采样点数为他，方 位角采样点数为N2,将其表示为向量形式为：
[0038]
[0039] 其中，f为单个散射中心回波歹(/，的的向量形式
表示观测矩阵中的原子，ypyUuyO为目标场景中第i个散射中心的值；
[0040] a3.利用步骤a2中单个散射中心的回波和场景散射函数的投影关系，构建目标场 景中有K个散射中心时的聚束SAR成像模型：
[0041] g = Θ μ+η,
[0042] 其中，g为回波向量，μ是场景散射系数向量，Θ为观测矩阵，η表示噪声向量。
[0043] 本发明的有益效果为，本发明联合图像域的峰值区域分割与信号域的稀疏贝叶斯 学习进行特征提取，首先通过PRS算法分割强散射区域，以避免强散射点对SBL算法的影响， 然后利用SBL算法自主学习提取弱散射特征点，弥补了基于图像域的PRS算法对弱目标分割 的不足。该方法充分利用了PRS算法和SBL算法各自的优势，对弱目标特征提取表现出了很 好的效果，具有可调参数少，鲁棒性强等优点;并且在PRS算法分割强散射区域后，测量矩阵 的维数也相应地下降，减小了运算量，提高了运算效率，从而可以广泛应用于SAR目标特征 提取领域，为后续自动目标识别(ATR)提供基础。
【附图说明】
[0044] 图1是SAR几何结构；
[0045] 图2是本发明的流程图；
[0046] 图3是仿真试验中用不同方法得到的特征提取效果图。
【具体实施方式】
[0047] 下面结合附图，详细描述本发明的技术方案：
[0048] 本发明使用的SAR几何结构如图1所示，其中雷达沿着移动路径移动的过程中不断 向目标场景发射与接收脉冲，Θ为雷达脉冲与场景y轴形成的夹角。
[0049] 参照图2,本发明的具体实施步骤如下：
[0050] 步骤1，构建聚束SAR成像模型：
[0051] 1.1)根据电磁理论，当目标散射体的电尺寸远大于波长时，目标的高频电磁散射 特性可以由局部散射特性的合成来表示，这些局部散射通常被称为等效多散射中心，由此 得到高频区的回波表达式：
[0052]
[0053] 其中K为目标场景中等效散射中心个数，μ(Xi，yi)表示场景中点（Xi，yi)的散射函 数，f为雷达发射频率，Θ为雷达脉冲与场景y轴形成的夹角，c为光速，j为虚数单位；
[0054] 1.2)将单个散射中心的回波进行频率和方位的离散化采样，频率采样点数为仏， 方位角米样点数为N2,将其表不为向量形式：
[0055]
12345678 其中，i为单个散射中心回波或/，的的向量形式 表示观测矩阵中的原子，yi = μ(Xi，yi)为目标场景中第i个散射中心的值；
2 1.3)利用1.2)表达式中单个散射中心的回波和场景散射函数的投影关系，构建目 标场景中有K个散射中心时的聚束SAR成像模型： 3 g = Θ μ+η, 4 其中，g为回波向量，μ是场景散射系数向量，Θ为观测矩阵，η表示噪声向量。 5 步骤2，由极坐标成像算法得到目标场景的成像结果； 6 步骤3,根据所述步骤2得到的图像，利用基于图像域的PRS算法对强散射点进行特 征提取： 7 3.1)设置强点门限τι，利用分水岭算法对强散射区域进行初始分割，得到办个初始 分割区域； 8 3.2)设置合并门限τ2，对3.1)得到的初始分割区域进行合并，当任意两个相邻区 域边界上像素点的幅值低于合并门限时，认为这两个相邻区域是由噪声和纹理杂波引起的 过分割导致的，因此合并这两个相邻区域，最后得到R2个合并后的分割区域，R2彡R1;
[0064]步骤4,分别利用加权位置估计和最小二乘获取所述步骤3中得到的分割区域所对 应的强散射点位置和幅度估计，并将位置索引放入集合Ω中；
[0065]步骤5,获取残余散射区域回波：
[0066] 5.1)原始回波信号的矩阵形式为：
[0067] g= Θ (Ai+A2)y+n = gi+g2+n，其中，A2为区域选择矩阵，其形式为将单位矩阵I中与 集合Ω中的原子所对应的列置为零，Αι+Α2=Ι，由此得到维数较小的测量矩阵，gi= ? Αιμ为 强散射区域回波，g2= ?―为残余散射区域回波；
[0068] 5.2)将残余散射区域回波方程等价描述为：
[0069] g2=?QW，其中，Θ〇表示残余散射区域的测量矩阵， w为相应的散射函数；
[0070]步骤6,利用基于信号域的SBL算法对残余信号的弱散射点进行特征提取，等价于 求解如下问题：
[0071]
，其中，ε为残差；
[0072] 6.1)基于高斯先验假设可得，所述残余散射区域的散射函数w的分布为：
[0073] ja(w:;.V) ~ Ν(0, ν)，
[0074] 其中，为目标场景中残余散射区域等效散射中心个数，％为一个非 负的参数，控制着w的稀疏度，当Vi = 0时，表示Wi为零；
[0075] 6.2)假设噪声也满足高斯分布，结合6.1)可得，w的后验分布可以表示成均值为
.方差为 _ . " 的高斯分布，其中，β是噪声的方差；
[0076] 6.3)当超参数ν和β被估计出来后，w的点估计由其后验分布均值确定；
[0077] 6.4)超参数ν和β由第二类最大似然估计获得，即等效于计算以下代价函数的极 值：
[0078]
车中，I为单位矩阵；
[0079] 6.5)采用期望最大化算法，得到参数ν和邱勺学习规则：
[0080]
[0081]
[0082] 其中，Ν是残余散射区域回波向量g2的维数；
[0083] 步骤7,结合所提取的强散射点和弱散射点，得到SAR目标特征提取结果。
[0084] 本发明的效果可以通过下述仿真实验加以说明：
[0085] 1.仿真条件
[0086]运行环境为CPU:Inter Core i5-2430M，内存 10.0GB，32位Windows操作系统，软件 采用MATLAB R2012b，系统参数设置如表1所示：
[0087] 表1系统参数
[0088]
[0089] 2.仿真内容与结果
[0090] 在上述条件下，采用本发明对电磁仿真数据进行特征提取，结果如图3所示。底图 为利用PFA算法得到的成像结果。由图3 (a)和图3 (b)可以看出，SBL和PRS算法对旁瓣和杂波 很敏感，而PRS-SBL算法克服了 SBL算法和PRS算法各自的缺点，对杂波环境中的弱散射特征 点提取表现出了很好的效果，如图3(c)所示。表格2列出了不同信噪比下各算法特征提取结 果的性能指标，其中，正确的位置估计定义为位置估计偏差在两个分辨单元之内，目标标准 散射中心位置由几何光学方法获得;残差的计算公式为:llg-111 2/|41124表示重构的回波 向量。设置强点门限为15dB、合并门限为3dB以及稀疏贝叶斯学习次数为4。可以看出，与单 独利用PRS和SBL算法相比，联合利用两种算法得到的位置估计准确率和残差都要更好。 [0091 ]表2不同信噪比下各算法特征提取结果的性能指标
[0092]
【主权项】
1. 一种用于合成孔径雷达的目标特征提取方法，其特征在于，包括W下步骤： a. 构建聚束SAR的成像模型如下公式1所示： g= Θμ+η(公式 1) 公式1中，g为回波向量，μ是场景散射系数向量，Θ为观测矩阵，η表示噪声向量； b. 由极坐标成像算法得到目标场景的成像图像； C.根据获得的成像图像，对强散射点进行特征提取，具体方法为： cl.对获得的成像图像进行分割，具体为： cl 1.设置强散射点口限τι，利用分水岭算法对强散射区域进行初始分割，得到Ri个初始 分割区域； cl2.设置合并口限T2,对步骤cll得到的初始分割区域进行合并，当任意两个相邻区域 边界上像素点的幅值低于合并口限时，认为运两个相邻区域是由噪声和纹理杂波引起的过 分割导致的，因此合并运两个相邻区域，最后得至帖个合并后的分割区域，R2《化； c2.获取分割区域所对应的强散射点位置和幅度估计，并将位置索引放入集合Ω中； d. 获取成像图像残余散射区域回波，具体为： dl.设原始回波信号的矩阵形式如下公式2: g= Θ (Ai+A2)y+n=gi+g2+n(公式2) 公式2中，A2为区域选择矩阵，其形式为将单位矩阵I中与集合Ω中的原子所对应的列置 为零，Ai+A2 = I，由此得到维数较小的测量矩阵，gi=0Aw为强散射区域回波，g2=0A2y为 残余散射区域回波； d2.将残余散射区域回波方程等价描述为如下公式3: 邑2= 0OW(公式3) 公式3中，00表示残余散射区域的测量矩阵，W为相应的散射函数； e. 对残余信号的弱散射点进行特征提取，具体为： el.基于高斯先验假设可得，所述残余散射区域的散射函数W的分布如下公式4:公式4中，v = A'Ay〇',)二.Μ为目标场景中残余散射区域等效散射中屯、个数，VI为一个非 负的参数，控制着W的稀疏度，当Vi = 0时，表示Wi为零； e2.假设噪声也满足高斯分布，结合公式4可得，W的后验分布可W表示成均值为向高斯分布，其中，β是噪声的方差；当超参数V 和β被估计出来后，W的点估计由其后验分布均值确定;超参数V和β由第二类最大似然估计 获得，即等效于计算公式5所示函数的极值：公式5中，I为单位矩阵； f. 采用期望最大化算法，得到参数V的学习规则如下公式6和β的学习规则如下公式7: 其中，N是残余散射区域回波向量g2的维数； g.根本步骤C中获得的强散射点和步骤e-f中获得的弱散射点，得到SAR目标特征提取 结果。2.根据权利要求1所述的一种用于合成孔径雷达的目标特征提取方法，其特征在于，所 述步骤a的具体方法为： al.根据电磁理论，当目标散射体的电尺寸远大于波长时，目标的高频电磁散射特性可 W由局部散射特性的合成来表示，运些局部散射通常被称为等效多散射中屯、，由此得到高 频区的回波表达式如下公式8所示：公式8中，K为目标场景中等效散射中屯、个数，y(xi，y〇表示场景中点（xi，y〇的散射函 数，f为雷达发射频率，Θ为雷达脉冲与场景y轴形成的夹角，C为光速，j为虚数单位； a2.将单个散射中屯、的回波进行频率和方位的离散化采样，频率采样点数为化，方位角 采样点数为化，将其表示为向量形式为：其中，忌为单个散射中屯、回波巧/，巧的向量形式表示观测矩阵中的原子，μι = μ(χι，7ι)为目标场景中第i个散射中屯、的值； a3.利用步骤曰2中单个散射中屯、的回波和场景散射函数的投影关系，构建目标场景中 有K个散射中屯、时的聚束SAR成像模型： 邑=Θμ+η， 其中，g为回波向量，μ是场景散射系数向量，Θ为观测矩阵，η表示噪声向量。
【文档编号】G01S13/90GK105974412SQ201610397199
【公开日】2016年9月28日
【申请日】2016年6月7日
【发明人】杨悦, 万群, 丛迅超, 张庆, 龙柯宇, 邹麟, 殷吉昊
【申请人】电子科技大学