[发明专利]基于谱聚类的极化SAR图像分类方法

权利要求书

1.一种基于谱聚类的极化SAR图像分类方法，包括如下步骤：

(1)对待分类的大小为R×Q的极化SAR图像进行滤波，去除斑点噪声；

(2)对滤波后的极化SAR图像每个像素点的相干矩阵T进行Cloude分解，提取每个像素的散射熵H特征，得到整幅图像的关于散射熵H的特征空间；

(3)用Mean Shift算法对获取到的特征空间进行分割，得到M个区域；

(4)在已获得的M个区域上，选取每个区域的典型代表点作为新的像素点Y_δ,δ＝1，...,M，得到M个新像素点，将这M个新像素点映射为具有M个节点的全连接图，并对这个全连接图进行谱聚类。

(4a)根据相似性准则，构造新像素点Y_δ的邻接矩阵A，将邻接矩阵A中每一个元素A_ij定义为：

Aij=exp(-d2(Ti,Tj)2σ2)i≠j0i=j]]>

其中A_ij是邻接矩阵A中第i行第j列的元素，i=1,...,M，j=1,…,M，σ为尺度参数，d(T_i,T_j)为新像素点Y_i和Y_j的相似性度量，T_i和T_j分别是新像素点Y_i和Y_j的相干矩阵，d(T_i,T_j)定义如下：

d(Ti,Tj)=12tr(Ti×Tj-1+Tj×Ti-1)-q]]>

其中q为常数，取值为q=3，tr(·)是矩阵的迹；

(4b)根据邻接矩阵A，构造规范的拉普拉斯矩阵L：

L=D-12AD-12]]>

其中D是对角矩阵，对角矩阵D对角线上的每一个元素为

(4c)对拉普拉斯矩阵L进行特征值分解，得到前k个最大特征值对应的特征向量，形成特征向量矩阵X＝[x₁,x₂,…,x_k]，k为分类类别数；

(4d)根据特征向量矩阵X得到规范化矩阵V：

规范化矩阵V中每一个元素V_θ,Z定义为：

Vθ,Z=Xθ,Z(ΣZ=1kXθ,Z2)-12]]>

其中X_θ,Z是特征向量矩阵X第θ行第Z列的元素，θ=1,…，M，Z=1,…,k；

(4e)用k-means算法对规范化矩阵V的行向量进行聚类；

(4f)对规范化矩阵V的行向量聚类后，当规范化矩阵V的第θ行为第c类时，将新像素点Y_δ标记为第c类，δ=1,…，M，c=1,…,k。

(5)在M个区域上，将由新像素点Yδ所代表的区域标记为与新像素点Y_δ相同的类别，完成对整幅图像的预分类。

(6)对预分类得到的整幅图像用能反映极化SAR分布特性的Wishart分类器进行迭代分类，得到更为准确的分类结果。