CN111999734A - 一种宽带斜视聚束sar两步成像方法及*** - Google Patents
一种宽带斜视聚束sar两步成像方法及*** Download PDFInfo
- Publication number
- CN111999734A CN111999734A CN202010888239.2A CN202010888239A CN111999734A CN 111999734 A CN111999734 A CN 111999734A CN 202010888239 A CN202010888239 A CN 202010888239A CN 111999734 A CN111999734 A CN 111999734A
- Authority
- CN
- China
- Prior art keywords
- azimuth
- signal
- sar
- spectrum
- imaging
- 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 80
- 208000004350 Strabismus Diseases 0.000 title claims abstract description 37
- 238000001228 spectrum Methods 0.000 claims abstract description 79
- 238000012545 processing Methods 0.000 claims abstract description 38
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 36
- 230000009466 transformation Effects 0.000 claims abstract description 36
- 230000006835 compression Effects 0.000 claims abstract description 29
- 238000007906 compression Methods 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000011084 recovery Methods 0.000 claims abstract description 13
- 230000001131 transforming effect Effects 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims description 7
- 238000013519 translation Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 8
- 238000007781 pre-processing Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 9
- 230000004044 response Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
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/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
- 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/904—SAR modes
- G01S13/9052—Spotlight mode
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)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种宽带斜视聚束SAR两步成像方法,包括:对原始回波数据进行距离向傅里叶变换,将信号变换到距离频域进行距离向脉冲压缩;对距离向脉冲压缩后的信号进行相位相乘;进行Chip‑Z变换完成变标傅里叶变换;得到方位频谱无模糊的信号;通过对方位频谱无模糊的信号进行相位相乘恢复得到原始回波信号的二维频谱;采用ω‑k算法进行SAR成像。本发明将宽带和斜视因素完美地融合在聚束SAR方位向预处理阶段,并创新地引入了Chirp‑Z变换完成信号的积分运算,从而整个处理过程与传统聚束SAR“两步”成像的算法构架完全相同,保证了不同SAR工作模式成像处理的一致性。
Description
技术领域
本发明涉及雷达技术领域,尤其是一种宽带斜视聚束SAR两步成像方法及***。
背景技术
星载高分辨率SAR通常都是通过聚束模式实现。对于复杂的成像几何关系、小场景的SAR成像,首选应是时域成像算法。但是,对于中小斜视角、较大场景的聚束SAR成像,仍然可以采用频域成像算法,除了计算效率方面的考虑以外,还有以下两个方面的因素:首先,星载SAR***一般都需要实现多种分辨率模式的SAR成像,低分辨SAR无一例外都是采用频域成像算法,因此,高分辨率SAR同样采用频域算法成像具有算法的继承性和一致性;其次,星载高分辨SAR成像多数都需要经过方位自聚焦处理,以补偿轨道、传播路径等误差造成的方位向散焦,虽然近些年随着并行计算技术的发展,SAR 研究者们也在逐渐使用时域成像的手段,但基于时域算法的自聚焦处理技术的研究仍然很少,相反,基于频域成像算法的自聚焦技术研究非常广泛且成熟。
对于星载聚束SAR,雷达***的脉冲重复频率(Pulse Repeat Frequency,PRF)只要求高于方位向的瞬时频率,不可能高于方位全多普勒带宽,因而,方位频谱混叠在所难免。对于聚束SAR频域成像算法,如何有效地去除方位频谱混叠是首先需要解决的一个问题。文献中有关方位频谱去混叠的方法主要有以下两种:一、方位向子孔径处理;二、方位向全孔径处理。前者通过在方位向选择较少的脉冲进行成像以使所处理信号的方位频谱小于***的瞬时带宽来避免方位频谱的混叠,得到的子孔径数据再融合成高分辨率SAR图像,但子孔径算法不可避免地会存在孔径划分、拼接及重叠部分的选择等问题,增加了信号处理的复杂性;后者先通过方位预处理,获得无模糊的方位频谱,再利用传统条带模式的方法进行成像,即所谓的“两步成像”算法(参考“Riccardo Lanari, Manlio Tesauro andEugenio Sansosti et al,Spotlight SAR Data Focusing based on a Two-StepProcessing Approach,IEEE Transactions on Geoscience and Remote sensing,2001,39(9):1993-2004.”)。
相对于子孔径处理繁琐的拼接过程,“两步成像”算法在去除方位频谱混叠后可以直接使用已有的条带SAR成像算法,处理过程简洁明了。但是,传统的“两步成像”算法是基于正侧视模型提出的,对于较大斜视角成像,该算法并不能完全去除场景内目标的方位频谱混叠现象。虽然可以通过在方位粗压缩时增加与doppler中心相关的线性相位消除由斜视带来的方位频谱平移(参考“韩冰等,斜视滑动聚束SAR成像的ECS算法,中国科学院研究生院学报,2012,29(5):674-680.”),但无法从根本上解决较大斜视角下全场景内目标频谱混叠的问题。因此,Daoxiang An等人提出了扩展的“两步成像”算法,用于斜视聚束SAR几何下进行频谱去混叠和成像(参考“Daoxiang An,Xiaotao Huang and Tian Jin et al,Extended Two-Step Focusing Approach for Squinted Spotlight SAR Imaging,IEEETransactions on Geoscience and Remote sensing,2012, 50(7):2889-2900.”)。该方法首先通过相位相乘去除了由距离频率引起的Doppler中心的变化,从而消除了二维频谱的平移和倾斜特性,同时校正了线性走动效应,然后采用传统的去斜函数进行去斜处理获得了无模糊的二维频谱,随后修正距离徙动算法完成了聚焦,但由于忽略了线性走动校正引起的方位空变问题,该方法只适用于小场景成像。韩晓磊分析了滑动聚束SAR方位带宽的组成,指出对于斜视情况距离方位耦合附加的带宽非常大,因此,作者在距离压缩的同时去除了线性走动,从而消除了频谱的二维倾斜,也消除了距离方位的耦合;传统去斜操作后,再在参考函数相乘之后、Stolt插值之前的方位时域补偿该线性走动,从而无需修改Stolt插值就可以完成聚焦(参考“韩晓磊等,斜视滑动聚束模式SAR成像算法研究,电子与信息学报,2013,35(12):2843-2849.”)。以上这两种方法类似,都是先在方位时域去除线性走动,然后利用传统去斜操作去除频谱模糊,不同点在于是否进行线性走动校正处理,导致了是否进行Stolt插值函数的校正。刘燕针对的不是斜视情况,而是正侧视下大距离带宽造成的频谱模糊,她分析了去斜后***带宽的大小,指出对于超高分辨率聚束SAR成像,方位预处理时如果采用传统的一维去斜函数,为保证方位频谱不混叠,SAR***所需的PRF要比采用二维去斜函数所需的PRF大得多,而对于星载SAR***,测绘带宽和***数据率的限制使得PRF的取值不可能很大,因此采用二维去斜函数具有非常明显的优势(参考“刘燕,高分辨率及新模式SAR成像算法研究,博士学位论文,2012.04,西安电子科技大学.”)。在正侧视的情况下,她提出方位预处理采用与距离频率相关的二维去斜函数代替传统的一维去斜函数,从而降低***PRF的需求,随后,通过不同距离频点上的补零处理,消除不同距离频点上输出采样间隔不一致的问题,其它步骤与传统“两步成像”方法类似。
文献中对于斜视聚束SAR成像的讨论很多集中在高频段,而在实际情况中,低频段宽带斜视聚束SAR有很多独特的问题。首先,低频段要达到很高的分辨率,方位向波束带宽要比高频段SAR高得多,造成合成孔径的长度非常大;其次,由于中心工作频率低,距离分辨率高,导致低频段发射信号的相对带宽非常大,甚至超过50%。通过试验发现,对于低频段宽带斜视聚束SAR频域一次性成像,任何文献中的方法单独使用都无法达到满意的聚焦效果。
发明内容
本发明的首要目的在于提供一种将斜视因素完美地融合在聚束SAR方位向预处理阶段,并创新地引入了Chirp-Z变换完成信号的积分运算,从而使整个处理过程与传统聚束SAR“两步”成像的算法构架完全相同,保证了不同SAR工作模式成像处理的一致性的宽带斜视聚束SAR两步成像方法。
为实现上述目的,本发明采用了以下技术方案:一种宽带斜视聚束SAR两步成像方法,包括下列顺序的步骤:
(1)对原始回波数据进行距离向傅里叶变换,将信号变换到距离频域进行距离向脉冲压缩;
(2)将距离向脉冲压缩后的信号与相位Φ1(fr,t)相乘,得到S1(fr,t)信号;
(3)对S1(fr,t)信号进行Chip-Z变换完成变标傅里叶变换,得到S2(fr,t′)信号;
(4)将S2(fr,t′)信号与相位补偿函数Φ2(fr,t)相乘,得到方位频谱无模糊的信号Sde(fr,t′);
(5)通过对方位频谱无模糊的信号Sde(fr,t′)进行相位相乘恢复得到原始回波信号的二维频谱S(fr,fa);
(6)采用ω-k算法进行SAR成像。
所述步骤(1)具体是指:
斜视条件下,聚束SAR的回波信号为:
其中,TP为脉冲宽度,为距离向时间,t为方位向时间,tc为聚束SAR成像的中心时刻;L为聚束SAR合成孔径长度;τ0是目标回波的延迟时间,且τ0=2R(t)/C,C为光速,r为tc时刻目标与雷达的斜距,为斜视角,v是平台的速度,Kr为距离向线性调频信号的调频率,fc为载波的中心频率, wr(.)、wa(.)分别是距离和方位向的窗函数;不失一般性,这里认为中心时刻tc为0,傅里叶变换后,式(1)变为:
傅里叶变换后,式(1)变为:
其中,fr为距离向频率;距离向脉冲压缩的函数为:
S(fr,t)与式(3)相乘得到S1(fr,t):
所述步骤(2)具体是指:
将式(4)中的R(t)在0处进行四阶泰勒展开,可得:
将上式代入式(4)中,则相位部分有:
利用考虑多普勒频率的二维去斜函数进行去斜:
距离压缩后的信号与去斜函数进行卷积后的信号为:
其中,t′为方位向新的时间,将Φ1(fr,t)与距离脉冲压缩后的信号S(fr,t)相乘:
得到S1(fr,t):
所述步骤(3)具体是指:
在完成步骤(2)后,式(8)变为:
这里,Δt、Δt′分别是原始和去斜后的方位向采样时间间隔,P是去斜后的方位向采样点数,N指自然数集;
式(11)中的积分利用Chirp-Z变换实现,Chip-Z变换后的信号为:
所述步骤(4)具体是指:
将S2(fr,t′)与相位补偿函数Φ2(fr,t)相乘:
根据式(11)得到方位频谱无模糊的信号Sde(fr,t′)。
所述步骤(5)具体是指:
信号Sde(fr,t′)方位傅里叶变换后的二维频谱为:
其中,S(fr,fa)是在新脉冲重复频率下原始SAR信号的二维频谱,FT(·)表示傅里叶变换;根据式(15),聚束SAR信号的二维频谱S(fr,fa)通过去斜后信号的二维频谱 Sde(fr,fa)乘上相位Φ3(fr,t)恢复;fa是方位向频率;
在所述步骤(2)中,采用式(7)中的去斜函数进行方位频谱去模糊时,去斜处理后信号的带宽见式(17):
在步骤(3)中,Chirp-Z变换的具体方法为:
(3a)将S1(fr,t)与变标核φ4相乘:
其中,n为表征所处方位点数的变量,实际处理时n=0,1,…,P-1;
(3b)对步骤(3a)的处理结果进行方位向快速傅里叶变换FFT;
(3c)对共轭变标核φ5进行FFT处理;
(3d)将步骤(3b)的结果与步骤(3c)的结果相乘;
(3e)将步骤(3d)的结果进行方位向逆快速傅里叶变换IFFT;
(3f)将步骤(3e)所得的结果与φ5进行相乘,完成Chirp-Z变换。
本发明的另一目的在于提供一种宽带斜视聚束SAR两步成像方法的***,包括:
距离压缩模块,完成对原始回波数据的距离向脉冲压缩;
考虑多普勒频率的二维去斜函数相乘模块,即方位去斜处理,在去除方位向调频率的同时,也去除信号在二维频谱上的平移和倾斜;
Chirp-Z变换模块,用Chirp-Z变换完成变标傅里叶变换,同时Chirp-Z变换转换为卷积运算,而卷积通过频域乘积实现;
考虑多普勒频率的二维残余相位补偿模块,用于补偿由于去斜操作所引入的残余相位;
考虑多普勒频率的二维频谱恢复模块,通过对去斜处理后的方位频谱无模糊的信号进行相位相乘,恢复得到原始回波信号的二维频谱;
传统条带SAR成像模块,完成聚束SAR方位向频谱去混叠和二维频谱恢复后,按照传统条带ω-k算法进行SAR成像。
由上述技术方案可知,本发明的有益效果为:第一,本发明在传统聚束SAR“两步”成像算法的基础上,考虑宽带和斜视两方面的综合因素对传统算法进行改进,打破了去斜处理只是针对方位相位的二次项进行操作的固定思维,创造性地将由于斜视造成的方位频谱展宽和由于距离向宽带造成的方位频谱展宽都融入到去斜函数中,统一进行频谱去混叠处理;第二,本发明无需单独考虑斜视引起的走动的去除和恢复的问题,从而不需要如文献中所示的修正传统距离徙动算法;第三,本发明将斜视因素完美地融合在聚束SAR方位向预处理阶段,即方位频谱去混叠步骤和频谱恢复步骤,并创新地引入了 Chirp-Z变换完成信号的积分运算,从而整个处理过程与传统聚束SAR“两步”成像的算法构架完全相同,保证了不同SAR工作模式成像处理的一致性。
附图说明
图1是本发明的方法流程图;
图2是Chirp-Z变换的方法流程图;
图3是星载聚束SAR***的仿真参数;
图4是场景中点目标的分布图;
图5是在PRF=2200Hz、K=1时,场景右下角目标压缩后的幅度图;
图6是在PRF=2200Hz、K=1时,场景右下角目标压缩后的等值线图;
图7是在PRF=2200Hz、K=1时,场景右下角目标压缩后的方位向点目标响应图;
图8是在PRF=2200Hz、K=1时,场景右下角目标压缩后的距离向点目标响应图;
图9是在PRF=2200Hz、K≠1时,场景右下角目标压缩后的幅度图;
图10是在PRF=2200Hz、K≠1时,场景右下角目标压缩后的等值线图;
图11是在PRF=2200Hz、K≠1时,场景右下角目标压缩后的方位向点目标响应图;
图12是在PRF=2200Hz、K≠1时,场景右下角目标压缩后的距离向点目标响应图;
图13是在PRF=3200Hz、K≠1时,场景右下角目标压缩后的幅度图;
图14是在PRF=3200Hz、K≠1时,场景右下角目标压缩后的等值线图;
图15是在PRF=3200Hz、K≠1时,场景右下角目标压缩后的方位向点目标响应图;
图16是在PRF=3200Hz、K≠1时,场景右下角目标压缩后的距离向点目标响应图。
具体实施方式
如图1所示,一种宽带斜视聚束SAR两步成像方法,包括下列顺序的步骤:
(1)对原始回波数据进行距离向傅里叶变换,将信号变换到距离频域进行距离向脉冲压缩;
(2)将距离向脉冲压缩后的信号与相位Φ1(fr,t)相乘,得到S1(fr,t)信号;
(3)对S1(fr,t)信号进行Chip-Z变换完成变标傅里叶变换,得到S2(fr,t′)信号;
(4)将S2(fr,t′)信号与相位补偿函数Φ2(fr,t)相乘,得到方位频谱无模糊的信号Sde(fr,t′);
(5)通过对方位频谱无模糊的信号Sde(fr,t′)进行相位相乘恢复得到原始回波信号的二维频谱S(fr,fa);
(6)采用ω-k算法进行SAR成像。
所述步骤(1)具体是指:
斜视条件下,聚束SAR的回波信号为:
其中,TP为脉冲宽度,为距离向时间,t为方位向时间,tc为聚束SAR成像的中心时刻;L为聚束SAR合成孔径长度;τ0是目标回波的延迟时间,且τ0=2R(t)/C,C为光速,r为tc时刻目标与雷达的斜距,为斜视角,v是平台的速度,Kr为距离向线性调频信号的调频率,fc为载波的中心频率, wr(.)、wa(.)分别是距离和方位向的窗函数;不失一般性,这里认为中心时刻tc为0,傅里叶变换后,式(1)变为:
其中,fr为距离向频率;距离向脉冲压缩的函数为:
S(fr,t)与式(3)相乘得到S1(fr,t):
所述步骤(2)具体是指:
将式(4)中的R(t)在0处进行四阶泰勒展开,可得:
将上式代入式(4)中,则相位部分有:
聚束SAR方位去斜时,传统的方法是去除雷达中心频率处方位向二次相位造成的调频率,去斜函数为:
在距离向带宽较大时,同时去除二次项调频率与距离频率的关系,此时的去斜函数为:
而对于斜视成像,斜视角会造成Doppler频率的偏移并随着距离频率倾斜,使得总的Doppler带宽变大,此时,通常的方法是在去斜处理时去除上式中的一次项以去除该带宽:
同时去除中心频率处的二次项,此时的去斜函数为:
在窄带假设下上述处理能够满足斜视聚束SAR方位向频谱去混叠的要求。但是,当存在斜视角并且***相对带宽达到35%或更高时,上述三种去斜函数都无法取得满意的压缩结果。在本申请中提出利用考虑多普勒频率的二维去斜函数进行去斜:
h(fr,t)=h1(fr,t)·h2(fr,t)
在去除方位向调频率的同时,也去除了信号在二维频谱的平移和倾斜。
距离压缩后的信号与去斜函数进行卷积后的信号为:
其中,t′为方位向新的时间,将Φ1(fr,t)与距离脉冲压缩后的信号S(fr,t)相乘:
得到S1(fr,t):
所述步骤(3)具体是指:
在完成步骤(2)后,式(8)变为:
这里,Δt、Δt′分别是原始和去斜后的方位向采样时间间隔,P是去斜后的方位向采样点数,N指自然数集;当存在系数K时,上述积分成为变标傅里叶变换(SCFT)。这里,提出利用Chirp-Z变换完成SCFT计算。Chirp-Z变换可以转换为卷积运算,而卷积又可以通过频域乘积快速实现,Chirp-Z变换后的信号为:
所述步骤(4)具体是指:
将S2(fr,t′)与相位补偿函数Φ2(fr,t)相乘:
以补偿由于去斜操作所造成的残余相位,就可以消除距离向大带宽造成的方位向频谱混叠,根据式(11)得到方位频谱无模糊的信号Sde(fr,t′)。
所述步骤(5)具体是指:
信号Sde(fr,t′)方位傅里叶变换后的二维频谱为:
其中,S(fr,fa)是在新脉冲重复频率下原始SAR信号的二维频谱,FT(·)表示傅里叶变换;根据式(15),聚束SAR信号的二维频谱S(fr,fa)通过去斜后信号的二维频谱 Sde(fr,fa)乘上相位Φ3(fr,t)恢复;fa是方位向频率;
完成步骤(1)至步骤(5)就完成了聚束SAR方位向频谱去混叠和二维频谱恢复,再按照传统条带SAR成像算法进行成像即可。对于宽带高分辨SAR成像,宜采用ω-k 算法。
需要说明的是,采用式(7)的去斜函数进行方位频谱去模糊时,虽然可以去除方位向相位的一次项和二次项,但对于相对带宽很大的SAR***来说,此时三次项和四次项引起的Doppler调频率仍将难以忽略。如果在PRF设计的时候不考虑这两项的带宽,仍然可能造成频谱的轻微混叠从而影响成像质量。
在所述步骤(2)中,采用式(7)中的去斜函数进行方位频谱去模糊时,去斜处理后信号的带宽见式(17):
它由方位向场景的瞬时带宽和整个合成孔径时间内三次项、四次项造成的剩余带宽组成,其中,Ls为成像的方位向场景,Br为距离向带宽,|t|max表示|t|的最大值,这里为采用式(7)的去斜函数处理时,雷达***所需的最小PRF必须大于BDe,否则将会影响成像质量。
在步骤(3)中,Chirp-Z变换的具体方法为:
(3a)将S1(fr,t)与变标核φ4相乘:
其中,n为表征所处方位点数的变量,实际处理时n=0,1,…,P-1;
(3b)对步骤(3a)的处理结果进行方位向快速傅里叶变换FFT;
(3c)对共轭变标核φ5进行FFT处理;
(3d)将步骤(3b)的结果与步骤(3c)的结果相乘;
(3e)将步骤(3d)的结果进行方位向逆快速傅里叶变换IFFT;
(3f)将步骤(3e)所得的结果与φ5进行相乘,完成Chirp-Z变换。
本***包括:
距离压缩模块,完成对原始回波数据的距离向脉冲压缩;
考虑多普勒频率的二维去斜函数相乘模块,即方位去斜处理,在去除方位向调频率的同时,也去除信号在二维频谱上的平移和倾斜;
Chirp-Z变换模块,用Chirp-Z变换完成变标傅里叶变换,同时Chirp-Z变换转换为卷积运算,而卷积通过频域乘积实现;
考虑多普勒频率的二维残余相位补偿模块,用于补偿由于去斜操作所引入的残余相位;
考虑多普勒频率的二维频谱恢复模块,通过对去斜处理后的方位频谱无模糊的信号进行相位相乘,恢复得到原始回波信号的二维频谱;
传统条带SAR成像模块,完成聚束SAR方位向频谱去混叠和二维频谱恢复后,按照传统条带ω-k算法进行SAR成像。
图1、图2给出了基于Chirp-z变换的宽带高分辨率聚束SAR成像的技术流程。为了验证该技术流程的有效性,利用仿真的P波段星载聚束SAR回波进行成像和分析。图 3是用于仿真的星载聚束SAR的***参数,中心斜视角为-12°。利用该参数仿真得到的距离向斜距分辨率约为0.6m,方位分辨率约为1m。成像的场景大小为10km×10km, 3行3列共九个点目标呈“田”字形均匀分布在场景中,见图4。为了排除轨道曲率的影响,这里假设卫星飞行的轨道为直线。根据方位向天线长度计算瞬时Doppler带宽, PRF=2200Hz可以达到***设计的要求;而根据式(17),在中心斜视角为-12°度时,三次项、四次项造成的去斜后剩余带宽最大可达到770Hz和70Hz。因此,这里分别采用 PRF=2200Hz和3200Hz进行成像仿真。
首先,PRF=2200Hz时,成像的去斜函数为h10(fr,t),亦即去斜函数的二次项没有考虑距离频率的影响,此时K=1。对于该场景成像,图像右下角目标的压缩结果见图5 至图8。可以看到,由于该去斜函数没有考虑距离频率的影响,不仅方位向压缩恶化,距离向压缩也受到了影响。
其次,PRF=2200Hz,去斜函数为h(fr,t),此时去斜函数的二次项考虑了距离频率的影响,此时K≠1。该区域右下角目标成像的结果见图9至图12。可以看到目标压缩的质量明显提高,方位向旁瓣明显不对称的现象也消失了,距离向压缩接近理想情况。但是方位向一侧的旁瓣有明显抬升现象,导致峰值旁瓣比(PSLR)达不到***要求,目标的主瓣也有轻微的展宽。
最后,PRF=3200Hz,去斜函数仍然为h(fr,t),右下角目标的压缩结果见图13至图16。此时,方位向旁瓣抬升和主瓣展宽现象消除了,两维分辨率也有所提高,图像压缩质量明显改善。
综上所述,本发明在传统聚束SAR“两步”成像算法的基础上,考虑宽带和斜视两方面的综合因素对传统算法进行改进,打破了去斜处理只是针对方位相位的二次项进行操作的固定思维,创造性地将由于斜视造成的方位频谱展宽和由于距离向宽带造成的方位频谱展宽都融入到去斜函数中,统一进行频谱去混叠处理。
Claims (6)
1.一种宽带斜视聚束SAR两步成像方法,其特征在于:包括下列顺序的步骤:
(1)对原始回波数据进行距离向傅里叶变换,将信号变换到距离频域进行距离向脉冲压缩;
(2)将距离向脉冲压缩后的信号与相位Φ1(fr,t)相乘,得到S1(fr,t)信号;
(3)对S1(fr,t)信号进行Chip-Z变换完成变标傅里叶变换,得到S2(fr,t′)信号;
(4)将S2(fr,t′)信号与相位补偿函数Φ2(fr,t)相乘,得到方位频谱无模糊的信号Sde(fr,t′);
(5)通过对方位频谱无模糊的信号Sde(fr,t′)进行相位相乘恢复得到原始回波信号的二维频谱S(fr,fa);
(6)采用ω-k算法进行SAR成像。
2.根据权利要求1所述的宽带斜视聚束SAR两步成像方法,其特征在于:所述步骤(1)具体是指:
斜视条件下,聚束SAR的回波信号为:
其中,TP为脉冲宽度,为距离向时间,t为方位向时间,tc为聚束SAR成像的中心时刻;L为聚束SAR合成孔径长度;τ0是目标回波的延迟时间,且τ0=2R(t)/C,C为光速,r为tc时刻目标与雷达的斜距,为斜视角,v是平台的速度,Kr为距离向线性调频信号的调频率,fc为载波的中心频率,wr(.)、wa(.)分别是距离和方位向的窗函数;不失一般性,这里认为中心时刻tc为0,傅里叶变换后,式(1)变为:
其中,fr为距离向频率;距离向脉冲压缩的函数为:
S(fr,t)与式(3)相乘得到S1(fr,t):
所述步骤(2)具体是指:
将式(4)中的R(t)在0处进行四阶泰勒展开,可得:
将上式代入式(4)中,则相位部分有:
利用考虑多普勒频率的二维去斜函数进行去斜:
距离压缩后的信号与去斜函数进行卷积后的信号为:
其中,t′为方位向新的时间,将Φ1(fr,t)与距离脉冲压缩后的信号S(fr,t)相乘:
得到S1(fr,t):
所述步骤(3)具体是指:
在完成步骤(2)后,式(8)变为:
这里,Δt、Δt′分别是原始和去斜后的方位向采样时间间隔,P是去斜后的方位向采样点数,N指自然数集;
式(11)中的积分利用Chirp-Z变换实现,Chip-Z变换后的信号为:
所述步骤(4)具体是指:
将S2(fr,t′)与相位补偿函数Φ2(fr,t)相乘:
根据式(11)得到方位频谱无模糊的信号Sde(fr,t′)。
6.实施权利要求1至5中任一项所述宽带斜视聚束SAR两步成像方法的***,其特征在于:包括:
距离压缩模块,完成对原始回波数据的距离向脉冲压缩;
考虑多普勒频率的二维去斜函数相乘模块,即方位去斜处理,在去除方位向调频率的同时,也去除信号在二维频谱上的平移和倾斜;
Chirp-Z变换模块,用Chirp-Z变换完成变标傅里叶变换,同时Chirp-Z变换转换为卷积运算,而卷积通过频域乘积实现;
考虑多普勒频率的二维残余相位补偿模块,用于补偿由于去斜操作所引入的残余相位;
考虑多普勒频率的二维频谱恢复模块,通过对去斜处理后的方位频谱无模糊的信号进行相位相乘,恢复得到原始回波信号的二维频谱;
传统条带SAR成像模块,完成聚束SAR方位向频谱去混叠和二维频谱恢复后,按照传统条带ω-k算法进行SAR成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010888239.2A CN111999734B (zh) | 2020-08-28 | 2020-08-28 | 一种宽带斜视聚束sar两步成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010888239.2A CN111999734B (zh) | 2020-08-28 | 2020-08-28 | 一种宽带斜视聚束sar两步成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111999734A true CN111999734A (zh) | 2020-11-27 |
CN111999734B CN111999734B (zh) | 2022-02-08 |
Family
ID=73465977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010888239.2A Active CN111999734B (zh) | 2020-08-28 | 2020-08-28 | 一种宽带斜视聚束sar两步成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111999734B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113341413A (zh) * | 2021-06-28 | 2021-09-03 | 南京航空航天大学 | 一种高分辨率斜视sar回波多普勒频谱恢复方法及*** |
CN113359132A (zh) * | 2021-04-30 | 2021-09-07 | 西安电子科技大学 | 星载斜视合成孔径雷达实时成像方法及装置 |
CN113900100A (zh) * | 2021-12-09 | 2022-01-07 | 中国科学院空天信息创新研究院 | 一种基于时域解混叠的聚束sar全孔径处理方法 |
CN114942443A (zh) * | 2022-07-25 | 2022-08-26 | 中国人民解放军国防科技大学 | 基于mimo-sar的介质目标快速成像方法和装置 |
CN116106852A (zh) * | 2023-04-12 | 2023-05-12 | 中国人民解放军63921部队 | 机载主杂波频道的确定方法、装置及电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102788972A (zh) * | 2011-05-18 | 2012-11-21 | 南京航空航天大学 | 一种适用于超高分辨率sar成像的自聚焦方法 |
CN103576147A (zh) * | 2012-08-02 | 2014-02-12 | 中国科学院电子学研究所 | 合成孔径雷达大斜视模式下成像方法 |
CN103792534A (zh) * | 2012-11-01 | 2014-05-14 | 南京航空航天大学 | 一种基于先验相位结构知识的sar两维自聚焦方法 |
CN105842694A (zh) * | 2016-03-23 | 2016-08-10 | 中国电子科技集团公司第三十八研究所 | 一种基于ffbp sar成像的自聚焦方法 |
CN109782277A (zh) * | 2017-11-14 | 2019-05-21 | 中电科海洋信息技术研究院有限公司 | 变pri的斜视聚束sar成像方法、装置、设备及存储介质 |
CN109870686A (zh) * | 2019-01-03 | 2019-06-11 | 西安电子科技大学 | 一种基于改进斜距模型的曲线轨迹sar斜视成像方法 |
-
2020
- 2020-08-28 CN CN202010888239.2A patent/CN111999734B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102788972A (zh) * | 2011-05-18 | 2012-11-21 | 南京航空航天大学 | 一种适用于超高分辨率sar成像的自聚焦方法 |
CN103576147A (zh) * | 2012-08-02 | 2014-02-12 | 中国科学院电子学研究所 | 合成孔径雷达大斜视模式下成像方法 |
CN103792534A (zh) * | 2012-11-01 | 2014-05-14 | 南京航空航天大学 | 一种基于先验相位结构知识的sar两维自聚焦方法 |
CN105842694A (zh) * | 2016-03-23 | 2016-08-10 | 中国电子科技集团公司第三十八研究所 | 一种基于ffbp sar成像的自聚焦方法 |
CN109782277A (zh) * | 2017-11-14 | 2019-05-21 | 中电科海洋信息技术研究院有限公司 | 变pri的斜视聚束sar成像方法、装置、设备及存储介质 |
CN109870686A (zh) * | 2019-01-03 | 2019-06-11 | 西安电子科技大学 | 一种基于改进斜距模型的曲线轨迹sar斜视成像方法 |
Non-Patent Citations (3)
Title |
---|
RICCARDO LANARI ET AL.: ""Spotlight SAR Data Focusing Based on a Two-Step Processing Approach"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
怀园园 等: ""一种基于方位谱重采样的大斜视子孔径SAR成像改进Omega-K算法"", 《电子与信息学报》 * |
鲁加国 等: """一步一停"假设不成立时的星载聚束SAR成像"", 《雷达科学与技术》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113359132A (zh) * | 2021-04-30 | 2021-09-07 | 西安电子科技大学 | 星载斜视合成孔径雷达实时成像方法及装置 |
CN113341413A (zh) * | 2021-06-28 | 2021-09-03 | 南京航空航天大学 | 一种高分辨率斜视sar回波多普勒频谱恢复方法及*** |
CN113341413B (zh) * | 2021-06-28 | 2024-04-09 | 南京航空航天大学 | 一种高分辨率斜视sar回波多普勒频谱恢复方法及*** |
CN113900100A (zh) * | 2021-12-09 | 2022-01-07 | 中国科学院空天信息创新研究院 | 一种基于时域解混叠的聚束sar全孔径处理方法 |
CN114942443A (zh) * | 2022-07-25 | 2022-08-26 | 中国人民解放军国防科技大学 | 基于mimo-sar的介质目标快速成像方法和装置 |
CN114942443B (zh) * | 2022-07-25 | 2022-10-21 | 中国人民解放军国防科技大学 | 基于mimo-sar的介质目标快速成像方法和装置 |
CN116106852A (zh) * | 2023-04-12 | 2023-05-12 | 中国人民解放军63921部队 | 机载主杂波频道的确定方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111999734B (zh) | 2022-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111999734B (zh) | 一种宽带斜视聚束sar两步成像方法 | |
CN105842694B (zh) | 一种基于ffbp sar成像的自聚焦方法 | |
EP1503223B1 (en) | Estimation and correction of phase for focusing search mode SAR images formed by range migration algorithm | |
CN104007440B (zh) | 一种加速分解后向投影聚束合成孔径雷达成像方法 | |
CN104316924B (zh) | 机载超高分辨率sar反投影图像的自聚焦运动补偿方法 | |
CN102680974B (zh) | 一种星载滑动聚束合成孔径雷达的信号处理方法 | |
CN111856461B (zh) | 基于改进pfa的聚束sar成像方法及其dsp实现 | |
CN108710111B (zh) | 一种机载双基前视sar方位相位的二维空变校正方法 | |
US7551119B1 (en) | Flight path-driven mitigation of wavefront curvature effects in SAR images | |
CN107918124A (zh) | 带有方位空变校正的机载大斜视高分辨sar成像方法 | |
CN109870686A (zh) | 一种基于改进斜距模型的曲线轨迹sar斜视成像方法 | |
Yang et al. | A novel three-step image formation scheme for unified focusing on spaceborne SAR data | |
CN108061890B (zh) | 一种sar成像方法 | |
CN102608576A (zh) | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 | |
CN103630905A (zh) | 阵列天线sar极坐标交叠子孔径成像方法 | |
CN114545411A (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
CN110361733B (zh) | 一种基于时频联合重采样的中轨sar大斜视成像方法 | |
CN113702974A (zh) | 一种机载/弹载合成孔径雷达图像快速优化方法 | |
CN102788978B (zh) | 一种斜视星机双基地合成孔径雷达成像方法 | |
CN112649808B (zh) | 基于移变构型的双基前视sar波数域成像方法 | |
CN112415512B (zh) | 基于进退法和黄金分割法的sar运动目标聚焦方法 | |
CN106646471B (zh) | 基于方位空变误差补偿的机载高分辨sar成像方法 | |
CN115015920A (zh) | 一种基于距离空变频谱矫正的快速后向投影成像方法 | |
CN116719027A (zh) | 一种双向滑动聚束模式下的星机双基sar成像方法 | |
CN111638518B (zh) | 基于去调频的滑动聚束模式sar成像处理方法及装置 |
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 |