CN104122552A - 一种滑动式双站圆周合成孔径雷达成像方法 - Google Patents

Abstract

本发明提供了一种滑动式双站圆周合成孔径雷达成像方法，包括：步骤1，计算滑动式双站圆周合成孔径雷达成像时的雷达接收天线俯仰向和航迹向-3dB波束宽度、雷达发射天线发射信号的脉冲重复频率和场景回波信号接收天线所在平台的圆周运动中心滑动间距；步骤2：根据场景回波信号接收天线所在平台自身惯性导航测量系统输出的姿态参数或雷达惯性导航测量系统输出场景回波信号接收天线和直达波信号接收天线的姿态参数和平台位置信息，并控制场景回波信号接收天线对观测场景的回波信号进行接收；步骤S3：根据步骤S2获取的滑动式双站圆周合成孔径雷达成像回波数据和直达波数据，通过相干叠加处理获取观测场景的二维或三维图像。

Description

一种滑动式双站圆周合成孔径雷达成像方法

技术领域

本发明涉及微波成像技术对地观测领域，是一种滑动式双站圆周合成孔径雷达成像方法。

背景技术

20世纪90年代中后期，Falconer等最早提出了圆周合成孔径雷达(Circular Synthetic Aperture Radar,简称CSAR)三维成像方法(参见Falconer D G,Moussally G J.Tomographic imaging of radar sata gathered on acircular flight path about a three-dimensional target zone.Proceedings of SPIE1995,Apr.,1995.)，主要通过雷达传感器随平台进行360°圆周运动、波束始终指向同一场景区域进行观测区域，实现对观测区域的三维成像数据获取；随后，Soumekh等为实现更为精细的三维成像能力，提出了E-CSAR(Elevation Circular SAR,简称E-CSAR)成像方法(参见Soumekh M.Synthetic aperture radar signal processing with matlab algorithms.New York；Wiley.1999.)，但这两种成像方法主要以实现圆周飞行轨迹中心区域一定范围观测场景的全方位成像能力，其成像区域完全取决于雷达天线俯仰向和航迹向-3dB波束宽度，对单次飞行而言，不能实现大条带式高分辨率全方位成像观测；同时，也不能实现对观测目标非后向散射特性的获取。

国内方面，唐智等针对CSAR成像进行了不同圆周SAR合成孔径方式进行分析，给出了不同平台飞行方式对成像性能的影响，但其隐含的成像区域仍然取决于雷达天线波束宽度(参见唐智,李景文,周荫清,等.曲线合成孔径雷达信号模型与孔径形状研究.系统工程与电子技术,2006,28(8):1115-1119.)；谭维贤等结合双站成像模式的特点给出了一种双站CSAR(Bi-static Circular Synthetic Aperture Radar,简称BiCSAR)成像方法，由于考虑以非后向散射特性获取为主，因此，并没有解决CSAR成像时观测区域较小的问题(参见谭维贤,合成孔径雷达三维成像理论与方法研究.中国科学院电子学研究所博士论文,2009.)；林赟等开展了多种CSAR成像方式的研究，给出了相应成像处理方法，同样，也尚未解决CSAR成像区域小的问题(参见林赟,圆迹合成孔径雷达成像算法研究.中国科学院电子学研究所博士论文,2011.)，也即尚未成像区域受雷达天线俯仰向和航迹向-3dB波束宽度限制的问题。

总的来讲，相对于常规SAR成像模式而言，尽管CSAR成像具有全方位、高分辨和三维观测能力，且通过机载试验获取了高质量的微波图像(参见Octavio Ponce,Pau Prats,Marc Rodriguez-Cassola,et al..Processingof circular SAR trajectories with fast factorized back-projection.IEEEInternational Geoscience&Remote Sensing Symposium,Vancouver,Canada,Jul.24-29,2011,3692-3695.以及参见文献：洪文,圆迹SAR成像技术研究进展.雷达学报,2013,1(2):124-135.)，但由于成像区域的限制，其应用方式一直倍受国内外科研和工业部门的广泛关注。滑动式双站圆周合成孔径雷达成像能够提高雷达系统的噪声等效后向散射系数，可以通过高空定点发射、低空无源被动接收来实现，对军事目标侦察和立体测绘等具有重要实用价值。目前，国内外均尚未给出滑动式双站圆周合成孔径雷达成像方法，因此，如何在获取观测区域非后向散射特性图像的同时，进一步扩大成像区域范围是CSAR发展的重要方向之一。

现有圆周合成孔径雷达(Circular Synthetic Aperture Radar,简称CSAR)获取观测区域范围较小，且无法实现非后向散射信息获取；

现有双站圆周合成孔径雷达(Bi-static Circular Synthetic Aperture Radar,简称BiCSAR)尽管可以获取目标的非后向散射信息，但其成像区域主要仍旧取决于雷达天线俯仰-方位3dB波束宽度，使其观测区域长期驻留在飞行轨迹范围内的中心区域，不利于获取观测区域的大面积非后向高分辨率图像。

发明内容

为了克服现有技术存在的缺陷，本发明提出了一种滑动式双站圆周合成孔径雷达成像方法。

根据本发明提出的滑动式双站圆周合成孔径雷达成像方法，包括步骤：步骤1，计算滑动式双站圆周合成孔径雷达成像时的雷达接收天线俯仰向和航迹向-3dB波束宽度、雷达发射天线发射信号的脉冲重复频率和场景回波信号接收天线所在平台的圆周运动中心滑动间距；步骤2：根据场景回波信号接收天线所在平台自身惯性导航测量系统输出的姿态参数或雷达惯性导航测量系统输出场景回波信号接收天线和直达波信号接收天线的姿态参数和平台位置信息，并控制场景回波信号接收天线对观测场景的回波信号进行接收，雷达发射天线发射具有一定带宽的信号；步骤S3：根据步骤S2获取的滑动式双站圆周合成孔径雷达成像回波数据和直达波数据，通过相干叠加处理获取观测场景的二维或三维图像。

