Abstract:The block-term tensor decomposition model with multilinear rank- $(L_{r},L_{r},1)$ terms (or the “ ${\mathsf{LL1}}$ tensor decomposition” in short) offers a valuable alternative formulation for hyperspectral unmixing (HU), which ensures the identifiability of the endmembers/abundances in cases where classic matrix factorization (MF) approaches cannot provide such guarantees. However, the existing ${\mathsf{LL1}}$ -tensor-decomposition-based HU algorithms use a three-factor parameterization of the tensor (i.e., the hyperspectral image cube), which causes difficulties in incorporating structural prior information arising in HU. Consequently, their algorithms often exhibit high per-iteration complexity and slow convergence. This article focuses on ${\mathsf{LL1}}$ tensor decomposition under structural constraints and regularization terms in HU. Our algorithm uses a two-factor reparameterization of the tensor model. Like in the MF-based approaches, the factors correspond to the endmembers and abundances in the context of HU. Thus, the proposed framework is natural to incorporate physics-motivated priors in HU. To tackle the formulated optimization problem, a two-block alternating gradient projection (GP)-based algorithm is proposed. Carefully designed projection solvers are proposed to implement the GP algorithm with a relatively low per-iteration complexity. An extrapolation-based acceleration strategy is proposed to expedite the GP algorithm. Such an extrapolated multiblock algorithm only had asymptotic convergence assurances in the literature. Our analysis shows that the algorithm converges to the vicinity of a stationary point within finite iterations, under reasonable conditions. Empirical study shows that the proposed algorithm often attains orders-of-magnitude speedup and substantial HU performance gains compared with the existing ${\mathsf{LL1}}$ -decomposition-based HU algorithms.

Constrained Block-Term Tensor Decomposition-Based Hyperspectral Unmixing Via Alternating Gradient Projection

Fast and Structured Block-Term Tensor Decomposition for Hyperspectral Unmixing.

Bilinear Hyperspectral Unmixing Via Tensor Decomposition

Blind Spectral Unmixing Considering the Adjacent Effect

Parallel Adaptive Sparsity-Constrained NMF Algorithm for Hyperspectral Unmixing.

Hyperspectral Unmixing Via Projected Mini-Batch Gradient Descent

Unmixing of Large-Scale Hyperspectral Data Based on Projected Mini-Batch Gradient Descent.

A Fast Fully Constrained Geometric Unmixing of Hyperspectral Images

Hypergraph-Regularized Sparse NMF for Hyperspectral Unmixing.

MultiHU-TD: Multifeature Hyperspectral Unmixing Based on Tensor Decomposition

Hyperspectral Super-Resolution: A Coupled Nonnegative Block-Term Tensor Decomposition Approach

Hyperspectral Unmixing via Latent Multiheterogeneous Subspace

Hyperspectral Unmixing Via Total Variation Regularized Nonnegative Tensor Factorization.

Nonnegative Matrix Factorization with Data-Guided Constraints for Hyperspectral Unmixing.

Factor-Regularized Nonnegative Tensor Decomposition for Blind Hyperspectral Unmixing.

Minimum Variance Block-Based Nonnegative Matrix Factorization Algorithm for Hyperspectral Unmixing.

Half-quadratic Based Robust Hyperspectral Unmixing Framework

Alternating direction method of multipliers for linear hyperspectral unmixing

Superpixel-Based Nonnegative Tensor Factorization for Hyperspectral Unmixing

Blind Hyperspectral Unmixing Considering the Adjacency Effect

Adaptive $textit{L}_{bf 1/2}$ Sparsity-Constrained NMF With Half-Thresholding Algorithm for Hyperspectral Unmixing