CN111175747A - 一种基于多通道复图像空间的相位误差估计方法 - Google Patents

Abstract

本发明属于通信技术领域，具体涉及一种基于多通道复图像空间的相位误差估计方法，包括：获取回波数据；对回波数据进行置零操作，得到置零回波数据；对置零回波数据进行方位频谱重构得到信号频谱；通过线频调变标算法对信号频谱进行成像得到第一图像集；对第一图像集进行低通滤波得到第二图像集；对第二图像集进行相位误差估计得到相位误差；根据相位误差对第一图像集进行误差补偿得到目标图像。本发明对多通道相位误差的估计是在图像空间中进行，避免多次迭代处理造成的繁重的计算量，提高了多通道相位误差估计的速度和效率；先进行重构和成像，再进行多通道相位误差的估计，由于成像后信号的能量相对于成像前更能集中，能得到更好的误差估计效果。

Description

一种基于多通道复图像空间的相位误差估计方法

技术领域

本发明属于通信技术领域，具体涉及一种基于多通道复图像空间的相位误差估计方 法。

背景技术

由于星载SAR具有免受国界区域和自然条件等方面限制，以及SAR全天时、全天候工作的特点，星载SAR成为近年来世界各地学者的热点研究对象。针对星载SAR先后发 展出多种成像算法，并取得了较好的成像结果。随着对SAR成像需求的提高，高方位分 辨率和宽测绘带成为近年来发展的重点目标。然而，由于低PRF和高分辨率之间是一组 基本矛盾，高方位分辨率和宽测绘带在传统的SAR系统中不能同时实现。

能够获得高方位分辨率和宽测绘带成像的方位多通道SAR系统吸引了越来越多的关 注，而方位多通道SAR系统的提出，有效地解决了低PRF和高分辨率之间的这一基本矛盾。该系统运行时PRF选用低于单通道奈奎斯特采样率的频率，通过方位信号重构消除 多通道回波在方位上存在多普勒模糊，因而常规的单通道的成像算法得以应用于多通道 数据的成像。

实际操作中，由于非理想因素存在于方位多通道SAR系统中，通道间存在相位增益误差，位置不确定和时间不确定，所以通道间不可避免地存在通道不匹配。这些通道间 的不匹配将会严重影响重构的性能，造成成像结果中存在方位模糊，严重降低成像质量。 因此，针对多通道SAR系统中通道间不匹配的估计和校正成为在实际操作中的关键问题。 针对方位多通道SAR系统通道间的相位误差，现有方法大多是在频域上进行估计补偿， 也有部分方法是在时域上采用最小熵的方法，通过迭代对通道间的相位误差进行迭代估 计。

针对方位多通道SAR系统通道间的不匹配，不少学者对这一问题进行了深入地研究， 其中西安电子科技大学李真芳教授针对分布式小卫星系统，提出的基于杂波的通道不匹 配的估计和校正的正交子空间法，这一方法已经成功应用于方位多通道SAR系统。其他研究者也提出了基于相邻通道相关特性的时域通道误差校正方法。上述研究所提出的方法是在图像成像前对通道误差进行估计，可将其归为一类信号空间方法。有部分研究者 提出在图像空间上采用最小熵的方法，对通道间的相位误差进行估计的方法，先后发展 出最小熵自聚焦方法、加权最小熵方法和基于最大似然估计的加权最小熵方法。由于无 法获得最小熵方法的自闭解，上述方法需要通过多次迭代来实现。

发明内容

为了解决现有技术中存在的上述问题，本发明提供了一种基于多通道复图像空间的 相位误差估计方法。本发明要解决的技术问题通过以下技术方案实现：

一种基于多通道复图像空间的相位误差估计方法，包括：

获取回波数据；

对所述回波数据进行置零操作，得到置零回波数据；

对所述置零回波数据进行方位频谱重构得到信号频谱；

通过线频调变标算法对所述信号频谱进行成像得到第一图像集；

对所述第一图像集进行低通滤波得到第二图像集；

对所述第二图像集进行相位误差估计得到相位误差；

