CN104597447B - 一种子孔径SAR大斜视改进Omega‑K成像方法 - Google Patents
一种子孔径SAR大斜视改进Omega‑K成像方法 Download PDFInfo
- Publication number
- CN104597447B CN104597447B CN201510051776.0A CN201510051776A CN104597447B CN 104597447 B CN104597447 B CN 104597447B CN 201510051776 A CN201510051776 A CN 201510051776A CN 104597447 B CN104597447 B CN 104597447B
- Authority
- CN
- China
- Prior art keywords
- prime
- signal
- azimuth
- distance
- theta
- 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 80
- 238000012937 correction Methods 0.000 claims abstract description 43
- 238000001228 spectrum Methods 0.000 claims abstract description 20
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 238000013507 mapping Methods 0.000 claims abstract description 3
- 238000012545 processing Methods 0.000 claims description 42
- 208000004350 Strabismus Diseases 0.000 claims description 39
- 230000035485 pulse pressure Effects 0.000 claims description 27
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 239000013256 coordination polymer Substances 0.000 claims description 4
- 238000010408 sweeping Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 abstract description 5
- 238000012952 Resampling Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 54
- 238000004088 simulation Methods 0.000 description 43
- 230000006870 function Effects 0.000 description 36
- 238000004422 calculation algorithm Methods 0.000 description 35
- 230000000875 corresponding effect Effects 0.000 description 10
- 230000008859 change Effects 0.000 description 8
- 230000004044 response Effects 0.000 description 4
- 230000001629 suppression Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000005526 G1 to G0 transition Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000012292 cell migration Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000001131 transforming effect Effects 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
- G01S13/904—SAR modes
-
- 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/9041—Squint 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成像技术领域,特别涉及一种子孔径SAR大斜视改进Omega‑K成像方法,其具体步骤为:(1)对原始回波数据进行距离向傅里叶变换;(2)对信号进行距离匹配滤波处理,并乘以旋转校正函数,实现波数谱支撑区的正侧化处理;(3)对信号进行方位向傅里叶变换,在二维波数域进行方位重采样插值;(4)对信号乘以统一相位补偿函数,并进行扩展Stolt插值,然后通过距离逆傅里叶变换实现距离向的空间位置域聚焦;(5)对信号乘以去斜校正函数,并进行方位逆傅里叶变换,在方位空间位置域乘以方位位置校正函数,最后再通过方位傅里叶变换,实现方位向的波数域聚焦。本发明应用范围广泛,可实现较高分辨率成像,可用于地图测绘,目标识别等领域。
Description
技术领域
本发明属于SAR成像技术领域,特别涉及一种子孔径SAR大斜视改进Omega-K成像方法,可用于机载、星载平台的SAR成像处理。
背景技术
SAR(合成孔径雷达)作为一种主动探测工具能够对观测场景进行高分辨微波成像,在遥感领域得到广泛应用。其通过发射宽频带的信号获得距离维的高分辨特性,而方位维则利用平台运动产生的多普勒信息实现高分辨。
近年来随着各种成像算法的成熟和发展,兼顾成像质量的同时,实时性成为很多成像系统的必备要求。子孔径概念的提出一定程度解决了成像质量与成像效率这对矛盾,由于子孔径数据成像具有录取时间短、积累的相位误差小等特点,可在适当损失分辨率的代价下简化处理流程,减小了运动补偿复杂度、计算量和存储量,实现实时成像。
而对于大斜视(斜视角大于45度)成像,最大的困难在于解决由大的斜视角所带来的距离单元徙动以及距离方位向的二维耦合问题。传统的斜视SAR成像方法大都是通过一定的近似消除这种耦合,包括距离多普勒算法(RDA),调频变标算法(CSA),SPECAN算法等,但这种近似在一定程度上使其使用范围得到限制。而Omega-K算法则可以通过Stolt插值实现无近似的RCM校正解耦合,是一种较为理想的成像方法。
现有Omega-K算法主要针对于全孔径成像的研究,并且忽略大斜视角所带来的波数谱支撑区的斜拉特性,较大的影响了成像质量。另外,传统的Omega-K算法直接应用于子孔径数据则会导致方位位置的混叠错位,无法反应真实方位位置,同时无法进行统一方位的加窗抑制旁瓣。
发明内容
本发明的目的在于提出一种子孔径SAR大斜视改进Omega-K成像方法,本发明将斜视二维波数谱正侧化为正侧视情况,增大波数域支撑区的利用率,并通过方位重采样解决方位调频率的空变,进而改善在大斜视情况下的成像质量。同时该发明能够适用于子孔径数据成像,解决了由于子孔径方位位置支撑区较小而出现的方位混叠问题,并实现了方位的统一加窗处理。
为实现上述技术目的,本发明采用如下技术方案予以实现。
一种子孔径SAR大斜视改进Omega-K成像方法包括以下步骤:
步骤1,利用机载合成孔径雷达发射线性调频信号,利用机载合成孔径雷达接收对应的回波信号;对回波信号进行解调,得到解调后的基带回波信号其中,为机载合成孔径雷达距离向位置变量,X为机载合成孔径雷达方位向位置变量,R0为机载合成孔径雷达波束中心扫过场景中心点时的斜距;对基带回波信号进行距离向傅里叶变换,得出距离波数域方位空间位置域信号S1(Kr,X),Kr为机载合成孔径雷达距离波数变量;
步骤2,将距离波数域方位空间位置域信号S1(Kr,X)进行距离向匹配滤波处理,得出距离匹配滤波后信号S2(Kr,X);
步骤3,将距离匹配滤波后信号S2(Kr,X)乘以旋转校正函数HLRWC(Kr,X),得出旋转校正后信号S3(Kr,X),其中,旋转校正函数HLRWC(Kr,X)=exp[-jKrXsinθ0],θ0为机载合成孔径雷达的波束中心斜视角;
步骤4,对旋转校正后信号S3(Kr,X)进行方位向傅里叶变换,得到二维波数域信号S4(Kr,Kx),Kx为机载合成孔径雷达方位波数变量;
步骤5,将R0=R0'-Xnsinθ0代入二维波数域信号S4(Kr,Kx)的表达式中,得出第一次形式变换后的二维波数域信号S5(Kr,Kx);令
将Kx'的表达式代入第一次形式变换后的二维波数域信号S5(Kr,Kx)的表达式中,得出第二次形式变换后的二维波数域信号S6(Kr,Kx');
步骤6,对二维波数域信号S6(Kr,Kx')依次进行相位补偿、Stolt插值处理、距离向逆傅里叶变换、去斜处理、方位向逆傅里叶变换、方位波数域聚焦成像,得出最终SAR聚焦成像结果。
本发明的有益效果为:
1)本发明采用Omega-K算法的Stolt插值实现RCM和耦合量的无近似校正,改善了成像效果。
2)传统Omega-K算法在应用时忽略了大斜视角所带来的二维波数谱的斜拉特点,在完成插值后进行二维逆傅里叶变换时,选取矩形区支撑区受限,导致利用率较低,影响成像质量。本发明通过波数谱的正侧化,实现了波数谱支撑区的“扳正”处理,进而扩大了支撑区的利用率,大幅改善成像分辨率。
3)对于子孔径数据成像,由于受限制于较小的子孔径位置支撑区,传统Omega-K算法在方位位置域聚焦将可能导方位聚焦位置出现混叠甚至反折,无法真实反映目标几何位置,需要进行方位补零扩展,大大加剧算法运算量。本发明通过改进Omega-K算法,将距离向与方位向分开处理,实现距离向空间位置域聚焦,而方位向实现波数域无混叠聚焦,避免了大量的补零操作,降低了运算量。
4)由于子孔径数据成像时,同一距离单元内不同方位位置的点目标,在完成距离压缩后,其方位向谱线处于同一距离单元但在方位向各自分离,无法进行统一的加窗处理。本发明通过“去斜”处理,通过在二维波数域乘以校正函数实现了各点目标数据在方位位置域的对齐,进而可通过统一加窗进行旁瓣抑制处理。
附图说明
图1a为斜视SAR成像几何模型示意图;图1b为机与场景散射点斜平面几何关系示意图;
图2为本发明的一种子孔径SAR大斜视改进Omega-K成像方法的流程图;
图3a为大斜视SAR的二维斜视成像点目标二维波数谱支撑区示意图,图3b为二维波数谱支撑区被扳正后的点目标二维波数谱支撑区示意图;
图4为本发明实施例中时域校正线性走动示意图;
图5a是本发明实施例中距离聚焦处理后相位-波数变化率分布线示意图;图5b是本发明实施例中统一“去斜”校正函数的相位-波数变化率分布线示意图;图5c为本发明实施例中“去斜”处理后相位-波数变化率分布线示意图;图5d为本发明实施例中“去斜”处理后进行方位IFFT后的相位-位置变换率分布线示意图;图5e为本发明实施例中方位位置域校正函数的相位-位置变换率分布线示意图;图5f为本发明实施例中经过方位位置域校正函数扳正后的相位-位置变换率分布线示意图;图5g为本发明实施例中最终经方位FFT在方位波数域聚焦时的相位-波数变化率分布线示意图;
图6为仿真实验中点目标仿真几何示意图;
图7a为仿真实验1中标准Omega-K算法的最终SAR聚焦成像结果示意图,图7b为仿真实验1中本发明得出的最终SAR聚焦成像结果示意图,图7c,为仿真实验1中利用本发明进行SAR聚焦成像时的成像结果插值示意图;
图8a为仿真实验2中点A通过本发明进行成像时方位脉压剖面图在加窗处理前的结果示意图,图8b为仿真实验2中点B通过本发明进行成像时方位脉压剖面图在加窗处理前的结果示意图,图8c为仿真实验2中点C通过本发明进行成像时方位脉压剖面图在加窗处理前的结果示意图;图8d为仿真实验2中点A通过本发明进行成像时方位脉压剖面图在加窗处理后的结果示意图,图8e为仿真实验2中点B通过本发明进行成像时方位脉压剖面图在加窗处理后的结果示意图,图8f为仿真实验2中点C通过本发明进行成像时方位脉压剖面图在加窗处理后的结果示意图;
图9a为仿真实验2中点A通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理前的结果示意图,图9b为仿真实验2中点B通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理前的结果示意图,图9c为仿真实验2中点C通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理前的结果示意图;图9d为仿真实验2中点A通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理后的结果示意图,图9e为仿真实验2中点B通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理后的结果示意图,图9f为仿真实验2中点C通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理后的结果示意图。
具体实施方式
下面结合附图对本发明作进一步说明:
参照图1a,为斜视SAR成像几何模型示意图;机载合成孔径雷达工作于条带模式,载机以水平速度v沿X方向(雷达的方位向)飞行。θ0为波束中心斜视角,θBW为方位向波束宽度,H为飞行高度,目标与载机飞行轨迹的最近距离为RB,R0为波束中心扫过目标时的斜距,点B为场景中心点,Xn为任一场景散射点P点沿方位向与点B的距离。参照图1b,为载机与场景散射点斜平面几何关系示意图,图1b中,A点为机载合成孔径雷达运动的方位零时刻点,此时机载合成孔径雷达的方位位置为零(X=0),波束中心指向点B,即场景中心点。P点为沿方位向与点B间距为Xn的目标,当经过任意时间后,雷达由A点运动到A'点,此时雷达与P点的瞬时斜距为A'P,瞬时斜视角为θ。由几何关系可知,点P的瞬时斜距R(X;R0)为:
其中X为方位向位置变量,这里我们使用空间位置域斜距表达式以与后续波数域表述相对应。由瞬时斜距表达式,可以容易得到目标的回波信号表达式,本发明中,机载合成孔径雷达发射信号为线性调频(LFM)信号,则经过解调后的基带回波信号可写为:
其中距离向即波束指向方向的位置变量,C为光速,γ为线性调频信号的调频率,λ为机载合成孔径雷达发射信号的载波波长,和wa(X)分别为线性调频信号的距离窗函数和方位窗函数。
参照图2,为本发明的一种子孔径SAR大斜视改进Omega-K成像方法的流程图。该子孔径SAR大斜视改进Omega-K成像方法的具体步骤如下:
步骤1,利用机载合成孔径雷达发射线性调频信号,利用机载合成孔径雷达接收对应的回波信号;对回波信号进行解调,得到解调后的基带回波信号其中,为机载合成孔径雷达距离向位置变量,X为机载合成孔径雷达方位向位置变量,R0为机载合成孔径雷达波束中心扫过场景中心点时的斜距;
采用驻定相位原理,将基带回波信号进行距离向傅里叶变换,得出距离波数域方位空间位置域信号S1(Kr,X),Kr为机载合成孔径雷达距离波数。
在步骤1中,距离波数域方位空间位置域信号S1(Kr,X)的具体表达式为:
其中Kr为机载合成孔径雷达距离波数,Wr(Kr)为距离窗函数的波数域形式,wa(X)表示方位窗函数,C为光速,Krc=4πfc/C,fc为线性调频信号的载波频率,Kr=Krc+ΔKr,ΔKr∈[-2πγTp/C,2πγTp/C],Tp为机载合成孔径雷达发射信号的脉冲宽度,γ为线性调频信号的调频率,R(X;R0)表示任一场景散射点P的瞬时斜距。
步骤2,将距离波数域方位空间位置域信号S1(Kr,X)进行距离匹配滤波处理,得出距离匹配滤波后信号S2(Kr,X)。
具体地说,将距离波数域方位空间位置域信号S1(Kr,X)乘以距离匹配滤波函数,消除距离二次项,得到距离匹配滤波后信号S2(Kr,X);
距离匹配滤波后信号S2(Kr,X)的表达式为:
S2(Kr,X)=S1(Kr,X)·HRMF(Kr)=Wr(Kr)wa(X)exp[-jKrR(X;R0)]
其中,HRMF(Kr)表示距离匹配滤波函数,其具体表达式为:
步骤3,将距离匹配滤波后信号S2(Kr,X)乘以旋转校正函数HLRWC(Kr,X),得出旋转校正后信号S3(Kr,X),其中,旋转校正函数HLRWC(Kr,X)=exp[-jKrX sinθ0],θ0为波束中心斜视角。
具体地说,大斜视SAR的二维波数谱的支撑区呈现斜拉特点。参照图3a,为大斜视SAR的二维斜视成像点目标二维波数谱支撑区示意图;图3a中,Kr为机载合成孔径雷达距离波数,Kx为机载合成孔径雷达方位波数,Krc=4πfc/C,当Kr=Krc时,Kxc=Kxc0。
此时需要将图3a的点目标二维波数谱正侧化以扩大支撑区可用面积。引入旋转校正函数HLRWC(Kr,X)与S2(Kr,X)相乘,得出旋转校正后信号S3(Kr,X),S3(Kr,X)=S2(Kr,X)·HLRWC(Kr,X)=Wr(Kr)wa(X)exp{-jKr[R(X;R0)+X sinθ0]}其中,HLRWC(Kr,X)表示旋转校正函数,HLRWC(Kr,X)=exp[-jKrX sinθ0],Wr(Kr)为距离窗的波数域形式,其自变量为Kr;wa(X)为方位窗函数,θ0为波束中心斜视角,R0为目标波束中心距离,R(X;R0)为瞬时斜距。
步骤4,对旋转校正后信号S3(Kr,X)进行方位向傅里叶变换,将其变换到二维波数域,得到二维波数域信号S4(Kr,Kx)。
具体地,在步骤4中,二维波数域信号S4(Kr,Kx)的表达式为:
S4(Kr,Kx)=Wr(Kr)Wa(Kx)exp{D}
其中,Wa(Kx)为方位窗函数的波数域形式,其自变量为Kx;Xn表示场景散射点相对场景中心点的方位位置,θ0为波束中心斜视角,R0为机载合成孔径雷达波束中心扫过场景中心点时的斜距(即目标波束中心距离)。此时其二维波数谱支撑区被扳正,参照图3b,为二维波数谱支撑区被扳正后的点目标二维波数谱支撑区示意图。
步骤5,将R0=R0'-Xnsinθ0代入二维波数域信号S4(Kr,Kx)的表达式中,得出第一次形式变换后的二维波数域信号S5(Kr,Kx);令
将Kx'的表达式代入第一次形式变换后的二维波数域信号S5(Kr,Kx)的表达式中,得出第二次形式变换后的二维波数域信号S6(Kr,Kx')。
具体地说,二维波数谱支撑区被扳正后,目标波束中心距离由原来的R0变为R0',R0'=R0+Xnsinθ0;参照图4,为本发明实施例中时域校正线性走动示意图;其中,横轴X为方位向位置变量,纵轴R为距离单元对应的最短距离,三条平行的实线为走动校正前的空间响应曲线,斜实线为距离走动校正线,虚线为走动校正后的空间响应曲线。对于点目标B,其方位位置坐标为0;对于点目标A,方位位置坐标X-n为负值;对于点目标C,方位位置坐标Xn为正值。
在图4中,以A、B、C三点分析为例说明:假设三点位于同一距离单元,其中B点位于场景中心,A点、C点位于左右两侧,则这三点的空间响应曲线具有相同的斜率,如图4中虚线所示,经过距离走动校正后,其空间响应曲线得到了扳正,并且A点向靠近平台(载机)飞行方向移动,距离向坐标变为R0+X-nsinθ0(X-n为A点方位位置坐标),C点向远离平台飞行方向移动,距离向坐标变为R0+Xnsinθ0(Xn为A点方位位置坐标),如图4中实线所示,相应的,校正后与B点落入同一距离单元(距离坐标为R0),不同方位位置的点目标对应不同的原始距离,这会造成校正后同一距离单元内的多普勒调频率沿方位空变。
所以,此时将R0=R0'-Xnsinθ0代入二维波数域信号S4(Kr,Kx)的表达式中,得出第一次形式变换后的二维波数域信号S5(Kr,Kx),S5(Kr,Kx)=S4(Kr,Kx),R0'表示二维波数谱支撑区被扳正后机载合成孔径雷达波束中心扫过场景中心点时的斜距。
本发明实施例中,第一次形式变换后的二维波数域信号S5(Kr,Kx)的表达式为:
S5(Kr,Kx)=Wr(Kr)Wa(Kx)exp{D1}
为了消除方位聚焦空变,达到方位统一聚焦处理,采用波数重采样的方法,通过方位向插值引入新的方位波数变量Kx',方位波数变量Kx'的表达式为:
将方位波数变量Kx'的表达式代入第一次形式变换后的二维波数域信号S5(Kr,Kx)的表达式中,得出第二次形式变换后的二维波数域信号S6(Kr,Kx'),显然,S6(Kr,Kx')=S5(Kr,Kx)=S4(Kr,Kx)。
第二次形式变换后的二维波数域信号S6(Kr,Kx')的表达式为:
其中,Kx'为新的方位波数变量,Wr(Kr)为距离窗的波数域形式,其自变量为Kr;Wa(Kx')为方位窗函数的波数域形式,其自变量为Kx';Xn为点目标相对场景中心点的方位位置,θ0为波束中心斜视角。
经过该步骤,原始斜视成像的二维波数谱已等效为正侧视的情况,二维波数谱的支撑区得以扳正,其矩形选取可用面积得到较大扩展,提高成像质量。
步骤6,对二维波数域信号S6(Kr,Kx')依次进行相位补偿、Stolt插值处理、距离向逆傅里叶变换、去斜处理、方位向逆傅里叶变换、方位波数域聚焦成像,得出最终SAR聚焦成像结果。
在步骤6中,应用改进Omega-K算法,对子孔径数据进行无模糊混叠的距离向空间位置域聚焦,方位向波数域聚焦,同时可实现方位向的统一加窗处理。
参照图5为,是本发明相位-波数变化率分布线变化示意图。参照图5a,是本发明实施例中距离聚焦处理后相位-波数变化率分布线示意图;参照图5b,是本发明实施例中统一“去斜”校正函数的相位-波数变化率分布线示意图;参照图5c,为本发明实施例中“去斜”处理后相位-波数变化率分布线示意图;图5a至图5c中,横轴表示机载合成孔径雷达方位波数变量,纵轴表示相位-位置变换率。参照图5d,为本发明实施例中“去斜”处理后进行方位IFFT后的相位-位置变换率分布线示意图;参照图5e,为本发明实施例中方位位置域校正函数的相位-位置变换率分布线示意图;参照图5f,为本发明实施例中经过方位位置域校正函数扳正后的相位-位置变换率分布线示意图;图5d至图5f中,横轴表示方位向位置变量,纵轴表示相位-位置变换率;参照图5g,为本发明实施例中最终经方位FFT在方位波数域聚焦时的相位-波数变化率分布线示意图;图5g中,横轴表示机载合成孔径雷达方位波数变量,纵轴表示相位-位置变换率。
下面结合图5a至图5g对步骤6进行具体说明。步骤6的具体子步骤为:
(6.1)对第二次形式变换后的二维波数域信号S6(Kr,Kx')进行相位补偿,得出相位补偿后信号S7(Kr,Kx')。
具体地说,对第二次形式变换后的二维波数域信号S6(Kr,Kx')乘以统一相位补偿函数HBulk(Kr,Kx'),得出相位补偿后信号S7(Kr,Kx')。统一相位补偿函数HBulk(Kr,Kx')的表达式为:
其中,Rs表示场景中心对应的参考距离,Rs为设定值。则相位补偿后信号S7(Kr,Kx')的表达式为:
其中,Xn为点目标相对场景中心点的方位位置,θ0为波束中心斜视角。在子步骤(6.2)中,相位补偿(一致补偿)处理将距离向位置以及距离徙动量均调整到以场景中心为参考。
(6.2)对相位补偿后信号S7(Kr,Kx')进行Stolt插值处理,得到Stolt插值处理后信号S8(Ky,Kx'),Ky表示Stolt插值处理时引入的距离波数变量。
具体地说,对相位补偿后信号S7(Kr,Kx')进行扩展Stolt插值处理,扩展Stolt插值处理的映射表达式为
则Stolt插值处理后信号S8(Ky,Kx')的表达式为:
其中,Ky为新的距离波数变量,Wr(Ky)是距离窗函数的波数域形式,其自变量为Ky;Wa(Kx')为方位窗函数的波数域形式,其自变量为Kx';θ0为波束中心斜视角,exp(·)为指数函数,Rs为场景中心对应的参考距离。该子步骤解耦合将距离向与方位向分离开,使得方位向可单独处理并使其适用于子孔径数据成像。
(6.3)对Stolt插值处理后信号S8(Ky,Kx')进行距离向逆傅里叶变换实现距离向聚焦,得到距离向脉压后信号为机载合成孔径雷达距离向位置变量。
具体地说,距离向脉压后信号的表达式为:
其中,为机载合成孔径雷达距离向位置变量,A'为设定的常数,Kx'为新的新的方位波数变量,Wa(Kx')为方位窗函数的波数域形式,其自变量为Kx';Xn为点目标相对场景中心点的方位位置,θ0为波束中心斜视角,Krc=4πfc/C,exp(i)为指数函数,Rs为场景中心对应的参考距离,
(6.4)对距离向脉压后信号进行去斜处理,得出去斜处理后信号对去斜处理后信号进行方位向逆傅里叶变换,得出二维位置域信号
具体地说,去斜处理的作用是消除相位-波数变化率分布线沿相位波数变化率轴的错位。距离向脉压后信号对应的相位-波数变化率分布线如图5a所示,X-n、O及Xn分别对应点A、B、C的方位位置,Xsub为子孔径数据对应的位置支撑区,超出位置支撑区的分布线将出现混叠如虚线所示,混叠后的位置如实线所示。同一距离单元内(R0'相同)的三个点目标A、B、C的三条相位-波数变化率分布线不仅在波数轴的投影是错开的,在相位波数变化率轴(纵轴)的投影也是错开的,所以在去斜处理时,引入去斜校正函数HCP(Kx'),将距离向脉压后信号乘以统一的“去斜”校正函数HCP(Kx'),得出去斜处理后信号去斜处理后信号的表达式为:
其中,“去斜”校正函数的相位-波数变化率分布线如图5b所示,去斜处理后信号的相位-波数变化率分布线如图5c所示。
对做方位IFFT处理,得出二维位置域信号二维位置域信号的表达式为:
二维位置域信号的相位-位置变换率分布线如图5d所示,此时,方位位置域的相位-位置变化率分布线在位置轴的投影是对齐的,并且具有相同的斜率
(6.5)对二维位置域信号进行方位波数域聚焦成像,得出最终SAR聚焦成像结果。
具体地说,对二维位置域信号乘以方位位置校正函数HCf(X),得出方位位置校正信号该操作“扳平”图5d所示相位-位置变化率分布线。本发明实施例中,方位位置校正函数HCf(X)的表达式为:
则方位位置校正信号的表达式为:
其中,wa(X)方位窗函数,其自变量为X;Xn为点目标相对场景中心点的方位位置,
方位位置校正函数的相位-位置变化率分布线如图5e所示,方位位置校正信号的相位-位置变换率分布线如图5f所示,为一平行于位置轴的直线,且数据在位置域已经对齐,可在该子步骤对方位向进行统一加窗进行旁瓣抑制操作。
在得出方位位置校正信号之后,对方位位置校正信号进行方位向傅里叶变换(方位FFT操作),得出最终SAR聚焦成像结果
本发明实施例中,最终SAR聚焦成像结果的相位-波数变化率分布线如图5g所示,其表达式为:
其中,B为设定的常数。
下面通过点目标仿真成像实验进一步说明本发明的正确性和有效性。
1)点目标仿真成像仿真条件
表1雷达参数
雷达系统仿真参数如表1所示,参照图6,为仿真实验中点目标仿真几何示意图;其中,X轴为SAR平台(载机)运动方向,Y轴垂直于雷达平台运动方向,Z轴为平台高度方向,SAR平台高度为H,以速度v沿X轴匀速直线飞行,设定Rs为场景中心对应的作用距离,P为场景中心点,点A、B、C均为场景边缘点,且正侧化后处于同一距离单元内。图6给出了点目标分布的空间几何模型,在地面场景中沿雷达视线投影方向布置3×3均匀点阵,沿视线投影方向和正交向的相邻点间隔为500m,波束中心线扫过P点时对应的作用距离Rs为12km,A、B、C为进行性能指标分析的远距点目标。
2)仿真内容
仿真实验1:为验证本发明的性能和有效性,这里通过本发明和标准Stolt插值的Omega-K算法(作者:Lei Zhang,公开日期:2011年11月8日,名称:Wavenumber-DomainAutofocusing for Highly Squinted UAV SAR Imagery,来源:IEEE SENSORS JOURNAL,VOL.12,NO.5)分别进行仿真实验,对比成像结果。
仿真实验2:为验证本发明的加窗特性,这里通过本发明和标准Omega-K算法分别进行仿真实验,通过在方位向加窗实验说明本发明适用于子孔径数据的加窗抑制旁瓣。
3)仿真结果分析
参照图7a,为仿真实验1中标准Omega-K算法的最终SAR聚焦成像结果示意图,参照图7b,为仿真实验1中本发明得出的最终SAR聚焦成像结果示意图,参照图7c,为仿真实验1中利用本发明进行SAR聚焦成像时的成像结果插值示意图。图7a和图7b中,横轴表示方位波数域单元,纵轴表示距离位置域单元;图7c中,横轴表示方位向,纵轴表示距离向。
对比图7a与图7b可以明显看出,标准Omega-K算法各点目标方位聚焦位置关系发生紊乱,甚至沿方位向产生反折,使得成像结果无法反应方位几何位置关系;与之对应的本发明的成像结果则能正确反应A、B、C三点的方位相互位置关系,这里B点相对于P点的方位位置偏移是由于成像几何模型引入的形变造成的,可以很容易的通过几何校正进行消除。图7c给出了本文算法各点成像结果插值图,从图中可以看出主瓣、副瓣明显分开,且呈现标准的“十字架”状,说明聚焦效果良好。
参照图8a,为仿真实验2中点A通过本发明进行成像时方位脉压剖面图在加窗处理前的结果示意图,参照图8b,为仿真实验2中点B通过本发明进行成像时方位脉压剖面图在加窗处理前的结果示意图,参照图8c,为仿真实验2中点C通过本发明进行成像时方位脉压剖面图在加窗处理前的结果示意图;参照图8d,为仿真实验2中点A通过本发明进行成像时方位脉压剖面图在加窗处理后的结果示意图,参照图8e,为仿真实验2中点B通过本发明进行成像时方位脉压剖面图在加窗处理后的结果示意图,参照图8f,为仿真实验2中点C通过本发明进行成像时方位脉压剖面图在加窗处理后的结果示意图;图8a至图8f中,横轴表示方位采样单元,纵轴表示各点的归一化幅度,单位为dB。
参照图9a,为仿真实验2中点A通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理前的结果示意图,参照图9b,为仿真实验2中点B通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理前的结果示意图,参照图9c,为仿真实验2中点C通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理前的结果示意图;参照图9d,为仿真实验2中点A通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理后的结果示意图,参照图9e,为仿真实验2中点B通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理后的结果示意图,参照图9f,为仿真实验2中点C通过标准Omega-K算法进行成像时方位脉压剖面图在加窗处理后的结果示意图;图9a至图9f中,横轴表示方位采样单元,纵轴表示各点的归一化幅度,单位为dB。
从图8a至图8f可以看出,通过本发明进行成像时加窗后各点旁瓣得到了良好的抑制,均处于-30dB以下,实现了同一距离单元内的方位统一加窗处理抑制旁瓣。从图9a至图9f可以看出,通过标准Omega-K算法进行成像时加窗后旁瓣及其他副瓣没有明显的降低,可知统一加窗无效。
表2进一步给出了A、B、C三点成像的指标参数计算结果,对比可知,本发明的成像指标与理论值基本吻合,且在边界点处仍保持良好的聚焦性能,证明了该算法的有效性,而参考算法由于经过“正侧化”处理,具有一定的聚焦能力,但方位位置错位,且严重损失方位分辨率。
表2成像算法性能指标对比(方位向)
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (5)
1.一种子孔径SAR大斜视改进Omega-K成像方法,其特征在于,包括以下步骤:
步骤1,利用机载合成孔径雷达发射线性调频信号,利用机载合成孔径雷达接收对应的回波信号;对回波信号进行解调,得到解调后的基带回波信号其中,为机载合成孔径雷达距离向位置变量,X为机载合成孔径雷达方位向位置变量,R0为机载合成孔径雷达波束中心扫过场景中心点时的斜距;对基带回波信号进行距离向傅里叶变换,得出距离波数域方位空间位置域信号S1(Kr,X),Kr为机载合成孔径雷达距离波数变量;
步骤2,将距离波数域方位空间位置域信号S1(Kr,X)进行距离向匹配滤波处理,得出距离匹配滤波后信号S2(Kr,X);
步骤3,将距离匹配滤波后信号S2(Kr,X)乘以旋转校正函数HLRWC(Kr,X),得出旋转校正后信号S3(Kr,X),其中,旋转校正函数HLRWC(Kr,X)=exp[-jKrXsinθ0],θ0为机载合成孔径雷达的波束中心斜视角;
步骤4,对旋转校正后信号S3(Kr,X)进行方位向傅里叶变换,得到二维波数域信号S4(Kr,Kx),Kx为机载合成孔径雷达方位波数变量;
步骤5,将R0=R0′-Xnsinθ0代入二维波数域信号S4(Kr,Kx)的表达式中,得出第一次形式变换后的二维波数域信号S5(Kr,Kx);令
将Kx′的表达式代入第一次形式变换后的二维波数域信号S5(Kr,Kx)的表达式中,得出第二次形式变换后的二维波数域信号S6(Kr,Kx′);R0′表示二维波数谱支撑区被扳正后机载合成孔径雷达波束中心扫过场景中心点时的斜距,Xn为点目标相对场景中心点的方位位置;
步骤6,对二维波数域信号S6(Kr,Kx′)依次进行相位补偿、Stolt插值处理、距离向逆傅里叶变换、去斜处理、方位向逆傅里叶变换、方位波数域聚焦成像,得出最终SAR聚焦成像结果。
2.如权利要求1所述的一种子孔径SAR大斜视改进Omega-K成像方法,其特征在于,在步骤1中,距离波数域方位空间位置域信号S1(Kr,X)的具体表达式为:
其中,Wr(Kr)为距离窗函数的波数域形式,其自变量为Kr;wa(X)表示方位窗函数,其自变量为X;C为光速,Krc=4πfc/C,fc为线性调频信号的载波频率,γ为线性调频信号的调频率,R(X;R0)表示任一场景散射点的瞬时斜距。
3.如权利要求1所述的一种子孔径SAR大斜视改进Omega-K成像方法,其特征在于,在步骤2中,距离匹配滤波后信号S2(Kr,X)的表达式为:
S2(Kr,X)=S1(Kr,X)·HRMF(Kr)=Wr(Kr)wa(X)exp[-jKrR(X;R0)]
其中,Wr(Kr)为距离窗函数的波数域形式,其自变量为Kr;wa(X)表示方位窗函数,其自变量为X;R(X;R0)表示任一场景散射点的瞬时斜距,HRMF(Kr)表示距离匹配滤波函数,其具体表达式为:
其中,C为光速,Krc=4πfc/C,fc为线性调频信号的载波频率,γ为线性调频信号的调频率。
4.如权利要求1所述的一种子孔径SAR大斜视改进Omega-K成像方法,其特征在于,在步骤5中,第一次形式变换后的二维波数域信号S5(Kr,Kx)的表达式为:
S5(Kr,Kx)=Wr(Kr)Wa(Kx)exp{D1}
其中,Wr(Kr)为距离窗函数的波数域形式,其自变量为Kr;Wa(Kx)为方位窗函数的波数域形式,其自变量为Kx;θ0为机载合成孔径雷达的波束中心斜视角,Xn为点目标相对场景中心点的方位位置;
在步骤5中,第二次形式变换后的二维波数域信号S6(Kr,Kx′)的表达式为:
其中,Wr(Kr)为距离窗的波数域形式,其自变量为Kr;Wa(Kx′)为方位窗函数的波数域形式,其自变量为Kx′。
5.如权利要求1所述的一种子孔径SAR大斜视改进Omega-K成像方法,其特征在于,所述步骤6的具体子步骤为:
(6.1)对第二次形式变换后的二维波数域信号S6(Kr,Kx′)进行相位补偿,得出相位补偿后信号S7(Kr,Kx′);所述相位补偿后信号S7(Kr,Kx′)的表达式为:
其中,Wr(Kr)为距离窗的波数域形式,其自变量为Kr;Wa(Kx′)为方位窗函数的波数域形式,其自变量为Kx′;Xn为点目标相对场景中心点的方位位置,θ0为机载合成孔径雷达的波束中心斜视角,Rs表示场景中心对应的参考距离,Rs为设定值;片Bulk(Kr,Kx′)为:
(6.2)对相位补偿后信号S7(Kr,Kx′)进行Stolt插值处理,得到Stolt插值处理后信号S8(Ky,Kx′),Ky表示Stolt插值处理时引入的距离波数变量;
在对相位补偿后信号S7(Kr,Kx′)进行Stolt插值处理时,Stolt插值处理的映射表达式为
Stolt插值处理后信号S8(Ky,Kx′)的表达式为:
其中,Wr(Ky)是距离窗函数的波数域形式,其自变量为Ky;Wa(Kx′)为方位窗函数的波数域形式,其自变量为Kx′;
(6.3)对Stolt插值处理后信号S8(Ky,Kx′)进行距离向逆傅里叶变换,得到距离向脉压后信号
(6.4)对距离向脉压后信号进行去斜处理,得出去斜处理后信号对去斜处理后信号进行方位向逆傅里叶变换,得出二维位置域信号
所述去斜处理后信号的表达式为:
其中,A′为设定的常数;HCP(Kx′)为:
其中,Krc=4πfc/C,fc为线性调频信号的载波频率,C为光速;
(6.5)对二维位置域信号乘以方位位置校正函数HCf(X),得出方位位置校正信号方位位置校正函数HCf(X)的表达式为:
其中,θ为瞬时斜视角,λ为机载合成孔径雷达发射信号的载波波长;则方位位置校正信号的表达式为:
其中,wa(X)表示方位窗函数,其自变量为X;
最终SAR聚焦成像结果表示为最终SAR聚焦成像结果的表达式为:
其中,B为设定的常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510051776.0A CN104597447B (zh) | 2015-01-30 | 2015-01-30 | 一种子孔径SAR大斜视改进Omega‑K成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510051776.0A CN104597447B (zh) | 2015-01-30 | 2015-01-30 | 一种子孔径SAR大斜视改进Omega‑K成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104597447A CN104597447A (zh) | 2015-05-06 |
CN104597447B true CN104597447B (zh) | 2017-03-08 |
Family
ID=53123353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510051776.0A Active CN104597447B (zh) | 2015-01-30 | 2015-01-30 | 一种子孔径SAR大斜视改进Omega‑K成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104597447B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105759267B (zh) * | 2016-03-11 | 2019-07-09 | 南京航空航天大学 | 一种大斜视角SAR的改进Omega-K成像方法 |
CN106054187B (zh) * | 2016-05-23 | 2018-11-09 | 西安电子科技大学 | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 |
CN106054152B (zh) * | 2016-05-23 | 2018-11-09 | 西安电子科技大学 | 基于逆扩展Omega-K算法的非理想轨迹SAR回波获取方法 |
CN106249237B (zh) * | 2016-07-19 | 2019-05-21 | 西安电子科技大学 | 一种曲线轨迹下大斜视sar频域成像方法 |
CN106908789B (zh) * | 2017-02-28 | 2019-08-02 | 西安电子科技大学 | 一种基于空间域支撑区融合的sar成像方法 |
CN107219524B (zh) * | 2017-05-11 | 2020-01-07 | 西北工业大学 | 一种基于全局最小相位近似的sar成像优化方法 |
CN107290747B (zh) * | 2017-06-12 | 2019-02-19 | 北京理工大学 | 一种大前斜视合成孔径雷达成像方法 |
CN112578383B (zh) * | 2020-11-19 | 2023-06-30 | 西安电子科技大学 | 一种基于扩展波数谱重构的高机动平台tops sar成像方法 |
CN114325704B (zh) * | 2021-12-31 | 2023-07-25 | 电子科技大学 | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 |
CN116736304B (zh) * | 2023-08-10 | 2023-10-20 | 西安电子科技大学 | 一种基于线性调频信号的随机跳频信号成像方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197317B (zh) * | 2013-03-26 | 2014-09-17 | 西安电子科技大学 | 基于fpga的sar成像方法 |
CN103675815B (zh) * | 2013-09-27 | 2017-06-20 | 西安电子科技大学 | 一种在大斜视sar成像模式下对多普勒调频率进行精确估计的方法 |
CN103760561A (zh) * | 2014-01-26 | 2014-04-30 | 西安电子科技大学 | 大斜视角sar成像模式下提高方位向非散焦长度的方法 |
-
2015
- 2015-01-30 CN CN201510051776.0A patent/CN104597447B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104597447A (zh) | 2015-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104597447B (zh) | 一种子孔径SAR大斜视改进Omega‑K成像方法 | |
CN111142105B (zh) | 复杂运动目标isar成像方法 | |
CN104316924B (zh) | 机载超高分辨率sar反投影图像的自聚焦运动补偿方法 | |
CN102288963B (zh) | 基于子孔径参数估计的双基地isar图像融合方法 | |
CN108459321B (zh) | 基于距离-方位圆模型的大斜视高分辨率sar成像方法 | |
CN105974414A (zh) | 基于二维自聚焦的高分辨聚束sar自聚焦成像方法 | |
CN109738894B (zh) | 一种用于大视场合成孔径雷达的高斜视多角度成像方法 | |
CN106249237A (zh) | 一种曲线轨迹下大斜视sar频域成像方法 | |
CN107229048A (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
CN102323581B (zh) | 斜视聚束合成孔径雷达的成像方法 | |
CN106324597A (zh) | 基于pfa的大转角isar雷达的平动补偿和成像方法 | |
CN104777479B (zh) | 基于多核dsp的前侧视sar实时成像方法 | |
CN108710111B (zh) | 一种机载双基前视sar方位相位的二维空变校正方法 | |
CN110346798B (zh) | 一种双基合成孔径雷达波数域高效成像处理方法 | |
CN106610492A (zh) | 基于rd算法的时频域混合校正距离徙动的sar成像方法 | |
CN110412587B (zh) | 一种基于解卷积的下视合成孔径三维成像方法及系统 | |
CN102608576B (zh) | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 | |
CN112859018B (zh) | 一种基于图像几何校正的视频sar成像方法 | |
CN105093224A (zh) | 一种大斜视合成孔径雷达成像处理方法 | |
CN102788978B (zh) | 一种斜视星机双基地合成孔径雷达成像方法 | |
CN103454638A (zh) | 一种圆迹合成孔径雷达三维层析成像方法 | |
CN115685200A (zh) | 一种高精度大前斜视sar成像运动补偿与几何校正方法 | |
CN106054187A (zh) | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 | |
CN109444882A (zh) | 基于变斜视椭圆波束同步模型的双站sar成像方法 | |
CN114545411A (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 |
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 |