CN107102328B - 基于fpga的实时成像信号处理方法及fpga - Google Patents
基于fpga的实时成像信号处理方法及fpga Download PDFInfo
- Publication number
- CN107102328B CN107102328B CN201710248893.5A CN201710248893A CN107102328B CN 107102328 B CN107102328 B CN 107102328B CN 201710248893 A CN201710248893 A CN 201710248893A CN 107102328 B CN107102328 B CN 107102328B
- Authority
- CN
- China
- Prior art keywords
- azimuth
- channel
- distance
- frequency
- fpga
- 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.)
- Active
Links
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
- G01S13/9004—SAR image acquisition 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
Abstract
一种基于FPGA的调频连续波方位向多通道FMCW SAR的实时成像信号处理方法及FPGA,包括:步骤一、计算重构矩阵;步骤二、重构方位向多通道数据,包括:步骤2.1、对各个通道的回波数据沿方位向分别间隔补零,并进行方位向傅里叶变换;步骤2.2、将方位向傅里叶变换之后各个通道方位向相同位置的点组合为一个向量并与重构矩阵相乘,得到重构完成的方位向数据;步骤2.3、对不同距离门的数据重复步骤2.3完成方位向重构,得到完成的重构后回波数据;步骤三、对重构后回波数据进行距离徙动校正及距离压缩处理;步骤四、对重构后回波数据进行方位向逆傅里叶变化,得到SAR图像。通过本发明可以提高图像处理质量。
Description
技术领域
本发明属于合成孔径雷达成像技术领域,尤其涉及一种基于FPGA的调频连续波方位多通道FMCW SAR的实时成像信号处理方法及FPGA。
背景技术
合成孔径雷达是一种主动微波遥感设备,能够全天时、全天候地实现对地观测,在国民经济与军事应用的诸多领域都发挥着重要的作用。
FMCW SAR由于其体积小、重量轻、分辨率高和成本低而在军用和民用领域备受青睐,世界各国纷纷开始了FMCW SAR系统的研制。多套机载FMCW SAR系统得到了成功研制并进行了试验验证。2009年,荷兰代尔夫特理工大学(TU Delft)、TNO(荷兰国家应用科学研究所)、SSBV公司等单位率先提出了星载FMCW SAR系统Panel SAR。
在星载条件下,更大的作用距离给系统提出了更高的要求。为了保证系统性能,在机载FMCW SAR的基础上,星载FMCW SAR需要对系统参数进行调整。提高发射功率、增大天线面积都是保证FMCW SAR系统信噪比的有效手段。但是提高发射功率对器件的要求大为增加,而且增加了系统的成本和功耗。增加天线面积可以提高天线的放大增益,但这样一方面会限制方位向的分辨率,另外一方面也会限制成像的测绘带宽度。同时,星载条件下,平台速度较快,需要较高的PRF以避免高方位向分辨率情况下的频谱混叠。方位多通道技术可以有效的解决上述问题,在避免方位频谱混叠的情况下,降低系统PRF并增加系统接收增益。但是方位多通道系统往往存在方位向非均匀采样的问题,需要进行方位多通道重构以保证成像质量。
为了提高星载SAR系统获取信息的时效性,需要在满足系统性能的前提下实现实时的成像处理。
发明内容
本发明解决的问题是现有方位向多通道系统存在方位向非均匀采样的问题;为解决所述问题本发明提供一种基于FPGA的方位多通道调频连续波SAR成像方法及FPGA。
本发明提供的基于FPGA的调频连续波方位向多通道FMCW SAR的实时成像信号处理方法包括:
步骤一、计算重构矩阵;
步骤二、重构方位向多通道数据,包括:
步骤2.1、对各个通道的回波数据沿方位向分别间隔补零,并进行方位向傅里叶变换;
步骤2.2、将方位向傅里叶变换之后各个通道方位向相同位置的点组合为一个向量并与重构矩阵相乘,得到重构完成的方位向数据;
步骤2.3、对不同距离门的数据重复步骤2.3完成方位向重构,得到完成的重构后回波数据;
步骤三、对重构后回波数据进行距离徙动校正及距离压缩处理;
步骤四、对重构后回波数据进行方位向逆傅里叶变化,得到SAR图像。
进一步,所述重构矩阵为:
其中:
f为方位向采样频率,di为第i+1通道与发射天线相位中心的距离,v为卫星等效速度,fp为脉冲重复频率,N为通道数。
进一步,所述步骤2.1包括:
间隔补零时,每个通道读入Na/N点,间隔补零至Na点,其中Na根据用户对每景图像大小的要求确定。每个通道数据补零方式示意如下:
通道1:a1,1 0…0|a1,N+1 0…0|a1,2N+1 0…0|…
通道2:0b2,2 0…0|0b2,N+2 0…0|0b2,2N+2 0…0|…
通道3:0 0c3,3 0…0|0 0c3,N+3 0…0|0 0c3,2N+3 0…0|…
…。
其中,第i个通道补零后的向量的第(k-1)*N+i个点为第i个通道补零前向量第k个点的值。
进一步,所述步骤2.1各个通道使用FPGA中的FFT核并行进行傅里叶变换。
进一步,所述步骤三依次包括:步骤3.1、对重构后回波数据与多普勒频移及频率变标函数相乘后做距离向傅里叶变换;步骤3.2、与残留视频相位函数相乘,并进行距离向傅里叶逆变换;;步骤3.3、与逆频率变标&二次距离压缩&距离徙动校正函数相乘,并进行距离向傅里叶变换;步骤3.4、与残留相位补偿函数&方位向匹配滤波相乘。
本发明还提供实现所述的实时成像信号处理方法的FPGA,包括:数据交互模块,方位向处理模块、距离向处理模块;所述数据交互模块接收方位向多通道回波数据后传输到方位向处理模块进行方位向处理,之后输出到距离向处理模块进行距离向处理。
本发明的优点包括:
提供一种基于FPGA的方位多通道调频连续波SAR成像方法,在利用多个方位接收通道突破方位向分辨率对雷达接收天线孔径的限制条件、保证系统高分辨并获得高天线增益的前提下,通过合理的算法设计结合FPGA实现实时的成像处理。
附图说明
图1为本发明实施例提供的基于FPGA的调频连续波方位向多通道FMCW SAR的实时成像信号处理方法的实现流程示意图。
图2为实现本发明实施例提供的基于FPGA的调频连续波方位向多通道FMCW SAR的实时成像信号处理方法的FPGA结构示意图。
图3为通过仿真得到的利用本发明提供的基于FPGA的调频连续波方位多通道FMCWSAR的实时成像信号处理方法得到的成像图像。
具体实施方式
下文中,结合附图和实施例对本发明做进一步阐述。
由背景技术可知,现有方位多通道FMCW SAR成像质量不够好,申请人分析认为是由于多方位通道系统存在方位向非均匀采样的问题,未解决所述问题,需要对方位向多通道进行重构,以解决所述问题,申请人经过进一步研究在本发明中提供基于FPGA的方位多通道调频连续波SAR成像方法,及实现所述方法的FPGA。
本发明提供的基于FPGA的调频连续波方位多通道FMCW SAR的实时成像信号处理方法包括:
步骤一、在开机后接收到回波数据之前完成重构矩阵及参考函数的计算;
本实施例中,在接收到回波数据之前完成重构矩阵及参考函数的计算,可以在后续数据处理中直接使用,提高成像的实时性,在本发明的其它实施例中,可以根据实际需要,在合适的步骤计算重构矩阵及参考函数。
步骤1.1、计算重构矩阵包括:
其中:
f为方位向采样频率,di为第i通道与发射天线相位中心的距离,v为卫星等效速度,fp为脉冲重复频率。在本发明的其它实施例中可以建立不同的重构矩阵,满足消除多方位通道系统方位向非均匀采样即可。
步骤1.2、参考函数的计算,如下:
1)多普勒频移及频率变标函数
其中:τ为距离向快时间,Kr为距离向调频率,
v为卫星等效速度;λ为波长;fη为方位向频率
2)残留视频相位函数
其中fτ为距离向频率。
3)逆频率变标&二次距离压缩&距离徙动校正函数
其中Rref为调频接收时的参考斜距;Rc为场景中心斜距,
4)残留相位补偿函数&方位向匹配滤波
步骤二、对各个通道的回波数据沿方位向分别间隔补零并使用不同的FFT核并行进行方位向傅里叶变换;
结合参考图2,所述间隔补零和方位向傅里叶变换在FPGA的方位向处理模块中进行,所述FFT核由FPGA生成。并行进行可以提高处理速度。
间隔补零时,每个通道读入Na/N点,间隔补零至Na点,每个通道数据补零方式示意如下:
通道1:a1,1 0…0|a1,M+1 0…0|a1,2M+1 0…0|…
通道2:0b2,2 0…0|0b2,M+2 0…0|0b2,2M+2 0…0|…
通道3:0 0c3,3 0…0|0 0c3,M+3 0…0|0 0c3,2M+3 0…0|…
…
通过算法公式推导结合FPGA的工作方式发现,间隔补零可以提高使用FPGA实现成像处理的效率。
步骤三、将方位向傅里叶变换之后每个通道方位向相同位置的点组合为一个向量并与重构矩阵相乘;之后将与重构矩阵相乘之后的序列依次排列得到重构完成的方位向数据;
[b1+(i-1)*M,b2+(i-1)*M,b3+(i-1)*M,…bM+(i-1)*M]=[a1,i,a2,i,a3,i,…aM,i]·G(f)
其中am,i分别为第m个通道方位向第i个点的值。
步骤四、对不同距离门的数据完成方位向重构得到完成的重构后回波数据;如图2所示,步骤二至步骤四在FPGA的方位向处理模块中进行。重构后回波数据输入到距离向处理模块,开始距离向处理。
步骤五、对所述多通道重构后回波数据进行距离徙动校正及距离压缩处理;包括:
步骤5.1、对重构后回波数据与多普勒频移及频率变标函数相乘后做距离向傅里叶变换;
步骤5.2、与残留视频相位函数相乘,并进行距离向傅里叶逆变换;
步骤5.3、与逆频率变标&二次距离压缩&距离徙动校正函数相乘,并进行距离向傅里叶变换;
步骤5.4、与残留相位补偿函数&方位向匹配滤波相乘。
步骤六、对完成距离向压缩的回波信号进行方位向逆傅里叶变化及聚焦处理,得到SAR图像。对回波信号进行逆傅里叶变化及聚焦处理的方法已为本领域技术人员所熟知,在此不再详述。
图3为通过仿真得到的利用本发明提供的基于FPGA的调频连续波方位多通道FMCWSAR的实时成像信号处理方法得到的成像图像。由图3可知通过本发明所提供的方案可以得到高质量的图像。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
Claims (3)
1.一种基于FPGA的调频连续波方位向多通道FMCW SAR的实时成像信号处理方法,其特征在于,包括:
步骤一、计算重构矩阵及参考函数;所述重构矩阵为:
其中:
f为方位向采样频率,di为第i+1通道与发射天线相位中心的距离,v为卫星等效速度,fp为脉冲重复频率,N为通道数;
所述参考函数包括:
1)多普勒频移及频率变标函数
其中:τ为距离向快时间,Kr为距离向调频率,
v为卫星等效速度;λ为波长;fη为方位向频率;
2)残留视频相位函数
fτ为距离向频率;
3)逆频率变标&二次距离压缩&距离徙动校正函数
其中Rref为调频接收时的参考斜距;Rc为场景中心斜距,
4)残留相位补偿函数&方位向匹配滤波
步骤二、重构方位向多通道数据,包括:
步骤2.1、对各个通道的回波数据沿方位向分别间隔补零,并进行方位向傅里叶变换;间隔补零时,每个通道读入Na/N点,间隔补零至Na点,其中Na根据用户对每景图像大小的要求确定;
步骤2.2、将方位向傅里叶变换之后各个通道方位向相同位置的点组合为一个向量并与重构矩阵相乘,得到重构完成的方位向数据;
步骤2.3、对不同距离门的数据重复步骤2.1~2.2完成方位向重构,得到完成的重构后回波数据;
步骤三、对重构后回波数据进行距离徙动校正及距离压缩处理;包括:步骤3.1、对重构后回波数据与多普勒频移及频率变标函数相乘后做距离向傅里叶变换;步骤3.2、与残留视频相位函数相乘,并进行距离向傅里叶逆变换;步骤3.3、与逆频率变标&二次距离压缩&距离徙动校正函数相乘,并进行距离向傅里叶变换;步骤3.4、与残留相位补偿函数&方位向匹配滤波相乘;
步骤四、对重构后回波数据进行方位向逆傅里叶变化,得到SAR图像。
2.依据权利要求1所述的基于FPGA的调频连续波方位向多通道FMCW SAR的实时成像信号处理方法,其特征在于,所述步骤2.1各个通道使用FPGA中的FFT核并行进行傅里叶变换。
3.实现权利要求1至2中任意一项所提供的实时成像信号处理方法的FPGA,其特征在于,包括:数据交互模块,方位向处理模块、距离向处理模块;所述数据交互模块接收方位向多通道回波数据后传输到方位向处理模块进行方位向处理,之后输出到距离向处理模块进行距离向处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710248893.5A CN107102328B (zh) | 2017-04-17 | 2017-04-17 | 基于fpga的实时成像信号处理方法及fpga |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710248893.5A CN107102328B (zh) | 2017-04-17 | 2017-04-17 | 基于fpga的实时成像信号处理方法及fpga |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107102328A CN107102328A (zh) | 2017-08-29 |
CN107102328B true CN107102328B (zh) | 2019-12-17 |
Family
ID=59675070
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710248893.5A Active CN107102328B (zh) | 2017-04-17 | 2017-04-17 | 基于fpga的实时成像信号处理方法及fpga |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107102328B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108983238A (zh) * | 2018-07-27 | 2018-12-11 | 山东航天电子技术研究所 | 一种星载Ka波段SAR实时成像系统 |
CN110161468B (zh) * | 2019-04-11 | 2023-05-26 | 上海卫星工程研究所 | 星载sar传输通道数据转成像接收通道功能模块设计方法 |
CN112505698B (zh) * | 2020-11-20 | 2023-09-12 | 内蒙古工业大学 | 多通道滑聚sar方位信号预处理方法、装置及存储介质 |
CN112135062B (zh) * | 2020-11-24 | 2021-02-05 | 宁波兰茜生物科技有限公司 | 一种基于fpga的相机聚焦控制的方法和装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323842A (zh) * | 2012-09-03 | 2013-09-25 | 中国科学院电子学研究所 | 一种调频连续波合成孔径雷达中的成像方法及装置 |
CN103969644B (zh) * | 2013-01-31 | 2016-06-01 | 中国人民解放军国防科学技术大学 | 一种多通道调频连续波sar成像方法 |
CN104833972B (zh) * | 2015-05-08 | 2017-09-19 | 电子科技大学 | 一种双基地调频连续波合成孔径雷达频率变标成像方法 |
CN105954751A (zh) * | 2016-06-03 | 2016-09-21 | 王辉 | Ka FMCW SAR的运动目标检测方法及装置 |
CN106443671A (zh) * | 2016-08-30 | 2017-02-22 | 西安电子科技大学 | 基于调频连续波的sar雷达动目标检测与成像方法 |
-
2017
- 2017-04-17 CN CN201710248893.5A patent/CN107102328B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107102328A (zh) | 2017-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111443339B (zh) | 一种双基sar空变校正成像方法、装置、设备及存储介质 | |
CN107102328B (zh) | 基于fpga的实时成像信号处理方法及fpga | |
CN105785327B (zh) | 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法 | |
CN103454632B (zh) | 一站固定式调频连续波双基地sar成像方法 | |
CN104698459B (zh) | 一种应用于缺失数据的条带sar压缩感知成像方法 | |
CN103323818B (zh) | 多通道合成孔径雷达系统非均匀采样奇异点的方法和装置 | |
CN104898107B (zh) | 一种多发多收合成孔径激光雷达信号处理方法 | |
CN108037497B (zh) | 多子阵合成孔径声纳数据的收发合置转换方法 | |
CN103336280B (zh) | 临近空间慢速平台合成孔径雷达大场景成像方法 | |
CN103472449B (zh) | 一种基于mimo图像域的bp宽带合成方法 | |
CN110187347B (zh) | 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法 | |
CN107544068A (zh) | 一种基于频域bp的图像域宽带合成方法 | |
CN103969644A (zh) | 一种多通道调频连续波sar成像方法 | |
CN103105610B (zh) | 基于非均匀采样的dpc-mab sar成像方法 | |
CN103885052A (zh) | 一种宽幅全极化星载sar的极化回波分离方法 | |
CN108008387B (zh) | 一种机载阵列天线下视三维成像方法 | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
CN103576151B (zh) | 基于压缩感知的方位向多通道sar成像方法及系统 | |
CN103576153A (zh) | 一种方位向多波束合成孔径雷达及其实现方法和装置 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN108845318B (zh) | 基于Relax算法的星载高分宽幅成像方法 | |
CN102928839A (zh) | 多通道波束指向sar的全孔径成像方法 | |
CN113189547B (zh) | 一种基于合成带宽频率变标的sar成像的方法及系统 | |
CN103217677A (zh) | 一种基于联合检测量的单通道sar动目标检测方法 | |
CN110308445A (zh) | 一种基于车载数字阵调频连续波雷达的成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |