一种适于多被试fMRI数据分析的快速移不变CPD方法技术

一种适于多被试fMRI数据分析的快速移不变CPD方法，属于医学信号处理领域。在移不变CPD算法的基础上，利用交替最小二乘法对被试共享SM成分、共享TC成分和各被试强度更新估计；在不影响时延估计性能的情况下，将维数高的共享SM成分和原始多被试fMRI数据通过矩阵相乘和转换成维数低的数据，从而加快算法运算速度以及降低运算内存。本发明专利技术所需内存明显降低，而且能够快速且有效地对任务态多被试fMRI数据的任务相关成分进行估计。

A Fast Shift-invariant CPD Method for Multi-subject fMRI Data Analysis

A fast shift invariant CPD method suitable for multi-subject fMRI data analysis belongs to the field of medical signal processing. On the basis of shift invariant CPD algorithm, alternant least squares method is used to update the estimation of shared SM components, shared TC components and strength of all subjects. Without affecting the performance of time delay estimation, high-dimensional shared SM components and original multi-subjects fMRI data are multiplied and converted into low-dimensional data by matrix, thus speeding up the algorithm and reducing the computational memory. The memory required by the invention is significantly reduced, and the task-related components of task-state multi-subject fMRI data can be quickly and effectively estimated.

【技术实现步骤摘要】
一种适于多被试fMRI数据分析的快速移不变CPD方法
本专利技术涉及医学信号处理领域，特别是涉及一种多被试功能磁共振成像fMRI(functionalmagneticresonanceimaging)数据的分析方法。
技术介绍
采用磁共振扫描仪对多个被试的大脑进行扫描，获取的脑功能数据称为多被试fMRI数据。fMRI凭借其无损伤且空间分辨率高等优势，成为了目前脑科学研究的一大重要技术。多被试fMRI数据一般看成为一个三维张量，包括空间维、时间维和被试维，可以采用张量分解算法来处理。CPD(canonicalpolyadicdecomposition)是一种典型的张量分解算法。CPD在多被试fMRI数据分析中物理意义明确，将多被试fMRI数据分解为各被试间共享脑空间激活成分(spatialmap,SM)和共享时间过程成分(timecourse,TC)，以及各被试的强度差异信息。这些信息可以为脑功能研究和脑疾病诊断提供重要依据。多被试fMRI数据不可避免地存在着被试间空时差异性。其中时间差异性问题主要是由于各被试响应时间和血液动力学时延方面存在差异导致的。针对该问题，等人在2008年文章“M.,Hansen,L.K.,Arnfred,S.M.,Lim,L.H.,Madsen,K.H.,Shift-invariantmultilineardecompositionofneuroimagingdata.NeuroImage42,2018,1439-1450.”中提出了一种移不变CPD算法，在无约束CPD模型中加入了被试TC时延估计，容许被试TC时延差异，分离性能优于无约束CPD算法。近几年，大规模的多被试fMRI数据(被试数目从几十到上万)的研究越来越受到人们的关注，而且fMRI数据的空间维数通常很高，比如5万多个脑内体素，因此当移不变CPD算法应用于多被试fMRI数据时，存在着运行速度慢且内存需求高的问题，特别是在被试时延的更新估计过程中，需要存储的数据量比较大。
技术实现思路
本专利技术提供了一种适于多被试fMRI数据的快速移不变CPD方法，通过矩阵相乘和变换，让被试时延估计速度变快、运行内存降低，且性能不会下降。在移不变CPD算法的基础上，利用交替最小二乘法对被试共享SM成分、共享TC成分和各被试强度更新估计；在不影响时延估计性能的情况下，将维数高的共享SM成分和原始多被试fMRI数据通过矩阵相乘和转换成维数低的数据，从而加快算法运算速度以及降低运算内存。本专利技术的技术方案步骤如下：第一步：输入多被试fMRI数据其中I表示脑内体素数目(空间维)，J表示扫描次数(时间维)，K表示被试数目。第二步：初始化。设成分个数为D(D为正整数且0<D≤J)。随机初始化共享SM成分共享TC成分和被试强度初始化被试时延为零矩阵。令迭代次数iter＝0，相对误差Δεiter＝1，计算迭代误差εiter：其中，bd(j-τk,d)表示对bj,d时移了τk,d(为整数)个点，τk,d为第k被试第d个成分的时延。若τk,d>0，第k被试第d个TC成分相对共享TC成分bd循环左移τk,d个点，若τk,d<0，则循环右移|τk,d|个点。式(1)也是移不变CPD算法模型。第三步：更新共享TC成分B。采用文章“M.,Hansen,L.K.,Arnfred,S.M.,Lim,L.H.,Madsen,K.H.Shift-invariantmultilineardecompositionofneuroimagingdata.NeuroImage42,2018,1439-1450.”中共享TC成分B的更新方法(交替最小二乘法)对B进行更新。第四步：对共享SM成分A和X进行降维。对A进行如下运算：其中由于D(一般小于100)通常远远小于I(一般5×105)，因此用维数低的替代维数高的A(D×I)，从而加快被试时延的估计，并且降低计算内存。进一步令其中为X的1模展开形式，将张量化为其维数也远小于X。第五步：更新被试时延将按3模展开，得到矩阵采用文章“M.,Hansen,L.K.,Arnfred,S.M.,Lim,L.H.,Madsen,K.H.Shift-invariantmultilineardecompositionofneuroimagingdata.NeuroImage42,1439-1450.”中被试时延估计方法，并用第四步中的和替换文章中的A和X(3)(为X的3模展开形式)，实现被试时延的快速估计。第六步：更新共享SM成分A。采用文章“M.,Hansen,L.K.,Arnfred,S.M.,Lim,L.H.,Madsen,K.H.Shift-invariantmultilineardecompositionofneuroimagingdata.NeuroImage42,2018,1439-1450.”中更新A的方法对A更新。第七步：更新被试强度C。采用文章“M.,Hansen,L.K.,Arnfred,S.M.,Lim,L.H.,Madsen,K.H.Shift-invariantmultilineardecompositionofneuroimagingdata.NeuroImage42,2018,1439-1450.”中更新C的方法对C更新。第八步：计算误差。令iter＝iter+1；根据式(1)，计算本次迭代误差εiter，以及相对误差Δεiter：Δεiter＝|(εiter-1-εiter)/εiter-1|。(4)第九步：若εiter小于预设误差阈值εiter_min，跳转到第十二步，否则执行第十步。第十步：若Δεiter小于预设相对误差阈值Δεiter_min，跳转到第十二步，否则执行第十一步。第十一步：若iter大于预设最大迭代次数itermax，跳转到第十二步，否则执行第三步。第十二步：输出共享SM成分A、共享TC成分B、被试时延与被试强度C。本专利技术所达到的效果和益处是，能够快速且有效地对任务态多被试fMRI数据的任务相关成分进行估计。在对10被试敲击手指任务的fMRI数据分析中，设置相同初始值的情况下，本专利技术所需的收敛次数跟移不变CPD算法相同，但每次迭代所需运行速度快约9.33倍，而且被试时延估计所需内存约是移不变CPD算法的2/25。本专利技术和移不变CPD算法估计的共享任务相关成分与先验参考信号的相关系数均值之比约为1，分离性能相似。因此，本专利技术能快速且有效地提取出多被试共享的脑功能信息，未来可与独立成分分析和稀疏表示等方法相结合，提高分离性能，在今后大规模fMRI数据研究和智慧医疗等方面也具有良好的发展前景。附图说明图1是本专利技术分析多被试fMRI数据的工作流程图。具体实施方式下面结合技术方案和附图，详细叙述本专利技术的一个具体实施例。现有10被试执行敲击手指任务下采集的fMRI数据，即K＝10。每个被试都进行了J＝165次扫描，每次扫描采集了53×63×46的全脑数据，去掉脑外数据体素，保留脑内数据体素I＝59610。假设共享SM和TC成分的成分个数D＝30，采用本专利技术进行多被试fMRI数据分析的步骤如附图所示。第一步：输入多被试fMRI数据第二步：初始化。随机初始化共享TC成分共享SM成分被试强本文档来自技高网...

