[发明专利]一种基于压缩感知的高光谱混合像元分解的方法

摘要

一种基于压缩感知的高光谱混合像元分解的方法，涉及压缩感知领域与高光谱遥感领域。解决了现有采用传统的奈奎斯特采样定理对高光谱图像数据采集时，混合像元分解速度慢的问题。该方法首先输入观测矩阵Φ和压缩观测矩阵Y，利用压缩感知理论建立光谱混合模型Y＝ΦXT＝Φ(AS)T，其次对端元丰度矩阵S的估计值和端元光谱矩阵A的估计值进行迭代处理，如果相邻两次获得的端元光谱矩阵A的估计值中对应的每个元素的绝对值之差小于0.1时，则停止迭代，输出端元丰度矩阵，完成对高光谱混合像元的分解，否则继续进行迭代处理。主要用于对高光谱混合像元分解。

主权项

一种基于压缩感知的高光谱混合像元分解的方法，其特征在于，该方法为，步骤一，输入观测矩阵Φ和压缩观测矩阵Y，利用压缩感知理论，建立光谱混合模型：Y＝ΦXT＝Φ(AS)T      (1)Φ∈RM×N为M×N的观测矩阵，R为实数，X∈RL×N为L×N的混合像元光谱矩阵，Y∈RM×L为M×L的压缩观测矩阵，S∈RP×N为P×N的端元丰度矩阵，A∈RL×P为L×P的端元光谱矩阵，步骤二，初始化，随机选取一个端元光谱矩阵A作为端元光谱矩阵A的估计值且为L×P的矩阵，令端元丰度矩阵S的估计值：S^=0...0...00P×N---(2)]]>其中，为P×N的矩阵，步骤三，令迭代次数变量t的初始值为1；步骤四，对端元丰度矩阵S的估计值和端元光谱矩阵A的估计值进行迭代处理；步骤五，在公式(1)两边同时乘以的伪逆，则公式(1)变形为Y1＝ΦST      (3)，步骤六，设小波基为Θ，Σ为端元丰度的稀疏系数，则公式(3)变形为Y1＝Φ(ΘΣ)T     (4)，采用BP算法求解下式：arg min||Σ||1 s.t. Y1＝Φ(ΘΣ)T     (5)，获得端元丰度矩阵S的估计值S^=ΘΣ---(6),]]>步骤七，在概率密度函数为P(Sk)∝αexp(α|Sk|)的条件下，归一化端元光谱矩阵A的估计值的列向量，更新端元光谱矩阵A的估计值A^=A1+λ{-A1(BS^T+I)}---(7),]]>其中，A1为上一次迭代得到的端元光谱矩阵的估计值，λ为迭代步长，I为单位矩阵，B为向量集合，且B＝{B1，B2，B3，......Bk}，为第k个源信号的先验概率分布，k为整数，α为大于零的实数，Sk表示第k个端元丰度矩阵S，表示第k个端元丰度矩阵S的估计值，步骤八，如果相邻两次获得的端元光谱矩阵A的估计值中对应的每个元素的绝对值之差小于0.1时，则停止迭代，执行步骤九，否则，令t＝t+1，且返回步骤四；步骤九，输出端元丰度矩阵完成对高光谱混合像元的分解。