技术领域
本发明涉及合成孔径雷达成像方法,尤其涉及一种基于抽样序列长度约束条件下的SAR成像方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)广泛应用于军事和民用领域。目前距离多普勒(Range Doppler,RD)方法虽然在许多模式的SAR成像处理中,尤其是正侧视SAR成像处理中仍然广为使用,但是其较低精度的SAR图像质量越来越不能满足当前实际应用的需要。
分数阶Fourier变换(Fractional Fourier Transform,FrFT)是一种广义的Fourier变换,它是一种统一的时频变换,随着变换阶数从0连续增长到1而展示出信号从时域逐步变化到频域的所有特征。近年来,将分数阶Fourier变换应用于SAR成像处理中,已有文献通过定义旋转中心距离,利用分数阶Fourier变换得到基于参数的聚焦方法,该方法适用于条带SAR、聚束SAR和扫描SAR成像,三种模式的SAR数据都可以通过选择合适的参数以及分数阶Fourier变换旋转角来进行处理。已有文献利用分数阶Fourier变换在旋转的时频平面内获得单视复图像(single-look complex,SLC)的特征描述符,并与传统多尺度方法如Gabor滤波、二阶统计量和谱分析法作对比分析,测试结果表明该方法更适合用来对地面目标进行SAR图像分类。已有文献将分数阶Fourier变换应用于传统的RD方法,虽然可获得SAR成像性能的提高,但同时计算复杂度也相应增加。已有文献通过局部最优处理来测量SAR回波信号的调频率并以此计算FrFT的最优阶数,所研究方法在提高弹载SAR成像性能方面是有效的但不具有适用的普遍性。
为解决弹载SAR大距离徙动及空间时变性,已有文献提出一种新的基于分数阶Fourier变换的CS方法,该方法通过近似选择距离向和方位向分数阶Fourier变换的旋转角,获得比传统CS方法更优的SAR图像。已有文献将分数阶Fourier变换应用于星载SAR的数据处理,利用分数阶Fourier变换在旋转时频平面上表征SAR信号,达到最优处理和分析动目标在方位向上的残余Chirp信号。为获得地面动目标清晰的SAR图像,已有文献提出分数阶Fourier变换与自适应迭代模糊数方法相结合的多普勒参数估计方法,已有文献提出联合利用Wigner-Ville分布和分数阶Fourier变换实时估计多普勒参数的方法,并且观测信号的WVD处理决定分数阶Fourier变换的旋转角。已有文献针对地面合成孔径雷达(GB-SAR)或其它由物理或合成线性孔径构成的雷达系统,研究了基于分数傅里叶变换(FrFT)进行方位压缩的成像方法。同年,已有文献利用距离-频率变量的等间隔采样,研究了基于距离-频率反变换分数傅里叶变换(RFRT-FrFT)的合成孔径雷达地面运动目标成像新方法。在满足SAR成像基本要求的抽样序列长度条件下,已有文献提出适用于方位向的SAR局部最优成像方法。
尽管众多学者在分数阶Fourier变换聚焦特性方面做了大量研究工作,针对分数阶Fourier变换的SAR成像方法也进行了较为深入的研究,但是为解决传统距离多普勒方法成像精度低的问题,将分数阶Fourier变换应用于该方法以获得全局最优效果的SAR成像处理方法一直是一个亟待解决的技术问题。
发明内容
发明目的:本发明提出一种基于抽样序列长度约束条件下的SAR成像方法,即根据距离向和方位向最佳抽样长度的约束条件构建高分辨率SAR成像方法,旨在解决传统距离多普勒方法成像精度低的问题。
技术方案:本发明基于抽样序列长度约束条件下的SAR成像方法,包括以下步骤:
(1)SAR回波信号应用分数阶Fourier变换的阶数计算;
(1.1)信号模型构建
假定SAR在任意位置发射的LFM(Linear Frequency Modulation)信号如式(1)所示:
其中τ为快时间变量,T
(1.2)分数阶Fourier变换定义式
连续信号f(x)的分数阶Fourier变换定义式为
其中K
(1.3)方位向最优阶数计算
SAR方位向回波信号近似为LFM信号,如式(3)所示:
s
其中f
其中arctan(·)为反正切函数,N
(1.4)距离向最优阶数计算
考虑回波信号的延迟,根据式(1)和式(2)进行变换分析计算SAR距离向回波信号,运用分数阶Fourier变换时的最优阶数为
其中N
(2)设定抽样序列长度约束条件
(2.1)初始序列长度设定
设PRF表示脉冲重复频率(Pulse Repetition Frequency),则距离向初始序列长度
其中INT[·]表示取整函数,γ、δ表示变常数,在机载SAR成像时取γ=1.2、δ=1.2。对于星载SAR,由于卫星平台运行速度很快,SAR成像时合成孔径时间T
(2.2)最佳抽样序列长度确定
分别以初始序列
其中argmin[·]表示使目标函数取最小值时的变量值;
(3)高分辨SAR成像方法的构建
本发明主要包括距离向及方位向信号处理两个部分:
(3.1)距离向最佳抽样序列长度的确定
根据距离向成像参数进行最优阶数分析,利用式(5)计算得到距离向最优阶数μ
(3.2)方位向最佳抽样序列长度的确定
根据方位向成像参数进行最优阶数分析,利用式(4)计算得到方位向最优阶数ν
工作原理:本发明根据SAR距离向和方位向上的相关成像参数,设计初始抽样序列长度的计算公式;利用获得的初始抽样序列长度以及SAR图像质量性能评估的相应准则,在合理区间变化范围内结合曲线极值确定距离向和方位向最佳抽样序列长度,旨在为全局最优成像提供最佳抽样数据长度支撑,以此为基础构建高分辨SAR成像方法。
有益效果:与现有技术相比,本发明具有以下优点:
(1)本发明显著提高了实测线性调频信号的聚焦性能。
(2)本发明在合成孔径雷达成像中克服了传统距离多普勒方法成像质量低的问题,显著提高了SAR图像距离向和方位向分辨率,地物目标可区分度高。
(3)本发明得到的SAR图像纹理清晰且相干斑噪声颗粒小,显著增强了目标图像的判读和解译性能。
附图说明
图1为本发明高分辨率SAR成像方法构建流程图;
图2为本发明距离向抽样序列长度对SAR成像性能影响图;
其中图(2a)为距离向抽样序列长度对距离向分辨率影响图;图(2b)为距离向抽样序列长度对距离向PSLR和ISLR影响图;
图3为本发明方位向抽样序列长度对SAR成像性能影响图;
其中图(3a)为方位向抽样序列长度对方位向分辨率影响图;图(3b)为方位向抽样序列长度对方位向PSLR和ISLR影响图;
图4为本发明SAR实测数据成像图(局部放大图);
其中图(4a)为传统RD方法图;图(4b)为本发明方法图;
具体实施方式
如图1所示,本发明基于抽样序列长度约束条件下的SAR成像方法包括以下步骤:
步骤(1),SAR回波信号应用分数阶Fourier变换的阶数计算;具体包括以下内容:
(1.1)信号模型构建
假定SAR在任意位置发射的LFM(Linear Frequency Modulation)信号如式(1)所示:
其中τ为快时间变量,T
(1.2)分数阶Fourier变换定义式
连续信号f(x)的分数阶Fourier变换定义式为
其中K
(1.3)方位向最优阶数计算
SAR方位向回波信号近似为LFM信号,如式(3)所示:
s
其中f
其中arctan(·)为反正切函数,N
(1.4)距离向最优阶数计算
考虑回波信号的延迟,根据式(1)和式(2)进行类似变换分析计算SAR距离向回波信号运用分数阶Fourier变换时的最优阶数为
其中N
步骤(2),抽样序列长度约束条件;具体包括以下内容:
(2.1)初始序列长度设定
设PRF表示脉冲重复频率(Pulse Repetition Frequency),则距离向初始序列长度
其中INT[·]表示取整函数,γ、δ表示变常数,在机载SAR成像时可取γ=1.2、δ=1.2。对于星载SAR,由于卫星平台运行速度很快,SAR成像时合成孔径时间T
(2.2)最佳抽样序列长度确定
分别以初始序列
其中argmin[·]表示使目标函数取最小值时的变量值;
步骤(3),高分辨SAR成像方法的构建,具体包括以下内容:
如图1所示,本发明包括距离向及方位向信号处理两个部分:
(3.1)距离向最佳抽样序列长度的确定
根据距离向成像参数进行最优阶数分析,利用式(5)计算得到距离向最优阶数μ
(3.2)方位向最佳抽样序列长度的确定
根据方位向成像参数进行最优阶数分析,利用式(4)计算得到方位向最优阶数ν
自此,完成了基于抽样序列长度约束条件下的SAR成像方法。
实验结果及分析
以正侧视点目标成像为例,设定机载SAR的仿真参数:LFM信号带宽为120MHz,方位向天线尺寸为3m,光速为2.998×10
根据设置参数,由式(7)计算得到方位向抽样序列长度初始值为
实际SAR成像处理系统中距离向和方位向最佳抽样序列长度将受SAR目标场景能量反射的波动程度以及抽样序列长度变化区间中步进长度的影响。理论上步进长度越小,选取的最佳抽样序列长度精度越高。对于加拿大RADARSAT-1相同场景(能量反射的波动程度相当)的原始回波数据,其相关参数为:距离向抽样率F
由图4中的图(4a)和图(4b)可知,经传统RD方法得到的SAR图像目标模糊不清,较低的分辨率严重影响目标判读和解译;而经本发明得到的SAR图像中,边界湾机场设施轮廓明显,跑道清晰可见;布莱恩高速公路边缘分界线易辨认;城市建筑反射信号强、目标聚焦性好,“丁”字形马路可见;其他不同地物目标可区分度高,纹理清晰,相干斑噪声颗粒小,因而显著增强了SAR图像目标的解译和识别性能。故对于实测数据成像,本发明在距离向和方位向通过应用最佳抽样序列长度,在局部最优成像性能基础上进一步提升了SAR成像质量,并获得全局最优成像。与传统的RD方法相比,其成像性能得到显著提高,但实际运算速度慢一些。
机译: 磁共振成像设备,SAR计算设备,磁共振成像设备的操作方法和SAR计算方法
机译: 磁共振成像系统,计算SAR的装置,磁共振成像系统的工作方法和计算SAR的方法
机译: 将二进制序列以最大转换密度约束将映射以2/3(1,K)游程长度限制的编码序列进行映射的方法和方法