根据所述相位误差对所述第一图像集进行误差补偿得到目标图像。

在本发明的一个实施例中，所述置零操作为：

取含有M个通道的所述回波数据中的第m通道数据，再将其余通道数据置零，每个通道数据只进行一次置零，经过M次置零后得到M组置零数据，其中m≤M。

在本发明的一个实施例中，对所述置零回波数据进行方位频谱重构得到信号频谱， 包括：

对所述M组置零数据进行逆滤波得到M组逆滤波信号；

对所述M组逆滤波信号进行重新排列得到信号频谱。

在本发明的一个实施例中，对所述第二图像集进行相位误差估计得到相位误差，包 括：

对所述第二图像集中的基频分量成像结果的协方差矩阵进行分析得到图像空间；

通过最小化代价函数对所述图像空间进行误差估计得到相位误差。

在本发明的一个实施例中，所述最小代价函数的表达式为：

其中，ψ为代价函数，Γ为相位误差，η＝[Γ₁₁,Γ₂₂,…Γ_MM]^T，Γ_MM为矩阵Γ的 第m行m列元素，Vⁱ＝diag P_Thⁱ，diag P_Thⁱ为一个对角元素为向量P_Thⁱ的对角 矩阵，L为正整数，hⁱ为H的第i列，H为导向矢量矩阵，pt为矩阵P的第t行，i＞0。

在本发明的一个实施例中，所述代价函数为：

其中，

为B的估计量。

本发明的有益效果：

本发明对多通道相位误差的估计是在图像空间中进行，且避免了多次迭代处理造成 的繁重的计算量，提高了多通道相位误差估计的速度和效率；还采用基于基频滤波的方式，将重构后的信号通过低通滤波器，降低了多通道信号在图像空间中的自由度，增加 了本发明的相位误差估计方法所适用的场景；先进行重构和成像，再进行多通道相位误 差的估计，由于成像后信号的能量相对于成像前更能集中，所以当选取信噪比高的子图 像进行通道相位误差估计时，能得到更好的误差估计效果。

以下将结合附图及实施例对本发明做进一步详细说明。

附图说明

图1是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法的流程 示意图；

图2是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法的流程 框图；

图3是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法未通过 低通滤波器进行自由度压缩的成像结果；

图4是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法通过低 通滤波器进行自由度压缩后的成像结果；

图5是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法估计得 到的通道相位误差补偿给图像后的成像结果；

图6是本发明实施例提供的使用传统的信号空间方法成像且未经过通道相位误差补 偿的成像结果。

具体实施方式

下面结合具体实施例对本发明做进一步详细的描述，但本发明的实施方式不限于此。

请参见图1，图1是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法的流程示意图，包括：

获取回波数据；

对所述回波数据进行置零操作，得到置零回波数据；

对所述置零回波数据进行方位频谱重构得到信号频谱；

通过线频调变标算法对所述信号频谱进行成像得到第一图像集；

对所述第一图像集进行低通滤波得到第二图像集；

对所述第二图像集进行相位误差估计得到相位误差；

根据所述相位误差对所述第一图像集进行误差补偿得到目标图像。

具体的，取原始多通道信号的第一个通道的数据，将余下的通道置零，对其余的M-1 通道进行相同的重复的操作，共得到M组置零的数据，由此得到置零回波数据；将方位 多通道SAR系统对原始多通道信号的采样过程看成是一个滤波过程，所以对多通道数据 进行逆滤波，然后将逆滤波后得到的频谱分量进行重新排列，即可得到完整的信号频谱； 对M组置零后的数据进行信号重构，得到M组重构完的数据，由于每个通道的数据都不 是完整的数据，所以M组重构完的数据仍存在方位模糊；对M组重构后的数据通过线频 调变标算法进行成像，可得到M组存在方位模糊或重影的图像I_e,I_e＝[I₁，I₂,…,I_M]^T， I_M表示第_m组的图像,[·]^T为转置操作；成像结果分为两部分,一部分用于最终的成 像，另一部分用于通道相位误差估计；将匹配滤波得到的M幅有方位模糊或重影的图像 I_e通过低通滤波器，将高频分量滤除后即可降低图像空间相位误差估计时的自由度，通 过低通滤波器滤波得到的I_e,subT可用于通道相位误差的估计,I_e,subT为每个通道的基频 分量的成像结果,I_{e, subT}F_lowpassI_e，

