[发明专利]基于几何空谱结构信息的高光谱混合像元分解方法

权利要求书

1.基于几何空谱结构信息的高光谱混合像元分解方法，其特征在于：包括以下步骤：

(1)输入高光谱数据，图像大小为a×b像素,并预处理得到高光谱数据L为处理的波段个数，N为像素点个数，N＝a×b；

(2)用VD估计法，计算纯端元个数p；

(3)提取高光谱图像的轮廓信息，通过PCA降维提取出第一主成分量提出对Y用Robert算子提取场景轮廓图N＝a×b；

(4)根据边缘信息、位置信息计算高光谱数据之间的几何空间距离度量提出此空间距离的公式；

大小为m×m像素的局部空间窗LW内若中心点为(m_i,n_i)，在局部窗LW内，通过轮廓线判定中心点(m_i,n_i)与周围任一点(m_j,n_j)是否处在同一区域内，局部空间窗LW内任一点(m_j,n_j)与中心点的几何空间距离度量为：

此处为两个光谱向量之间的差异；

(5)根据数据光谱统计信息计算谱间距离度量提出此谱间距离的公式；

在光谱维上，对数据建立最近邻图，则两个光谱向量x_i,x_j之间的几何谱间距离度量定义为：

di,je=e-||xi-xj||22,ifxi∈NB(xj)e-γ||xi-xj||22,ifxi∉NB(xj)]]>

此处γ＝[(m_i-m_j)+(n_i-n_j)]²为两光谱向量位置距离，NB(x_j)为向量x_j的k个最近邻点集合；

(6)构造新的NMF算法框架

由步骤(4)和步骤(5)中的几何空谱距离度量和根据图论中如果原数据x_i,x_j之间相似则其分解后的丰度分布向量之间必定相似的原则构造空谱拉普拉斯流形正则：

||f||Ma=Σi,j=1N||si-sj||2/di,ja=STLaS]]>和||f||Me=Σi,j=1N||si-sj||2/di,je=STLeS]]>

此处的L^a，L^e邻接图的拉普拉斯矩阵，和是每个节点相关的对角矩阵,d~iie=Σj=1n1/dije,d~iia=Σj=1n1/dija;]]>

由此构造新的非负矩阵分解框架式，为解混后得出的纯端元特征矩阵，则对应的丰度分布比例矩阵为

f(M,S)=fe(M,S)+||f||Ma+||f||Me=12||X-MS||22+tr(SLaST)+tr(SLeST)]]>

(7)通过新的非负矩阵分解模型输出端元光谱曲线矩阵和对应的每个端元分布丰度图矩阵。