(完整word)各种SAR成像算法总结,推荐文档.docVIP

各种 SAR 成像算法总结 1 SAR 成像原理 SAR 成像处理的目的是要得到目标区域散射系数的二维分布，它是一个二 维相关处理过程， 通常可以分成距离向处理和方位向处理两个部分。 在处理过程中，各算法的区别在于如何定义雷达与目标的距离模型以及如何解决距离－方位耦合问题，这些问题直接导致了各种算法在成像质量和运算量方面的差异。 一般来说，忽略多普勒频移所引起的距离向相位变化， 距离向处理变为一维的移不变过程且相关核已知， 即退化为一般的脉冲压缩处理； 同时将雷达与目标的距离按 2 阶 Taylor 展开并忽略高次项，则方位向处理也是一个一维的移不变过 程，并退化为一般的脉冲压缩处理，这就是经典的距离多普勒（ Range-Doppler RD）算法的实质。 若考虑多普勒频移对距离向相位的影响， 同时精确的建立雷达与目标的距离模型，则不论距离向处理还是方位向处理都变为二维的移变相关过程。 线性调频尺度变换（ Chirp-Scaling CS）算法即在此基础之上将二维数据变换到频域， 利用 Chirp Scaling 原理及频域的相位校正方法，对二维数据进行距离徙动校正处理、 距离向及方位向的聚焦处理，最终完成二维成像处理。 当方位向数据积累延迟小于全孔径时间 （即方位向为子孔径数据） 的情况下， 方位向处理必须使用去斜（ dechirp）处理及频谱分析的方法。在 RD 和 CS 算法 的基础之上，采用 dechirp 处理及频谱分析的方法完成方位向处理的算法分别称 为频谱分析（ SPECAN）算法和扩展 CS（Extended Chirp Scaling ECS）算法。 1.1 SAR 成像原理 本节以基本的正侧视条带工作模式为例，对 SAR 的成像原理进行分析和讨 论。 正侧视条带 SAR 的空间几何关系如下图所示。 图中，αoβ平面为地平面， oγ 垂直于 αoβ平面。 SAR 运动平台位于 S 点，其在地面的投影为 G 点。SAR 运动 平台的运动方向 Sx 平行于 oβ，速度大小为 va 。SAR 天线波束中心与地面的交点 为 C，CG 与运动方向 Sx 垂直； S 与 C 的距离为 Rs ， B1SB2 称为天线波束的方 位向宽度，大小为 a 。 P 为测绘带内的某一点，一般情况下取斜距平面 CSP 进 行分析，称 SAR 运动的方向 Sx 为方位向（或方位维），称天线波束指向方向 SC 为距离向（或距离维）。 x γ va S βa R(t) β G 测绘带 Rs o P B1 vat C B2 X α 正侧视条带 SAR 几何关系示意图 假定 P 的方位向坐标为 X ；在 t 时刻， SAR 运动平台 S 与 P 的距离为 R t 。 若当 t 0 时刻， SAR 运动平台位于方位向 0 点，则当 t 时刻， R(t) 的表达式为： R(t ) Rs2 (vat X ) 2 (1.1) 将式 (1.1)在 t X /Va 附近进行 2 阶 Taylor 展开，有： R X 1 R 2 R(t ) R X t X X t X Va Va Va 2 Va Va (vat X )2 (1.2) Rs 2Rs 假设雷达发射连续的正弦波，即发射信号 st (t ) 为： st (t ) Re[ Ae j ct ] (1.3) 其中， A 为发射正弦波的幅度， c 为发射信号的载频。 发射信号 st (t) 经点目标 P 散射后，雷达接收机收到的信号 sr (t ) 为： sr (t ) Re KAe j c (t ) F ( x) (1.4) 其中： c 为光速， K 为复常数， 为回波信号相对于发射信号的时间延迟： 2R(t) / c (1.5) F ( x) 为考虑雷达水平方向增益变化而引入的加权函数。若不考虑雷达天线 的加权作用，即令 F (x) 1，则式 (1.4)变为： s (t ) Re KAej c ( t ) (1.6) r 根据式 (1.6)，雷达运动平台相对于点目标的运动将造成回波信号的相位随时间不断变化，从而引起回波瞬时频率的变化，产生多普勒频移。多普勒频移量 fd (t ) 为： 1 d( c ) 1 d 2 c R(t ) f d (t ) dt 2 dt (1.7) 2 c 将式 (1.2)内的 R(t ) 代入可得： f d (t ) 1 d 2 c R(t ) 2 dt c 2va2 (1.8) 1 d 2 c Rs 1 (vat X )2 / 2Rs2 (t t0 ) 2 dt c Rs 其中： 为雷达工作波长，且 2 c / c ，t 0 X / va 为雷达波束中心通过 P 点的时间。 回波信号的瞬时频率 fr (