利用本发明提供的滑动式双站圆周合成孔径雷达成像方法，能够实现大范围观测区域成像，不仅能够实现雷达随平台的运动对目标进行全方位观测，而且还能实现保留常规CSAR成像的优势；解决了现有CSAR成像区域受雷达天线俯仰向和航迹向-3dB波束宽度限制的问题，实现了大范围高分辨率全方位圆周成像数据获取与成像处理；同时还能实现非后向散射特性获取。

附图说明：

图1常规圆周合成孔径雷达成像几何示意图；

图2滑动式双站圆周合成孔径雷达成像几何示意图；

图3滑动式双站圆周合成孔径雷达成像刈宽与入射角之间的关系示意图；

图4滑动式双站圆周合成孔径雷达成像飞行轨迹示意图；

图5a是滑动式双站圆周合成孔径雷达成像时场景回波信号接收天线所在平台飞行轨迹与照射区域示意图，图5b是场景回波信号接收天线所在平台单个完整圆周飞行时的飞行轨迹，图5c是场景回波信号接收天线所在平台单个完整圆周飞行时的成像区域；

图6滑动式双站圆周合成孔径雷达成像方法流程图。

具体实施方式：

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。

如图1所示，本发明给出了一种滑动式双站圆周合成孔径雷达成像方法，同常规圆周合成孔径雷达成像(如图1所示)相比，滑动式双站圆周合成孔径雷达成像(Sliding Bi-static CSAR,简称S-BiCSAR)通过高空相对地面静止或缓慢移动的平台携带雷达发射机分系统发射信号，机载平台携带接收机分系统“滑动式”前移，如图2所示，既能实现非后向散射特性获取，还能实现对观测区域的全方位观测，同时还能实现大面积高分辨率图像获取。

滑动式双站圆周合成孔径雷达由1副发射天线、2副接收天线(分别为场景回波信号接收天线和直达波信号接收天线)、单通道发射机、双通道接收机、定时与同步模块和位置测量模块组成。直达波信号接收天线用于接收来自发射天线的雷达信号，场景回波信号接收天线用于接收来自地面场景反射的雷达信号。

为叙述方便，结合滑动式双站圆周合成孔径雷达成像几何对其进行说明。如图2所示，设平台运动时的场景回波信号接收天线的三维位置坐标为(x,y,z)，直达波信号接收天线的三维位置坐标为(x_ref,y_ref,z_ref)，雷达发射天线的三维位置坐标为(x_T,y_T,z_T)，相对地面静止或缓慢移动，且z_T>z，ρ为平台携带包含2副接收天线的接收机分系统的平台进行圆周飞行半径，场景回波信号接收天线俯仰向和航迹向-3dB波束宽度分别为φ_-3dB和θ_-3dB，雷达发射天线俯仰向和方位向-3dB波束宽度分别为φ_-3dBT和θ_-3dBT，目标P_n对应的三维位置坐标为P_n(x_n,y_n,z_n)。Tr表示雷达发射天线，相应三维位置坐标为(x_T,y_T,z_T)，相对地面静止或缓慢移动，雷达发射天线所在平台可以为卫星，或者飞艇；Re表示接收天线，包括场景回波信号接收天线Re(x,y,z)和直达波信号接收天线Re(x_ref,y_ref,z_ref)，场景回波信号接收天线通过携带接收天线和直达波信号接收天线进行圆周飞行，每个圆周飞行对应不同的观测场景区域，ρ_n为观测场景区域的半径，通过场景回波信号接收天线所在平台的多个完整和不完整圆周飞行实现多个观测区域的连续覆盖，从而扩大成像观测范围，这多个完整和不完整的形成过程即为“滑动式”运动。

考虑到直达波接收天线主要直接接收来自发射天线的雷达信号，对直达波接收天线的波束宽度没有严格限制。

本发明的滑动式双站圆周合成孔径雷达成像方法，包括：

步骤S1：滑动式双站圆周合成孔径雷达成像雷达参数计算，根据系统三维分辨率、发射天线所在平台高度、场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线所在平台飞行速度、场景回波信号接收天线入射角和成像刈宽计算滑动式双站圆周合成孔径雷达成像时的雷达接收天线俯仰向和航迹向-3dB波束宽度、雷达发射天线发射信号的脉冲重复频率和场景回波信号接收天线所在平台的圆周运动中心滑动间距，具体地：

步骤S11：根据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线入射角和成像刈宽计算俯仰向-3dB波束宽度，如下所示：

${\mathrm{\φ}}_{-3\mathrm{dB}}=\arctan \left(\frac{2H_{\mathrm{Re}}\×\tan {\mathrm{\φ}}_{\mathrm{in}}+W_{\mathrm{\φ}}}{2H_{\mathrm{Re}}}\right)-\arctan \left(\frac{2H_{\mathrm{Re}}\×\tan {\mathrm{\φ}}_{\mathrm{in}}-W_{\mathrm{\φ}}}{2H_{\mathrm{Re}}}\right)---\left(1\right)$

其中，φ_-3dB为场景回波信号接收天线入射角俯仰向-3dB波束宽度，H_Re为场景回波信号接收天线所在平台飞行高度，φ_in为场景回波信号接收天线入射角，W_φ为成像刈宽，与场景回波信号接收天线入射角以及接收天线所在平台飞行高度的关系如下：

$W_{\mathrm{\φ}}=H_{\mathrm{Re}}\×[\tan ({\mathrm{\φ}}_{\mathrm{in}}+\frac{{\mathrm{\φ}}_{-3\mathrm{dB}}}{2})-\tan ({\mathrm{\φ}}_{\mathrm{in}}-\frac{{\mathrm{\φ}}_{-3\mathrm{dB}}}{2})],$

其中，H_Re∈[1000,1000000]米，(φ_in+φ_-3dB/2)和(φ_in-φ_-3dB/2)分别表示场景回波信号接收天线的最大观测角度和最小观测角度，且0<(φ_in-φ_-3dB/2)<(φ_in+φ_-3dB/2)<π/2，如图3所示；

步骤S11：根据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线入射角和成像刈宽计算航迹向-3dB波束宽度，如下所示：

${\mathrm{\&theta;}}_{-3\mathrm{dB}}\&GreaterEqual;{\mathrm{\&phi;}}_{-3\mathrm{dB}}=\arctan \left(\frac{2H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}+W_{\mathrm{\&phi;}}}{2H_{\mathrm{Re}}}\right)-\arctan \left(\frac{2H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}-W_{\mathrm{\&phi;}}}{2H_{\mathrm{Re}}}\right)---\left(2\right)$

其中，θ_-3dB为场景回波信号接收天线航迹向-3dB波束宽度，φ_-3dB为场景回波信号接收天线俯仰向-3dB波束宽度，H_Re为场景回波信号接收天线所在平台飞行高度，φ_in为场景回波信号接收天线入射角，W_φ为成像刈宽；W_φ为成像刈宽，且

$W_{\mathrm{\&phi;}}=H_{\mathrm{Re}}\&times;[\tan ({\mathrm{\&phi;}}_{\mathrm{in}}+\frac{{\mathrm{\&phi;}}_{-3\mathrm{dB}}}{2})-\tan ({\mathrm{\&phi;}}_{\mathrm{in}}-\frac{{\mathrm{\&phi;}}_{-3\mathrm{dB}}}{2})],$

其中，H_Re∈[1000,1000000]米，(φ_in+φ_-3dB/2)和(φ_in-φ_-3dB/2)分别表示场景回波信号接收天线的最大观测角度和最小观测角度，且0<(φ_in-φ_-3dB/2)<(φ_in+φ_-3dB/2)<π/2；

θ_-3dB≥φ_-3dB表示系统工作时，保证雷达场景回波信号接收天线航迹向-3dB波束宽度大于等于雷达接收天线俯仰向-3dB波束宽度，从而保证滑动式双站圆周合成孔径雷达成像时成像刈宽最大限度地由场景回波信号接收天线俯仰向-3dB波束宽度决定，不降低雷达系统的成像刈宽；若θ_-3dB<φ_-3dB，则滑动式双站圆周合成孔径雷达成像时成像刈宽由场景回波信号接收天线航迹向-3dB波束宽度决定，因此，在进行雷达系统设计时，当成像刈宽确定后，则优先保证θ_-3dB≥φ_-3dB。

步骤S13：根据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线所在平台飞行速度、场景回波信号接收天线入射角和成像刈宽计算雷达发射天线的脉冲重复频率，如下所示：

$\left\{\begin{array}{c}\mathrm{PRF}\&GreaterEqual;\frac{V_{\mathrm{\&theta;}}(\sqrt{A+{\mathrm{\&rho;W}}_{\mathrm{\&phi;}}}-\sqrt{A-{\mathrm{\&rho;W}}_{\mathrm{\&phi;}}})}{\mathrm{\&lambda;\&rho;}}\\ \mathrm{PRF}\&le;\frac{1}{2[\frac{\sqrt{W_{0.5}^2+H_{\mathrm{Tr}}^2}+\sqrt{{(\mathrm{\&rho;}+W_{0.5})}^2+H_{\mathrm{Re}}^2}}{C}-\frac{H_{\mathrm{Tr}}+\sqrt{{\mathrm{\&rho;}}^2+H_{\mathrm{Re}}^2}}{C}]}\end{array}\right.---\left(3\right)$

其中，PRF为雷达发射天线的脉冲重复频率，V_θ表示场景回波信号接收天线所在平台的飞行速度，ρ表示场景回波信号接收天线所在平台飞行半径，且ρ＝H_Re×tanφ_in(H_Re为场景回波信号接收天线所在平台飞行高度，φ_in为场景回波信号接收天线入射角)，W_0.5＝W_φ/2，W_φ为成像刈宽，，λ为雷达工作波长，C表示电磁波传播速度，H_Tr为发射天线所在平台飞行高度，且H_Tr>H_Re；

步骤S14：根据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线所在平台飞行速度、场景回波信号接收天线入射角和成像刈宽计算圆周运动中心滑动间距，如下所示：

$\frac{\mathrm{\&lambda;}\sqrt{{\mathrm{\&rho;}}^2+H_{\mathrm{Re}}^2}}{2\mathrm{\&rho;}}\&le;\mathrm{\&Delta;\&rho;}\&le;\frac{H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}}{2}---\left(4\right)$

其中，Δρ为圆周运动中心滑动间距，H_Re为场景回波信号接收天线所在平台飞行高度，ρ表示场景回波信号接收天线所在平台飞行半径，且ρ＝H_Re×tanφ_in(φ_in为场景回波信号接收天线入射角)；

步骤S2：滑动式双站圆周合成孔径雷达成像回波数据获取，根据场景回波信号接收天线所在平台自身惯性导航测量系统输出的姿态参数或雷达惯性导航测量系统输出场景回波信号接收天线和直达波信号接收天线的姿态参数和平台位置信息，并控制场景回波信号接收天线对观测场景的回波信号进行接收、雷达发射天线发射具有一定带宽的信号，具体地：

步骤S21：设置雷达数据获取参数，主要包括场景回波信号接收天线航迹向-3dB波束宽度φ_-3dB、航迹向-3dB波束宽度θ_-3dB、雷达发射天线的脉冲重复频率PRF和圆周运动中心滑动间距Δρ；

步骤S22：设置场景回波信号接收天线所在平台飞行轨迹坐标，根据飞行轨迹坐标控制场景回波信号接收天线所在平台进行滑动式圆周飞行，场景回波信号接收天线所在平台飞行轨迹坐标为：

$\mathrm{Track}(x,y,z)\&Element;\left\{\begin{array}{c}x=\mathrm{\&rho;}\&times;\cos \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)+\mathrm{\&Delta;x}+(n-1)\&times;\mathrm{\&Delta;\&rho;}\\ y=\mathrm{\&rho;}\&times;\sin \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)+\mathrm{\&Delta;y}+(n-1)\&times;\mathrm{\&Delta;\&rho;}\\ z=H+\mathrm{\&Delta;z}\end{array}\right.---\left(5\right)$

其中，Track(x,y,z)为场景回波信号接收天线所在平台飞行轨迹函数，(x,y,z)为场景回波信号接收天线所在平台的飞行轨迹坐标，ρ表示场景回波信号接收天线所在平台飞行半径，Δρ为圆周运动中心滑动间距，n表示滑动圆周数，且n为正整数(n＝1,2,3,…,N)，N为总的滑动圆周数，

