一种SAR成像方法、装置及系统[发明专利]

(10)申请公布号 (43)申请公布日 2013.12.25C N 103472448 A (21)申请号 201310296193.5
(22)申请日 2013.07.15
G01S 13/90(2006.01)
(71)申请人中国科学院电子学研究所
地址100190 北京市海淀区北四环西路19
号
(72)发明人刘亚波 宋明聪 郑世超 王宇
赵凤军 邓云凯
(74)专利代理机构北京派特恩知识产权代理事
务所(普通合伙) 11270
代理人张振伟
王黎延
(54)发明名称
一种SAR 成像方法、装置及系统
(57)摘要
本发明公开了一种SAR 成像方法，包括：将接
收到的第一全孔径SAR 回波数据沿方位向划分为
两段以上子孔径SAR 回波数据，对所述两段以上
子孔径SAR 回波数据进行距离向处理；将距离向
处理后的所述两段以上子孔径SAR 回波数据合并
为第二全孔径SAR 回波数据；将所述第二全孔径
SAR 回波数据沿距离向划分为两个以上SAR 回波
数据块，对所述两个以上SAR 回波数据块进行方
位向处理；将方位向处理后的所述两个以上SAR
回波数据块合并为SAR 图像。

本发明还同时公开
了一种SAR 成像装置及系统。

采用本发明的方案，
实现了对大规模SAR 回波数据进行快速、高质量
的成像处理，具有良好的成像效果和GPU 加速比。

(51)Int.Cl.
权利要求书3页 说明书9页 附图5页
(19)中华人民共和国国家知识产权局(12)发明专利申请权利要求书3页 说明书9页 附图5页(10)申请公布号CN 103472448 A
*CN103472448A*
1.一种合成孔径雷达SAR成像方法，其特征在于，所述方法包括：
将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据，对所述两段以上子孔径SAR回波数据进行距离向处理；
将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据；
将所述第二全孔径SAR回波数据沿距离向划分为两个以上SAR回波数据块，对所述两个以上SAR回波数据块进行方位向处理；
将方位向处理后的所述两个以上SAR回波数据块合并为SAR图像。

2.根据权利要求1所述的方法，其特征在于，所述将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据，包括：
根据预设的时间周期T
sub
，将接收到的第一全孔径SAR回波数据按接收时间先后顺序
沿方位向划分为两段以上子孔径SAR回波数据；其中，所述T
sub 满足T
sub
×B
sub
＞＞1，所述
B
sub
为所述子孔径SAR回波数据对应的带宽，所述＞＞为远大于。

3.根据权利要求1所述的方法，其特征在于，所述对所述两段以上子孔径SAR回波数据进行距离向处理，包括：
按照接收时间先后顺序分别对所述两段以上子孔径SAR回波数据进行方位向傅里叶变换；
在方位向傅里叶变换后的所述子孔径SAR回波数据中插入线性调频尺度CS因子；
对插入CS因子后的所述子孔径SAR回波数据进行距离向傅里叶变换；
对距离向傅里叶变后的所述子孔径SAR回波数据进行距离压缩；
对距离压缩后的所述子孔径SAR回波数据中的残余距离二次距离压缩SRC项和一致距离徙动修正RCMC项进行矫正；
对SRC项和一致RCMC项进行矫正后的所述子孔径SAR回波数据进行距离向傅里叶逆变换；
对距离向傅里叶逆变换的所述子孔径SAR回波数据进行方位向傅里叶逆变换。

4.根据权利要求1所述的方法，其特征在于，所述将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据，包括：
获取相邻两段所述子孔径SAR回波数据的重叠率；
根据所述重叠率将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据。

5.根据权利要求1所述的方法，其特征在于，所述对所述两个以上SAR回波数据块进行方位向处理之前，所述方法还包括：
分别对所述两个以上SAR回波数据块进行转置。

6.根据权利要求1至5任一项所述的方法，其特征在于，所述对所述两个以上SAR回波数据块进行方位向处理，包括：
按照接收的时间先后顺序分别对所述两个以上SAR回波数据块进行方位向傅里叶变换；
对方位向傅里叶变换后的所述SAR回波数据块进行方位压缩和相位校正；
对方位压缩和相位校正后的所述SAR回波数据块进行方位向傅里叶逆变换。

