Functional magnetic resonance imaging (fMRI) is one of the most popular methods in the study of the human brain. It is a technique that uses magnetic resonance imaging (MRI) to access brain activity by detecting local changes in blood oxygen levels.


FMRI is used both for brain function research and as a diagnostic tool. Through fMRI, it is possible to study the neural networks formed in the human brain in order to perform brain functions at the level of topology and structure. Conclusions on brain function, such as the structure and morphology of neural networks, are achieved by examining multiple healthy subjects that perform common activities. The use of fMRI as a diagnostic tool has been suggested for many diseases through a large number of tasks. Any deviations and alterations in the characteristics of neural networks can be detected in a timely manner, even at a pre-symptomatic stage, resulting in early detection of neurodegenerative diseases (Alzeimer's disease, Parkinson's disease, etc.) or mental disorders (schizophrenia, bipolar disorder, etc.).


Analysis of fMRI data is particularly complex, mainly due to two factors. The first factor is associated with the need for high spatial and temporal clarity, which leads to the creation of large volume data. The second factor concerns the nature of the data. The fMRI signal is not directly related to neuronal activity but to the changes in oxygen levels observed locally and associated with the metabolic requirements of neurons (associated and unrelated to a particular brain function). In addition, the fMRI signal is weak and overlaps by many noise sources associated with the signal acquisition process. Therefore, the processing of fMRI is composed by the removal of the noise (pre-processing), as well as the components unrelated to the intended brain function, as well as the extraction of information about the neural function.


The fMRI signal processing methods, after the pre-processing step, fall into two categories. Univariate analysis methods, based on the analysis of the signal generated by each brain area separately (voxel-wise signal analysis), appear in the first category. In the second category, there are multivariate analysis methods that analyze the signals that arise from each area of the brain and aim to draw conclusions on how a brain function is performed in a human brain (finding neural networks associated with a brain function and temporal behavior of these).


A new trend in the multivariate analysis methods is the use of tensor models (PARAFAC and TUCKER). In contrast to the most popular methods for analyzing multiple variables (PCA, ICA, etc.), the use of tensor models allows data analysis to maintain its multidimensional structure and offers unique data decompositions (PARAFAC model) to components, each expressed by a set of one-dimensional signatures, making it easier to interpret the results. When the data are not characterized by the PARAFAC model's linear structure, degeneration phenomena occur during the decomposition process, resulting in the extraction of components that are highly concordant and are useless to interpret information.


The fMRI data do not follow the multilinear structure of the PARAFAC model, as the signals corresponding to adjacent regions exhibit delayed correlation phenomena due to the hemodynamic phenomena that develop in the interior of the brain. In order to enable the processing of fMRI data, variants of the PARAFAC model have been proposed in the literature, where delays between signals originating from adjacent regions of the brain have been predicted. In particular, models have been porposed that incorporate the unique delays for each voxel per component (Shift Invariant PARAFAC), as well as the existence of a parameterizable array of delays for each voxel per component (Convolutive PARAFAC). However, the decomposition process for both models is particularly time consuming and no parallel disintegration pattern has been proposed, as is a method of estimating the appropriate number of components.