F_lowpass表示低通滤波器；对基频 分量的成像结果I_e,subT的协方差矩阵进行特征分解，得到特征值，并根据特征值的大小 得到图像的信号空间和图像的噪声空间，最后定义代价函数ψ，在线性约束条件下，通 过最小化代价函数ψ来估计通道相位误差Γ，并将误差补偿给最终的成像。

本发明对多通道相位误差的估计是在图像空间中进行，且避免了多次迭代处理造成 的繁重的计算量，提高了多通道相位误差估计的速度和效率；还采用基于基频滤波的方式，将重构后的信号通过低通滤波器，降低了多通道信号在图像空间中的自由度，增加 了本发明的相位误差估计方法所适用的场景；先进行重构和成像，再进行多通道相位误 差的估计，由于成像后信号的能量相对于成像前更能集中，所以当选取信噪比高的子图 像进行通道相位误差估计时，能得到更好的误差估计效果。

在本发明的一个实施例中，所述置零操作为：

取含有M个通道的所述回波数据中的第m通道数据，再将其余通道数据置零，每个通道数据只进行一次置零，经过M次置零后得到M组置零数据，其中m≤M。

在本发明的一个实施例中，对所述置零回波数据进行方位频谱重构得到信号频谱， 包括：

对所述M组置零数据进行逆滤波得到M组逆滤波信号；

对所述M组逆滤波信号进行重新排列得到信号频谱。

在本发明的一个实施例中，对所述第二图像集进行相位误差估计得到相位误差，包 括：

对所述第二图像集中的基频分量成像结果的协方差矩阵进行分析得到图像空间；

通过最小化代价函数对所述图像空间进行误差估计得到相位误差。

在本发明的一个实施例中，所述最小代价函数的表达式为：

其中，ψ为代价函数，Γ为相位误差，η＝[Γ₁₁,Γ₂₂,…Γ_MM]^T，Γ_MM为矩阵Γ的 第m行m列元素，Vⁱ＝diag P_Thⁱ，diag P_Thⁱ为一个对角元素为向量P_Thⁱ的对角 矩阵，L为正整数，hⁱ为H的第i列，H为导向矢量矩阵，pt为矩阵P的第t行，i＞0。

在本发明的一个实施例中，所述代价函数为：

其中，

为B的估计量。

对于通道数为M，方位模糊数为N的SAR系统，对完整信号重构的原理如下：

有方位模糊的信号可表示为：S_amb f_aHS f_a

H为导向矢量矩阵，H的表达式为：

H＝[κ_-L,κ_-L+1,…,0，…,κ_L-1，κ_L]

其中，L为正整数，2L+1＝N，N为SAR系统的方位模糊数。

κ_i的表达式为：

其中，exp·为指数函数，j为虚数单位，f_a为方位频率，-f_p/2f_a≤f_p/2， i＝[-L，-L+1，…，0，…,L-1,L],L为正整数，2L+1＝N，N为SAR系统的方位模 糊数，f_p为脉冲重复频率，f_p的取值范围为B/Nf_p≤B，B为信号带宽，v为雷 达平台运动的速度，X_M为第m个通道的相位中心的位置，

m＝1,2,…,M-1,M,M为SAR系统的通道数,d为天线两个相邻通道之间的距离， [·]^T为转置操作。

S f_a为无方位模糊的信号，S f_a的表达式为：

S f_a＝[S(f_a-L·f_p) … S f_p+L·f_p]^T

其中，f_a为方位频率，-f_p/2f_a≤f_p/2，L为正整数，2L1＝N，N为SAR 系统的方位模糊数，f_p为脉冲重复频率，[·]^T为转置操作。

