Application of the sequential matrix diagonalization algorithm to high-dimensional functional MRI data. (English) Zbl 1482.62006

Summary: This paper introduces an adaptation of the sequential matrix diagonalization (SMD) method to high-dimensional functional magnetic resonance imaging (fMRI) data. SMD is currently the most efficient statistical method to perform polynomial eigenvalue decomposition. Unfortunately, with current implementations based on dense polynomial matrices, the algorithmic complexity of SMD is intractable and it cannot be applied as such to high-dimensional problems. However, for certain applications, these polynomial matrices are mostly filled with null values and we have consequently developed an efficient implementation of SMD based on a GPU-parallel representation of sparse polynomial matrices. We then apply our “sparse SMD” to fMRI data, i.e. the temporal evolution of a large grid of voxels (as many as 29,328 voxels). Because of the energy compaction property of SMD, practically all the information is concentrated by SMD on the first polynomial principal component. Brain regions that are activated over time are thus reconstructed with great fidelity through analysis-synthesis based on the first principal component only, itself in the form of a sequence of polynomial parameters.


62-08 Computational methods for problems pertaining to statistics
62R10 Functional data analysis
65F15 Numerical computation of eigenvalues and eigenvectors of matrices
94A08 Image processing (compression, reconstruction, etc.) in information and communication theory


Full Text: DOI


