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.

Alternating direction method of multipliers for linear hyperspectral unmixing

Blind Spectral Unmixing Considering the Adjacent Effect

Mixing Support Detection-Based Alternating Direction Method of Multipliers for Sparse Hyperspectral Image Unmixing

Hyperspectral Unmixing with Scaled and Perturbed Linear Mixing Model to Address Spectral Variability

A Dual Symmetric Gauss-Seidel Alternating Direction Method of Multipliers for Hyperspectral Sparse Unmixing

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

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

Correntropy Maximization Via ADMM: Application to Robust Hyperspectral Unmixing

Dual Alternating Direction Method of Multipliers for Inverse Imaging

Parallel Adaptive Sparsity-Constrained NMF Algorithm for Hyperspectral Unmixing.

Unsupervised Nonlinear Hyperspectral Unmixing Based on an Extended Multilinear Mixing Model-Inspired Dual-Stream Network

Alternating Direction Method of Multipliers for Linear Inverse Problems

MultiHU-TD: Multifeature Hyperspectral Unmixing Based on Tensor Decomposition

Alternating Direction Iterative Nonnegative Matrix Factorization Unmixing for Multispectral and Hyperspectral Data Fusion

A Plug-and-Play Priors Framework for Hyperspectral Unmixing

Bilinear Hyperspectral Unmixing Via Tensor Decomposition

Half-quadratic Based Robust Hyperspectral Unmixing Framework

Hypergraph-Regularized Sparse NMF for Hyperspectral Unmixing.

Hyperspectral Unmixing Using Robust Deep Nonnegative Matrix Factorization

Relaxing Alternating Direction Method of Multipliers (ADMM) for Linear Inverse Problems

Hyperspectral Unmixing Via Projected Mini-Batch Gradient Descent