CN107255800B - 多极化雷达成像的双层次块稀疏匹配追踪方法及装置 - Google Patents

Abstract

本发明提供多极化雷达成像的双层次块稀疏匹配追踪方法及装置，方法包括：获取多个极化通道下雷达回波并向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K、参数δ、基信号集和各极化通道回波残余初始值；计算各极化通道的稀疏解并找出其中K个最大值对应的位置作为支撑集；对各极化下支撑集投票选出票数最高的支撑集组成各极化通道的共同支撑集；用概率图模型开发目标的双层次块稀疏性，获得新的共同支撑集；根据新的共同支撑集获取新的稀疏解，根据新的稀疏解获取各极化通道新的回波残余；判断新的回波残余是否满足预设迭代停止条件，若满足则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。本发明计算量较小，能保证图像质量。

Description

多极化雷达成像的双层次块稀疏匹配追踪方法及装置

技术领域

本发明涉及雷达成像技术领域，尤其涉及一种多极化雷达成像的双层次块稀疏匹配追踪方法及装置。

背景技术

为了获取高分辨率的雷达图像，合成孔径雷达利用雷达与目标的相对运动，是一种切实可行的雷达成像技术。微波能够穿透云团、障碍物等，并且不会受到光照条件的影响，合成孔径雷达能够对观测场景进行全天候的观测。

雷达有不同的极化发射或接收方式。一个目标，例如飞机或者舰船，对不同的极化波形有着不同的散射特性。采用多极化雷达观测能够提供更丰富的关于目标的信息，例如目标的位置、方向、大小等信息。

压缩感知是最近几年兴起的一种的信号处理理论。在压缩感知的框架下，可以采用少量的雷达观测数据，获得更高质量的雷达图像。在场景稀疏的高分辨雷达图像中，一个目标往往占据了多个相邻的像素，即成簇出现，这种图像特性可以称之为块稀疏性。针对多极化合成孔径雷达，考虑到各个极化通道下对目标场景的划分是一致的，则不同极化下的图像呈现联合稀疏性。联合稀疏性在某种程度上也可以称之为块稀疏性。因此，在多极化合成孔径雷达成像中，可以开发目标的成簇特性和联合稀疏特性，即双层次块稀疏性。开发雷达图像中双层次块稀疏性能够更好的提高成像质量，使得目标更加清晰，更好的抑制杂点。现有的针对双层次块稀疏性的压缩感知方法为贝叶斯类的方法。贝叶斯类的方法计算量较大，且存在大量的未知参数需要手动调节，不利于进行快速的雷达成像。

鉴于此，如何针对多极化合成孔径雷达成像中双层次块稀疏性，提供一种计算量较小且能够保证图像质量的压缩感知方法及装置成为目前需要解决的技术问题。

发明内容

为解决上述的技术问题，本发明提供一种多极化雷达成像的双层次块稀疏匹配追踪方法及装置，是针对多极化合成孔径雷达成像中双层次块稀疏性的压缩感知方法，计算量较小，能保证图像质量。

第一方面，本发明提供一种多极化雷达成像的双层次块稀疏匹配追踪方法，包括：

获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，设置基信号集和各极化通道的回波残余的初始值；

计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；

对各极化下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；

利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；

根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；

判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

可选地，所述根据所述新的共同支撑集，获取各极化通道的新的稀疏解，包括：

根据所述新的共同支撑集，利用最小二乘计算得到各极化通道的新的稀疏解。

可选地，所述双层次块稀疏性包括：成簇特性和联合稀疏性。

可选地，所述预设迭代停止条件，包括：

其中，Q为极化通道的通道总数量，q＝1,2,…,Q，r^(q,old)为根据所述新的稀疏解获取各极化通道新的回波残余之前第q个极化通道的回波残余，r^(q,new)为根据所述新的稀疏解所获取的第q个极化通道的新的回波残余。

