CN104215962A - 基于单片FPGA的Chirp Scaling成像方法 - Google Patents
基于单片FPGA的Chirp Scaling成像方法 Download PDFInfo
- Publication number
- CN104215962A CN104215962A CN201410379642.7A CN201410379642A CN104215962A CN 104215962 A CN104215962 A CN 104215962A CN 201410379642 A CN201410379642 A CN 201410379642A CN 104215962 A CN104215962 A CN 104215962A
- Authority
- CN
- China
- Prior art keywords
- module
- data
- mrow
- doppler
- complex multiplication
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000013139 quantization Methods 0.000 claims abstract description 16
- 230000017105 transposition Effects 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 15
- 238000005452 bending Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 102100023006 Basic leucine zipper transcriptional factor ATF-like 2 Human genes 0.000 description 1
- 101000903615 Homo sapiens Basic leucine zipper transcriptional factor ATF-like 2 Proteins 0.000 description 1
- 208000004350 Strabismus Diseases 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- YOBAEOGBNPPUQV-UHFFFAOYSA-N iron;trihydrate Chemical compound O.O.O.[Fe].[Fe] YOBAEOGBNPPUQV-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于单片FPGA的Chirp Scaling成像方法,具体过程为:在一片FPGA上构建FFT运算模块、多普勒中心频率估计和拟合模块、多普勒调频斜率估计和拟合模块、参数解算模块、Chirp Scaling因子生成模块、距离向补偿因子生成模块、方位向补偿因子生成模块、复乘模块、数据转置模块、核心计算模块和图像量化模块;基于所构建的模块执行Chirp Scaling成像的过程。本发明对Chirp Scaling成像实现的架构进行了分析,最大限度的提高成像算法的并行度,能够充分发挥FPGA的优势,进一步提高成像的实时性。
Description
技术领域
本发明属于信号处理技术领域,具体涉及一种基于单片FPGA的ChirpScaling成像方法。
背景技术
随着SAR成像技术的发展,对星载SAR分辨率的要求越来越高,原始数据的数据率急剧增大。由于图像数据的相关性远高于原始回波数据,如果能在星上进行实时成像,通过图像压缩算法可以大大降低需要传输的数据量。
目前常用的SAR成像平台有DSP、FPGA+DSP、FPGA+ASIC+DSP,其中FPGA的并行处理能力较强,适于处理不太复杂的大批量数据,DSP灵活性强,适于处理比较复杂的运算,ASIC集成度高、比较稳定,适于处理相对独立的核心模块。这三种工具协作可以完成SAR成像任务,FPGA作为主处理器,DSP和ASIC作为从处理器。然而虽然DSP能够完成复杂的运算,但是它只能进行串行操作,运行速度慢,不利于SAR实时成像。并且多板卡的结构设计比较复杂,体积较大,而且多块处理板之间通过接插件链接,传输数据很慢而且容易产生误码,对散热和稳定性造成很大的挑战,不利于星上实时成像的实现。
西安电子科技大学提出的发明名称为“基于FPGA的SAR成像方法”(专利申请号201310101045.3,公开号103197317)的专利申请。这种成像方法主要解决现有系统结构复杂、数据处理速度慢的问题,采用的是SPECAN算法,简化了SAR成像系统结构,提高了处理速度,适用于弹载模式下的快速成像。但是星载SAR成像要求比较高的分辨率,一般会采用Chirp Scaling成像算法。这种算法的运算比较复杂,三个补偿因子的生成及其参数解算包含三角函数、开根号等超越运算,需要占用很大的硬件逻辑资源。目前CS成像算法的实现多是利用DSP来进行参数解算和补偿因子生成,FPGA完成FFT运算、数据转置以及数据与因子的复乘操作。但是DSP运行速度慢,板卡之间交互存在误码,会影响成像的实时性和稳定性。
发明内容
本发明的目的在于克服上述已有技术的不足,提出一种基于单片FPGA的Chirp Scaling成像方法,该方法能够在一片FPGA上实现Chirp Scaling成像,并且具有很好的实时性和稳定性。
实现本发明的技术方案如下:
一种基于单片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,
其中,λ表示回波波长,f表示方位向回波频率,V(r)表示载体的运动速度,表示波束前斜角,K表示回波的调频斜率,r表示雷达与目标之间的斜距,c表示光速。
进一步地,本发明在多普勒调频斜率估计完成后进行了剔除野值的处理,首先剔除几个最大值和最小值,然后剔除偏离均值较大的一些值,再利用最小二乘法对其进行拟合,从而得到精确的多普勒调频斜率值。
进一步地,本发明所述Chirp Scaling因子生成模块、距离向补偿因子生成模块和方位向补偿因子生成模块在每个聚焦深度内实现相应因子的生成。
进一步地,本发明所述量化模块存有针对不同场景进行图像量化时所需的最大值和最小值,选择与当前目标所处场景对应的最大值和最小值对二维时域的数据信号进行量化。
有益效果
第一,本发明使用单片FPGA实现Chirp Scaling成像,利用FPGA的高速并行处理特性,提高了SAR成像的实时性,而且避免了多板交互,减小了板间数据传输的误码率,提高了系统的稳定性。本发明对Chirp Scaling成像实现的架构进行了分析,最大限度的提高成像算法的并行度,能够充分发挥FPGA的优势,进一步提高成像的实时性。
第二,本发明利用低阶泰勒展开对现有计算弯曲因子Cs和等效距离向调频常数Ks进行改进,避免相差很大的两个数据相减,提高运算的精度,保证单片FPGA实现Chirp Scaling成像的精度,同时减少了运算量。
第三,本发明FFT运算模块均采用定点FFT运算,这样大大减小了硬件资源,从而保证整个Chirp Scaling成像可以在单片FPGA上实现。
第四,本发明中FFT运算模块、复乘模块、核心计算模块均采用时分复用的形式,从而减少所需的硬件资源,降低系统功耗。
第五,量化模块采用经验值存储的方法,能够大大减小硬件资源,提高成像实时性。
附图说明
图1是基于单片FPGA的Chirp Scaling成像方法流程图。
图2是基于单片FPGA的SAR成像算法模块结构图。
具体实施方式
下面结合附图和具体实例对本发明进行详细说明。
本发明一种基于单片FPGA的Chirp Scaling成像方法,如图1所示,具体步骤为:
步骤一:在一片FPGA上构建FFT运算模块、多普勒中心频率估计和拟合模块、多普勒调频斜率估计和拟合模块、参数解算模块、Chirp Scaling因子生成模块、距离向补偿因子生成模块、方位向补偿因子生成模块、复乘模块、数据转置模块、核心计算模块和图像量化模块,如图2所示。
步骤二:并行执行数据处理和Chirp Scaling因子生成:
数据处理为:FFT运算模块从回波数据中读取方位向回波数据,然后进行方位向定点FFT运算得到距离多普勒域的数据信号。
Chirp Scaling因子生成的过程为:
S101、多普勒中心频率估计和拟合模块接收回波数据,对回波数据进行自相关估计和最小二乘拟合,得到多普勒中心频率值;
S102、参数解算模块计算补偿因子生成所需的参数;
参数解算模块首先根据多普勒中心频率值计算回波方位向频率,其次根据事先输入的回波数据对应的初始参数以及回波方位向频率,通过加、减、乘、除、开根号等操作计算补偿因子生成所需要的参数;其中在执行占用硬件资源较大的(除法和开根号)运算时,所述参数解算模块调用核心计算模块中的运算资源实现;
补偿因子生成所需要的参数主要包括弯曲因子Cs和等效距离向调频常数Ks;
设 则有
其中,λ表示回波波长,f表示方位向回波频率,V(r)表示载体的运动速度,表示波束前斜角,K表示回波的调频斜率,r表示雷达与目标之间的斜距;
为了提高补偿因子计算的精度和运算实时性,本发明中对补偿因子的参数计算中用到的一些运算算法进行了改进。在计算弯曲因子Cs和等效距离向调频常数Ks时,跟据Chirp Scaling成像的特性可知,xpowr的值很小,远远小于1。一般硬件系统中都采用单精度浮点数据格式,当xpowr的值小于10-5时,(1-xpowr)得到的结果在单精度格式下表示为1,即有一大批接近于0的xpowr值被舍为0,则
本发明利用低阶泰勒展开来避免两个相差特别大的数据相减带来的误差。
针对弯曲因子Cs,利用一阶泰勒展开可得:
由于星载SAR的波束前斜角接近于90度,故是一个接近于1的值,所以是一个很小的数值,而也是一个很小的数值,两者数值相差不太大,在进行相加操作时损失的精度较小。与传统的计算方法比,本发明的计算方法可以修正的运算误差。
针对等效距离向调频常数Ks,利用一阶泰勒展开可以做以下简化:
则, 可简化为
同样可以修正计算误差,同时减小运算量。
S103、Chirp Scaling因子生成模块根据补偿因子生成所需要的参数,计算Chirp Scaling因子;其中生成Chirp Scaling因子需要用到的除法和三角函数运算调用核心计算模块中的运算资源。
传统成像方法是针对每条数据都生成一条相应的补偿因子与该条数据相乘,由于Chirp Scaling因子的运算比较复杂,其中除法、开根号和三角函数运算所需要的运算时间比较长,这样会导致Chirp Scaling因子生成的速度跟不上数据传输的速度,可能会造成数据流等待补偿因子计算的情况,从而会增加成像时间,降低成像的实时性。为了提高成像的实时性,本发明在每个聚焦深度内更新一次Chirp Scaling因子。即本发明是令相邻的几条数据与同一条补偿因子相乘,这样能保证补偿因子在数据到来之前就已经计算完成,则数据就可以进行流水处理,大大提高成像的实时性。
同时,本步骤中FFT运算模块对输入的方位向回波数据进行方位向定点FFT运算,完成方位向信号的时频转换,得到距离多普勒域的数据信号。一般的SAR成像方法中多采用单精度浮点的FFT运算,不仅运算时间长,而且占用的硬件资源很大。本发明采用定点FFT运算,在能够满足系统精度要求的前提下使用最小的数据字长,从而能够大大减小硬件资源。
步骤三:复乘模块对距离多普勒域的数据信号与Chirp Scaling因子进行复乘操作,完成第一次复乘操作;
步骤四:并行计算二维频域的信号数据和距离向补偿因子;
数据转置模块对第一次复乘操作得到的数据进行转置,然后利用FFT运算模块对转置后的数据进行距离向定点FFT运算,完成距离向信号的时频转换,得到二维频域的信号数据;
这里的FFT操作也是采用的定点FFT运算模块,而且与步骤五利用的是同一个运算模块,在不同的阶段进行时分复用,大大节约了硬件资源。
距离向补偿因子生成模块根据所述补偿因子生成所需要的参数,计算得到距离向补偿因子;其中生成距离向补偿因子需要用到的除法和三角函数运算调用核心计算模块中的运算资源;
这一步骤中距离向补偿因子的生成也是采用聚焦深度内更新的方法,令相邻的几条数据与同一条距离向补偿因子相乘,大大提高成像的实时性。
步骤五:复乘模块对二维频域的信号数据与距离向补偿因子进行复乘操作,完成第二次复乘操作。
步骤六:FFT运算模块对第二次复乘操作之后的数据进行距离向IFFT操作,完成距离向信号的时频转换,得到距离多普勒域的数据信号;
这里的FFT操作也是采用的定点FFT运算模块,而且与步骤五和步骤七利用的是同一个运算模块,在不同的阶段进行时分复用。
步骤七:数据转置模块对距离多普勒域的数据信号进行转置,并将转置后的数据传输给多普勒调频斜率估计和拟合模块和复乘模块,多普勒调频斜率估计和拟合模块针对转置后的数据信号,利用杂波锁定自聚焦法对多普勒调频斜率进行估计。
在Chirp Scaling成像算法中,多普勒调频斜率的值直接影响成像质量,如果多普勒调频斜率估计不准确,则会造成图像散焦。为了提高成像质量,本发明在多普勒调频斜率估计完成后进行了剔除野值的处理,首先剔除几个最大值和最小值,然后剔除偏离均值较大的一些值,再利用最小二乘法对其进行拟合,从而得到精确的多普勒调频斜率值。
步骤八:方位向补偿因子生成模块根据所述补偿因子生成所需要的参数和多普勒调频斜率,计算得到方位向补偿因子;其中生成方位向补偿因子需要用到的除法和三角函数运算调用核心计算模块中的运算资源。
这一步骤中方位向补偿因子的生成也是采用聚焦深度内更新的方法,令相邻的几条数据与同一条方位向补偿因子相乘,大大提高成像的实时性。
步骤九:复乘模块将得到的转置后的数据信号与方位向补偿因子进行复乘操作,完成第三次复乘操作。
步骤十:FFT运算模块对第三次复乘操作之后的数据进行方位向IFFT操作,完成方位向信号的时频转换,得到二维时域的数据信号;
这里的FFT操作也是采用的定点FFT运算模块,而且与之前步骤中利用的是同一个运算模块,在不同的阶段进行时分复用。
量化模块存有不同场景的图像量化时所需要的最大值和最小值,操作人员根据当前目标所处的场景选择相应的最大值和最小值,根据所选取的最大值和最小值对二维时域的数据信号进行量化,得到二维量化的SAR图像。
传统的SAR成像方法中,量化模块需要计算所有数据的均值和方差,需要的运算量很大,而且需要至少遍历两次数据,所需的运算时间很长。
本发明为了减小运算量,提高成像实时性,对量化模块进行了改进,根据大批量的不同场景下的数据仿真测试,得到各个不同场景下图像量化时所需要的最大值和最小值的经验值,存储在存储器中,当数据进行量化时,根据不同场景信息去选择不同的量化最大值和最小值,经过大量数据验证,这种方法的量化效果可以满足指标要求,而且能够大大减小硬件资源,提高成像实时性。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于单片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、步骤四及步骤八的计算过程中,在执行需占用硬件资源较大的运算时,调用核心计算模块来实现。
2.根据权利要求1所述基于单片FPGA的Chirp Scaling成像方法,其特征在于,所述参数包括弯曲因子Cs和等效距离向调频常数Ks,
其中,λ表示回波波长,f表示方位向回波频率,V(r)表示载体的运动速度,表示波束前斜角,K表示回波的调频斜率,r表示雷达与目标之间的斜距,c表示光速。
3.根据权利要求1所述基于单片FPGA的Chirp Scaling成像方法,其特征在于,在步骤七多普勒调频斜率估计完成后进行了剔除野值的处理,首先剔除最大值和最小值,然后剔除偏离均值较大的一些值,再利用最小二乘法对其进行拟合,得到精确的多普勒调频斜率。
4.根据权利要求1所述基于单片FPGA的Chirp Scaling成像方法,其特征在于,所述Chirp Scaling因子生成模块、距离向补偿因子生成模块和方位向补偿因子生成模块在每个聚焦深度内实现因子的生成。
5.根据权利要求1所述基于单片FPGA的Chirp Scaling成像方法,其特征在于,所述量化模块存有针对不同场景进行图像量化时所需的最大值和最小值,选择与当前目标所处场景对应的最大值和最小值,对二维时域的数据信号进行量化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410379642.7A CN104215962B (zh) | 2014-08-04 | 2014-08-04 | 基于单片FPGA的Chirp Scaling成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410379642.7A CN104215962B (zh) | 2014-08-04 | 2014-08-04 | 基于单片FPGA的Chirp Scaling成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104215962A true CN104215962A (zh) | 2014-12-17 |
CN104215962B CN104215962B (zh) | 2016-08-17 |
Family
ID=52097659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410379642.7A Active CN104215962B (zh) | 2014-08-04 | 2014-08-04 | 基于单片FPGA的Chirp Scaling成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104215962B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108614266A (zh) * | 2018-03-13 | 2018-10-02 | 南京航空航天大学 | 一种视频sar高速处理技术的fpga的实现方法 |
CN108627830A (zh) * | 2017-03-22 | 2018-10-09 | 北京行易道科技有限公司 | 车辆成像方法和系统 |
CN113885029A (zh) * | 2021-09-28 | 2022-01-04 | 理工雷科电子(西安)有限公司 | 一种基于DSP实现Sar成像中提高精度和时效性的方法 |
CN114626006A (zh) * | 2022-03-21 | 2022-06-14 | 电子科技大学 | 雷达成像中cs算法补偿因子实时生成的fpga实现方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7511656B2 (en) * | 2006-02-10 | 2009-03-31 | Raytheon Company | Spotlight synthetic aperture radar (SAR) system and method for generating a SAR map in real-time using a modified polar format algorithm (PFA) approach |
CN101441271A (zh) * | 2008-12-05 | 2009-05-27 | 航天恒星科技有限公司 | 基于gpu的sar实时成像处理设备 |
US7876257B2 (en) * | 2008-04-28 | 2011-01-25 | Mitsubishi Electric Research Laboratories, Inc. | Method and apparatus for compressing SAR signals |
CN102298139B (zh) * | 2011-05-18 | 2013-03-20 | 中国科学院计算技术研究所 | 基于fpga的sar成像系统的二维加窗方法 |
-
2014
- 2014-08-04 CN CN201410379642.7A patent/CN104215962B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7511656B2 (en) * | 2006-02-10 | 2009-03-31 | Raytheon Company | Spotlight synthetic aperture radar (SAR) system and method for generating a SAR map in real-time using a modified polar format algorithm (PFA) approach |
US7876257B2 (en) * | 2008-04-28 | 2011-01-25 | Mitsubishi Electric Research Laboratories, Inc. | Method and apparatus for compressing SAR signals |
CN101441271A (zh) * | 2008-12-05 | 2009-05-27 | 航天恒星科技有限公司 | 基于gpu的sar实时成像处理设备 |
CN102298139B (zh) * | 2011-05-18 | 2013-03-20 | 中国科学院计算技术研究所 | 基于fpga的sar成像系统的二维加窗方法 |
Non-Patent Citations (5)
Title |
---|
CHI-LI YU等: "TRANSPOSE-FREE SAR IMAGINE ON FPGA PLATFORM", 《CIRCUIT AND SYSTEMS(ISCAS),2012 IEEE INTERNATIONAL SYMPOSIUM》 * |
曾海彬等: "聚束式SAR斜视模型子孔径CS算法", 《北京理工大学学报》 * |
杨磊等: "高效机载SAR实时成像处理系统设计", 《电路与系统学报》 * |
熊君君等: "用FPGA实现星载SAR实时成像处理器的工程方法", 《计算机工程》 * |
闫雯等: "一种补偿因子区域不变的CS成像算法", 《北京理工大学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627830A (zh) * | 2017-03-22 | 2018-10-09 | 北京行易道科技有限公司 | 车辆成像方法和系统 |
CN108614266A (zh) * | 2018-03-13 | 2018-10-02 | 南京航空航天大学 | 一种视频sar高速处理技术的fpga的实现方法 |
CN113885029A (zh) * | 2021-09-28 | 2022-01-04 | 理工雷科电子(西安)有限公司 | 一种基于DSP实现Sar成像中提高精度和时效性的方法 |
CN113885029B (zh) * | 2021-09-28 | 2024-08-16 | 理工雷科电子(西安)有限公司 | 一种基于DSP实现Sar成像中提高精度和时效性的方法 |
CN114626006A (zh) * | 2022-03-21 | 2022-06-14 | 电子科技大学 | 雷达成像中cs算法补偿因子实时生成的fpga实现方法 |
CN114626006B (zh) * | 2022-03-21 | 2023-03-14 | 电子科技大学 | 雷达成像中cs算法补偿因子实时生成的fpga实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104215962B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109212526B (zh) | 用于高频地波雷达的分布式阵列目标角度测量方法 | |
CN104215962B (zh) | 基于单片FPGA的Chirp Scaling成像方法 | |
Patton et al. | Efficient design of radar waveforms for optimised detection in coloured noise | |
CN103197317B (zh) | 基于fpga的sar成像方法 | |
CN105842694B (zh) | 一种基于ffbp sar成像的自聚焦方法 | |
US9063220B2 (en) | Object information acquiring apparatus | |
CN111736131A (zh) | 一种剔除一比特信号谐波虚假目标的方法及相关组件 | |
US9022937B2 (en) | Ultrasound device and method for real-time computation of point-by-point apodization coefficients with trigonometric functions | |
CN109597072B (zh) | 一种双基合成孔径雷达sar系统的成像处理方法及装置 | |
Shi et al. | Accelerating parallel Jacobi method for matrix eigenvalue computation in DOA estimation algorithm | |
CN107329110A (zh) | 基于稀疏阵列直接内插的波达方向估计方法 | |
CN104188687A (zh) | 基于超声回波射频信号的多普勒血流速度估测方法和系统 | |
CN104007437A (zh) | 基于fpga和多dsp的sar实时成像处理方法 | |
CN109765526A (zh) | 一种基于空间谱的目标搜索方法及装置 | |
Tan et al. | An efficient range-Doppler domain ISAR imaging approach for rapidly spinning targets | |
US9664782B2 (en) | Object information acquiring apparatus | |
CN109541567B (zh) | 基于深度学习的高速机动目标检测方法 | |
Yang et al. | GF-3 data real-time processing method based on multi-satellite distributed data processing system | |
CN114626006B (zh) | 雷达成像中cs算法补偿因子实时生成的fpga实现方法 | |
CN116243313A (zh) | 一种基于距离分区的sar快速智能稀疏自聚焦技术 | |
CN113702971B (zh) | 用于合成孔径无源定位系统中的雷达波束宽度设计方法 | |
CN114706078A (zh) | 基于数字信号处理的太赫兹合成孔径雷达成像系统及方法 | |
CN115856782A (zh) | 一种基于fpga实现框架的车载雷达定位方法 | |
CN113203997A (zh) | 基于硬浮点fpga的雷达超分辨测向方法、系统及数据处理设备 | |
KR20210115194A (ko) | 레이더 영상 생성 방법 및 이를 수행하는 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CB03 | Change of inventor or designer information |
Inventor after: Chen He Inventor after: Yang Zhu Inventor after: Yan Wen Inventor after: Zeng Tao Inventor after: Long Teng Inventor before: Chen He Inventor before: Yan Wen Inventor before: Zeng Tao Inventor before: Long Teng |
|
CB03 | Change of inventor or designer information |