基于深度学习的自闭症谱系障碍fMRI数据分类研究.pdf
《基于深度学习的自闭症谱系障碍fMRI数据分类研究.pdf》由会员分享,可在线阅读,更多相关《基于深度学习的自闭症谱系障碍fMRI数据分类研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、现代电子技术Modern Electronics Technique2023年11月1日第46卷第21期Nov.2023Vol.46 No.210 引 言自闭症谱系障碍(Autistic Spectrum Disorder,ASD)是一种发育障碍,可导致重大的社交、沟通和行为挑战1。2014年,自闭症总体患病率估计为每1 000名8岁儿童有 16.8 人患病,部分社区 ASD 患病率2接近 3%。近年来,人们一直担心儿童 ASD 患病率的增加34。然而,目前对自闭症的诊断高度依赖于传统的行为症状,这些行为症状通常具有主观性,容易导致忽视早期症状而误诊56。因此,寻找一种客观的生物标志物用于自闭
2、症的早期诊断和及时干预治疗越来越受到精神病学和神经科学领域的关注。静息态功能核磁共振成像(RestingState Functional Magnetic Resonance Image,rsfMRI)是一种无需执行任何 特 定 任 务 即 可 测 量 受 试 者 血 氧 水 平 依 赖(Blood OxygenationLevel Dependent,BOLD)信号的技术,已广泛应用于神经影像学分析78。使用rsfMRI通常需要构建 功 能 连 接 网 络(Functional Connectivity Networks,FCNs)来表示每个主题,每个元素表示大脑 ROI(Region of
3、 Interest)之间的成对关系910。因此,FCNs 倾向于捕捉大脑成对 ROI的 BOLD信号之间的依赖关系,并已被用于识别潜在的神经影像生物标志物以诊断神经系统基于深度学习的自闭症谱系障碍fMRI数据分类研究陈茂洲,刘化东,许博俊,李梦琪(昆明理工大学 信息工程与自动化学院,云南 昆明 650500)摘 要:与健康儿童对照相比,自闭症谱系障碍(ASD)患者的脑结构和功能存在显著异常,因此文中采用脑影像的方法实现ASD的辅助诊断。由于基于传统的脑图谱构建脑网络节点依赖于人为的先验知识和假设,因此将由受试者的数据通过数据驱动的方法即独立成分分析(ICA)进行全脑独立成分(ICs)分析,并将
4、提取的 ICs作为感兴趣区域(ROI)提取时间序列,随后用计算得到的相关性矩阵构建脑网络并通过Node2vec方法将脑网络的节点向量化以生成节点特征,最后通过图卷积网络(GCN)刻画脑网络连接水平状态,并通过多层感知机(MLP)实现对ASD的识别。所用方法ICA+GCN在79名ASD患者和105名典型对照组中取得了86.1%的分类精度。综上,文中所用方法有望用于ASD辅助诊断。关键词:自闭症谱系障碍;独立成分分析;脑影像;脑网络;Node2vec;图卷积网络;多层感知机;辅助诊断中图分类号:TN911.134;TP391 文献标识码:A 文章编号:1004373X(2023)21004807S
5、tudy on ASD fMRI data classification based on deep learningCHEN Maozhou,LIU Huadong,XU Bojun,LI Mengqi(Faculty of Information Engineering and Automation,Kunming University of Science and Technology,Kunming 650500,China)Abstract:In comparison with healthy children,individuals with autistic spectrum d
6、isorder(ASD)exhibit significant abnormalities in brain structures and functions.Therefore,the brain imaging method is used to assist in the diagnosis of ASD.Since the brain network construction based on the traditional brain atlases depends on prior knowledge and assumptions,the participant data are
7、 subjected to wholebrain independent components(ICs)analysis by a datadriven method called independent component analysis(ICA).The extracted ICs are then used as regions of interest(ROI)to extract time series,and the resulting correlation matrix is used to construct a brain network.The method of Nod
8、e2vec is used to vectorize the brain network nodes to generate node features,and the status of brain network connectivity is characterized by a graph convolutional network(GCN).The recognition of ASD is realized by multilayer perceptron(MLP).The method ICA+GCN achieves classification accuracy of 86.
9、1%in a sample of 79 ASD patients and 105 typical controls.Therefore,this method has the potential to be used as an aid in the diagnosis of ASD.Keywords:ASD;ICA;brain imaging;brain network;Node2vec;GCN;MLP;assisted diagnosisDOI:10.16652/j.issn.1004373x.2023.21.010引用格式:陈茂洲,刘化东,许博俊,等.基于深度学习的自闭症谱系障碍fMRI
10、数据分类研究J.现代电子技术,2023,46(21):4854.收稿日期:20230307 修回日期:202303294848第21期疾病1112。它可以帮助了解大脑组织模式和诊断神经系统疾病,例如 ASD1314、阿尔兹海默症及其前期症状(即轻度认知障碍)15、帕金森16。然而,以往的研究往往首先从 FCNs中提取手工制作的网络特征(如节点的度和聚类系数),然后进行ASD诊断的模型预测17,这两个步骤被分开处理并高度依赖专家知识。随着深度学习的发展,特别是图神经网络(Graph Convolution Network,GCN)已被用于识别大脑 FCNs 中潜在的fMRI生物标志物,以用于疾病
11、诊断1819。一般来说,每个大脑网络都可以看作是一个由不规则数据组成的复杂图结构20,不仅包含节点特征,还包含不同节点之间的拓扑信息。GCN提供了一个端到端的深度学习架构,可以自动学习节点特征和节点间的拓扑信息。目前的研究表明,将 GCN用于 fMRI分析有助于通过卷积和池化等操作自动捕获脑网络的高级拓扑信息,从而显著提高脑疾病的诊断性能21。1 材料与方法1.1 数据集本研究使用的数据来自美国自闭症脑成像数据交换项目(Autism Brain Imaging Data Exchange,ABIDE)数据集中纽约大学提供的184名受试者的fMRI数据,其中包括 79 名 ASD 和 105 名
12、典型对照(Typical Controls,TC)。该站点是该数据集中具有最大样本量的站点。所有受试者经筛选,受试者必须没有癫痫发作史,没有超过 5 min 的意识丧失,没有抽动障碍或不自主运动,也没有任何已知的遗传、神经或精神障碍。其中典型对照组中的受试者必须没有自闭症谱系障碍家族史。典型对照组在群体水平上与自闭症谱系障碍受试者的年龄、性别、惯用手性和全面智商相匹配。所有的fMRI数据均为静息态数据,受试者被要求在扫描仪中躺着不动、闭上眼睛、保持清醒。屏幕上没有显示任何刺激物,扫描结束后,所有参与者都确认他们在整个过程中保持清醒。1.2 研究流程本研究首先使用数据驱动的方法即独立成分分析(I
13、ndependent Component Correlation Algorithm,ICA)对fMRI数据进行全脑独立成分(Independent Components,ICs)分析,并将提取的 ICs作为感兴趣区域(ROI),提取ROI的时间序列;随后通过对比不同方法量化序列间的相关性程度得到相关性矩阵,根据阈值选取策略构建脑网络;构建的脑网络需要进行两步操作,一是将其经过对称化、二值化等策略转换为无向脑网络,二是通过Node2vec方法将脑网络的节点向量化以生成节点特征;然后将得到的无向脑网络和节点特征输入 GCN网络以刻画脑网络连接水平状态;最后通过一层多层感知机(Multilayer
14、 Perceptron,MLP)实现对 ASD的识别,如图 1所示。图1 研究流程图1.3 实验设计由本研究的流程,分别设计了如下实验:1)ICA对fMRI数据进行全脑ICs分析,并将提取的ICs作为ROI,提取ROI的时间序列。2)脑网络阈值选取及属性分析实验。通过分析阈值、脑网络密度和小世界属性三者之间的关系确定阈值;然后计算不同脑网络密度下 ASD和 TC的小世界属性值。3)对比多个方法量化 ASDTC 功能子网络间相关程度。4)根据阈值选取策略构建脑网络。5)将构建的带权有向脑网络经过对称化、二值化等策略转换为无向脑网络且输入 GCN网络来刻画脑网络连接水平状态,并通过MLP实现ASD
15、的识别。1.4 方 法1.4.1 独立成分分析独立成分分析(ICA)是一种完全数据驱动的分析方法,适用于更广泛的领域。之前对功能核磁共振成像数据进行的研究都是在时间序列预先知道的情况下,通过相关分析得到大脑的激活区。如果不知道大脑收集数据的时间序列以及激活大脑的因素,那么就无法确定信号源和混合矩阵。这种情况下通常会采用盲源分离方法来解决问题,其中最常用的方法是独立成分分析。通过这种方法可以对功能核磁共振成像数据进行混合分离,并提取与任务相关的成分。线性分解方法是一种替代聚类来定义节点的方法,其中 ICA 方法较为常用。与大多数体素聚类方法不同的是,ICA等方法不考虑空间邻近信息(也就是说,它们
16、不区分体素在空间上相邻还是相隔较远。因此,这类方法通常会造成不连续的分割)。简而言之,ICA 旨在将数据分解成一组成分,这些成分由独立的空间位置图和描述这些成分随时间变化的时间序列图构成。ICA 对陈茂洲,等:基于深度学习的自闭症谱系障碍fMRI数据分类研究49现代电子技术2023年第46卷fMRI数据进行全脑 ICs分析,并将提取的 ICs作为 ROI,提取ROI的时间序列。图2是组独立成分分析流程。图2 组独立成分分析流程1.4.2 功能子网络间相关程度分析计算各 ICs 间的相关性,计算各 ICs 内所有体素的平均时间序列,作为该ICs的时间序列xi()t。1)Graphical Las
17、son CV稀疏逆协方差(Sparse Inverse Covariance)也称为精度逆协方差(Precision Matrix)或精确相关系数矩阵的逆矩阵(Inverse of the Precise Correlation Matrix)。稀疏逆协方差是指对于一个给定的高维数据集,通过分析数据集中各个变量之间的线性关系,得到一个矩阵,可以用来衡量不同变量之间的相关性,同时也可以用于数据降维和特征选择。稀疏逆协方差公式由式(1)计算:()ij=-1ijSij+12iij(1)式中:为稀疏逆协方差矩阵;i和j为第i和第j个变量 的 标 准 差;S为 协 方 差 矩 阵;ij为 克 罗 内 克
18、(Kronecker)函数,即当i=j时,取值为1,否则为0。2)Spearman斯皮尔曼相关系数也被称为等级相关系数或顺序相关系数,它是一种计算两个变量之间关联性的统计量,使用一种基于样本值大小顺序排列的方法,将实际数据转化为等级形式进行计算。具体而言,该系数的计算方式如下:给定两个向量:X=(X1,X2,Xn),Y=(Y1,Y2,Yn),首先对这两个向量进行排序,得到元素的排行x、y,以xi、yi分别表示Xi在X中的排行以及Yi在Y中的排行,再将集合x、y中的元素对应相减得到一个排行差分集合d,则斯皮尔曼相关系数计算公式为:r=1-6i=1NdiN(N2-1)(2)3)Pearson采 用
19、 式(3)的 Pearson 相 关 系 数 计 算 各 ICs 的 相关性:rij=t=1Txi()t-x ixj()t-x jt=1Txi()t-x i2t=1Txj()t-x j2(3)式中:xi()t和xj()t分别为独立成分i和独立成分j在t时刻的血氧水平依赖(BOLD)值,i,j 1,N,N为独立成分数量;T为时长;x i和x j为脑区i和脑区j在T内的平均BOLD 值;脑区i和脑区j之间的皮尔森相关系数值用rij表示,最终得到大小为N N的相关系数矩阵。1.4.3 脑网络rsfMRI能够较好地反映大脑的自发神经元活动22,已有研究表明 ASD患者的静息态数据反映的功能状态与 TC
20、的不同23。因此,对 rsfMRI构建的脑网络进行量化可以用来识别 ASD。构建脑网络主要有如下三个步骤:1)使用 ICA 对 fMRI 数据进行全脑 ICs 分析,将大脑划分成多个功能区域,把每一个 ICs 定义成网络的一个节点。2)计算各ICs间的相关性。3)制定阈值策略将脑区相关性矩阵转换成邻接矩阵。脑功能网络是一种复杂而稀疏的网络,其主要由大量短程边和少量长程边组成。脑功能网络需要制定合适的阈值将脑网络间的相关系数矩阵转换为二值化邻接矩阵。当相关系数大于设定的阈值时,可以认为两个脑区之间存在功能连接边。反之,如果相关系数小于该阈值,则认为两个脑区之间不存在功能连接边。研究表明,脑网络通
21、常具有低消耗特性,其密度往往不超过50%。当密度在这个范围内逐渐增加时,脑网络的全局和局部效率呈现非线性增长趋势。当脑网络呈现出大于 50%的稀疏区间时,其网络结构将会变得随机,这种随机化可能导致网络性能的下降,因为脑网络的稳定性和信息传递能力受到了影响24。因此,在本研究中为确保脑网络密度小于50%,由式(4)计算网络密度:S=1N()N-1i=1Nki(4)式中ki为节点i的度。人类的大脑在信息传递过程中采用了小世界现象作为一种优化模式。研究表明,ASD患者的脑功能网络小世界属性低于普通人。小世界属性可通过式(5)计算得出,当结果中的 1时,说明该脑功能网络具有小世界拓扑结构25。=(5)
22、式中:=C Crandom,C为脑功能网络的聚类系数,由式(6)计算,Crandom=kN,为随机网络下的聚类系数,k为节点的平均度;=L Lrandom,Lrandom=ln N ln k,为随机网络下的平均路径长度,L为脑功能网络的平均路径长度,由式(7)计算,k为节点的度。C=1Ni=1NCi=1Ni=1N2Eiki()ki-1(6)50第21期L=1Ni=1NLi=1Ni=1Nj=1,j iNdijN-1(7)式中:节点i的聚类系数记为Ci;Ei为节点i与相邻节点间存在的连接边数;Li为节点i到其他所有节点的平均最短路径长度;dij为节点i和j的最短路径长度。1.4.4 Node2ve
23、c算法Node2vec是图嵌入节点转换为向量的算法,输入为网络结构,输出为每个节点对应的低维向量,而非高维onehot编码向量,如图3所示。该算法不仅最大程度地保留脑网络原有的非欧几里德空间结构特性,也便于深度学习模型提取其内在特征26。图3 Node2vec算法示意图Node2vec 是 介 于 深 度 优 先 搜 索(Depth First Search,DFS)和 宽 度 优 先 搜 索(Breadth First Search,BFS)之间的搜索算法,通过参数调整其搜索方式偏向DFS 或偏向 BFS,既能由 DFS 体现网络的结构性,也能由 BFS 体现网络的同质性。Node2vec
24、的核心算法为node2vecwalk,从节点v开始,按固定步长1游走到节点x的概论由式(8)计算:P()ci=|x ci-1=v=pq()v,x wvxZ(8)式中:ci表示当前节点;ci-1表示前一个节点;Z为归一化常数;wvx为有向图中边权重(本研究采用无向图,wvx=1);pq()v,x为元转移概率,由式(9)计算。pq()v,x=1p,dvx=01,dvx=11q,dvx=2(9)式中:dvx表示三种候选节点的类型,dvx=0表示候选节点为前一个时刻节点y,代表回头路,dvx=1代表候选节点x到v的距离与前一时刻节点y到v的距离相等,即BFS,dvx=2代表DFS;p和q代表方向参数。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 深度 学习 自闭症 谱系 障碍 fMRI 数据 分类 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。