基于BP算法的高分辨率P波段SAR空变运动误差补偿算法制造技术

本发明专利技术公开了一种基于BP算法的高分辨率P波段SAR空变运动误差补偿算法，本发明专利技术首先利用常规方法进行方位空不变误差补偿，然后按照BP算法进行剩余的空变运动误差补偿。剩余的空变运动误差补偿时利用二次相位扰动技术降低了运算量，节省了时间；并且由于BP算法自身的优点，在二维时域对单个目标单独进行补偿，避免了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图像的二维频域数据；步骤二，选择一个时间带宽积值Α，在二维频域上构造一个调频率随距离频率变化的方位线性调频信号:步骤三，利用步骤二得到的方位线性调频信号对步骤一得到的二维频域数据进行相位扰动，并返回二维时域，得到二维时域数据:步骤四，计算二维时域中的合成孔径覆盖点数L，进而得到以每个点为中心的能量范围；步骤五，构造SAR图像像素网格I (mTm，nTn)，其中m，η分别为距离像素位置和方位向像素位置，Tm, Tn分别为距离、方位向像素空间间隔，τ..:=γ , m= Ι,.,.,Μ, η =I,..., N, Μ, N分别为SAR图像距离、方位向像素个数，C为光速，Fs为脉冲采样率，Prf为脉冲重复间隔，V为载机参考速度；步骤六，对于步骤五构造的SAR图像像素网格I (mTm，nTn)中每一个像素点，根据步骤四确定的能量范围，计算每个像素点的回波在二维时域中的位置曲线；步骤七，针对步骤六得到的位置曲线，利用合成孔径雷达空间几何关系，计算得到所述位置曲线上的方位空变运动误差；步骤八，针对步骤六得到的位置曲线，进而得知所述位置曲线上的能量，根据步骤七得到的方位空变运动误差利用BP方位空变误差补偿算法先进行相位校正然后再将能量叠加，最后得到SAR图像像素点的复数灰度值。本专利技术的有益效果:1.与现有空变运动误差补偿算法相比，本专利技术首先利用常规方法进行方位空不变误差补偿，然后按照BP算法进行剩余的空变运动误差补偿。剩余的空变运动误差补偿时利用二次相位扰动技术降低了运算量，节省了时间。2.通过步骤三的相位扰动，并返回二维时域，在二维时域中进行合成孔径，不随距离向变化，因此不同距离单元目标具有相同的合成孔径长度，为后续的误差补偿及积累带来便利。3.并且由于BP算法自身的优点，在二维时域对单个目标单独进行补偿，避免了SAR数据的二维耦合性，可以进行距离迁移误差的校正，因此能够适应更大的载机运动误差，具有较高的补偿精度。附图说明图1为本专利技术的流程图；图2为ω-Κ算法+常规运动补偿后的距离迁移校正效果；图3为原始合成孔径长度与线性调频信号扰动后合成孔径长度对比；图4为对成像结果点目标进行指标测量；具体实施例方式为了更好地理解本专利技术的技术方案，下面结合附图及具体实施例对本专利技术做进一步详细描述。本专利技术提出了一种基于BP算法的高分辨率P波段SAR空变运动误差补偿算法——相位扰动BP空变误差补偿算法，与已有方位空变误差补偿算法相比，由于新方法能够补偿方位空变误差导致的距离迁移校正误差，因此能够适应更大的方位空变误差。如图1，该方法具体为:步骤一，对接收到的SAR原始回 波数据利用结合运动补偿的ω-Κ算法进行处理，得到SAR图像的二维频域数据；具体为:将SAR原始回波数据进行距离压缩、距离迁移校正、相位补偿、Stolt插值之后得到二维频域数据S(f τ，f n)，其中f\，f n分别为距离频率坐标和方位频率坐标；步骤一先利用ω-Κ算法以及常规运动补偿算法进行基本的成像处理以及空不变、距离空变运动误差补偿，能够补偿掉大部分的运动误差，使得后续补偿步骤的推导成立。步骤二，选择一个时间带宽积值A (—般选择100)，在步骤一的得到的二维频域数据范围内，构造一个调频率随距离频率变化的方位线性调频信号p(f\，fn):本文档来自技高网...

【技术保护点】
一种基于BP算法的高分辨率P波段SAR空变运动误差补偿算法，其特征在于，步骤一，对接收到的SAR原始回波数据利用结合运动补偿的ω?κ算法进行处理，得到SAR图像的二维频域数据；步骤二，选择一个时间带宽积值A，在二维频域上构造一个调频率随距离频率变化的方位线性调频信号：步骤三，利用步骤二得到的方位线性调频信号对步骤一得到的二维频域数据进行相位扰动，并返回二维时域，得到二维时域数据：步骤四，计算二维时域中的合成孔径覆盖点数L，进而得到以每个点为中心的能量范围；步骤五，构造SAR图像像素网格I(mTm,nTn)，其中m,n分别为距离像素位置和方位向像素位置，Tm,Tn分别为距离、方位向像素空间间隔，m＝1,...,M，n＝1,...,N，M,N分别为SAR图像距离、方位向像素个数，C为光速，Fs为脉冲采样率，Prf为脉冲重复间隔，v为载机参考速度；步骤六，对于步骤五构造的SAR图像像素网格I(mTm,nTn)中每一个像素点，根据步骤四确定的能量范围，计算每个像素点的回波在二维时域中的位置曲线；步骤七，针对步骤六得到的位置曲线，利用合成孔径雷达空间几何关系，计算得到所述位置曲线上的方位空变运动误差；步骤八，针对步骤六得到的位置曲线，进而得知所述位置曲线上的能量，根据步骤七得到的方位空变运动误差利用BP方位空变误差补偿算法先进行相位校正然后再将能量叠加，最后得到SAR图像像素点的复数灰度值。FDA00003061911800011.jpg,FDA00003061911800012.jpg...

【技术特征摘要】
1.一种基于BP算法的高分辨率P波段SAR空变运动误差补偿算法，其特征在于， 步骤一，对接收到的SAR原始回波数据利用结合运动补偿的ω-Κ算法进行处理，得到SAR图像的二维频域数据； 步骤二，选择一个时间带宽积值Α，在二维频域上构造一个调频率随距离频率变化的方位线性调频信号: 步骤三，利用步骤二得到的方位线性调频信号对步骤一得到的二维频域数据进行相位扰动，并返回二维时域，得到二维时域数据: 步骤四，计算二维时域中的合成孔径覆盖点数L，进而得到以每个点为中心的能量范围； 步骤五，构造SAR图像像素网格I (mTm，nTn)，其中m，η分别为距离像素位置和方位向像素位置，Tffl, Tn 分别为距离、方位向像素空间间隔，I ，& =I，m = 1，...，Μ，η =I,..., N, ...

【专利技术属性】
技术研发人员：孟大地，林雪，胡玉新，丁赤飚，
申请(专利权)人：中国科学院电子学研究所，
类型：发明
国别省市：