CN106908789B - 一种基于空间域支撑区融合的sar成像方法 - Google Patents
一种基于空间域支撑区融合的sar成像方法 Download PDFInfo
- Publication number
- CN106908789B CN106908789B CN201710111407.5A CN201710111407A CN106908789B CN 106908789 B CN106908789 B CN 106908789B CN 201710111407 A CN201710111407 A CN 201710111407A CN 106908789 B CN106908789 B CN 106908789B
- Authority
- CN
- China
- Prior art keywords
- sub
- aperture
- data
- azimuth
- domain
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 86
- 230000004927 fusion Effects 0.000 title claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 112
- 238000001914 filtration Methods 0.000 claims abstract description 32
- 230000035485 pulse pressure Effects 0.000 claims abstract description 27
- 230000008859 change Effects 0.000 claims description 24
- 238000001514 detection method Methods 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 19
- 230000001427 coherent effect Effects 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 34
- 238000010586 diagram Methods 0.000 description 18
- 238000004088 simulation Methods 0.000 description 13
- 238000004422 calculation algorithm Methods 0.000 description 11
- 230000001419 dependent effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000005526 G1 to G0 transition Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004576 sand Substances 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
- 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
-
- 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
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成像方法,其主要思路为:获取子孔径回波数据;对子孔径回波数据依次进行2‑维快速傅里叶变换、距离匹配滤波处理、相位补偿处理、扩展stolt插值处理、距离向脉压处理,进而计算得到距离向脉压处理后的距离位置域‑方位波数域子孔径数据;对距离向脉压处理后的距离位置域‑方位波数域子孔径数据进行去斜处理,得到去斜处理后的距离位置域‑方位波数域子孔径数据,然后对去斜处理后的距离位置域‑方位波数域子孔径数据进行方位向逆快速傅里叶变换,得到2‑维位置域子孔径数据,然后计算每一个子孔径的成像数据;进而得到方位向聚焦的全孔径高分辨率SAR成像。
Description
技术领域
本发明属于SAR成像技术领域,特别涉及一种基于空间域支撑区融合的SAR成像方法,适用于机载或星载平台的SAR成像处理。
背景技术
SAR(合成孔径雷达)作为一种主动探测工具能够对观测场景进行高分辨微波成像,在遥感领域得到广泛应用,其通过发射宽频带的信号获得距离维的高分辨特性,而方位维则利用平台运动产生的多普勒信息实现高分辨。
近年来随着各种成像算法的成熟和发展,兼顾成像质量的同时,实时性成为很多成像系统的必备要求,子孔径概念的提出一定程度解决了成像质量与成像效率这对矛盾;由于子孔径数据成像具有录取时间短、积累的相位误差小等特点,可以减小运动补偿复杂度、计算量和存储量,并实现实时成像;因此,通过子孔径数据成像,首先快速扫描场景,得到低分辨的场景图像,而当获取到感兴趣的目标场景时,通过相关子孔径数据间的融合,实现感兴趣目标的高分辨成像;这样在节省时间的同时,又实现了特点场景的高分辨成像。
针对SAR数据的子孔径成像,人们提出了很多算法,但这些算法存在或多或少的近似;但其中的Omega-K算法通过Stolt插值能够实现距离及方位向的无近似解耦合,并由二维快速逆傅里叶变换(IFFT)完成距离及方位向的聚焦,因此具有更高的精准度。
然而传统的Omega-K算法由于其特殊的插值方式,使得方位向聚焦与解耦合同时完成;该传统的Omega-K算法一方面令成像算法难以与现有运动补偿方法相结合,另一方面,各子孔径间将无法进行方位向的相干融合以实现高分辨成像;同时,现有Omega-K算法主要针对于全孔径成像,由于方位向数据聚焦在空间域,使得必须对子孔径数据进行方位向补零操作以扩展其支撑区,大大增加了运算量,降低了算法的实时性。
发明内容
针对上述现有技术存在的不足,本发明的目的在于提出一种基于空间域支撑区融合的SAR成像方法,该种基于空间域支撑区融合的SAR成像方法能够在实现距离空间域和方位波数域的子孔径成像的同时,对各子孔径方位向数据进行相干性处理,实现方位向支撑区在空间域的融合,以满足高分辨SAR成像的需求。
为达到上述技术目的,本发明采用如下技术方案予以实现。
一种基于空间域支撑区融合的SAR成像方法,包括以下步骤:
步骤1,确定机载合成孔径雷达,所述机载合成孔径雷达发射信号并接收回波信号,并将回波信号记为子孔径回波数据;
步骤2,对子孔径回波数据进行2-维快速傅里叶变换,得到2-维波数域子孔径数据;
步骤3,对2-维波数域子孔径数据进行距离匹配滤波处理,得到距离匹配滤波后的2-维波数域子孔径数据;
步骤4,对距离匹配滤波后的2-维波数域子孔径数据进行相位补偿处理,得到相位补偿处理后的2-维波数域子孔径数据。
步骤5,对相位补偿处理后的2-维波数域子孔径数据进行扩展stolt插值处理,得到扩展stolt插值处理后的2-维波数域子孔径数据;
步骤6,对扩展stolt插值处理后的2-维波数域子孔径数据进行距离向脉压处理,即进行距离向逆快速傅里叶变换,进而计算得到距离向脉压处理后的距离位置域-方位波数域子孔径数据;
步骤7,对距离向脉压处理后的距离位置域-方位波数域子孔径数据进行去斜处理,得到去斜处理后的距离位置域-方位波数域子孔径数据,然后对去斜处理后的距离位置域-方位波数域子孔径数据进行方位向逆快速傅里叶变换,得到2-维位置域子孔径数据;
步骤8,初始化:令i∈{1,2,…,N},N为子孔径个数,i的初始值为1;
对第i个子孔径,取2-维位置域子孔径数据S7(Rr,X)的相位项,记为第i个子孔径数据的相位项,然后对第i个子孔径数据的相位项进行相干叠加处理及方位向空间域支撑区融合,即对第i个子孔径数据的相位项依次进行方位匹配滤波处理、常数相位补偿、方位位置补偿、方位向快速傅里叶变换、全局化处理、方位向逆快速傅里叶变换、相干叠加、方位波数域聚焦成像,进而得到i个子孔径的成像数据;
步骤9,分别令i取1至N,分别得到第1个子孔径的成像数据至第N个子孔径的成像数据,然后对第1个子孔径的成像数据至第N个子孔径的成像数据进行叠加,得到全孔径雷达成像数据,最后对全孔径雷达成像数据进行方位向快速傅里叶变换,进而得到方位向聚焦的全孔径高分辨率SAR成像。
本发明的有益效果为:
1)本发明结合改进的Omega-K算法,一方面进行了距离空间域和方位波数域的子孔径成像,实现子孔径快速实时成像;另一方面对各子孔径方位向数据进行相干性处理,实现了向方位向支撑区在空间域的融合,以满足特定场景的高分辨SAR成像需求。
2)传统的Omega-K算法通过映射关系进行stolt插值实现距离方位解耦合操作,同时方位向匹配滤波一并完成,最后通过2维逆快速傅里叶变换(IFFT)完成距离和方位向的同时聚焦;由于子孔径较短的方位空间域支撑区,将会使得方位向目标位置发生混叠;而本发明方法使用扩展stolt插值实现解耦合操作和方位调制项的独立化,并通过去斜操作对方位向进行波数域聚焦,能够解决子孔径数据方位聚焦位置混叠的现象。
3)本发明方法考虑到子孔径间的差异情况,在进行不同子孔径图像相干叠加时,进行相应的相干性处理,消除了不同子孔径中相位-位置变化率分布线的错位现象,使得方位位置一致,并且方位快速傅里叶变换(FFT)后,同一点目标在不同子孔径中将会聚焦到相同的位置。
4)本发明方法能够使得各子孔径数据从各自的局部坐标系统一到全孔径的全局坐标系中,且子孔径数据的空间域支撑区处于各自正确的区间位置,融合为全孔径数据的完整支撑区,并且在所有子孔径数据的相干叠加后进行方位快速傅里叶变换(FFT),即可得到全孔径高分辨的SAR图像。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1为本发明的一种基于空间域支撑区融合的SAR成像方法流程图;
图2为本发明的三维空间直角坐标系OXYZ中的SAR成像几何模型示意图;
图3a为本发明实施例中距离向脉压处理后的距离位置域-方位波数域子孔径数据对应的相位-波数变化率分布线示意图;
图3b为本发明实施例中去斜校正函数对应的相位-波数变化率分布线示意图;
图3c为本发明实施例中去斜处理后的距离位置域-方位波数域子孔径数据的相位-波数变化率分布线示意图
图3d为本发明实施例中2-维位置域子孔径数据的相位-位置变化率分布线示意图;
图4为本发明中对2-维位置域子孔径数据进行相干叠加处理及方位向空间域支撑区融合的流程图;
图5a为本发明实施例中距离聚焦处理后,处于不同方位位置的点目标A、点目标B、点目标各自的相位-位置变化率分布线示意图;
图5b为本发明实施例中方位位置补偿后2-维位置域子孔径数据中第i个子孔径数据的二维位置域信号的相位-位置变化率分布线示意图;
图5c为本发明实施例中把各子孔径数据从各自的局部坐标系统一到全孔径的全局坐标系后的相位-位置变化率分布线示意图;
图6a为仿真实验中使用本发明得出的单个点目标的二维升采样后的单子孔径结果图;
图6b为图6a对应的方位向剖面图;
图7a为仿真实验中使用本发明得出的单个点目标的二维升采样后的两子孔径融合结果图;
图7b为图7a对应的方位向剖面图;
图8a为仿真实验中使用本发明得出的单个点目标的二维升采样后的四子孔径融合结果图;
图8b为图8a对应的方位向剖面图;
图9a为仿真实验中使用本发明得出的单个点目标的二维升采样后的八子孔径融合结果图;
图9b为图9a对应的方位向剖面图。
具体实施方式
参照图1,为本发明的一种基于空间域支撑区融合的SAR成像方法流程图;其中基于空间域支撑区融合的SAR成像方法,包括以下步骤:
步骤1,确定机载合成孔径雷达,所述机载合成孔径雷达发射信号并接收回波信号,并将回波信号记为子孔径回波数据。
具体地,参照图2,为本发明的三维空间直角坐标系XOYZ中的SAR成像几何模型示意图;其中在三维坐标系XOYZ中,机载合成孔径雷达工作于条带模式,机载合成孔径雷达的载机以水平速度v沿X方向(即机载合成孔径雷达的方位向)飞行,H为机载合成孔径雷达的载机飞行高度,点P为机载合成孔径雷达照射场景中的任意一个散射点,记为点目标P;Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,点B为机载合成孔径雷达照射检测范围内场景的中心点,XP为点目标P沿方位向与点B之间的距离;机载合成孔径雷达以重复周期T、发射时刻tm依次发射信号,tm=mT,m为正整数变量。
当经过任意时间后,由几何关系计算得到机载合成孔径雷达与点P之间的瞬时斜距R(X;Rb)为:
其中,X为机载合成孔径雷达的方位向位置变量,X=vtm,tm=mT,T为机载合成孔径雷达发射信号的重复周期,m为正整数;此处使用空间位置域斜距表达式以与后续波数域表述相对应;根据机载合成孔径雷达与点目标P之间的瞬时斜距R(X;Rb),可以容易得到点目标P的回波信号表达式;机载合成孔径雷达发射信号并接收回波信号;本发明中,机载合成孔径雷达发射线性调频(LFM)信号并接收回波信号,然后对接收到的回波信号进行相干检波,即乘以exp(-j2πfct),fc为机载合成孔径雷达发射信号的载波频率,t为时间变量;进而得到机载合成孔径雷达接收的基带回波信号,所述机载合成孔径雷达接收的基带回波信号为子孔径回波数据S0(Rr,X;Rb),其表达式为:
其中,Rr为机载合成孔径雷达的距离向位置变量,C为光速,t为时间变量;R(X;Rb)为机载合成孔径雷达与点目标P之间的瞬时斜距,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,X为机载合成孔径雷达的方位向位置变量,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,γ为机载合成孔径雷达发射信号的调频率,λ为机载合成孔径雷达发射信号的载波波长,wr(Rr)为机载合成孔径雷达发射信号的距离窗函数,wa(X)为机载合成孔径雷达发射信号的方位窗函数,exp表示指数函数,j为虚数单位。
步骤2,采用驻定相位原理,对子孔径回波数据S0(Rr,X;Rb)进行2-维快速傅里叶变换(FFT),得到2-维波数域子孔径数据S1(Kr,Kx)。
其中,所述2-维波数域子孔径数据S1(Kr,Kx)的表达式为:
其中,Kr为机载合成孔径雷达的距离波数变量,Kr=4πf/C,f表示机载合成孔径雷达发射信号的频率在点目标P上的频率分量,f=fc+fb,fb为机载合成孔径雷达发射信号的频率在点目标P上的基频,
Kr=Krc+△Kr,Krc为机载合成孔径雷达的距离波数中心值,
Krc=4πfc/C,fc为机载合成孔径雷达发射信号的载波频率,△Kr为差值,△Kr∈[-2πγTp/C,2πγTp/C],C为光速,γ为机载合成孔径雷达发射信号的调频率,Tp为机载合成孔径雷达发射信号的脉冲宽度,∈表示属于;W(Kr)为机载合成孔径雷达发射信号的距离窗函数的波数域形式,Kx为机载合成孔径雷达的方位波数变量,Kx=2πfd/v,v为机载合成孔径雷达的载机水平速度,fd为机载合成孔径雷达接收到的回波信号的多普勒频率;W(Kx)为机载合成孔径雷达发射信号的方位窗函数的波数域形式,C为光速,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,点B为机载合成孔径雷达照射检测范围内场景的中心点,X为机载合成孔径雷达的方位向位置变量,exp表示指数函数,j为虚数单位。
步骤3,对2-维波数域子孔径数据S1(Kr,Kx)进行距离匹配滤波处理,得到距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx)。
具体地,确定距离匹配滤波函数HRMF(Kr),其表达式为:
将2-维波数域子孔径数据S1(Kr,Kx)乘以距离匹配滤波函数HRMF(Kr),消除距离二次项,得到距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx)。
所述距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx)的表达式为:
其中,W(Kr)为机载合成孔径雷达发射信号的距离窗函数的波数域形式,W(Kx)为机载合成孔径雷达发射信号的方位窗函数的波数域形式,Kx为机载合成孔径雷达的方位波数变量,Kr为机载合成孔径雷达的距离波数变量,Krc为机载合成孔径雷达的距离波数中心值,C为光速,γ为机载合成孔径雷达发射信号的调频率,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,点B为机载合成孔径雷达照射检测范围内场景的中心点,exp表示指数函数,j为虚数单位。
步骤4,对距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx)进行相位补偿处理,得到相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)。
具体地,相位补偿处理是将点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离Rb以及距离徙动量均调整到以机载合成孔径雷达照射检测范围内场景的中心点B为参考;其中将距离徙动量记为Rq,Rq=Re-Rb,Re为机载合成孔径雷达照射检测范围内场景中的边缘散射点与机载合成孔径雷达的载机飞行轨迹之间的距离,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,点目标P为机载合成孔径雷达照射场景中的任意一个散射点。
确定相位补偿函数HBulk(Kr,Kx),其表达式为:
其中,Rs表示机载合成孔径雷达照射检测范围内场景的中心点B对应的参考距离,Kx为机载合成孔径雷达的方位波数变量,Kr为机载合成孔径雷达的距离波数变量,exp表示指数函数,j为虚数单位。
将距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx)乘以相位补偿函数HBulk(Kr,Kx),得到相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)。
所述相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)的表达式为:
其中,W(Kr)为机载合成孔径雷达发射信号的距离窗函数的波数域形式,W(Kx)为机载合成孔径雷达发射信号的方位窗函数的波数域形式,Kx为机载合成孔径雷达的方位波数变量,Kr为机载合成孔径雷达的距离波数变量,Rs表示机载合成孔径雷达照射检测范围内场景的中心点B对应的参考距离,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,点B为机载合成孔径雷达照射检测范围内场景的中心点,exp表示指数函数,j为虚数单位。
步骤5,对相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)进行扩展stolt插值处理,得到扩展stolt插值处理后的2-维波数域子孔径数据S4(Ky,Kx);Ky表示stolt插值处理时引入的机载合成孔径雷达的距离波数变量,
具体地,对相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)进行扩展stolt插值处理,其中扩展stolt插值处理的映射表达式为进而计算得到扩展stolt插值处理后的2-维波数域子孔径数据S4(Ky,Kx),其表达式为:
S4(Ky,Kx)=W(Ky)W(Kx)exp(-jKxXP)
其中,Ky表示stolt插值处理时引入的机载合成孔径雷达的距离波数变量,
该步骤进行扩展stolt插值处理的实质是对相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)进行距离方位解耦合操作,同时将相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)中的机载合成孔径雷达的距离波数变量Kr和机载合成孔径雷达的方位波数变量Kx分离开,使得扩展stolt插值处理后的2-维波数域子孔径数据S4(Ky,Kx)的方位向能够单独处理并适用于子孔径数据成像。
步骤6,对扩展stolt插值处理后的2-维波数域子孔径数据进行距离向脉压处理,即进行距离向逆快速傅里叶变换(IFFT),实现距离向聚焦,进而计算得到距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)。
具体地,所述距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)的表达式为:
S5(Rr,Kx)=W(Kx)sinc{A[Rr-(Rb-Rs)]}
其中,A为距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)的常数项,W(Kx)表示机载合成孔径雷达发射信号的方位窗函数的波数域形式,exp表示指数函数,j为虚数单位,a=A[Rr-(Rb-Rs)];Rs表示机载合成孔径雷达照射检测范围内场景的中心点B对应的参考距离。
步骤7,对距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)进行去斜处理,得到去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx),然后对去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx)进行方位向逆快速傅里叶变换(IFFT),得到2-维位置域子孔径数据S7(Rr,X)。
具体地,对距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)进行去斜处理的作用是消除距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)对应的相位-波数变化率分布线沿相位波数变化率轴的错位,其中距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)对应的相位-波数变化率分布线如图3a所示。
在图3a所示的平面直角坐标系中,横轴为机载合成孔径雷达的方位波数变量Kx,纵轴为距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)对应的相位-波数变化率 X-p对应点目标A'的方位位置,O对应点目标B'的方位位置,Xp对应点目标C'的方位位置,Xsub为距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)对应的位置支撑区,即距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)在波数域的取值范围;超出距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)对应的位置支撑区的分布线将出现混叠,如图3a中的虚线所示,混叠后的位置如图3a中的实线所示。
对混叠后的点目标进行方位向逆快速傅里叶变换(IFFT),会造成对应点目标方位位置在纵坐标的投影出现在错误位置上,无法反映对应点目标的真实方位位置关系,其中点目标为点目标A'、点目标B'或点目标C';本发明方法能够有效解决这一问题。
点目标A'、点目标B'和点目标C'为距离向位置分别相同的三个点目标,即点目标A'、点目标B'和点目标C'与机载合成孔径雷达的载机飞行轨迹之间的最近距离分别相同;点目标A'、点目标B'和点目标C'各自的相位-波数变化率分布线不仅在方位波束轴的投影都是错开的,而且在相位-波数变化率轴的投影也都是错开的,所以在去斜处理时,引入去斜校正函数H3(Kx),
然后,将距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)乘以去斜校正函数H3(Kx),进而得到去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx);所述去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx)的表达式为:
其中,a=A[Rr-(Rb-Rs)],A为距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)的常数项;去斜校正函数H3(Kx)对应的相位-波数变化率分布线如图3b所示,去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx)的相位-波数变化率分布线如图3c所示,图3c中的k表示斜率,
对去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx)进行方位向逆快速傅里叶变换(IFFT)处理,得到2-维位置域子孔径数据S7(Rr,X);所述2-维位置域子孔径数据S7(Rr,X)的表达式为:
其中,wa(X)为机载合成孔径雷达发射信号的方位窗函数。
2-维位置域子孔径数据S7(Rr,X)的相位-位置变化率分布线如图3d所示,在图3d的平面直角坐标系中,横轴为机载合成孔径雷达的方位向位置变量X,纵轴为2-维位置域子孔径数据S7(Rr,X)的相位-位置变化率 此时,点目标A'、点目标B'和点目标C'各自的方位位置域的相位-位置变化率分布线在位置轴的投影都是对齐的,并且具有相同的斜率
步骤8,初始化:令i∈{1,2,…,N},N为子孔径个数,i的初始值为1。
对第i个子孔径,取2-维位置域子孔径数据S7(Rr,X)的相位项,记为第i个子孔径数据的相位项,然后对第i个子孔径数据的相位项进行相干叠加处理及方位向空间域支撑区融合,即对第i个子孔径数据的相位项依次进行方位匹配滤波处理、常数相位补偿、方位位置补偿、方位向快速傅里叶变换、全局化处理、方位向逆快速傅里叶变换、相干叠加、方位波数域聚焦成像,进而得到i个子孔径的成像数据Φi4(X)。
具体地,参照图4,为本发明中对2-维位置域子孔径数据S7(Rr,X)进行相干叠加处理及方位向空间域支撑区融合的流程图;为了区分全孔径数据与子孔径数据之间的差异,此处重新定义几个方位向变量:
参照图2所示的本发明的三维空间直角坐标系OXYZ中的SAR成像几何模型示意图,其中机载合成孔径雷达载机从左向右飞行,且机载合成孔径雷达以重复周期T、发射时刻为tm依次发射信号;本实施例中机载合成孔径雷达发射的信号分别为圆锥体波束,tm=mT,m为正整数。
第m个发射时刻tm时机载合成孔径雷达发射的信号照射到点目标P并产生回波数据,记为第m个子孔径,且第m个子孔径对应的回波数据为第m个子孔径数据。
令m分别取1至N,进而分别得到第1个子孔径至第N个子孔径,以及第1个子孔径数据至第N个子孔径数据;然后分别将第1个子孔径至第N个子孔径记为全孔径,将第1个子孔径数据至第N个子孔径数据记为全孔径数据。
建立平面直角坐标系,其横轴为机载合成孔径雷达的方位向位置变量X,纵轴为相位-位置变化率 λ为机载合成孔径雷达发射信号的载波波长,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离;并将该平面直角坐标系作为全孔径坐标系,其中全孔径坐标系中对应包含N个子孔径坐标系,将第i个子孔径中心在全孔径坐标系下的位置记为Ci,i∈{1,2,…,N},N为子孔径个数。
Xf为全孔径方位向采样位置,Xf∈[-L/2,L/2],L为合成孔径长度,Xs为每个子孔径方位向采样位置,Xs∈[-Ls/2,Ls/2],Ls为每个子孔径长度,Ci为第i个子孔径中心在全孔径坐标系下的位置,i∈{1,2,…,N},N为子孔径个数,L=NLs。
因此,对第i个子孔径,取2-维位置域子孔径数据S7(Rr,X)的相位项,记为第i个子孔径数据的相位项Φi0(X),其表达式为:
Φi0(X)=exp{-j[α(Xs+Ci-XP)2]}
其中,i∈{1,2,…,N},N为子孔径个数,α=2π/λRb;展开第i个子孔径数据的相位项Φi0(X),得到第i个子孔径数据的相位项表达式Φ'i0(X):
其中第一个相位项为方位调制项,决定方位向聚焦;
第二项exp{j2αXs(Xp-Ci)}为Xs的线性项,决定方位聚焦位置;
第三项为常数相位项,不影响方位向聚焦。
需要注意的是,在完成方位聚焦后,同一点目标在不同子孔径中将会聚焦在不同的方位位置,Ki为点目标P在第i个子孔径数据中的方位位置,且Ki=2k(XP-Ci),该方位位置与第i个子孔径中心在全孔径坐标系下的位置Ci有关;同时,在常数相位项中,后两项相位均与Ci相关,α=2π/λRb,这种子孔径位置的相关性将会对子孔径数据的相干叠加造成困难。
下面结合图4以及图5a至图5c对步骤8进行具体说明,步骤8的具体子步骤为:
(8.1)对第i个子孔径数据的相位项表达式Φ'i0(X)进行方位匹配滤波处理,得出方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X)。
具体地,将第i个子孔径数据的相位项表达式Φ'i0(X)乘以方位匹配滤波函数Hquad,得出方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X);所述方位匹配滤波函数Hquad的表达式为:
其中,exp表示指数函数,j为虚数单位;所述方位匹配滤波函数Hquad的相位-位置变化率分布线如图5a中的短实线所示,图5a为本发明实施例中距离聚焦处理后,处于不同方位位置的点目标A'、点目标B'和点目标C'各自的相位-位置变化率分布线示意图;所述方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X)的表达式为:
其中,i∈{1,2,…,N},N为子孔径个数,α=2π/λRb,
(8.2)对方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X)进行常数相位补偿,得到常数相位补偿后的第i个子孔径数据的二维位置域信号Φi2(X)。
具体地,将方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X)乘以对应第i个子孔径数据的常数相位补偿函数 补偿第i个子孔径数据中的常数相位项,进而得到常数相位补偿后的第i个子孔径数据Φi2(X),其表达式为:
由图5a可以看出,点目标A'或点目标B'或点目标C'各自的相位-位置变化率分布线在不同子孔径中分别沿纵坐标错开分布,如图5a中虚线所示,并且在纵坐标轴上的投影分别处于不同的位置;对常数相位补偿后的第i个子孔径数据Φi2(X)做方位向快速傅里叶变换(FFT)后,点目标A'或点目标B'或点目标C'在不同的子孔径中将分别出现在各自不同的方位位置处,使得不同子孔径中的同一点目标无法相干叠加。
(8.3)对常数相位补偿后的第i个子孔径数据Φi2(X)进行方位位置补偿,得到方位位置补偿后的第i个子孔径数据Φi3(X)。
具体地,将常数相位补偿后的第i个子孔径数据Φi2(X)乘以对应第i个子孔径数据的方位位置补偿函数Hi,得到方位位置补偿后的第i个子孔径数据Φi3(X);所述对应第i个子孔径数据的方位位置补偿函数Hi的表达式为:
Hi=exp{j2αXsCi}
所述方位位置补偿后的第i个子孔径数据Φi3(X)的表达式为:
其中,exp表示指数函数,j为虚数单位;该操作消除了同一点目标(点目标A'或点目标B'或点目标C')在不同子孔径中相位-位置变化率分布线的错位现象,使得同一点目标(点目标A'或点目标B'或点目标C')在不同子孔径中的方位位置一致;方位位置补偿后的第i个子孔径数据Φi3(X)的相位-位置变化率分布线如图5b所示,可以看出,同一点目标(点目标A'或点目标B'或点目标C')分别在各子孔径中的分布线在纵轴上的投影相重合;方位向快速傅里叶变换(FFT)后,同一点目标(点目标A'或点目标B'或点目标C')在不同子孔径中将会聚焦到相同的位置。
需要注意的是,方位位置补偿后的第i个子孔径数据Φi3(X)的常数相位中存在子孔径位置相关项exp{j2αXpCi},该子孔径位置相关项exp{j2αXpCi}破坏了各子孔径中的同一点目标(点目标A'或点目标B'或点目标C')相位间的相干性;同时,在各子孔径中,同一点目标在方位向分别存在于各自独立的局部坐标系中;同一点目标(点目标A'或点目标B'或点目标C')的空间域支撑区(即点目标A'或点目标B'或点目标C'各自的相位-位置变化率分布线分别在横轴上的投影)在方位空间轴(Xs轴)重合,均为-Ls/2~Ls/2。因此,直接进行子孔径数据的叠加无法使得空间域支撑区得到扩展,不能得到全孔径数据的高分辨成像。
(8.4)对方位位置补偿后的第i个子孔径数据Φi3(X)做方位向快速傅里叶变换(FFT)变换到方位波数域即成像域,进而得到聚焦的子孔径图像,并将聚焦的子孔径图像乘以支撑区全局化函数Hglb,Hglb=exp{-jKxCi},然后再做方位向逆快速傅里叶变换(IFFT),得到第i个子孔径的成像数据Φi4(X)。
所述第i个子孔径的成像数据Φi4(X)的表达式为:
第i个子孔径的成像数据Φi4(X)的相位-位置变化率分布线如图5c所示,可以看出,对聚焦的子孔径图像乘以支撑区全局化函数,能够使得同一点目标(点目标A'或点目标B'或点目标C')从各自的局部坐标系统一到全孔径的全局坐标系,其中点目标A'有三个子孔径,对应分别为第1个子孔径中心在全孔径坐标系下的位置C1、第2个子孔径中心在全孔径坐标系下的位置C2和第3个子孔径中心在全孔径坐标系下的位置C3;点目标B'有三个子孔径,对应分别为第1个子孔径中心在全孔径坐标系下的位置C1、第2个子孔径中心在全孔径坐标系下的位置C2和第3个子孔径中心在全孔径坐标系下的位置C3;点目标C'有三个子孔径,对应分别为第1个子孔径中心在全孔径坐标系下的位置C1、第2个子孔径中心在全孔径坐标系下的位置C2和第3个子孔径中心在全孔径坐标系下的位置C3;同样地,第1个子孔径中心在全孔径坐标系下的位置C1包含三个点目标,分别为点目标A'、点目标B'和点目标C';第2个子孔径中心在全孔径坐标系下的位置C2包含三个点目标,分别为点目标A'、点目标B'和点目标C';第3个子孔径中心在全孔径坐标系下的位置C3包含三个点目标,分别为点目标A'、点目标B'和点目标C'。
并且各子孔径中的同一点目标(点目标A'或点目标B'或点目标C')的空间域支撑区处于各自正确的区间位置,融合为全孔径数据的完整支撑区;与此同时,第i个子孔径的成像数据Φi4(X)中子孔径位置相关的常数项exp{j2αXPCi}得以消除,仅余下一个统一的常数相位项
步骤9,分别令i取1至N,分别得到第1个子孔径的成像数据Φ14(X)至第N个子孔径的成像数据ΦN4(X),然后对第1个子孔径的成像数据Φ14(X)至第N个子孔径的成像数据ΦN4(X)进行叠加,得到全孔径雷达成像数据Φ4(X),最后对全孔径雷达成像数据Φ4(X)进行方位向快速傅里叶变换(FFT),进而得到方位向聚焦的全孔径高分辨率SAR成像。
下面通过点目标的仿真成像实验进一步说明本发明方法的正确性和有效性。
(一)点目标仿真成像仿真条件
表1雷达参数
(二)仿真内容
为验证本发明的性能和有效性,此处取4096方位点数作为处理数据,子孔径点数设为512点,子孔径数为8;使用本发明方法对各子孔径数据进行成像及相干性处理,得到子图像后进行方位向空间域支撑区的融合,形成高分辨率的SAR图像。
3)仿真结果分析
参照图6a,为仿真实验中使用本发明得出的单个点目标的二维升采样后的结果图,横轴表示方位向,纵轴表示距离向,这里选取了单子孔径结果;图6b为其相对应的方位向剖面图,横轴表示方位采样单元,纵轴表示各点的归一化幅度,单位为dB,3dB分辨率如图左上角所示。
参照图7a,为仿真实验中使用本发明得出的单个点目标的二维升采样后的结果图,横轴表示方位向,纵轴表示距离向,这里选取了两子孔径融合结果;图7b为其相对应的方位向剖面图,横轴表示方位采样单元,纵轴表示各点的归一化幅度,单位为dB,3dB分辨率如图左上角所示。
参照图8a,为仿真实验中使用本发明得出的单个点目标的二维升采样后的结果图,横轴表示方位向,纵轴表示距离向,这里选取了四子孔径融合结果;图8b为其相对应的方位向剖面图,横轴表示方位采样单元,纵轴表示各点的归一化幅度,单位为dB,3dB分辨率如图左上角所示。
参照图9a,为仿真实验中使用本发明得出的单个点目标的二维升采样后的结果图,横轴表示方位向,纵轴表示距离向,这里选取了八子孔径融合结果;图9b为其相对应的方位向剖面图,横轴表示方位采样单元,纵轴表示各点的归一化幅度,单位为dB,3dB分辨率如图左上角所示。
从图6a至图9b可以看出,随着子孔径数据的不断融合叠加,方位向聚焦明显改善,方位向分辨率成倍增长。当为八子孔径融合时,其成像质量达到了实际4096点数据一次性处理时的分辨率(0.27m)。该实验验证了方位向空间域支撑区融合的有效性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种基于空间域支撑区融合的SAR成像方法,其特征在于,包括以下步骤:
步骤1,确定机载合成孔径雷达,所述机载合成孔径雷达发射信号并接收回波信号,并将回波信号记为子孔径回波数据;
步骤2,对子孔径回波数据进行2-维快速傅里叶变换,得到2-维波数域子孔径数据;
步骤3,对2-维波数域子孔径数据进行距离匹配滤波处理,得到距离匹配滤波后的2-维波数域子孔径数据;
步骤4,对距离匹配滤波后的2-维波数域子孔径数据进行相位补偿处理,得到相位补偿处理后的2-维波数域子孔径数据;
步骤5,对相位补偿处理后的2-维波数域子孔径数据进行扩展stolt插值处理,得到扩展stolt插值处理后的2-维波数域子孔径数据;
步骤6,对扩展stolt插值处理后的2-维波数域子孔径数据进行距离向脉压处理,即进行距离向逆快速傅里叶变换,进而计算得到距离向脉压处理后的距离位置域-方位波数域子孔径数据;
步骤7,对距离向脉压处理后的距离位置域-方位波数域子孔径数据进行去斜处理,得到去斜处理后的距离位置域-方位波数域子孔径数据,然后对去斜处理后的距离位置域-方位波数域子孔径数据进行方位向逆快速傅里叶变换,得到2-维位置域子孔径数据;
步骤8,初始化:令i∈{1,2,…,N},N为子孔径个数,i的初始值为1;
对第i个子孔径,取2-维位置域子孔径数据的相位项,记为第i个子孔径数据的相位项,然后对第i个子孔径数据的相位项进行相干叠加处理及方位向空间域支撑区融合,即对第i个子孔径数据的相位项依次进行方位匹配滤波处理、常数相位补偿、方位位置补偿、方位向快速傅里叶变换、全局化处理、方位向逆快速傅里叶变换、相干叠加、方位波数域聚焦成像,进而得到第i个子孔径的成像数据;
步骤9,分别令i取1至N,分别得到第1个子孔径的成像数据至第N个子孔径的成像数据,然后对第1个子孔径的成像数据至第N个子孔径的成像数据进行叠加,得到全孔径雷达成像数据,最后对全孔径雷达成像数据进行方位向快速傅里叶变换,进而得到方位向聚焦的全孔径高分辨率SAR成像。
2.如权利要求1所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤1中,所述子孔径回波数据,其得到过程为:
建立三维坐标系XOYZ,其中机载合成孔径雷达工作于条带模式,机载合成孔径雷达的载机以水平速度v沿机载合成孔径雷达的方位向飞行,机载合成孔径雷达的载机飞行高度为H,选取机载合成孔径雷达照射场景中的任意一个散射点,记为点目标P;点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离为Rb,机载合成孔径雷达照射检测范围内场景的中心点为点B,点目标P沿方位向与点B之间的距离为XP;机载合成孔径雷达以重复周期T、发射时刻tm依次发射信号,tm=mT,m为正整数变量;
当经过任意时间后,由几何关系计算得到机载合成孔径雷达与点目标P之间的瞬时斜距R(X;Rb)为:
其中,X为机载合成孔径雷达的方位向位置变量,X=vtm,tm=mT,T为机载合成孔径雷达发射信号的重复周期,m为正整数变量;
机载合成孔径雷达发射信号并接收回波信号,然后对接收到的回波信号进行相干检波,即乘以exp(-j2πfct),fc为机载合成孔径雷达发射信号的载波频率,t为时间变量;进而得到机载合成孔径雷达接收的基带回波信号,所述机载合成孔径雷达接收的基带回波信号为子孔径回波数据S0(Rr,X;Rb),其表达式为:
其中,Rr为机载合成孔径雷达的距离向位置变量,C为光速,t为时间变量;R(X;Rb)为机载合成孔径雷达与点目标P之间的瞬时斜距,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,X为机载合成孔径雷达的方位向位置变量,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,γ为机载合成孔径雷达发射信号的调频率,λ为机载合成孔径雷达发射信号的载波波长,wr(Rr)为机载合成孔径雷达发射信号的距离窗函数,wa(X)为机载合成孔径雷达发射信号的方位窗函数,exp表示指数函数,j为虚数单位。
3.如权利要求1所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤2中,所述2-维波数域子孔径数据为S1(Kr,Kx),其表达式为:
其中,Kr为机载合成孔径雷达的距离波数变量,Kr=4πf/C,f表示机载合成孔径雷达发射信号的频率在点目标P上的频率分量,f=fc+fb,fb为机载合成孔径雷达发射信号的频率在点目标P上的基频,
Kr=Krc+△Kr,Krc为机载合成孔径雷达的距离波数中心值,
Krc=4πfc/C,fc为机载合成孔径雷达发射信号的载波频率,△Kr为差值,△Kr∈[-2πγTp/C,2πγTp/C],C为光速,γ为机载合成孔径雷达发射信号的调频率,Tp为机载合成孔径雷达发射信号的脉冲宽度,∈表示属于;W(Kr)为机载合成孔径雷达发射信号的距离窗函数的波数域形式,Kx为机载合成孔径雷达的方位波数变量,Kx=2πfd/v,v为机载合成孔径雷达的载机水平速度,fd为机载合成孔径雷达接收到的回波信号的多普勒频率;W(Kx)为机载合成孔径雷达发射信号的方位窗函数的波数域形式,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,点B为机载合成孔径雷达照射检测范围内场景的中心点,X为机载合成孔径雷达的方位向位置变量,exp表示指数函数,j为虚数单位。
4.如权利要求1所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤3中,所述距离匹配滤波后的2-维波数域子孔径数据,其得到过程为:
确定距离匹配滤波函数HRMF(Kr),其表达式为:
将2-维波数域子孔径数据S1(Kr,Kx)乘以距离匹配滤波函数HRMF(Kr),消除距离二次项,得到距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx);
所述距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx)的表达式为:
其中,W(Kr)为机载合成孔径雷达发射信号的距离窗函数的波数域形式,W(Kx)为机载合成孔径雷达发射信号的方位窗函数的波数域形式,Kx为机载合成孔径雷达的方位波数变量,Kr为机载合成孔径雷达的距离波数变量,Krc为机载合成孔径雷达的距离波数中心值,C为光速,γ为机载合成孔径雷达发射信号的调频率,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,点B为机载合成孔径雷达照射检测范围内场景的中心点,exp表示指数函数,j为虚数单位。
5.如权利要求1所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤4中,所述相位补偿处理后的2-维波数域子孔径数据,其得到过程为:
确定相位补偿函数HBulk(Kr,Kx),其表达式为:
其中,Rs表示机载合成孔径雷达照射检测范围内场景的中心点B对应的参考距离,Kx为机载合成孔径雷达的方位波数变量,Kr为机载合成孔径雷达的距离波数变量,exp表示指数函数,j为虚数单位;
将距离匹配滤波后的2-维波数域子孔径数据S2(Kr,Kx)乘以相位补偿函数HBulk(Kr,Kx),得到相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx);
所述相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)的表达式为:
其中,W(Kr)为机载合成孔径雷达发射信号的距离窗函数的波数域形式,W(Kx)为机载合成孔径雷达发射信号的方位窗函数的波数域形式,Kx为机载合成孔径雷达的方位波数变量,Kr为机载合成孔径雷达的距离波数变量,Rs表示机载合成孔径雷达照射检测范围内场景的中心点B对应的参考距离,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离,点目标P为机载合成孔径雷达照射场景中的任意一个散射点,点B为机载合成孔径雷达照射检测范围内场景的中心点,exp表示指数函数,j为虚数单位。
6.如权利要求5所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤5中,所述扩展stolt插值处理后的2-维波数域子孔径数据,其得到过程为:
对相位补偿处理后的2-维波数域子孔径数据S3(Kr,Kx)进行扩展stolt插值处理,其中扩展stolt插值处理的映射表达式为
进而计算得到扩展stolt插值处理后的2-维波数域子孔径数据S4(Ky,Kx),其表达式为:
S4(Ky,Kx)=W(Ky)W(Kx)exp(-jKxXP)
其中,Ky表示stolt插值处理时引入的机载合成孔径雷达的距离波数变量,Krc为机载合成孔径雷达的距离波数中心值。
7.如权利要求5所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤6中,所述距离向脉压处理后的距离位置域-方位波数域子孔径数据为S5(Rr,Kx),其表达式为:
S5(Rr,Kx)=W(Kx)sinc{A[Rr-(Rb-Rs)]}
其中,A为距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)的常数项,W(Kx)表示机载合成孔径雷达发射信号的方位窗函数的波数域形式,exp表示指数函数,j为虚数单位,
a=A[Rr-(Rb-Rs)];Rs表示机载合成孔径雷达照射检测范围内场景的中心点B对应的参考距离;Rr为机载合成孔径雷达的距离向位置变量;Krc为机载合成孔径雷达的距离波数中心值。
8.如权利要求5所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤7中,所述2-维位置域子孔径数据,其得到过程为:
确定去斜校正函数H3(Kx),
然后,将距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)乘以去斜校正函数H3(Kx),进而得到去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx);所述去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx)的表达式为:
其中,a=A[Rr-(Rb-Rs)],A为距离向脉压处理后的距离位置域-方位波数域子孔径数据S5(Rr,Kx)的常数项;Rr为机载合成孔径雷达的距离向位置变量;Krc为机载合成孔径雷达的距离波数中心值;
对去斜处理后的距离位置域-方位波数域子孔径数据S6(Rr,Kx)进行方位向逆快速傅里叶变换处理,得到2-维位置域子孔径数据S7(Rr,X);所述2-维位置域子孔径数据S7(Rr,X)的表达式为:
其中,wa(X)为机载合成孔径雷达发射信号的方位窗函数;X为机载合成孔径雷达的方位向位置变量,X=vtm,tm=mT,T为机载合成孔径雷达发射信号的重复周期,m为正整数变量;λ为机载合成孔径雷达发射信号的载波波长;v为机载合成孔径雷达的载机以沿机载合成孔径雷达的方位向飞行的水平速度。
9.如权利要求1或5所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤8中,所述N为子孔径个数,其得到过程为:
在三维坐标系XOYZ中机载合成孔径雷达载机从左向右飞行,且机载合成孔径雷达以重复周期T、发射时刻为tm依次发射信号,tm=mT,m为正整数变量;
第m个发射时刻tm时机载合成孔径雷达发射的信号照射到点目标P并产生回波数据,记为第m个子孔径,且第m个子孔径对应的回波数据为第m个子孔径数据;
令m分别取1至N,进而分别得到第1个子孔径至第N个子孔径,以及第1个子孔径数据至第N个子孔径数据;然后分别将第1个子孔径至第N个子孔径记为全孔径,将第1个子孔径数据至第N个子孔径数据记为全孔径数据;
建立平面直角坐标系,其横轴为机载合成孔径雷达的方位向位置变量X,纵轴为相位-位置变化率 λ为机载合成孔径雷达发射信号的载波波长,Rb为点目标P与机载合成孔径雷达的载机飞行轨迹之间的最近距离,XP为点目标P沿方位向与点B之间的距离;并将该平面直角坐标系作为全孔径坐标系,其中全孔径坐标系中对应包含N个子孔径坐标系,将第i个子孔径中心在全孔径坐标系下的位置记为Ci,i∈{1,2,…,N},N为子孔径个数。
10.如权利要求9所述的一种基于空间域支撑区融合的SAR成像方法,其特征在于,在步骤8中,所述得到第i个子孔径的成像数据,其过程为:
对第i个子孔径,取2-维位置域子孔径数据S7(Rr,X)的相位项,记为第i个子孔径数据的相位项Φi0(X),其表达式为:
Φi0(X)=exp{-j[α(Xs+Ci-XP)2]}
其中,i∈{1,2,…,N},N为子孔径个数,α=2π/λRb;展开第i个子孔径数据的相位项Φi0(X),得到第i个子孔径数据的相位项表达式Φ'i0(X):Xs为每个子孔径方位向采样位置,Xs∈[-Ls/2,Ls/2];
(8.1)对第i个子孔径数据的相位项表达式Φ'i0(X)进行方位匹配滤波处理,得出方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X);
将第i个子孔径数据的相位项表达式Φ'i0(X)乘以方位匹配滤波函数Hquad,得出方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X);所述方位匹配滤波函数Hquad的表达式为:
其中,exp表示指数函数,j为虚数单位;所述方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X)的表达式为:
(8.2)对方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X)进行常数相位补偿,得到常数相位补偿后的第i个子孔径数据的二维位置域信号Φi2(X);
将方位匹配滤波处理后的第i个子孔径数据的二维位置域信号Φi1(X)乘以对应第i个子孔径数据的常数相位补偿函数 补偿第i个子孔径数据中的常数相位项,进而得到常数相位补偿后的第i个子孔径数据的二维位置域信号Φi2(X),其表达式为:
(8.3)对常数相位补偿后的第i个子孔径数据的二维位置域信号Φi2(X)进行方位位置补偿,得到方位位置补偿后的第i个子孔径数据的二维位置域信号Φi3(X);
将常数相位补偿后的第i个子孔径数据的二维位置域信号Φi2(X)乘以对应第i个子孔径数据的方位位置补偿函数Hi,得到方位位置补偿后的第i个子孔径数据的二维位置域信号Φi3(X);所述对应第i个子孔径数据的方位位置补偿函数Hi的表达式为:
Hi=exp{j2αXsCi}
所述方位位置补偿后的第i个子孔径数据的二维位置域信号Φi3(X)的表达式为:
其中,exp表示指数函数,j为虚数单位;
(8.4)对方位位置补偿后的第i个子孔径数据的二维位置域信号Φi3(X)做方位向快速傅里叶变换变换到方位波数域即成像域,进而得到聚焦的子孔径图像,并将聚焦的子孔径图像乘以支撑区全局化函数Hglb,Hglb=exp{-jKxCi},其中,Kx为机载合成孔径雷达的方位波数变量;然后再做方位向逆快速傅里叶变换,得到第i个子孔径的成像数据Φi4(X);所述第i个子孔径的成像数据Φi4(X)的表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111407.5A CN106908789B (zh) | 2017-02-28 | 2017-02-28 | 一种基于空间域支撑区融合的sar成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710111407.5A CN106908789B (zh) | 2017-02-28 | 2017-02-28 | 一种基于空间域支撑区融合的sar成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106908789A CN106908789A (zh) | 2017-06-30 |
CN106908789B true CN106908789B (zh) | 2019-08-02 |
Family
ID=59208610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710111407.5A Active CN106908789B (zh) | 2017-02-28 | 2017-02-28 | 一种基于空间域支撑区融合的sar成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106908789B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110515080B (zh) * | 2019-08-30 | 2022-12-02 | 西安电子科技大学 | 基于半径角度插值的高超声速机动平台sar成像方法 |
CN112395546B (zh) * | 2020-11-27 | 2022-07-01 | 北京理工大学 | 一种基于线性正则域的有限新息率信号降采样与重构方法 |
CN112946649B (zh) * | 2021-04-08 | 2022-08-26 | 电子科技大学 | 一种适用于任意子孔径长度的pfa成像方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6670907B2 (en) * | 2002-01-30 | 2003-12-30 | Raytheon Company | Efficient phase correction scheme for range migration algorithm |
CN101833095A (zh) * | 2010-04-14 | 2010-09-15 | 电子科技大学 | 一种基于空域展开的星机联合sar二维频域成像方法 |
CN104597447A (zh) * | 2015-01-30 | 2015-05-06 | 西安电子科技大学 | 一种子孔径SAR大斜视改进Omega-K成像方法 |
CN104865571A (zh) * | 2015-04-02 | 2015-08-26 | 西安电子科技大学 | 多通道多子带滑动聚束模式sar成像方法 |
CN105652273A (zh) * | 2016-03-17 | 2016-06-08 | 哈尔滨工程大学 | 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法 |
CN105974409A (zh) * | 2016-01-13 | 2016-09-28 | 合肥工业大学 | 基于多频子带并发的星载滑动聚束mimo-sar成像方法 |
CN106054152A (zh) * | 2016-05-23 | 2016-10-26 | 西安电子科技大学 | 基于逆扩展Omega‑K算法的非理想轨迹SAR回波获取方法 |
-
2017
- 2017-02-28 CN CN201710111407.5A patent/CN106908789B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6670907B2 (en) * | 2002-01-30 | 2003-12-30 | Raytheon Company | Efficient phase correction scheme for range migration algorithm |
CN101833095A (zh) * | 2010-04-14 | 2010-09-15 | 电子科技大学 | 一种基于空域展开的星机联合sar二维频域成像方法 |
CN104597447A (zh) * | 2015-01-30 | 2015-05-06 | 西安电子科技大学 | 一种子孔径SAR大斜视改进Omega-K成像方法 |
CN104865571A (zh) * | 2015-04-02 | 2015-08-26 | 西安电子科技大学 | 多通道多子带滑动聚束模式sar成像方法 |
CN105974409A (zh) * | 2016-01-13 | 2016-09-28 | 合肥工业大学 | 基于多频子带并发的星载滑动聚束mimo-sar成像方法 |
CN105652273A (zh) * | 2016-03-17 | 2016-06-08 | 哈尔滨工程大学 | 一种基于混合匹配追踪算法的mimo雷达稀疏成像算法 |
CN106054152A (zh) * | 2016-05-23 | 2016-10-26 | 西安电子科技大学 | 基于逆扩展Omega‑K算法的非理想轨迹SAR回波获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106908789A (zh) | 2017-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101430380B (zh) | 基于非均匀采样的大斜视角机载sar聚束模式成像方法 | |
CN105974414B (zh) | 基于二维自聚焦的高分辨聚束sar自聚焦成像方法 | |
Zhang et al. | Spaceborne/stationary bistatic SAR imaging with TerraSAR-X as an illuminator in staring-spotlight mode | |
CN110146889B (zh) | 基于最优回波子区域选择的大转角isar成像方法 | |
Fan et al. | A high-precision method of phase-derived velocity measurement and its application in motion compensation of ISAR imaging | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN104898119B (zh) | 一种基于相关函数的动目标参数估计方法 | |
CN108490443B (zh) | 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法 | |
CN104597447B (zh) | 一种子孔径SAR大斜视改进Omega‑K成像方法 | |
CN106908789B (zh) | 一种基于空间域支撑区融合的sar成像方法 | |
CN109143237A (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN108152820B (zh) | 一种基于层析原理的双基地雷达成像方法 | |
CN112859018B (zh) | 一种基于图像几何校正的视频sar成像方法 | |
CN109444882B (zh) | 基于变斜视椭圆波束同步模型的双站sar成像方法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
CN103576130A (zh) | 一种进动锥体的三维成像方法 | |
Pailhas et al. | Impact of temporal Doppler on synthetic aperture sonar imagery | |
CN101620272A (zh) | 一种逆合成孔径雷达的目标转速估计方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN108594196B (zh) | 一种目标散射中心参数提取方法及装置 | |
CN114325704A (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN110632616B (zh) | 一种稀疏采样下机载逆合成孔径激光雷达微动成像方法 | |
CN115453530B (zh) | 一种基于参数化模型的双基sar滤波反投影两维自聚焦方法 | |
Khwaja et al. | SAR raw data generation using inverse SAR image formation algorithms | |
Saeedi et al. | Improved navigation-based motion compensation for LFMCW synthetic aperture radar imaging |
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 |