CN106772372B - 一种Ka波段机载SAR***的实时成像方法和*** - Google Patents
一种Ka波段机载SAR***的实时成像方法和*** Download PDFInfo
- Publication number
- CN106772372B CN106772372B CN201611068315.5A CN201611068315A CN106772372B CN 106772372 B CN106772372 B CN 106772372B CN 201611068315 A CN201611068315 A CN 201611068315A CN 106772372 B CN106772372 B CN 106772372B
- Authority
- CN
- China
- Prior art keywords
- compensation
- data
- sar system
- doppler
- real time
- 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
- G01S13/9019—Auto-focussing of the SAR signals
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达技术领域,涉及一种Ka波段机载SAR***的实时成像方法和***。所述方法包括如下步骤:S1,对机载SAR***的原始成像数据进行预处理;S2,进行多普勒中心频率补偿;S3,进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理;S4,根据S3处理后的原始成像数据计算多普勒调频率,根据所述多普勒调频率进行相位补偿;S5,进行方位脉冲压缩处理,及横向位移补偿处理,得到输出数据;S6,根据所述输出数据得到单视实图像,对单视实图像进行量化、多视处理及拼接得到实时成像数据。本发明实现Ka波段机载SAR***的实时成像技术,满足ka波段机载SAR***高分辨成像需求,实时成像效果好,聚焦效果好。
Description
技术领域
本发明属于雷达技术领域,尤其涉及一种Ka波段机载SAR***的实时成像方法和***。
背景技术
高分辨率机载SAR***要求载机匀速直线飞行,但是由于受气流和载机自身等因素影响,其飞行轨迹常常偏离设定航迹,因此会引入运动误差。运动误差的大小对机载SAR实时成像的效果影响不同,轻则造成图像的信噪比下降、分辨率降低,重则造成图像散焦无法使用。不同频段SAR***对运动误差的敏感程度也不同,频段越高的雷达对运动误差越敏感,因此ka波段机载SAR***对运动误差的敏感程度相对于常见的X波段机载SAR***要高很多。
机载SAR***的实时成像通常将机载惯导数据和基于回波数据的多普勒参数估计相结合的方法。由于以往我国机载SAR***使用的惯导***精度不高,所以一般只用惯导数据进行粗补偿,然后再利用基于回波数据的多普勒参数估计进行精补偿,以达到高分辨率成像的要求。近几年虽然国外的一些高精度惯导设备的引进和国产惯导精度的提高,但是由于ka波段机载SAR***的特点,Ka波段机载SAR***的实时成像还未实现。
发明内容
本发明所要解决的技术问题是针对现有技术的不足,提供一种Ka波段机载SAR***的实时成像方法和***。
本发明解决上述技术问题的技术方案如下:一种Ka波段机载SAR***的实时成像方法,包括如下步骤:
S1,对机载SAR***的原始成像数据进行预处理;
S2,根据机载SAR***得到的惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿;
S3,对经多普勒中心频率补偿后的原始成像数据依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理;
S4,根据S3处理后的原始成像数据计算多普勒调频率,依据所述多普勒调频率对经距离徙动校正处理的原始成像数据进行相位补偿;
S5,对经相位补偿的原始成像数据进行方位脉冲压缩处理,及横向位移补偿处理,得到输出数据;
S6,根据所述输出数据得到单视实图像,对单视实图像进行量化和多视处理得到每个单幅实图像,对相邻单幅实图像进行拼接得到实时成像数据。
本发明的有益效果是:对机载SAR***的原始成像数据进行预处理,根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,再依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理,实现在实时成像时间内相邻图像中心位置和相位一致,依据多普勒调频率进行相位补偿,实时微调,方位脉冲压缩处理后进行横向位移补偿,得到输出数据,再经过量化和多视处理得到单幅实图像,拼接得到实时成像数据,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述S2包括:根据机载SAR***测量得到的惯导数据计算多普勒中心频率fdc,fdc的表达式是:
其中,λ表示雷达工作波长,R表示雷达瞬时斜距,VN表示载机北向速度,VE表示载机东向速度,VD表示载机地向速度,H表示载机相对目标场景的高度,Θ表示载机俯仰角,α表示载机航迹角,所述λ是机载SAR***参数,R、VN、VE、VD、H、Θ和α均是所述惯导数据;;
根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿。
采用上述进一步方案的有益效果是:根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
进一步,所述S1包括:对机载SAR***的原始成像数据进行距离向的脉冲压缩处理,所述原始成像数据包括回波信号数据。
进一步,所述运动误差补偿处理包括:根据惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿。
采用上述进一步方案的有益效果是:通过惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿实现在实时成像时间内相邻图像中心位置和相位一致,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
进一步,所述S5包括:对经相位补偿的原始成像数据进行方位脉冲压缩处理,再进行横向位移补偿处理,得到输出数据,所述横向位移Δx的表达式是:
其中,Δx表示横向位移,fdc表示所述多普勒中心频率,v表示载机飞行速度,R表示雷达瞬时斜距,θ表示载机俯仰角,所述λ是机载SAR***参数,v、R和θ均是所述惯导数据;
所述输出数据Xa的表达式是:
xa=ifft(S(f)×H(f)×exp(j2πf×Δx/v))
其中,所述S(f)表示回波信号频谱,H(f)表示滤波器频谱,f表示多普勒频率,所述S(f)由回波信号数据计算得出,H(f)由机载SAR***参数和惯导数据计算得出,所述多普勒频率是将时域信号转化为频率域信号的变量。
采用上述进一步方案的有益效果是:通过方位脉冲压缩处理后再进行横向位移补偿,得到输出数据,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
进一步,所述S4中所述多普勒调频率是所述多普勒中心频率的变化量。
本发明解决上述技术问题的另一个技术方案如下:一种Ka波段机载SAR***的实时成像***,包括:
预处理模块,用于对机载SAR***的原始成像数据进行预处理;
多普勒中心频率补偿模块,用于根据机载SAR***得到的惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿;
第一处理模块,用于对经多普勒中心频率补偿后的原始成像数据依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理;
相位补偿模块,根据第一处理模块处理后的原始成像数据计算多普勒调频率,依据所述多普勒调频率对经距离徙动校正处理的原始成像数据进行相位补偿;
方位脉冲压缩及位移补偿处理模块,用于对经相位补偿的原始成像数据进行方位脉冲压缩处理,及横向位移补偿处理,得到输出数据;
第二处理模块,用于根据所述输出数据得到单视实图像,对单视实图像进行量化和多视处理得到每个单幅实图像,对相邻单幅实图像进行拼接得到实时成像数据。
本发明的有益效果是:预处理模块对机载SAR***的原始成像数据进行预处理,多普勒中心频率补偿模块根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,第一处理模块依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理,实现在实时成像时间内相邻图像中心位置和相位一致,相位补偿模块依据多普勒调频率进行相位补偿,实时微调,方位脉冲压缩及位移补偿处理模块进行方位脉冲压缩处理后进行横向位移补偿,得到输出数据,第二处理模块进行量化和多视处理得到单幅实图像,拼接得到实时成像数据,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述多普勒中心频率补偿模块包括:用于根据机载SAR***测量得到的惯导数据计算多普勒中心频率fdc,fdc的表达式是:
其中,λ表示雷达工作波长,R表示雷达瞬时斜距,VN表示载机北向速度,VE表示载机东向速度,VD表示载机地向速度,H表示载机相对目标场景的高度,Θ表示载机俯仰角,α表示载机航迹角,所述λ是机载SAR***参数,R、VN、VE、VD、H、Θ和α均是所述惯导数据;
根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿。
采用上述进一步方案的有益效果是:通过多普勒中心频率补偿模块根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
进一步,所述预处理模块包括:对机载SAR***的原始成像数据进行距离向的脉冲压缩处理,所述原始成像数据包括回波信号数据。
进一步,所述运动误差补偿包括:根据惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿。
采用上述进一步方案的有益效果是:通过惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿实现在实时成像时间内相邻图像中心位置和相位一致,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
进一步,所述方位脉冲压缩及位移补偿处理模块包括:用于对经相位补偿的原始成像数据进行方位脉冲压缩处理,再进行横向位移补偿处理,得到输出数据,所述横向位移Δx的表达式是:
其中,Δx表示横向位移,fdc表示所述多普勒中心频率,v表示载机飞行速度,R表示雷达瞬时斜距,θ表示载机俯仰角,所述λ是机载SAR***参数,v、R和θ均是所述惯导数据;
所述输出数据Xa的表达式是:
xa=ifft(S(f)×H(f)×exp(j2πf×Δx/v))
其中,所述S(f)表示回波信号频谱,G(f)表示滤波器频谱,f表示多普勒频率,所述S(f)由回波信号数据计算得出,H(f)由机载SAR***参数和惯导数据计算得出,所述多普勒频率是将时域信号转化为频率域信号的变量。
采用上述进一步方案的有益效果是:通过方位脉冲压缩处理后再进行横向位移补偿,得到输出数据,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
附图说明
图1为现有技术中载机的实际航线和理想航线示意图;
图2为本发明实施例提供的一种Ka波段机载SAR***的实时成像方法流程示意图;
图3为本发明实施例提供的一种Ka波段机载SAR***的实时成像***结构示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
由于载机受气流和载机自身因素等影响,其飞行轨迹常常偏离设定轨迹,造成运动误差,由于运动误差使得载机的实时成像效果差。
如图1所示,实际飞行轨迹偏离设定轨迹,造成运动误差ΔR,其表达式为:
ΔrR=Y(tm)sinβ+(H-Z(tm))cosβ
其中,β表示载机上雷达天线侧视角,v表示载机飞行速度,H表示载机相对目标场景的高度,Rs表示载机设定轨迹相对于地面成像条带中心的距离,Xn表示目标在X轴上的坐标,tm表示时间,X(tm)、Y(tm)和Z(tm)表示载机上雷达天线相位中心在tm时刻的坐标,ΔrR表示视线误差,ΔrN表示航线误差,R0和R1分别是载机上雷达天线相位中心参考位置和实际位置相对于散射点A的斜距。
本发明实施例通过消除视线误差解决运动误差的问题,进而实现实时成像效果的改善。
本发明实施例提供一种Ka波段机载SAR***的实时成像方法和***,如图2所示,Ka波段机载SAR***的实时成像方法,包括如下步骤:
S1,对机载SAR***的原始成像数据进行预处理;
S2,根据机载SAR***得到的惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿;
S3,对经多普勒中心频率补偿后的原始成像数据依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理;
S4,根据S3处理后的原始成像数据计算多普勒调频率,依据所述多普勒调频率对经距离徙动校正处理的原始成像数据进行相位补偿;
S5,对经相位补偿的原始成像数据进行方位脉冲压缩处理,及横向位移补偿处理,得到输出数据;
S6,根据所述输出数据得到单视实图像,对单视实图像进行量化和多视处理得到每个单幅实图像,对相邻单幅实图像进行拼接得到实时成像数据。
对机载SAR***的原始成像数据进行预处理,根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,再依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理,实现在实时成像时间内相邻图像中心位置和相位一致,依据多普勒调频率进行相位补偿,实时微调,方位脉冲压缩处理后进行横向位移补偿,得到输出数据,再经过量化和多视处理得到单幅实图像,拼接得到实时成像数据,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述S2包括:根据机载SAR***测量得到的惯导数据计算多普勒中心频率fdc,fdc的表达式是:
其中,λ表示雷达工作波长,R表示雷达瞬时斜距,VN表示载机北向速度,VE表示载机东向速度,VD表示载机地向速度,H表示载机相对目标场景的高度,Θ表示载机俯仰角,α表示载机航迹角,所述λ是机载SAR***参数,,R、VN、VE、VD、H、Θ和α均是所述惯导数据;
根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿。
采用现有技术中的基于回波数据的多普勒中心频率fdc估计方法,对地面特殊场景采用现有技术中的估计方法和本发明实施例提供的方法分别计算fdc,该地面特殊场景采用基于回波数据的多普勒中心频率fdc估计方法可以得到准确的fdc,本发明实施例提供的fdc的计算方法与现有技术的fdc值一致。
根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述S1包括:对机载SAR***的原始成像数据进行距离向的脉冲压缩处理,所述原始成像数据包括回波信号数据。所述原始成像数据由机载雷达得到。
具体地,该实施例中,所述运动误差补偿处理包括:根据惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿。
通过惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿实现在实时成像时间内相邻图像中心位置和相位一致,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述S5包括:对经相位补偿的原始成像数据进行方位脉冲压缩处理,再进行横向位移补偿处理,得到输出数据,所述横向位移Δx的表达式是:
其中,Δx表示横向位移,fdc表示所述多普勒中心频率,v表示载机飞行速度,R表示雷达瞬时斜距,θ表示载机俯仰角,所述λ是机载SAR***参数,v、R和θ均是所述惯导数据;
方位向位移即横向位移,方位向信号可以等效为线性调频信号,若线性调频信号调制一个多普勒中心频率,对其进行脉冲压缩会产生横向位移,横向位移的表达式是:
其中:fdc是多普勒中心频率,μ是多普勒调频率,v是载机飞行速度。其多普勒调频率的表达式是:
其中:v是载机飞行速度,θ表示载机俯仰角,R表示雷达瞬时斜距,所述λ是机载SAR***参数。
将多普勒调频率的表达式带入横向位移的表达式中可得:
所述输出数据Xa的表达式是:
xa=ifft(S(f)×H(f)×exp(j2πf×Δx/v))
其中,所述S(f)表示回波信号频谱,H(f)表示滤波器频谱,f表示多普勒频率,所述S(f)由回波信号数据计算得出,H(f)由机载SAR***参数和惯导数据计算得出,所述多普勒频率是将时域信号转化为频率域信号的变量,所述时域信号是回波信号等。
为了消除由于多普勒中心频率造成的横向位移,可以采用在多普勒中心频率域,给信号频谱S(f)和滤波器频谱H(f)乘积再乘以一个线性相位项然后再经过逆傅立叶变换ifft()将信号转换成时域信号即可消除。
通过方位脉冲压缩处理后再进行横向位移补偿,得到输出数据,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述S4中所述多普勒调频率是所述多普勒中心频率的变化量。
本发明实施例提供一种Ka波段机载SAR***的实时成像***,如图3所示,Ka波段机载SAR***的实时成像***,包括:
预处理模块,用于对机载SAR***的原始成像数据进行预处理;
多普勒中心频率补偿模块,用于根据机载SAR***得到的惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿;
第一处理模块,用于对经多普勒中心频率补偿后的原始成像数据依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理;
相位补偿模块,用于根据第一处理模块处理后的原始成像数据计算多普勒调频率,依据所述多普勒调频率对经距离徙动校正处理的原始成像数据进行相位补偿;
方位脉冲压缩及位移补偿处理模块,用于对经相位补偿的原始成像数据进行方位脉冲压缩处理,及横向位移补偿处理,得到输出数据;
第二处理模块,用于根据所述输出数据得到单视实图像,对单视实图像进行量化和多视处理得到每个单幅实图像,对相邻单幅实图像进行拼接得到实时成像数据。
预处理模块对机载SAR***的原始成像数据进行预处理,多普勒中心频率补偿模块根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,第一处理模块依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理,实现在实时成像时间内相邻图像中心位置和相位一致,相位补偿模块依据多普勒调频率进行相位补偿,实时微调,方位脉冲压缩及位移补偿处理模块进行方位脉冲压缩处理后进行横向位移补偿,得到输出数据,第二处理模块进行量化和多视处理得到单幅实图像,拼接得到实时成像,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述多普勒中心频率补偿模块包括:用于根据机载SAR***测量得到的惯导数据计算多普勒中心频率fdc,fdc的表达式是:
其中,λ表示雷达工作波长,R表示雷达瞬时斜距,VN表示载机北向速度,VE表示载机东向速度,VD表示载机地向速度,H表示载机相对目标场景的高度,Θ表示载机俯仰角,α表示载机航迹角,所述λ是机载SAR***参数,R、VN、VE、VD、H、Θ和α均是所述惯导数据;
根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿。
通过多普勒中心频率补偿模块根据惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿,实现实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述预处理模块包括:对机载SAR***的原始成像数据进行距离向的脉冲压缩处理,所述原始成像数据包括回波信号数据。
具体地,该实施例中,所述运动误差补偿包括:根据惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿。
通过惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿实现在实时成像时间内相邻图像中心位置和相位一致,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述方位脉冲压缩及位移补偿处理模块包括:用于对经相位补偿的原始成像数据进行方位脉冲压缩处理,再进行横向位移补偿处理,得到输出数据,所述横向位移Δx的表达式是:
其中,Δx表示横向位移,fdc表示所述多普勒中心频率,v表示载机飞行速度,R表示雷达瞬时斜距,θ表示载机俯仰角,所述λ是机载SAR***参数,v、R和θ均是所述惯导数据;
所述输出数据Xa的表达式是:
xa=ifft(S(f)×H(f)×exp(j2πf×Δx/v))
其中,所述S(f)表示回波信号频谱,H(f)表示滤波器频谱,f表示多普勒频率,所述s(f)由回波信号数据计算得出,H(f)由机载SAR***参数和惯导数据计算得出,所述多普勒频率是将时域信号转化为频率域信号的变量,所述时域信号是回波信号等。
通过方位脉冲压缩处理后再进行横向位移补偿,得到输出数据,实时成像效果好,聚焦效果好,满足ka波段机载SAR***高分辨成像需求。
具体地,该实施例中,所述多普勒调频率是所述多普勒频率的变化量。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种Ka波段机载SAR***的实时成像方法,其特征在于,包括如下步骤:
S1,对机载SAR***的原始成像数据进行预处理;
S2,根据机载SAR***得到的惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿;
S3,对经多普勒中心频率补偿后的原始成像数据依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理;
S4,根据S3处理后的原始成像数据计算多普勒调频率,依据所述多普勒调频率对经距离徙动校正处理的原始成像数据进行相位补偿;
S5,对经相位补偿的原始成像数据进行方位脉冲压缩处理,及横向位移补偿处理,得到输出数据;
S6,根据所述输出数据得到单视实图像,对单视实图像进行量化和多视处理得到每个单幅实图像,对相邻单幅实图像进行拼接得到实时成像数据。
2.根据权利要求1所述一种Ka波段机载SAR***的实时成像方法,其特征在于,所述S2包括:根据机载SAR***测量得到的惯导数据计算多普勒中心频率fdc,fdc的表达式是:
其中,λ表示雷达工作波长,R表示雷达瞬时斜距,VN表示载机北向速度,VE表示载机东向速度,VD表示载机地向速度,H表示载机相对目标场景的高度,Θ表示载机俯仰角,α表示载机航迹角,所述λ是机载SAR***参数,R、VN、VE、VD、H、Θ和α均是所述惯导数据;
根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿。
3.根据权利要求2所述一种Ka波段机载SAR***的实时成像方法,其特征在于,所述S1包括:对机载SAR***的原始成像数据进行距离向的脉冲压缩处理,所述原始成像数据包括回波信号数据。
4.根据权利要求3所述一种Ka波段机载SAR***的实时成像方法,其特征在于,所述运动误差补偿处理包括:根据惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿。
5.根据权利要求4所述一种Ka波段机载SAR***的实时成像方法,其特征在于,所述S5包括:对经相位补偿的原始成像数据进行方位脉冲压缩处理,再进行横向位移补偿处理,得到输出数据,所述横向位移Δx的表达式是:
其中,Δx表示横向位移,fdc表示所述多普勒中心频率,ν表示载机飞行速度,R表示雷达瞬时斜距,Θ表示载机俯仰角,所述λ是机载SAR***参数,ν、R和Θ均是所述惯导数据;
所述输出数据xa的表达式是:
xa=ifft(S(f)×H(f)×exp(j2πf×Δx/ν))
其中,所述S(f)表示回波信号频谱,H(f)表示滤波器频谱,f表示多普勒频率,所述S(f)由回波信号数据计算得出,H(f)由机载SAR***参数和惯导数据计算得出,所述多普勒频率是将时域信号转化为频率域信号的变量。
6.根据权利要求1-5任意一项所述一种Ka波段机载SAR***的实时成像方法,其特征在于,所述S4中所述多普勒调频率是所述多普勒频率的变化量。
7.一种Ka波段机载SAR***的实时成像***,其特征在于,包括:
预处理模块,用于对机载SAR***的原始成像数据进行预处理;
多普勒中心频率补偿模块,用于根据机载SAR***得到的惯导数据计算多普勒中心频率,根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿;
第一处理模块,用于对经多普勒中心频率补偿后的原始成像数据依次进行方位预滤波、抽取、运动误差补偿和距离徙动校正处理;
相位补偿模块,用于根据第一处理模块处理后的原始成像数据计算多普勒调频率,依据所述多普勒调频率对经距离徙动校正处理的原始成像数据进行相位补偿;
方位脉冲压缩及位移补偿处理模块,用于对经相位补偿的原始成像数据进行方位脉冲压缩处理,及横向位移补偿处理,得到输出数据;
第二处理模块,用于根据所述输出数据得到单视实图像,对单视实图像进行量化和多视处理得到每个单幅实图像,对相邻单幅实图像进行拼接得到实时成像数据。
8.根据权利要求7所述的一种Ka波段机载SAR***的实时成像***,其特征在于,所述多普勒中心频率补偿模块包括:用于根据机载SAR***测量得到的惯导数据计算多普勒中心频率fdc,fdc的表达式是:
其中,λ表示雷达工作波长,R表示雷达瞬时斜距,VN表示载机北向速度,VE表示载机东向速度,VD表示载机地向速度,H表示载机相对目标场景的高度,Θ表示载机俯仰角,α表示载机航迹角,所述λ是机载SAR***参数,R、VN、VE、VD、H、Θ和α均是所述惯导数据;
根据多普勒中心频率对经预处理的原始成像数据进行多普勒中心频率补偿。
9.根据权利要求8所述的一种Ka波段机载SAR***的实时成像***,其特征在于,所述运动误差补偿包括:根据惯导数据解算出的一次运动误差和二次运动误差对经距离向的脉冲压缩、多普勒中心频率补偿、方位预滤波和抽取处理后的原始成像数据进行运动误差补偿。
10.根据权利要求9所述的一种Ka波段机载SAR***的实时成像***,其特征在于,所述方位脉冲压缩及位移补偿处理模块包括:用于对经相位补偿的原始成像数据进行方位脉冲压缩处理,再进行横向位移补偿处理,得到输出数据,所述横向位移Δx的表达式是:
其中,Δx表示横向位移,fdc表示所述多普勒中心频率,ν表示载机飞行速度,R表示雷达瞬时斜距,Θ表示载机俯仰角,所述λ是机载SAR***参数,ν、R和Θ均是所述惯导数据;
所述输出数据xa的表达式是:
xa=ifft(S(f)×H(f)×exp(j2πf×Δx/ν))
其中,所述S(f)表示回波信号频谱,H(f)表示滤波器频谱,f表示多普勒频率,所述S(f)由回波信号数据计算得出,H(f)由机载SAR***参数和惯导数据计算得出,所述多普勒频率是将时域信号转化为频率域信号的变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611068315.5A CN106772372B (zh) | 2016-11-29 | 2016-11-29 | 一种Ka波段机载SAR***的实时成像方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611068315.5A CN106772372B (zh) | 2016-11-29 | 2016-11-29 | 一种Ka波段机载SAR***的实时成像方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772372A CN106772372A (zh) | 2017-05-31 |
CN106772372B true CN106772372B (zh) | 2019-05-17 |
Family
ID=58904939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611068315.5A Active CN106772372B (zh) | 2016-11-29 | 2016-11-29 | 一种Ka波段机载SAR***的实时成像方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772372B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387836A (zh) * | 2018-12-11 | 2019-02-26 | 北京无线电测量研究所 | 一种逆合成孔径雷达的成像方法 |
CN111722221B (zh) * | 2020-06-10 | 2023-06-06 | 中国极地研究中心 | 机载雷达数据处理方法、装置、计算机设备及存储介质 |
CN112684447B (zh) * | 2020-12-10 | 2023-11-10 | 北京无线电测量研究所 | 一种毫米波机载sar实时成像优化方法及*** |
CN113376632B (zh) * | 2021-05-18 | 2023-12-15 | 南京航空航天大学 | 一种基于预处理和改进pfa的大斜视机载sar成像方法 |
CN114624656B (zh) * | 2022-03-08 | 2023-03-28 | 北京卫星信息工程研究所 | 星载多基线双波段雷达***及其空间目标探测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4617567A (en) * | 1984-09-28 | 1986-10-14 | The Boeing Company | Automatic motion compensation for residual motion in a synthetic aperture radar system |
CN102141611A (zh) * | 2010-12-29 | 2011-08-03 | 电子科技大学 | 斜视合成孔径雷达多普勒模糊数快速测定方法 |
CN102590812A (zh) * | 2012-01-13 | 2012-07-18 | 西安电子科技大学 | 基于调频连续波的sar实时成像方法 |
CN103197317A (zh) * | 2013-03-26 | 2013-07-10 | 西安电子科技大学 | 基于fpga的sar成像方法 |
CN103777201A (zh) * | 2012-10-19 | 2014-05-07 | 中国航天科工集团第二研究院二〇七所 | 基于gps数据的机载sar运动补偿方法 |
CN105403887A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子工程研究所 | 基于ins的机载sar实时运动补偿方法 |
-
2016
- 2016-11-29 CN CN201611068315.5A patent/CN106772372B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4617567A (en) * | 1984-09-28 | 1986-10-14 | The Boeing Company | Automatic motion compensation for residual motion in a synthetic aperture radar system |
CN102141611A (zh) * | 2010-12-29 | 2011-08-03 | 电子科技大学 | 斜视合成孔径雷达多普勒模糊数快速测定方法 |
CN102590812A (zh) * | 2012-01-13 | 2012-07-18 | 西安电子科技大学 | 基于调频连续波的sar实时成像方法 |
CN103777201A (zh) * | 2012-10-19 | 2014-05-07 | 中国航天科工集团第二研究院二〇七所 | 基于gps数据的机载sar运动补偿方法 |
CN103197317A (zh) * | 2013-03-26 | 2013-07-10 | 西安电子科技大学 | 基于fpga的sar成像方法 |
CN105403887A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子工程研究所 | 基于ins的机载sar实时运动补偿方法 |
Non-Patent Citations (1)
Title |
---|
机载SAR多普勒中心估计和调频率估计算法的实现;王凡等;《设计与研发》;20070930;56-60 |
Also Published As
Publication number | Publication date |
---|---|
CN106772372A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772372B (zh) | 一种Ka波段机载SAR***的实时成像方法和*** | |
CN109856635B (zh) | 一种csar地面动目标重聚焦成像方法 | |
CN106597437B (zh) | 高频机载调频连续波sar运动补偿及成像处理方法 | |
CN103487809B (zh) | 一种基于BP算法和时变基线的机载InSAR数据处理方法 | |
CN104597447B (zh) | 一种子孔径SAR大斜视改进Omega‑K成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN110488288B (zh) | 一种机载sar高分辨层析成像方法 | |
CN102176018B (zh) | 机械扫描雷达的多普勒波束锐化快速成像方法 | |
CN104898119B (zh) | 一种基于相关函数的动目标参数估计方法 | |
CN104316924A (zh) | 机载超高分辨率sar反投影图像的自聚焦运动补偿方法 | |
CN105372657A (zh) | 基于回波数据的视频合成孔径雷达运动补偿成像方法 | |
CN106054188A (zh) | 无人机合成孔径雷达成像的图像偏移自聚焦方法 | |
CN105607055A (zh) | 一种基于天线方向图的机载雷达单脉冲前视成像方法 | |
CN107271997B (zh) | 机载多通道cssar地面运动目标运动参数估计方法 | |
CN103777188A (zh) | 一种基于双频共轭处理的InISAR成像方法 | |
CN105572635A (zh) | 基于最小二乘法的单站无源快速定位方法 | |
CN101609152B (zh) | 一种基于单通道sar子图像对消的运动目标检测方法 | |
CN108710111A (zh) | 一种机载双基前视sar方位相位的二维空变校正方法 | |
CN103744081A (zh) | 一种机载圆迹合成孔径雷达高精度三维成像补偿方法 | |
CN109387836A (zh) | 一种逆合成孔径雷达的成像方法 | |
CN114325704B (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN104793196A (zh) | 基于距离徙动改进算法的实时sar成像方法 | |
CN103675815A (zh) | 一种在大斜视sar成像模式下对多普勒调频率进行精确估计的方法 | |
CN106772308A (zh) | 基于脉内干涉的太赫兹宽带雷达微多普勒解模糊方法 | |
CN109799502B (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 |