基于FPGA和多个多核DSP的SAR前侧视成像方法技术

本发明专利技术提出一种基于FPGA和多个多核DSP的SAR前侧视成像方法，用于解决现有高速飞行器SAR成像平台中存在的成像算法实时性差的技术问题，实现步骤为：1.初始化信号处理机及SAR参数；2.FPGA获取回波信号矩阵块并发送；3.DSPn对回波信号矩阵块

【技术实现步骤摘要】
基于FPGA和多个多核DSP的SAR前侧视成像方法
本专利技术属于数字信号处理
，涉及一种SAR前侧视成像方法，具体是一种基于FPGA和多个多核DSP的SAR前侧视成像方法，可应用于高速飞行器快速成像处理等领域。
技术介绍
合成孔径雷达(SyntheticApertureRadar，SAR)是雷达成像的主体，其应用范围最广。合成孔径雷达成像技术通过对雷达回波信号矩阵进行距离维和方位维两维处理，获得距离维和方位维两维高分辨率的SAR图像，将目标的形状和精细结构特征清晰地呈现出来，大大提高了对目标的检测和鉴别能力。在实际应用中，由于要保证高速飞行器具有一定的机动能力，进行例如转弯变向等动作，因此需要雷达能够对目标场景进行提前观测，这就要求SAR成像算法工作于前侧视模式下，相比工作于传统的正侧视模式，前者的运算量要大的多，成像的实时性难以保证，限制了SAR成像技术在高速飞行器中的应用。现有SAR前侧视成像算法可以分为两类：一类是基于傅里叶变换的频域SAR成像算法，另一类是基于逐像素点插值和相干积累的时域SAR成像算法。然而，大多数频域SAR成像算法对目标传输函数进行了近似，这些成像算法中的近似条件对SAR参数较为敏感，在大场景成像中带来较为严重的空变误差，导致SAR图像质量变差，从而限制了SAR成像技术在高速飞行器中的应用范围。为了解决上述问题，现有的研究现状是，一方面，使用频域SAR成像算法，并控制SAR成像的场景范围。该种方法降低了SAR成像算法的运算复杂度，保证了SAR成像的实时性。但是该种方法仅适用于成像场景范围不大的情况，在成像场景范围较大的情况下，该种成像方法运算复杂度增加，SAR成像的实时性降低；另一方面，采用时域SAR成像算法，并增加信号处理机上数字信号处理芯片的个数。时域成像算法观测几何灵活、聚焦精度高、图像分辨率可控，并且时域成像算法具有较高的并行度，但时域SAR成像算法的计算复杂度十分巨大，通过增加信号处理机上数字信号处理芯片的个数，可以有效解决时域成像算法因计算复杂度巨大而导致的实时性差的问题。赵浩浩在2019年3月发表的名为大斜视时域SAR成像算法的多核DSP并行架构实现技术一文中，公开了一种利用四个多核DSP，两两分组，并采用组内流水，组间乒乓的方式，实现对大斜视时域SAR成像算法的处理方法，该方法利用时域SAR成像算法解决了频域成像算法在大斜视角下成像质量较差的问题。但是，该方法的不足之处是：第一，该方法采用的时域SAR成像算法利用两维插值的方式实现子孔径图像融合，导致时域SAR成像算法计算复杂度增加，该方法成出一幅SAR前侧视图像的时间需要2058ms，时域SAR成像算法实时性差的问题依然存在。第二，该方法设计的多核DSP并行架构未能充分利用多核DSP芯片的信号处理能力，使用较多硬件资源，信号处理机功耗较大。
技术实现思路
本专利技术的目的在于针对上述技术中存在的不足，提出一种基于FPGA和多个多核DSP的SAR前侧视成像方法，用于解决现有高速飞行器SAR成像平台中存在的成像算法实时性差的技术问题。为实现上述目的，本专利技术采取的技术方案，包括如下步骤：(1)初始化信号处理机及SAR参数：初始化信号处理机上的FPGA获取SAR向目标场景发送载频为fc的线性调频信号所产生的回波采样矩阵为SM×N；初始化信号处理机上的g个并行排列的多核DSP芯片为{DSP1,DSP2,…,DSPn,…,DSPg}；初始化包括合成孔径长度L、子孔径长度lsub以及斜距Rs的SAR参数；其中，M为距离维采样点数，N为方位维采样点数，DSPn为第n个包含Q个核的DSP芯片，g为多核DSP芯片的个数，Q≥8，g≥3，n∈[1,g]；(2)FPGA获取回波信号矩阵块并发送：FPGA对回波采样矩阵SM×N进行距离维脉冲压缩，得到距离维脉冲压缩后回波信号矩阵并将依方位维划分为与多核DSP芯片个数相等的g个回波信号矩阵块后，将第n个回波信号矩阵块发送至对应的DSPn内存中，其中，gn为的列数，P为距离维脉冲压缩；(3)DSPn对回波信号矩阵块进行距离维FFT插值：DSPn将回波信号矩阵块依方位维划分为Q个回波信号矩阵块并通过第q核以J为插值倍数对进行距离维FFT插值，得到距离维FFT插值后回波信号矩阵块表示距离维FFT插值后第q+1个回波信号矩阵块，其中，gn(q+1)为的列数，K为进行距离维FFT插值后距离维点数，K＝M×J，I为FFT插值，q∈{0,1,…,Q-1}；(4)DSPn对距离维FFT插值后回波信号矩阵块进行BP积分：(4a)DSPn在以SAR的合成孔径中心为极点、SAR的航迹为极轴的统一极坐标系(r,θ)中，计算SAR到任意点目标U(rU,θU)的瞬时距离R(D；rU,θU)，其中，r和θ分别为极径和极角，rU为任意点目标U在(r,θ)中的极径，θU为任意点目标U在(r,θ)中的极角，D为SAR在(r,θ)中的极径，D∈[-L/2,L/2)；(4b)DSPn依方位维将FFT插值后回波信号矩阵块划分为个回波信号矩阵块并利用R(D；rU,θU)对进行BP积分，得到低分辨子孔径图像其中，为第x个回波信号矩阵块，gnx为的列数，Sx为第x幅低分辨子孔径图像，(5)DSPn对进行子孔径图像融合；(5a)DSPn根据线性调频信号载频fc、合成孔径长度L、子孔径长度lsub、斜距Rs，计算Sx的真实波数谱中心Kx，并通过Kx对Sx的波数谱中心进行校正，得到波数谱中心校正后的低分辨子孔径图像Sx′为波数谱中心校正后的第x幅低分辨子孔径图像；(5b)DSPn对波数谱中心校正后的低分辨子孔径图像在方位维进行IFFT，得到幅时域子孔径图像为波数谱中心校正后的第x幅时域子孔径图像，其中，T为时域；(5c)DSPn对幅时域子孔径图像在方位维进行叠加，得到全孔径波数谱SAR图像；(5d)DSPn对全孔径波数谱SAR图像在方位维进行FFT，得到SAR前侧视图像。本方明与现有技术相比，具有如下优点：1.本专利技术通过对BP积分后得到的低分辨子孔径图像进行波数谱中心校正，方位维IFFT处理，最终通过子孔径图像叠加的方法实现子孔径图像融合，与现有采用在距离维与方位维进行二维插值实现子孔径图像融合的技术相比，本专利技术的成像算法计算复杂度大幅降低，因此多核DSP芯片的负载与开销得以降低，时域SAR成像算法的实时性得到大幅提升，本专利技术成出一幅SAR前侧视图像的时间仅需要789ms。2.本专利技术对FFT插值后的回波数据进行BP积分时，采用逐列的对回波数据进行BP积分的方法，相比于现有的对回波数据方位维进行分块BP积分的方法，解决了因回波数据量大而BP积分必须在计算速度较低的DSP内存中运行的问题，因此使用较少的硬件资源，即可保证成像算法的实时性，并且信号处理机的功耗更低。附图说明图1为本专利技术采用的信号处理机结构图；图2为本专利技术的实现流程图；图3为本专利技术回波采样矩阵SM×N的DSP片间本文档来自技高网...

【技术保护点】
1.一种基于FPGA和多个多核DSP的SAR前侧视成像方法，其特征在于，包括如下步骤：/n(1)初始化信号处理机及SAR参数：/n初始化信号处理机上的FPGA获取SAR向目标场景发送载频为f

【技术特征摘要】
1.一种基于FPGA和多个多核DSP的SAR前侧视成像方法，其特征在于，包括如下步骤：
(1)初始化信号处理机及SAR参数：
初始化信号处理机上的FPGA获取SAR向目标场景发送载频为fc的线性调频信号所产生的回波采样矩阵为SM×N；初始化信号处理机上的g个并行排列的多核DSP芯片为{DSP1,DSP2,…,DSPn,…,DSPg}；初始化包括合成孔径长度L、子孔径长度lsub以及斜距Rs的SAR参数；其中，M为距离维采样点数，N为方位维采样点数，DSPn为第n个包含Q个核的DSP芯片，Q≥8，g≥3，n∈[1,g]；
(2)FPGA获取回波信号矩阵块并发送：
FPGA对回波采样矩阵SM×N进行距离维脉冲压缩，得到距离维脉冲压缩后回波信号矩阵并将依方位维划分为与多核DSP芯片个数相等的g个回波信号矩阵块后，将第n个回波信号矩阵块发送至对应的DSPn内存中，其中，gn为的列数，P为距离维脉冲压缩；
(3)DSPn对回波信号矩阵块进行距离维FFT插值：
DSPn将回波信号矩阵块依方位维划分为Q个回波信号矩阵块并通过第q核以J为插值倍数对进行距离维FFT插值，得到距离维FFT插值后回波信号矩阵块其中，表示距离维FFT插值后第q+1个回波信号矩阵块，gn(q+1)为的列数，K为进行距离维FFT插值后距离维点数，K＝M×J，I为FFT插值，q∈{0,1,…,Q-1}；
(4)DSPn对距离维距离维FFT插值后回波信号矩阵块进行BP积分：
(4a)DSPn在以SAR的合成孔径中心为极点、SAR的航迹为极轴的统一极坐标系(r,θ)中，计算SAR到任意点目标U(rU,θU)的瞬时距离R(D；rU,θU)，其中，r和θ分别为极径和极角，rU为任意点目标U在(r,θ)中的极径，θU为任意点目标U在(r,θ)中的极角，D为SAR在(r,θ)中的极径，D∈[-L/2,L/2)；
(4b)DSPn依方位维将FFT插值后回波信号矩阵块划分为个回波信号矩阵块并利用R(D；rU,θU)对进行BP积分，得到低分辨子孔径图像其中，为第x个回波信号矩阵块，gnx为的列数，Sx为第x幅低分辨子孔径图像，
(5)DSPn对进行子孔径图像融合；
(5a)DSPn根据线性调频信号载频fc、合成孔径长度L、子孔径长度lsub、斜距Rs，计算Sx的真实波数谱中心K...

【专利技术属性】
技术研发人员：谢荣，张政，刘峥，何志超，冉磊，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61