发明名称 基于单片FPGA的Chirp Scaling成像方法
摘要 本发明提供一种基于单片FPGA的Chirp Scaling成像方法,具体过程为:在一片FPGA上构建FFT运算模块、多普勒中心频率估计和拟合模块、多普勒调频斜率估计和拟合模块、参数解算模块、Chirp Scaling因子生成模块、距离向补偿因子生成模块、方位向补偿因子生成模块、复乘模块、数据转置模块、核心计算模块和图像量化模块;基于所构建的模块执行Chirp Scaling成像的过程。本发明对Chirp Scaling成像实现的架构进行了分析,最大限度的提高成像算法的并行度,能够充分发挥FPGA的优势,进一步提高成像的实时性。
申请公布号 CN104215962B 申请公布日期 2016.08.17
申请号 CN201410379642.7 申请日期 2014.08.04
申请人 北京理工大学 发明人 陈禾;闫雯;曾涛;龙腾
分类号 G01S13/90(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 仇蕾安;李爱英
主权项 一种基于单片FPGA的Chirp Scaling成像方法,其特征在于,具体步骤为:步骤一:在一片FPGA上构建FFT运算模块、多普勒中心频率估计和拟合模块、多普勒调频斜率估计和拟合模块、参数解算模块、Chirp Scaling因子生成模块、距离向补偿因子生成模块、方位向补偿因子生成模块、复乘模块、数据转置模块、核心计算模块和图像量化模块;步骤二:并行执行数据处理和Chirp Scaling因子生成:数据处理为:FFT运算模块从回波数据中读取方位向回波数据,然后进行方位向定点FFT运算得到距离多普勒域的数据信号;Chirp Scaling因子生成的过程为:S101、多普勒中心频率估计和拟合模块接收回波数据,对回波数据进行自相关估计和最小二乘拟合,得到多普勒中心频率值;S102、参数解算模块根据多普勒中心频率值计算补偿因子生成所需的参数;S103、Chirp Scaling因子生成模块根据所述参数,计算Chirp Scaling因子;步骤三、复乘模块对距离多普勒域的数据信号与Chirp Scaling因子进行复乘操作,完成第一次复乘操作;步骤四:并行计算二维频域的信号数据和距离向补偿因子;数据转置模块对第一次复乘操作得到的数据进行转置,FFT运算模块对转置后的数据进行距离向定点FFT运算,得到二维频域的信号数据;距离向补偿因子生成模块根据所述参数,计算得到距离向补偿因子;步骤五:复乘模块对二维频域的信号数据与距离向补偿因子进行复乘操作,完成第二次复乘操作;步骤六:FFT运算模块对第二次复乘操作之后的数据进行距离向IFFT操作,得到距离多普勒域的数据信号;步骤七:数据转置模块对距离多普勒域的数据信号进行转置,并将转置后的数据传输给所述多普勒调频斜率估计和拟合模块及所述复乘模块;多普勒调频斜率估计和拟合模块针对转置后的数据信号,利用杂波锁定自聚焦法估计多普勒调频斜率;步骤八:方位向补偿因子生成模块根据所述参数和多普勒调频斜率,计算得到方位向补偿因子;步骤九:复乘模块将得到的转置后的数据信号与方位向补偿因子进行复乘操作,完成第三次复乘操作;步骤十:FFT运算模块对第三次复乘操作之后的数据进行方位向IFFT操作,得到二维时域的数据信号;量化模块对二维时域的数据信号进行量化,得到二维量化的SAR图像;其中所述S102、S103、步骤四及步骤八的计算过程中,在执行需占用硬件资源较大的运算时,调用核心计算模块来实现;所述参数包括弯曲因子Cs和等效距离向调频常数Ks,<img file="FDA0000947725540000021.GIF" wi="718" he="117" /><maths num="0001"><math><![CDATA[<mrow><msub><mi>K</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>;</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>K</mi><mrow><mn>1</mn><mo>+</mo><mi>K</mi><mo>&CenterDot;</mo><mi>r</mi><mo>&CenterDot;</mo><mfrac><mrow><mn>2</mn><mi>&lambda;</mi></mrow><msup><mi>c</mi><mn>2</mn></msup></mfrac><mi>x</mi><mi>p</mi><mi>o</mi><mi>w</mi><mi>r</mi><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>x</mi><mi>p</mi><mi>o</mi><mi>w</mi><mi>r</mi><mo>)</mo></mrow><mn>3</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0000947725540000022.GIF" wi="875" he="178" /></maths>其中,<img file="FDA0000947725540000023.GIF" wi="390" he="147" />λ表示回波波长,f表示方位向回波频率,V(r)表示载体的运动速度,<img file="FDA0000947725540000024.GIF" wi="37" he="45" />表示波束前斜角,K表示回波的调频斜率,r表示雷达与目标之间的斜距,c表示光速。
地址 100081 北京市海淀区中关村南大街5号