[发明专利]基于BP算法的高分辨率P波段SAR空变运动误差补偿算法

说明书

技术领域

本发明涉及一种基于BP算法的高分辨率P波段SAR空变运动误差补偿算法，涉及机载SAR成像处理及运动补偿领域。

背景技术

机载合成孔径雷达（SAR）需要利用成像处理算法对SAR原始回波数据进行匹配滤波处理，得到SAR图像。由于时域处理算法非常耗时，一般利用频域处理算法（如WK算法、Chirp Scaling算法、距离多普勒算法）进行成像处理。但频域处理算法只能适应载机做理想水平匀速直线运动的情况，需要通过运动补偿算法补偿载机的非理想运动。由于SAR信号的二维耦合性，出于运算效率考虑，运动补偿算法一般只能补偿运动误差的空不变分量以及距离空变分量。在SAR系统分辨率较低、合成孔径时间较短（如X波段）、载机运动误差较小的情况下，上述处理方式能够满足各种应用要求。但对于高分辨率P波段SAR系统，由于其距离单元较小、合成孔径时间较长，在载机运动误差较大时，常规运动补偿后的剩余方位空变运动误差不但导致方位向相位误差，同时导致距离迁移校正误差，因此上述方位空变相位误差补偿不能对目标回波信号进行完全积累，从而导致SAR图像分辨率降低，成像质量下降。因此需要提供一种空变运动误差补偿方法，能够解决高分辨率P波段SAR系统在运动误差较大时的方位空变运动误差补偿，尤其是所导致的距离迁移校正误差，从而使得高分辨率P波段SAR系统即使在较大的运动误差时也能得到良好聚焦的SAR图像产品。

反投影（Back Projection,BP）算法是一种用于SAR成像处理的时域逐点处理算法。与频域处理算法对载机运动误差较差的适应性相比，BP算法由于其时域逐点处理特性，能够在任何的载机运动状态下根据载机运动轨迹精确恢复出SAR图像每个像素点的灰度值。与频域处理算法相比，BP算法虽然能够适应大运动误差，但却需要数百倍的运算量，在实际应用中是无法忍受的，大大限制了BP算法的应用。

发明内容

本发明提出了一种基于BP算法的高分辨率P波段SAR空变运动误差补偿算法，通过相位扰动BP空变误差补偿，能够补偿方位空变误差导致的距离迁移校正误差，因此能够适应更大的方位空变误差。

本发明的目的是通过下述技术方案实现的：

步骤一，对接收到的SAR原始回波数据利用结合运动补偿的ω-κ算法进行处理，得到SAR图像的二维频域数据；

步骤二，选择一个时间带宽积值A，在二维频域上构造一个调频率随距离频率变化的方位线性调频信号：

步骤三，利用步骤二得到的方位线性调频信号对步骤一得到的二维频域数据进行相位扰动，并返回二维时域，得到二维时域数据：

步骤四，计算二维时域中的合成孔径覆盖点数L，进而得到以每个点为中心的能量范围；

步骤五，构造SAR图像像素网格I(mT_m,nT_n)，其中m,n分别为距离像素位置和方位向像素位置，T_m,T_n分别为距离、方位向像素空间间隔，m＝1,...,M，n＝1,...,N，M,N分别为SAR图像距离、方位向像素个数，C为光速，Fs为脉冲采样率，P_rf为脉冲重复间隔，v为载机参考速度；

步骤六，对于步骤五构造的SAR图像像素网格I(mT_m,nT_n)中每一个像素点，根据步骤四确定的能量范围，计算每个像素点的回波在二维时域中的位置曲线；

步骤七，针对步骤六得到的位置曲线，利用合成孔径雷达空间几何关系，计算得到所述位置曲线上的方位空变运动误差；

步骤八，针对步骤六得到的位置曲线，进而得知所述位置曲线上的能量，根据步骤七得到的方位空变运动误差利用BP方位空变误差补偿算法先进行相位校正然后再将能量叠加，最后得到SAR图像像素点的复数灰度值。

本发明的有益效果：

1.与现有空变运动误差补偿算法相比，本发明首先利用常规方法进行方位空不变误差补偿，然后按照BP算法进行剩余的空变运动误差补偿。剩余的空变运动误差补偿时利用二次相位扰动技术降低了运算量，节省了时间。

2.通过步骤三的相位扰动，并返回二维时域，在二维时域中进行合成孔径，不随距离向变化，因此不同距离单元目标具有相同的合成孔径长度，为后续的误差补偿及积累带来便利。

3.并且由于BP算法自身的优点，在二维时域对单个目标单独进行补偿，避免了SAR数据的二维耦合性，可以进行距离迁移误差的校正，因此能够适应更大的载机运动误差，具有较高的补偿精度。

附图说明

图1为本发明的流程图；