本发明涉及通信
技术领域：
，具体地，涉及一种方位多通道sar系统通道间相位偏差校正方法和系统。尤其地，涉及一种基于子孔径图像最大对比度的方位多通道高分辨率宽幅(hrws)合成孔径雷达(sar)系统通道间相位偏差校正方法。
背景技术：
：星载方位多通道sar体制可有效克服传统单通道sar面临的“最小天线面积”约束，实现高分辨率宽幅sar成像，目前已广泛应用于国内外星载sar系统中，在大面积国土普查、农作物监测等领域得到了广泛的应用。多普勒谱重构是方位多通道sar成像处理的关键步骤，而通道间幅相误差将造成多普勒谱重构性能下降，进而造成系统成像性能下降。为有效校正通道间幅相误差，通常星载sar系统通过设计内定标链路对系统收发通道的幅相特性进行标定。但这里需要指出的是，天线安装误差、在轨天线热形变等将导致各通道天线相位中心位置存在测量误差。其中，垂直航向位置测量误差可等效为相位误差。因此，在方位多通道hrwssar成像处理前还需对通道间幅相误差进行自适应校正。在方位多通道hrwssar系统中，如何实现通道相位偏差的校正是研制人员必须考虑、无法回避的问题。目前，该问题已成为国内外研究的热点问题。本发明与国内外现有技术对比如下：王志斌等人，phasebiasestimationformulti-channelhrwssarbasedondopplerspectrumoptimisation,ieeelectronicletters,2016，52(21):1805-1807.该文提出了基于多普勒频谱最大对比度的相位偏差估计算法，但多普勒频谱对比度对通道间相位偏差不敏感，这导致该算法估计精度相对较差，而本发明给出的通道间相位偏差校正算法采用子孔径图像对比度作为优化准则，对通道间相位偏差更加敏感，可取得更高的估计精度。张双喜等人，asedonminimum-entropyhrwssarchannel-calibrationmethod.ietinternationalradarconference,2013:1-5.该文提出了利用子孔径图像熵作为优化准则对通道间相位偏差进行估计，但图像熵对通道间相位偏差的敏感度在非均匀场景中较图像对比度低，本发明在非均匀场景将具备更高的估计精度。李真芳等人，performanceimprovementforconstellationsarusingsignalprocessingtechniques.ieeetransactionsonaerospaceandelectronicsystems,2006,42(2):436-452.该文提出利用子空间估计理论给出了相位偏差校正算法，而本发明提出的利用子孔径sar图像的对比度对通道间相位偏差进行估计，两者原理方法区别明显。专利文献cn105158759b(申请号：cn201510466324.9)公开了一种基于杂波相消的hrwssar通道相位偏差校正方法，主要解决了方位多通道sar通道成像处理面临的相位偏差校正问题，包括以下步骤：(1)输入多通道hrwssar回波数据；(2)各通道回波数据方位向傅里叶变换；(3)选择参考通道；(4)计算各通道用于重构参考通道回波的信号分量；(5)构造杂波相消优化代价函数(6)将杂波相消问题转为恒模优化问题；(7)求解恒模优化问题；(8)输出通道间相位偏差。该发明提出的最大化杂波相消的相位偏差估计算法通过最大化杂波相消性能实现通道间相位误差估计，而本发明利用子孔径图像的对比度进行通道间相位偏差估计，两者原理方法区别明显。技术实现要素：针对现有技术中的缺陷，本发明的目的是提供一种方位多通道sar系统通道间相位偏差校正方法和系统。根据本发明提供的方位多通道sar系统通道间相位偏差校正方法，包括：步骤1，输入方位多通道hrwssar回波数据，通道数为m；步骤2，对各通道回波数据进行距离压缩，实现距离向成像；步骤3，构造参考函数对各通道回波数据进行方位向去斜；步骤4，构造计算子孔径无模糊图像重构函数，用于重构子孔径无模糊sar图像；步骤5，计算子孔径无模糊sar图像的对比度；步骤6，以最大化子孔径无模糊sar图像的对比度为优化目标，求解优化问题，得到优化结果；步骤7，根据优化结果输出通道间相位偏差。优选的，所述步骤3对各通道回波数据进行方位去斜处理，计算公式为：sd,m(fa)＝ft{sm(ta)exp(jπka(ta δta,m)2)}下标d代表去斜操作，ft{·}为傅里叶变换，ka为多普勒调频率；下标m表示第m个通道；下标a表示方位向；sd,m(·)表示第m个通道去斜后的sar图像；fa表示多普勒频率；sm表示第m个通道距离压缩后的信号；ta表示方位时间；j表示虚数单位；δta,m表示第m通道相对参考通道的方位时延。优选的，所述步骤4中构造子孔径无模糊图像重构函数，包括如下步骤：步骤4.1，构造方位多通道hrwssar子孔径图像间传递函数矩阵h(fa)；步骤4.2，对传递函数矩阵h(fa)求逆，得到子孔径无模糊图像重构矩阵p(fa)。优选的，所述步骤5中构造多通道子孔径图像对比度优化代价函数，包括如下步骤：步骤5.1，自变量选为通道间相位偏差φm，构造相位偏差校正矩阵：γ＝diag(1,e-jφ2,…,e-jφm)，j为虚数单位；步骤5.2，计算在步骤5.1中得到的相位偏差下子孔径无模糊图像：s(fa)＝p(fa)γse(fa)，其中，se(fa)为各通道去斜后子孔径图像，fa为多普勒频率；步骤5.3，计算步骤5.2得到的子孔径无模糊图像的图像对比度：优选的，所述步骤6中以最大化子孔径图像对比度为优化目标，采用牛顿迭代法求解优化问题。根据本发明提供的方位多通道sar系统通道间相位偏差校正系统，包括：模块m1，输入方位多通道hrwssar回波数据，通道数为m；模块m2，对各通道回波数据进行距离压缩，实现距离向成像；模块m3，构造参考函数对各通道回波数据进行方位向去斜；模块m4，构造计算子孔径无模糊图像重构函数，用于重构子孔径无模糊sar图像；模块m5，计算子孔径无模糊sar图像的对比度；模块m6，以最大化子孔径无模糊sar图像的对比度为优化目标，求解优化问题，得到优化结果；模块m7，根据优化结果输出通道间相位偏差。优选的，所述模块m3对各通道回波数据进行方位去斜处理，计算公式为：sd,m(fa)＝ft{sm(ta)exp(jπka(ta δta,m)2)}下标d代表去斜操作，ft{·}为傅里叶变换，ka为多普勒调频率；下标m表示第m个通道；下标a表示方位向；sd,m(·)表示第m个通道去斜后的sar图像；fa表示多普勒频率；sm表示第m个通道距离压缩后的信号；ta表示方位时间；j表示虚数单位；δta,m表示第m通道相对参考通道的方位时延。优选的，所述模块m4中构造子孔径无模糊图像重构函数，包括如下模块：模块m4.1，构造方位多通道hrwssar子孔径图像间传递函数矩阵h(fa)；模块m4.2，对传递函数矩阵h(fa)求逆，得到子孔径无模糊图像重构矩阵p(fa)。优选的，所述模块m5中构造多通道子孔径图像对比度优化代价函数，包括如下模块：模块m5.1，自变量选为通道间相位偏差φm，构造相位偏差校正矩阵：γ＝diag(1,e-jφ2,…,e-jφm)，j为虚数单位；模块m5.2，计算在模块m5.1中得到的相位偏差下子孔径无模糊图像：s(fa)＝p(fa)γse(fa)，其中，se(fa)为各通道去斜后子孔径图像，fa为多普勒频率；模块m5.3，计算模块m5.2得到的子孔径无模糊图像的图像对比度：优选的，所述模块m6中以最大化子孔径图像对比度为优化目标，采用牛顿迭代法求解优化问题。与现有技术相比，本发明具有如下的有益效果：1.本发明采用从最大图像对比度出发进行的相位偏差校正方法，在国内外尚数首次，填补了现有技术的空白；2.本发明有效解决了方位多通道hrwssar系统相位偏差校正问题，可有效保障sar成像质量。附图说明通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：图1是星载多通道hrwssar观测几何示意图；图2是本发明的数据处理流程；图3是本发明对实测数据试验结果。具体实施方式下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。实施例：如图2，本实例提供了一种基于子孔径图像最大对比度的方位多通道高分辨率宽幅(hrws)合成孔径雷达(sar)系统通道间相位偏差校正方法，包括以下步骤：步骤1，输入方位多通道hrwssar回波数据，通道数为m；步骤2，对各通道回波数据进行距离压缩，实现sar距离向成像；步骤3，构造参考函数对各通道回波数据进行方位向去斜，得到子孔径方位去斜sar图像se(fa)；步骤4，构造子孔径无模糊图像重构函数矩阵p，实现无模糊sar图像重构；步骤5，计算多通道子孔径图像对比度，作为通道间相位偏差估计所需的优化代价函数；步骤6，以最大化子孔径图像对比度为优化目标，求解优化问题；步骤7，根据步骤7得到的优化结果输出通道间相位偏差。进一步地，所述步骤4中构造子孔径无模糊图像重构函数，包括步骤如下：步骤4.1，构造方位多通道hrwssar子孔径图像间传递函数矩阵h；步骤4.2，传递函数矩阵h求逆，得到子孔径无模糊图像重构矩阵p。进一步地，所述步骤5中构造多通道子孔径图像对比度优化代价函数，包括步骤如下：步骤5.1，自变量选为通道间相位偏差φm，构造相位偏差校正矩阵：γ＝diag(1,e-jφ2,…,e-jφm)，j为虚数单位；步骤5.2，计算在步骤5.1中得到的相位偏差下子孔径无模糊图像：s(fa)＝pγse(fa)，其中，se(fa)为各通道去斜后子孔径图像，fa为多普勒频率；步骤5.3，计算步骤5.2得到的子孔径无模糊图像的图像对比度：下面结合附图对本实施实例进一步描述。参照图1，本实施案例的理论基础简介描述如下：假定方位多通道hrwssar共有m个接收通道，经距离压缩后，对第m个通道信号补偿一常数相位可将其回波信号可视为参考通道s0的方位延迟：sm(ta)＝s0(ta δta,m)其中，sm(ta)为第m个通道接收信号获取的信号，τ为距离向时间，ta为方位时间，δta,m为第m个通道相对参考通道的时间延迟。对方位多通道hrwssar距离压缩后信号进行方位去斜：sd,m(fa)＝ft{sm(ta)exp(jπka(ta δta,m)2)}下标d代表去斜操作，ft{·}为傅里叶变换，ka为多普勒调频率。假定去斜后参考通道子孔径无模糊信号可表示为sd,0(fa)，由于方位多通道hrwssar系统发射脉冲频率小于信号带宽，去斜后sar图像将存在模糊。第m个通道去斜后模糊信号可表示为：式中，上标p为模糊标识，prf为脉冲重复频率。不失一般性，这里假设多普勒中心为零，忽略多普勒频率|fa|＞m·prf/2的信号，将各通道去斜后信号进一步写为矩阵形式：式中，s0(fa)＝[s0(fa k0·prf)s0(fa (k0 1)·prf)…s0(fa (k0 m-1)·prf)]th(fd)＝[h1(fa)h2(fa)…hm(fa)]thm(fd)＝[h1,mh2,m…hm,m]thn,m＝exp(-j2π(fa (k0 n-1)prf)δta,m)其中，为向下取整运算符。由式可知，定义p(fd)＝h-1(fd)，考虑通道间相位偏差校正：式中，γ＝diag(1,e-jφ2,…,e-jφm)，φm为第m个通道相对参考通道的相位偏差。通道间相位偏差将导致子孔径多通道重构图像方位模糊度指标下降，导致子孔径sar图像对比度下降。因此，可利用子孔径图像对比度作为优化目标，通过最大化子孔径图像对比度实现通道间相位偏差估计。多通道sar子孔径sar图像的对比度可表示为：通道间相位偏差可通过求解如下所示的优化问题获得：上述优化问题，可利用牛顿迭代法求解。下面结合实测数据对本发明的效果做进一步的说明。这里采用某机载x波段方位多通道sar试验系统回波数据验证本文通道相位偏差估计算法的性能，主要系统参数如表1所示。该系统第1通道的sar图像如图3a所示。表1x波段方位四通道sar试验系统参数参数名取值载频9.6ghz发射信号带宽20mhz距离向采样率24mhz工作prf100hz多普勒带宽280hz载机速度148m/s通道数4等效相位中心间距0.4m采用本实施例提供的通道相位偏差校正方法估计得到的系统各通道相对于第1通道的相位偏差分别为0、-19.86°、37.80°及-32.32°。利用本发明估计得到的通道相位偏差对原始回波进行校正前后聚焦得到的sar图像如图3b和图3c所示。通过比较可知，经本发明进行通道相位偏差校正后sar图像质量得到了显著的提升，由通道间相位偏差引起的sar图像方位模糊得到了有效的抑制(如图中圆圈标注区域)。上述试验有效验证了本发明对于方位多通道hrwssar系统通道间相位偏差估计校正的有效性。根据本发明提供的方位多通道sar系统通道间相位偏差校正系统，包括：模块m1，输入方位多通道hrwssar回波数据，通道数为m；模块m2，对各通道回波数据进行距离压缩，实现距离向成像；模块m3，构造参考函数对各通道回波数据进行方位向去斜；模块m4，构造计算子孔径无模糊图像重构函数，用于重构子孔径无模糊sar图像；模块m5，计算子孔径无模糊sar图像的对比度；模块m6，以最大化子孔径无模糊sar图像的对比度为优化目标，求解优化问题，得到优化结果；模块m7，根据优化结果输出通道间相位偏差。优选的，所述模块m3对各通道回波数据进行方位去斜处理，计算公式为：sd,m(fa)＝ft{sm(ta)exp(jπka(ta δta,m)2)}下标d代表去斜操作，ft{·}为傅里叶变换，ka为多普勒调频率；下标m表示第m个通道；下标a表示方位向；sd,m(·)表示第m个通道去斜后的sar图像；fa表示多普勒频率；sm表示第m个通道距离压缩后的信号；ta表示方位时间；j表示虚数单位；δta,m表示第m通道相对参考通道的方位时延。优选的，所述模块m4中构造子孔径无模糊图像重构函数，包括如下模块：模块m4.1，构造方位多通道hrwssar子孔径图像间传递函数矩阵h(fa)；模块m4.2，对传递函数矩阵h(fa)求逆，得到子孔径无模糊图像重构矩阵p(fa)。优选的，所述模块m5中构造多通道子孔径图像对比度优化代价函数，包括如下模块：模块m5.1，自变量选为通道间相位偏差φm，构造相位偏差校正矩阵：γ＝diag(1,e-jφ2,…,e-jφm)，j为虚数单位；模块m5.2，计算在模块m5.1中得到的相位偏差下子孔径无模糊图像：s(fa)＝p(fa)γse(fa)，其中，se(fa)为各通道去斜后子孔径图像，fa为多普勒频率；模块m5.3，计算模块m5.2得到的子孔径无模糊图像的图像对比度：优选的，所述模块m6中以最大化子孔径图像对比度为优化目标，采用牛顿迭代法求解优化问题。本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以，本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件，而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构；也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。当前第1页1 2 3
技术特征：1.一种方位多通道sar系统通道间相位偏差校正方法，其特征在于，包括：
步骤1，输入方位多通道hrwssar回波数据，通道数为m；
步骤2，对各通道回波数据进行距离压缩，实现距离向成像；
步骤3，构造参考函数对各通道回波数据进行方位向去斜；
步骤4，构造计算子孔径无模糊图像重构函数，用于重构子孔径无模糊sar图像；
步骤5，计算子孔径无模糊sar图像的对比度；
步骤6，以最大化子孔径无模糊sar图像的对比度为优化目标，求解优化问题，得到优化结果；
步骤7，根据优化结果输出通道间相位偏差。
2.根据权利要求1所述的方位多通道sar系统通道间相位偏差校正方法，其特征在于，所述步骤3对各通道回波数据进行方位去斜处理，计算公式为：
sd,m(fa)＝ft{sm(ta)exp(jπka(ta δta,m)2)}
下标d代表去斜操作，ft{·}为傅里叶变换，ka为多普勒调频率；
下标m表示第m个通道；下标a表示方位向；sd,m(·)表示第m个通道去斜后的sar图像；fa表示多普勒频率；sm表示第m个通道距离压缩后的信号；ta表示方位时间；j表示虚数单位；δta,m表示第m通道相对参考通道的方位时延。
3.根据权利要求1所述的方位多通道sar系统通道间相位偏差校正方法，其特征在于，所述步骤4中构造子孔径无模糊图像重构函数，包括如下步骤：
步骤4.1，构造方位多通道hrwssar子孔径图像间传递函数矩阵h(fa)；
步骤4.2，对传递函数矩阵h(fa)求逆，得到子孔径无模糊图像重构矩阵p(fa)。
4.根据权利要求1所述的方位多通道sar系统通道间相位偏差校正方法，其特征在于，所述步骤5中构造多通道子孔径图像对比度优化代价函数，包括如下步骤：
步骤5.1，自变量选为通道间相位偏差φm，构造相位偏差校正矩阵：
j为虚数单位；
步骤5.2，计算在步骤5.1中得到的相位偏差下子孔径无模糊图像：
s(fa)＝p(fa)γse(fa)，其中，se(fa)为各通道去斜后子孔径图像，fa为多普勒频率；
步骤5.3，计算步骤5.2得到的子孔径无模糊图像的图像对比度：
5.根据权利要求1所述的方位多通道sar系统通道间相位偏差校正方法，其特征在于，所述步骤6中以最大化子孔径图像对比度为优化目标，采用牛顿迭代法求解优化问题。
6.一种方位多通道sar系统通道间相位偏差校正系统，其特征在于，包括：
模块m1，输入方位多通道hrwssar回波数据，通道数为m；
模块m2，对各通道回波数据进行距离压缩，实现距离向成像；
模块m3，构造参考函数对各通道回波数据进行方位向去斜；
模块m4，构造计算子孔径无模糊图像重构函数，用于重构子孔径无模糊sar图像；
模块m5，计算子孔径无模糊sar图像的对比度；
模块m6，以最大化子孔径无模糊sar图像的对比度为优化目标，求解优化问题，得到优化结果；
模块m7，根据优化结果输出通道间相位偏差。
7.根据权利要求6所述的方位多通道sar系统通道间相位偏差校正系统，其特征在于，所述模块m3对各通道回波数据进行方位去斜处理，计算公式为：
sd,m(fa)＝ft{sm(ta)exp(jπka(ta δta,m)2)}
下标d代表去斜操作，ft{·}为傅里叶变换，ka为多普勒调频率；
下标m表示第m个通道；下标a表示方位向；sd,m(·)表示第m个通道去斜后的sar图像；fa表示多普勒频率；sm表示第m个通道距离压缩后的信号；ta表示方位时间；j表示虚数单位；δta,m表示第m通道相对参考通道的方位时延。
8.根据权利要求6所述的方位多通道sar系统通道间相位偏差校正系统，其特征在于，所述模块m4中构造子孔径无模糊图像重构函数，包括如下模块：
模块m4.1，构造方位多通道hrwssar子孔径图像间传递函数矩阵h(fa)；
模块m4.2，对传递函数矩阵h(fa)求逆，得到子孔径无模糊图像重构矩阵p(fa)。
9.根据权利要求6所述的方位多通道sar系统通道间相位偏差校正系统，其特征在于，所述模块m5中构造多通道子孔径图像对比度优化代价函数，包括如下模块：
模块m5.1，自变量选为通道间相位偏差φm，构造相位偏差校正矩阵：
j为虚数单位；
模块m5.2，计算在模块m5.1中得到的相位偏差下子孔径无模糊图像：
s(fa)＝p(fa)γse(fa)，其中，se(fa)为各通道去斜后子孔径图像，fa为多普勒频率；
模块m5.3，计算模块m5.2得到的子孔径无模糊图像的图像对比度：
10.根据权利要求6所述的方位多通道sar系统通道间相位偏差校正系统，其特征在于，所述模块m6中以最大化子孔径图像对比度为优化目标，采用牛顿迭代法求解优化问题。
技术总结本发明提供了一种方位多通道SAR系统通道间相位偏差校正方法和系统，包括：步骤1，输入方位多通道HRWS SAR回波数据，通道数为M；步骤2，对各通道回波数据进行距离压缩，实现距离向成像；步骤3，构造参考函数对各通道回波数据进行方位向去斜；步骤4，构造计算子孔径无模糊图像重构函数，用于重构子孔径无模糊SAR图像；步骤5，计算子孔径无模糊SAR图像的对比度；步骤6，以最大化子孔径无模糊SAR图像的对比度为优化目标，求解优化问题，得到优化结果；步骤7，根据优化结果输出通道间相位偏差。本发明可有效解决方位多通道HRWS SAR成像处理中面临的通道间向为偏差影响问题，有效保障系统成像性能。
技术研发人员：陈筠力;刘艳阳;陈红飞;陈重华;路瑞峰;侯雨生;赵迪;范军
受保护的技术使用者：上海卫星工程研究所
技术研发日：2020.11.11
技术公布日：2021.03.12
