CN104991251B - 基于匀加速建模的超高分辨率星载sar成像方法 - Google Patents
基于匀加速建模的超高分辨率星载sar成像方法 Download PDFInfo
- Publication number
- CN104991251B CN104991251B CN201510363102.4A CN201510363102A CN104991251B CN 104991251 B CN104991251 B CN 104991251B CN 201510363102 A CN201510363102 A CN 201510363102A CN 104991251 B CN104991251 B CN 104991251B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- orientation
- 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
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
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于匀加速建模的超高分辨率星载SAR成像方法,包括以下步骤:(1)将卫星的等效速度沿方位向的变化建模为匀加速运动;(2)接收星载SAR回波信号,并对其进行去载频处理,得到去载频后的回波信号建立点目标(r,tc)的精确斜距历程R(tm);(3)对去载频后的回波信号进行距离傅立叶变换,并进行距离压缩,得到距离压缩后的回波信号S(fr,tm);(4)将距离压缩后的回波信号S(fr,tm)依次进行方位向插值、方位快速傅立叶变换和方位向相位误差补偿,得到相位误差补偿后的回波信号S(fr,fa);(5)消除相位误差补偿后的回波信号S(fr,fa)中的卫星的等效速度沿距离向的变化,得到最终的距离‑多普勒域回波信号并对其进行方位逆傅立叶变换,得到星载SAR图像。
Description
技术领域
本发明属于雷达信号处理领域,涉及星载合成孔径雷达(SAR)成像技术,特别涉及一种基于匀加速建模的超高分辨率星载SAR成像方法,可用于星载平台的超高分辨SAR成像。
背景技术
随着合成孔径雷达在对地观测方面的应用越来越广泛,大场景高分辨率成像已在国内外引起了高度的关注。距离高分辨率星载SAR技术近年来得到了迅猛地发展,国际上一些先进的距离高分辨率星载SAR能达到优于1米的分辨率,但是方位高分辨率星载SAR技术仍有待发展。由于卫星轨道以及地球旋转分量的时变性,卫星的等效速度不仅沿距离向变化,还沿方位向缓慢改变。在低分辨率或小场景的情况下,只需要考虑卫星的等效速度沿距离向的变化,但在超高分辨率及宽方位测绘带的情况下,卫星的等效速度沿方位向的变化也要予以考虑。
对星载SAR进行高分辨成像需要对斜距历程进行精确地建模,但由于星载SAR的几何关系较为复杂,现有星载SAR的斜距历程模型大都为近似模型,这些模型都没有考虑等效速度为非匀速的性质。传统的双曲线斜距模型具有便于推导成像算法这一优势,但由于它没有考虑等效速度沿方位向的变化,在高分辨率的情况下,尤其是超高分辨的情况下传统的双曲线斜距模型的近似误差较大,会导致点目标在方位向上产生散焦;现有的一些算法中,比较典型的算法是基于距离徙动算法(RMA)的大场景超高分辨星载SAR成像方法,该方法针对卫星等效速度沿距离向的变化问题对距离徙动算法进行了改进,能够很好地处理卫星等效速度沿距离向的变化对聚焦的影响,但是无法处理方位测绘带较宽时卫星等效速度沿方位向的变化问题对聚焦的影响。
发明内容
针对上述现有技术的不足,本发明的目的在于提出了一种基于匀加速建模的超高分辨率星载SAR成像方法,该方法首先将卫星的等效速度沿方位向的变化建模为匀加速运动,然后对方位向回波信号进行插值,以克服卫星的等效速度沿方位向的变化,最后运用距离徙动算法,得到超高分辨率的星载SAR图像;本发明能够解决宽方位测绘带时卫星等效速度沿方位向的变化导致的点目标在方位向散焦的问题。
为实现上述技术目的,本发明采用以下技术方案予以实现。
一种基于匀加速建模的超高分辨率星载SAR成像方法,其特征在于,包括以下步骤:
步骤1,将卫星的等效速度沿方位向的变化建模为匀加速运动,并将点目标(r,tc)的精确斜距历程R(tm)表示为近似斜距历程R′(tm)与斜距历程误差ΔR(tm)的和,r为卫星距点目标(r,tc)的最短斜距,tc为点目标(r,tc)的方位零多普勒时刻,tm为慢时间;
步骤2,接收星载SAR回波信号,将星载SAR回波信号进行去载频处理,得到去载频后的回波信号其中,为快时间;
步骤3,对去载频后的回波信号进行距离快速傅立叶变换(FFT),并进行距离压缩,得到距离压缩后的回波信号S(fr,tm),其中,fr为距离频率;
步骤4,将距离压缩后的回波信号S(fr,tm)进行方位插值,得到去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′);对去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′)进行方位快速傅立叶变换(FFT),得到去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa);对去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)进行方位向相位误差补偿,得到方位向相位误差补偿后的回波信号S(fr,fa);其中,tm′为方位插值慢时间,fa为多普勒频率;
步骤5,运用距离徙动算法(RMA)消除方位向相位误差补偿后的回波信号S(fr,fa)中的卫星的等效速度沿距离向的变化,得到最终的距离-多普勒域回波信号对最终的距离-多普勒域回波信号进行方位逆傅立叶变换,得到星载SAR图像。
本发明的有益效果为:针对现有算法只克服了卫星的等效速度沿距离向的变化,且在超高分辨率宽方位测绘带的情况下,点目标在方位向上散焦严重的问题,本发明提出的基于匀加速建模的超高分辨率星载SAR成像方法,能够较好地处理卫星的等效速度沿方位向的变化对聚焦的影响,有效提升了超高分辨率星载SAR的成像质量。
附图说明
下面结合附图说明和具体实施方式对本发明作进一步详细说明。
图1是本发明的流程图。
图2是仿真实验的布点示意图,横坐标为方位向长度,单位为千米(Km),纵坐标为距离向长度,单位为千米(Km)。
图3是使用传统双曲线斜距模型进行点目标成像的结果图,横坐标为方位向采样单元,纵坐标为距离向采样单元,其中:
图3a为点目标1的成像结果图;
图3b为点目标5的成像结果图;
图3c为点目标9的成像结果图。
图4是使用本发明方法进行点目标成像的结果图,横坐标为方位向采样单元,纵坐标为距离向采样单元,其中:
图4a为点目标1的成像结果图;
图4b为点目标5的成像结果图;
图4c为点目标9的成像结果图。
具体实施方式
参照图1,本发明的基于匀加速建模的超高分辨率星载SAR成像方法,包括以下步骤:
步骤1,将卫星的等效速度沿方位向的变化建模为匀加速运动,并将点目标(r,tc)的精确斜距历程R(tm)表示为近似斜距历程R′(tm)与斜距历程误差ΔR(tm)的和,点目标(r,tc)的精确斜距历程R(tm)的表达式为:
R(tm)=R′(tm)+ΔR(tm)
其中,tm为慢时间,r为卫星距点目标(r,tc)的最短斜距,tc为点目标(r,tc)的方位零多普勒时刻,R′(tm)为近似斜距历程,其表达式为:
其中,为参考目标(r,tref)的方位向等效速度,tref为选定的参考时间,αr为点目标(r,tc)的方位向等效加速度;
ΔR(tm)为斜距历程误差,其表达式为:
ΔR(tm)=β(tm-tc)3+γ(tm-tc)4
其中,和分别为斜距历程误差ΔR(tm)的三次项系数和四次项系数,β和γ在整个场景中的变化很小,认定β和γ为恒定不变;为点目标(r,tc)的方位向等效速度;c3=-12fd3/λ,c4=-48fd4/λ,fd3为方位调频率的变化率,fd4为方位调频率的变化率fd3的变化率,fd3和fd4的表达式分别为:
其中,Rsat、Vsat以及Asat分别为卫星的位置、速度以及加速度矢量;A′sat和A″sat分别为卫星加速度矢量的导数和二次导数;Rtar、vtar以及Atar分别为地面目标的位置、速度以及加速度矢量;A′tar和A″tar分别表示地面目标的加速度矢量的导数和二次导数;fd为多普勒中心频率,fd2为方位调频率,fd和fd2的表达式分别为:
步骤2,接收星载SAR回波信号,将星载SAR回波信号进行去载频处理,得到去载频后的回波信号为:
其中,为快时间,c为光速,λ为星载SAR的中心波长,fr为距离频率,Wr(·)和Wa(·)分别表示星载SAR的距离窗函数和方位窗函数,R(tm)为点目标(r,tc)的精确斜距历程。
步骤3,对去载频后的回波信号进行距离快速傅立叶变换(FFT),并构造距离压缩函数进行距离压缩,κ为距离向调频率,得到距离压缩后的回波信号S(fr,tm)为:
其中,fc为星载SAR发射信号载频。
步骤4,将距离压缩后的回波信号S(fr,tm)进行方位向插值,得到去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′);对去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′)进行方位快速傅立叶变换(FFT),得到去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa);对去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)进行方位向相位误差补偿,得到方位向相位误差补偿后的回波信号S(fr,fa)。
步骤4的具体子步骤为:
4.1将距离压缩后的回波信号S(fr,tm)进行方位向插值,即将慢时间tm表示为方位向插值慢时间tm′, 为插值系数,为场景中心(rS,tc)的方位向等效加速度,为场景中心的方位向等效速度,rS为场景中心的最短斜距;再将点目标(r,tc)的方位零多普勒时刻tc表示为点目标(r,tc)的方位向插值零多普勒时刻tc′,得到去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′)为:
其中,R(tm′)为将慢时间tm进行方位向插值后的精确斜距历程,其表达式为:
4.2对去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′)进行方位快速傅立叶变换(FFT),得到去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)为:
其中,A0为星载SAR回波信号幅度,和分别为去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)的三次相位误差和四次相位误差,和的表达式分别为:
4.3构造方位相位误差补偿函数H2,对去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)的三次相位误差和四次相位误差进行方位向相位误差补偿,方位相位误差补偿函数H2为:
得到方位向相位误差补偿后的回波信号S(fr,fa)为:
步骤5,运用距离徙动算法(RMA)消除方位向相位误差补偿后的回波信号S(fr,fa)中的卫星的等效速度沿距离向的变化,得到最终的距离-多普勒域回波信号对最终的距离-多普勒域回波信号进行方位逆傅立叶变换,得到星载SAR图像。
步骤5的具体子步骤为:
5.1构造距离徙动相位误差补偿函数H3(fr,fa)为:
对方位向相位误差补偿后的回波信号S(fr,fa)进行距离徙动相位误差补偿,得到剩余相位误差为
5.2对剩余相位误差进行进行Stolt变量替换,即令为Stolt变量替换后的距离频率,得到变量替换后的剩余相位误差为:
5.3对变量替换后的剩余相位误差进行二阶泰勒因式分解,则其中,φ0为变量替换后的剩余相位误差的常数项,φ1为变量替换后的剩余相位误差的一次项系数,φ0和φ1的表达式分别为:
其中,为场景中心点的最大多普勒频率,faM;r为点目标(r,tc)的最大多普勒频率,
5.4根据变量替换后的剩余相位误差的二阶泰勒分解因式将方位向相位误差补偿后的回波信号S(fr,fa)变换到距离-多普勒域,得到距离-多普勒域回波信号为:
其中,fd为多普勒中心频率;
在距离-多普勒域,根据插值公式对距离-多普勒域回波信号进行插值处理;并将距离-多普勒域回波信号乘以剩余相位误差补偿函数补偿变量替换后的剩余相位误差的常数项φ0,得到最终的距离-多普勒域回波信号为:
其中,B为发射信号带宽,为距离插值快时间;
5.5对最终的距离-多普勒域回波信号进行方位逆傅立叶变换,得到星载SAR图像。
本发明的有效性可通过以下仿真和实测数据作进一步说明。
1)点目标仿真数据
在10千米×10千米的仿真场景中,在方位向和距离向放置9个点目标,如图2所示;仿真参数如表1所示;
表1仿真参数
2)仿真内容和结果分析
分别用传统的双曲线斜距模型和本发明方法对图2中的点目标进行仿真成像。图3是使用传统双曲线斜距模型进行点目标成像的结果图,其中:图3a为点目标1的成像结果图;图3b为点目标5的成像结果图;图3c为点目标9的成像结果图。图4是使用本发明方法进行点目标成像的结果图,其中:图4a为点目标1的成像结果图;图4b为点目标5的成像结果图;图4c为点目标9的成像结果图。
从图3a-图3c中可以看出,运用传统的双曲线斜距模型对点目标成像时,由于点目标5处于方位向中间,所以点目标5得到了很好的聚焦,而点目标1和点目标9由于位于方位向的两边,点目标1和点目标9的方位向等效速度不同于中心点目标5的方位向等效速度,造成点目标1和点目标9沿方位向散焦非常严重。
从图4a-图4c中可以看出,运用本发明方法对点目标成像时,三个点目标都得到了很好的聚焦,说明本发明方法能够很好地解决卫星的等效速度沿方位向的变化对点目标聚焦效果的影响,说明了本发明的有效性。
分别用传统的双曲线斜距模型(简称传统模型)和本发明方法对图2中的点目标进行仿真成像的成像质量指标如表2所示,其中ρa为方位向分辨率,PSLR为峰值旁瓣比,ISLR为积分旁瓣比,ρr为距离向分辨率。
表2点目标的成像质量指标
从表2中可以看出,使用本发明方法得到的所有的点目标的方位向峰值旁瓣比均约为-13.2,方位向积分旁瓣比均约为-10.5,距离向峰值旁瓣比均约为-13.3,距离向积分旁瓣比均约为-9.95,上述指标值非常接近理想点目标的指标值,远远优于传统的双曲线斜距模型的成像指标;从表2中可以看出,本发明对点目标的成像结果的分辨率也要优于传统的双曲线斜距模型,进一步证明了本发明的有效性。
Claims (3)
1.一种基于匀加速建模的超高分辨率星载SAR成像方法,其特征在于,包括以下步骤:
步骤1,将卫星的等效速度沿方位向的变化建模为匀加速运动,并将点目标(r,tc)的精确斜距历程R(tm)表示为近似斜距历程R′(tm)与斜距历程误差ΔR(tm)的和,r为卫星距点目标(r,tc)的最短斜距,tc为点目标(r,tc)的方位零多普勒时刻,tm为慢时间;
其中,所述点目标(r,tc)的精确斜距历程R(tm),其表达式为:
R(tm)=R′(tm)+ΔR(tm)
其中,tm为慢时间,r为卫星距点目标(r,tc)的最短斜距,tc为点目标(r,tc)的方位零多普勒时刻,R′(tm)为近似斜距历程,其表达式为:
<mrow>
<msup>
<mi>R</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&alpha;</mi>
<mi>r</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&alpha;</mi>
<mi>r</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
其中,为参考目标(r,tref)的方位向等效速度,tref为选定的参考时间,αr为点目标(r,tc)的方位向等效加速度;
ΔR(tm)为斜距历程误差,其表达式为:
ΔR(tm)=β(tm-tc)3+γ(tm-tc)4
其中,和分别为斜距历程误差ΔR(tm)的三次项系数和四次项系数,β和γ在整个场景中的变化很小,认定β和γ为恒定不变;为点目标(r,tc)的方位向等效速度;c3=-12fd3/λ,c4=-48fd4/λ,λ为星载SAR的中心波长,fd3为方位调频率的变化率,fd4为方位调频率的变化率fd3的变化率,fd3和fd4的表达式分别为:
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>6</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msub>
<mi>&lambda;f</mi>
<mi>d</mi>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>8</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>6</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3</mn>
<msubsup>
<mi>&lambda;f</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&lambda;f</mi>
<mi>d</mi>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mn>3</mn>
</mrow>
</msub>
</mrow>
<mi>r</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Rsat、Vsat以及Asat分别为卫星的位置、速度以及加速度矢量;A′sat和A″sat分别为卫星加速度矢量的导数和二次导数;Rtar、Vtar以及Atar分别为地面目标的位置、速度以及加速度矢量;A′tar和A″tar分别表示地面目标的加速度矢量的导数和二次导数;fd为多普勒中心频率,fd2为方位调频率,fd和fd2的表达式分别为:
<mrow>
<msub>
<mi>f</mi>
<mi>d</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&lambda;</mi>
<mi>r</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&lambda;f</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<mi>r</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
步骤2,接收星载SAR回波信号,将星载SAR回波信号进行去载频处理,得到去载频后的回波信号其中,为快时间;
步骤3,对去载频后的回波信号进行距离快速傅立叶变换,并进行距离压缩,得到距离压缩后的回波信号S(fr,tm),其中,fr为距离频率;
步骤4,将距离压缩后的回波信号S(fr,tm)进行方位向插值,得到去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′);对去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′)进行方位快速傅立叶变换,得到去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa);对去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)进行方位向相位误差补偿,得到方位向相位误差补偿后的回波信号S(fr,fa);其中,tm′为方位插值慢时间,fa为多普勒频率;
步骤5,消除方位向相位误差补偿后的回波信号S(fr,fa)中的卫星的等效速度沿距离向的变化,得到最终的距离-多普勒域回波信号对最终的距离-多普勒域回波信号进行方位逆傅立叶变换,得到星载SAR图像。
2.如权利要求1所述的基于匀加速建模的超高分辨率星载SAR成像方法,其特征在于,步骤4的具体子步骤为:
4.1将距离压缩后的回波信号S(fr,tm)进行方位向插值,即将慢时间tm表示为方位向插值慢时间tm′,为插值系数,为场景中心(rS,tc)的方位向等效加速度,为场景中心的方位向等效速度,rS为场景中心的最短斜距,tref为选定的参考时间;再将点目标(r,tc)的方位零多普勒时刻tc表示为点目标(r,tc)的方位向插值零多普勒时刻tc′,得到去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′)为:
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>W</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mi>c</mi>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
<mi>R</mi>
<mo>(</mo>
<mrow>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
其中,Wr(·)和Wa(·)分别表示星载SAR的距离窗函数和方位窗函数,fc为星载SAR发射信号载频,c为光速;R(tm′)为将慢时间tm进行方位向插值后的精确斜距历程,表达式为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>&beta;</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>4</mn>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>4</mn>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>&gamma;</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>4</mn>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>4</mn>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&ap;</mo>
<msqrt>
<mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<mi>&beta;</mi>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>+</mo>
<mi>&gamma;</mi>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>t</mi>
<mi>m</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,为参考目标(r,tref)的方位向等效速度,β和γ分别为斜距历程误差ΔR(tm)的三次项系数和四次项系数;
4.2对去除卫星的等效速度沿方位向的变化的回波信号S(fr,tm′)进行方位快速傅立叶变换,得到去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)为:
其中,A0为星载SAR回波信号幅度,和分别为去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)的三次相位误差和四次相位误差,和的表达式分别为:
4.3构造方位相位误差补偿函数H2,对去除卫星的等效速度沿方位向的变化的二维频域回波信号Sech(fr,fa)的三次相位误差和四次相位误差进行方位向相位误差补偿,方位相位误差补偿函数H2为:
得到方位向相位误差补偿后的回波信号S(fr,fa)为:
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<msub>
<mi>W</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>W</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>j</mi>
<mn>4</mn>
<mi>&pi;</mi>
<mi>r</mi>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>f</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>a</mi>
</msub>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3.如权利要求1所述的基于匀加速建模的超高分辨率星载SAR成像方法,其特征在于,步骤5的具体子步骤为:
5.1构造距离徙动相位误差补偿函数H3(fr,fa)为:
<mrow>
<msub>
<mi>H</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>j</mi>
<mn>4</mn>
<msub>
<mi>&pi;r</mi>
<mi>S</mi>
</msub>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>f</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>v</mi>
<mrow>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,为场景中心的方位向等效速度,rS为场景中心的最短斜距,tref为选定的参考时间,c为光速;
对方位向相位误差补偿后的回波信号S(fr,fa)进行距离徙动相位误差补偿,得到剩余相位误差为
其中,fc为星载SAR发射信号载频;
5.2对剩余相位误差进行进行Stolt变量替换,即令fr′为Stolt变量替换后的距离频率,得到变量替换后的剩余相位误差为:
5.3对变量替换后的剩余相位误差进行二阶泰勒因式分解,则其中,φ0为变量替换后的剩余相位误差的常数项,φ1为变量替换后的剩余相位误差的一次项系数,φ0和φ1的表达式分别为:
<mrow>
<msub>
<mi>&phi;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mi>r</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>a</mi>
<mi>M</mi>
<mo>;</mo>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>a</mi>
<mi>M</mi>
<mo>;</mo>
<mi>r</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&phi;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>-</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>r</mi>
</mrow>
<mi>c</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>a</mi>
<mi>M</mi>
<mo>;</mo>
<msub>
<mi>r</mi>
<mi>S</mi>
</msub>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>a</mi>
</msub>
<msub>
<mi>f</mi>
<mrow>
<mi>a</mi>
<mi>M</mi>
<mo>;</mo>
<mi>r</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为场景中心点的最大多普勒频率,faM;r为点目标(r,tc)的最大多普勒频率, 为参考目标(r,tref)的方位向等效速度,λ为星载SAR的中心波长;
5.4根据变量替换后的剩余相位误差的二阶泰勒分解因式将方位向相位误差补偿后的回波信号S(fr,fa)变换到距离-多普勒域,得到距离-多普勒域回波信号为:
其中,A0为星载SAR回波信号幅度,Wr(·)和Wa(·)分别表示星载SAR的距离窗函数和方位窗函数,B为发射信号带宽,fd为多普勒中心频率;
在距离-多普勒域,根据插值公式对距离-多普勒域回波信号进行插值处理,为距离插值快时间;并将距离-多普勒域回波信号乘以剩余相位误差补偿函数补偿变量替换后的剩余相位误差的常数项φ0,得到最终的距离-多普勒域回波信号为:
5.5对最终的距离-多普勒域回波信号进行方位逆傅立叶变换,得到星载SAR图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510363102.4A CN104991251B (zh) | 2015-06-26 | 2015-06-26 | 基于匀加速建模的超高分辨率星载sar成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510363102.4A CN104991251B (zh) | 2015-06-26 | 2015-06-26 | 基于匀加速建模的超高分辨率星载sar成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104991251A CN104991251A (zh) | 2015-10-21 |
CN104991251B true CN104991251B (zh) | 2017-08-25 |
Family
ID=54303087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510363102.4A Active CN104991251B (zh) | 2015-06-26 | 2015-06-26 | 基于匀加速建模的超高分辨率星载sar成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104991251B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443670B (zh) * | 2016-08-29 | 2018-10-19 | 西北工业大学 | 一种非匀速大斜视cs成像方法 |
CN106872980B (zh) * | 2017-04-12 | 2019-10-18 | 中国科学院电子学研究所 | 基于匀加速曲线运动模型的收发分离斜距确定方法 |
CN109270507B (zh) * | 2018-11-05 | 2023-08-15 | 金华航大北斗应用技术有限公司 | Gnss-r双基sar成像机载运动误差补偿模型 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102759734B (zh) * | 2012-05-23 | 2013-09-11 | 北京理工大学 | 一种用于高分辨率星载sar的改进距离徙动校正算法的成像方法 |
CN103323842A (zh) * | 2012-09-03 | 2013-09-25 | 中国科学院电子学研究所 | 一种调频连续波合成孔径雷达中的成像方法及装置 |
CN103323828B (zh) * | 2013-05-29 | 2015-08-19 | 中国科学院电子学研究所 | 超高分辨率星载sar成像处理方法及装置 |
-
2015
- 2015-06-26 CN CN201510363102.4A patent/CN104991251B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104991251A (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104931967B (zh) | 一种改进的高分辨率sar成像自聚焦方法 | |
CN108490441B (zh) | 基于两级滤波的俯冲段大斜视sar子孔径成像空变校正方法 | |
CN102176016B (zh) | 一种大斜视滑动聚束sar成像处理方法 | |
CN103885061B (zh) | 基于改进相位梯度自聚焦的双基sar运动补偿方法 | |
CN105759263B (zh) | 一种高分辨率大场景下的星载斜视sar雷达成像方法 | |
CN103543453B (zh) | 一种地球同步轨道合成孔径雷达干涉的高程反演方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN106405552B (zh) | 基于wvd—pga算法的sar雷达目标聚焦方法 | |
CN107918124A (zh) | 带有方位空变校正的机载大斜视高分辨sar成像方法 | |
CN109655830B (zh) | 一种针对星上去走动与二维空变校正的地球同步轨道sar斜视成像处理装置 | |
CN110596701B (zh) | 基于二次椭圆模型的非平飞双站sar频域fenlcs成像方法 | |
CN109270528B (zh) | 基于全解析距离模型的一站固定式双站sar成像方法 | |
CN103091674A (zh) | 基于hrrp序列的空间目标高分辨成像方法 | |
CN109444882B (zh) | 基于变斜视椭圆波束同步模型的双站sar成像方法 | |
CN104991251B (zh) | 基于匀加速建模的超高分辨率星载sar成像方法 | |
CN104597447A (zh) | 一种子孔径SAR大斜视改进Omega-K成像方法 | |
CN103809180B (zh) | 用于InSAR地形测量的方位向预滤波处理方法 | |
CN106054187A (zh) | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 | |
CN111208515B (zh) | 一种基于二维非线性映射的sar运动补偿方法 | |
CN103576148B (zh) | 模拟星载sar距离模糊噪声图像的方法 | |
CN110361733B (zh) | 一种基于时频联合重采样的中轨sar大斜视成像方法 | |
CN103837874B (zh) | 用于地球同步轨道sar成像的二维非线性变调频方法 | |
CN114325704B (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN103823210B (zh) | 一种非合作式星地双基地sar时频同步方法 | |
CN105022060A (zh) | 针对快速空天目标的步进isar成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |