CN102890277A - 移不变双基地合成孔径雷达距离徙动成像方法 - Google Patents
移不变双基地合成孔径雷达距离徙动成像方法 Download PDFInfo
- Publication number
- CN102890277A CN102890277A CN2012103476472A CN201210347647A CN102890277A CN 102890277 A CN102890277 A CN 102890277A CN 2012103476472 A CN2012103476472 A CN 2012103476472A CN 201210347647 A CN201210347647 A CN 201210347647A CN 102890277 A CN102890277 A CN 102890277A
- Authority
- CN
- China
- Prior art keywords
- eta
- frequency
- tan
- receiving station
- phi
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种移不变双基地SAR距离徙动成像方法,本发明的方法基于广义Loffeld变换的双基地SAR二维频谱,对频谱相位沿接收站最短斜距进行线性展开,导出Stolt频率变换表达式,通过该变换完成了残余距离单元徙动矫正、残余二次距离压缩和残余方位压缩,实现残余相位的空域线性化和频域线性化,进而实现了移不变双基地SAR的精确聚焦。本发明的方法与现有移不变双基地SAR距离徙动算法相比,形式简单,精度较高,同时运算效率较高。
Description
技术领域
本发明属于雷达技术领域,特别涉及合成孔径雷达成像技术中的移不变双基地合成孔径雷达(Synthetic Aperture Radar,SAR)的成像方法。
背景技术
SAR是一种全天时、全天候的现代高分辨率微波遥感成像雷达,在军事侦察、地形测绘、植被分析、海洋及水文观测、环境及灾害监视、资源勘探以及地壳微变检测等领域,SAR发挥了越来越重要的作用。
双基地SAR由于收发分置而有着很多突出的优点,它能获取目标的非后向散射信息,具有作用距离远、隐蔽性和抗干扰性强等特点。另外,由于双基地SAR接收机不含大功率器件,其功耗低、体积小、重量轻,便于多种类型的飞机携带,造价较低。总之,双基地SAR作为一种空间对地观测的新手段,在诸多领域都有着广阔的发展空间。
相对于传统的成像算法,如距离多普勒算法和Chirp Scaling算法,距离徙动算法通常被认为是最精确的成像算法。针对移不变双基地SAR,在文献“I.Walterscheid,J.Ender,A.Brenner,and O.Loffeld:Bistatic SAR processing and experiments,IEEE Trans.Geosci.RemoteSens.,vol.44,no.10,pp.2710–2717,2006”和文献“X.Qiu,D.Hu,and C.Ding:An Omega-Kalgorithm with phase error compensation for bistatic SAR of a translational invariant case,IEEETrans.Geosci.Remote Sens.,vol.46,no.8,pp.2224–2232,2008”中,提出了一种距离徙动成像算法,但是该算法采用的点目标响应二维频谱是基于数值计算方法得到的,运算效率较低。在文献“B.Liu,T.Wang,Q.Wu,and Z.Bao:Bistatic SAR data focusing using an Omega-Kalgorithm based on method of series reversion,IEEE Trans.Geosci.Remote Sens.,vol.47,no.8,pp.2899-2912,2009”中,提出了一种解析的距离徙动成像算法,但是该方法采用的点目标响应二维频谱形式复杂,并且该方法是基于双基地距离和对二维频谱进行线性化,精度较低。
发明内容
本发明的目的是为了解决现有的移不变双基地SAR成像算法存在的上述问题,提出了一种移不变双基地合成孔径雷达距离徙动成像方法。
为了方便描述本发明的内容,首先对以下术语进行解释:
术语1:双基地SAR,双基地SAR是指系统发射站和接收站分置于不同平台上的SAR系统,其中至少有一个平台为运动平台,在概念上属于双基地雷达。
术语2:移不变双基地SAR,移不变双基地SAR是双基地SAR的一种,收发站平行飞行,且速度相同。
本发明的技术方案为:一种移不变双基地SAR距离徙动成像方法,具体包括如下步骤:
步骤一:对原始回波数据进行二维傅立叶变换:
这里的步骤一具体可以通过如下过程实现:
双基地距离和为Rb(η;x,y)=RT(η;x,y)+RR(η;x,y),其中,η为方位时间,RT(η;x,y),RR(η;x,y)分别为发射站和接收站的距离历程:
其中,θsT和θsR分别为发射站和接收站的斜视角,V为接收站平台速度,rT和rR分别为发射站和接收站的最短斜距, hT和hR分别表示发射站和接收站的飞行高度;
原始回波数据在距离频域、方位时域的表达式为:
其中,f为距离频率, Br为发射信号带宽,Kr为发射信号的调频斜率,f0为系统载频,c为光速;
再对上式S(f,η;x,y)进行方位向傅里叶变换,基于广义Loffeld变换,得到原始回波在二维频域的表达式为:
S2df(f,fη;x,y)=S0(f)cxp{-jΦG(f,fη;x,y)}
其中,fη为方位频率,
fηT(fη)和fηR(fη)分别为发射站和接收站的多普勒频率:
其中,fηcT,fηcR分别为发射站和接收站对应的多普勒质心;fηrT,fηrR分别为发射站和接收站对应的多普勒调频斜率;fηc和fηr为系统总的多普勒质心和多普勒调频斜率;
步骤二:选取参考点,对步骤一中得到的二维频域数据进行参考函数匹配,完成粗聚焦;
这里的步骤二具体可以通过如下过程实现:
设选取的参考点为(x0,y0),该点回波的二维频谱为:
S2df(f,fη;x0,y0)=S0(f)cxp{-jΦG(f,fη;x0,y0)}
其中,
根据步骤一中rT和rR的表达式,将rT用rR进行表示:
然后将rT(rR)在rR0处关于rR进行线性泰勒展开,得到:rT(rR)≈rT0+aRΔr,其中,rT0=rT(rR0),Δr=rR-rR0,
则有:
其中,Δy=y-y0。
步骤三:对步骤二中的匹配结果进行Stolt频率变换:
这里的Stolt频率变换可以通过如下过程实现:
所述的Stolt频率变换的具体通过表达式aRFT(f,fη)+FR(f,fη)=f'+f0完成,其中,f′为变换后的距离频率,完成该变换后,则有:
变换后结果为:
步骤四:对步骤三中频率变换的结果SRES(f',fη)进行二维傅立叶反变换,得到聚焦SAR图像,变换结果为:Simage(rR,y)≈sinc(rR-Δr)sinc(y-Δy)。
本发明的有益效果:本发明的方法基于广义Loffeld变换的双基地SAR二维频谱,对频谱相位沿接收站最短斜距进行线性展开,导出Stolt频率变换表达式,通过该变换完成了残余距离单元徙动矫正、残余二次距离压缩和残余方位压缩,实现残余相位的空域线性化和频域线性化,进而实现了移不变双基地SAR的精确聚焦。本发明的方法与现有移不变双基地SAR距离徙动算法相比,形式简单,精度较高,同时运算效率较高,能够满足移不变双基地SAR成像处理的要求,可以应用于地球遥感等领域。
附图说明
图1是本发明具体实施例采用的移不变双基地SAR系统结构图。
图2是本发明具体实施例采用的移不变双基地SAR系统参数表。
图3是本发明方法的流程示意图。
图4是本发明具体实施例中采用的目标场景布置图,其中,黑色圆点为布置于地面上的共5个点目标,这5个点沿x方向(切航迹)间隔500米,沿y方向(沿航迹)间隔200米,平台沿y轴运动。
图5是本发明具体实施例中回波二维频谱示意图。
图6是本发明具体实施例中经过Stolt变换后的二维频谱示意图。
图7是本发明具体实施例中对图4中5个点目标进行成像的结果示意图。
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2010上验证正确。下面就具体实施方式对本发明作进一步的详细描述。
本发明具体实施方式采用的固定站双基地SAR系统结构图如图1所示,在直角坐标系中设发射平台位置为(xT ,yT,hT),接收站零时刻位置记为(xR,yR,hR),接收站速度记为υ,任意成像点坐标记为P(x,y),系统坐标系以成像中心点目标O位坐标原点,平台沿y轴运动,x轴为切航迹方向,z轴为垂直地面方向。在介绍本发明方法之前,先对系统参数进行成像参数初始化,生成目标回波,采用图2中给出的系统参数,对图4中的点目标逐点生成回波数据,累加得到目标回波,记为S(τ,η)。本发明方法的具体流程如图3所示,具体步骤如下:
步骤一:对回波数据S(τ,η)分别进行方位向和距离向FFT(Fast FourierTransformation),得到移不变双基地SAR二维频谱,记为S2df(f,fη),其中,τ为距离时间。
步骤二:选取成像区域中心点即O点为参考点,其二维频谱记为S2df(f,fη;x=0,y=0)。利用参考点二维频谱对移不变双基地SAR进行粗匹配聚焦,即用步骤二得到的回波S2df(f,fη)与参考点二维频谱S2df(f,fη;x=0,y=0)共轭相乘:得到匹配后的数据,记为SRES(f,fη)。回波二维频谱示意图如图5所示。
步骤三:利用8点sinc插值,对步骤二得到的二维频域数据SRES(f,fη)进行aRFT(f,fη)+FR(f,fη)=f'+f0距离向频率stolt变换,从而实现了距离频率的线性化,得到变换后的数据,记为SRES(f′,fη)。经过Stolt变换后的二维频谱示意图如图6所示。
由于距离向频率stolt变换为非线性变换,所以需要进行一维插值来实现,该频率变换可以完成残余距离单元徙动矫正,残余二次距离压缩和残余方位压缩。
步骤四:将步骤三得到频率变换后SRES(f,fη)进行二维IFFT操作,得到聚焦SAR图像,图像记为Simage(rR,y)。
至此,完成移不变双基地斜视合成孔径雷达频域成像处理,成像结果如图7所示。通过本发明具体实施方式可以看出,本发明可以实现对移不变双基地SAR回波的精确聚焦。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (4)
1.一种移不变双基地SAR距离徙动成像方法,具体包括如下步骤:
步骤一:对原始回波数据进行二维傅立叶变换:
步骤二:选取参考点,对步骤一中得到的二维频域数据进行参考函数匹配,完成粗聚焦;
步骤三:对步骤二中的匹配结果进行Stolt频率变换:
步骤四:对步骤三中频率变换的结果进行二维傅立叶反变换,得到聚焦SAR图像。
2.根据权利要求1所述的移不变双基地SAR距离徙动成像方法,其特征在于,步骤一的具体过程如下:
双基地距离和为Rb(η;x,y)=RT(η;x,y)+RR(η;x,y),其中,η为方位时间,RT(η;x,y),RR(η;x,y)分别为发射站和接收站的距离历程:
其中,θsT和θsR分别为发射站和接收站的斜视角,V为接收站平台速度,rT和rR分别为发射站和接收站的最短斜距, hT和hR分别表示发射站和接收站的飞行高度;
原始回波数据在距离频域、方位时域的表达式为:
再对上式S(f,η;x,y)进行方位向傅里叶变换,基于广义Loffeld变换,得到原始回波在二维频域的表达式为:
S2df(f,fη;x,y)=S0(f)cxp{-jΦG(f,fη;x,y)}
其中,fη为方位频率,
fηT(fη)和fηR(fη)分别为发射站和接收站的多普勒频率:
其中,fηcT,fηcR分别为发射站和接收站对应的多普勒质心;fηrT,fηrR分别为发射站和接收站对应的多普勒调频斜率;fηc和fηr为系统总的多普勒质心和多普勒调频斜率。
3.根据权利要求2所述的移不变双基地SAR距离徙动成像方法,其特征在于,步骤二的具体过程如下:设选取的参考点为(x0,y0),该点回波的二维频谱为:
S2df(f,fη-x0,y0)=S0(f)cxp{-jΦG(f,fη;x0,y0)}
其中,
根据步骤一中rT和rR的表达式,将rT用rR进行表示:
然后将rT(rR)在rR0处关于rR进行线性泰勒展开,得到:rT(rR)≈rT0+aRΔr,其中,rT0=rT(rR0),Δr=rR-rR0,
则有:
其中,Δy=y-y0。
4.根据权利要求2或3所述的的移不变双基地SAR距离徙动成像方法,其特征在于,步骤三所述的Stolt频率变换的具体过程如下:
所述的Stolt频率变换的具体通过表达式aRFT(f,fη)+FR(f,fη)=f'+f0完成,其中,f′为变换后的距离频率,完成该变换后,则有:
变换后结果为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210347647.2A CN102890277B (zh) | 2012-09-19 | 2012-09-19 | 移不变双基地合成孔径雷达距离徙动成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210347647.2A CN102890277B (zh) | 2012-09-19 | 2012-09-19 | 移不变双基地合成孔径雷达距离徙动成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102890277A true CN102890277A (zh) | 2013-01-23 |
CN102890277B CN102890277B (zh) | 2014-07-09 |
Family
ID=47533832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210347647.2A Active CN102890277B (zh) | 2012-09-19 | 2012-09-19 | 移不变双基地合成孔径雷达距离徙动成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102890277B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412309A (zh) * | 2013-08-26 | 2013-11-27 | 电子科技大学 | 移不变双基地前视合成孔径雷达nlcs成像方法 |
CN104749570A (zh) * | 2015-04-01 | 2015-07-01 | 电子科技大学 | 一种移不变机载双基合成孔径雷达目标定位方法 |
CN105204020A (zh) * | 2015-10-09 | 2015-12-30 | 电子科技大学 | 基于粒子群优化的移不变双基前视sar飞行模式设计方法 |
CN108490443A (zh) * | 2018-02-24 | 2018-09-04 | 中国人民解放军海军工程大学 | 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法 |
-
2012
- 2012-09-19 CN CN201210347647.2A patent/CN102890277B/zh active Active
Non-Patent Citations (5)
Title |
---|
《电子科技大学学报》 20071031 武俊杰 等 双基地斜视SAR距离徙动成像算法 1106-1109 2 第36卷, 第5期 * |
INGO WALTERSCHEID 等: "Bistatic SAR Processing and Experiments", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
VINCENT GIROUX 等: "An Omega-K algorithm for SAR bistatic systems", 《IEEE》 * |
武俊杰 等: "双基地斜视SAR距离徙动成像算法", 《电子科技大学学报》 * |
武俊杰: "双基地斜视合成孔径雷达成像技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412309A (zh) * | 2013-08-26 | 2013-11-27 | 电子科技大学 | 移不变双基地前视合成孔径雷达nlcs成像方法 |
CN104749570A (zh) * | 2015-04-01 | 2015-07-01 | 电子科技大学 | 一种移不变机载双基合成孔径雷达目标定位方法 |
CN104749570B (zh) * | 2015-04-01 | 2018-06-01 | 电子科技大学 | 一种移不变机载双基合成孔径雷达目标定位方法 |
CN105204020A (zh) * | 2015-10-09 | 2015-12-30 | 电子科技大学 | 基于粒子群优化的移不变双基前视sar飞行模式设计方法 |
CN105204020B (zh) * | 2015-10-09 | 2017-07-14 | 电子科技大学 | 基于粒子群优化的移不变双基前视sar飞行模式设计方法 |
CN108490443A (zh) * | 2018-02-24 | 2018-09-04 | 中国人民解放军海军工程大学 | 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法 |
Also Published As
Publication number | Publication date |
---|---|
CN102890277B (zh) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
Wang et al. | Bistatic SAR system and signal processing technology | |
CN102147469B (zh) | 一种双基地前视合成孔径雷达成像方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN102778681B (zh) | 固定发射站双基地前视合成孔径雷达成像方法 | |
CN101587500B (zh) | 双站合成孔径雷达海面成像的计算机仿真方法 | |
CN103412310B (zh) | 双基地前视合成孔径雷达地面动目标检测方法与成像方法 | |
CN101887122B (zh) | 可消除地面高程误差的星载sar图像目标定位方法 | |
CN103983974B (zh) | 双站调频连续波合成孔径雷达成像方法 | |
CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN104330795B (zh) | 一种基于Keystone变换的地基合成孔径雷达快速成像方法 | |
CN102004250B (zh) | 基于频域展开的星机联合双基地合成孔径雷达成像方法 | |
CN105487074B (zh) | 一种双基地合成孔径雷达数值距离多普勒成像方法 | |
CN102914773B (zh) | 一种多航过圆周sar三维成像方法 | |
CN104749570A (zh) | 一种移不变机载双基合成孔径雷达目标定位方法 | |
CN104251990A (zh) | 合成孔径雷达自聚焦方法 | |
CN104898119A (zh) | 一种基于相关函数的动目标参数估计方法 | |
CN103869314A (zh) | 双飞移变双基地前视合成孔径雷达成像方法 | |
CN102967862B (zh) | 双飞移变模式双基地合成孔径雷达成像方法 | |
CN105223572A (zh) | 一种基于pfa算法的正前视双基sar成像处理方法 | |
CN105158745A (zh) | 移变双基前视合成孔径雷达距离徙动校正方法 | |
CN102890277B (zh) | 移不变双基地合成孔径雷达距离徙动成像方法 | |
CN103412309A (zh) | 移不变双基地前视合成孔径雷达nlcs成像方法 | |
CN103630903B (zh) | 基于顺轨干涉sar测量海面流场径向速度的方法 | |
CN103630899B (zh) | 地面运动目标高分辨雷达压缩感知成像的方法 |
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 |