Song Mei-ping, Xu Xing-wei, Chang Chein-I, An Ju-bai, Yao Li
Guang Pu Xue Yu Guang Pu Fen Xi. 2015 Dec;35(12):3465-70.
Spectrum unmixing is an important part of hyperspectral technologies, which is essential for material quantity analysis in hyperspectral imagery. Most linear unmixing algorithms require computations of matrix multiplication and matrix inversion or matrix determination. These are difficult for programming, especially hard for realization on hardware. At the same time, the computation costs of the algorithms increase significantly as the number of endmembers grows. Here, based on the traditional algorithm Orthogonal Subspace Projection, a new method called. Orthogonal Vector Projection is prompted using orthogonal principle. It simplifies this process by avoiding matrix multiplication and inversion. It firstly computes the final orthogonal vector via Gram-Schmidt process for each endmember spectrum. And then, these orthogonal vectors are used as projection vector for the pixel signature. The unconstrained abundance can be obtained directly by projecting the signature to the projection vectors, and computing the ratio of projected vector length and orthogonal vector length. Compared to the Orthogonal Subspace Projection and Least Squares Error algorithms, this method does not need matrix inversion, which is much computation costing and hard to implement on hardware. It just completes the orthogonalization process by repeated vector operations, easy for application on both parallel computation and hardware. The reasonability of the algorithm is proved by its relationship with Orthogonal Sub-space Projection and Least Squares Error algorithms. And its computational complexity is also compared with the other two algorithms', which is the lowest one. At last, the experimental results on synthetic image and real image are also provided, giving another evidence for effectiveness of the method.
光谱分解是高光谱技术的重要组成部分,对于高光谱图像中的物质定量分析至关重要。大多数线性分解算法需要进行矩阵乘法、求逆或行列式运算。这些运算编程困难,尤其难以在硬件上实现。同时,随着端元数量的增加,算法的计算成本会显著提高。在此,基于传统的正交子空间投影算法,利用正交原理提出了一种名为正交向量投影的新方法。该方法通过避免矩阵乘法和求逆简化了这一过程。它首先通过Gram-Schmidt过程为每个端元光谱计算最终的正交向量。然后,将这些正交向量用作像素特征的投影向量。通过将特征投影到投影向量上并计算投影向量长度与正交向量长度之比,可以直接获得无约束丰度。与正交子空间投影算法和最小二乘误差算法相比,该方法无需矩阵求逆,而矩阵求逆计算成本高且难以在硬件上实现。它只需通过重复向量运算完成正交化过程,便于在并行计算和硬件上应用。通过与正交子空间投影算法和最小二乘误差算法的关系证明了该算法的合理性。并且还将其计算复杂度与其他两种算法进行了比较,结果表明其计算复杂度最低。最后,给出了合成图像和真实图像上的实验结果,进一步证明了该方法的有效性。