CN103412309B - 移不变双基地前视合成孔径雷达nlcs成像方法 - Google Patents

移不变双基地前视合成孔径雷达nlcs成像方法 Download PDF

Info

Publication number
CN103412309B
CN103412309B CN201310375801.1A CN201310375801A CN103412309B CN 103412309 B CN103412309 B CN 103412309B CN 201310375801 A CN201310375801 A CN 201310375801A CN 103412309 B CN103412309 B CN 103412309B
Authority
CN
China
Prior art keywords
mrow
msub
eta
msup
mfrac
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
Application number
CN201310375801.1A
Other languages
English (en)
Other versions
CN103412309A (zh
Inventor
武俊杰
李中余
黄钰林
孙稚超
杨建宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201310375801.1A priority Critical patent/CN103412309B/zh
Publication of CN103412309A publication Critical patent/CN103412309A/zh
Application granted granted Critical
Publication of CN103412309B publication Critical patent/CN103412309B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种移不变双基地前视合成孔径雷达NLCS成像方法,具体将四次滤波和距离向四阶非线性Chirp-Scaling相结合,解决了移不变双基地前视SAR距离徙动(RCM)随距离向的非线性空变及二次距离压缩(SRC)随距离向的线性空变问题。本发明方法的特点:首先在二维频域对回波进行四次滤波,增加了两维自由度,其次在方位频域距离时域进行距离向四阶非线性Chirp-Scaling操作,从而消除了RCM随距离向的非线性空变及SRC随距离向的线性空变,可以实现移不变双基地前视SAR距离徙动的精确校正和SRC空变性的去除,完成双基地前视SAR的精确聚焦。

Description

移不变双基地前视合成孔径雷达NLCS成像方法
技术领域
本发明属于雷达技术领域,具体涉及合成孔径雷达成像技术中的移不变双基前视SAR的成像方法。
背景技术
合成孔径雷达(SyntheticApertureRadar,SAR)是一种全天时、全天候的现代高分辨率微波遥感成像雷达,它利用雷达天线和目标区域间的相对运动来获得空间的高分辨率。在地形测绘、植被分析、海洋及水文观测、环境及灾害监视、资源勘探以及地壳微变检测等领域,合成孔径雷达发挥了越来越重要的作用。但是由于本身工作体制的限制,现有单基地SAR并不能实现飞行器前视区域的高分辨成像,从而使SAR技术在飞行器前视对地、自主着陆、物资空投等方面不能充分的发挥作用。
双基地前视SAR(BFSAR)是一种新的雷达体制,系统发射站和接收站分置于不同平台上,通过合理的几何配置,可以实现在接收站正前方进行高分辨成像。另外,收发分置的特点使其具备了许多突出的优点,如获取目标信息丰富、作用距离远、安全性好、抗干扰能力强等。
移不变BFSAR是BFSAR的一种模式,是指发射站与接收站以相同的速度、沿着平行的航迹飞行的模式。与传统SAR相比,在移不变BFSAR中,由于接收站对前视区域的对称性,距离徙动沿着距离向具有非线性空变,SRC沿着距离向具有较大的线性空变等特点。
在文献“R.K.Raney,H.Runge,R.Bamler,I.G.Cumming,andF.H.Wong:PrecisionSARprocessingusingchirpscaling,IEEETrans.Geosci.RemoteSens.,vol.32,no.4,pp.786-799,Jul.1994”中,应用了传统线性调频变标(Chirp-Scaling,CS)方法对BFSAR进行成像处理,但是线性Chirp—Scaling方法忽略了二次距离压缩的空变性,而在BFSAR中,距离徙动沿着距离向具有非线性空变特点,将直接影响最终的成像效果。
在文献“X.Qiu,D.Hu,andC.Ding:SomereflectionsonbistaticSARofforward-lookingconfiguration,IEEEGeosci.RemoteSens.Lett.,vol.5,no.4,pp.735-739,2008”中,应用了距离多普勒方法对BFSAR进行成像处理,但是该方法没有考虑SRC的空变性,而在BFSAR中,SRC沿着距离向具有较大的线性空变性,将直接影响最终的成像效果。
在文献“Z.HuaandL.Xingzhao:AnExtendedNonlinearChirp-ScalingAlgorithmforFocusingLarge-BaselineAzimuth-InvariantBistaticSARData,IEEEGeosci.RemoteSens.Lett.,vol.6,no.3,pp.548-552,2009”中,应用了非线性CS方法,但也忽略了BFSAR中SRC沿着距离向较大的线性空变问题。
发明内容
本发明的目的是针对现有技术存在的上述缺陷,研究设计一种移不变BFSAR非线性Chirp—Scaling(NLCS)成像方法。
为了方便描述本发明的内容,首先对以下术语进行解释:
术语1:双基地SAR
双基地SAR是指系统发射站和接收站分置于不同平台上的SAR系统,其中至少有一个平台为运动平台,在概念上属于双基地雷达。
术语2:移不变双基地前视SAR
移不变双基地前视SAR是双基地SAR的一种,收发站平行飞行,且速度相同,发射站正侧视,接收站前视。
本发明的技术方案为:一种移不变双基地前视SARNLCS成像方法,具体包括如下步骤:
步骤一:获取回波数据的二维频谱,
在直角坐标系中设发射平台位置为(xT,yT,hT),接收站零时刻位置记为(xR,yR,hR),发射站和接收站速度均为v,并沿y轴平行飞行,任意成像点坐标记为P(x,y),双基地距离和为Rb(η;x,y)=RT(η;x,y)+RR(η;x,y),其中,η为方位时间,RT(η;x,y),RR(η;x,y)分别为发射站和接收站的距离历程:
R T ( η ; x , y ) = ( r T cos θ sT ) 2 + v 2 ( η - y v ) 2 - 2 r T v ( η - y v ) tan θ sT
R R ( η ; x , y ) = ( r R cos θ sR ) 2 + v 2 ( η - y v ) 2 - 2 r R v ( η - y v ) tan θ sR
其中,θsT和θsR分别为发射站和接收站的斜视角,rT和rR分别为发射站和接收站的最短斜距,表达式分别为 r T = ( x - x T ) 2 + h T 2 , r R = ( x - x R ) 2 + h R 2 ;
原始回波数据在距离频域、方位时域的表达式为:
S ( f , η ; x , y )
= S 0 ( f ) exp { - j 2 π ( f + f 0 ) R T ( η ; x , y ) + R R ( η ; x , y ) c }
其中,f为距离频率变量,f0为发射信号载频,c为光速,Br为发射信号带宽,Kr为发射信号的调频斜率,rect[·]代表距离时间窗;
基于广义Loffeld模型,得到原始回波在二维频域的表达式为:
s2df(f,fη;x,y)=exp{jΦG(f,fη;x,y)}
其中,fη为方位频率变量,二维频谱相位为:
Φ G ( f , f η ; x , y ) = - πf 2 K r - φ T ( η PT ) - φ R ( η PR )
= - π f 2 K r - 2 π c [ r T F T ( f , f η ) + r R F R ( f , f η ) ]
- 2 π v [ r T f ηT ( f η ) tan θ sT + r R f ηR ( f η ) ta n - 1 θ dR ]
- 2 π f η y v
F T ( f , f η ) = ( f + f 0 ) 2 - ( c f ηT ( f η ) v ) 2
F R ( f , f η ) = ( f + f 0 ) 2 - ( c f ηR ( f η ) v ) 2
其中,fηT(fη)和fηR(fη)分别为发射站和接收站的对总的多普勒频率的贡献,表达式如下:
f ηT ( f η ) = f ηcT + f ηrT f ηr ( f η - f ηc ) - f ηrT f η 3 - f η 3 T f ηr f ηr 3 ( f η - f ηc ) 2 · · ·
f ηR ( f η ) = f ηcR + f ηrR f ηr ( f η - f ηc ) - f ηrR f η 3 - f η 3 R f ηr f ηr 3 ( f η - f ηc ) 2 · · · ;
其中,fηcT,fηcR分别为发射站和接收站对应的多普勒质心;fηrT,fηrR分别为发射站和接收站对应的多普勒调频斜率;fη3T,fη3R分别为发射站和接收站对应的多普勒三阶调频率;fηc,fηr和fη3为系统总的多普勒质心、多普勒调频斜率和多普勒三阶调频率;
然后ΦG(f,fη;x,y)被分解为:
ΦG(f,fη;x,y)=ΦRCM(f,fη;x)+ΦRC(f,fη;x)
3rd(f,fη;x)+Φ4th(f,fη;x)
AC(fη;x)+ΦAL(fη;x,y)
其中,φRCM(f,fη;x)是随距离频率变化的线性项,代表距离徙动因子和沿距离方向的位置;ΦRC(f,fη;x)是距离频率的二阶项,代表沿距离方向的调制;Φ3rd(f,fη;x)、Φ4th(f,fη;x)是距离频率的三阶项和四阶项,代表距离和方位之间的耦合;ΦAC(fη;x)是方位压缩因子,该因子随距离变化而变化,具体表达式如下:
Φ RCM ( f , f η ; x ) = - 2 π c ( r T D ( f ηT ( f η ) ) + r R D ( f ηR ( f η ) ) ) f
Φ RC ( f , f η ; x ) = - π f 2 K r + π ( r T f ηT 2 ( f η ) D 3 ( f ηT ( f η ) ) + r R f ηR 2 ( f η ) D 3 ( f ηR ( f η ) ) ) cf 2 v 2 f 0 3 = - π f 2 K m ( f η ; x )
Φ 3 rd ( f , f η ; x ) = - π ( r T f ηT 2 ( f η ) D 5 ( f ηT ( f η ) ) + r R f ηR 2 ( f η ) D 5 ( f ηR ( f η ) ) ) cf 3 v 2 f 0 4 = π K c ( f η ; x ) f 3
Φ 4 rd ( f , f η ; x ) = π r T f ηT 2 ( f η ) ( c 2 f ηT 2 ( f η ) + 4 f 0 2 v 2 ) D 7 ( f ηT ( f η ) ) + r R f ηR 2 ( f η ) ( c 2 f ηR 2 ( f η ) + 4 f 0 2 v 2 ) D 7 ( f ηR ( f η ) ) cf 4 16 v 4 f 0 7
= π K 4 th ( f η ; x ) f 4
Φ 4 rd ( f , f η ; x ) = - 2 π f 0 c [ r T D ( f ηT ( f η ) ) + r R D ( f ηR ( f η ) ) ]
Φ AL ( f η ; x , y ) = - 2 π υ [ r T f ηT ( f η ) tan θ sT + r R f ηR ( f η ) tan - 1 θ dR ] - 2 π f η y υ
其中, D ( f ηT ( f η ) ) = 1 - ( cf ηT ( f η ) / vf 0 ) 2 , D ( f ηR ( f η ) ) = 1 - ( cf ηR ( f η ) / vf 0 ) 2 .
将Km(fη;x)沿着距离向泰勒展开到二次,可得:
K m ( f η ; x ) = K m ref + k 1 Δτ + k 2 Δ τ 2
其中,为参考点的调频率,k1、k2分别为Km(fη;x)沿距离向泰勒展开的一次项和二次项系数,τd(fη)为任意目标点的距离徙动量,为参考目标点的距离徙动量;
将Kc(fη;x)沿着距离向泰勒展开到一次,可得:
K c ( f η ; x ) = K c ref + k c Δτ
其中,为参考点的三阶调频率,kc为Kc(fη;x)沿距离向泰勒展开的一次项系数;
步骤二:四次相位滤波;
在步骤一对回波数据进行二维傅立叶变换之后,利用一个四次滤波器对其进行滤波,其中滤波器的表达式为:
H 1 ( f ) = exp { j 2 π 3 Y 1 ( f η ) f 3 + j 2 π 4 Y 2 ( f η ) f 4 }
其中,Y1(fη)和Y2(fη)是H1(f)的系数;
假设
2 π 3 Y 1 m ( f η ) = 2 π 3 Y 1 ( f η ) + π K c ref + π k c Δτ
2 π 4 Y 2 m ( f η ) = 2 π 4 Y 2 ( f η ) + π K 4 th ( f η , x )
可得:
Y 1 m ( f η ) = k 1 ( a - 0.5 ) - K m ref a 2 β ( K m ref ) 3 ( a - 1 ) Y 2 m ( f η ) = 3 ( a - 0.5 ) ( K m ref ) 3 k c + 6 ( a - 0.5 ) Y 1 m ( f η ) ( K m ref ) 2 k 1 - ak 2 - 2 a Y 1 m ( f η ) β ( K m ref ) 3 - 2 a q 3 β 3 ( a - 1 ) ( K m ref ) 4
其中, α = m 1 ( f η ) m 1 ( f η 0 ) , β = 1 m 1 2 ( f η ) [ m 2 ( f η 0 ) - m 1 ( f η 0 ) m 1 ( f η ) m 2 ( f η ) ] ,
m 1 ( f η ) = 1 cD ( f ηT ( f η ) ) - r T 0 k T 1 ( f ηT ref ( f η ) ) c D 2 ( f ηT ( f η ) ) + a T 1 cD ( f ηR ( f η ) ) - r R 0 k R 1 ( f ηR ref ( f η ) ) c D 2 ( f ηR ( f η ) )
m 2 ( f η ) = a T 2 2 cD ( f ηR ( f η ) ) - 2 k T 1 ( f ηT ref ( f η ) ) + r T 0 k T 2 ( f ηT ref ( f η ) ) 2 c D 2 ( f ηT ( f η ) )
+ r T 0 k T 1 2 ( f ηT ref ( f η ) ) cD 3 ( f ηT ( f η ) ) - 2 a R 1 k R 1 ( f ηR ref ( f η ) ) + r R 0 k R 2 ( f ηR ref ( f η ) ) 2 cD 2 ( f ηR ( f η ) )
+ r R 0 k R 1 2 ( f ηR ref ( f η ) ) cD 3 ( f ηR ( f η ) )
m 1 ( f η 0 ) = m 1 ( f η ) | f η = f η 0 , m 1 ( f η 0 ) = m 1 ( f η ) | f η = f η 0 ,
其中, 是D(fηT)对rT的一阶和二阶导数; 是D(fηR)对rT的一阶和二阶导数;aT1、aT2为rR对rT的一阶和二阶导数,η0是预先设定的参考值,fη0为一个方位参考频率;
利用上述四次滤波器H1(f)对回波二维频谱进行滤波,并将滤波后的二维频谱变换到距离多普勒域,其表达式为:
Sfilter(τ,fη)=exp{jΦRD(τ,fη)}
其中,τ为距离时间变量,相位ΦRD(τ,fη)表达为式:
Φ RD ( τ , f η ) = π K m ( f η ; x ) [ τ - τ d ( f η ) ] 2 + 2 π 3 Y 1 m ( f η ) K m 3 ( f η ; x ) [ τ - τ d ( f η ) ] 3
+ 2 π 4 Y 2 m ( f η ) K m 4 ( f η ; x ) [ τ - τ d ( f η ) ] 4 + Φ AC ( f η ; x ) + Φ AL ( f η ; x , y )
步骤三:距离向四阶非线性Chirp-Scaling处理;
将四次相位滤波后的距离多普勒域信号乘以四阶非线性CS因子,非线性CS因子表达式为:
SCS(τ,fη)=exp{jΦCS(τ,fη)}
其相位ΦCS(τ,fη)为:
Φ CS ( τ , f η ) = π q 2 ( f η ) [ τ - τ d ref ( f η ) ] 2 + 2 π 3 q 3 ( f η ) [ τ - τ d ref ( f η ) ] 3
+ 2 π 4 q 4 ( f η ) [ τ - τ d ref ( f η ) ] 4
其中,q2(fη),q3(fη)和q4(fη)是二阶、三阶和四阶项系数,表达式为:
q 2 ( f η ) = K m ref ( a - 1 ) q 3 ( f η ) = k 1 2 ( a - 1 ) - K m ref a 2 β q 4 ( f η ) = - 1 3 k 2 a + ( a - 1 ) ( K m ref ) 3 k c 2 + ( a - 1 ) Y 1 m ( f η ) ( K m ref ) 2 k 1 - 2 3 Y 1 m ( f η ) aβ ( K m ref ) 3 - 2 3 q 3 aβ
步骤四:RCM校正,距离压缩和高阶项处理;
经过非线性Chirp-Scaling处理后距离多普勒信号,距离徙动、二次距离压缩和三阶项的距离向空变性已经去除,因此,距离徙动、二次距离压缩和三阶项可以直接在二维频域来处理,校正函数为:
S jz ( f , f η ) = exp { j Φ 2 D CS ( f , f η ) }
其相位为:
Φ RD CS ( f , f η ) = 2 πΔ τ ref f - π f 2 α K m ref + 2 π [ Y 1 m ref ( f η ) ( K m ref ) 3 + q 3 ] 3 α 3 ( K m ref ) 3 f 3
+ 2 π [ Y 2 m ( f η ) ( K m ref ) 4 + q 4 ] 4 α 4 ( K m ref ) 4 f 4 + Φ AC ( f η ; x ) + Φ AL ( f η ; x , y )
步骤五:方位向压缩;
将上述RCM校正、距离压缩和高阶项处理后的数据变换到距离时域方位频域,接下来将进行方位压缩,方位压缩函数为:
SAC(τ,fη)=exp{j[ΦAC(fη;x)+ΦAL(fη;x,y)]}
完成方位压缩后,进行方位向傅里叶反变换得到时域聚焦图像,从而完成移不变双基地前视SAR的聚焦成像。
本发明的有益效果:本发明的方法首先在二维频域对回波进行四次滤波,增加了两维自由度,其次在方位频域距离时域进行距离向四阶非线性Chirp-Scaling操作,从而同时消除距离徙动随距离向的非线性空变及二次距离压缩随距离向的线性空变,实现移不变双基地前视SAR距离徙动的精确校正和SRC空变的去除,完成双基地前视SAR的精确聚焦。本发明的方法解决了现有BFSAR成像算法无法同时去除距离徙动(RCM)和二次距离压缩SRC的空变问题,实现BFSAR的精确聚焦成像。本发明的方法采用将四次滤波和距离向四阶非线性Chirp-Scaling相结合的思想,有效解决了移不变双基地前视SAR距离徙动(RCM)随距离向的非线性空变及二次距离压缩(SRC)随距离向的线性空变问题。本方法的特点是首先在二维频域对回波进行四次滤波,增加了两维自由度,其次在方位频域距离时域进行距离向四阶非线性Chirp-Scaling操作,从而消除了RCM随距离向的非线性空变及SRC随距离向的线性空变,可以实现移不变双基地前视SAR距离徙动的精确校正和SRC空变性的去除,完成双基地前视SAR的精确聚焦;同时本发明方法的精度高,运算效率较高,还能够满足BFSAR实时成像处理的需求。
附图说明
图1为本发明实施例方法的流程示意图。
图2为本发明实施例的系统结构图。
图3为本发明实施例的所使用的参数表图。
图4为本发明实施例的成像场景示意图。
图5为本发明实施例的成像结果图。
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2010上验证正确。下面就具体实施方式对本发明作进一步的详细描述。
本实例中采用的系统结构如图2所示,系统坐标系以成像中心点目标O位坐标原点,双平台沿着y轴以相同的速度平飞,x轴为切航迹方向,z轴为垂直地面方向。
本发明的移不变双基地前视合成孔径雷达NLCS成像算法的流程示意图如图1所示,
具体过程如下:
步骤一:成像系统参数初始化
系统参数列表如图3所示。发射站的位置坐标为(-14000,-10000,12000)m,接收站零时刻位置坐标为(0,-14000,10000)m,波速中心位于场景坐标原点处时记为零时刻,平台速度为200m/s,系统载频为9.6GHz,发射信号带宽为80MHz,脉冲重复频率为800Hz。
本发明实施采用的目标场景如图4所示,图中的黑色圆点为布置于地面上的15个点目标,每相邻两点沿着x轴的距离为750m,沿着y轴的距离为300m,中心点O位于坐标原点。
步骤二:获取目标回波
采用图1所示的系统采用图3中给出的系统参数,对图4中的点目标逐点生成回波数据,从而生成目标回波,将生成的目标回波分别做方位向傅里叶变换和距离向傅里叶变换,得到回波的二维频谱,记为S2df(f,fη;x,y)。
步骤三:四次相位滤波
对步骤二得到的回波二维频谱S2df(f,fη;x,y)进行四次相位滤波,四次滤波器表达式为:
H 1 ( f ) = exp { j 2 π 3 Y 1 ( f η ) f 3 + j 2 π 4 Y 2 ( f η ) f 4 }
将滤波后的二维频谱变换到距离多普勒域,其表达式为
Sfilter(τ,fη)=exp{jΦRD(τ,fη)}
相位ΦRD(τ,fη)表达为式:
Φ RD ( τ , f η ) = π K m ( f η ; x ) [ τ - τ d ( f η ) ] 2
+ 2 π 3 Y 1 m ( f η ) K m 3 ( f η ; x ) [ τ - τ d ( f η ) ] 3
+ 2 π 4 Y 2 m ( f η ) K m 4 ( f η ; x ) [ τ - τ d ( f η ) ] 4
+ Φ AC ( f η ; x ) + Φ AL ( f η ; x , y )
步骤四:距离向四阶非线性Chirp-Scaling处理
对四次滤波后的距离多普勒域信号Sfilter(τ,fη)进行四阶非线性Chirp-Scaling处理,非线性Chirp-Scaling因子表达式为:
S cs ( τ , f η ) = exp { j Φ CS ( τ , f η ) }
其相位Φcs(τ,fη)为:
Φ CS ( τ , f η ) = π q 2 ( f η ) [ τ - τ d ref ( f η ) ] 2 + 2 π 3 q 3 ( f η ) [ τ - τ d ref ( f η ) ] 3
+ 2 π 4 q 4 ( f η ) [ τ - τ d ref ( f η ) ] 4
距离向四阶非线性Chirp-Scaling处理处理后信号RCM、SRC和三次项的距离向空变都得到消除。
步骤五:RCM校正,距离压缩和高阶项处理
经过非线性CS处理后距离多普勒信号,RCM,SRC和三阶项的距离向空变性已经去除,因此一致RCM,SRC和三阶项可以直接在二维频域来处理。校正函数为
S jz ( f , f η ) = exp { j Φ 2 D CS ( f , f η ) }
其相位为:
Φ RD CS ( f , f η ) = 2 πΔ τ ref f - π f 2 α K m ref
+ 2 π [ Y 1 m ref ( f η ) ( K m ref ) 3 + q 3 ] 3 α 3 ( K m ref ) 3 f 3
+ 2 π [ Y 2 m ( f η ) ( K m ref ) 4 + q 4 ] 4 α 4 ( K m ref ) 4 f 4
+ Φ AC ( f η ; x ) + Φ AL ( f η ; x , y )
步骤六:方位向压缩
将上述RCM校正、距离压缩和高阶项处理后的数据变换到距离时域方位频域,接下来将进行方位压缩,方位压缩函数为
SAC(τ,fη)=exp{j[ΦAC(fη;x)+ΦAL(fη;x,y)]}
完成方位压缩后,进行方位向傅里叶反变换得到时域聚焦图像,从而完成移不变双基地前视SAR的聚焦成像,其成像结果如图5所示。从图中可以看出,本发明提供的方法可以很好的实现移不变双基地前视SAR成像数据处理。
通过本发明具体实施方式可以看出,本发明解决了移不变双基地前视SAR距离徙动随距离向的非线性空变、二次距离压缩随距离向的线性空变问题,可以实现移不变双基地前视SAR目标回波较好的聚焦成像。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (1)

1.一种移不变双基地前视SARNLCS成像方法,具体包括如下步骤:
步骤一:获取回波数据的二维频谱,
在直角坐标系中设发射平台位置为(xT,yT,hT),接收站零时刻位置记为(xR,yR,hR),发射站和接收站速度均为v,并沿y轴平行飞行,任意成像点坐标记为P(x,y),双基地距离和为Rb(η;x,y)=RT(η;x,y)+RR(η;x,y),其中,η为方位时间,RT(η;x,y),RR(η;x,y)分别为发射站和接收站的距离历程:
R T ( η ; x , y ) = ( r T cos θ sT ) 2 + v 2 ( η - y v ) 2 - 2 r T v ( η - y v ) tan θ sT
R R ( η ; x , y ) = ( r R cos θ sR ) 2 + v 2 ( η - y v ) 2 - 2 r R v ( η - y v ) tan θ sR
其中,θsT和θsR分别为发射站和接收站的斜视角,rT和rR分别为发射站和接收站的最短斜距,表达式分别为 r T = ( x - x T ) 2 + h T 2 , r R = ( x - x R ) 2 + h R 2 ;
原始回波数据在距离频域、方位时域的表达式为:
S ( f , η ; x , y ) = S 0 ( f ) exp { - j 2 π ( f + f 0 ) R T ( η ; x , y ) + R R ( η ; x , y ) c }
其中,f为距离频率变量,f0为发射信号载频,c为光速,Br为发射信号带宽,Kr为发射信号的调频斜率,rect[·]代表距离时间窗;
基于广义Loffeld模型,得到原始回波在二维频域的表达式为:
S2df(f,fη;x,y)=exp{jΦG(f,fη;x,y)}
其中,fη为方位频率变量,二维频谱相位为:
Φ G ( f , f η ; x , y ) = - πf 2 K r - 2 π c [ r T F T ( f , f η ) + r R F R ( f , f η ) ] - 2 π v [ r T f ηT ( f η ) tan θ sT + r R r ηR ( f η ) tan - 1 θ dR - 2 π f η y v
其中,θdR=90°-θsR
F T ( f , f η ) = ( f + f 0 ) 2 - ( cf ηT ( f η ) v ) 2
F R ( f , f η ) = ( f + f 0 ) 2 - ( cf ηR ( f η ) v ) 2
其中,fηT(fη)和fηR(fη)分别为发射站和接收站的对总的多普勒频率的贡献,表达式如下:
f ηT ( f η ) = f ηcT + f ηrT f ηr ( f η - f ηc ) - f ηrT f η 3 - f η 3 T f ηr f ηr 3 ( f η - f ηc ) 2 · · ·
f ηR ( f η ) = f ηcR + f ηrR f ηr ( f η - f ηc ) - f ηrR f η 3 - f η 3 R f ηr f ηr 3 ( f η - f ηc ) 2 · · · ;
其中,fηcT,fηcR分别为发射站和接收站对应的多普勒质心;fηrT,fηrR分别为发射站和接收站对应的多普勒调频斜率;fη3T,fη3R分别为发射站和接收站对应的多普勒三阶调频率;fηc,fηr和fη3为系统总的多普勒质心、多普勒调频斜率和多普勒三阶调频率;
然后ΦG(f,fη;x,y)被分解为:
Φ G ( f , f η ; x , y ) = Φ RCM ( f , f η ; x ) + Φ RC ( f , f η ; x ) + Φ 3 rd ( f , f η ; x ) + Φ 4 th ( f , f η ; x ) + Φ AC ( f η ; x ) + Φ AL ( f η ; x , y )
其中,ΦRCM(f,fη;x)是随距离频率变化的线性项,代表距离徙动因子和沿距离方向的位置;ΦRC(f,fη;x)是距离频率的二阶项,代表沿距离方向的调制;Φ3rd(f,fη;x)、Φ4th(f,fη;x)是距离频率的三阶项和四阶项,代表距离和方位之间的耦合;ΦAC(fη;x)是方位压缩因子,该因子随距离变化而变化,具体表达式如下:
Φ RCM ( f , f η ; x ) = - 2 π c ( r T D ( f ηT ( f η ) ) + r R D ( f ηR ( f η ) ) ) f
Φ RC ( f , f η ; x ) = - π f 2 K r + π ( r T f ηT 2 ( f η ) D 2 ( f ηT ( f η ) ) + r R f ηR 2 ( f η ) D 3 ( f ηR ( f η ) ) ) cf 2 v 2 f 0 3 = - πf 2 K m ( f η ; x )
Φ 3 rd ( f , f η ; x ) = - π ( r T f ηT 2 ( f η ) D 5 ( f ηT ( f η ) ) + r R f ηR 2 ( f η ) D 5 ( f ηR ( f η ) ) ) cf 3 v 2 f 0 4 = πK c ( f η ; x ) f 3
Φ 4 th ( f , f η ; x ) = π r T f ηT 2 ( f η ) ( c 2 f ηT 2 ( f η ) + 4 f 0 2 v 2 ) D 7 ( f ηT ( f η ) ) + r R f ηR 2 ( f η ) ( c 2 f ηR 2 ( f η ) + 4 f 0 2 v 2 ) D 7 ( f ηR ( f η ) ) cf 4 16 v 4 f 0 7 = π K 4 th ( f η ; x ) f 4
Φ AC ( f η ; x ) = - 2 π f 0 c [ r T D ( f ηT ( f η ) ) + r R D ( f ηR ( f η ) ) ]
Φ AL ( f η ; x , y ) = - 2 π v [ r T f ηT ( f η ) tan θ sT + r R f ηR ( f η ) tan - 1 θ dR ] - 2 π f η y v
其中, D ( f ηT ( f η ) ) = 1 - ( cf ηT ( f η ) / vf 0 ) 2 , D ( f ηR ( f η ) ) = 1 - ( cf ηR ( f η ) / vf 0 ) 2 .
将Km(fη;x)沿着距离向泰勒展开到二次,可得:
K m ( f η ; x ) = K m ref + k 1 Δτ + k 2 Δτ 2
其中,为参考点的调频率,k1、k2分别为Km(fη;x)沿距离向泰勒展开的一次项和二次项系数,τd(fη)为任意目标点的距离徙动量,为参考目标点的距离徙动量;
将Kc(fη;x)沿着距离向泰勒展开到一次,可得:
K c ( f η ; x ) = K c ref + k c Δτ
其中,为参考点的三阶调频率,kc为Kc(fη;x)沿距离向泰勒展开的一次项系数;
步骤二:四次相位滤波;
在步骤一对回波数据进行二维傅立叶变换之后,利用一个四次滤波器对其进行滤波,其中滤波器的表达式为:
H 1 ( f ) = exp { j 2 π 3 Y 1 ( f η ) f 3 + j 2 π 4 Y 2 ( f η ) f 4 }
其中,Y1(fη)和Y2(fη)是H1(f)的系数;
假设
3 π 3 Y 1 m ( f η ) = 2 π 3 Y 1 ( f η ) + π K c ref + πk c Δτ
2 π 4 Y 2 m ( f η ) = 2 π 4 Y 2 ( f η ) + π K 4 th ( f η , x )
可得:
Y 1 m ( f η ) = k 1 ( a - 0.5 ) - K m ref a 2 β ( K m ref ) 3 ( a - 1 ) Y 2 m ( f η ) = 3 ( a - 0.5 ) ( K m ref ) 3 k c + 6 ( a - 0.5 ) Y 1 m ( f η ) ( K m ref ) 2 k 1 - ak 2 - 2 a Y 1 m ( f η ) β ( K m ref ) 2 - 2 aq 3 β 3 ( a - 1 ) ( K m ref ) 4
其中, α = m 1 ( f η ) m 1 ( f η 0 ) , β = 1 m 1 2 ( f η ) [ m 2 ( f η 0 ) - m 1 ( f η 0 ) m 1 ( f η ) m 2 ( f η ) ] ,
m 1 ( f η ) = 1 cD ( f ηT ( f η ) ) - r T 0 k T 1 ( f ηT ref ( f η ) ) c D 2 ( f ηT ( f η ) ) + a T 1 cD ( f ηR ( f η ) ) - r R 0 k R 1 ( f ηR ref ( f η ) ) cD 2 ( f ηR ( f η ) )
m 2 ( f η ) = a T 2 2 cD ( f ηR ( f η ) ) - 2 k T 1 ( f ηT ref ( f η ) ) + r T 0 k T 2 ( f ηT ref ( f η ) ) 2 cD 2 ( f ηT ( f η ) ) + r T 0 k T 1 2 ( f ηT ref ( f η ) ) cD 3 ( f ηT ( f η ) ) - 2 a T 1 k R 1 ( f ηR ref ( f η ) ) + r R 0 k R 2 ( f ηR ref ( f η ) ) 2 cD 2 ( f ηR ( f η ) ) + r R 0 k R 1 2 ( f ηR ref ( f η ) ) c D 3 ( f ηR ( f η ) )
m 1 ( f η 0 ) = m 1 ( f η ) | f η = f η 0 , m 2 ( f η 0 ) = m 2 ( f η ) | f η = f η 0 ,
其中,是D(fηT)对rT的一阶和二阶导数;是D(fηR)对rT的一阶和二阶导数;aT1、aT2为rR对rT的一阶和二阶导数,η0是预先设定的参考值,fη0为一个方位参考频率;
利用上述四次滤波器H1(f)对回波二维频谱进行滤波,并将滤波后的二维频谱变换到距离多普勒域,其表达式为:
Sfilter(τ,fη)=exp{jΦRD(τ,fη)}
其中,τ为距离时间变量,相位ΦRD(τ,fη)表达为式:
Φ RD ( τ , f η ) = π K m ( f η ; x ) [ τ - τ d ( f η ) ] 2 + 2 π 3 Y 1 m ( f η ) K m 3 ( f η ; x ) [ τ - τ d ( f η ) ] 3 + 2 π 4 Y 2 m ( f η ) K m 4 ( f η ; x ) [ τ - τ d ( f η ) ] 4 + Φ AC ( f η ; x ) + Φ AL ( f η ; x , y )
步骤三:距离向四阶非线性Chirp-Scaling处理;
将四次相位滤波后的距离多普勒域信号乘以四阶非线性CS因子,非线性CS因子表达式为:
Scs(τ,fη)=exp{jΦCS(τ,fη)}
其相位ΦCS(τ,fη)为:
Φ CS ( τ , f η ) = π q 2 ( f η ) [ τ - τ d ref ( f η ) ] 2 + 2 π 3 q 3 ( f η ) [ τ - τ f ref ( f η ) ] 3 + 2 π 4 q 4 ( f η ) [ τ - τ f ref ( f η ) ] 4
其中,q2(fη),q3(fη)和q4(fη)是二阶、三阶和四阶项系数,表达式为:
q 2 ( f η ) = K m ref ( a - 1 ) q 3 ( f η ) = k 1 2 ( a - 1 ) - K m ref a 2 β q 4 ( f η ) = - 1 3 k 2 a + ( a - 1 ) ( K m ref ) 3 k c 2 + ( a - 1 ) Y 1 m ( f η ) ( K m ref ) 2 k 1 - 2 3 Y 1 m ( f η ) aβ ( K m ref ) 3 - 2 3 q 3 aβ
步骤四:RCM校正,距离压缩和高阶项处理;
经过非线性Chirp-Scaling处理后距离多普勒信号,距离徙动、二次距离压缩和三阶项的距离向空变性已经去除,因此,距离徙动、二次距离压缩和三阶项直接在二维频域来处理,校正函数为:
S jz ( f , f η ) = exp { j Φ 2 D CS ( f , f η ) }
其相位为:
Φ 2 D CS ( f , f η ) = 2 π Δτ ref f - πf 2 aK m ref + 2 π [ Y 1 m ( f η ) ( K m ref ) 3 + q 3 ] 3 α 3 ( K m ref ) 3 f 3 + 2 π [ Y 2 m ( f η ) ( K rm ref ) 4 + q 4 ] 4 α 4 ( K m ref ) 4 f 4 + Φ AC ( f η ; x ) + Φ AL ( f η ; x , y )
步骤五:方位向压缩;
将上述RCM校正、距离压缩和高阶项处理后的数据变换到距离时域方位频域,接下来将进行方位压缩,方位压缩函数为:
SAC(τ,fη)=exp{j[ΦAC(fη;x)+ΦAL(fη;x,y)]}
完成方位压缩后,进行方位向傅里叶反变换得到时域聚焦图像,从而完成移不变双基地前视SAR的聚焦成像。
CN201310375801.1A 2013-08-26 2013-08-26 移不变双基地前视合成孔径雷达nlcs成像方法 Active CN103412309B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310375801.1A CN103412309B (zh) 2013-08-26 2013-08-26 移不变双基地前视合成孔径雷达nlcs成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310375801.1A CN103412309B (zh) 2013-08-26 2013-08-26 移不变双基地前视合成孔径雷达nlcs成像方法

Publications (2)

Publication Number Publication Date
CN103412309A CN103412309A (zh) 2013-11-27
CN103412309B true CN103412309B (zh) 2015-12-23

Family

ID=49605336

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310375801.1A Active CN103412309B (zh) 2013-08-26 2013-08-26 移不变双基地前视合成孔径雷达nlcs成像方法

Country Status (1)

Country Link
CN (1) CN103412309B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103728617B (zh) * 2014-01-13 2016-01-20 电子科技大学 双基地合成孔径雷达时域快速成像方法
CN103869314B (zh) * 2014-03-18 2016-05-25 电子科技大学 双飞移变双基地前视合成孔径雷达成像方法
CN105204020B (zh) * 2015-10-09 2017-07-14 电子科技大学 基于粒子群优化的移不变双基前视sar飞行模式设计方法
CN106249237B (zh) * 2016-07-19 2019-05-21 西安电子科技大学 一种曲线轨迹下大斜视sar频域成像方法
CN107102330A (zh) * 2017-04-24 2017-08-29 清华大学 一种面向星载地理参考条带sar的wnlcs成像方法
CN106990397B (zh) * 2017-06-07 2020-06-30 电子科技大学 双基地前视sar非系统距离徙动校正方法
CN109444882B (zh) * 2018-11-05 2020-10-23 杭州电子科技大学 基于变斜视椭圆波束同步模型的双站sar成像方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6614386B1 (en) * 2002-02-28 2003-09-02 Raytheon Company Bistatic radar system using transmitters in mid-earth orbit
CN102147469B (zh) * 2010-12-29 2012-11-07 电子科技大学 一种双基地前视合成孔径雷达成像方法
CN102890277B (zh) * 2012-09-19 2014-07-09 电子科技大学 移不变双基地合成孔径雷达距离徙动成像方法
CN102914775B (zh) * 2012-10-10 2014-04-02 北京理工大学 一种改进的双基地合成孔径雷达nlcs成像算法

Also Published As

Publication number Publication date
CN103412309A (zh) 2013-11-27

Similar Documents

Publication Publication Date Title
CN103412309B (zh) 移不变双基地前视合成孔径雷达nlcs成像方法
CN103983974B (zh) 双站调频连续波合成孔径雷达成像方法
CN111142105B (zh) 复杂运动目标isar成像方法
CN103278820B (zh) 临近空间慢速平台sar动目标检测方法及成像方法
CN104833972B (zh) 一种双基地调频连续波合成孔径雷达频率变标成像方法
CN102288963B (zh) 基于子孔径参数估计的双基地isar图像融合方法
CN103412310B (zh) 双基地前视合成孔径雷达地面动目标检测方法与成像方法
CN102426360B (zh) 空中微动旋转目标的二维isar成像方法
CN102288961B (zh) 一种合成孔径雷达非线性调频变标的成像方法
CN103454635B (zh) 基于高超声速飞行器平飞段的前斜视sar成像方法
CN108427115B (zh) 合成孔径雷达对运动目标参数的快速估计方法
CN102226841B (zh) 基于高阶多项式距离方程的同步轨道sar成像方法
CN105487074B (zh) 一种双基地合成孔径雷达数值距离多普勒成像方法
CN102778681A (zh) 固定发射站双基地前视合成孔径雷达成像方法
CN110988878A (zh) 一种基于rd算法的sar海浪成像仿真方法
CN102749621A (zh) 一种双基地合成孔径雷达频域成像方法
CN109031299B (zh) 低信噪比条件下基于相位差分的isar平动补偿方法
CN102928841A (zh) 一种基于级数反演的机载圆迹环扫sar成像方法
CN103885062B (zh) 双基前视sar动目标成像方法与动目标速度估计方法
CN102788978A (zh) 一种斜视星机双基地合成孔径雷达成像方法
CN109143236B (zh) 适用于复杂飞行轨迹的双基聚束sar大场景成像方法
CN103630899A (zh) 地面运动目标高分辨雷达压缩感知成像的方法
CN104931965A (zh) 一种固定发射站双基正侧视sar成像方法
Xu et al. A variable PRF imaging method for high squint diving SAR
CN111624602A (zh) 一种旋转FrFT域构建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