7.一种SAR成像装置，其特征在于，所述装置包括：第一划分单元、第一合并单元、第二划分单元以及第二合并单元；其中，
所述第一划分单元，用于将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据并发送至图形处理器GPU；
所述第一合并单元，用于接收所述GPU发送的距离向处理后的所述两段以上子孔径SAR回波数据，并将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据；
所述第二划分单元，用于将所述第二全孔径SAR回波数据沿距离向划分为两个以上SAR回波数据块并发送至所述GPU；
所述第二合并单元，用于接收所述GPU发送的方位向处理后的所述两个以上SAR回波数据块，并将方位向处理后的所述两个以上SAR回波数据块合并为SAR图像。

8.根据权利要求7所述的装置，其特征在于，所述第一划分单元，还用于根据预设的时间周期T
sub
，将接收到的第一全孔径SAR回波数据按接收先后顺序沿方位向划分为两段以
上子孔径SAR回波数据；其中，所述T
sub 满足T
sub
×B
sub
＞＞1，所述B
sub
为所述子孔径SAR回
波数据对应的带宽，所述＞＞为远大于。

9.根据权利要求7所述的装置，其特征在于，所述第一合并单元包括：获取子单元以及合并子单元；其中，
所述获取子单元，用于获取相邻两段所述子孔径SAR回波数据的重叠率；
所述合并子单元，用于根据所述重叠率将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据。

10.根据权利要求7至9任一项所述的装置，其特征在于，所述装置还包括转置单元，用于分别对所述两个以上SAR回波数据块进行转置。

11.一种SAR成像装置，其特征在于，所述装置包括：第一处理单元以及第二处理单元；其中，
所述第一处理单，用于接收中央处理器CPU发送的两段以上子孔径SAR回波数据，并对所述两段以上子孔径SAR回波数据进行距离向处理；将距离向处理后的所述两段以上子孔径SAR回波数据发送至所述CPU；
所述第二处理单元，用于接收CPU发送的两个以上SAR回波数据块，并对所述两个以上SAR回波数据块进行方位向处理；将方位向处理后的所述两个以上SAR回波数据块发送至所述CPU。

12.根据权利要求11所述的装置，其特征在于，所述第一处理单元包括：第一处理子单元、第二处理子单元、第三处理子单元、第四处理子单元、第五处理子单元、第六处理子单元以及第七处理子单元；其中，
所述第一处理子单元，用于按接收先后顺序分别对所述两段以上子孔径SAR回波数据进行方位向傅里叶变换；
所述第二处理子单元，用于在方位向傅里叶变换后的所述子孔径SAR回波数据中插入CS因子；
所述第三处理子单元，用于对插入CS因子后的所述子孔径SAR回波数据进行距离向傅里叶变换；
所述第四处理子单元，用于对距离向傅里叶变后的所述子孔径SAR回波数据进行距离压缩；
所述第五处理子单元，用于对距离压缩后的所述子孔径SAR回波数据中的SRC项和一致RCMC项进行矫正；
所述第六处理子单元，用于对SRC项和一致RCMC项进行矫正后的所述子孔径SAR回波数据进行距离向傅里叶逆变换；
所述第七处理子单元，用于对距离向傅里叶逆变换的所述子孔径SAR回波数据进行方位向傅里叶逆变换。

13.根据权利要求11或12所述的装置，其特征在于，所述第二处理单元包括：第八处理子单元、第九处理子单元以及第十处理子单元；其中，
所述第八处理子单元，用于按接收先后顺序分别对所述两个以上SAR回波数据块进行方位向傅里叶变换；
所述第九处理子单元，用于对方位向傅里叶变换后的所述SAR回波数据块进行方位压缩和相位校正；
所述第十处理子单元，用于对方位压缩和相位校正后的所述SAR回波数据块进行方位向傅里叶逆变换。

14.一种SAR成像系统，其特征在于，所述系统包括权利要求7至10任一项所述的装置、以及权利要求11至13任一项所述的装置。

一种SAR成像方法、装置及系统
技术领域
[0001] 本发明涉及合成孔径雷达（SAR，Synthetic Aperture Radar）成像技术，尤其涉及一种SAR成像方法、装置及系统。

背景技术
[0002] SAR是一种主动式的对地观测技术，与传统的光学传感器相比，SAR可以实现全天时、全天候对地实时观测，并具有一定的地表穿透能力。

SAR技术的发展带来了数据量和计算量的急剧膨胀，而且随着实时成像需求的增强，对大规模SAR回波数据的实时成像处理越发重要。

[0003] 近年来，随着图形处理器（GPU，Graphic Processing Unit）并行处理能力的不断提高，GPU正被越来越多地应用于科学计算领域，例如图像分割，SAR回波仿真等。

在SAR成像处理领域，由于SAR成像处理算法主要由向量相乘组成，非常利于并行实现，所以GPU加速可以发挥巨大作用。

对于小于GPU内存的SAR回波数据，GPU加速可以取得非常理想的加速比。

但随着SAR回波数据的增大，受限于GPU显存的容量，GPU不能一次完成整个场景的SAR回波处理。

如果运用传统的成像算法，需要不停地在中央处理器（CPU，Central Processing Unit）和GPU之间进行大规模数据的传输，而新一代总线接口（PCI-E，Peripheral Component Interconnect-Express）的传输带宽相对于GPU的内存读取带宽来说非常小，因此这将严重制约GPU的处理速度，所以对于大规模的SAR回波处理，仅仅将原有成像算法在GPU上并行实现是远远不够的。

[0004] 综上所述，如何对大规模的SAR回波数据进行快速、高质量的成像处理是亟待解决的问题。

发明内容
[0005] 有鉴于此，本发明的主要目的在于提供一种SAR成像方法、装置及系统，能够实现对大规模SAR回波数据进行快速、高质量的成像处理，具有良好的成像效果和GPU加速比。

[0006] 为达到上述目的，本发明的技术方案是这样实现的：
[0007] 一种合成孔径雷达SAR成像方法，所述方法包括：
[0008] 将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据，对所述两段以上子孔径SAR回波数据进行距离向处理；
[0009] 将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据；
[0010] 将所述第二全孔径SAR回波数据沿距离向划分为两个以上SAR回波数据块，对所述两个以上SAR回波数据块进行方位向处理；
[0011] 将方位向处理后的所述两个以上SAR回波数据块合并为SAR图像。

[0012] 其中，所述将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据，包括：
[0013] 根据预设的时间周期T sub，将接收到的第一全孔径SAR回波数据按接收时间先后
顺序沿方位向划分为两段以上子孔径SAR回波数据；其中，所述T
sub 满足T
sub
×B
sub
＞＞1，
所述B
sub
为所述子孔径SAR回波数据对应的带宽，所述＞＞为远大于。

[0014] 其中，所述对所述两段以上子孔径SAR回波数据进行距离向处理，包括：[0015] 按照接收时间先后顺序分别对所述两段以上子孔径SAR回波数据进行方位向傅里叶变换；
[0016] 在方位向傅里叶变换后的所述子孔径SAR回波数据中插入线性调频尺度CS因子；
[0017] 对插入CS因子后的所述子孔径SAR回波数据进行距离向傅里叶变换；[0018] 对距离向傅里叶变后的所述子孔径SAR回波数据进行距离压缩；
[0019] 对距离压缩后的所述子孔径SAR回波数据中的残余距离二次距离压缩SRC项和一致距离徙动修正RCMC项进行矫正；
[0020] 对SRC项和一致RCMC项进行矫正后的所述子孔径SAR回波数据进行距离向傅里叶逆变换；
[0021] 对距离向傅里叶逆变换的所述子孔径SAR回波数据进行方位向傅里叶逆变换。

[0022] 其中，所述将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据，包括：
[0023] 获取相邻两段所述子孔径SAR回波数据的重叠率；
[0024] 根据所述重叠率将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据。

[0025] 其中，所述对所述两个以上SAR回波数据块进行方位向处理之前，所述方法还包括：
[0026] 分别对所述两个以上SAR回波数据块进行转置。

