基于最优分数域Gabor谱特征的SAR目标检测方法
【专利摘要】本发明公开了一种基于最优分数域Gabor谱特征的SAR目标检测方法，属于图像处理【技术领域】，其包括如下步骤：读入SAR图像信号；将所述SAR图像信号按行、列方向展开为和的信号；对两个方向的信号设计最优窗函数，并做GT；对两个方向得到的空间-频率谱做FrFT；对空间-频率谱进行能量衰减梯度特征提取；将两个方向的特征表示对应空间位置做乘积，从而得到原始SAR图像的特征空间，完成检测。本发明将FrFT与GT结合，因为比常规算法多出了分数域参数，使得算法抗干扰性更强，且能够检测出微弱目标区域，比常规算检测算法精度更高，也适用于各种场景，具有较好的通用性。
【专利说明】基于最优分数域Gabor谱特征的SAR目标检测方法
【技术领域】
[0001]本发明涉及图像处理【技术领域】，具体涉及一种基于最优分数域Gabor谱特征的SAR目标检测方法，可用于SAR目标的检测，其基于分数阶傅里叶变换(FractionalFourier Transform, FrFT)与时频分析方法(Time Frequency Analysis, TFA),提取图像的分数域Gabor谱特征，从而高精度检测出目标。
【背景技术】
[0002]SAR具有全天时、全天候、穿透能力强、分辨率高等特点，已被广泛运用于军事侦察、地形测绘、地质勘探、海洋应用、农林监测、水资源利用和科学研究等领域。SAR发展迅速，SAR图像系统的信息获取技术明显超前于SAR图像成像后的信息处理技术，因此各个国家都对 SAR ATR (Automation Target Recognition, ATR)展开了大量研究工作。SAR 目标检测是SAR ATR的关键环节。
[0003]FrFT实际上就是Fourier变换的广义形式。FrFT的概念其实最早在1929年，就由N.Wiener提出，但是并没有受到太多的瞩目。1980年Namias从特征值和特征函数的角度，以纯数学的方式重新提出它的概念，并揭示了 FRFT的几个特性以后，FrFT变换逐渐受到各国学者的重视。1993年Almeida指出分数阶傅里叶变换可以理解为时频平面的旋转。1996年Ozaktas等提出了一种计算量与FFT相当的离散算法后，开始出现大量相关研究文章。目前FrFT已广泛应用于雷达、声纳、通信、信息安全等众多领域。
[0004]时频分析的思想始于上世纪四十年代年D.Gabor提出的Gabor变换。常用的时频分析方法包括短时傅里叶变换(Short Time Fourier Transform, STFT)、Gabor变换(GaborTransform, GT)、小波变换、Winger-Vile 分布(Wigner-Vile Distributioan, WVD)、S 变换等。Gabor变换是通过时域窗函数将时间信号进行截断，对截断后的局部信号做Fourier变换，并表示为时间和频率的联合分布，来揭示信号的时变谱特征，满足瞬时信号分析的要求，更适用于非平稳信号的研究。对图像信号进行时频分析实际上就是进行图像信号的空间-频率分析。
[0005]近年来兴起的分数域时频分析方法正是在FrFT的基础上发展起来的。与传统时频分析方法相比，在某些条件下，分数域时频分析能够取得更好的效果，这是因为分数域时频分析方法比传统时频分析多了一个分数域参量。目前，分数域时频分析方法已经被广泛用在科学研究与工程应用的很多相关领域，如相位恢复、信号重构、信号检测、参数估计、分数域滤波器设计、神经网络、图像数字水印、图像加密、油气勘测等。
[0006]由于SAR图像的成像机理，所以其不能很好的描述轮廓和细节；其次，是SAR目标对方位角敏感，因此，不同目标及不同照射方位下所形成的SAR图像，其背景散射强度与目标散射强度也不相同；而且，SAR图像存在特有的斑点噪声，且在成像过程中辐射特征以及几何特征发生畸变，常规的图像检测方法不适用于SAR目标检测当中，因此需要一个高精度，抗噪性强的检测方法。
【发明内容】
[0007]针对上述现有技术，本发明的目的在于如何提供一种基于最优分数域Gabor谱特征的SAR目标检测方法，其旨在抑制SAR图像相干斑噪声，提升传统Gabor变换时频的聚集性，改善时频分布，并能更精确的完成SAR目标检测。
[0008]为达到上述目的，本发明采用如下技术方案:
[0009]一种基于最优分数域Gabor谱特征的SAR目标检测方法，其特片在于，包括如下步骤:
[0010]步骤一:读入mXn SAR图像信号；
[0011]步骤二:将所述SAR图像信号按行、列方向展开为I X丽和丽X I的信号；
[0012]步骤三:对两个方向的信号设计最优窗函数，并分别对它们做GT ；
[0013]步骤四:对两个方向得到的空间-频率谱做FrFT，并搜索最优变换阶次，改善空间_频率分布；
[0014]步骤五:对经过步骤三和步骤四优化后的空间-频率谱进行能量衰减梯度特征提取，分别得到两个方向的特征表示；
[0015]步骤六:将两个方向的特征表示对应空间位置做乘积，从而得到原始SAR图像的特征空间，完成检测。
[0016]所述步骤三具体又分解为如下步骤:
[0017]步骤3.1:对行方向信号进行GT后，得到丽X丽的空间-频率谱，利用FrFT的旋转性，并结合广义时频带宽积准则(Generalized Time-Frequency band product, GTBP)设计出最优窗；
[0018]步骤3.2:对于列方向也进行与步骤3.1相同的处理。
[0019]所述步骤四具体包括如下步骤:
[0020]步骤4.1:以经过步骤3处理过的行方向的空间-频率谱为基础，固定空间位置，抽取对应位置以频率表示的MN道信号，分别对每一道信号进行不同阶次的FrFT，并通过能量峰值法对每一道信号进行最优阶次的FrFT ；
[0021]步骤4.2:对于列方向也进行与步骤3.1相同的处理。
[0022]所述步骤五具体包括如下步骤:
[0023]步骤5.1:以行方向最优空间-频率谱为基础，固定空间位置，抽取对应位置以频率表示的丽道信号，找到每一道信号峰值Max的位置，记为fmax，在fmax之后找到总能量的65%对应的位置，记为& ;在fmax之后找到总能量的85%对应的位置，记为;计算&，f1； fmax的拟合斜率得到IXMN的行方向空间-频率谱特征表示。
[0024]步骤5.2:对于列方向也进行与步骤3.1相同的处理。
[0025]步骤5.3:对行方向IXMN按对应位置重置为MXN，得到原始SAR图像的行方向的特征表示；对行方向丽X I按对应位置重置为MXN，得到原始SAR图像的列方向的特征表示；
[0026]与现有技术相比，本发明具有以下有益效果:
[0027]—、本发明采用分数域时频分析方法，并成功将其应用于SAR图像处理，对非平稳信号和信号瞬态分析提供了一种新的思路；本发明抗噪性强、检测精度高，时频特征剖面加大了目标区域和背景区域的对比度，同时抑制了相干斑噪声，克服了常规检测方法对噪声敏感的缺点。
[0028]二、本发明将FrFT引入窗函数的设计，利用分数阶旋转性，并搜索最优旋转因子，找到满足GTBP的最优阶次，完成窗函数设计，显著提升空间-频率聚集性。
[0029]三、本发明将FrFT与GT结合，因为比常规算法多出了分数域参数，使得算法抗干扰性更强，且能够检测出微弱目标区域，比常规算检测算法精度更高，也适用于各种场景，具有较好的通用性。
【专利附图】
【附图说明】
[0030]图1为本发明方法流程示意图；
[0031]图2为FrFT的对时频面的旋转性示意图；
[0032]图3为最优旋转因子搜索示意图；
[0033]图4为能量衰减梯度示意图。
【具体实施方式】
[0034]下面将结合附图及【具体实施方式】对本发明作进一步的描述。
[0035]一种基于最优分数域Gabor谱特征的SAR目标检测方法，该方法包括以下步骤:
[0036]步骤一:读入MXN维SAR图像信号；
[0037]步骤二:将读入的图像信号按行、列方向展开为I X丽和丽X I的信号；
[0038]步骤三:对两个方向的信号设计最优窗函数，并分别对它们做GT，使得GT的空间-频率聚集性得到显著提升，从而抑制杂波。
[0039]步骤四:对两个方向得到的空间-频率谱做FrFT，并搜索最优变换阶次，改善空间-频率分布，从而对微弱目标敏感，提高检测率；
[0040]步骤五:对经过步骤三和步骤4优化后的空间-频率谱进行能量衰减梯度特征提取，分别得到两个方向的特征表示；
[0041]步骤六:将两个方向的特征表示对应空间位置做乘积，从而得到原始SAR图像的特征空间，完成检测。
[0042]所述步骤三具体包括如下步骤:
[0043]步骤3.1:对行方向信号进行GT后，得到丽X丽的空间-频率谱。利用FrFT的旋转性，并结合广义时频带宽积准则(Generalized Time-Frequency band product, GTBP),使得空间宽和频率宽都达到最小，从而设计出最优窗，使得GT拥有很好的空频聚集性。
[0044]如图2所示FrFT对时频面的旋转,相当于把信号的时频轴在时频面上逆时针旋转α = P Ji /2，P为变换阶次。
[0045]步骤3.1.1:改变旋转阶次使得其在分数域上持续时间最短，从而拥有最大的模值，该阶次既为最优旋转因子。
[0046]如图3所示，图3(a)所示为信号的原始时频图，图3(b)所示，当时频图的长轴恰好和分数域频率轴重合，此时时频图的短轴将与时间轴重合时，则信号持续时间最短，时频图的模值最大，据此可以通过遍历不同阶次FrFT，来确定最优旋转因子。
[0047]步骤3.1.2:求得最优旋转因子后，计算出最优窗函数；
[0048]步骤3.2:对于列方向也进行与步骤3.1相同的处理。[0049]如图4所示，对理论chirp信号分别用最优窗进行GT和普通窗进行GT，时频聚集性上得到了显著的提升，能够很好的抑制因GT本身而产生的杂波，为高精度检测打下基础。
[0050]所述步骤四具体包括如下步骤:
[0051]步骤4.1:以经过步骤3处理过的行方向的空间-频率谱为基础。固定空间位置，抽取对应位置以频率表示的MN道信号，分别对每一道信号进行不同阶次的FrFT，并通过能量峰值法对对每个空间位置的频率信号都进行最优FrFT，最优阶由下式确定:
[0052]Fopt = arg max( Energyλ>)
[0053]其中，￡mrgf *表示不同阶次下最优窗Gabor变换得到的空间-频率谱所抽取的每个空间位置的频率信号的能量值，从而改善时频分布，使得本方法对微弱目标更敏感，提高了检测精度。
[0054]步骤4.2:对于列方向也进行与步骤3.1相同的处理。
[0055]对两个方向得到的空间-频率谱做FrFT，并搜索最优变换阶次，改善空间-频率分布，从而对微弱目标敏感，提闻检测率。
[0056]所述步骤五具体包括如下步骤:
[0057]通过步骤三与步骤四的处理，我们得到了经过最优窗设计和最优FrFT后的最优空间-频率谱。
[0058]步骤5.1:以行方向最优空间-频率谱为基础，固定空间位置，抽取对应位置以频率表示的丽道信号，找到每一道信号峰值Ma`x的位置，记为fmax，在fmax之后找到总能量的65%对应的位置，记为& ;在fmax之后找到总能量的85%对应的位置，记为;计算&，f1； fmax的拟合斜率，实现了空间-频率谱的降维，得到I XMN的行方向空间-频率谱特征表示；
[0059]步骤5.2:对于列方向也进行与步骤3.1相同的处理；
[0060]步骤5.3:对行方向IXMN按对应位置重置为MXN，得到原始SAR图像的行方向的特征表示；对行方向丽X I按对应位置重置为MXN，得到原始SAR图像的列方向的特征表
/Jn ο
[0061]步骤六:将两个方向的特征表示对应空间位置做乘积，从而得到原始SAR图像的特征空间，完成检测。
[0062]传统FT是一种全局性变换，能够得到信号的整体频谱，因此，FT不能用来有效处理非平稳信号。为了克服FT在信号处理中的局限性，更好的分析信号的局部特性，提出了FrFT,时间信号X (t)的P阶FrFT变换定义为:
[0063]Xp(U) = ? χ(τ)ΚΡ(υ,τ)Λ
m—iXs
[0064]其中，p为FrFT的阶次，其周期为4，α = ρ π/2为分数域坐标轴u相对时间轴t逆时针旋转的角度，当P = I时实际上就是FT,其中Kp(t，u)表示FrFT的核函数,可以表示为:
【权利要求】
1.一种基于最优分数域Gabor谱特征的SAR目标检测方法，包括如下步骤: 步骤一:读入mXn SAR图像信号； 步骤二:将所述SAR图像信号按行、列方向展开为I X丽和丽X I的信号； 步骤三:对两个方向的信号设计最优窗函数，并分别对它们做GT ； 步骤四:对两个方向得到的空间-频率谱做FrFT，并搜索最优变换阶次，改善空间-频率分布； 步骤五:对经过步骤三和步骤四优化后的空间-频率谱进行能量衰减梯度特征提取，分别得到两个方向的特征表示； 步骤六:将两个方向的特征表示对应空间位置做乘积，从而得到原始SAR图像的特征空间，完成检测。
2.根据权利要求1所述的基于最优分数域Gabor谱特征的SAR目标检测方法，其特片在于，所述步骤三具体又分解为如下步骤: 步骤3.1:对行方向信号进行GT后，得到MNXMN的空间-频率谱，利用FrFT的旋转性，并结合广义时频带宽积准则设计出最优窗； 步骤3.2:对于列方向也进行与步骤3.1相同的处理。
3.根据权利要求1所述的基于最优分数域Gabor谱特征的SAR目标检测方法，其特片在于，所述步骤四具体包括如下步骤: 步骤4.1:以经过步骤3处理过的行方向的空间-频率谱为基础，固定空间位置，抽取对应位置以频率表示的MN道信号，分别对每一道信号进行不同阶次的FrFT，并通过能量峰值法对每一道信号进行最优阶次的FrFT ； 步骤4.2:对于列方向也进行与步骤3.1相同的处理。
4.根据权利要求1所述的基于最优分数域Gabor谱特征的SAR目标检测方法，其特片在于，所述步骤五具体包括如下步骤: 步骤5.1:以行方向最优空间-频率谱为基础，固定空间位置，抽取对应位置以频率表示的MN道信号，找到每一道信号峰值Max的位置，记为fmax，在fmax之后找到总能量的65%对应的位置，记为& ;在fmax之后找到总能量的85%对应的位置，记为;计算&，f1； fmax的拟合斜率得到IXMN的行方向空间-频率谱特征表示。 步骤5.2:对于列方向也进行与步骤3.1相同的处理。 步骤5.3:对行方向IXMN按对应位置重置为MXN，得到原始SAR图像的行方向的特征表示；对行方向MNX I按对应位置重置为MXN，得到原始SAR图像的列方向的特征表示。
【文档编号】G06T5/00GK103456015SQ201310403403
【公开日】2013年12月18日 申请日期:2013年9月6日 优先权日:2013年9月6日
【发明者】彭真明, 杨俊涛, 黄振星, 魏瑞鹏, 何艳敏, 张萍 申请人:电子科技大学