S_amb f_a的表达式为：

令P＝H^-1，H^-1为H的逆矩阵，H为导向矢量矩阵，P左乘以S_amb f_a可得到 S f_a＝PS_ambf_a，所以将有方位模糊的信号S_amb f_a通过逆滤波器P后，可得到无方 位模糊的信号S f_a。

进一步地说明本发明的方法步骤，请参见图2，图2是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法的流程框图：

1.取数据和置零

1a)取原始多通道信号的第一个通道的数据，将余下的通道置零，对于第m个通道的数据，置零后的结果为S_cm f_a＝[0,…,S_cm(f_a)，…，0]^T，[·]^T为装置操作；

1b)对其余的M-1个通道进行相同的重复的操作，共得到M组置零的数据，置零后每组数据的大小与原始的M个通道相同；

1c)得到的数据可用于通道相位误差的估计和最终的成像。

2.方位信号重构

本步骤是针对M组置零后的数据进行信号重构，可得到M组重构完的数据，由于每个通道的数据都不是完整的数据，所以M组重构完的数据仍存在方位模糊。

2a)分别对M个通道的数据进行逆滤波，对第m个通道的数据进行逆滤波的表达式为：

S_rm f_a＝PS_cm f_a

其中，P＝H^-1，H^-1为H的逆矩阵，H为导向矢量矩阵， S_rm f_a＝[S_rm(f_a-L·f_p) … S_rmf_a+L·f_p]^T，f_a为 方位频率，-fp/2≤fa≤fp/2， L为正整数，2L1＝N，N为SAR系统的方位模糊数，f_p为脉冲重复频率，[·]^T为 转置操作；

2b)对第m个通道滤波得到的信号频谱S_rm f_a进行重新排列，可得到完整的信号 频谱

f_p/2≤f_a≤f_p/2，

3.多通道数据成像

将得到完整信号频谱的多通道数据通过线频调变标(CS,Chirp Scaling)算法进行成像，在图像空间中可得到M幅图像。

对方位信号成像可等效于对方位信号进行匹配滤波，对步骤2中重构得到的完整的 信号频谱

进行匹配滤波，将第m个通道信号匹配滤波后的结果表示为I_Mx_n， 上述过程的表达式为:

其中，x_n为零时刻SAR平台在X轴的初始位置，

为完整的信号频谱， exp·为指数函数，j为虚数单位，v为雷达平台运动的速度，R_B为场景中心到雷达运 动轨迹的最近距离，f_am为最大多普勒频率，

λ为信号波长，

将匹配滤波函数

和傅里叶变换因子合成的新函数表示为

的表达式为：

其中，匹配滤波函数表达式为：

所有通道信号匹配滤波后的表达式为：

其中，I_M x_n为第m个通道信号匹配滤波后的结果,x_n为零时刻SAR平台在X轴 的初始位置；

为对第m个通道滤波得到的信号频谱进行重新排列后得到的完整 的离散化信号频谱，

为[-L·f_p/2,L·f_p/2]频率范围内均匀采样nan个点组成的列 向量,nan为方位向的采样点数，L为正整数，2L+1＝N，N为SAR系统的方位模糊 数，f_p为脉冲重复频率；

为函数

的一种离散化形式；

进一步可将所有通道信号匹配滤波后的结果表示为I，I的表达式为：

其中，L为正整数，2L+1＝N，N为SAR系统的方位模糊数,diag{p_i+L+1}表 示为一个对角元素为向量p_i+L+1的对角矩阵，p_i+L+1是矩阵P的第i+L+1行， P＝H^-1，H^-1为H的逆矩阵，H为导向矢量矩阵，i＝[-L,-L+1,…L-1，L]；S为无 方位模糊信号S(_fa)的离散化形式，S＝[S(f_a-L·f_p) … S(f_a+L·f_p)]^T，f_a为 [-f_p/2,f_p/2]频率范围内均匀采样nan个点组成的列向量,nan为方位向的采样点数， f_p为脉冲重复频率；M(f_a+i·f_p)为函数

