[发明专利]一种高光谱图像端元个数自动估计的端元提取方法

权利要求书

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。