技术领域
本发明涉及图像处理、模式识别等领域，尤其涉及合成孔径雷达 图像处理与目标识别领域。
背景技术
合成孔径雷达（Synthetic Aperture Radar，SAR）图像目标识别是 SAR图像解译和处理的关键步骤，具有重要的应用价值。SAR图像目 标识别是对获取的含有目标的SAR图像进行目标检测、定位并判别出 目标的种类和型号。由于阴影效应，信号与背景的交互，三维场景在 斜平面的投影以及其他雷达横截面对方位角敏感的因素，SAR图像对 目标方位角的变化十分敏感。因而，SAR图像中目标的识别能力也会 随着方位角的变化而改变。所以，可利用同目标不同方位角下的多视 SAR图像来提高目标识别的性能。利用多视SAR图像进行目标识别 可在目标识别框架的某些层次进行。如在数据层，可使用主成分分析、 离散小波分解等数据融合算法将多视SAR图像进行数据层融合；如在 决策层，可使用贝叶斯决策、DS决策将多视SAR图像的类别决策进 行决策层的融合。这些方法均获得了较单幅SAR图像目标识别更高的 识别率。目前，绝大多数利用同目标不同方位角下的多视SAR图像进 行目标识别的方法都仅在目标识别框架的一个层次进行，这使得目标 识别率的提高受到了限制。
为了进一步提高SAR图像目标识别性能，利用多视SAR图像进 行目标识别可同时在目标识别框架的两个层次——数据层和决策层进 行，并将两个层次获得的识别率收益相互结合，生成最终类别决策， 从而提高目标识别的正确率。
发明内容
本发明利用同目标不同方位角下的多视合成孔径雷达目标图像进 行目标识别，在数据层用凸集投影超分辨率重构算法将多视合成孔径 雷达图像重构出一幅高分辨率雷达图像，在决策层用带权重的贝叶斯 决策融合算法进行决策融合，将数据层和决策层两个层次获得的识别 率收益相互结合，提高合成孔径雷达图像目标识别的正确率。
本发明解决其技术问题所采用的技术方案是：
多层次结合的多视合成孔径雷达图像目标识别方法，其特征在于： 包括以下步骤：
步骤1，对采集到的N幅未知目标的同类不同方位角下的多视合 成孔径雷达目标图像，k=1,2,…N进行图像预处理；
步骤2，对预处理后的每一幅合成孔径雷达图像用小波分解和主 成分分析进行特征提取，提取得到代表该幅图像的多维特征矢量；
步骤3，用经过训练的支持多目标分类的支持向量机对步骤2得 到特征矢量进行分类，得到该幅图像xk，k=1,2,…N属于某个目标类 别q的后验概率p(q|xk)；
步骤4，对预处理后的多视合成孔径雷达目标图像xk，k=1,2,…N 在数据层用凸集投影超分辨率重构算法重构一幅高分辨率雷达图像 xc，xc的分辨率是xk，k=1,2,…N的2倍；
步骤5，对重构后的高分辨率图像用小波分解和主成分分析进行 特征提取，提取得到代表该重构图像的多维特征矢量；
步骤6，用经过训练的支持多目标分类的支持向量机对步骤5得 到的特征矢量进行分类，得到该重构图像xc属于某个目标类别q的后 验概率p(q|xc)；
步骤7，将步骤3得到的每幅图像的后验概率p(q|xk)，k=1,2,…N 和步骤6得到的重构图像的后验概率p(q|xc)用带权重的贝叶斯决策 融合方法进行决策层融合，得到该多视合成孔径雷达目标图像的目标 类别。
进一步，所述步骤1中，预处理包括图像截取，噪声处理，幅度 值归一化处理，以及合成孔径雷达图像中目标方位角估计。
进一步，所述步骤4中，在数据层使用的凸集投影超分辨率重构 算法包括：选择一幅合成孔径雷达图像作为参考图像，采用插值算法 将其构造成为分辨率为原分辨率2倍的参考帧，作为超分辨率图像的 初始估计；设置迭代次数，对每一幅低分辨率合成孔径雷达目标图像 进行运动补偿，残差计算和像素值修正，得到一幅分辨率为原分辨率 2倍的凸集投影超分辨率重构图像。
进一步，所述步骤7中，带权重的贝叶斯决策融合方法包括：计 算
Q为目标类别总数。
本发明的技术构思为：该方法的一个递推过程包括以下基本步 骤：
1）、对采集到的不同方位角下的多视合成孔径雷达目标图像进行 图像预处理；
2）、用小波分解和主成分分析方法对预处理后的合成孔径雷达图 像进行特征提取；
3）、用支持向量机对特征分类，得到图像属于某个类别的后验概 率；
4）、对预处理后的多视合成孔径雷达目标图像在数据层用凸集投 影超分辨率重构算法重构一幅高分辨率雷达图像；
5）、对重构后的高分辨率图像用小波分解和主成分分析进行特征 提取；
6）、用支持向量机对上一步骤提取得到的特征进行分类，得到重 构图像属于某个类别的后验概率；
7）、用带权重的贝叶斯决策融合方法对单图像的后验概率和重构 图像的后验概率进行决策层融合，得到不同方位角下的多视合成孔径 雷达目标图像的类别。
本发明的有益效果主要表现在：同时在目标识别框架的两个层次 ——数据层和决策层利用同目标不同方位角下的多视合成孔径雷达图 像进行目标识别，并将两个层次获得的识别率收益相互结合，提高了 合成孔径雷达图像目标识别的正确率，有利于提高合成孔径雷达图像 解译的正确率。
附图说明
图1为本发明的多层次结合的多视合成孔径雷达图像目标识别方 法的流程图。
图2为凸集投影超分辨率重构过程示意图。
图3为本发明方法和其他对比方法在不同方位角间隔下得到的目 标正确识别率。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
参照图1，多层次结合的多视合成孔径雷达图像目标识别方法， 包括以下步骤：
步骤1，对采集到的N幅未知目标的同类不同方位角下的多视合 成孔径雷达目标图像xk，k=1,2,…N进行图像预处理；
步骤2，对预处理后的每一幅合成孔径雷达图像用小波分解和主 成分分析进行特征提取，提取得到代表该幅图像的多维特征矢量；
步骤3，用经过训练的支持多目标分类的支持向量机对步骤2得 到特征矢量进行分类，得到该幅图像xk，k=1,2,…N属于某个目标类 别q的后验概率p(q|xk)；
步骤4，对预处理后的多视合成孔径雷达目标图像xk，k=1,2,…N 在数据层用凸集投影超分辨率重构算法重构一幅高分辨率雷达图像 xc，xc的分辨率是xk，k=1,2,…N的2倍；
步骤5，对重构后的高分辨率图像用小波分解和主成分分析进行 特征提取，提取得到代表该重构图像的多维特征矢量；
步骤6，用经过训练的支持多目标分类的支持向量机对步骤5得 到的特征矢量进行分类，得到该重构图像xc属于某个目标类别q的后 验概率p(q|xc)；
步骤7，将步骤3得到的每幅图像的后验概率p(q|xk)，k=1,2,…N 和步骤6得到的重构图像的后验概率p(q|xc)用带权重的贝叶斯决策 融合方法进行决策层融合，得到该多视合成孔径雷达目标图像的目标 类别。
进一步，所述步骤1中，预处理包括图像截取，噪声处理，幅度 值归一化处理，以及合成孔径雷达图像中目标方位角估计。
进一步，所述步骤4中，在数据层使用的凸集投影超分辨率重构 算法包括：选择一幅合成孔径雷达图像作为参考图像，采用插值算法 将其构造成为分辨率为原分辨率2倍的参考帧，作为超分辨率图像的 初始估计；设置迭代次数，对每一幅低分辨率合成孔径雷达目标图像 进行运动补偿，残差计算和像素值修正，得到一幅分辨率为原分辨率 2倍的凸集投影超分辨率重构图像。
进一步，所述步骤7中，带权重的贝叶斯决策融合方法包括：计 算
Q为目标类别总数。
本实施例采用美国DARPA/AFRL MSTAR项目组提供的实测 SAR地面静止军用目标数据，将本发明方法和单幅图像的SAR目标 识别方法以及仅使用数据层方法，仅使用决策层方法进行比较。
MSTAR目标数据包含俯视角分别为15°和17°的SAR目标图像数 据，本实施例选取俯视角为17°的图像数据为训练样本，俯视角为15° 的图像数据为测试样本，因而，本实施例的训练样本数为2747，测试 样本数为3203，具体训练样本和测试样本可见表1。
表1
对采集到的N幅未知目标的同类不同方位角下的多视合成孔径 雷达目标图像xk，k=1,2,…N进行图像预处理。本实施例中N=4。 MSTAR数据库中的SAR目标图像大于或等于128×128像素，目标主 要集中在小于64×64像素的区域中，其余区域均为杂波，杂波区域过 大将不必要地增加后续处理过程的计算量，因此，需要对原始SAR目 标图像进行截取，形成尺寸为64×64像素的新目标图像。原始SAR 图像中，存在着乘性噪声，需要在识别前将目标图像的幅度值作对数 变换，使图像中的乘性噪声转换为加性噪声。本实施例针对MSTAR 数据库的十类地面静止军用目标数据进行目标识别，在原始SAR图像 中，不同类别的目标图像的幅度值分布差异较大，因此需要对图像幅 度值进行归一化处理，使归一化后图像中各像素的幅度值的均值为0， 标准差为1。本实施例中SAR图像目标的方位角由MSTAR数据提供。
用小波分解和主成分分析（Principal Component Analysis,PCA） 提取每一单幅SAR图像的特征。对每一幅图像，首先用小波分解进行 三层二维离散小波分解，提取三层分解后的低频子带图像，并将其表 示为一个64维的列向量。将10类目标的训练样本得到的所有64维列 向量按列方式排列组成数据矩阵X，求矩阵X的自相关阵C=E[X·XT]， 将C作特征值分解，得到特征值λi，i=1,2,…,64。对λi进行归一化，即 
将
从大到小依次排序，选择
中最大30个值对应的特征 向量为基向量构造投影子空间W=[v1,v2,…,v30]T，将目标图像x对投影 子空间投影e=W·x，可得到表征该目标图像的30维特征向量e。
用经过训练的支持10种种类目标分类的支持向量机（Support Vector Machine，SVM）对该30维的特征矢量e分类，得到该幅图像xk， k=1,2,3,4属于某个目标类别q，q=1,2,…10的后验概率p(q|xk)。
对预处理后的多视合成孔径雷达目标图像xk，k=1,2,3,4在数据层 用凸集投影（Projections onto Convex Sets,POCS）超分辨率重构算法 重构，生成分辨率为原分辨率2倍的超分辨率图像xc。选择该4幅合 成孔径雷达目标图像中的其中一幅作为参考图像，采用插值算法将其 构造成为分辨率为原分辨率2倍的参考帧，作为超分辨率图像的初始 估计；设置凸集投影超分辨率重构算法的迭代次数；对每一幅低分辨 率合成孔径雷达目标图像进行运动补偿、残差计算和像素值修正，得 到一幅分辨率为原分辨率2倍的凸集投影超分辨率重构图像。凸集投 影超分辨率重构过程如图2所示，其中图2(a)、2(b)、2(c)和2(d) 分别是目标BMP2_c21在俯视角17°和方位角101°、103°、105°和106° 下的SAR图像，图2(e)是通过双线性插值法得到的高分辨率图像，图 2(f)是经过20次迭代后得到的高分辨率图像。该4幅合成孔径雷达目 标图像的相邻方位角间隔应小于30°，才可获得较未进行超分辨率重 构的目标识别方法更高的目标正确识别率，因而，本实施例的方位角 间隔均选择小于30°。
对重构后的高分辨率图像用小波分解和主成分分析进行特征提 取，提取过程与单图像特征提取过程类似，这里不再赘述。最终提取 得到代表该重构图像的200维特征矢量。
用经过训练的支持10种种类目标分类的支持向量机对该200维 特征矢量进行分类，得到该重构图像xc属于某个目标类别q， q=1,2,…10的后验概率p(q|xc)。
将每幅图像经过支持向量机得到的后验概率p(q|xk)和重构图像 经过支持向量机得到的后验概率p(q|xc)用带权重的贝叶斯决策融合 方法进行决策层融合。计算
其中， p(q|xk)是单幅图像经过支持向量机得到的后验概率，p(q|xc)是重构 图像经过支持向量机得到的后验概率，wk是单幅图像后验概率的权 重，这里wk=1,k=1,2,3,4，wc是重构图像后验概率的权重，由于重构 图像包含了4幅SAR图像的信息，所以1≤wc≤4，本实施例中wc=2； 最终的类别决策为：
Q为目标类别数，这里Q=10。从而 识别出该多视合成孔径雷达目标图像的目标类别。
将本发明方法得到的目标正确识别率与单幅图像的SAR目标识 别方法以及仅使用数据层方法，仅使用决策层方法进行比较，比较结 果如图3所示。由图3可见，本发明方法得到的目标正确识别率高于 单幅图像的目标正确识别率，同时高于仅使用数据层方法和仅使用决 策层方法得到的目标正确识别率；且本发明方法得到的识别率对于多 视SAR图像的相邻方位角间隔的敏感性更低，也就是说，本发明方法 对于方位角间隔的变化更具有鲁棒性。
综上所述，本发明的多层次结合的多视合成孔径雷达图像目标识 别方法，在目标识别框架的两个层次——数据层和决策层同时利用同 目标不同方位角下的多视合成孔径雷达图像进行目标识别，并将两个 层次获得的识别率收益相互结合，提高了合成孔径雷达图像目标识别 的正确率，有利于提高合成孔径雷达图像解译的正确率。
显而易见，在不偏离本发明的真实精神和范围的前提下，在此描 述的本发明可以有许多变化。因此，所有对于本领域技术人员来说显 而易见的改变，都应包括在本权利要求书所涵盖的范围之内。本发明 所要求保护的范围仅由所述的权利要求书进行限定。