[0027] 其中，所述对所述两个以上SAR回波数据块进行方位向处理，包括：
[0028] 按照接收的时间先后顺序分别对所述两个以上SAR回波数据块进行方位向傅里叶变换；
[0029] 对方位向傅里叶变换后的所述SAR回波数据块进行方位压缩和相位校正；[0030] 对方位压缩和相位校正后的所述SAR回波数据块进行方位向傅里叶逆变换。

[0031] 一种SAR成像装置，所述装置包括：第一划分单元、第一合并单元、第二划分单元以及第二合并单元；其中，
[0032] 所述第一划分单元，用于将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据并发送至图形处理器GPU；
[0033] 所述第一合并单元，用于接收所述GPU发送的距离向处理后的所述两段以上子孔径SAR回波数据，并将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据；
[0034] 所述第二划分单元，用于将所述第二全孔径SAR回波数据沿距离向划分为两个以上SAR回波数据块并发送至所述GPU；
[0035] 所述第二合并单元，用于接收所述GPU发送的方位向处理后的所述两个以上SAR 回波数据块，并将方位向处理后的所述两个以上SAR回波数据块合并为SAR图像。

[0036] 其中，所述第一划分单元，还用于根据预设的时间周期T sub，将接收到的第一全孔径SAR回波数据按接收先后顺序沿方位向划分为两段以上子孔径SAR回波数据；其中，所述
T sub 满足T
sub
×B
sub
＞＞1，所述B
sub
为所述子孔径SAR回波数据对应的带宽，所述＞＞为远
大于。

[0037] 其中，所述第一合并单元包括：获取子单元以及合并子单元；其中，
[0038] 所述获取子单元，用于获取相邻两段所述子孔径SAR回波数据的重叠率；[0039] 所述合并子单元，用于根据所述重叠率将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据。

[0040] 其中，所述装置还包括转置单元，用于分别对所述两个以上SAR回波数据块进行转置。

[0041] 一种SAR成像装置，所述装置包括：第一处理单元以及第二处理单元；其中，[0042] 所述第一处理单，用于接收中央处理器CPU发送的两段以上子孔径SAR回波数据，并对所述两段以上子孔径SAR回波数据进行距离向处理；将距离向处理后的所述两段以上子孔径SAR回波数据发送至所述CPU；
[0043] 所述第二处理单元，用于接收CPU发送的两个以上SAR回波数据块，并对所述两个以上SAR回波数据块进行方位向处理；将方位向处理后的所述两个以上SAR回波数据块发送至所述CPU。

[0044] 其中，所述第一处理单元包括：第一处理子单元、第二处理子单元、第三处理子单元、第四处理子单元、第五处理子单元、第六处理子单元以及第七处理子单元；其中，[0045] 所述第一处理子单元，用于按接收先后顺序分别对所述两段以上子孔径SAR回波数据进行方位向傅里叶变换；
[0046] 所述第二处理子单元，用于在方位向傅里叶变换后的所述子孔径SAR回波数据中插入CS因子；
[0047] 所述第三处理子单元，用于对插入CS因子后的所述子孔径SAR回波数据进行距离向傅里叶变换；
[0048] 所述第四处理子单元，用于对距离向傅里叶变后的所述子孔径SAR回波数据进行距离压缩；
[0049] 所述第五处理子单元，用于对距离压缩后的所述子孔径SAR回波数据中的SRC项和一致RCMC项进行矫正；
[0050] 所述第六处理子单元，用于对SRC项和一致RCMC项进行矫正后的所述子孔径SAR 回波数据进行距离向傅里叶逆变换；
[0051] 所述第七处理子单元，用于对距离向傅里叶逆变换的所述子孔径SAR回波数据进行方位向傅里叶逆变换。

[0052] 其中，所述第二处理单元包括：第八处理子单元、第九处理子单元以及第十处理子单元；其中，
[0053] 所述第八处理子单元，用于按接收先后顺序分别对所述两个以上SAR回波数据块进行方位向傅里叶变换；
[0054] 所述第九处理子单元，用于对方位向傅里叶变换后的所述SAR回波数据块进行方位压缩和相位校正；
[0055] 所述第十处理子单元，用于对方位压缩和相位校正后的所述SAR回波数据块进行方位向傅里叶逆变换。

[0056] 一种SAR成像系统，所述系统包括上述任意SAR成像装置。