其中，L_Y为拟观测区域沿Y方向观测范围的有效距离，W_φ为成像刈宽，表示上取整函数；Δx、Δy和Δz分别为沿X、Y和Z方向上所允许最大轨迹偏离值(当Δx、Δy和Δz均不为0时，其运动轨迹与理想圆周运动轨迹存在差异，这里仍然称之为圆周运动)，H为平台飞行高度，α为轨迹角度变量，且α∈[0,2π]，f(α)为滑动圆周函数，

$f\left(\mathrm{\&alpha;}\right)=\left\{\begin{array}{c}1,n=1,\mathrm{\&alpha;}\&Element;\left[\mathrm{0,2}\mathrm{\&pi;}\right]\\ 1,n=2,\&CenterDot;\&CenterDot;\&CenterDot;,N,\mathrm{\&alpha;}\&Element;[0,\frac{\mathrm{\&pi;}}{2}-a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right)]\&cup;[\frac{\mathrm{\&pi;}}{2}+a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right),2\mathrm{\&pi;}]\\ 0,n=2,\&CenterDot;\&CenterDot;\&CenterDot;,N,\mathrm{\&alpha;}\&Element;(\frac{\mathrm{\&pi;}}{2}-a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right),\frac{\mathrm{\&pi;}}{2}+a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right))\end{array}\right.---\left(7\right)$

