端元初始化

端元估计在高光谱图像处理中的重要性

高光谱图像处理领域经常会遇到端元矩阵初始化的问题,而端元矩阵的准确估计对于最终处理结果的性能有着重要的影响。在处理凸优化问题时,初始值的选择对最终结果的影响相对较小。然而,面对非凸优化问题,一个接近全局最优解的初始值能够显著提升算法的性能。

凸优化与非凸优化

端元数量的重要性

在高光谱领域的多种模型中,端元数量扮演着关键角色。例如,给定一个高光谱图像 \(\mathcal{Z} \in \mathbb{R}^{M\times N\times B}\),我们假设其端元数量为 \(R\)。在非负矩阵分解模型中,我们有:

\[ \mathbf{Z} \in \mathbb{R}^{B\times MN} = (\mathbf{A} \in \mathbb{R}^{B\times R})(\mathbf{S} \in \mathbb{R}^{R\times MN}) \]

式中,\(\mathbb{R}^{B\times MN}\),表示端元矩阵。

在张量CP分解模型中,表达式变为:

\[ \mathcal{Z} = \sum_{r=1}^{R} \omega_{r}(\mathbf{a}_{r} \circ \mathbf{b}_{r} \circ \mathbf{c}_{r}) \]

式中\(\omega_{r}\)表示第\(r\)个分量的权重,将其中的\(a\)\(b\)进行矩阵化可以得到秩1近似的空间信息矩阵,而\(c\)矩阵就是包含光谱信息的矩阵

对于Tucker分解模型,表达式为:

\[ \mathcal{Z} \in \mathbb{R}^{M\times N\times B} = (\mathcal{G} \in \mathbb{R}^{m\times n\times R}) \times_{1} \mathbf{A} \times_{2} \mathbf{B} \times_{3} (\mathbf{C} \in \mathbb{R}^{B\times R}) \] 式中,\(\mathbf{C} \in \mathbb{R}^{B\times R}\)表示光谱因子矩阵,包含了光谱的成分信息。将它初始化为端元矩阵往往会提升算法的性能以及收敛速度。

最后,在BTD(Block Term Decomposition)分解中,我们得到:

\[ \mathcal{Z} = \sum_{r=1}^{R} \mathcal{G}_{r} \times_{1} \mathbf{A}_{r} \times_{2} \mathbf{B}_{r} \times_{3} \mathbf{C}_{r} \approx \sum_{r=1}^{R} \mathbf{A}_{r} \cdot \mathbf{B}_{r}^{T} \circ \mathbf{c}_{r} = \sum_{r=1}^{R} \mathbf{E}_{r} \circ \mathbf{c}_{r} \]

在BTD分解中进行地址假设对标准分解形式进行华建,最终得到空间因子矩阵和光谱因子矩阵的乘积,矩阵\(\mathbf{C}\)表示高光谱图像中的光谱特征矩阵。同样将它初始化为端元矩阵,可以提升算法的性能和收敛速度。

在这些模型中,端元数量 \(R\) 的选择对于算法的性能和结果的准确性至关重要。正确估计端元数量不仅可以提高处理效率,还能确保最终结果的准确性和可靠性。


端元初始化
http://jingmengzhiyue.top/2024/03/19/endmember/
作者
Jingmengzhiyue
发布于
2024年3月19日
许可协议