Shift-Invariant Canonical Polyadic Decomposition of Complex-Valued Multi-Subject fMRI Data with a Phase Sparsity Constraint

Li Dan Kuang, Qiu Hua Lin, Xiao Feng Gong, Fengyu Cong, Yu Ping Wang, Vince D. Calhoun

Research output: Contribution to journalArticlepeer-review


Canonical polyadic decomposition (CPD) of multi-subject complex-valued fMRI data can be used to provide spatially and temporally shared components among groups with both magnitude and phase information. However, the CPD model is not well formulated due to the large subject variability in the spatial and temporal modalities, as well as the high noise level in complex-valued fMRI data. Considering that the shift-invariant CPD can model temporal variability across subjects, we propose to further impose a phase sparsity constraint on the shared spatial maps to denoise the complex-valued components and to model the inter-subject spatial variability as well. More precisely, subject-specific time delays are first estimated for the complex-valued shared time courses in the framework of real-valued shift-invariant CPD. Source phase sparsity is then imposed on the complex-valued shared spatial maps. A smoothed $\ell _{\mathbf {{0}}}$ norm is specifically used to reduce voxels with large phase values after phase de-ambiguity based on the small phase characteristic of BOLD-related voxels. The results from both the simulated and experimental fMRI data demonstrate improvements of the proposed method over three complex-valued algorithms, namely, tensor-based spatial ICA, shift-invariant CPD and CPD without spatiotemporal constraints. When comparing with a real-valued algorithm combining shift-invariant CPD and ICA, the proposed method detects 178.7% more contiguous task-related activations.

Original languageEnglish (US)
Article number8805171
Pages (from-to)844-853
Number of pages10
JournalIEEE transactions on medical imaging
Issue number4
StatePublished - Apr 2020


  • Canonical polyadic decomposition (CPD)
  • complex-valued fMRI data
  • shift-invariant
  • source phase sparsity
  • spatiotemporal constraints

ASJC Scopus subject areas

  • Software
  • Radiological and Ultrasound Technology
  • Computer Science Applications
  • Electrical and Electronic Engineering


Dive into the research topics of 'Shift-Invariant Canonical Polyadic Decomposition of Complex-Valued Multi-Subject fMRI Data with a Phase Sparsity Constraint'. Together they form a unique fingerprint.

Cite this