[发明专利]一种基于双线性混合模型的高光谱图像非线性解混方法

权利要求书

1.一种基于双线性混合模型的高光谱图像非线性解混方法，其特征在于：在双线性混合模型的基础上，引入一个非线性端点使高光谱遥感图像的非线性解混问题投影转换为简单的线性解混，然后采用线性解混算法在修正投影位置的同时迭代地求解该问题；具体步骤如下：

(1)确定双线性混合模型的非线性端点

基于不同的光谱混合模型，高光谱遥感图像中的每个像元表示为：

式(1)到式(4)分别对应于线性混合模型LMM和三种双线性混合模型BMM：FM，GBM及PPNM；其中，表示端元矩阵，r是端元数，s_i,j是端元a_i在像元x_j中的丰度，s_k,j是端元a_k在像元x_j中的丰度，ε_j是噪声误差；丰度皆满足非负及“和为1”约束：

非线性参数γ_i,k,j，0≤γ_i,k,j≤1，与b_j用于调控像元中的非线性影响程度；

首先计算双线性混合模型中的r个非线性混合中点，然后对数据降维以确定r个超平面，最后计算这些超平面的交点，即得到非线性端点，具体流程为：

输入高光谱数据和端元矩阵

步骤1.1：FM和GBM按式(6)而PPNM按式(7)分别计算各自对应的r个非线性混合中点ω₁,...,ω_r-1,ω_r，即各端元丰度值相等：

步骤1.2：利用PCA算法将数据降到r维，确定顶点集{ω₁,a₂,...,a_r}，{a₁,ω₂,...,a_r}，…，{a₁,a₂,…,ω_r}在该空间中所张成的r个r-1维超平面H₁…H_r；

步骤1.3：求解线性方程组(8)，非线性端点p即为其唯一解：

其中，U^T为输入的高光谱数据的数据前r个特征向量构成的矩阵，det表示行列式；

输出非线性端点p；

(2)修正投影位置并迭代估计丰度

先利用p点计算投影重心坐标s′_j，然后代入模型计算所对应的非线性混合部分，再确定使像元到其所在直线距离最短的系数，从而得到修正后的新投影位置，最后用FCLS算法求解新投影位置的丰度，并以此迭代进行直到收敛，具体流程为：

输入：高光谱数据和端元矩阵非线性端点p；

输出：丰度矩阵

步骤2.1：以最小二乘法求解优化问题(9)得投影重心坐标s′_j

这里，

步骤2.2：对于每一个观测像元x_j,j＝1,...,m执行如下操作：

2.2a)初始化迭代次数t＝0和丰度

2.2b)计算对应投影点的BMM非线性混合向量

对于GBM和FM，定义为

对于PPNM，定义为：

2.2c)求解最优化问题(10)得系数并以式(11)计算新投影位置

2.2d)用线性丰度估计算法求解的丰度s^(t+1)；

2.2e)如达到最大迭代次数或满足收敛精度则输出丰度矩阵S＝S^(t+1)，否则回到2.2b)；

步骤2.3：输出结果S＝(s₁,s₂,…,s_m)，并用线性规划求得非线性参数γ_i,k,j与b_j。

2.根据权利要求1所述的解混方法，其特征在于：步骤2.2c)中的系数为最小二乘解：

3.根据权利要求1所述的解混方法，其特征在于：步骤2.2d)中采用FCLS算法求解s^(t+1)。