CN117116434A - 人脑白质结构连接组的个体差异评估方法、应用及装置 - Google Patents

人脑白质结构连接组的个体差异评估方法、应用及装置 Download PDF

Info

Publication number
CN117116434A
CN117116434A CN202311385895.0A CN202311385895A CN117116434A CN 117116434 A CN117116434 A CN 117116434A CN 202311385895 A CN202311385895 A CN 202311385895A CN 117116434 A CN117116434 A CN 117116434A
Authority
CN
China
Prior art keywords
individual
brain
brain regions
individual difference
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202311385895.0A
Other languages
English (en)
Other versions
CN117116434B (zh
Inventor
舒妮
黄伟杰
刘祯钊
冯国政
陈豪杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Normal University
Original Assignee
Beijing Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Normal University filed Critical Beijing Normal University
Priority to CN202311385895.0A priority Critical patent/CN117116434B/zh
Publication of CN117116434A publication Critical patent/CN117116434A/zh
Application granted granted Critical
Publication of CN117116434B publication Critical patent/CN117116434B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/048Activation functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30016Brain
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30204Marker

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Molecular Biology (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Quality & Reliability (AREA)
  • Public Health (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

人脑白质结构连接组的个体差异评估方法、应用及装置,本发明通过图像采集、预处理、构建脑网络、去除中心效应、定义个体差异、获取结构连接个体差异与大脑其他结构属性的关系、获取结构连接个体差异与基因表达关系、基于激活函数的岭回归模型并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能,从而通过探索大脑结构连接的个体差异的时空模式,能够获得形成该时空模式可能的原因,是否受到其他组织结构的约束,是否受基因调控。

Description

人脑白质结构连接组的个体差异评估方法、应用及装置
技术领域
本发明涉及医学图像处理的技术领域,尤其涉及一种人脑白质结构连接组的个体差异评估方法,该方法的应用,以及人脑白质结构连接组的个体差异评估装置。
背景技术
每个个体都是独一无二的,他们的思考,行为,行动以及表达自我的方式都是不同的,而这些可能正是大脑功能与结构的个体差异的体现。最早通对尸检后大脑的肉眼观察,发现大脑皮层上一些小的沟回结构存在着明显的个体差异。随后,基于细胞染色技术发现在布罗德曼44区的细胞构筑特征也存在着明显的个体差异。而功能上,通过任务态fMRI(Functional magnetic resonance imaging,功能性磁共振成像)发现在重认测试的情景检索过程中,大脑激活脑区也存在显著的个体差异。
随后,基于全脑的观察研究发现大脑的结构和功能的个体差异在大脑中分布是异质的,即支持高阶关联和整合过程的神经系统在人群存在着明显的个体差异,而处理单模态信息的神经系统在人群中则较为一致。例如语言区域在细胞构筑定义的体积以及基于fMRI衍生的功能定位都表现出较高的个体差异。宏观尺度上,联络皮层沟回的个体差异明显高于运动皮层。此外,通过对白质纤维束的染色研究发现,长联合白质纤维束的拓扑地形特征比视辐射和皮质脊髓束的个体差异更大。
人脑连接组的提出为认识大脑神经系统提供了一个全新的视角,且为将大脑作为一个完整系统研究提供了有效的工具。大量研究表明人脑网络的个体差异可以解释健康被试在不同认知域上行为的变异,并且与神经系统和精神疾病患者的症状严重程度相关。近年来,一些学者开始探究大脑功能连接的个体差异模式。基于功能连接的研究发现,人脑功能连接的个体差异在大脑中也是非均匀分布的。个体差异在包括外侧前额叶和颞顶叶交界处的多模态联合皮层中最大,而在单模态感觉和运动皮层中个体差异最小。基于Yeo的7个特定的功能网络发现额顶控制和注意力网络表现出高度的功能变异性,而感觉运动和视觉系统的变异性最低,默认网络表现出中等水平的变异性,低于额顶和注意力网络,但高于感觉运动和视觉网络的变异性。通过对早产儿以及正常出生的婴儿研究发现,在发育早期(约33周),全脑都表现出非常高的功能连接个体差异,随着大脑逐渐发育,全脑都表现出个体差异下降的趋势,但是单模态区域如感觉运动皮层、视觉皮层个体差异下降最快,而多模态联合皮层的个体差异下降速率则较慢。这样的趋势一直延续到产后一周岁,才形成和成年人相似的个体差异空间模式。而Ma等人发现,老化过程中功能连接个体差异的变化趋势与发育过程正好相反,全脑的功能连接个体差异普遍上升,且在发育过程中个体差异下降最快的脑区上升速率最快。对于结构连接的个体差异的时空模式则研究较少,主要发现结构连接的个体差异远小于功能连接的个体差异,且它们的空间模式也并不对应。
因此,本研究旨在探索大脑结构连接的个体差异的时空模式,以及形成该时空模式可能的原因,是否受到其他组织结构的约束,是否受基因调控。最后通过关联个体差异的时空模式与认知的个体差异来探究大脑结构连接的个体差异的时空模式对于个体与群体的意义。
发明内容
为克服现有技术的缺陷,本发明要解决的技术问题是提供了一种人脑白质结构连接组的个体差异评估方法,其通过探索大脑结构连接的个体差异的时空模式,能够获得形成该时空模式可能的原因,是否受到其他组织结构的约束,是否受基因调控。
本发明的技术方案是:这种人脑白质结构连接组的个体差异评估方法,其包括以下步骤:
(1)图像采集:使用矢状三维磁化制备的快速梯度回波序列,获得覆盖全脑的高分辨率T1-MRI;使用磁共振扫描仪获得dMRI数据;
(2)对T1-MRI、dMRI数据进行预处理;
(3)构建脑网络:使用人脑连接组分区图谱将大脑分割成246个感兴趣区,这些区域代表结构网络中的节点,为了将图谱映射到个体空间,将在T1配准到ICBM模板时生成的变换取逆,将逆变换应用到在标准空间上的人脑连接组分区图谱得到个体空间的大脑分区,这一步使用最近邻插值方法来保留离散的标记值;对于每一对节点,如果至少有一条流线的两个端点位于这对节点中,则认为这对节点在结构上是相连的,将连接两个节点之间的纤维数目作为边的权重,最后在使用球形反卷积信息过滤的方法对边的权重进行校正,使获得的纤维连接更有生物学意义;
(4)使用ComBat算法来去除数据集的中心效应;
(5)定义个体差异:结构连接的被试内差异定义如下:
(3)
其中,表示第/>个脑区的结构连接的被试内差异,/>表示的是第/>个在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个在第二个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
然后定义被试间的个体差异如下:
(4)
其中,表示第/>个脑区的结构连接的被试间差异,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
最后定义结构连接的个体差异为被试间差异回归掉被试内差异的残差,定义如下:
(5)
其中和/>是需要拟合的线性模型参数;
(6)获取结构连接个体差异与大脑其他结构属性的关系,使用斯皮尔曼秩相关评估图谱之间的相关性,使用基于莫兰指数的空间置换来评价相关系数的显著性;
(7)获取结构连接个体差异与基因表达关系,包括:不同脑区基因表达估计、结构连接个体差异与基因表达的空间相关、基因富集分析、基于BrainSpan数据集的鲁棒性分析;
(8)基于激活函数的岭回归模型,并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能。
本发明通过图像采集、预处理、构建脑网络、去除中心效应、定义个体差异、获取结构连接个体差异与大脑其他结构属性的关系、获取结构连接个体差异与基因表达关系、基于激活函数的岭回归模型并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能,从而通过探索大脑结构连接的个体差异的时空模式,能够获得形成该时空模式可能的原因,是否受到其他组织结构的约束,是否受基因调控。
还提供了这种人脑白质结构连接组的个体差异评估方法应用于预测个体认知。
还提供了人脑白质结构连接组的个体差异评估装置,其包括:
图像采集模块,其配置来使用矢状三维磁化制备的快速梯度回波序列,获得覆盖全脑的高分辨率T1-MRI;使用磁共振扫描仪dMRI数据;
数据预处理模块,其配置来对T1-MRI、dMRI预处理;
脑网络构建模块,其配置来使用人脑连接组分区图谱将大脑分割成246个感兴趣区,这些区域代表结构网络中的节点,为了将图谱映射到个体空间,将在T1配准到ICBM模板时生成的变换取逆,将逆变换应用到在标准空间上的人脑连接组分区图谱得到个体空间的大脑分区,这一步使用最近邻插值方法来保留离散的标记值;对于每一对节点,如果至少有一条流线的两个端点位于这对节点中,则认为这对节点在结构上是相连的,将连接两个节点之间的纤维数目作为边的权重,最后在使用球形反卷积信息过滤的方法对边的权重进行校正,使获得的纤维连接更有生物学意义;
中心效应去除模块,其配置来使用ComBat算法来去除数据集的中心效应;
定义个体差异模块,其配置来定义个体差异:结构连接的被试内差异定义如下:
(3)
其中,表示第/>个脑区的结构连接的被试内差异,/>表示的是第/>个在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个在第二个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
然后定义被试间的个体差异如下:
(4)
其中,表示第/>个脑区的结构连接的被试间差异,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
最后定义结构连接的个体差异为被试间差异回归掉被试内差异的残差,定义如下:
(5)
其中和/>是需要拟合的线性模型参数;
第一获取模块,其配置来获取结构连接个体差异与大脑其他结构属性的关系,使用斯皮尔曼秩相关评估图谱之间的相关性,使用基于莫兰指数的空间置换来评价相关系数的显著性;
第二获取模块,其配置来获取结构连接个体差异与基因表达关系,包括:不同脑区基因表达估计、结构连接个体差异与基因表达的空间相关、基因富集分析、基于BrainSpan数据集的鲁棒性分析;
评估模块,其配置来基于激活函数的岭回归模型,并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能。
附图说明
图1是根据本发明的人脑白质结构连接组的个体差异评估方法的流程图。
具体实施方式
如图1所示,这种人脑白质结构连接组的个体差异评估方法,其包括以下步骤:
(1)图像采集:使用矢状三维磁化制备的快速梯度回波序列,获得覆盖全脑的高分辨率T1-MRI;使用磁共振扫描仪获得dMRI数据;
(2)对T1-MRI、dMRI数据进行预处理;
(3)构建脑网络:使用人脑连接组分区图谱将大脑分割成246个感兴趣区,这些区域代表结构网络中的节点,为了将图谱映射到个体空间,将在T1配准到ICBM模板时生成的变换取逆,将逆变换应用到在标准空间上的人脑连接组分区图谱得到个体空间的大脑分区,这一步使用最近邻插值方法来保留离散的标记值;对于每一对节点,如果至少有一条流线的两个端点位于这对节点中,则认为这对节点在结构上是相连的,将连接两个节点之间的纤维数目作为边的权重,最后在使用球形反卷积信息过滤的方法对边的权重进行校正,使获得的纤维连接更有生物学意义;
(4)使用ComBat算法来去除数据集的中心效应;
(5)定义个体差异:结构连接的被试内差异定义如下:
(3)
其中,表示第/>个脑区的结构连接的被试内差异,/>表示的是第/>个在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个在第二个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
然后定义被试间的个体差异如下:
(4)
其中,表示第/>个脑区的结构连接的被试间差异,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
最后定义结构连接的个体差异为被试间差异回归掉被试内差异的残差,定义如下:
(5)
其中和/>是需要拟合的线性模型参数;
(6)获取结构连接个体差异与大脑其他结构属性的关系,使用斯皮尔曼秩相关评估图谱之间的相关性,使用基于莫兰指数的空间置换来评价相关系数的显著性;
(7)获取结构连接个体差异与基因表达关系,包括:不同脑区基因表达估计、结构连接个体差异与基因表达的空间相关、基因富集分析、基于BrainSpan数据集的鲁棒性分析;
(8)基于激活函数的岭回归模型,并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能。
本发明通过图像采集、预处理、构建脑网络、去除中心效应、定义个体差异、获取结构连接个体差异与大脑其他结构属性的关系、获取结构连接个体差异与基因表达关系、基于激活函数的岭回归模型并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能,从而通过探索大脑结构连接的个体差异的时空模式,能够获得形成该时空模式可能的原因,是否受到其他组织结构的约束,是否受基因调控。
优选地,所述步骤(2)中,通过qsiprep方法处理T1-MRI和dMRI数据;对于T1-MRI数据,首先使用ANTs中的antsBrainExtraction命令移除T1-MRI中非脑组织部分并创建大脑的掩模图像,然后使用FSL中的fast命令将大脑组织分割为灰质、白质和脑脊液三部分,最后使用ANTs中的antsRegistration将T1-MRI图像配准到ICBM模板上;对于单壳dMRI数据,使用FSL中的eddy_correct进行头动和涡流校正;对于多壳dMRI数据,首先对图像进行降噪处理,然后使用FSL中的topup对图像进行磁敏感校正和eddy进行头动和涡流校正,最后基于多个b0图像创建平均的b0模板并与T1图像配准。
优选地,所述步骤(3)中,基于dMRI数据,首先使用MRtrix3估计每个体素的纤维方向分布函数,使用单壳的dMRI数据球形谐波去卷积算法估计纤维方向分布,而对于多壳的dMRI数据则使用多壳的球形谐波去卷积算法;然后使用概率追踪算法重建全脑纤维束;在追踪时,增加解剖结构限制,限制代表纤维束的流线只经过白质区域,且两个端点只落在灰质区域。
优选地,所述步骤(4)中,中心效应包含加性和乘性效应,对于每个特征建立了一个线性模型来拟合这些效应,该线性模型为:
(1)
其中是第/>个中心第/>个被试第/>个连接的连接强度,/>是第/>个连接的平均值,/>是可能会影响结构连接的协变量,/>为其对应的系数,/>是表示中心的变量,/>是对应的相关系数;假设残差项/>均值为0,但是每个中心残差项的标准差不同,/>则是第/>个中心第/>个连接对应的残差的尺度,使用经验贝叶斯估计模型参数/>与/>,调和后的连接定义如下:
(2)
这一步通过调用Python包neuroCombat实现。
优选地,所述步骤(6)中,使用欧式距离度量脑区之间的相似性,根据相似矩阵对大脑图谱进行10000次随机置换且使得置换后的图谱具有和原图谱相似的莫兰指数,保持了原始图谱的空间自相关性;然后使用斯皮尔曼秩相关计算结构连接个体差异图谱与置换后的图谱的相关系数得到一个零分布,通过计算零分布中大于真实相关系数的比例得到代表显著性的p值,这一过程使用Python包Brainspace实现。
优选地,所述步骤(7)中,
针对不同脑区基因表达估计,首先根据最新且最准确的探针基因对应关系将探针分配到每个基因上,选择测得的基因表达值与RNA序列技术测得的值最相关的探针作为代表探针;然后将采样样本分配到距离样本最近的脑区,如果一个样本距离所有脑区的距离都大于2 mm则舍弃该样本;最后对每个被试内所有样本使用缩放的且稳定的Sigmoid标准化方法控制个体差异,然后对每个脑区计算每个基因在6个被试上表达的平均值,得到最终的基因表达图谱,选择差异稳定性大于0.1的基因,最终得到11261个基因在左脑123个脑区的表达水平,差异稳定性定义如下:
(6)
其中是探针/>在第/>个被试所有脑区测得基因表达值与其在第/>个被试所有脑区测得的基因表达值的斯皮尔曼秩相关系数,/>是被试数量;
选择AHBA筛选后的基因与基于GTEx数据库鉴别基因表达显著高于身体的其他部位的2979个基因的交集,得到2113个基因在全脑不同脑区的表达;
针对结构连接个体差异与基因表达的空间相关,使用偏最小二乘分析探究结构连接个体差异与基因表达关系;
针对基因富集分析,分别对正向与负向权重最大的前50%基因进行基因富集分析,对于每一个生物过程类别,将该生物过程涉及的基因的平均权重作为富集分数;然后,使用基于Moran指数的空间置换对结构连接个体差异图谱进行10000次置换,对每一次置换得到的结构连接个体差异图谱进行PLS和基因富集分析;将真实的生物过程类别的富集分数与10000次置换得到的富集分数的零分布进行比较,显著性p值为零分布中高于真实富集分数的比例;
针对基于BrainSpan数据集的鲁棒性分析,估计16个一致的大脑区域在五个生命阶段的基因分数,为了与来自AHBA的PLS分析结果进行比较,通过将基因表达矩阵和PLS分析衍生的基因权重相乘来估计16个区域的基因分数;为了将估计的基因分数与结构连接个体差异联系起来,定义一个从16节点分区到123节点分区的区域对应关系,对于这16个脑区,计算其对应的脑区的结构连接的个体差异的平均值作为其对应的个体差异。
优选地,所述步骤(8)中,基于激活函数的岭回归模型为两层神经网络,第一层为特征的线性变换,然后将变换后的特征输入到线性整流激活函数中,通过激活函数的特征再被输入到第二层的岭回归模型中,该模型表示为:
(11)
其中和/>为神经网络中第一层特征线性变换参数,/>和/>分别为第二层岭回归模型的回归系数与截距,/>是线性整流函数,训练时,该模型的代价函数与岭回归的代价函数相同,但是使用反向传播更新模型参数。
优选地,所述步骤(8)中,针对嵌套交叉验证,将数据划分成五个子集,在每一折外层的交叉验证中,选择一个子集作为测试集,剩余四个子集作为训练集,训练集特征的均值和标准差被用来标准化训练集和测试集的特征,然后基于内层的四折交叉验证选择的最优超参数,然后我其应用到整个训练集训练一个最终的模型,并对测试集进行预测,基于预测结果,预测的认知分数与实际认知分数的平均绝对值误差和皮尔逊相关系数被用来评价模型准确性;
在内层的交叉验证中,首先从[2-5,2-4,…,24,25]选择一个值作为超参数,然后基于该超参数进行四折交叉验证,每一折中,计算测试集上预测的认知分数和实际认知分数的皮尔逊相关系数评价模型准确性,然后计算四折的皮尔逊相关系数的平均值用以评估给定超参数/>的模型准确性,选择使模型准确性最高的/>作为最优的/>
使用置换检验来评价模型表现是否显著优于随机模型,按照预测框架执行1000次模型训练与预测,在每次执行时,将训练集的标签随机置换,对每一个认知功能的每一折交叉验证中都获得1000次准确性评价;通过比较真实模型的准确性评价与置换检验得到的1000次模型评价组成的分布获得显著性p值,皮尔逊相关系数的p值是置换检验获得的相关系数中高于真实模型的比例,平均绝对值误差的p值是置换检验获得的平均绝对值误差中低于真实模型的比例;
对于每一个被试,定义被试特异的特征激活图谱度量该被试每一个特征对于预测的贡献,特征激活图谱定义如下:
(12)
其中是第/>个被试第/>个特征对于预测模型的贡献值,对于每一个被试,将五折交叉验证得到的五个模型的特征激活图谱进行平均,得到最终的平均特征激活图谱评价该被试的每一个特征对于预测模型的最终贡献;最后对于三个不同认知域,使用结构连接个体差异计算方法计算特征激活图谱的个体间差异,然后使用斯皮尔曼秩相关评价结构连接个体差异与认知预测模型权重的关系。
还提供了这种人脑白质结构连接组的个体差异评估方法应用于预测个体认知。本发明研发了一个稳定且可重复的白质结构连接个体差异时空模式。并发现白质连接个体差异的空间结构与大脑的可塑性有关。该可塑性一方面与大脑的微观结构有关,另一方面与基因的选择性表达有关。最后基于大脑白质连接个体差异的特定模式,创新性提出了一种基于激活函数的岭回归模型,能够使用白质连接准确的预测个体认知。该模型一方面表明结构连接的个体差异可能是受个体不同的认知和行为塑造的结果,另一方面为未来的大脑认知关联研究提供启示,应当考虑到脑网络中同一连接在不同个体中对于认知的参与程度可能存在差异。
本领域普通技术人员可以理解,实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,包括上述实施例方法的各步骤,而所述的存储介质可以是:ROM/RAM、磁碟、光盘、存储卡等。因此,与本发明的方法相对应的,本发明还同时包括一种人脑白质结构连接组的个体差异评估装置,该装置通常以与方法各步骤相对应的功能模块的形式表示。该装置包括:
图像采集模块,其配置来使用矢状三维磁化制备的快速梯度回波序列,获得覆盖全脑的高分辨率T1-MRI与dMRI数据;
数据预处理模块,其配置来对T1-MRI、dMRI数据进行预处理;
脑网络构建模块,其配置来使用人脑连接组分区图谱将大脑分割成246个感兴趣区,这些区域代表结构网络中的节点,为了将图谱映射到个体空间,将在T1配准到ICBM模板时生成的变换取逆,将逆变换应用到在标准空间上的人脑连接组分区图谱得到个体空间的大脑分区,这一步使用最近邻插值方法来保留离散的标记值;对于每一对节点,如果至少有一条流线的两个端点位于这对节点中,则认为这对节点在结构上是相连的,将连接两个节点之间的纤维数目作为边的权重,最后在使用球形反卷积信息过滤的方法对边的权重进行校正,使获得的纤维连接更有生物学意义;
中心效应去除模块,其配置来使用ComBat算法来去除数据集的中心效应;
定义个体差异模块,其配置来定义个体差异:结构连接的被试内差异定义如下:
(3)
其中,表示第/>个脑区的结构连接的被试内差异,/>表示的是第/>个在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个在第二个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
然后定义被试间的个体差异如下:
(4)
其中,表示第/>个脑区的结构连接的被试间差异,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
最后定义结构连接的个体差异为被试间差异回归掉被试内差异的残差,定义如下:
(5)
其中和/>是需要拟合的线性模型参数;
第一获取模块,其配置来获取结构连接个体差异与大脑其他结构属性的关系,使用斯皮尔曼秩相关评估图谱之间的相关性,使用基于莫兰指数的空间置换来评价相关系数的显著性;
第二获取模块,其配置来获取结构连接个体差异与基因表达关系,包括:不同脑区基因表达估计、结构连接个体差异与基因表达的空间相关、基因富集分析、基于BrainSpan数据集的鲁棒性分析;
评估模块,其配置来基于激活函数的岭回归模型,并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属本发明技术方案的保护范围。

Claims (10)

1.人脑白质结构连接组的个体差异评估方法,其特征在于:其包括以下步骤:
(1)图像采集:使用矢状三维磁化制备的快速梯度回波序列,获得覆盖全脑的高分辨率T1-MRI;使用磁共振扫描仪获得dMRI;
(2)对T1-MRI、dMRI数据进行预处理进行预处理;
(3)构建脑网络:使用人脑连接组分区图谱将大脑分割成246个感兴趣区,这些区域代表结构网络中的节点,为了将图谱映射到个体空间,将在T1配准到ICBM模板时生成的变换取逆,将逆变换应用到在标准空间上的人脑连接组分区图谱得到个体空间的大脑分区,这一步使用最近邻插值方法来保留离散的标记值;对于每一对节点,如果至少有一条流线的两个端点位于这对节点中,则认为这对节点在结构上是相连的,将连接两个节点之间的纤维数目作为边的权重,最后在使用球形反卷积信息过滤的方法对边的权重进行校正,使获得的纤维连接更有生物学意义;
(4)使用ComBat算法来去除数据集的中心效应;
(5)定义个体差异:结构连接的被试内差异定义如下:
(3)
其中,表示第/>个脑区的结构连接的被试内差异,/>表示的是第/>个在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个在第二个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
然后定义被试间的个体差异如下:
(4)
其中,表示第/>个脑区的结构连接的被试间差异,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
最后定义结构连接的个体差异为被试间差异回归掉被试内差异的残差,定义如下:
(5)
其中和/>是需要拟合的线性模型参数;
(6)获取结构连接个体差异与大脑其他结构属性的关系,使用斯皮尔曼秩相关评估图谱之间的相关性,使用基于莫兰指数的空间置换来评价相关系数的显著性;
(7)获取结构连接个体差异与基因表达关系,包括:不同脑区基因表达估计、结构连接个体差异与基因表达的空间相关、基因富集分析、基于BrainSpan数据集的鲁棒性分析;
(8)基于激活函数的岭回归模型,并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能。
2.根据权利要求1所述的人脑白质结构连接组的个体差异评估方法,其特征在于:所述步骤(2)中,通过qsiprep方法处理T1-MRI和dMRI数据;对于T1-MRI数据,首先使用ANTs中的antsBrainExtraction命令移除T1-MRI中非脑组织部分并创建大脑的掩模图像,然后使用FSL中的fast命令将大脑组织分割为灰质、白质和脑脊液三部分,最后使用ANTs中的antsRegistration将T1-MRI图像配准到ICBM模板上;对于dMRI数据,首先对图像进行降噪处理,然后使用FSL中的topup对图像进行磁敏感校正和eddy进行头动和涡流校正,最后基于多个b0图像创建平均的b0模板并与T1图像配准。
3.根据权利要求2所述的人脑白质结构连接组的个体差异评估方法,其特征在于:所述步骤(3)中,对于单壳dMRI数据,首先使用MRtrix3估计每个体素的纤维方向分布函数,使用单壳的球形谐波去卷积算法估计纤维方向分布,而对于多壳dMRI数据则使用多壳的球形谐波去卷积算法;然后使用概率追踪算法重建全脑纤维束;在追踪时,增加解剖结构限制,限制代表纤维束的流线只经过白质区域,且两个端点只落在灰质区域。
4.根据权利要求3所述的人脑白质结构连接组的个体差异评估方法,其特征在于:所述步骤(4)中,中心效应包含加性和乘性效应,对于每个特征建立了一个线性模型来拟合这些效应,该线性模型为:
(1)
其中是第/>个中心第/>个被试第/>个连接的连接强度,/>是第/>个连接的平均值,/>是可能会影响结构连接的协变量,/>为其对应的系数,/>是表示中心的变量,/>是对应的相关系数;假设残差项/>均值为0,但是每个中心残差项的标准差不同,/>则是第/>个中心第/>个连接对应的残差的尺度,使用经验贝叶斯估计模型参数/>与/>,调和后的连接定义如下:
(2)
这一步通过调用Python包neuroCombat实现。
5.根据权利要求4所述的人脑白质结构连接组的个体差异评估方法,其特征在于:所述步骤(6)中,使用欧式距离度量脑区之间的相似性,根据相似矩阵对大脑图谱进行10000次随机置换且使得置换后的图谱具有和原图谱相似的莫兰指数,保持了原始图谱的空间自相关性;然后使用斯皮尔曼秩相关计算结构连接个体差异图谱与置换后的图谱的相关系数得到一个零分布,通过计算零分布中大于真实相关系数的比例得到代表显著性的p值,这一过程使用Python包Brainspace实现。
6.根据权利要求5所述的人脑白质结构连接组的个体差异评估方法,其特征在于:所述步骤(7)中,
针对不同脑区基因表达估计,首先根据最新且最准确的探针基因对应关系将探针分配到每个基因上,选择测得的基因表达值与RNA序列技术测得的值最相关的探针作为代表探针;然后将采样样本分配到距离样本最近的脑区,如果一个样本距离所有脑区的距离都大于2 mm则舍弃该样本;最后对每个被试内所有样本使用缩放的且稳定的Sigmoid标准化方法控制个体差异,然后对每个脑区计算每个基因在6个被试上表达的平均值,得到最终的基因表达图谱,选择差异稳定性大于0.1的基因,最终得到11261个基因在左脑123个脑区的表达水平,差异稳定性定义如下:
(6)
其中是探针/>在第/>个被试所有脑区测得基因表达值与其在第/>个被试所有脑区测得的基因表达值的斯皮尔曼秩相关系数,/>是被试数量;
选择AHBA筛选后的基因与基于GTEx数据库鉴别基因表达显著高于身体的其他部位的2979个基因的交集,得到2113个基因在全脑不同脑区的表达;
针对结构连接个体差异与基因表达的空间相关,使用偏最小二乘分析探究结构连接个体差异与基因表达关系;
针对基因富集分析,分别对正向与负向权重最大的前50%基因进行基因富集分析,对于每一个生物过程类别,将该生物过程涉及的基因的平均权重作为富集分数;然后,使用基于Moran指数的空间置换对结构连接个体差异图谱进行10000次置换,对每一次置换得到的结构连接个体差异图谱进行PLS和基因富集分析;将真实的生物过程类别的富集分数与10000次置换得到的富集分数的零分布进行比较,显著性p值为零分布中高于真实富集分数的比例;
针对基于BrainSpan数据集的鲁棒性分析,估计16个一致的大脑区域在五个生命阶段的基因分数,为了与来自AHBA的PLS分析结果进行比较,通过将基因表达矩阵和PLS分析衍生的基因权重相乘来估计16个区域的基因分数;为了将估计的基因分数与结构连接个体差异联系起来,定义一个从16节点分区到123节点分区的区域对应关系,对于这16个脑区,计算其对应的脑区的结构连接的个体差异的平均值作为其对应的个体差异。
7.根据权利要求6所述的人脑白质结构连接组的个体差异评估方法,其特征在于:所述步骤(8)中,基于激活函数的岭回归模型为两层神经网络,第一层为特征的线性变换,然后将变换后的特征输入到线性整流激活函数中,通过激活函数的特征再被输入到第二层的岭回归模型中,该模型表示为:
(11)
其中和/>为神经网络中第一层特征线性变换参数,/>和/>分别为第二层岭回归模型的回归系数与截距,/>是线性整流函数,训练时,该模型的代价函数与岭回归的代价函数相同,但是使用反向传播更新模型参数。
8.根据权利要求7所述的人脑白质结构连接组的个体差异评估方法,其特征在于:所述步骤(8)中,针对嵌套交叉验证,将数据划分成五个子集,在每一折外层的交叉验证中,选择一个子集作为测试集,剩余四个子集作为训练集,训练集特征的均值和标准差被用来标准化训练集和测试集的特征,然后基于内层的四折交叉验证选择的最优超参数,然后我其应用到整个训练集训练一个最终的模型,并对测试集进行预测,基于预测结果,预测的认知分数与实际认知分数的平均绝对值误差和皮尔逊相关系数被用来评价模型准确性;
在内层的交叉验证中,首先从[2-5,2-4,…,24,25]选择一个值作为超参数,然后基于该超参数进行四折交叉验证,每一折中,计算测试集上预测的认知分数和实际认知分数的皮尔逊相关系数评价模型准确性,然后计算四折的皮尔逊相关系数的平均值用以评估给定超参数/>的模型准确性,选择使模型准确性最高的/>作为最优的/>
使用置换检验来评价模型表现是否显著优于随机模型,按照预测框架执行1000次模型训练与预测,在每次执行时,将训练集的标签随机置换,对每一个认知功能的每一折交叉验证中都获得1000次准确性评价;通过比较真实模型的准确性评价与置换检验得到的1000次模型评价组成的分布获得显著性p值,皮尔逊相关系数的p值是置换检验获得的相关系数中高于真实模型的比例,平均绝对值误差的p值是置换检验获得的平均绝对值误差中低于真实模型的比例;
对于每一个被试,定义被试特异的特征激活图谱度量该被试每一个特征对于预测的贡献,特征激活图谱定义如下:
(12)
其中是第/>个被试第/>个特征对于预测模型的贡献值,对于每一个被试,将五折交叉验证得到的五个模型的特征激活图谱进行平均,得到最终的平均特征激活图谱评价该被试的每一个特征对于预测模型的最终贡献;最后对于三个不同认知域,使用结构连接个体差异计算方法计算特征激活图谱的个体间差异,然后使用斯皮尔曼秩相关评价结构连接个体差异与认知预测模型权重的关系。
9.根据权利要求8所述的人脑白质结构连接组的个体差异评估方法应用于预测个体认知。
10.人脑白质结构连接组的个体差异评估装置,其特征在于:其包括:
图像采集模块,其配置来使用矢状三维磁化制备的快速梯度回波序列,获得覆盖全脑的高分辨率T1-MRI;使用磁共振扫描仪获得dMRI数据;
数据预处理模块,其配置来对T1-MRI、dMRI数据进行预处理;
脑网络构建模块,其配置来使用人脑连接组分区图谱将大脑分割成246个感兴趣区,这些区域代表结构网络中的节点,为了将图谱映射到个体空间,将在T1配准到ICBM模板时生成的变换取逆,将逆变换应用到在标准空间上的人脑连接组分区图谱得到个体空间的大脑分区,这一步使用最近邻插值方法来保留离散的标记值;对于每一对节点,如果至少有一条流线的两个端点位于这对节点中,则认为这对节点在结构上是相连的,将连接两个节点之间的纤维数目作为边的权重,最后在使用球形反卷积信息过滤的方法对边的权重进行校正,使获得的纤维连接更有生物学意义;
中心效应去除模块,其配置来使用ComBat算法来去除数据集的中心效应;
定义个体差异模块,其配置来定义个体差异:结构连接的被试内差异定义如下:
(3)
其中,表示第/>个脑区的结构连接的被试内差异,/>表示的是第/>个在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个在第二个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
然后定义被试间的个体差异如下:
(4)
其中,表示第/>个脑区的结构连接的被试间差异,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>表示的是第/>个被试在第一个时间点的第/>个脑区与其他脑区的结构连接,/>代表被试数量;
最后定义结构连接的个体差异为被试间差异回归掉被试内差异的残差,定义如下:
(5)
其中和/>是需要拟合的线性模型参数;
第一获取模块,其配置来获取结构连接个体差异与大脑其他结构属性的关系,使用斯皮尔曼秩相关评估图谱之间的相关性,使用基于莫兰指数的空间置换来评价相关系数的显著性;
第二获取模块,其配置来获取结构连接个体差异与基因表达关系,包括:不同脑区基因表达估计、结构连接个体差异与基因表达的空间相关、基因富集分析、基于BrainSpan数据集的鲁棒性分析;
评估模块,其配置来基于激活函数的岭回归模型,并使用嵌套交叉验证评估人脑结构连接模式能否被用来预测个体的执行、注意和记忆功能。
CN202311385895.0A 2023-10-25 2023-10-25 人脑白质结构连接组的个体差异评估方法、应用及装置 Active CN117116434B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311385895.0A CN117116434B (zh) 2023-10-25 2023-10-25 人脑白质结构连接组的个体差异评估方法、应用及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311385895.0A CN117116434B (zh) 2023-10-25 2023-10-25 人脑白质结构连接组的个体差异评估方法、应用及装置

Publications (2)

Publication Number Publication Date
CN117116434A true CN117116434A (zh) 2023-11-24
CN117116434B CN117116434B (zh) 2024-03-01

Family

ID=88797048

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311385895.0A Active CN117116434B (zh) 2023-10-25 2023-10-25 人脑白质结构连接组的个体差异评估方法、应用及装置

Country Status (1)

Country Link
CN (1) CN117116434B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117316293A (zh) * 2023-11-30 2023-12-29 北京师范大学 基于神经影像的脑网络结构与功能耦合方法及装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140058241A1 (en) * 2012-08-22 2014-02-27 Ross Apparies Method and Apparatus for Assessing Neurocognitive Status
US20150088024A1 (en) * 2012-03-19 2015-03-26 University Of Florida Research Foundation, Inc. Methods and systems for brain function analysis
CN106021949A (zh) * 2016-05-31 2016-10-12 常州大学 一种大脑默认网络的功能连接分析方法
US20190206057A1 (en) * 2016-09-13 2019-07-04 Ohio State Innovation Foundation Systems and methods for modeling neural architecture
CN113610808A (zh) * 2021-08-09 2021-11-05 中国科学院自动化研究所 基于个体脑连接图的群体脑图谱个体化方法、系统和设备
US20220020154A1 (en) * 2020-07-15 2022-01-20 Siemens Healthcare Gmbh Method and system for characterizing an impact of brain lesions on brain connectivity using mri
CN115829054A (zh) * 2022-12-01 2023-03-21 华南理工大学 基于机器学习模型对脑网络异常边连接的预测方法
CN116784820A (zh) * 2022-11-12 2023-09-22 上海大学 一种基于种子点连接的脑功能网络构建方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150088024A1 (en) * 2012-03-19 2015-03-26 University Of Florida Research Foundation, Inc. Methods and systems for brain function analysis
US20140058241A1 (en) * 2012-08-22 2014-02-27 Ross Apparies Method and Apparatus for Assessing Neurocognitive Status
CN106021949A (zh) * 2016-05-31 2016-10-12 常州大学 一种大脑默认网络的功能连接分析方法
US20190206057A1 (en) * 2016-09-13 2019-07-04 Ohio State Innovation Foundation Systems and methods for modeling neural architecture
US20220020154A1 (en) * 2020-07-15 2022-01-20 Siemens Healthcare Gmbh Method and system for characterizing an impact of brain lesions on brain connectivity using mri
CN113610808A (zh) * 2021-08-09 2021-11-05 中国科学院自动化研究所 基于个体脑连接图的群体脑图谱个体化方法、系统和设备
CN116784820A (zh) * 2022-11-12 2023-09-22 上海大学 一种基于种子点连接的脑功能网络构建方法及系统
CN115829054A (zh) * 2022-12-01 2023-03-21 华南理工大学 基于机器学习模型对脑网络异常边连接的预测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
何李;李?;庄恺祥;陈群林;孙江洲;杨文静;位东涛;邱江;: "创造性的大脑网络连接特征与研究展望", 科学通报, no. 01, 10 January 2020 (2020-01-10), pages 31 - 42 *
何李等: "创造性的大脑网络连接特征与研究展望", 科学通报, no. 01, pages 31 - 42 *
梁夏;王金辉;贺永;: "人脑连接组研究:脑结构网络和脑功能网络", 科学通报, no. 16, 5 June 2010 (2010-06-05), pages 41 - 59 *
梁夏;王金辉;贺永;: "人脑连接组研究:脑结构网络和脑功能网络", 科学通报, no. 16, pages 41 - 59 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117316293A (zh) * 2023-11-30 2023-12-29 北京师范大学 基于神经影像的脑网络结构与功能耦合方法及装置
CN117316293B (zh) * 2023-11-30 2024-04-19 北京师范大学 基于神经影像的脑网络结构与功能耦合方法及装置

Also Published As

Publication number Publication date
CN117116434B (zh) 2024-03-01

Similar Documents

Publication Publication Date Title
CN113040715B (zh) 一种基于卷积神经网络的人脑功能网络分类方法
CN109376751B (zh) 一种基于卷积神经网络的人脑功能网络分类方法
CN109528197B (zh) 基于脑功能图谱进行精神疾病的个体化预测方法和系统
CN110097968B (zh) 基于静息态功能磁共振影像的婴儿脑龄预测方法、系统
CN109589092B (zh) 基于集成学习的阿尔茨海默症确定方法及系统
CN113255728B (zh) 一种基于图嵌入和多模态脑网络的抑郁症分类方法
CN111009324B (zh) 脑网络多特征分析的轻度认知障碍辅助诊断系统及方法
CN117116434B (zh) 人脑白质结构连接组的个体差异评估方法、应用及装置
CN111863244B (zh) 基于稀疏池化图卷积的功能连接精神疾病分类方法和系统
CN112418337B (zh) 基于脑功能超网络模型的多特征融合数据分类方法
CN111728590A (zh) 基于动态功能连接的个体认知能力预测方法和系统
CN110598793A (zh) 一种大脑功能网络特征分类方法
CN111063442A (zh) 基于弱监督多任务矩阵补全的脑疾病进程预测方法及系统
CN111539509B (zh) 一种虚拟数字脑构建方法与系统及智能机器人控制系统
CN110148145A (zh) 一种融合边界信息的图像目标区提取方法及应用
CN115294054A (zh) 一种基于大脑功能网络特征的认知功能评分预测方法
CN112614126A (zh) 基于机器学习的磁共振图像脑区划分方法、系统和装置
CN111513717A (zh) 一种大脑功能状态提取方法
CN112861879A (zh) 基于双重迁移的认知障碍辅助决策支持方法和系统
CN115496953A (zh) 基于时空图卷积的脑网络分类方法
CN112233805A (zh) 基于多图谱神经影像学数据进行生物标志物的挖掘方法
CN111488934A (zh) 脑影像数据处理方法、存储介质、计算机设备和装置
CN113197578A (zh) 一种基于多中心模型的精神分裂症分类方法及系统
CN116369891A (zh) 轻度认知障碍发展进程预测方法、装置和计算机设备
CN113283518B (zh) 一种基于聚类的多模态大脑网络特征选择方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant