[发明专利]联合解混及自适应端元提取的高光谱图像分类方法

权利要求书

1.一种联合解混及自适应端元提取的高光谱图像分类方法，其特征在于包括以下步骤：

第一步，采用基于统计理论的最大似然分类算法对原始高光谱图像进行粗分类；当样本数量N满足N>1000，每类随机选择5‰的像元作为训练样本，100<N<1000则每类随机选择5%像元作为训练样本，否则随机选择5个像元作为训练样本；对原始高光谱图像数据采用主成分分析方法实现数据降维；

步骤二、采用混淆矩阵作为筛选类别端元集的依据，通过迭代不断更新端元集；将GIS数据计算类别平均光谱作为端元光谱；设图像共有N个像元，像元集合为X=[x₁,x₂,...,x_N]，图像有K个波段像元的类标签集合为Y=[y₁,y₂,...,y_N]，共有M个类别，因此像元x_i的类标签y_i满足图像中任意一个类别m_i的端元集提取过程如下：

（1）计算各个类别的平均光谱作为端元光谱e_i，获得图像的端元集Ε=[e₁,e₂,...,e_M]；其中，N_i表示第i类的像元个数，表示第i类的第j个像元；

ei=1NiΣj=1Nixj(i)---(1)]]>

（2）利用混淆矩阵统计第i类地物m_i错分到其他类别的像元个数；设定阈值λ，满足.0＜λ＜＜1，对于第j类地物，若第i类像元错分到第j类的像元个数C_ij与第i类分对的像元个数C_ii之比大于阈值λ，则将第j类的端元加入到第i类地物的端元集中，Β_i为图像端元集Ε的子集，S_i是第i类地物可能包含的地物种类的个数；

T={t|C_it≠0,C_it≥λC_ii,t∈Μ}                         （2）

Β_i=Ε_T,i∈Μ                       （3）

其中，T表示混淆矩阵中满足条件的类别号的集合；Β_i表示由相应类别号从图像端元集Ε中选择对应的端元组成第i类地物的端元集；λ是阈值；

（3）根据以上过程获得各个类别相应的端元集Β_k，针对第k类地物中的混合像元，利用其端元集合Β_k进行线性光谱解混，获得混合像元的丰度值信息；

步骤三、采用线性光谱解混模型估计混合像元的丰度值信息；实际的高光谱数据丰度值系数A=[a₁,a₂,...,a_M]需要满足采用全约束的最小二乘方法条件的混合像元丰度值信息；

步骤四、利用基于多元逻辑回归的概率分类器估计像元的概率类标签；其中，混合像元采用其丰度值信息优化MLR的回归系数ω；具体的回归参数估计过程如下：

（1）x_i的类标签为当y_i中仅含有一个值时，x_i为纯净像元，否则为混合像元；MLR模型的回归系数为利用最大后验概率法估计回归系数ω，如公式（4）和公式（5）；使用标记样本集D_L={(y₁,x₁),(y₂,x₂),...,(y_L,x_L)}作为训练样本，其中，包含的混合像元的类标签是其丰度值；

ω^=argmaxωl(ω)---(4)]]>

l(ω)=logΠi=1Lp(yi|xi,ω)=Σi=1L(Σt=1Myi(t)ω(k)Txi-logΣt=1Mexp(ω(k)Txi))---(5)]]>

（2）根据回归系数建立多元逻辑回归模型，用于计算像元概率类标签；

P(yit=1|xi,ω^)=exp(ω^(i)Txi)Σj=1Mexp(ω^(j)Txj)---(6)]]>

步骤五、利用MLR分类器对待分类像元分配概率类标签，具体分类步骤如下：

（1）将待分类像元x_j输入MLR分类器，估计像元属于第i类地物的概率P(yji=1|xj,ω);]]>

（2）选择概率最大的类别作为像元的所属类别；

xj∈mt,s.t.P(yjt=1|xj,ω^)=argmax{P(yji=1|xj,ω^)},i=1,2,...,M---(7)]]>

步骤六、根据以上分类结果获得混淆矩阵，重复步骤二到步骤六；通过更新各类端元集，实现MLR分类器的优化。