Fast and Structured Block-Term Tensor Decomposition for Hyperspectral Unmixing

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.

Keywords

• Hyperspectral unmixing
• structured block-term tensor decomposition
• alternating gradient projection