CN111624602A - 一种旋转FrFT域构建SAR方位向分辨率成像方法 - Google Patents
一种旋转FrFT域构建SAR方位向分辨率成像方法 Download PDFInfo
- Publication number
- CN111624602A CN111624602A CN202010668076.7A CN202010668076A CN111624602A CN 111624602 A CN111624602 A CN 111624602A CN 202010668076 A CN202010668076 A CN 202010668076A CN 111624602 A CN111624602 A CN 111624602A
- Authority
- CN
- China
- Prior art keywords
- azimuth
- sar
- fourier transform
- signal
- fractional
- 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.)
- Pending
Links
Images
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
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种旋转分数阶Fourier变换(Fractional Fourier Transform,FrFT)域构建SAR(Synthetic Aperture Radar)方位向高分辨率成像方法,该方法通过对SAR方位向信号运用分数阶Fourier变换,得出方位向分数阶阶数的计算表达式,该阶数具有唯一性且取决于方位向调频率、方位向抽样序列长度和方位向抽样频率;进而以获得的分数阶阶数为基础完成旋转分数阶Fourier变换域和方位向高分辨率SAR成像方法的构建。本发明在合成孔径雷达成像中克服了传统方法方位向低精度成像的问题,提高了SAR方位向成像分辨率,满足了当前实际应用需求。
Description
技术领域
本发明涉及合成孔径雷达成像方法,尤其涉及一种旋转FrFT域构建SAR方位向分辨率成像方法。
背景技术
作为一种主动式微波传感器,合成孔径雷达具有不受光照和气候条件限制而实现全天时、全天候对地观测的特点,可穿透地表或植被获取其掩盖的信息,在民用和军事领域具有广阔的应用前景。
SAR通过把小孔径天线雷达装载在运动的载体上,利用雷达与地面测绘带内各种目标的相对运动进而利用相干处理实现距离向和方位向成像。常见的SAR成像方法有距离多普勒(Range Doppler,RD)方法、波数域(ωK)方法和CS(Chirp Scaling)方法,其中距离多普勒方法是SAR成像处理中最常见、最经典的方法。目前距离多普勒方法虽然在许多模式的SAR尤其是正侧视SAR成像处理中仍然广为使用,但是其较低质量的SAR图像越来越不能满足实际应用需求。
众所周知,Fourier变换在传统距离多普勒方法成像处理中起着至关重要的作用,而同时Fourier变换也是分数阶Fourier变换(Fractional Fourier Transform,FrFT)阶数为1时的特例。分数阶Fourier变换是一种广义的Fourier变换,它是一种统一的时频变换,随着变换阶数从0连续增长到1而展示出信号从时域逐步变化到频域的所有特征。将分数阶Fourier变换应用于SAR成像处理是近年以来的研究热点。已有文献针对线性调频(LinearFrequency Modulation,LFM)信号的参数估计问题提出在分数阶Fourier域进行二维谱峰搜索来确定分数阶Fourier变换的最佳变换阶数,该方法稳定性较好但搜索分数阶谱峰计算量大,同时有限数据样本会导致参数估计精度的降低。已有文献利用几何变换关系得到适用于线性调频信号的分数阶Fourier变换分数阶变换阶数,但相应的分数阶Fourier变换无法代替传统距离多普勒方法中的Fourier变换进而实现信号重构。已有文献将分数阶Fourier变换应用于传统的RD方法,虽然获得了SAR成像性能的提高,但同时计算复杂度也相应增加。已有文献通过局部最优处理来测量SAR回波信号的调频率并以此计算FrFT的分数阶阶数,所研究方法在提高弹载SAR成像性能方面是有效的但不具有适用的普遍性。为获得地面动目标清晰的SAR图像,已有文献提出分数阶Fourier变换与自适应迭代模糊数方法相结合的多普勒参数估计方法,已有文献提出联合利用Wigner-Ville分布和分数阶Fourier变换实时估计多普勒参数的方法,并且观测信号的WVD处理决定分数阶Fourier变换的旋转角。
发明内容
发明目的:本发明提出一种旋转分数阶Fourier变换FrFT域构建SAR方位向分辨率成像方法,该方法提高了SAR方位向成像分辨率,解决了传统方法方位向成像质量低的问题。
技术方案:本发明旋转分数阶Fourier变换FrFT域构建SAR方位向分辨率成像方法,该方法包括以下步骤:
(1)SAR方位向信号变换分析;具体包括以下内容:
(1.1)构建方位向信号模型;
点目标回波在方位向为线性调频信号的形式,如公式(1)所示:
sa(t)=Wa(t)exp(j2πfdct+jπκat2) (1)
其中fdc为多普勒中心频率,κa为方位向调频率,Wa(t)=rect(t/Ta),Ta为合成孔径时间,t为慢时间变量;
其中连续信号f(x)的分数阶Fourier变换定义式为
其中Kβ(u,x)为分数阶Fourier变换的核函数,如公式(3)所示
(1.2)利用分数阶傅里叶变换处理方位向信号;
(1.3)获取信号的驻定相位点;
(1.4)获取泰勒展开级数并代入分析;
(1.5)变量变换分析;
(1.6)菲涅耳积分计算;
(1.7)变量替换;
(1.8)方位向分数阶阶数计算;
(1.9)能量聚焦轴计算;
(2)旋转分数阶Fourier变换域的构建;
(3)方位向高分辨率成像方法的构建;包括以下内容:
(3.1)利用已知的SAR成像参数,计算出方位向上的分数阶阶数νopt;
(3.3)对完成距离向信号处理的SAR回波方位向信号、方位压缩参考函数均作阶数为ν′opt的分数阶Fourier变换,完成方位向脉冲压缩;
(3.4)对方位向信号作阶数为1的分数阶Fourier逆变换,完成方位向信号重构。
步骤(2)中,对于SAR接收的原始回波信号作分数阶Fourier变换,采用公式(12)
其中,κa为方位向调频率、Na为抽样序列长度,Fa为抽样频率,计算得到方位向分数阶阶数νopt。
步骤(12)中,将公式(1)代入公式(2)得
步骤(13)中,通过公式(5)
步骤(1.7)中,假设SAR回波方位向信号抽样序列长度为Na,抽样频率为Fa,对方位向信号量纲归一化后调频率为κ′a,旋转角度为β′,代入公式(10)
得
步骤(1.8)中,根据公式(11)得到SAR回波方位向信号进行分数阶Fourier变换时的分数阶阶数
其中arctan(·)为反正切函数。
工作原理:本发明通过对SAR方位向信号运用分数阶Fourier变换,得出方位向分数阶阶数的计算表达式,该阶数具有唯一性且取决于方位向调频率、方位向抽样序列长度和方位向抽样频率;进而以获得的分数阶阶数为基础完成旋转分数阶Fourier变换域和方位向高分辨率SAR成像方法的构建。本发明在合成孔径雷达成像中克服了传统方法方位向低精度成像的问题,提高了SAR方位向成像分辨率,满足当前实际应用需求。
有益效果:与现有技术相比,本发明具有以下优点:
(1)本发明构建的旋转分数阶傅里叶变换域,显著提高了线性调频信号的方位向成像聚焦性能。
(2)本发明显著提高了机载SAR仿真数据和星载SAR实测数据的方位向成像分辨率。
(3)本发明应用于SAR数据成像时,明显降低峰值点目标方位向主瓣宽度,明显降低旁瓣电平,尤其是第一旁瓣能量显著降低。
附图说明
图1为本发明方位向信号时频域旋转图;
图2为本发明方位向高分辨率成像方法的构建流程图;
图3为本发明插值后点目标的轮廓图;
其中,图(3a)为RD方法图;图(3b)为本发明方法图;
图4为插值后峰值点目标成像方位向剖面图;
图5为SAR实测数据成像;
其中,图(5a)为RD方法图;图(5b)为本发明方法图;
图6为插值后峰值点目标成像轮廓图;
其中,图(6a)为RD方法图;图(6b)为本发明方法图;
图7为插值后峰值点目标成像方位向剖面图。
具体实施方式
如图1所示,本发明旋转分数阶傅里叶变换FrFT域构建SAR方位向分辨率成像方法包括以下步骤:
步骤(1),SAR方位向信号变换分析;具体包括以下内容:
①构建方位向信号模型
在SAR成像系统中,理想的点目标回波在方位向近似为线性调频信号的形式,如公式(1)所示:
sa(t)=Wa(t)exp(j2πfdct+jπκat2) (1)
其中fdc为多普勒中心频率,κa为方位向调频率,Wa(t)=rect(t/Ta),Ta为合成孔径时间,t为慢时间变量。连续信号f(x)的分数阶Fourier变换定义式为:
其中Kβ(u,x)为分数阶Fourier变换的核函数,如公式(3)所示,β为旋转角度且ν为分数阶Fourier变换的阶数。公式(3)中β≠2nπ,若β=2nπ,Kβ(u,x)=δ(u-x);若β=(2n+1)π,Kβ(u,x)=δ(u+x)。
②利用分数阶傅里叶变换处理方位向信号;
将公式(1)代入公式(2)得
③获取信号的驻定相位点;
为简化计算,仅考虑SAR条带模式下正侧视成像情况,此时fdc=0。令U(t)=Wa(t),V(t)=exp{jπ[(t2+u2)cotβ-2utcscβ+κat2]},根据驻定相位定理(Principal ofStationary Phase),通过公式(5)
④获取泰勒(Taylor)展开级数并代入分析;
利用公式(5)计算得到的驻定相位点tk,将公式(4)改写为
将公式(7)代入公式(6)得
⑤变量变换分析;
⑥菲涅耳积分计算;
⑦变量替换
假设SAR回波方位向信号抽样序列长度为Na,抽样频率为Fa,对方位向信号量纲归一化后调频率为κ′a,旋转角度为β′,代入公式(10)得
⑧方位向分数阶阶数计算;
根据公式(11),若使分数阶能量谱Fβ′[sa(u)]在β′旋转角度下高度聚焦在某个分数阶Fourier变换域轴上,必然满足κ′a+cotβ′=0,此时Fβ′[sa(u)]为冲激函数。由于β′=β,因此得到SAR回波方位向信号进行分数阶Fourier变换时的分数阶阶数
其中arctan(·)为反正切函数。对于实际给定的SAR回波抽样信号,方位向调频率κa、抽样序列长度Na和抽样频率Fa均是已知的,因此根据公式(12)直接计算出相应的分数阶阶数νopt。
⑨能量聚焦轴计算;
在获得分数阶阶数νopt情况下,根据冲激函数的Fourier变换定义,令公式(4)中慢时间变量t的系数为零,得到Fβ′[sa(u)]的能量谱高度聚焦在分数阶Fourier变换域轴上,如公式(13)所示:
步骤(2),旋转分数阶Fourier变换域的构建;
对于SAR接收的原始回波信号作分数阶Fourier变换,采用公式(12)计算得到方位向分数阶阶数νopt,相对应的时频域旋转角为β1=-νopt×π/2,如图1所示,其中负号表示逆时针旋转,tν和fν分别表示方位向时间域和方位向频率域。图1中,表示对方位向时域信号作旋转角为β1时的FrFT相对应的分数阶Fourier变换域,表示对分数阶Fourier变换域中信号作旋转角为β2(β2=π/2,对应阶数为1)时的FrFT相对应的分数阶Fourier变换域,t′ν表示tν的平行轴,显然tν⊥fν,方位向回波信号在时频域(tν,fν)中分解形式与其在分数阶Fourier变换域中的分解形式是等价的,且后者更利于方位向回波信号能量的聚焦。若域中分解信号等同于fν域中的频率信号,则域中分解信号等同于tν域中的时间信号。
步骤(3),方位向高分辨率成像方法的构建:
(31)利用已知的SAR成像参数,采用公式(12)计算出方位向上的分数阶阶数νopt。
(32)为提高方位向信号在分数阶Fourier变换域中的成像聚焦效果,将相应的时频角旋转π/2变换至获得方位向分数阶阶数为ν′opt=1-νopt。当νopt=0即ν′opt=1时,本发明退化为传统的距离多普勒方法。
(33)对完成距离向信号处理的SAR回波方位向信号、方位压缩参考函数均作阶数为ν′opt的分数阶Fourier变换,完成方位向脉冲压缩。
(34)对方位向信号作阶数为1的分数阶Fourier逆变换,完成方位向信号重构。
自此,完成了基于分数阶Fourier变换的高分辨率SAR成像方法的构建。
实验结果及分析
表1给出了机载SAR正侧视点目标成像仿真参数。仿真实验中所选脉冲重复频率(PRF)的方位向过抽样率系数为1.5,3个点目标幅度系数分别为0.8、1和1.2。对于采用传统的RD方法、本发明的3个点目标成像,由表2数据可知,在峰值旁瓣比(Peak Side LobeRatio,PSLR)和积分旁瓣比(Integrated Side lobe Ratio,ISLR)相当的情况下,本发明获得的方位向分辨率绝对值降低0.42,相应的分辨率提高比值为31.3%。本发明的方位向高分辨率成像在图3和图4中也得到了明显体现。图3所示的插值后点目标成像轮廓图表明,经本发明成像得到的点目标其方位向主瓣能量聚焦性能更优,相应的旁瓣能量明显减弱。图4所示的插值后峰值点目标成像剖面图表明,经本发明成像的点目标其方位向波形主瓣宽度明显变窄,旁(副)瓣电平也相应降低,因而其成像的清晰度得到显著提高。
图5给出了星载SAR实测数据成像结果,其原始数据来自于加拿大RADARSAT-1精细模式下温哥华场景(截取)。图5成像结果表明:与传统RD方法相比,经本发明得到的SAR图像分辨率高,所获得的图像中自然景观公园、马路和高速公路清晰可见,斑点噪声小,纹理和边界线明显,强点目标突出,强弱目标对比度高。图6和图7分别给出了与图5相对应插值后的峰值点目标成像轮廓图、峰值点目标成像方位向剖面图。表3给出了方位向成像性能对比。由图6和图7可知,由本发明成像的峰值点目标方位向主瓣宽度明显变窄,旁瓣电平明显降低,尤其是第一旁瓣能量显著降低。由表3中实验结果可知,在PSLR和ISLR相当的情况下,本发明成像方位向分辨率绝对值比传统RD方法降低12.1,相应的方位向分辨率提高比值为41.8%。
以上机载SAR仿真数据和星载SAR实测数据成像实验表明,由本发明得到的SAR图像均获得比传统RD方法更优异的成像性能,这是因为本发明在方位向依据最优阶数获得旋转分数阶Fourier变换域具有比传统Fourier域更好的目标成像聚焦效果。
表1机载SAR成像仿真参数
表2方位向成像性能对比
表3方位向成像性能对比
Claims (7)
1.一种旋转FrFT域构建SAR方位向分辨率成像方法,其特征在于:所述方法包括以下步骤:
(1)SAR方位向信号变换分析;具体包括以下内容:
(1.1)构建方位向信号模型;
点目标回波在方位向为线性调频信号的形式,如公式(1)所示:
sa(t)=Wa(t)exp(j2πfdct+jπκat2) (1)
其中fdc为多普勒中心频率,κa为方位向调频率,Wa(t)=rect(t/Ta),Ta为合成孔径时间,t为慢时间变量;
其中连续信号f(x)的分数阶Fourier变换定义式为
其中Kβ(u,x)为分数阶Fourier变换的核函数,如公式(3)所示
(1.2)利用分数阶傅里叶变换处理方位向信号;
(1.3)获取信号的驻定相位点;
(1.4)获取泰勒展开级数并代入分析;
(1.5)变量变换分析;
(1.6)菲涅耳积分计算;
(1.7)变量替换;
(1.8)方位向分数阶阶数计算;
(1.9)能量聚焦轴计算;
(2)旋转分数阶Fourier变换域的构建;
(3)方位向高分辨率成像方法的构建;包括以下内容:
(3.1)利用已知的SAR成像参数,计算出方位向上的分数阶阶数νopt;
(3.3)对完成距离向信号处理的SAR回波方位向信号、方位压缩参考函数均作阶数为ν′opt的分数阶Fourier变换,完成方位向脉冲压缩;
(3.4)对方位向信号作阶数为1的分数阶Fourier逆变换,完成方位向信号重构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010668076.7A CN111624602A (zh) | 2020-07-13 | 2020-07-13 | 一种旋转FrFT域构建SAR方位向分辨率成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010668076.7A CN111624602A (zh) | 2020-07-13 | 2020-07-13 | 一种旋转FrFT域构建SAR方位向分辨率成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111624602A true CN111624602A (zh) | 2020-09-04 |
Family
ID=72271464
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010668076.7A Pending CN111624602A (zh) | 2020-07-13 | 2020-07-13 | 一种旋转FrFT域构建SAR方位向分辨率成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111624602A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112559973A (zh) * | 2021-02-25 | 2021-03-26 | 中国人民解放军空军预警学院 | 基于STFrFT的自适应多分量线性调频信号参数估计方法 |
CN114866165A (zh) * | 2021-06-29 | 2022-08-05 | 哈尔滨工业大学 | 一种多频段室内信号分布场的快速测量获取方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103456015A (zh) * | 2013-09-06 | 2013-12-18 | 电子科技大学 | 基于最优分数域Gabor谱特征的SAR目标检测方法 |
EP2743727A2 (en) * | 2014-01-16 | 2014-06-18 | Institute of Electronics, Chinese Academy of Sciences | Method for implementing high-resolution wide-swath spaceborne SAR system |
CN108459321A (zh) * | 2018-02-07 | 2018-08-28 | 杭州电子科技大学 | 基于距离-方位圆模型的大斜视高分辨率sar成像方法 |
-
2020
- 2020-07-13 CN CN202010668076.7A patent/CN111624602A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103456015A (zh) * | 2013-09-06 | 2013-12-18 | 电子科技大学 | 基于最优分数域Gabor谱特征的SAR目标检测方法 |
EP2743727A2 (en) * | 2014-01-16 | 2014-06-18 | Institute of Electronics, Chinese Academy of Sciences | Method for implementing high-resolution wide-swath spaceborne SAR system |
CN108459321A (zh) * | 2018-02-07 | 2018-08-28 | 杭州电子科技大学 | 基于距离-方位圆模型的大斜视高分辨率sar成像方法 |
Non-Patent Citations (1)
Title |
---|
ZHENLI WANG等: "Application of optimal FrFT order for improving the azimuth resolution of range Doppler imaging algorithm", 《IET IMAGE PROCESSING》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112559973A (zh) * | 2021-02-25 | 2021-03-26 | 中国人民解放军空军预警学院 | 基于STFrFT的自适应多分量线性调频信号参数估计方法 |
CN112559973B (zh) * | 2021-02-25 | 2021-05-14 | 中国人民解放军空军预警学院 | 自适应多分量线性调频信号参数估计方法 |
CN114866165A (zh) * | 2021-06-29 | 2022-08-05 | 哈尔滨工业大学 | 一种多频段室内信号分布场的快速测量获取方法 |
CN114866165B (zh) * | 2021-06-29 | 2024-04-26 | 哈尔滨工业大学 | 一种多频段室内信号分布场的快速测量获取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Motion compensation/autofocus in airborne synthetic aperture radar: A review | |
CN111142105B (zh) | 复杂运动目标isar成像方法 | |
Li et al. | Scaled Radon-Wigner transform imaging and scaling of maneuvering target | |
CN109407070B (zh) | 一种高轨平台地面运动目标检测方法 | |
CN111624602A (zh) | 一种旋转FrFT域构建SAR方位向分辨率成像方法 | |
Zhou et al. | Analysis of wide-angle radar imaging | |
CN103278819B (zh) | 基于滑动接收窗的机载高分辨率斜视聚束sar成像方法 | |
CN113484859B (zh) | 一种基于融合技术的二维超分辨雷达成像方法 | |
CN111060909A (zh) | 一种机载雷达斜前视超分辨成像方法 | |
Xingyu et al. | Approach for ISAR imaging of near-field targets based on coordinate conversion and image interpolation | |
CN112305541A (zh) | 一种基于抽样序列长度约束条件下的sar成像方法 | |
Wang et al. | Velocity estimation of moving targets in SAR imaging | |
Tran et al. | The slow-time k-space of radar tomography and applications to high-resolution target imaging | |
Xu et al. | A variable PRF imaging method for high squint diving SAR | |
CN116719027A (zh) | 一种双向滑动聚束模式下的星机双基sar成像方法 | |
CN115755046A (zh) | 一种用于条带sar大方位幅宽成像的扩展pfa算法 | |
CN112731392B (zh) | 高效的太赫兹步进频率合成孔径雷达成像方法 | |
Gorham | Large scene SAR image formation | |
Chen et al. | An improved range-Doppler algorithm for SAR imaging at high squint angles | |
Yang et al. | A subaperture imaging algorithm to highly squinted TOPS SAR based on SPECAN and deramping | |
CN109633641B (zh) | 一种太赫兹频段的旋翼叶片逆合成孔径雷达成像算法 | |
CN107422320A (zh) | 一种消除降雨对x波段雷达观测海浪的影响的方法 | |
CN110045374B (zh) | 一种基于Chirp Scaling的多角度超长聚束SAR超分辨率成像算法 | |
Wang et al. | High resolution SAR image algorithm with sample length constraints for the range direction | |
CN111965641B (zh) | 一种基于分数阶Fourier变换的SAR成像方法 |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Wang Zhenli Inventor after: Wang Qun Inventor after: Ma Rupo Inventor after: Cai Mengsheng Inventor before: Wang Zhenli Inventor before: Wang Qun Inventor before: Ma Rupo |
|
CB03 | Change of inventor or designer information |