[发明专利]一种补偿调频连续波圆周SAR脉内运动的成像方法

摘要

本发明属于合成孔径雷达成像技术领域，公开了一种补偿调频连续波圆周SAR脉内运动的成像方法，本发明通过在成像过程中增加一个仅和调频连续波圆周SAR成像系统参数有关的判决过程来快速判定此系统在时域方法成像过程中是否可以忽略脉内运动，并指导选取合适的成像方法，若可以忽略脉内运动，采用脉冲体制圆周SAR时域方法来成像，反之，对于不能忽略脉内运动的情况，采用本发明的通过在距离域补偿脉内运动引起的距离偏移误差的时域方法来成像。本发明解决了调频连续波圆周SAR成像系统在时域方法成像过程中停‑走‑停近似不再有效以及是否可以忽略脉内运动的问题，实现了调频连续波圆周SAR时域方法高效高分辨成像。

主权项

1.一种补偿调频连续波圆周SAR脉内运动的成像方法，其特征在于，包括下列步骤：步骤1获取调频连续波圆周SAR回波信号，包括调频连续波圆周SAR成像系统参数初始化及调频连续波圆周SAR回波信号建模，基于以下步骤：(1.1)调频连续波圆周SAR成像系统参数初始化：雷达平台圆形运动轨迹半径为R，高度为H，方位向观测角为θ，雷达天线相位中心位置A为(x_r,y_r,H)，其中x_r、y_r和H分别为位置A的x轴、y轴和z轴坐标，A的坐标也可以表示为(Rcosθ,Rsinθ,H)；地平面的圆形观测场景半径为R_a，场景中任意一个点目标P的坐标为(x,y,0)，也可表示为其中r和分别为目标P的径向半径和方位角；雷达运动角速度为ω，载机速度V＝Rω，t时刻对应的方位观测角为θ＝ωt，全时间其中t_m表示慢时间，为快时间；(1.2)调频连续波圆周SAR回波信号建模：雷达位置A与目标P之间的瞬时斜距R_r为：计算R_r在快时间处的一阶泰勒展开为：其中，R_cen为停‑走‑停近似下的瞬时斜距，即脉冲体制圆周SAR的瞬时斜距，为调频连续波圆周SAR脉内运动引入的距离徙动误差，也是脉冲圆周SAR忽略的误差；雷达发射信号为线性调频信号时，目标P的回波信号为：其中，rect(·)为矩形窗函数，c为光速，T_p为脉冲持续时间,且f_c为载频，K_r调频率，且K_r＝B/T_p，B为雷达发射信号带宽；步骤2解线频调处理以及残余视频相位消除和基频回波信号获取，基于以下步骤：(2.1)解线频调处理：回波信号与参考信号的共轭相乘得到差频信号为：其中，参考信号为：R_ref为参考距离，取为场景中心处目标与雷达之间的瞬时斜距，R_Δ＝R_r‑R_ref；差频信号中的最后一个指数项为残余视频相位项；(2.2)残余视频相位消除和基频回波信号获取：对差频信号进行距离向傅里叶变换并与指数项exp(‑jπf_i²/K_r)相乘消除残余视频相位，其中f_i＝‑2K_rR_Δ/c为差频频率，再对差频频率做逆傅里叶变换得到基频回波信号为：其中，λ＝f_c/c为载波波长，R_Δm＝R_cen‑R_ref；步骤3距离压缩并分析脉内平台运动引起的距离偏移误差，基于以下步骤：(3.1)距离压缩：基频回波信号最后一个指数项中的系数为其最大值为4πB/2c，在毫米波段f_c＞＞B/2，则4π/λ＝4πf_c/c＞＞4πB/2c，即基频回波信号中指数项对应的相位远小于对应的相位，此时，忽略基频信号中的最后一个指数项对成像的影响，然后沿距离向进行傅里叶变换得到距离压缩信号为：(3.2)分析脉内运动引起的距离偏移误差：调频连续波圆周SAR脉内平台运动在距离域引起的频率偏移误差为Δf_i＝‑2k₁/λ，将其转换为距离偏移误差为：其中，θ_m为慢时间t_m对应的方位观测角，即θ_m＝ωt_m；计算距离偏移误差ΔR₁的最大值，即调频连续波圆周SAR脉内运动对场景中某一目标引起的最大距离偏移误差为：上式中令r＝R_a，即得到调频连续波圆周SAR脉内运动对整个观测场景引起的最大距离偏移误差，并将其与调频连续波圆周SAR场景中心处的理论平面分辨率Δρ′₀相除得到脉内运动引起的最大距离偏移误差相对于分辨率的比值为ΔC，其中Δρ′₀＝2.4λ/(4πcosθ_d0)，为调频连续波圆周SAR对场景中心目标的俯仰角，ΔC为：步骤4判定调频连续波圆周SAR时域方法成像过程中是否可以忽略脉内运动并指导选取合适的成像方法：将调频连续波圆周SAR成像系统参数带入ΔC的表达式中并计算ΔC的值，当调频连续波圆周SAR脉内运动引起的距离偏移误差小于一个距离分辨单元，即ΔC＜1时，采用调频连续波圆周SAR时域方法成像过程中脉内运动引起的距离偏移误差可忽略不计，此时可以忽略脉内运动，并采用脉冲体制圆周SAR时域方法来成像，即采用停‑走‑停近似的时域成像方法来得到最终的调频连续波圆周SAR图像；反之，当调频连续波圆周SAR脉内运动引起的距离偏移误差大于或等于一个距离分辨单元，即ΔC≥1时，采用调频连续波圆周SAR时域方法成像过程中不能忽略脉内运动，此时需要对脉内运动进行额外的补偿，并采用本发明的补偿调频连续波圆周SAR脉内运动的成像方法来成像，继续下面的步骤；步骤5成像场景网格划分，将雷达观测场景划分为大小一致的若干个单位网格，单位网格的大小要求小于场景中心处的理论平面分辨率Δρ′₀，即Δx′＜Δρ′₀，Δy′＜Δρ′₀，其中，Δx′和Δy′分别为场景网格的距离向间隔和方位向间隔；步骤6计算网格点在某一方位慢时刻包含了脉内平台运动引入的距离偏移误差的新的瞬时斜距R′，具体步骤为：计算网格点在该方位时刻停‑走‑停近似下的瞬时斜距R_cen为：计算脉内平台运动在距离域引入的距离偏移误差ΔR₁为：其中，Δf_i＝‑2k₁/λ为调频连续波圆周SAR脉内平台运动在距离域引起的频率偏移误差，θ_m为慢时间t_m对应的方位观测角，即θ_m＝ωt_m；计算新的瞬时斜距R′＝R_cen+ΔR₁；步骤7根据新的瞬时斜距计算出该方位时刻的距离压缩信号在成像网格中相对应的位置，即将该方位时刻的距离压缩信号按照新的瞬时斜距投影到成像网格中，并记为S′_if(f_i,t_m)；步骤8方位向相位补偿，得到某一方位时刻的子图像为：其中，x′和y′为场景中的网格坐标，exp(j4πR_Δm/λ)为方位向相位补偿函数；步骤9得到最终的调频连续波圆周SAR图像，判断是否得到了所有方位时刻的子图像，即是否得到了方位向所有采样点处的子图像，若没有，进入下一方位时刻并重复上述步骤6～8，若得到了所有方位时刻的子图像，进行子图像相干叠加得到最终的调频连续波圆周SAR图像为：其中，M为整个360°孔径内方位向的采样点数。