基于耦合张量分解的高维多被试复数fMRI数据分析

功能磁共振成像（fMRI）是脑科学研究的重要手段之一。多被试fMRI本质上是高维、复数数据，但分析时一般将其低维化且只分析幅值数据，导致分析精度损失及脑信息利用不全。充分利用复数fMRI数据的高维结构，有助于提取更完整、更切实际的脑功能特性。然而，高维多被试复数fMRI数据存在脑外无用数据多（约占全部扫描数据的2/3），脑内各切片体素数目不同、多被试空时差异、复数数据噪声高等问题。为此，本项目拟基于耦合张量分解给出解决方案。耦合张量分解模型灵活，既善于刻画高维数据结构，又能对不同张量进行联合分解，容许不同切片张量存在差异。本项目拟以高维多被试复数fMRI脑内数据为研究对象，建立其耦合张量化模型，提出复数空时约束的耦合张量分解算法，提出同时考虑统计和结构特性的模型阶数估计方法，挖掘复数共享空时成分的脑功能信息，提升精神分裂症的分类准确率。本项目方法易于推广到其他脑疾病，助力我国脑科学发展。

In recent years, studies in brain sciences have attracted worldwide attention. Functional magnetic resonance imaging (fMRI) with the advantages of non-invasive nature and high spatial resolution, has been widely applied in brain science research. In fact, multi-subject fMRI data are intrinsically high-dimensional (i.e., five-ways tensor) and complex-valued data. However, the majority of studies reshape the data into three-ways tensor or matrix, and only focus on magnitude-only data analysis. Thus, this leads to the analysis accuracy loss and deficient brain functional information exploitation. Making full use of the high-dimensional structure of complex-valued fMRI data, can help to extract more complete and practical brain functional characteristics. Nevertheless, when we analyze the high-dimensional multi-subject complex-valued fMRI data, we will face the following three problems: 1) the scan-acquired whole fMRI data contain many useless out-brain voxels (about 2/3 of whole voxels), and the number of in-brain voxels of each slice is different; 2) large inter-subject variability; 3) high-noisy nature of complex-valued fMRI data. This proposal intends to utilize coupled tensor decomposition to give the solutions for the above problems. Coupled tensor decomposition model is very flexible. It can not only depict high-dimensional data structure, but also jointly decompose different tensors, allowing differences among different slice tensors. To this end, this proposal firstly aims to establish coupled high-dimensional tensorization models of high-dimensional multi-subject complex-valued fMRI in-brain data. Secondly, coupled tensor decomposition with spatiotemporal constraints are to be proposed. Thirdly, efficient model order estimation methods are to be proposed by jointly exploiting both the statistical and structural properties of the data. Finally, based on the resting-state fMRI data consisting of schizophrenia patients (SZs) and healthy controls (HCs), the proposal intends to explore the brain functional information of shared complex-valued spatiotemporal components, and improve the classification accuracy of SZs and HCs. Our proposed methods can be straightforwardly extended to other brain disease studies, and contribute to our Chinese brain science research.