[0041] 以下结合附图对本发明作进一步说明。
[0042] 如图1所示,本发明包括以下步骤:
[0043] 步骤(1).预处理
[0044] 1.1方位角的标记
[0045] 由于同一种目标在不同方位角条件下的性质有很大的差别,因此在目标识别中应该考虑目标的方向性。估计出目标方位角可极大地减少目标识别所需的特征个数,从而减少识别时间并提高识别性能。采用手工方法标记方位角,记录目标与标准方位角(90度)之间的夹角。并将标记的方位角与图像的索引建立一种映射关系,即将方位角和图像的索引一一对应,以便于通过索引来找到图像的方位角。
[0046] 1.2去噪、滤波
[0047] SAR图像中的相干斑噪声、背景、杂波和方位角变化等都为目标的识别带来困难,这就要求在提取特征之前进行预处理使得图像中的信息更接近真实值。主要包括:SAR图像去噪(抑制SAR图像相干斑噪声),SAR图像分辨率的改善,SAR图像分割特征。对图像进行了裁剪、小波变换和均值滤波处理,裁剪即读取以图像中心点为中心的61×61图像区域,小波变换中使用的小波是Haar小波,均值滤波的模板是[33]。
[0048] 步骤(2).分类器模型设计
[0049] 2.1建立压缩感知模型和基于压缩感知分类
[0050] 2.1.1模型建立
[0051] 在训练样本中,选取第i类目标的ni张图像,将第j(j=1,…,ni)张图像进行预处理,然后将图像的信息存放在一个列向量里来构成Ai(Ai表示由第i类目标所有图像的信息)的第j列vi,j(vi,j表示第i类目标第j张图像的信息),得到一个过完备字典
[0052] A=[A1;...;Ai;...;Ak]=[v1,1,...,v1,n1;....;vi,1,...,vi,j,...,vi,ni;...;vk,1,...,vk,nk](m×n)。
[0053] 采用MSTAR数据库中的图像,由训练样本直接构成字典时,一幅图像用961×1的列向量(vi,j)表示,字典A的大小为961×1622。
[0054] 对于第i类目标的测试样本y表示为
[0055] 式(1);
[0056] 即
[0057] y=Ax0 式(2);
[0058] 其中 标量 j=1,…,ni。x0是稀疏系数,x0中只有与第i类目标相对应的值为非零,其它值均为零。
[0059] 用l1范数求解x0,用数学公式表达如下:
[0060] (P1) min||x0||1 subject to y=Ax0 式(3);
[0061] 上式为压缩感知的优化重建问题,可求出稀疏解x0,求解出的x0的大小为1622×1。
[0062] 2.1.2基于压缩感知分类
[0063] 由于相干斑噪声、方位角等因素的影响,x0中非对应于测试图像类别的系数值也是非零的,这些值相对于x0中对应于测试图像类别的系数值是很小的,因此可采用寻找最小重构误差的方法进行目标类别判别,该方法也可视为最近邻法(Nearest Neighbor,NN)。对于第i类目标,定义函数 对于 是系数向量,其中
x0中仅与第i类目标相对应的系数保持不变,其它值均为零。定义残差为
[0064] ri(y)=||y-Aδi(x0)||2 式(4);
[0065] 残差值最小意味着测试图像与样本最接近,通过求解的残差最小值来判断类别。则测试图像y的所在类别s可由式子(5)判断:
[0066] 式(5);
[0067] 2.2选择支持向量机模型参数
[0068] 支持向量机的模型训练,最重要的部分是选定核函数。目前使用最普遍的核函数是径向基核函数,这是因为径向基核函数对应的特征空间是无穷维的,有限的样本在该特征空间中肯定是线性可分的。径向基核函数参数的选择主要是C和γ。C是惩罚因子,表示对错误的惩罚程度,我们可以通过C来控制训练精度,一般情况下取值为[0,50];γ主要影响着样本数据在高维特征空间中分布的复杂程度,γ的取值范围是[0,1]。C和γ的最优结果由交叉验证和网格参数寻优方法训练得到,最优结果为C取值32和γ取值0.5。
[0069] 步骤(3).姿态矫正
[0070] 用SVM方法进行SAR图像目标识别时,方位角的校正对识别结果很重要;但是对于测试图像来说,它的信息全部都是未知的,无法校正方位角。对于CS方法的优化求解数据,稀疏解x0不仅能够得到识别结果,还可以得到与测试图像最相似的训练样本图像的姿态信息,可用来估计测试图像的姿态信息和校正方位角。
[0071] 利用训练样本构造字典后,利用压缩感知模型进行优化求解时,可根据获得稀疏解x0判断与之最相近的训练样本,即稀疏解x0的最大值位置t处表明测试样本与该训练样本最相似:
[0072] 式(6);
[0073] 此时根据训练样本的姿态信息估计测试样本的姿态信息,从而对测试样本进行姿态矫正。
[0074] 步骤(4).决策级融合
[0075] 采用投票法(Majority Vote)对CS方法和SVM方法(校正方位角、未校正方位角)识别的三种结果进行决策级融合,即
[0076] r=MajorityVote(rcs,rsvm1,rsvm2) 式(7);
[0077] 式中rcs为未经过姿态校正的基于CS方法的识别结果,rsvm1为未经过姿态校正的基于SVM方法的识别结果,rsvm2为经过姿态校正的基于SVM方法的识别结果。当出现三种结果不一致时,根据每种算法对训练样本的识别准确率作为先验知识进行判别,即
[0078] 式(8);
[0079] 式中 为采用CS方法进行识别时识别结果为第i类的概率, 为基于未经校正的样本采用SVM方法进行识别时识别结果为第j类的概率, 为基于经校正的样本采用SVM方法进行识别时识别结果为第k类的概率。