的另一种离散化形式。

所以各通道的信号匹配滤波后，可将M幅有方位模糊或重影的图像表示为：

I_e＝ΓI

＝ΓI_sub1+ΓI_sub2+…+ΓI_subN

其中，

为第一个通道与其他通道之间的相位误差，

ΓI_{sub N}表示图像I_e的第n个分量，I_{sub N}＝P_NHSM_N， P_N＝diag{p_n}，diag{p_n}表示为一个对角元素为向量p_n的对角矩阵，p_n为矩阵P的 第n行，p_n＝[P_n1,…P_nM]，P_nM为矩阵P第n行第m列的元素，P＝H^-1，H^-1为H的 逆矩阵，H为导向矢量矩阵；S为无方位模糊信号S(f_a)的离散化形式，S＝[S(f_a-L·f_p) … S(f_a+L·f_p)]^T，f_a为[-f_p/2,f_p/2]频率范围内均匀采样nan 个点组成的列向量,nan为方位向的采样点数，f_p为脉冲重复频率，L为正整数， 2L+1＝N，N为SAR系统的方位模糊数，f_p为脉冲重复频率；M_N为

的第n 个分量，n＝1,2,…,N-1,N,

为函数

的一种离散化形式。

将成像结果分为两部分,一部分用于最终的成像，另一部分用于通道相位误差估计。 其目的在于使最终的成像结果既能得到通道相位误差补偿，又能保留高频分量，能得到 更清晰的图像。

4.自由度压缩

由于重构时需要将多个频谱进行重排，而各个频谱分量之间不能完全对齐，对方位 信号成像等效于一次匹配滤波，所以在通道数为M，方位模糊数为N的SAR系统中，成 像后图像空间中各个通道信号的自由度增加为2N-1。将图像I_e通过低通滤波器，将高 频分量滤除后即可降低图像空间相位误差估计时的自由度，通过低通滤波器滤波得到的 I_e,subT可用于通道相位误差的估计，这一过程的表达式为:

I_e,subT＝F_lowpassI_e

＝ΓP_THSM_T

其中，I_e,subT为每个通道的基频分量的成像结果，

N为SAR系统的方 位模糊数；F_lowpass表示低通滤波器；

为第一个通道与其他通道之间的相 位误差，

P_T＝diag p_t，diag p_t表示为一个对角元素为向量p_t的对角矩阵，pt为矩阵P的第t行，p_t＝[P_t1,…P_tM],P_tM为矩阵P第t行第m列的元 素，P＝H^-1，H^-1为H的逆矩阵，H为导向矢量矩阵；S为无方位模糊信号S f_a的离 散化形式，S＝[S(f_a-L·f_p) … S(f_a+L·f_p)]^T，f_a为[-f_p/2,f_p/2]频率范围内均 匀采样nan个点组成的列向量,nan为方位向的采样点数，L为正整数，2L+1＝N，N 为SAR系统的方位模糊数,f_p为脉冲重复频率；M_T为

的第t个分量，

为函 数

的一种离散化形式。

5.通道相位误差估计

5a)对步骤4中得到的基频分量的成像结果I_e,subT的协方差矩阵R I_e,subT进行分析，分析过程的表达式为：

其中，E·为求均值操作；

为第一个通道与其他通道之间的相位误差，

