SAR图像序列特征提取方法技术

一种圆迹SAR图像序列特征提取方法，包括：联合应用似然比和恒虚警检测实现对圆迹SAR初始图像的方位向多散射中心的检测、消除和非相干累加，得到新的圆迹SAR图像序列，并获得所述新的圆迹SAR图像的似然比曲线Λ(i)和子孔径消除曲线Ap(i)；通过对所述似然比曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向似然比特征；以及通过对子孔径消除曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向子孔径消除特征。

【技术实现步骤摘要】

本专利技术涉及雷达图像处理领域，尤其涉及一种圆迹合成孔径雷达(Synthetic Aperture Radar，SAR)图像序列的特征提取方法。
技术介绍
目前，在圆迹SAR成像过程中，通常对各子孔径进行全相干累加全孔径图像。若目标对所有孔径都满足方位向同散射中心假设，全相干累加可以获得分辨率和信噪比最优的图像。然而，由于实际地物的散射特征仅在有限的方位角内保持稳定，方位向多散射中心的存在使得圆迹SAR出现新的特征。传统圆迹SAR图像应用并没有有效利用方位向多散射中心的特征。
技术实现思路
本专利技术实施例提供了一种圆迹SAR图像序列特征提取方法，包括：步骤S1，联合应用似然比和恒虚警检测实现对圆迹SAR初始图像的方位向多散射中心的检测、消除和非相干累加，得到新的圆迹SAR图像序列，并获得所述新的圆迹SAR图像的似然比曲线和子孔径消除曲线；步骤S2，通过对所述似然比曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向似然比特征；以及步骤S3，通过对子孔径消除曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向子孔径消除特征。针对传统图像合成方法未完全利用由圆迹SAR成像中各子孔径视向不同而产生的丰富的和视向相关的信息，根据本专利技术的技术方案，利用通过似然比检测、恒虚警检测获得的圆迹SAR图像序列进行圆迹SAR图像的多视向特征提取。通过重组原始的子孔径图像序列获得一种物理意义明确且和观测对象所处方位无关的数据，从序列图像中提取出随散射中心变化的多视向特征，应用差异性似然比和恒虚警检测中生成的随视向变化的差异性似然比曲线和子孔径消除曲线，分别对差异性似然比曲线和子孔径消除曲线进行前向差分、编码和转码，并合并类别数，获得多视向似然比特征和多视向子孔径消除特征。附图说明图1示出了根据本专利技术的SAR图像序列特征提取方法的流程图；图2示出了根据本专利技术的提取SAR图像序列的差异性似然比特征的流程图；图3示出了根据本专利技术的提取圆迹SAR图像序列的子孔径消除特征的流程图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本专利技术作进一步的详细说明。图1示出了根据本专利技术的圆迹SAR图像序列特征提取方法的流程图。如图1所示，所述方法可以包括：步骤S1，联合应用似然比和恒虚警检测实现对圆迹SAR初始图像序列的方位向多散射中心的检测、消除和非相干累加，得到新的圆迹SAR图像序列，并获得新的圆迹SAR图像的似然比曲线和子孔径消除曲线；步骤S2，通过对似然比曲线进行前向差分、编码和转十进制码，并对编码后的特征进行合并类别数，得到多视向似然比特征；以及步骤S3，通过对子孔径消除曲线进行前向差分、编码和转十进制码，并对编码后的特征进行合并类别数，得到多视向子孔径消除特征。接下来将参考图2和图3来详细描述根据本专利技术的圆迹SAR图像序列特征提取方法。图2示出了根据本专利技术的提取SAR图像序列的差异性似然比特征的流程图，图3示出了根据本专利技术的提取圆迹SAR图像序列的子孔径消除特征的流程图。首先，本领域技术人员可以理解，可以根据本领域公知的各种方法联合应用似然比和恒虚警检测实现对圆迹SAR初始图像序列的方位向多散射中心的检测、消除和非相干累加，得到新的圆迹SAR图像序列，并获得新的圆迹SAR图像序列中每个像素点随序列号变化的似然比大小和子孔径序号，构成似然比曲线和子孔径消除曲线。例如可以，参见Ferro-Famil L，Reigber A，Pottier E，et al..Scene characterization using subaperture polarimetric SAR data[J].IEEE Transactions on Geoscience and Remote Sensing，2003，41(10)：2264-2276。如图2所示，步骤S2具体可以包括：在步骤S21，对似然比曲线进行前向差分。同一类散射目标的似然比曲线绝对值差异较大，这是由于似然比计算过程中主要依赖散射强度信息。散射强度容易受到多种因素的影响，例如天线的散射方向图、地形坡度以及目标本身的介电常数和几何结构。由此可见，似然比曲线处理应以曲线特征提取及归一化为主，似然比曲线的变化趋势(上升、降低或不变)可以反映目标的方位向散射结构。因此需结合似然比曲线的前向差分运算及状态编码来提取曲线特征并归一化。若Λ(i)代表似然比曲线，其前向差分值可以被定义为：ΔlogΛ(i)＝logΛ(i+1)-logΛ(i)其中1≤i≤N，N为图像序列长度。在步骤S22，对似然比曲线的前向差分值进行编码。由于同一类型地物的前向差分值在实际数据中的差异较大(和似然比相同，受到多种因素的共同作用)，但变化趋势表现出很强的模式性，且在模型分析中具有明确的含义。因此，将似然比曲线特征简化为变化趋势的三种状态：上升、降低或不变。为进一步应用于分类，可以将三种状态表达为三进制数Num，一组编码ki，i＝1，2，3…N-1，按照Num＝k1k2…组成三进制数，编码的位数可以根据应用需求来确定。其中，若ΔlogΛ(i)＝0，则编码ki＝0；若ΔlogΛ(i)＞0，则ki＝1；若ΔlogΛ(i)＜0，则ki＝2。编码的重要性及可信度从左至右依次降低的。因此，若在后期应用时发现过度分割的现象或希望减少类别数，可以在编码阶段调整尾数的位数，以此来舍弃不重要或可信度差的编码。该编码是不同类型的散射目标和其自身似然比曲线变化趋势唯一匹配的“身份”信息，如同提取了每个散射目标的指纹信息。然后，在步骤S23，对Num进行十进制转码。在得到指纹信息后，将多位Num进一步转化为一个十进制数字Numd。三进制转十进制从编码右侧开始，令似然比差异较小的类型所对应的十进制数值比较接近。变化次数多的和变化次数少的类型所对应的十进制数值相差较大。相同变化次数但在不同迭代周期中发生的也有差别。然后，在步骤S24，基于十进制编码Numd的直方图合并类别数目，从而获得多视向似然比特征。合并类别数目方法需根据应用需求决定。通过调整预设的直方图分辨间隔ΔL，可以调节特征合并类别过程中的敏感度。如果需要得到细小的特征差别，则令直方图间隔ΔL变小，敏感度调高。如果特征的类别合并结果有过度分类现象，则增大ΔL以降低敏感度。计算Numd和第l个类别之间的距离Lenl，l＝1，2…P，P是类别的数量。对距离进行排序得到最短距离Lenmin＝Lenx，则Numd被归类到第x个类别中。如图3所示，步骤S3可以具体包括：在步骤S31，对子孔径消除曲线进行前向差分。除似然比外，新的圆迹SAR图像序列还产生了消除子孔径标号序列。被消除的子孔径标号代表了该散射中心的后向散射方向图主要覆盖方向。由于试验中的消除子孔径标号在[0，9]之间，因此图中把方位向同散射中心的孔径差异统一设置为10。曲线的特征不是非常明显，很难从中直接分辨不同地物类型。基于前文对似然比曲线的分析，可以采用类似的思路来处理子孔径消除标号曲线。若假设子孔径的数量R是偶数，Ap(i)代表子孔径消除曲线，其前向差分值可以被定义为： Δ A p ( i 本文档来自技高网...

【技术保护点】
一种圆迹合成孔径雷达SAR图像序列特征提取方法，包括：步骤S1，联合应用似然比和恒虚警检测实现对圆迹SAR初始图像的方位向多散射中心的检测、消除和非相干累加，得到新的圆迹SAR图像序列，并获得所述新的圆迹SAR图像的似然比曲线Λ(i)和子孔径消除曲线Ap(i)；步骤S2，通过对所述似然比曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向似然比特征；以及步骤S3，通过对子孔径消除曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向子孔径消除特征。

【技术特征摘要】
1.一种圆迹合成孔径雷达SAR图像序列特征提取方法，包括：步骤S1，联合应用似然比和恒虚警检测实现对圆迹SAR初始图像的方位向多散射中心的检测、消除和非相干累加，得到新的圆迹SAR图像序列，并获得所述新的圆迹SAR图像的似然比曲线Λ(i)和子孔径消除曲线Ap(i)；步骤S2，通过对所述似然比曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向似然比特征；以及步骤S3，通过对子孔径消除曲线进行前向差分、编码和转码，并对编码后的特征进行合并类别数，得到多视向子孔径消除特征。2.根据权利要求1所述的方法，其中，所述步骤S2中对所述似然比曲线进行前向差分包括ΔlogΛ(i)＝logΛ(i+1)-logΛ(i)其中1≤i≤N，N为图像序列长度，ΔlogΛ(i)表示似然比曲线前向差分值。3.根据权利要求2所述的方法，其中，所述步骤S2中对所述似然比曲线进行编码包括将所述似然比曲线前向差分值编码为三进制数Num＝k1k2…，ki，i＝1，2，3…N-1，其中若ΔlogΛ(i)＝0，则编码ki＝0；若ΔlogΛ(i)＞0，则ki＝1；若ΔlogΛ(i)＜0，则ki＝2。4.根据权利要求3所述的方法，其中，所述步骤S2中对所述似然比曲线进行编码还包括：对Num进行十进制转码得到十进制数Numd，并基于十进制数Numd的直方图合并类别数目，从而获得多视向似然比特征。5.根据权利要求1所述的方法，其中，所述步骤S3中对所述子孔径消除曲线Ap(i)进行前向差分包括：在子孔径的数量R是偶数的情况下，将Ap(i)的前向差分值ΔAp(i)定义为： Δ A p ( j ) = | A p ( i + 1 ) - A p ( i ) | | A p ( i + 1 ...

【专利技术属性】
技术研发人员：洪文，李洋，薛斐腾，郭小洋，王建峰，郭胜龙，
申请(专利权)人：中国科学院电子学研究所，
类型：发明
国别省市：北京;11