[1] Alrmah MA, Weiss S, Redif S, Lambotharan S, McWhirter JG (2013) Angle of arrival estimation for broadband signals: a comparison. In: Proceedings on intelligent signal processing conference, IET, London, UK
[2] Bai, Z.; Demmel, J.; Dongarra, J.; Ruhe, A.; Van der Vorst, H., Templates for the solution of algebraic eigenvalue problems: a practical guide (2000), Philadelphia: SIAM, Philadelphia · Zbl 0965.65058
[3] Baumgartner, R.; Ryner, L.; Richter, W.; Summers, R.; Jarmasz, M.; Somorjai, R., Comparison of two exploratory data analysis methods for fMRI: fuzzy clustering vs. principal component analysis, Magn Reson Imaging, 18, 89-94 (2000)
[4] Beckmann, C.; Smith, S., Tensorial extensions of independent component analysis for multisubject fMRI analysis, NeuroImage, 25, 1, 294-311 (2005)
[5] Bell, AJ; Sejnowski, TJ, An information maximisation approach to blind separation and blind deconvolution, Neural Comput, 7, 6, 1129-1159 (1995)
[6] Carcenac, M.; Redif, S., A highly scalable modular bottleneck neural network for image dimensionality reduction and image transformation, Appl Intell, 44, 3, 557-610 (2016)
[7] Carcenac, M.; Redif, S.; Kasap, S., GPU parallelization of the sequential matrix diagonalization algorithm and its application to high-dimensional data, J Supercomput, 73, 8, 3603-3634 (2017)
[8] Comon, P., Independent component analysis—a new concept?, Signal Process, 36, 3, 287-314 (1994) · Zbl 0791.62004
[9] Corr J, Thomson K, Weiss S, McWhirter JG, Redif S, Proudler IK (2014) Multiple shift maximum element sequential matrix diagonalisation for parahermitian matrices. In: IEEE workshop on statistical signal processing, Gold Coast, Australia, pp 312-315
[10] Fisher M (2014) Marching cubes. https://graphics.stanford.edu/ mdfisher/MarchingCubes.html. Accessed 20 July 2018
[11] Friston, K.; Josephs, O.; Rees, G.; Tuner, R., Nonlinear event-related responses in fMRI, Magn Reson Med, 39, 41-52 (1998)
[12] Glover, G., Overview of functional magnetic resonance imaging, Neurosurg Clin N Am, 22, 2, 133-139 (2011)
[13] Golay, X.; Kollias, S.; Stoll, G.; Meier, D.; Valavanis, A.; Boesiger, PA, New correlation-based fuzzy logic clustering algorithm for fMRI, Magn Reson Imaging, 40, 249-260 (1998)
[14] Golub, GH; Van Loan, CF, Matrix computations (2013), Baltimore, MD: The Johns Hopkins Univ. Press, Baltimore, MD
[15] Hanke M, Dinga R, Hausler C, Guntupalli JS, Casey M, Kaule FR, Stadler J (2015) High-resolution 7-Tesla fMRI data on the perception of musical genres—an extension to the studyforrest dataset. https://f1000research.com/articles/4-174/v1. Accessed 20 July 2018
[16] Intel Corporation (2018) Developer reference for Intel Math Kernel Library 2018 - C. https://software.intel.com/en-us/mkl-reference-manual-for-c. Accessed 20 July 2018
[17] Kailath, T., Linear systems (1980), Englewood Cliffs, NJ: Prentice-Hall, Englewood Cliffs, NJ · Zbl 0458.93025
[18] Kasap, S.; Redif, S., Novel field-programmable gate array architecture for computing the eigenvalue decomposition of para-Hermitian polynomial matrices, IEEE Trans VLSI Syst, 22, 3, 522-536 (2014)
[19] Kim, SG; Ugurbil, K., Functional magnetic resonance imaging of the human brain, J Neurosci Methods, 74, 229-243 (1997)
[20] Krishnan, A.; Williams, LJ; McIntosh, AR; Abdi, H., Partial least squares methods for neuroimaging: a tutorial and review, NeuroImage, 56, 455-475 (2011)
[21] Lazar, NA; Luna, B.; Sweeney, JA; Eddy, WF, Combining brains: a survey of methods for statistical pooling of information, NeuroImage, 16, 538-550 (2002)
[22] Mandelkow, H.; De Zwart, JA; Duyn, JH, Linear discriminant analysis achieves high classification accuracy for the BOLD fMRI response to naturalistic movie stimuli, Front Hum Neurosci, 10 (2016)
[23] McLachlan, GJ, Discriminant analysis and statistical pattern recognition (1992), Hoboken: Wiley, Hoboken
[24] McWhirter, JG; Baxter, PD; Cooper, T.; Redif, S.; Foster, J., An EVD algorithm for para-Hermitian polynomial matrices, IEEE Trans Signal Process, 55, 5, 2158-2169 (2007) · Zbl 1388.94021
[25] Miller, K.; Luh, W.; Lie, T.; Martinez, A.; Obata, T.; Wong, E.; Frank, L.; Buxton, R., Nonlinear temporal dynamics the cerebral blood flow response, Hum Brain Mapp, 13, 1-12 (2001)
[26] Moret N, Tonello A, Weiss S (2011) MIMO precoding for filter bank modulation systems based on PSVD. In: Proceedings on IEEE 73rd vehicle technology conference, pp 1-95
[27] Nvidia Corporation (2018) CUDA Toolkit Documentation v9.1.85. http://docs.nvidia.com/cuda. Accessed 20 July 2018
[28] Poldrack Lab and Center for Reproducible Neuroscience at Stanford University (2018) OpenfMRI. https://openfmri.org. Accessed 20 July 2018 (recently superseded by https://openneuro.org)
[29] Polizzi, E., Density-matrix-based algorithm for solving eigenvalue problems, Phys Rev B, 79, 115112 (2009)
[30] Redif, S., Fetal electrocardiogram estimation using polynomial eigenvalue decomposition, Turk J Electr Eng Comput Sci (2015)
[31] Redif, S.; Kasap, S., Novel reconfigurable hardware architecture for polynomial matrix multiplications, IEEE Trans VLSI Syst, 23, 3, 454-465 (2015)
[32] Redif S, McWhirter JG, Baxter P, Cooper T (2006) Robust broadband adaptive beamforming via polynomial eigenvalues. In: Proceeding on IEEE OCEAN conference, pp 1-6
[33] Redif S, Weiss S, McWhirter JG (2011) An approximate polynomial matrix eigenvalue decomposition algorithm for para-Hermitian matrices. In: Proceedings on 11th IEEE international symposium on signal processing and information technology, Bilbao, Spain, pp 421-425
[34] Redif, S.; Weiss, S.; McWhirter, JG, Sequential matrix diagonalisation algorithms for polynomial EVD of parahermitian matrices, IEEE Trans Signal Process, 63, 1, 81-89 (2015) · Zbl 1394.94479
[35] Redif, S.; Weiss, S.; McWhirter, JG, Relevance of polynomial matrix decompositions to broadband blind signal separation, Signal Process, 134, 76-86 (2017)
[36] Ros, BP; Bijma, F.; De Gunst, MC; De Munck, JC, A three domain covariance framework for EEG/MEG data, NeuroImage, 119, 305-315 (2015)
[37] Shen, H.; Wang, LB; Liu, YD; Hu, DW, Discriminative analysis of resting-state functional connectivity patterns of schizophrenia using low dimensional embedding of fMRI, NeuroImage, 49, 3110-3121 (2014)
[38] Tkacenko A (2010) Approximate eigenvalue decomposition of para-Hermitian systems through successive FIR paraunitary transformations. In: Proceedings on IEEE international conference on acoustics, speech and signal processing, Dallas, TX, USA, pp 4074-4077
[39] Tohidian, M.; Amindavar, H.; Reza, AM, A DFT-based approximate eigenvalue and singular value decomposition of polynomial matrices, EURASIP J Adv Signal Process, 1, 1-16 (2013)
[40] Turner, B.; Paul, E.; Miller, M.; Barbey, A., Small sample sizes reduce the replicability of task-based fMRI studies, Commun Biol (2018)
[41] Vaidyanathan, PP, Multirate systems and filter banks (1993), Englewood Cliffs, NJ: Prentice-Hall, Englewood Cliffs, NJ
[42] Virta J, Li B, Nordhausen K, Oja H (2016) Independent component analysis for tensor-valued data. Preprint available as arXiv:1602.00879 · Zbl 1381.62107
[43] Weiss, S.; Pestana, J.; Proudler, IK, On the existence and uniqueness of the eigenvalue decomposition of a parahermitian matrix, Trans Signal Process, 66, 10, 2659-2672 (2018) · Zbl 1414.15019
[44] Wellcome Trust Centre for Neuroimaging (2018) Statistical parametric mapping (SPM). http://www.fil.ion.ucl.ac.uk/spm. Accessed 20 July 2018
[45] Zarahn, E.; Aquirre, GK; D’Esposito, M., Empirical analysis of BOLD fMRI statistics, NeuroImage, 5, 179-197 (1997)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.