【技术保护点】
1.一种适于多被试fMRI数据的快速移不变CPD方法，其特征包括以下步骤：第一步：输入多被试fMRI数据

【技术特征摘要】
1.一种适于多被试fMRI数据的快速移不变CPD方法，其特征包括以下步骤：第一步：输入多被试fMRI数据其中I表示脑内体素数目，J表示扫描次数，K表示被试数目；第二步：初始化：设成分个数为D，D为正整数且0<D≤J；随机初始化共享SM成分共享TC成分和被试强度初始化被试时延为零矩阵；令迭代次数iter＝0，相对误差Δεiter＝1，计算迭代误差εiter：其中，bd(j-τk,d)表示对bj,d时移了τk,d(为整数)个点，τk,d为第k被试第d个成分的时延。若τk,d>0，第k被试第d个TC成分相对共享TC成分bd循环左移τk,d个点，若τk,d<0，则循环右移|τk,d|个点；式(1)也是移不变CPD方法模型；第三步：更新共享TC成分B；第四步：对共享SM成分A和X进行降维；对A进行如下运算：其中由于D小于I，用维数低的替代维数高的A(D×I)，从而加快被试时延的估...

【专利技术属性】
技术研发人员：邝利丹，林秋华，龚晓峰，丛丰裕，
申请(专利权)人：长沙理工大学，
类型：发明
国别省市：湖南,43