可选地，在判断所述新的回波残余是否满足预设迭代停止条件之后，所述方法还包括：

若所述新的回波残余不满足预设迭代停止条件，则根据各极化通道的新的稀疏解和新的回波残余，返回执行所述计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集的步骤进行继续迭代。

第二方面，本发明提供一种多极化雷达成像的双层次块稀疏匹配追踪装置，包括：

第一获取模块，用于获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，设置基信号集和各极化通道的回波残余的初始值；

计算模块，用于计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；

投票模块，用于对各极化下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；

第二获取模块，用于利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；

第三获取模块，用于根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；

结果获取模块，用于判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

可选地，所述第三获取模块，具体用于

根据所述新的共同支撑集，利用最小二乘计算得到各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余。

可选地，所述双层次块稀疏性包括：成簇特性和联合稀疏性。

可选地，所述预设迭代停止条件，包括：

其中，Q为极化通道的通道总数量，q＝1,2,…,Q，r^(q,old)为根据所述新的稀疏解获取各极化通道新的回波残余之前第q个极化通道的回波残余，r^(q,new)为根据所述新的稀疏解所获取的第q个极化通道的新的回波残余。

可选地，所述装置还包括：

迭代模块，用于若所述新的回波残余不满足预设迭代停止条件，则根据各极化通道的新的稀疏解和新的回波残余，返回执行所述计算模块进行继续迭代。

由上述技术方案可知，本发明的多极化雷达成像的双层次块稀疏匹配追踪方法及装置，通过获取多个极化通道下雷达回波并向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K、参数δ、基信号集和各极化通道回波残余初始值；计算各极化通道的稀疏解并找出其中K个最大值对应的位置作为支撑集；对各极化下支撑集投票选出票数最高的支撑集组成各极化通道的共同支撑集；用概率图模型开发目标的双层次块稀疏性，获得新的共同支撑集；根据新的共同支撑集获取新的稀疏解，根据新的稀疏解获取各极化通道新的回波残余；判断新的回波残余是否满足预设迭代停止条件，若满足则停止迭代并根据当前迭代结果计算并融合得到最终成像结果，由此，能够实现针对多极化合成孔径雷达成像中双层次块稀疏性的压缩感知，计算量较小，能保证图像质量。

附图说明

图1为本发明一实施例提供的一种多极化雷达成像的双层次块稀疏匹配追踪方法的流程示意图；

图2为本发明实施例提供的一种多极化雷达成像的双层次块稀疏匹配追踪方法的实验设置图；

图3为本发明实施例提供的双层次块稀疏性的示例示意图；

图4为本发明实施例采用概率图模型的方式来对双层次块稀疏性进行建模后得到的概率图模型示意图；

图5(a)-图5(c)为图4建模后得到的概率图模型的一种概率密度估计示意图；

图6为本发明另一实施例提供的一种多极化雷达成像的双层次块稀疏匹配追踪方法的流程示意图；

图7为本发明一实施例提供的一种多极化雷达成像的双层次块稀疏匹配追踪方法的成像结果示意图；

图8为本发明一实施例提供的一种多极化雷达成像的双层次块稀疏匹配追踪装置的结构示意图；

图9为本发明实施例提供的一种电子设备的实体结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他的实施例，都属于本发明保护的范围。

图1示出了本发明一实施例提供的多极化雷达成像的双层次块稀疏匹配追踪方法的流程示意图，如图1所示，本实施例的多极化雷达成像的双层次块稀疏匹配追踪方法如下所述。

101、获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，设置基信号集和各极化通道的回波残余的初始值。

在具体应用中，可以参考图2所示的实验设置图，在多极化雷达成像中假定雷达所观测的场景V被分为N_x×N_y个像素点，N_x和N_y均为正整数，

表示像素点的位置，{1,…,q,…,Q}代表不同的极化通道，Q为极化通道的通道总数量，q＝1,2,…,Q。第q个极化通道下：合成孔径雷达的M^(q)个天线的位置分别是