$\left\{\begin{array}{c}0\&le;\mathrm{\&Delta;x}\&le;0.5\mathrm{\&rho;}\\ 0\&le;\mathrm{\&Delta;y}\&le;0.5\mathrm{\&rho;}\\ 0\&le;\mathrm{\&Delta;z}\&le;0.1H\end{array}\right.---\left(8\right)$

图4给出了滑动式双站圆周合成孔径雷达成像时，场景回波信号接收天线所在平台飞行轨迹示意图，通过控制场景回波信号接收天线所在平台按照平台飞行轨迹函数Track(x,y,z)勾画的轨迹进行飞行，飞行起始点设置在Pos₁(Pos₁为相邻两个圆周运动轨迹在第一或第二坐标象限上的交点)。

步骤S23：场景回波信号接收天线所在平台以观测区域不同中心点进行圆周运动，根据平台自身惯性导航测量系统输出的姿态参数或雷达惯性导航测量系统输出的姿态参数以及平台位置信息进行场景回波信号接收天线入射角调整，其中，惯性导航测量系统输出的姿态参数主要包括航向角、横滚角和俯仰角信息，通过雷达发射天线发射一定带宽的信号，场景回波信号接收天线接收观测场景反射回来的雷达回波信号；直达波信号接收天线需要始终对准雷达发射天线，接收来自雷达发射天线的信号，具体地，如公式(5)、(7)和(8)所示，若n＝1时，即场景回波信号接收天线所在平台以O₁(O₁与O重合)为中心、Pos₁为起始点进行完整的第1周圆周运动；n＝2,…,N时，场景回波信号接收天线所在平台以O_n为中心进行第n周滑动式不完整圆周运动，用于扩大成像观测区域，同时又保持系统性能和分辨率不受损失。

场景回波信号接收天线所在平台飞行轨迹坐标由Track(x,y,z)确定，且飞行过程中根据系统输出的航向角、横滚角和俯仰角信息控制场景回波信号接收天线始终指向以O_n为中心的观测场景m＝n区域，观测场景大小也是以O_n为中心、W_φ/2为半径的圆周区域，场景回波信号接收天线所在平台飞行中，雷达发射按照步骤S21已设置的发射天线的脉冲重复频率PRF以脉冲重复时间PRT为间隔均匀发射信号

s(t)＝exp{j(2πf_ct+πK_rt²)}g(t)            (9)

其中， $g\left(t\right)=\left\{\begin{array}{cc}0,& \left|t\right|>\frac{T}{2}\\ 1,& \left|t\right|\&le;\frac{T}{2}\end{array}\right.,$ K_r为调频率，g(t)＝T为脉冲持续时间，

f_c＝1/λ为雷达工作中心频率，λ为雷达工作波长，场景回波信号接收天线接收观测场景反射回来的雷达回波信号，经下变频后的回波数据E(t,x,y,z,PRT)_q为

$\left\{\begin{array}{c}E{(t,x,y,z,\mathrm{PRT})}_q\\ =\underset{V}{\&Integral;\&Integral;\&Integral;}{\mathrm{\&delta;}}_n(x_n,y_n,z_n)\exp \left\{j\right[-2\mathrm{\&pi;}{f}_c\frac{2R_{\mathrm{ntr}}}{C}+\mathrm{\&pi;}{K}_r{(t-\frac{2R_{\mathrm{ntr}}}{C})}^2\left]\right\}\\ \&times;g(t-\frac{2R_{\mathrm{ntr}}}{C})\mathrm{dxdydz}\end{array}\right.---\left(9\right)$

其中，R_ntr为收发信号传播等效单程距离，且

$\begin{array}{c}{R}_{\mathrm{ntr}}=\frac{\sqrt{{(x_T-x_n)}^2+{(y_T-y_n)}^2+{(z_T-z_n)}^2}}{2}\\ +\frac{\sqrt{{(x-x_n)}^2+{(y-y_n)}^2+{(z-z_n)}^2}}{2}\end{array}---\left(10\right)$

直达波信号接收天线接收来自雷达发射天线的信号，经下变频后的直达波数据E_ref(t,x_ref,y_ref,z_ref,PRT)为

$\begin{array}{c}{E}_{\mathrm{ref}}(t,x_{\mathrm{ref}},y_{\mathrm{ref}},z_{\mathrm{ref}},\mathrm{PRT})\\ =A_{\mathrm{ref}}\exp \left\{j(-2\mathrm{\&pi;}{f}_c\frac{R_{\mathrm{tref}}}{C}+\mathrm{\&pi;}{K}_r{(t-\frac{R_{\mathrm{tref}}}{C})}^2)\right\}g(t-\frac{R_{\mathrm{tref}}}{C})\end{array}---\left(11\right)$

其中，R_tref为直达波信号参考距离，且

$R_{\mathrm{tref}}=\sqrt{{(x_T-x_{\mathrm{ref}})}^2+{(y_T-y_{\mathrm{ref}})}^2+{(z_T-z_{\mathrm{ref}})}^2}---\left(12\right)$

其中，场景回波信号接收天线的三维位置坐标为(x,y,z)，直达波信号接收天线的三维位置坐标为(x_ref,y_ref,z_ref)，雷达发射天线的三维位置坐标为(x_T,y_T,z_T)，V表示单次收发时天线覆盖的观测场景，表示电磁波传播速度，在真空中为3×10⁸m/s。δ_n(x_n,y_n,z_n)为观测场景中目标P_n(x_n,y_n,z_n)的复散射系数，A_ref为直达波信号幅度，R_ntr表示处于(x,y,z)处的场景回波信号接收天线和处于(x_T,y_T,z_T)处的雷达发射天线与目标P_n(x_n,y_n,z_n)间的距离之和，R_tref表示处于(x_ref,y_ref,z_ref)处的直达波信号接收天线和处于(x_T,y_T,z_T)处的雷达发射天线之间的距离，PRT＝1/PRF为脉冲重复时间，q＝1,…,Q_n表示场景回波信号接收天线对应接收回波信号的序列号，第n周圆周运动运动时，整个观测与数据采集过程如图5a所示，场景回波信号接收天线所在平台飞行中的总接收回波数据序列数为Q＝Q₁+Q₂+…+Q_N，ρ表示场景回波信号接收天线所在平台飞行半径，n表示场景回波信号接收天线所在平台第n个飞行圆周轨迹，m表示m个观测场景，n＝1表示场景回波信号接收天线所在平台以O₁(O₁与O重合)为中心、Pos₁为起始点进行完整的第1周圆周运动，对应观测场景区域以O₁为中心；n＝2时，表示以Pos₁为起始点，O₂为中心进行第n周滑动式不完整圆周运动，对应观测场景区域以O₂为中心；n＝3时，表示以Pos₂为起始点，O₃为中心进行第n周双站滑动式不完整圆周运动，对应观测场景区域以O₃为中心；n＝4时，表示以Pos₃为起始点，O₄为中心进行第n周双站滑动式不完整圆周运动，对应观测场景区域以O₄为中心；依次类推。图5b示出了每一周双站滑动式圆周运动的飞行轨迹，所有n周圆周运动的半径均为ρ,当n＝1时，是一个完整的圆周，n>1时为不完整的圆周，图5c示出了每一周滑双站动式圆周运动的观测场景；整个观测场景在水平面(XOY平面)上的投影表示如下：

$\mathrm{Scene}(x_n,y_n)\&Element;\left\{\begin{array}{c}{x}_n=\frac{W_{\mathrm{\&phi;}}}{2}\&times;\cos \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)+(n-1)\&times;\mathrm{\&Delta;\&rho;}\\ y_n=\frac{W_{\mathrm{\&phi;}}}{2}\&times;\sin \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)+(n-1)\&times;\mathrm{\&Delta;\&rho;}\end{array}\right.---\left(13\right)$

相比于现有CSAR成像模式，其整个观测场景在水平面(XOY平面)上的投影表示如下：

$\mathrm{Scene}(x_n,y_n)\&Element;\left\{\begin{array}{c}{x}_n=\frac{W_{\mathrm{\&phi;}}}{2}\&times;\cos \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)\\ y_n=\frac{W_{\mathrm{\&phi;}}}{2}\&times;\sin \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)\end{array}\right.---\left(14\right)$

因此，现有CSAR成像模式与本发明的滑动式成像方法相比，不具备可扩展性，其成像范围也较小；同时，现有CSAR成像模式不能够实现对观测区域的非后向散射信息获取。本发明提出的滑动式双站圆周合成孔径雷达不仅可以实现对观测区域的全方位大范围成像，而且还能获取观测区域的非后向散射信息。

根据系统输出的航向角、横滚角和俯仰角信息控制雷达接收天线始终指向观测场景区域的方式主要有有源相控扫描方式或机械扫描方式两种，目前，这两种方式较为成熟，实际实现也较为容易，不是本发明内容，在此不再敖述，具体控制方式可以参考《雷达手册(第三版)》(Merrill l.Skolnik主编，南京电子技术研究所译，电子工业出版社，2010-7-1出版)。

在步骤S2中，滑动式双站圆周合成孔径雷达成像回波数据获取所需的场景回波信号接收天线所在平台可以为飞机，也可以为飞艇；雷达发射天线所在平台可以为卫星，也可以为飞艇。其中滑动式双站圆周合成孔径雷达成像回波数据获取中，场景回波信号接收天线和直达波信号接收天线之间可以采用刚性基线连接，保持相对静止，刚性基线为平台运动过程中场景回波信号接收天线和直达波信号接收天线之间的三维相对位置变化不超过1mm的金属杆；场景回波信号接收天线和直达波信号接收天线之间也可以直接安装在飞行平台上，并通过激光测量获得两副天线的精确的三维相对位置变化关系，其测量精度需要优于1mm。另外在滑动式双站圆周合成孔径雷达成像回波数据获取中，直达波信号接收天线需要始终对准雷达发射天线，并接收雷达直达波信号。

步骤S3：滑动式双站圆周合成孔径雷达成像数据处理，根据步骤S2获取的滑动式双站圆周合成孔径雷达成像回波数据和直达波数据，通过相干叠加处理获取观测场景的二维或三维图像，具体地：

步骤S31：观测场景对应的图像空间三维离散化，具体地：

以Δx、Δy和Δz的像素间隔分别沿X、Y和Z轴方向对观测场景区域进行三维离散化得到三维图像空间I(n_xΔx,n_yΔy,n_zΔz)，n_x＝1，…，N_x，n_y＝1,…,N_y，n_z＝1,…,N_z，N_x、N_y和N_z分别为沿X、Y和Z轴方向离散化后的像素点数，其中，

$\left\{\begin{array}{cc}\mathrm{\&Delta;x}=\frac{C}{2K_{r}T},& N_x=\frac{W_{\mathrm{\&phi;}}}{\mathrm{\&Delta;x}}\\ \mathrm{\&Delta;y}=\frac{C}{f_c},& N_y=\frac{L_Y}{\mathrm{\&Delta;y}}\\ \mathrm{\&Delta;z}=\frac{C}{f_c},& N_z=\frac{H_z}{\mathrm{\&Delta;z}}\end{array}\right.---\left(15\right)$

步骤S32：距离向傅里叶变换与匹配滤波，具体地：

对步骤S2中平台飞行中所采集到观测场景回波数据E(t,x,y,z,PRT)_q，其中q＝1,...,(Q₁+…+Q_N)，沿距离向进行傅里叶变换，而后进行匹配滤波得到信号E_M(K_ω,x,y,z,PRT)_q，

$\left\{\begin{array}{c}{E}_M{(f_t,x,y,z,\mathrm{PRT})}_q\\ ={\mathrm{FT}}_t\left\{E{(t,x,y,z,\mathrm{PRT})}_q\right\}H\left(f_t\right)\\ \&times;\exp \{-j\frac{2\mathrm{\&pi;}(f_c+f_t)R_{\mathrm{tref}}}{C}\}\end{array}\right.---\left(16\right)$

其中，FT_t表示沿距离向t进行傅里叶变换，H(f_t)为匹配滤波器，

H(f_t)＝conj{FT_t{E_ref(t,x_ref,y_ref,z_ref,PRT)}}      (17)

其中，conj表示复共轭操作。

则匹配滤波得到信号E_M(ft,x,y,z,PRT)_q表示为

$\left\{\begin{array}{c}{E}_M{(f_t,x,y,z,\mathrm{PRT})}_q\\ =\underset{V}{\&Integral;\&Integral;\&Integral;}{\mathrm{\&delta;}}_n(x_n,y_n,z_n)\exp \{-j\frac{4\mathrm{\&pi;}(f_c+f_t)R_{\mathrm{ntr}}}{C}\}\mathrm{rect}\left\{\frac{f_t}{\left|K_r\right|T}\right\}\mathrm{dxdydz}\end{array}\right.---\left(18\right)$

其中， $\mathrm{rect}\{*\}=\left\{\begin{array}{cc}1,& |*|\&le;0.5\\ 0,& |*|>0.5\end{array}\right.,$ V表示单次收发时天线覆盖的观测场景，δ_n(x_n,y_n,z_n)为观测场景中目标P_n(x_n,y_n,z_n)的复散射系数，t_n表示处于(x,y,z)处的场景回波信号接收天线和雷达发射天线距离目标P_n(x_n,y_n,z_n)之间的电磁波传播时间，R_n为对应(x,y,z)处的场景回波信号接收天线和(x_T,y_T,z_T)处的雷达发射天线与目标P_n(x_n,y_n,z_n)之间距离，(x,y,z)为场景回波信号接收天线所在平台运动时的飞行轨迹坐标；

记K_ω＝2π(f_c+f_t)/C，f_c为雷达工作中心频率，C为电磁波传播速度，式(18)记为

$\left\{\begin{array}{c}{E}_M{(K_{\mathrm{\&omega;}},x,y,z,\mathrm{PRT})}_q=E_M{(f_t,x,y,z,\mathrm{PRT})}_q\\ =\underset{V}{\&Integral;\&Integral;\&Integral;}{\mathrm{\&delta;}}_n(x_n,y_n,z_n)\exp \{-j2K_{\mathrm{\&omega;}}{R}_{\mathrm{ntr}}\}\mathrm{dxdydz}\end{array}\right.---\left(19\right)$

步骤S33：滑动圆周相干叠加，通过对步骤S32得到的距离向傅里叶变换与匹配滤波结果进行相干叠加循环求解图像的每个像素值，具体地：

步骤S330：取n_x＝1，n_y＝1，n_z＝1，其中n_x，n_y和n_z表示像素序号，且n_x＝1，…，N_x，n_y＝1,…,N_y，n_z＝1,…,N_z，N_x、N_y和N_z分别为沿X、Y和Z轴方向离散化后的像素点数；

步骤S331：计算图像I(n_xΔx,n_yΔy,n_zΔz)中第(n_x,n_y,n_z)像素对应的坐标位置(n_xΔx,n_yΔy,n_zΔz)到场景回波信号接收天线(x,y,z)和雷达发射天线(x_T,y_T,z_T)获取雷达回波数据时对应每个位置(x,y,z)的距离，并生成匹配函数，

$\begin{array}{c}{H}_M(K_{\mathrm{\&omega;}},x,y,z,n_x\mathrm{\&Delta;x},n_y\mathrm{\&Delta;y},n_z\mathrm{\&Delta;z})=\exp \left\{j2K_{\mathrm{\&omega;}}{R}_{I-(x,y,z)}\right\}\\ =\exp \left\{j2K_{\mathrm{\&omega;}}\right[\frac{\sqrt{{(x_T-x_n)}^2+{(y_T-y_n)}^2+{(z_T-z_n)}^2}}{2}\\ +\frac{\sqrt{{(x-n_x\mathrm{\&Delta;x})}^2+{(y-n_y\mathrm{\&Delta;y})}^2+{(z-n_z\mathrm{\&Delta;z})}^2}}{2}\left]\right\}\end{array}---\left(20\right)$

其中，Δx、Δy和Δz为待计算图像I(n_xΔx,n_yΔy,n_zΔz)的像素间隔，

n_x＝1，…，N_x，n_y＝1,…,N_y，n_z＝1,…,N_z，N_x、N_y和N_z分别为沿X、Y和Z轴方向离散化后的像素点数；

步骤S332：求解图像的第(n_x,n_y,n_z)个像素对应的数值，具体地：

$\begin{array}{c}I(n_x\mathrm{\&Delta;x},n_y\mathrm{\&Delta;y},n_z\mathrm{\&Delta;z})\\ =\underset{q=1}{\overset{Q}{\mathrm{\&Sigma;}}}\underset{K_{\mathrm{\&omega;}}}{\&Integral;}{H}_M(K_{\mathrm{\&omega;}},x,y,z,n_x,\mathrm{\&Delta;x},n_y\mathrm{\&Delta;y},n_z\mathrm{\&Delta;z})\&times;E_M{(K_{\mathrm{\&omega;}},x,y,z,\mathrm{PRT})}_q\end{array}---\left(21\right)$ 其中，Q＝Q₁+Q₂+…+Q_N为雷达总接收回波数据序列数，表示沿K_ω方向对*进行积分，表示对*进行相干求和；

步骤S333：若n_x<N_x，使n_x加1，并返回至步骤S331，若n_x≥N_x，继续执行步骤S334；

步骤S333：若n_x<N_x，使n_x加1，并返回至步骤S331，若n_x≥N_x，继续执行步骤S334；

步骤S334：若n_y<N_y，使n_y加1，并返回至步骤S331，若n_y≥N_y，继续执行步骤S335；

步骤S335：若n_z<N_z，使n_z加1，并返回至步骤S331，若n_z≥N_z，继续执行步骤S34。

步骤S335，也可以只针对1～N_z某一个像素n_z求解获得观测场景的二维图像，即

I₂(n_xΔx,n_yΔy)＝I(n_xΔx,n_yΔy,n_z＝N_nz)

其中，I₂(n_xΔx,n_yΔy)表示观测区域的二维图像，Δx、Δy和Δz为待计算图像I(n_xΔx,n_yΔy,n_zΔz)的像素间隔，n_x＝1，…，N_x，n_y＝1,…,N_y，n_z＝1,…,N_z，N_x、N_y和N_z分别为沿X、Y和Z轴方向离散化后的像素点数，N_nz取1～N_z某一个数值；

步骤S34：输出滑动圆周合成孔径雷达成像三维图像I(n_xΔx,n_yΔy,n_zΔz)。

综上，滑动式双站圆周合成孔径雷达成像方法的整个流程如图6所示。

本发明给出了一种滑动式双站圆周合成孔径雷达成像方法，实现了大范围观测区域成像，不仅能够实现雷达随平台的运动对目标进行全方位观测，而且还能实现保留常规CSAR成像的优势；解决了现有CSAR成像区域受雷达天线俯仰向和航迹向-3dB波束宽度限制的问题，实现了大范围高分辨率全方位圆周成像数据获取与成像处理；同时还能实现非后向散射特性获取。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种滑动式双站圆周合成孔径雷达成像方法，该方法包括步骤：

步骤1，计算滑动式双站圆周合成孔径雷达成像时的雷达接收天线俯仰向和航迹向-3dB波束宽度、雷达发射天线发射信号的脉冲重复频率和场景回波信号接收天线所在平台的圆周运动中心滑动间距；

步骤2：根据场景回波信号接收天线所在平台自身惯性导航测量系统输出的姿态参数或雷达惯性导航测量系统输出场景回波信号接收天线和直达波信号接收天线的姿态参数和平台位置信息，并控制场景回波信号接收天线对观测场景的回波信号进行接收，雷达发射天线发射具有一定带宽的信号；

步骤S3：根据步骤S2获取的滑动式双站圆周合成孔径雷达成像回波数据和直达波数据，通过相干叠加处理获取观测场景的二维或三维图像。

2.根据权利要求1所述的方法，其特征在于，步骤1进一步包括：

步骤S11，根据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线入射角和成像刈宽计算航迹向-3dB波束宽度；

步骤S12，根据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线入射角和成像刈宽计算航迹向-3dB波束宽度；

步骤S13，据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线所在平台飞行速度、场景回波信号接收天线入射角和成像刈宽计算雷达发射天线的脉冲重复频率；

步骤S14，根据场景回波信号接收天线所在平台飞行高度、场景回波信号接收天线所在平台飞行速度、场景回波信号接收天线入射角和成像刈宽计算圆周运动中心滑动间距。

3.根据权利要求2所述的方法，其特征在于，所述俯仰向-3dB波束宽度按照下式进行计算：

${\mathrm{\&phi;}}_{-3\mathrm{dB}}=\arctan \left(\frac{2H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}+W_{\mathrm{\&phi;}}}{2H_{\mathrm{Re}}}\right)-\arctan \left(\frac{2H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}-W_{\mathrm{\&phi;}}}{2H_{\mathrm{Re}}}\right)$

其中，φ_-3dB为场景回波信号接收天线入射角俯仰向-3dB波束宽度，H_Re为场景回波信号接收天线所在平台飞行高度，φ_in为场景回波信号接收天线入射角，W_φ为成像刈宽。

4.根据权利要求3所述的方法，其特征在于，所述航迹向-3dB波束宽度按照下式进行计算：

${\mathrm{\&theta;}}_{-3\mathrm{dB}}\&GreaterEqual;{\mathrm{\&phi;}}_{-3\mathrm{dB}}=\arctan \left(\frac{2H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}+W_{\mathrm{\&phi;}}}{2H_{\mathrm{Re}}}\right)-\arctan \left(\frac{2H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}-W_{\mathrm{\&phi;}}}{2H_{\mathrm{Re}}}\right)$

其中，θ_-3dB为场景回波信号接收天线航迹向-3dB波束宽度，φ_-3dB为场景回波信号接收天线俯仰向-3dB波束宽度，H_Re为场景回波信号接收天线所在平台飞行高度，φ_in为场景回波信号接收天线入射角，W_φ为成像刈宽，θ_-3dB≥φ_-3dB表示系统工作时，保证雷达场景回波信号接收天线航迹向-3dB波束宽度大于等于雷达接收天线俯仰向-3dB波束宽度，以保证滑动式双站圆周合成孔径雷达成像时成像刈宽最大限度地由场景回波信号接收天线俯仰向-3dB波束宽度决定，不降低雷达系统的成像刈宽。

5.根据权利要求4所述的方法，其特征在于，雷达发射天线的脉冲重复频率通过下式进行计算：

$\left\{\begin{array}{c}\mathrm{PRF}\&GreaterEqual;\frac{V_{\mathrm{\&theta;}}(\sqrt{A+{\mathrm{\&rho;W}}_{\mathrm{\&phi;}}}-\sqrt{A-{\mathrm{\&rho;W}}_{\mathrm{\&phi;}}})}{\mathrm{\&lambda;\&rho;}}\\ \mathrm{PRF}\&le;\frac{1}{2[\frac{\sqrt{W_{0.5}^2+H_{\mathrm{Tr}}^2}+\sqrt{{(\mathrm{\&rho;}+W_{0.5})}^2+H_{\mathrm{Re}}^2}}{C}-\frac{H_{\mathrm{Tr}}+\sqrt{{\mathrm{\&rho;}}^2+H_{\mathrm{Re}}^2}}{C}]}\end{array}\right.$

其中，PRF为雷达发射天线的脉冲重复频率，V_θ表示场景回波信号接收天线所在平台的飞行速度，ρ表示场景回波信号接收天线所在平台飞行半径，且ρ＝H_Re×tanφ_in(H_Re为场景回波信号接收天线所在平台飞行高度，φ_in为场景回波信号接收天线入射角)，W_0.5＝W_φ/2，W_φ为成像刈宽，，λ为雷达工作波长，C表示电磁波传播速度，H_Tr为发射天线所在平台飞行高度，且H_Tr>H_Re。

6.根据权利要求5所述的方法，其特征在于，圆周运动中心滑动间距通过下式进行计算：

$\frac{\mathrm{\&lambda;}\sqrt{{\mathrm{\&rho;}}^2+H_{\mathrm{Re}}^2}}{2\mathrm{\&rho;}}\&le;\mathrm{\&Delta;\&rho;}\&le;\frac{H_{\mathrm{Re}}\&times;\tan {\mathrm{\&phi;}}_{\mathrm{in}}}{2}$

其中，Δρ为圆周运动中心滑动间距，H_Re为场景回波信号接收天线所在平台飞行高度，ρ表示场景回波信号接收天线所在平台飞行半径，且ρ＝H_Re×tanφ_in，φ_in为场景回波信号接收天线入射角。

7.根据权利要求6所述的方法，其特征在于，步骤2进一步包括：

步骤S21，设置雷达数据获取参数，包括场景回波信号接收天线航迹向-3dB波束宽度φ_-3dB、航迹向-3dB波束宽度θ_-3dB、雷达发射天线的脉冲重复频率PRF和圆周运动中心滑动间距Δρ；

步骤S22，设置场景回波信号接收天线所在平台飞行轨迹坐标，根据飞行轨迹坐标控制场景回波信号接收天线所在平台进行滑动式圆周飞行；

步骤S23：场景回波信号接收天线所在平台以观测区域不同中心点进行圆周运动，根据平台自身惯性导航测量系统输出的姿态参数或雷达惯性导航测量系统输出的姿态参数以及平台位置信息进行场景回波信号接收天线入射角调整，其中，惯性导航测量系统输出的姿态参数包括航向角、横滚角和俯仰角信息，通过雷达发射天线发射一定带宽的信号，场景回波信号接收天线接收观测场景反射回来的雷达回波信号，直达波信号接收天线需要始终对准雷达发射天线，接收来自雷达发射天线的信号。

8.根据权利要求7所述的方法，其特征在于，所述场景回波信号接收天线所在平台飞行轨迹坐标为：

$\mathrm{Track}(x,y,z)\&Element;\left\{\begin{array}{c}x=\mathrm{\&rho;}\&times;\cos \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)+\mathrm{\&Delta;x}+(n-1)\&times;\mathrm{\&Delta;\&rho;}\\ y=\mathrm{\&rho;}\&times;\sin \mathrm{\&alpha;}\&times;f\left(\mathrm{\&alpha;}\right)+\mathrm{\&Delta;y}+(n-1)\&times;\mathrm{\&Delta;\&rho;}\\ z=H+\mathrm{\&Delta;z}\end{array}\right.$

其中，Track(x,y,z)为场景回波信号接收天线所在平台飞行轨迹函数，(x,y,z)为场景回波信号接收天线所在平台的飞行轨迹坐标，ρ表示场景回波信号接收天线所在平台飞行半径，Δρ为圆周运动中心滑动间距，n表示滑动圆周数，且n为正整数，n＝1,2,3,…,N，N为总的滑动圆周数，

其中，L_Y为拟观测区域沿Y方向观测范围的有效距离，W_φ为成像刈宽，表示上取整函数；Δx、Δy和Δz分别为沿X、Y和Z方向上所允许最大轨迹偏离值，H为飞行高度，α为轨迹角度变量，且α∈[0,2π]，f(α)为滑动圆周函数，

$f\left(\mathrm{\&alpha;}\right)=\left\{\begin{array}{c}1,n=1,\mathrm{\&alpha;}\&Element;\left[\mathrm{0,2}\mathrm{\&pi;}\right]\\ 1,n=2,\&CenterDot;\&CenterDot;\&CenterDot;,N,\mathrm{\&alpha;}\&Element;[0,\frac{\mathrm{\&pi;}}{2}-a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right)]\&cup;[\frac{\mathrm{\&pi;}}{2}+a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right),2\mathrm{\&pi;}]\\ 0,n=2,\&CenterDot;\&CenterDot;\&CenterDot;,N,\mathrm{\&alpha;}\&Element;(\frac{\mathrm{\&pi;}}{2}-a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right),\frac{\mathrm{\&pi;}}{2}+a\tan \left(\frac{\mathrm{\&Delta;\&rho;}}{2\mathrm{\&rho;}}\right))\end{array}\right.$

$\left\{\begin{array}{c}0\&le;\mathrm{\&Delta;x}\&le;0.5\mathrm{\&rho;}\\ 0\&le;\mathrm{\&Delta;y}\&le;0.5\mathrm{\&rho;}\\ 0\&le;\mathrm{\&Delta;z}\&le;0.1H\end{array}\right..$

9.根据权利要求8所述的方法，其特征在于，步骤S3进一步包括：

步骤S31：观测场景对应的图像空间三维离散化；

步骤S32：对步骤S2中平台飞行中所采集到观测场景回波数据，沿距离向进行傅里叶变换，而后进行匹配滤波；

步骤S33：滑动圆周相干叠加，通过对步骤S32得到的距离向傅里叶变换与匹配滤波结果进行相干叠加循环求解图像的每个像素值；

步骤S34：输出滑动圆周合成孔径雷达成像三维图像。