[發明專利]一種高光譜圖像端元個數自動估計的端元提取方法

有效

• 申請號: 201310601834.3
• 申請日: 2013-11-25
• 公開（公告）號: CN103617424A
• 公開（公告）日: 2014-03-05
• 發明（設計）人: 許毅平;田巖;譚志敏
• 申請（專利權）人: 華中科技大學
• 主分類號: G06K9/46
• 分類號: G06K9/46
• 代理公司: 華中科技大學專利中心 42201
• 代理人: 朱仁玲

1.一種高光譜圖像端元個數自動估計的端元提取方法，包括以下步驟：

S1、使用PPI算法計算每個像元的PPI值，將PPI值大于0的像元作為初始候選端元，并用矩陣E_L×N表示，其中L為數據維數，N為候選端元的個數；

S2、對初始候選端元按PPI值從大到小進行排序，根據候選端元間的相關性逐個剔除端元，并獲得候選端元集的光譜重構精度曲線R和獨立性曲線I；

S3、根據光譜重構精度曲線和獨立性曲線隨端元個數減少的變化率，確定端元個數的取值范圍；

S4、在步驟S3確定的端元個數取值范圍內，利用SGA算法提取不同端元個數所對應的候選端元集，并計算這些候選端元集的光譜重構精度曲線，選擇曲線變化最顯著的位置值作為最終的端元個數值，并將該個數值對應應的候選端元集作為最終的端元集。

2.根據權利要求1所述的端元提取，其中，步驟S2具體為：

（a）令p＝N，S_p＝E_L×p；

（b）將初始候選像元根據其PPI值從大到小排序，結果存入S_p中；

（c）計算p個候選端元兩兩間的光譜角距離，用距離矩陣C_p×p表示，該矩陣為下三角矩陣，對角線上的值設為0；

（d）根據式(1)得到I(p)，根據式(2)～(4)得到R(p)；

I(p)＝c_i,j??????????(1)

αp=(SpT·Sp)-1·SpT·γ---(2)]]>

γp^=Sp·αp---(3)]]>

R(p)=cos<γ,γ^p>=γ·γ^|γ||γ^|----(4)]]>

其中，c_i,j為距離矩陣C_p×p中的最大值，γ為圖像數據的均值向量；

（e）刪除矩陣C_p×p中的第i行和第j列；

（f）令p＝p-1；如果p＞0，則跳至（d）繼續，否則退出S2。