L^(q)个频点分别为

则合成孔径雷达在第q个极化通道下，第m个天线、第l个频点处收到的雷达回波信号

的表达式为：

其中，

表示点

处的散射强度，c表示电磁波的传播速度，。在对雷达回波信号的表达式(1)进行向量化之后，可以得到下述公式(2)：

y^(q)＝Φ^(q)σ^(q)+w^(q)(2)

其中，w^(q)代表加性噪声，

n_x＝1,2,…,N_x,n_y＝1,2,…,N_y，矩阵Φ^(q)称为基信号矩阵，矩阵Φ^(q)的每一列为基信号，σ^(q)是

的向量化表示，σ^(q)中非零元素的个数是稀疏度K。

在具体应用中，在所述步骤101之前，还可以预先设定多极化合成孔径雷达参数得到雷达所观测的场景，例如可以参考表1，表1为一种多极化合成孔径雷达参数设置表。

表1

步进频率范围	1GHz
		中心频率	2.5GHz
频点间隔	5MHz
		雷达孔径长度	1.51m
天线个数	69个
		极化方式	HH、VV、HV

102、计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数。

在具体应用中，可以将第q个极化通道的支撑集定义为Λ^(q)，表示σ^(q)中非零元素的位置(K个)。

103、对各极化下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集。

在具体应用中，可以参考图3所示的双层次块稀疏性的示例，图3中填充斜线的区域1表示目标区域，白色区域2表示非目标区域。可以定义支撑集区域为s_Λ(q)＝1，非支撑集区域为

其中，s是一个N_xN_y×1的向量，

代表Λ^(q)的补集。考虑到各个极化通道之间的联合稀疏性，可以定义各极化通道的共同支撑集为Λ＝Λ^(q)(q＝1,2,…,Q)。因此，共同支撑集区域可定义为s_Λ＝1。

104、利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集。

其中，所述双层次块稀疏性包括：成簇特性和联合稀疏性。

在具体应用中，所述步骤104可以用马尔科夫随机场对目标的成簇特性进行描述，得到：

其中，N_i表示与像素i相邻的像素区域，像素i′∈N_i，s_i表示支撑集区域s中第i个元素，

表示支撑集区域s中像素i的邻域，

表示归一化的系数，这里考虑了二阶邻域***；

然后，可以采用概率图模型的方式来对双层次块稀疏性进行建模，建模后得到的模型可以参考图4，由于目标存在成簇特性，因此

影响着s_i；又由于联合稀疏性，因此s_i决定着

是否为0。所以，建模后得到的模型中各项的联合分布可以表示为：

其中，J(s_i)是关于s_i的函数。将J(s_i)最大化，可以得到：

对公式(5)求解，可以得到公式(5)的解为：

其中，

在具体应用中，可以根据图5(a)来定义

根据图5(b)来定义

根据图5(c)来定义

在图5(a)-5(c)中，横坐标

代表

的绝对值，纵坐标代表概率密度取值，α是

的最大值，τ为用于区别“大值”与“小值”的阈值，0＜ε₂＜ε₁，0＜ε₄＜ε₃；可以将上述公式(7)表示为：

其中，

可以假设

δ一般取一个很小的值，例如，本实施例可以取δ＝10^-6。所以，公式(8)可以表示为：

105、根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余。

在具体应用中，可以用r^(q,new)表示根据所述新的稀疏解所获取的第q个极化通道的新的回波残余，用r^(q,old)表示根据所述新的稀疏解获取各极化通道新的回波残余之前第q个极化通道的回波残余。

在具体应用中，所述步骤105可以根据所述新的共同支撑集，利用最小二乘计算得到各极化通道的新的稀疏解。

106、判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

在具体应用中，所述预设迭代停止条件可以包括：