I_subT＝P_THSM_T，P_T＝diag p_t，diag p_t表示为一个对 角元素为向量pt的对角矩阵，pt为矩阵P的第t行，p_t＝[P_t1,…P_tM］,P_tM为矩阵P第 t行第m列的元素，P＝H^-1，H^-1为H的逆矩阵，H为导向矢量矩阵；Q＝SM_T，S为 无方位模糊信号S f_a的离散化形式，S＝[S(f_a-L·f_p) … S(f_a+Lf_p)^T，f_a为

[-f_p/2,fp/2]频率范围内均匀采样nan个点组成的列向量,nan为方位向的采样点数， L为正整数，2L+1＝N，N为SAR系统的方位模糊数,f_p为脉冲重复频率；M_T为

的第t个分量，

为函数

的一种离散化形式。将RQ的估计量表示

τ为时间，K为用于估计协 方差矩阵的服从独立同分布采样点的个数，2K+1≥2M+1,k0,1,2…2K，·^H为 对矩阵进行共轭转置操作。

5b)通过分解

的特征值可得到图像空间。

根据特征值的大小，可将得到的图像空间分为图像的信号空间和图像的噪声空间， 其中大的特征值对应的特征向量组成图像的信号空间，小的特征值对应的特征向量组成 图像的噪声空间，理论上信号空间与噪声空间正交。

特征值的分解满足如下两个条件：

1.特征值λ_m满足λ₁＞λ₂＞…＞λ_2L+1＞＞λ_2L+2＝…λ_M，特征值λ_m对应的特征向 量为b_m，其中m=1,2,…,M；

2.矩阵ΓP_TH的每列正交于矩阵B的每列，B＝[b_2L+2,b_2L+3，…，b_M]，B为特 征值λ_2L+2,λ_2L+3,…,λ_M对应的特征向量b_2L+2,b_2L+3,...,b_M组成的矩阵。

可将代价函数定义为:

其中,

为B的估计量,hⁱ为H的第i列，H为导向矢量矩阵。

5c)通过最小化代价函数ψ来估计通道相位误差Γ,当

等于B，Γ和P_T以及hⁱ取得真实值时，代价函数ψ达到最小值，这一过程的表达式为：

其中_，η＝[Γ₁₁,Γ₂₂,…Γ_MM]^T，Γ_MM表示矩阵Γ的第m行第m列的元素，

为第一个通道与其他通道之间的相位误差，

Vⁱ＝diag P_Thⁱ，diag P_Thⁱ表示为一个对角元素为向量P_Thⁱ的对角矩阵， P_T＝diag p_t，diag p_t表示为一个对角元素为向量p_t的对角矩阵，p_t为矩阵P的第 t行，p_t＝P_t1,…P_tM］，P_tM为矩阵P第t行第m列的元素，P＝H^-1，H^-1为H的逆矩 阵，H为导向矢量矩阵，hⁱ为H的第i列。

在线性约束条件η^HD＝1，D＝[1,0,0,…,0]^T下，求解代价函数ψ的最小值问题的 最优解为

其中，

可得到通道间的相位 误差估计值为

6.误差补偿

将估计得到的通道相位误差

补偿给步骤3中得到的M幅图像,再将补偿后的M幅图像进行叠加，最终可得到无方位模糊或重影的成像结果。

本发明的效果可以通过下述仿真实验加以说明：

将接收的实测数据在频域补零和时域采样处理后，进行仿真实验处理。多通道天线 阵列沿航迹方向均匀分布，发射通道位于天线阵列中心，天线阵列的任一通道都可作为接收回波信号的通道。相关参数如表1所示：

表1仿真参数表

参数	参数值	单位
			斜距	850	Km
载频	5.6	GHz
			斜视角	28.34	°
飞行平台速度	7024	m/s
			发射信号带宽	100	MHz
脉冲重复频率	2930	Hz

为了评估本发明提出的基于多通道复图像空间的相位误差估计方法的估计性能，在 存在相同通道相位误差的条件下，对比本发明提出的通道相位误差估计方法和在信号空 间的通道相位误差估计方法，仿真结果如表2所示：

表2两种通道相位误差估计方法的仿真结果

对比分析上述表格中的两种通道相位误差估计方法可知，使用本发明提出的基于多 通道复图像空间的相位误差估计方法，估计得到的通道相位误差更接近实际通道相位误 差。统计两种通道相位误差估计方法的方差，得到基于多通道复图像空间的相位误差估计方法的方差为0.0002，信号空间的通道相位误差估计方法的方差为0.9361，对比可 知本发明中提出的通道相位误差估计方法的方差更小。

经过上述分析，基于多通道复图像空间的相位误差估计方法的误差估计性能优于传 统的信号空间的通道相位误差估计方法。

