发明名称 基于对角子类判决分析的合成孔径雷达目标识别方法
摘要 本发明提出一种对角子类判决分析的合成孔径雷达目标识别方法,主要解决现有的合成孔径雷达目标识别性能差的问题。该方法的过程是:对原始图像进行自适应的阈值分割、形态学滤波、几何聚类操作和图像增强的预处理;对预处理后的每类目标采用二维快速全局k均值聚类算法进行最优的子类划分;用对角子类判决分析或者对角子类判决分析和二维子类判决分析找到最优的投影矩阵;将预处理后的训练和测试图像向投影矩阵投影,得到它们的特征矩阵;计算测试目标与每个训练目标的特征矩阵之间的欧氏距离,并采用最近邻准则确定测试目标的类别属性。仿真实验表明,本发明抑制背景杂波效果好、目标图像质量高、特征维数低的优点,可用于遥感系统中。
申请公布号 CN101526995B 申请公布日期 2011.06.29
申请号 CN200910020969.4 申请日期 2009.01.19
申请人 西安电子科技大学 发明人 刘宏伟;胡利平;苏洪涛;冯大政
分类号 G06K9/00(2006.01)I;G06K9/62(2006.01)I;G01S7/00(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;黎汉华
主权项 1.一种基于对角子类判决分析的合成孔径雷达目标识别方法,包括如下步骤:(1)对原始的训练和测试图像依次进行阈值分割、形态学滤波、几何聚类操作、图像增强和二维傅里叶变换的预处理;(2)采用二维快速全局k均值聚类算法对预处理后的每类目标图像进行最优的子类划分,并按如下步骤采用对角子类判决分析找到最优子类的最优投影矢量:(2a)将每类预处理后的合成孔径雷达训练目标图像采用二维全局k均值聚类算法找到最优的子类划分结果;(2b)将预处理后的合成孔径雷达训练图像沿对角方向将相应的图像像素取出来,重新排列后转化为相应的对角图像;(2c)根据上述的子类划分结果和相应的对角图像,定义如下公式:<img file="FSB00000474222600011.GIF" wi="1272" he="329" />其中,G<sub>B</sub>和G<sub>W</sub>分别表示所有训练图像对应对角图像的子类类间和子类类内散布矩阵,c为训练样本的类别数,d<sub>j</sub>表示第j(j=1,…,c)类目标包含的子类个数,d<sub>l</sub>表示第l(l=1,…,c)类目标包含的子类个数,<img file="FSB00000474222600012.GIF" wi="75" he="48" />表示第j类的第i个子类的均值,<img file="FSB00000474222600013.GIF" wi="77" he="47" />表示第l类的第h个子类的均值,<img file="FSB00000474222600014.GIF" wi="75" he="58" />表示<img file="FSB00000474222600015.GIF" wi="60" he="58" />的对角图像,<img file="FSB00000474222600016.GIF" wi="62" he="58" />表示第j类的第i个子类中的第k个样本,p<sub>j,i</sub>=N<sup>j,i</sup>/M为第j类的第i个子类的先验概率,p<sub>l,h</sub>=N<sup>l,h</sup>/M为第l类的第h个子类的先验概率,N<sup>j,i</sup>为第j类的第i个子类中的样本个数,N<sup>l,h</sup>为第l类的第h个子类中的样 本个数,M为所有训练样本的个数,R<sup>n×n</sup>表示子类类间和子类类内散布矩阵的大小,T表示矩阵转置;(2d)对矩阵G<sub>W</sub>求逆后与G<sub>B</sub>相乘,得到矩阵<img file="FSB00000474222600021.GIF" wi="145" he="60" />再对它进行特征分解,得到相应的特征值和特征向量,并取前d个较大特征值对应的特征向量w<sub>1</sub>,…,w<sub>d</sub>组成投影矩阵W=[w<sub>1</sub>,…,w<sub>d</sub>]∈R<sup>n×d</sup>,R<sup>m×d</sup>表示特征矩阵的维数;(3)将预处理后的训练和测试图像向最优投影矢量进行投影,得到训练和测试目标的特征矩阵:(3a)将第i个训练样本I<sub>i</sub>∈R<sup>m×n</sup>向投影矩阵投影,得到它的特征矩阵B<sub>i</sub>:<img file="FSB00000474222600022.GIF" wi="1501" he="101" />其中,<img file="FSB00000474222600023.GIF" wi="63" he="68" />表示I<sub>i</sub>∈R<sup>m×n</sup>向w<sub>k</sub>(k=1,2,…,d)投影得到的特征分量,<img file="FSB00000474222600024.GIF" wi="33" he="46" />为所有训练图像的均值图像,R<sup>m×d</sup>表示特征矩阵的维数,(3b)将测试样本I∈R<sup>m×n</sup>向投影矩阵投影,得到它的特征矩阵B为:<img file="FSB00000474222600025.GIF" wi="1176" he="126" />其中,y<sub>k</sub>表示I<sub>i</sub>∈R<sup>m×n</sup>向w<sub>k</sub>(k=1,2,…,d)投影得到的特征分量;(4)计算测试目标与每个训练目标的特征矩阵之间的欧氏距离,并采用最近邻分类器准则确定测试目标的类别属性。
地址 710071 陕西省西安市太白路2号