其中，Q表示极化通道的通道总数量，q＝1,2,…,Q，r^(q,old)表示根据所述新的稀疏解获取各极化通道新的回波残余之前第q个极化通道的回波残余，r^(q,new)表示根据所述新的稀疏解所获取的第q个极化通道的新的回波残余。

在具体应用中，如图6所示，所述步骤106在判断所述新的回波残余是否满足预设迭代停止条件之后，还可以包括：

若所述新的回波残余不满足预设迭代停止条件，则根据各极化通道的新的稀疏解和新的回波残余，返回执行所述计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集的步骤进行继续迭代。

进一步地，为了方便描述本实施例的多极化雷达成像的双层次块稀疏匹配追踪方法，可以预先定义本方法中所使用到的三个函数：

1、

2、

3、Λ_vote＝majority-vote({Λ(^q)|q＝1,2,…,Q})：(a)设置

(b)对每一个q∈{1,2,…,Q}，n＝add（n,Λ^(q））；(c)Λ_vote＝max_ind(n,K)；

然后本实施例所述方法的具体步骤可以包括：

输入:{y^(q),Φ^(q)|q＝1,2,…,Q},K,δ。

初始化:针对q＝1,2,…,Q，设置r^(q,old)＝y^(q)；针对q＝1,2,…,Q，设置

迭代步骤:

S1:针对q＝1,2,…,Q，

S2:针对q＝1,2,…,Q，

S3:Λ_vote＝majority-vote({Λ(^q)|q＝1,2,…,Q}),然后设置

S4:对于每一个i∈{1,2,…,N_xN_y},利用上述公式(9)和

和

计算Δ；如果Δ＞0,

否则

在公式(9)中，τ可以先由步骤S1中得到的稀疏解的第2K大的幅度的值确定。

S5:

S6:针对q＝1,2,…,Q，

然后设置

中最小的(P_xP_y-K)个系数为0。

S7:针对q＝1,2,…,Q，

S8:如果

设置r^(q,old)＝r^(q,new),

返回步骤S1；否则停止迭代。

输出:

本实施例的多极化雷达成像的双层次块稀疏匹配追踪方法，应用于处理器，通过获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，设置基信号集和各极化通道的回波残余的初始值；计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；对各极化下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果；由此，能够实现针对多极化合成孔径雷达成像中双层次块稀疏性的压缩感知，计算量较小，能保证图像质量。

本实施例所述方法，通过多极化合成孔径雷达采集雷达回波数据，再采用本发明所提出的方法重建并融合目标场景，得到雷达的成像结果。本实施例所述方法着重开发了雷达图像中的目标的成簇特性和多极化下图像的联合稀疏性。由于本实施例所述方法采用了匹配追踪的框架，其计算量

也显著低于贝叶斯类方法的计算量

(其中T是贝叶斯类算法的迭代次数)。经过验证，采用本实施例所述方法进行多极化合成孔径雷达成像，可参考图7，能够使目标更加清晰，提高成像的视觉效果，保留弱目标的能力也得到增强，杂波得到抑制，从而有力的证明了本实施例所述方法的有效性。

图8示出了本发明一实施例提供的一种多极化雷达成像的双层次块稀疏匹配追踪装置的结构示意图，如图8所示，本实施例的多极化雷达成像的双层次块稀疏匹配追踪装置，包括：第一获取模块81、计算模块82、投票模块83、第二获取模块84、第三获取模块85和结果获取模块86；其中：

第一获取模块81，用于获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，设置基信号集和各极化通道的回波残余的初始值；

计算模块82，用于计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；

投票模块83，用于对各极化下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；

第二获取模块84，用于利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；

第三获取模块85，用于根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；

结果获取模块86，用于判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

其中，所述双层次块稀疏性包括：成簇特性和联合稀疏性。

在具体应用中，所述第三获取模块85，可具体用于

根据所述新的共同支撑集，利用最小二乘计算得到各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余。

在具体应用中，所述预设迭代停止条件，可以包括：

其中，Q为极化通道的通道总数量，q＝1,2,…,Q，r^（q,old）为根据所述新的稀疏解获取各极化通道新的回波残余之前第q个极化通道的回波残余，r^（q,new)为根据所述新的稀疏解所获取的第q个极化通道的新的回波残余。