[0057] 本发明记载的SAR成像方法、装置及系统，将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据，对所述两段以上子孔径SAR回波数据进行距离向处理；将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据；将所述第二全孔径SAR回波数据沿距离向划分为两个以上SAR回波数据块，对所述两个以上SAR回波数据块进行方位向处理；将方位向处理后的所述两个以上SAR回波数据块合并为SAR图像。

如此，可以对大规模的SAR回波数据进行划分，分别对划分的SAR回波数据进行成像处理，提高了SAR回波数据处理的速度和质量，进而实现了良好的成像效果和GPU加速比。

附图说明
[0058] 图1为本发明实施例SAR成像方法的实现流程示意图一；
[0059] 图2为本发明实施例SAR成像装置的结构组成示意图一；
[0060] 图3为本发明实施例SAR成像装置的结构组成示意图二；
[0061] 图4为本发明实施例SAR成像系统的结构组成示意图；
[0062] 图5为本发明实施例SAR成像方法的实现流程示意图二；
[0063] 图6为本发明实施例SAR成像处理时序图；
[0064] 图7为本发明实施例三种算法的加速比对比示意图。

具体实施方式
[0065] 为了能够更加详尽地了解本发明的特点与技术内容，下面结合附图对本发明的实现进行详细阐述，所附附图仅供参考说明之用，并非用来限定本发明。

[0066] 本发明实施例记载了一种SAR成像方法，如图1所示，该方法包括以下步骤：[0067] 步骤101：将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据，对所述两段以上子孔径SAR回波数据进行距离向处理。

[0068] 这里，第一全孔径SAR回波数据是指由地面返回的全孔径SAR回波数据。

[0069] 优选地，所述将接收到的第一全孔径SAR回波数据沿方位向划分为两段以上子孔径SAR回波数据，包括：
[0070] 根据预设的时间周期T sub，将接收到的第一全孔径SAR回波数据按接收先后顺序
沿方位向划分为两段以上子孔径SAR回波数据；其中，所述T
sub 满足 T
sub
×B
sub
＞＞1，所述
B
sub
为所述子孔径SAR回波数据对应的带宽，所述＞＞为远大于。

[0071] 具体地，对地面返回的全孔径SAR回波数据进行子孔径划分时，需要满足每个子孔径的时间带宽积远大于1，子孔径的时间带宽积为（1）式，
[0072]
[0073] 因此，当划分子孔径时，时间带宽积需要满足（2）式，
[0074]
[0075] 其中，R0s是最短距离，V r是SAR系统沿飞行路径的速度，λ是SAR信号波长，PRF 是划分的子孔径长度。

是脉冲重复频率，N
sub
[0076] 如此，在接下来的距离向处理之前，第一全孔径SAR回波数据就已经沿方位向划分成不同的子孔径，如果划分的子孔径的大小小于GPU的显存，GPU就可以以较高的加速比处理子孔径对应的SAR回波数据；优选地，可以将多段子孔径SAR回波数据分别交给不同的GPU去做并行处理。

此外，由于对子孔径SAR回波数据进行处理，因此，不必等到所有的SAR 回波数据全部采集到再进行子孔径SAR回波数据的距离向处理，当一段子孔径SAR回波数据的采集完成后，就可以进行该子孔径SAR回波数据的距离向处理，如此，SAR回波数据的采集和子孔径SAR回波数据的距离向处理形成链接，在SAR回波数据的采集过程中也完成了子孔径SAR回波数据的距离向处理。

[0077] 优选地，所述对所述两段以上子孔径SAR回波数据进行距离向处理，包括：[0078] 按接收先后顺序分别对所述两段以上子孔径SAR回波数据进行方位向傅里叶变换；
[0079] 在方位向傅里叶变换后的所述子孔径SAR回波数据中插入线性调频尺度（CS，Chirp Scaling）因子；
[0080] 对插入CS因子后的所述子孔径SAR回波数据进行距离向傅里叶变换；[0081] 对距离向傅里叶变后的所述子孔径SAR回波数据进行距离压缩；
[0082] 对距离压缩后的所述子孔径SAR回波数据中的残余距离二次距离压缩（SRC，Secondary Range Compression）项和一致距离徙动修正（RCMC，Range Cell Migration Correction）项进行矫正；
[0083] 对SRC项和一致RCMC项进行矫正后的所述子孔径SAR回波数据进行距离向傅里叶逆变换；
[0084] 对距离向傅里叶逆变换的所述子孔径SAR回波数据进行方位向傅里叶逆变换。