为了验证多通道信号在信号重构后通过低通滤波器进行自由度压缩的有效性，请参 见图3和图4，图3是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法未通过低通滤波器进行自由度压缩的成像结果，图4是本发明实施例提供的一种基 于多通道复图像空间的相位误差估计方法通过低通滤波器进行自由度压缩后的成像结 果，对比分析图3和图4。在通道数为5，方位模糊数为3的SAR系统中，由于多通道 信号的信号重构需要将多个频谱进行重排，而不同的频谱分量之间不能完全对齐，所以 对于方位模糊数为N的信号，重构后每个通道信号的自由度为2N-1。当重构后的多通 道信号未通过低通滤波器进行自由度压缩时，成像结果如图3所示，对于方位模糊数为 3的方位信号，成像后图像中会存在5个模糊点；当重构后的多通道信号通过低通滤波 器进行自由度压缩时，成像结果如图4所示，成像后图像中模糊点数从5减少为3。对 比可知，将信号重构后的多通道信号，通过低通滤波器即可有效降低自由度。

请参见图5和图6，图5是本发明实施例提供的一种基于多通道复图像空间的相位误差估计方法估计得到的通道相位误差补偿给图像后的成像结果，图6是本发明实施例 提供的使用传统的信号空间方法成像且未经过通道相位误差补偿的成像结果，将本发明 提出的基于多通道复图像空间的相位误差估计方法估计得到的通道相位误差，补偿给图 像后的成像结果如图5所示，使用信号空间方法成像且未经通道相位误差补偿的成像结 果如图6所示，对比观察图5和图6发现，经过通道相位误差补偿的成像结果更加清晰。 进一步说明通过本发明提出的通道相位误差估计方法估计得到的相位误差，将其补偿给 图像后，能够提高成像质量，获得清晰度更高的成像结果。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本 发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明 的保护范围。

Claims (6)

1.基于多通道复图像空间的相位误差估计方法，其特征在于，包括：

获取回波数据；

对所述回波数据进行置零操作，得到置零回波数据；

对所述置零回波数据进行方位频谱重构得到信号频谱；

通过线频调变标算法对所述信号频谱进行成像得到第一图像集；

对所述第一图像集进行低通滤波得到第二图像集；

对所述第二图像集进行相位误差估计得到相位误差；

根据所述相位误差对所述第一图像集进行误差补偿得到目标图像。

2.根据权利要求1所述的基于多通道复图像空间的相位误差估计方法，其特征在于，所述置零操作为：

取含有M个通道的所述回波数据中的第m通道数据，再将其余通道数据置零，每个通道数据只进行一次置零，经过M次置零后得到M组置零数据，其中m≤M。

3.根据权利要求2所述的基于多通道复图像空间的相位误差估计方法，其特征在于，对所述置零回波数据进行方位频谱重构得到信号频谱，包括：

对所述M组置零数据进行逆滤波得到M组逆滤波信号；

对所述M组逆滤波信号进行重新排列得到信号频谱。

4.根据权利要求1所述的基于多通道复图像空间的相位误差估计方法，其特征在于，对所述第二图像集进行相位误差估计得到相位误差，包括：

对所述第二图像集中的基频分量成像结果的协方差矩阵进行分析得到图像空间；

通过最小化代价函数对所述图像空间进行误差估计得到相位误差。

5.根据权利要求4所述的基于多通道复图像空间的相位误差估计方法，其特征在于，所述最小代价函数的表达式为：

其中，Ψ为代价函数，Γ为相位误差，η＝[Γ₁₁,Γ₂₂,…Γ_MM]^T，Γ_MM为矩阵Γ的第m行m列元素，Vⁱ＝diag P_Thⁱ，diag P_Thⁱ为一个对角元素为向量P_Thⁱ的对角矩阵，L为正整数，hⁱ为H的第i列，H为导向矢量矩阵，p_t为矩阵P的第t行，i＞0。

6.根据权利要求5所述的基于多通道复图像空间的相位误差估计方法，其特征在于，所述代价函数为：

其中，

为B的估计量。

Images