在具体应用中，所述装置还可以包括图中未示出的：

迭代模块，用于若所述新的回波残余不满足预设迭代停止条件，则根据各极化通道的新的稀疏解和新的回波残余，返回执行所述计算模块进行继续迭代。

本实施例的多极化雷达成像的双层次块稀疏匹配追踪装置，可以用于执行前述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。

本实施例的多极化雷达成像的双层次块稀疏匹配追踪装置，应用于处理器中，能够实现针对多极化合成孔径雷达成像中双层次块稀疏性的压缩感知，计算量较小，能保证图像质量，能够使目标更加清晰，提高成像的视觉效果，保留弱目标的能力也得到增强，杂波得到抑制。

图9示出了本发明实施例提供的一种电子设备的实体结构示意图，如图9所示，该电子设备可以包括：处理器11、存储器12、总线13及存储在存储器12上并可在处理器11上运行的计算机程序；

其中，所述处理器11，存储器12通过所述总线13完成相互间的通信；

所述处理器11执行所述计算机程序时实现上述方法实施例所提供的方法，例如包括：获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，设置基信号集和各极化通道的回波残余的初始值；计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；对各极化下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

本发明实施例提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述方法实施例所提供的方法，例如包括：获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，设置基信号集和各极化通道的回波残余的初始值；计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；对各极化下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