[0085] 这里，所述CS因子根据所述SAR系统的各项参数设定，所述SRC项和一致RCMC项为根据CS算法得到的内容，均为现有技术，此处不再赘述。

[0086] 步骤102：将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据。

[0087] 这里，所述第二全孔径SAR回波数据是指将距离向处理后的所述两段以上子孔径SAR回波数据合并后的全孔径SAR回波数据。

优选地，合并时，需要保证第二全孔径SAR回波数据的频谱的光滑度，如此，最终的成像质量可以不受上述子孔径划分的影响。

[0088] 具体地，由于SAR斜视角和SAR天线波束宽度的影响，子孔径的合并时基于子孔径的成像算法需要一定的重叠率来保证合并后第二全孔径SAR回波数据的频谱的光滑度，进而使成像质量不受子孔径划分的影响。

[0089] 优选地，所述将距离向处理后的所述两段以上子孔径SAR回波数据合并为第二全孔径SAR回波数据，包括：
[0090] 获取相邻两段所述子孔径SAR回波数据的重叠率；
[0091] 根据所述重叠率将距离向处理后的所述两段以上子孔径SAR回波数据合并为第
二全孔径SAR回波数据。

[0092] 例如，极限情况下，相邻两个子孔径间需要的最大的重叠率为式（3）：
[0093]
[0094] 其中，β是SAR天线的波束宽度，F r是发射线性调频信号的带宽，θ是斜视角，f0是载波频率，L
满足式（4），
sub
[0095]
[0096] 由于子孔径时间带宽积的限制，划分子孔径的个数是有限的，式（3）得到的是极限重叠率，在实际的成像处理中，利用一个相对较小的重叠率就可以得到比较理想的成像质量。

[0097] 例如，采用两块英伟达（NVDIA）型号为Tesla C2050的GPU，型号为Intel Xeon E5620的CPU，且CPU和GPU的参数如表1所示：
[0098]
器件型号
CPU2Intel Xeon E56204Core2.4GHz
RAM96GB DDR31033MHz
GPU2Tesla C2050448Cores
GPU RAM3GB GDDR5
[0099] 表1
[0100] 若第一SAR回波数据的大小为65024x32768，计算得到相邻子孔径间的理论重叠点数（对应上述重叠率）为212，因此，子孔径间的重叠点数可设为256，共划分18个子孔径，每个子孔径的长度为4096。

[0101] 步骤103：将所述第二全孔径SAR回波数据沿距离向划分为两个以上SAR回波数据块，对所述两个以上SAR回波数据块进行方位向处理。

[0102] 优选地，所述对所述两个以上SAR回波数据块进行方位向处理之前，所述方法还包括：分别对所述两个以上SAR回波数据块进行转置。

[0103] 优选地，借助多线程技术，使得CPU在进行SAR回波数据块转置的同时，GPU进行SAR回波数据块方位向处理。

由于各个SAR回波数据块的处理并不相关，所以CPU完成第一块SAR回波数据块的转置后，在CPU转置第二SAR回波数据块的同时，GPU可以对完成转置的第一SAR回波数据块进行方位向处理。

依次类推，CPU和GPU协同处理完所有SAR回波数据块的转置和方位向处理。

[0104] 所述对所述两个以上SAR回波数据块进行方位向处理，包括：
[0105] 按接收先后顺序分别对所述两个以上SAR回波数据块进行方位向傅里叶变换；[0106] 对方位向傅里叶变换后的所述SAR回波数据块进行方位压缩和相位校正；[0107] 对方位压缩和相位校正后的所述SAR回波数据块进行方位向傅里叶逆变换。

[0108] 步骤104：将方位向处理后的所述两个以上SAR回波数据块合并为SAR图像。

[0109] 本发明实施例还记载了一种适用于上述SAR成像方法的装置，如图2所示，所述装。