本领域内的技术人员应明白，本申请的实施例可提供为方法、***、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置/***。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。术语“上”、“下”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。

本发明的说明书中，说明了大量具体细节。然而能够理解的是，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。类似地，应当理解，为了精简本发明公开并帮助理解各个发明方面中的一个或多个，在上面对本发明的示例性实施例的描述中，本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该公开的方法解释呈反映如下意图：即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如权利要求书所反映的那样，发明方面在于少于前面公开的单个实施例的所有特征。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明的单独实施例。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。本发明并不局限于任何单一的方面，也不局限于任何单一的实施例，也不局限于这些方面和/或实施例的任意组合和/或置换。而且，可以单独使用本发明的每个方面和/或实施例或者与一个或更多其他方面和/或其实施例结合使用。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。

Claims (10)

1.一种多极化雷达成像的双层次块稀疏匹配追踪方法，其特征在于，包括：

获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，0<δ<1，设置基信号集和各极化通道的回波残余的初始值；

计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；

对各极化通道下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；定义支撑集区域为

非支撑集区域为

其中，s是一个N_xN_y×1的向量，

代表Λ^(q)的补集；定义各极化通道的共同支撑集为Λ＝Λ^(q)(q＝1,2,…,Q)，共同支撑集区域定义为s_Λ＝1；

利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；其中，采用概率图模型的方式对双层次块稀疏性进行建模，建模后得到的模型中各项的联合分布表示为：

式中，

表示支撑集区域s中像素i的邻域，s_i表示支撑集区域s中第i个元素，{1,…,q,…,Q}代表不同的极化通道，Q为极化通道的通道总数量，

表示在第1个极化通道下第i个元素的散射强度，

表示在第q个极化通道下第i个元素的散射强度，

表示在第Q个极化通道下第i个元素的散射强度，J(s_i)是关于s_i的函数；将J(s_i)最大化，得到：

对公式(5)求解，得到公式(5)的解为：

其中，

N_i表示与像素i相邻的像素区域，像素i′∈N_i；

根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；

判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

2.根据权利要求1所述的方法，其特征在于，所述根据所述新的共同支撑集，获取各极化通道的新的稀疏解，包括：

根据所述新的共同支撑集，利用最小二乘计算得到各极化通道的新的稀疏解。

3.根据权利要求1所述的方法，其特征在于，所述双层次块稀疏性包括：成簇特性和联合稀疏性。

4.根据权利要求3所述的方法，其特征在于，所述预设迭代停止条件，包括：

其中，Q为极化通道的通道总数量，q＝1,2,…,Q，r^(q,old)为根据所述新的稀疏解获取各极化通道新的回波残余之前第q个极化通道的回波残余，r^(q,new)为根据所述新的稀疏解所获取的第q个极化通道的新的回波残余。

5.根据权利要求1-4中任一项所述的方法，其特征在于，在判断所述新的回波残余是否满足预设迭代停止条件之后，所述方法还包括：

若所述新的回波残余不满足预设迭代停止条件，则根据各极化通道的新的稀疏解和新的回波残余，返回执行所述计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集的步骤进行继续迭代。

6.一种多极化雷达成像的双层次块稀疏匹配追踪装置，其特征在于，包括：

第一获取模块，用于获取多个极化通道下的雷达回波，对所述雷达回波进行向量化处理，根据向量化的雷达回波计算基信号矩阵，设置稀疏度K和参数δ，0<δ<1，设置基信号集和各极化通道的回波残余的初始值；

计算模块，用于计算各极化通道的稀疏解，找出所有稀疏解中K个最大值所对应的位置作为支撑集，K为正整数；

投票模块，用于对各极化通道下的支撑集进行投票，选出票数最高的支撑集组成各极化通道的共同支撑集；定义支撑集区域为

非支撑集区域为

其中，s是一个N_xN_y×1的向量，

代表Λ^(q)的补集；定义各极化通道的共同支撑集为Λ＝Λ^(q)(q＝1,2,…,Q)，共同支撑集区域定义为s_Λ＝1；

第二获取模块，用于利用概率图模型，开发目标的双层次块稀疏性，获得新的共同支撑集；其中，采用概率图模型的方式对双层次块稀疏性进行建模，建模后得到的模型中各项的联合分布表示为：

式中，

表示支撑集区域s中像素i的邻域，s_i表示支撑集区域s中第i个元素，{1,…,q,…,Q}代表不同的极化通道，Q为极化通道的通道总数量，

表示在第1个极化通道下第i个元素的散射强度，

表示在第q个极化通道下第i个元素的散射强度，

表示在第Q个极化通道下第i个元素的散射强度，J(s_i)是关于s_i的函数；将J(s_i)最大化，得到：

对公式(5)求解，得到公式(5)的解为：

其中，

N_i表示与像素i相邻的像素区域，像素i′∈N_i；

第三获取模块，用于根据所述新的共同支撑集，获取各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余；

结果获取模块，用于判断所述新的回波残余是否满足预设迭代停止条件，若满足，则停止迭代并根据当前迭代结果计算并融合得到最终成像结果。

7.根据权利要求6所述的装置，其特征在于，所述第三获取模块，具体用于

根据所述新的共同支撑集，利用最小二乘计算得到各极化通道的新的稀疏解，以及根据所述新的稀疏解，获取各极化通道的新的回波残余。

8.根据权利要求6所述的装置，其特征在于，所述双层次块稀疏性包括：成簇特性和联合稀疏性。

9.根据权利要求8所述的装置，其特征在于，所述预设迭代停止条件，包括：

其中，Q为极化通道的通道总数量，q＝1,2,…,Q，r^(q,old)为根据所述新的稀疏解获取各极化通道新的回波残余之前第q个极化通道的回波残余，r^(q,new)为根据所述新的稀疏解所获取的第q个极化通道的新的回波残余。

10.根据权利要求6-9中任一项所述的装置，其特征在于，所述装置还包括：

迭代模块，用于若所述新的回波残余不满足预设迭代停止条件，则根据各极化通道的新的稀疏解和新的回波残余，返回执行所述计算模块进行继续迭代。

Images