CN110718301B - 基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法 - Google Patents
基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法 Download PDFInfo
- Publication number
- CN110718301B CN110718301B CN201910916563.8A CN201910916563A CN110718301B CN 110718301 B CN110718301 B CN 110718301B CN 201910916563 A CN201910916563 A CN 201910916563A CN 110718301 B CN110718301 B CN 110718301B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- dynamic
- functional nuclear
- brain function
- 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.)
- Active
Links
- 230000003925 brain function Effects 0.000 title claims abstract description 130
- 238000003745 diagnosis Methods 0.000 title claims abstract description 66
- 208000024827 Alzheimer disease Diseases 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 34
- 210000004556 brain Anatomy 0.000 claims abstract description 111
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 238000012549 training Methods 0.000 claims abstract description 36
- 238000012216 screening Methods 0.000 claims abstract description 24
- 238000007781 pre-processing Methods 0.000 claims abstract description 18
- 238000002599 functional magnetic resonance imaging Methods 0.000 claims abstract description 17
- 238000005481 NMR spectroscopy Methods 0.000 claims description 207
- 238000012937 correction Methods 0.000 claims description 43
- 238000001914 filtration Methods 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 239000003638 chemical reducing agent Substances 0.000 claims description 12
- 238000012706 support-vector machine Methods 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 4
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000003068 static effect Effects 0.000 abstract description 3
- 230000007177 brain activity Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000009499 grossing Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000007635 classification algorithm Methods 0.000 description 2
- 238000004195 computer-aided diagnosis Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005714 functional activity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 238000002610 neuroimaging Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 238000000700 time series analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/38—Registration of image sequences
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Epidemiology (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公布了一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法。诊断装置包括fMRI数据预处理单元、构建动态脑功能网络单元、生成用于训练的特征单元以及svm分类辅助诊断单元,该诊断装置的使用方法为:首先进行图像预处理,然后构建动态脑网络,其次计算分割后的脑网络的节点度量,并通过时间序列生成器将每个节点度量构成一个时间序列,随后通过特征提取器为构成的时间序列提取特征,再通过特征过滤器将过滤后的特征拼接成一个矩阵并通过特征筛选器筛选,最后通过数据训练器进行数据的分类训练,最终通过辅助诊断器实现对阿尔茨海默病的诊断。该方法克服了静态脑功能网络无法表示动态信息的缺陷,起到了更好的为医疗辅助诊断服务的效果。
Description
技术领域
本发明属于计算机辅助诊断技术领域,涉及一种基于支持向量机分类算法的阿尔茨海默病辅助诊断装置及方法,特别涉及一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法。
背景技术
近些年来,神经影像技术特别是功能影像学取得了飞速的进步,这为研究患有阿尔茨海默病病人各个脑区具有的功能提供了相应的技术。功能磁共振成像作为其中一种成熟的功能影像学检测技术,它的优势在于无创监测大脑功能和活动,时空分辨率较高。使用基于血氧水平依赖方法来测量大脑各个区域之间的相关性已被证明是研究大脑功能组织的有力工具。在静息态fMRI中,可以看出人脑的各个脑区在此时存在着有序的功能活动。
大脑是一个动态结构,随着时间的变化,大脑里的神经元之间的连接是变化的,基于动态网络的方法对脑功能连接网络进行分析,可以更好的分析出fMRI数据中每个时间段的瞬时特性。通过对动态脑功能网络的构建和分析,可以更好地描述大脑的活动状态以及各个神经元或脑区之间的交互。
发明内容
针对现有技术的不足,本发明的目的是提出基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法,利用支持向量机分类算法技术与构建动态脑功能网络技术对功能核磁共振图像(functional Magnetic Resonance Imaging,fMRI)进行脑网络的分类,进而能够有效的对阿尔茨海默病进行诊断,这种方法可以获取更多大脑动态活动的信息,并且可以准确地判断不同状态下脑网络之间的差异性,具体方案如下:
一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置包括fMRI数据预处理单元、构建动态脑功能网络单元、生成用于训练的特征单元以及svm分类辅助诊断单元,首先将获取的i个待测功能核磁共振图像通过所述fMRI数据预处理单元得到预处理后的的i个标准功能核磁共振图像,然后将预处理后的i个标准功能核磁共振图像通过所述构建动态脑功能网络单元进行时间维度上的分割,构建出i个功能核磁共振图像的i*s个动态脑功能网络,其次将构建出的i*s个动态脑功能网络通过所述生成用于训练的特征单元横向的提取出每个动态脑功能网络的动态特征,最后将从i*s个动态脑功能网络提取出的经过特征过滤器过滤后的动态特征拼接成一个动态特征矩阵,最后通过所述svm分类辅助诊断单元使用Fisher算法进行特征筛选,筛选出具有代表性的特征,并进行训练,用于辅助诊断。
所述的fMRI数据预处理单元包括时间片校正器、头动校正器、空间标准化器、平滑降噪器,首先将获取的i个待测功能核磁共振图像通过所述时间片校正器进行时间片校正,得到i个时间片校正后的功能核磁共振图像,然后将所述i个时间校正后的功能核磁共振图像通过所述头动校正器进行头动校正,得到i个头动校正后的功能核磁共振图像,其次将所述i个头动校正后的功能核磁共振图像通过所述空间标准化器进行空间标准化,得到i个空间标准化后的功能核磁共振图像,最后将所述i个空间标准化后的功能核磁共振图像通过所述平滑降噪器进行平滑降噪,得到i个标准功能核磁共振图像;
所述时间片校正器用于将输入的i个待测功能核磁图像进行时间片校正,得到i个时间片校正后的功能核磁共振图像(I-1,I-2,I-3,…,I-i),其中i表示选取的待测功能核磁共振图像的个数;
所述头动校正器用于将i个时间片校正后的功能核磁图像(I-1,I-2,I-3,…,I-i)进行头动校正,得到i个头动校正后的功能核磁共振图像(H-1,H-2,H-3,...,H-i);
所述空间标准化器用于将i个头动校正后的功能核磁图像(H-1,H-2,H-3,...,H-i)进行空间标准化,得到i个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,...,F-i);
所述平滑降噪器用于将i个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,…,F-i)进行平滑降噪,得到i个标准功能核磁共振图像(S-1,S-2,S-3,…,S-i)。
所述构建动态脑功能网络单元包括模板匹配器、时间序列分割器以及脑网络构建器,首先将i个标准功能核磁共振图像中的每个标准功能核磁共振图像通过所述模板匹配器与规格为m的标准匹配模板进行感兴趣区域的匹配,则每个匹配后的功能核磁共振图像包含m个脑区,然后将匹配后的功能核磁共振图像通过所述时间序列分割器将每个脑区的时间序列进行分割,得到s段瞬时时间序列,最后将分割得到的s段瞬时时间序列通过所述脑网络构建器构建出i个分割后的功能核磁共振图像的i*s个动态脑功能网络;
所述模板匹配器用于将每个标准功能核磁共振图像与规格为m的标准匹配模板进行匹配,匹配后的每个功能核磁共振图像包含m个脑区,i个匹配后的功能核磁共振图像表示为(A-1,A-2,A-3,...,A-i);
所述时间序列分割器用于将每个匹配后的功能核磁共振图像中的m个脑区的时间序列进行分割,将每个脑区的时间序列分割为s段瞬时时间序列,每段瞬时时间序列代表一个脑区的瞬时信息(T-1,T-2,...,T-s),i个分割后的功能核磁共振图像表示为其中s的取值范围根据预设的分割间隔确定;
所述脑网络构建器用于将每个分割后的功能核磁共振图像中的s段瞬时时间序列构建出s个动态脑功能网络,则i个分割后的功能核磁共振图像得到i*s个动态脑功能网络
所述生成用于训练的特征单元包括节点度量生成器、时间序列生成器、特征提取器以及特征过滤器,首先将构建动态脑功能网络单元输出的i个分割后的功能核磁共振图像的i*s个动态脑功能网络通过所述节点度量生成器计算i*s个动态脑功能网络的节点度量,然后通过时间序列生成器将i*s个动态脑功能网络中的每个节点度量构成一个时间序列,其次通过所述特征提取器为每个节点度量形成的时间序列提取新特征值,最后将提取到的所有的新特征值通过所述特征过滤器过滤后拼接为一个动态特征矩阵;
所述节点度量生成器用于生成i*s个动态脑功能网络的节点度量,并计算得到z个动态特征的特征值表示为所述z个动态特征包括i*s*x个全局特征和i*s*m*y个局部特征,即z=i*s*x+i*s*m*y,x表示每个动态脑功能网络计算的全局特征的个数,y表示每个动态脑功能网络计算的局部特征的个数;
所述时间序列生成器用于将得到的每组特征值生成一个时间序列,则i*s个动态脑功能网络构成的z/s个时间序列表示为所述每组特征值包括全局的每组特征值和局部的每组特征值,所述全局的每组特征值包括每个分割后的功能核磁共振图像中的每个全局特征的s个动态脑功能网络的特征值,所述局部的每组特征值包括从每个分割后的功能核磁共振图像中的s*m*y个局部特征中的每个局部特征按照m个脑区分组得到的s个动态脑功能网络的特征值;
所述特征提取器用于对每组特征生成的时间序列,基于小波的时间序列熵再次提取特征,得到z/s个特征的新特征值
所述特征过滤器用于对所述节点度量中的i*s*m*y个局部特征进行过滤,首先通过节点度量中的degree方法计算出每个匹配后的功能核磁共振图像中的每个脑区的degree,并计算出每个匹配后的功能核磁共振图像中的m个脑区degree的平均值然后计算出每个匹配后的功能核磁共振图像中的m个脑区的标准方差σ,其次将所述m个脑区中degree属于区间/>的脑区过滤为关键脑区,最后将i个匹配后的功能核磁共振图像过滤得到的N个脑区的N*y个局部特征生成的新特征值与i*x个全局特征生成的新特征值通过特征过滤器拼接为一个动态特征矩阵,/>其中nj表示第j个匹配后的功能核磁共振图像中的关键脑区的个数,i表示匹配后的功能核磁共振图像的个数。
所述的svm分类辅助诊断单元包括特征筛选器、数据训练器、辅助诊断器,首先将生成用于训练的特征单元输出的动态特征矩阵通过特征筛选器,使用Fisher算法对所述动态特征矩阵进行动态特征筛选,然后通过数据训练器进行分类训练,最终通过辅助诊断器实现对阿尔茨海默病的诊断;
所述特征筛选器用于将特征过滤器提取的动态特征矩阵,使用Fisher算法对所述动态特征矩阵中的N*y+i*x个特征进行打分,并按照分数从高到底的顺序排序,筛选分数高的前w个特征作为最具有代表性的特征,w根据实际情况确定;
所述数据训练器用于在支持向量机分类中,根据功能核磁共振图像对筛选出的前w个最具有代表性的特征进行训练得到一个分类器;
所述辅助诊断器用于根据训练得到的分类器进行阿尔茨海默病的辅助诊断。
一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置的使用方法,包括以下步骤:
步骤1:功能核磁共振图像的预处理;
步骤2:利用预处理后的图像构建动态脑功能网络;
步骤3:对动态脑功能网络计算节点度量,并提取每个动态脑功能网络的动态特征,然后将过滤后的动态特征的特征值拼接为一个动态特征矩阵;
步骤4:运用生成的动态特征矩阵对阿尔茨海默病进行辅助诊断。
所述的步骤1功能核磁共振图像的预处理包括以下步骤:
1.1)将获取的i个待测功能核磁共振图像通过所述时间片校正器进行时间片校正,得到i个时间片校正后的功能核磁共振图像(I-1,I-2,I-3,…,I-i),其中i表示选取的待测功能核磁共振图像的个数;
1.2)将所述i个时间校正后的功能核磁共振图像通过所述头动校正器进行头动校正,得到i个头动校正后的功能核磁共振图像(H-1,H-2,H-3,...,H-i);
1.3)将所述i个头动校正后的功能核磁共振图像通过所述空间标准化器进行空间标准化,得到i个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,...,F-i);
1.4)所述i个空间标准化后的功能核磁共振图像通过所述平滑降噪器进行平滑降噪,得到i个标准功能核磁共振图像(S-1,S-2,S-3,...,S-i)。
所述的步骤2利用预处理后的图像构建动态脑功能网络,包括以下步骤:
2.1)将i个标准功能核磁共振图像中的每个标准功能核磁共振图像通过所述模板匹配器与规格为m的标准匹配模板进行感兴趣区域的匹配,匹配后的每个功能核磁共振图像包含m个脑区,得到i个匹配后的功能核磁共振图像(A-1,A-2,A-3,...,A-i);
2.2)将i个匹配后的功能核磁共振图像通过所述时间序列分割器对每个匹配后的功能核磁共振图像中的每个脑区的时间序列进行分割,得到s段瞬时时间序列,每段瞬时时间序列代表一个脑区的瞬时信息(T-1,T-2,...,T-s),i个分割后的功能核磁共振图像表示为其中s的取值范围根据预设的分割间隔确定;
2.3)将分割得到的s段瞬时时间序列通过所述脑网络构建器构建出i个分割后的功能核磁共振图像的i*s个动态脑功能网络
所述的步骤3对动态脑功能网络计算节点度量,并提取每个动态脑功能网络的动态特征,然后将过滤后的动态特征的特征值拼接为一个动态特征矩阵,包括以下步骤:
3.1)将i个分割后的功能核磁共振图像的i*s个动态脑功能网络通过所述节点度量生成器生成i*s个动态脑功能网络的节点度量,并计算得到z个动态特征的特征值表示为所述z个动态特征包括i*s*x个全局特征和i*s*m*y个局部特征,即z=i*s*x+i*s*m*y,x表示每个动态脑功能网络计算的全局特征的个数,y表示每个动态脑功能网络计算的局部特征的个数;
3.2)将每组特征值通过时间序列生成器生成一个时间序列,则i*s个动态脑功能网络构成的z/s个时间序列表示为所述每组特征值包括全局的每组特征值和局部的每组特征值,所述全局的每组特征值包括每个分割后的功能核磁共振图像中的每个全局特征的s个动态脑功能网络的特征值,所述局部的每组特征值包括从每个分割后的功能核磁共振图像中的s*m*y个局部特征中的每个局部特征按照m个脑区分组得到的s个动态脑功能网络的特征值;
3.3)将每组特征生成的时间序列通过特征提取器中的基于小波的时间序列熵再次提取特征,得到z/s个特征的新特征值
3.4)将所述节点度量中的i*s*m*y个局部特征通过特征过滤器进行过滤,首先通过节点度量中的degree方法计算出每个匹配后的功能核磁共振图像中的每个脑区的degree,并计算出每个匹配后的功能核磁共振图像中的m个脑区degree的平均值然后计算出每个匹配后的功能核磁共振图像中的m个脑区的标准方差σ,其次将所述m个脑区中degree属于区间/>的脑区过滤为关键脑区,最后将i个匹配后的功能核磁共振图像过滤得到的N个脑区的N*y个局部特征生成的新特征值与i*x个全局特征生成的新特征值通过特征过滤器拼接为一个动态特征矩阵,/>其中nj表示第j个匹配后的功能核磁共振图像中的关键脑区的个数,i表示匹配后的功能核磁共振图像的个数;
所述步骤4运用生成的动态特征矩阵对阿尔茨海默病进行辅助诊断,包括以下步骤:
4.1)将特征过滤器提取的动态特征矩阵通过特征筛选器中的Fisher算法对所述动态特征矩阵中的N*y+i*x个特征进行打分,并按照分数从高到底的顺序排序,筛选分数高的前w个特征作为最具有代表性的特征,w根据实际情况确定;
4.2)将筛选出的前w个最具有代表性的特征通过数据训练器中的支持向量机进行训练得到分类器;
4.3)使用辅助诊断器通过得到的分类器进行阿尔茨海默病的辅助诊断。
本发明的有益效果是:
本发明是一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法,克服了以往的静态脑功能网络无法表示动态信息的缺陷,使得对大脑活动信息进行动态分析成为可能,令功能核磁共振图像的信号信息得到充分发挥,起到了更好的为医疗辅助诊断服务的效果。
附图说明
图1为本发明实施例中的基于动态脑功能网络的阿尔茨海默病辅助诊断装置结构框图。
图2为本发明实施例中的动态脑功能网络构建方法流程图。
图3为本发明实施例中的动态信息特征提取和过滤及辅助诊断的方法流程图。
具体实施方式
下面是结合附图对本发明的技术方案进行详细说明。
在fMRI数据中,不同时间下的信号包含了大脑活动信息,而现有的研究都是基于对复杂网络进行静态网络分析,忽略了大脑是一个动态结构,随着时间的变化,大脑里的神经元之间的连接是变化的,因而需要构建动态脑功能网络来对大脑活动信息进行更全面的描述,并计算每个动态脑功能网络的节点度量特征值的特征。因此,提出了基于动态脑功能网络的设计方法,既考虑了脑网络在多个时间段下的瞬时差异,又考虑到了脑网络的结构,并将其用于基于支持向量机进行阿尔茨海默病的辅助诊断。
如图1所示,基于动态脑功能网络的阿尔茨海默病辅助诊断装置包括fMRI数据预处理单元、构建动态脑功能网络单元、生成用于训练的特征单元以及svm分类辅助诊断单元,首先将获取的606个待测功能核磁共振图像通过所述fMRI数据预处理单元得到预处理后的606个标准功能核磁共振图像,然后将预处理后的606个标准功能核磁共振图像通过所述构建动态脑功能网络单元进行时间维度上的分割,构建出606个功能核磁共振图像的15150个动态脑功能网络,其次将构建出的15150个动态脑功能网络通过所述生成用于训练的特征单元横向的提取出每个动态脑功能网络的动态特征,最后将从15150个动态脑功能网络提取出的经过特征过滤器过滤后的动态特征拼接成一个动态特征矩阵,最后通过所述svm分类辅助诊断单元使用Fisher算法进行特征筛选,筛选出具有代表性的特征,并进行训练,用于辅助诊断。
所述fMRI数据预处理单元包括时间片校正器、头动校正器、空间标准化器、平滑降噪器,首先将获取的606个待测功能核磁共振图像通过所述时间片校正器进行时间片校正,得到606个时间片校正后的功能核磁共振图像,然后将所述606个时间校正后的功能核磁共振图像通过所述头动校正器进行头动校正,得到606个头动校正后的功能核磁共振图像,其次将所述606个头动校正后的功能核磁共振图像通过所述空间标准化器进行空间标准化,得到606个空间标准化后的功能核磁共振图像,最后将所述606个空间标准化后的功能核磁共振图像通过所述平滑降噪器进行平滑降噪,得到606个标准功能核磁共振图像;
所述时间片校正器用于将输入的606个待测功能核磁图像进行时间片校正,得到606个时间片校正后的功能核磁共振图像(I-1,I-2,I-3,...,I-606);
所述头动校正器用于将606个时间片校正后的功能核磁图像(I-1,I-2,I-3,...,I-606)进行头动校正,得到606个头动校正后的功能核磁共振图像(H-1,H-2,H-3,...,H-606);
所述空间标准化器用于将606个头动校正后的功能核磁图像(H-1,H-2,H-3,...,H-606)进行空间标准化,得到606个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,...,F-606);
所述平滑降噪器用于将606个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,...,F-606)进行平滑降噪,得到606个标准功能核磁共振图像(S-1,S-2,S-3,...,S-606)。
所述构建动态脑功能网络单元包括模板匹配器、时间序列分割器以及脑网络构建器,首先将606个标准功能核磁共振图像中的每个标准功能核磁共振图像通过所述模板匹配器与264个脑区的标准power-264模板进行感兴趣区域的匹配,则每个匹配后的功能核磁共振图像包含264个脑区,然后将匹配后的功能核磁共振图像通过所述时间序列分割器将每个脑区的时间序列进行分割,得到25段瞬时时间序列,最后将分割得到的25段瞬时时间序列通过所述脑网络构建器构建出606个分割后的功能核磁共振图像的15150个动态脑功能网络;
所述模板匹配器用于将每个标准功能核磁共振图像与264个脑区的标准power-264模板进行匹配,匹配后的每个功能核磁共振图像包含264个脑区,606个匹配后的功能核磁共振图像表示为(A-1,A-2,A-3,...,A-606);
所述时间序列分割器用于将每个匹配后的功能核磁共振图像中的264个脑区的时间序列进行分割,将每个脑区的时间序列分割为25段瞬时时间序列,每段瞬时时间序列代表一个脑区的瞬时信息(T-1,T-2,...,T-25),606个分割后的功能核磁共振图像表示为
所述脑网络构建器用于将每个分割后的功能核磁共振图像中的25段瞬时时间序列构建出25个动态脑功能网络,则606个分割后的功能核磁共振图像得到15150个动态脑功能网络
所述生成用于训练的特征单元包括节点度量生成器、时间序列生成器、特征提取器以及特征过滤器,首先将构建动态脑功能网络单元输出的606个分割后的功能核磁共振图像的15150个动态脑功能网络通过所述节点度量生成器计算15150个动态脑功能网络的节点度量,然后通过时间序列生成器将15150个动态脑功能网络中的每个节点度量构成一个时间序列,其次通过所述特征提取器为每个节点度量形成的时间序列提取新特征值,最后将提取到的所有的新特征值通过所述特征过滤器过滤后拼接为一个动态特征矩阵;
所述节点度量生成器用于生成15150个动态脑功能网络的节点度量,并计算得到28057800个动态特征的特征值表示为所述28057800个动态特征包括60600个全局特征和27997200个局部特征,每个动态脑功能网络计算的全局特征的个数为4个,每个动态脑功能网络计算的局部特征的个数为7个;
所述时间序列生成器用于将得到的每组特征值生成一个时间序列,则15150个动态脑功能网络构成的1122312个时间序列表示为所述每组特征值包括全局的每组特征值和局部的每组特征值,所述全局的每组特征值包括每个分割后的功能核磁共振图像中的每个全局特征的25个动态脑功能网络的特征值,所述局部的每组特征值包括从每个分割后的功能核磁共振图像中的46200个局部特征中的每个局部特征按照264个各脑区分组得到的25个动态脑功能网络的特征值;
所述特征提取器用于对每组特征生成的时间序列,基于小波的时间序列熵再次提取特征,得到1122312个特征的新特征值
所述特征过滤器用于对所述节点度量中的27997200个局部特征进行过滤,首先通过节点度量中的degree方法计算出每个匹配后的功能核磁共振图像中的每个脑区的degree(度),并计算出每个匹配后的功能核磁共振图像中的264个脑区degree的平均值然后计算出每个匹配后的功能核磁共振图像中的264个脑区的标准方差σ,其次将264个脑区中degree属于区间/>的脑区过滤为关键脑区,最后将606个匹配后的功能核磁共振图像过滤得到的8484个脑区的59388个局部特征生成的新特征值与2424个全局特征生成的新特征值通过特征过滤器拼接为一个动态特征矩阵。
所述的svm分类辅助诊断单元包括特征筛选器、数据训练器、辅助诊断器,首先将生成用于训练的特征单元输出的动态特征矩阵通过特征筛选器,使用Fisher算法对所述动态特征矩阵进行动态特征筛选,然后通过数据训练器进行分类训练,最终通过辅助诊断器实现对阿尔茨海默病的诊断;
所述特征筛选器用于将特征过滤器提取的动态特征矩阵,使用Fisher算法对所述动态特征矩阵中的61812个特征进行打分,并按照分数从高到底的顺序排序,筛选分数高的前80个特征作为最具有代表性的特征;
所述数据训练器用于在支持向量机分类中,根据功能核磁共振图像对筛选出的前80个最具有代表性的特征进行训练得到一个分类器;
所述辅助诊断器用于根据训练得到的分类器进行阿尔茨海默病的辅助诊断。
一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置的使用方法,包括以下步骤:
步骤1:功能核磁共振图像的预处理,包括以下步骤:
1.1)将获取的606个待测功能核磁共振图像通过所述时间片校正器进行时间片校正,得到606个时间片校正后的功能核磁共振图像(I-1,I-2,I-3,…,I-606);
1.2)将所述606个时间校正后的功能核磁共振图像通过所述头动校正器进行头动校正,得到606个头动校正后的功能核磁共振图像(H-1,H-2,H-3,…,H-606);
1.3)将所述606个头动校正后的功能核磁共振图像通过所述空间标准化器进行空间标准化,得到606个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,…,F-606);
1.4)所述606个空间标准化后的功能核磁共振图像通过所述平滑降噪器进行平滑降噪,得到606个标准功能核磁共振图像(S-1,S-2,S-3,...,S-606)。
步骤2:利用预处理后的图像构建动态脑功能网络,具体构建流程如图2所示,包括以下步骤:
2.1)将606个标准功能核磁共振图像中的每个标准功能核磁共振图像通过所述模板匹配器与264个脑区的标准power-264模板进行匹配,匹配后的每个功能核磁共振图像包含264个脑区,得到606个匹配后的功能核磁共振图像(A-1,A-2,A-3,...,A-606);
2.2)将606个匹配后的功能核磁共振图像通过所述时间序列分割器对每个匹配后的功能核磁共振图像中的每个脑区的时间序列进行分割,得到25段瞬时时间序列,每段瞬时时间序列代表一个脑区的瞬时信息(T-1,T-2,...,T-25),606个分割后的功能核磁共振图像表示为并将分割后的每个图像依次通过互信息计算器,对分割得到的图像的各个节点之间的匹配关系进行计算,确定每个动态脑功能网络中边的表示,采用互信息的方法来计算各个节点之间的相关关系,并将其构成关联矩阵;
2.3)将计算得到的互信息值通过所述脑网络构建器构建出606个分割后的功能核磁共振图像的15150个动态脑功能网络
计算出fMRI图像任意两个节点的互信息并得到关联矩阵后,设置互信息的阈值为0.18,将互信息值与设定的阈值0.18进行比较,当互信息值比设定的阈值0.18大时,将其边设为1,即两节点之间有关联,否则将其边设为0,即无关联,由此将关联矩阵转化为邻接矩阵,构建出动态脑功能网络。
本实施方式是对功能核磁共振图像进行计算机辅助诊断,帮助医生进行诊断,其中动态信息特征提取和过滤和筛选及辅助诊断的方法流程图如图3所示,具体步骤包括步骤3和步骤4两部分。
步骤3对动态脑功能网络计算节点度量,并提取每个动态脑功能网络的动态特征,然后将过滤后的动态特征的特征值拼接为一个动态特征矩阵,包括以下步骤:
3.1)将606个分割后的功能核磁共振图像的15150个动态脑功能网络通过所述节点度量生成器生成15150个动态脑功能网络的节点度量,并计算得到28057800个动态特征的特征值表示为所述28057800个动态特征包括60600个全局特征和27997200个局部特征,每个动态脑功能网络计算的全局特征的个数为4个,每个动态脑功能网络计算的局部特征的个数为7个;
4个全局特征分别为:clustering_coefficients(聚类系数)、characteristicpath length(特征路径长度)、Global efficiency(全局效率)、Transitivity(传递性);
7个局部特征分别为:Local efficiency(局部效率)、degree(度)、betweennesscentrality(中介中心性)、Pagerank centrality(pagerank中心性)、node strength(节点强度)、k-coreness centrality(k-核心中心性)、flow coefficient(流量系数)。
3.2)将每组特征值通过时间序列生成器生成一个时间序列,则15150个动态脑功能网络构成的1122312个时间序列表示为所述每组特征值包括全局的每组特征值和局部的每组特征值,所述全局的每组特征值包括每个分割后的功能核磁共振图像中的每个全局特征的25个动态脑功能网络的特征值,所述局部的每组特征值包括从每个分割后的功能核磁共振图像中的46200个局部特征中的每个局部特征按照264个脑区分组得到的25个动态脑功能网络的特征值;
对于一个功能核磁共振图像可以通过时间序列分割器再通过脑网络构建器,构建25个脑功能网络,每个构成的脑网络计算全局特征和局部特征,对于每种特征在25个脑功能网络下可以通过时间序列生成器生成具有25个时间点的时间序列。
3.3)将每组特征生成的时间序列通过特征提取器中的基于小波的时间序列熵再次提取特征,得到1122312个特征的新特征值
基于小波的时间序列熵再次提取特征。熵这个概念最早是用于热力学中,毕竟这个字是火字旁,用于衡量一个系统能量的不可用程度,熵越大,能量的不可用程度就越大;越小能量的不可用程度越低。它的物理意义是体系中混乱程度或者复杂程度的度量。关于熵的应用也在不断拓展,从热力学到生物学、物理学,以及在时间序列分析上都有应用。香农熵(信息熵)是一个数学上颇为抽象的概念,可以把信息熵理解成某种特定信息的出现概率(离散随机事件的出现概率)。一个系统越是有序,信息熵就越低;反之,一个系统越是混乱,信息熵就越高。信息熵也可以说是系统有序化程度的一个度量
3.4)将所述节点度量中的27997200个局部特征通过特征过滤器进行过滤,首先通过节点度量中的degree方法计算出每个匹配后的功能核磁共振图像中的每个脑区的degree(度),并计算出每个匹配后的功能核磁共振图像中的264个脑区degree的平均值然后计算出每个匹配后的功能核磁共振图像中的264个脑区的标准方差σ,其次将264个脑区中degree属于区间/>的脑区过滤为关键脑区,最后将606个匹配后的功能核磁共振图像过滤得到的8484个脑区的59388个局部特征生成的新特征值与2424个全局特征生成的新特征值通过特征过滤器拼接为一个动态特征矩阵。
步骤4运用生成的动态特征矩阵对阿尔茨海默病进行辅助诊断,包括以下步骤:
4.1)将特征过滤器提取的动态特征矩阵通过特征筛选器中的Fisher算法对所述动态特征矩阵中的61812个特征进行打分,并按照分数从高到底的顺序排序,筛选分数高的前80个特征作为最具有代表性的特征;
将所有功能核磁共振图像构建的动态脑功能网络,对于所有节点度量生成的时间序列提取的特征拼成一个矩阵,通过Fisher算法进行特征筛选,Fisher准则基本原理是找到一个最合适的投影轴,使两类样本在该轴上投影之间的距离尽可能远,而每一类样本的投影尽可能紧凑,从而使分类效果为最佳。
4.2)将筛选出的前80个最具有代表性的特征通过数据训练器中的支持向量机进行训练得到分类器;
4.3)使用辅助诊断器通过得到的分类器进行阿尔茨海默病的辅助诊断。
Claims (7)
1.一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置,其特征在于,包括fMRI数据预处理单元、构建动态脑功能网络单元、生成用于训练的特征单元以及svm分类辅助诊断单元,首先将获取的i个待测功能核磁共振图像通过所述fMRI数据预处理单元得到预处理后的的i个标准功能核磁共振图像,然后将预处理后的i个标准功能核磁共振图像通过所述构建动态脑功能网络单元进行时间维度上的分割,构建出i个功能核磁共振图像的i*s个动态脑功能网络,其次将构建出的i*s个动态脑功能网络通过所述生成用于训练的特征单元横向的提取出每个动态脑功能网络的动态特征,最后将从i*s个动态脑功能网络提取出的经过特征过滤器过滤后的动态特征拼接成一个矩阵,最后通过所述svm分类辅助诊断单元使用Fisher算法进行特征筛选,筛选出具有代表性的特征,并进行训练,用于辅助诊断;
所述的fMRI数据预处理单元包括时间片校正器、头动校正器、空间标准化器、平滑降噪器,首先将获取的i个待测功能核磁共振图像通过所述时间片校正器进行时间片校正,得到i个时间片校正后的功能核磁共振图像,然后将所述i个时间片校正后的功能核磁共振图像通过所述头动校正器进行头动校正,得到i个头动校正后的功能核磁共振图像,其次将所述i个头动校正后的功能核磁共振图像通过所述空间标准化器进行空间标准化,得到i个空间标准化后的功能核磁共振图像,最后将所述i个空间标准化后的功能核磁共振图像通过所述平滑降噪器进行平滑降噪,得到i个标准功能核磁共振图像;
所述时间片校正器用于将输入的i个待测功能核磁图像进行时间片校正,得到i个时间片校正后的功能核磁共振图像(I-1,I-2,I-3,...,I-i),其中i表示选取的待测功能核磁共振图像的个数;
所述头动校正器用于将i个时间片校正后的功能核磁图像(I-1,I-2,I-3,...,I-i)进行头动校正,得到i个头动校正后的功能核磁共振图像(H-1,H-2,H-3,...,H-i);
所述空间标准化器用于将i个头动校正后的功能核磁图像(H-1,H-2,H-3,...,H-i)进行空间标准化,得到i个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,...,F-i);
所述平滑降噪器用于将i个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,...,F-i)进行平滑降噪,得到i个标准功能核磁共振图像(S-1,S-2,S-3,...,S-i);
所述构建动态脑功能网络单元包括模板匹配器、时间序列分割器以及脑网络构建器,首先将i个标准功能核磁共振图像中的每个标准功能核磁共振图像通过所述模板匹配器与规格为m的标准匹配模板进行感兴趣区域的匹配,则每个匹配后的功能核磁共振图像包含m个脑区,然后将匹配后的功能核磁共振图像通过所述时间序列分割器将每个脑区的时间序列进行分割,得到s段瞬时时间序列,最后将分割得到的s段瞬时时间序列通过所述脑网络构建器构建出i个分割后的功能核磁共振图像的i*s个动态脑功能网络;
所述模板匹配器用于将每个标准功能核磁共振图像与规格为m的标准匹配模板进行匹配,匹配后的每个功能核磁共振图像包含m个脑区,i个匹配后的功能核磁共振图像表示为(A-1,A-2,A-3,...,A-i);
所述时间序列分割器用于将每个匹配后的功能核磁共振图像中的m个脑区的时间序列进行分割,将每个脑区的时间序列分割为s段瞬时时间序列,每段瞬时时间序列代表一个脑区的瞬时信息(T-1,T-2,...,T-s),i个分割后的功能核磁共振图像表示为其中s的取值范围根据预设的分割间隔确定;
所述脑网络构建器用于将每个分割后的功能核磁共振图像中的s段瞬时时间序列构建出s个动态脑功能网络,则i个分割后的功能核磁共振图像得到i*s个动态脑功能网络
所述生成用于训练的特征单元包括节点度量生成器、时间序列生成器、特征提取器以及特征过滤器,首先将构建动态脑功能网络单元输出的i个分割后的功能核磁共振图像的i*s个动态脑功能网络通过所述节点度量生成器计算i*s个动态脑功能网络的节点度量,然后通过时间序列生成器将i*s个动态脑功能网络中的每个节点度量构成一个时间序列,其次通过所述特征提取器为每个节点度量形成的时间序列提取新特征值,最后将提取到的所有的新特征值通过所述特征过滤器过滤后拼接为一个动态特征矩阵;
所述节点度量生成器用于生成i*s个动态脑功能网络的节点度量,并计算得到z个动态特征的特征值表示为所述z个动态特征包括i*s*x个全局特征和i*s*m*y个局部特征,即z=i*s*x+i*s*m*y,x表示每个动态脑功能网络计算的全局特征的个数,y表示每个动态脑功能网络计算的局部特征的个数;
所述时间序列生成器用于将得到的每组特征值生成一个时间序列,则i*s个动态脑功能网络构成的z/s个时间序列表示为所述每组特征值包括全局的每组特征值和局部的每组特征值,所述全局的每组特征值包括每个分割后的功能核磁共振图像中的每个全局特征的s个动态脑功能网络的特征值,所述局部的每组特征值包括从每个分割后的功能核磁共振图像中的s*m*y个局部特征中的每个局部特征按照m个脑区分组得到的s个动态脑功能网络的特征值;
所述特征提取器用于对每组特征生成的时间序列,基于小波的时间序列熵再次提取特征,得到z/s个特征的新特征值
所述特征过滤器用于对所述节点度量中的i*s*m*y个局部特征进行过滤,首先通过节点度量中的degree方法计算出每个匹配后的功能核磁共振图像中的每个脑区的degree,并计算出每个匹配后的功能核磁共振图像中的m个脑区degree的平均值然后计算出每个匹配后的功能核磁共振图像中的m个脑区的标准方差σ,其次将所述m个脑区中degree属于区间/>的脑区过滤为关键脑区,最后将i个匹配后的功能核磁共振图像过滤得到的N个脑区的N*y个局部特征生成的新特征值与i*x个全局特征生成的新特征值通过特征过滤器拼接为一个动态特征矩阵,/>其中nj表示第j个匹配后的功能核磁共振图像中的关键脑区的个数,i表示匹配后的功能核磁共振图像的个数。
2.根据权利要求1所述的一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置,其特征在于,所述的svm分类辅助诊断单元包括特征筛选器、数据训练器、辅助诊断器,首先将生成用于训练的特征单元输出的动态特征矩阵通过特征筛选器,使用Fisher算法对所述动态特征矩阵进行动态特征筛选,然后通过数据训练器进行分类训练,最终通过辅助诊断器实现对阿尔茨海默病的诊断;
所述特征筛选器用于将特征过滤器提取的动态特征矩阵,使用Fisher算法对所述动态特征矩阵中的N*y+i*x个特征进行打分,并按照分数从高到底的顺序排序,筛选分数高的前w个特征作为最具有代表性的特征,w根据实际情况确定;
所述数据训练器用于在支持向量机分类中,根据功能核磁共振图像对筛选出的前w个最具有代表性的特征进行训练得到一个分类器;
所述辅助诊断器用于根据训练得到的分类器进行阿尔茨海默病的辅助诊断。
3.根据权利要求1-2任意一项所述的一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置的使用方法,其特征在于,包括以下步骤:
步骤1:功能核磁共振图像的预处理;
步骤2:利用预处理后的图像构建动态脑功能网络;
步骤3:对动态脑功能网络计算节点度量,并提取每个动态脑功能网络的动态特征,然后将过滤后的动态特征的特征值拼接为一个动态特征矩阵;
步骤4:运用生成的动态特征矩阵对阿尔茨海默病进行辅助诊断。
4.根据权利要求3所述的一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置的使用方法,其特征在于,所述的步骤1功能核磁共振图像的预处理包括以下步骤:
1.1)将获取的i个待测功能核磁共振图像通过所述时间片校正器进行时间片校正,得到i个时间片校正后的功能核磁共振图像(I-1,I-2,I-3,...,I-i),其中i表示选取的待测功能核磁共振图像的个数;
1.2)将所述i个时间校正后的功能核磁共振图像通过所述头动校正器进行头动校正,得到i个头动校正后的功能核磁共振图像(H-1,H-2,H-3,...,H-i);
1.3)将所述i个头动校正后的功能核磁共振图像通过所述空间标准化器进行空间标准化,得到i个空间标准化后的功能核磁共振图像(F-1,F-2,F-3,...,F-i);
1.4)所述i个空间标准化后的功能核磁共振图像通过所述平滑降噪器进行平滑降噪,得到i个标准功能核磁共振图像(S-1,S-2,S-3,...,S-i)。
5.根据权利要求3所述的一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置的使用方法,其特征在于,所述的步骤2利用预处理后的图像构建动态脑功能网络,包括以下步骤:
2.1)将i个标准功能核磁共振图像中的每个标准功能核磁共振图像通过所述模板匹配器与规格为m的标准匹配模板进行感兴趣区域的匹配,匹配后的每个功能核磁共振图像包含m个脑区,得到i个匹配后的功能核磁共振图像(A-1,A-2,A-3,...,A-i);
2.2)将i个匹配后的功能核磁共振图像通过所述时间序列分割器对每个匹配后的功能核磁共振图像中的每个脑区的时间序列进行分割,得到s段瞬时时间序列,每段瞬时时间序列代表一个脑区的瞬时信息(T-1,T-2,...,T-s),i个分割后的功能核磁共振图像表示为其中s的取值范围根据预设的分割间隔确定;
2.3)将分割得到的s段瞬时时间序列通过所述脑网络构建器构建出i个分割后的功能核磁共振图像的i*s个动态脑功能网络
6.根据权利要求3所述的一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置的使用方法,其特征在于,所述的步骤3对动态脑功能网络计算节点度量,并提取每个动态脑功能网络的动态特征,然后将过滤后的动态特征的特征值拼接为一个动态特征矩阵,包括以下步骤:
3.1)将i个分割后的功能核磁共振图像的i*s个动态脑功能网络通过所述节点度量生成器生成i*s个动态脑功能网络的节点度量,并计算得到z个动态特征的特征值表示为所述z个动态特征包括i*s*x个全局特征和i*s*m*y个局部特征,即z=i*s*x+i*s*m*y,x表示每个动态脑功能网络计算的全局特征的个数,y表示每个动态脑功能网络计算的局部特征的个数;
3.2)将每组特征值通过时间序列生成器生成一个时间序列,则i*s个动态脑功能网络构成的z/s个时间序列表示为所述每组特征值包括全局的每组特征值和局部的每组特征值,所述全局的每组特征值包括每个分割后的功能核磁共振图像中的每个全局特征的s个动态脑功能网络的特征值,所述局部的每组特征值包括从每个分割后的功能核磁共振图像中的s*m*y个局部特征中的每个局部特征按照m个脑区分组得到的s个动态脑功能网络的特征值;
3.3)将每组特征生成的时间序列通过特征提取器中的基于小波的时间序列熵再次提取特征,得到z/s个特征的新特征值
3.4)将所述节点度量中的i*s*m*y个局部特征通过特征过滤器进行过滤,首先通过节点度量中的degree方法计算出每个匹配后的功能核磁共振图像中的每个脑区的degree,并计算出每个匹配后的功能核磁共振图像中的m个脑区degree的平均值然后计算出每个匹配后的功能核磁共振图像中的m个脑区的标准方差σ,其次将所述m个脑区中degree属于区间/>的脑区过滤为关键脑区,最后将i个匹配后的功能核磁共振图像过滤得到的N个脑区的N*y个局部特征生成的新特征值与i*x个全局特征生成的新特征值通过特征过滤器拼接为一个动态特征矩阵,/>其中nj表示第j个匹配后的功能核磁共振图像中的关键脑区的个数,i表示匹配后的功能核磁共振图像的个数。
7.根据权利要求3所述的一种基于动态脑功能网络的阿尔茨海默病辅助诊断装置的使用方法,其特征在于,所述步骤4运用生成的动态特征矩阵对阿尔茨海默病进行辅助诊断,包括以下步骤:
4.1)将特征过滤器提取的动态特征矩阵通过特征筛选器中的Fisher算法对所述动态特征矩阵中的N*y+i*x个特征进行打分,并按照分数从高到底的顺序排序,筛选分数高的前w个特征作为最具有代表性的特征,w根据实际情况确定;
4.2)将筛选出的前w个最具有代表性的特征通过数据训练器中的支持向量机进行训练得到分类器;
4.3)使用辅助诊断器通过得到的分类器进行阿尔茨海默病的辅助诊断。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910916563.8A CN110718301B (zh) | 2019-09-26 | 2019-09-26 | 基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910916563.8A CN110718301B (zh) | 2019-09-26 | 2019-09-26 | 基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110718301A CN110718301A (zh) | 2020-01-21 |
CN110718301B true CN110718301B (zh) | 2024-04-19 |
Family
ID=69210922
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910916563.8A Active CN110718301B (zh) | 2019-09-26 | 2019-09-26 | 基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110718301B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113177633B (zh) * | 2021-04-20 | 2023-04-25 | 浙江大学 | 一种深度解耦时间序列预测方法 |
CN112991335B (zh) * | 2021-04-23 | 2024-04-26 | 合肥中聚源智能科技有限公司 | 一种想像能力评估方法、系统及可读存储介质 |
CN116541751B (zh) * | 2023-07-03 | 2023-09-12 | 中国医学科学院生物医学工程研究所 | 一种基于脑功能连接网络特征的脑电信号分类方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104715261A (zh) * | 2015-03-23 | 2015-06-17 | 南京工业大学 | fMRI动态脑功能子网络构建及并联SVM加权识别方法 |
CN107256408A (zh) * | 2017-05-12 | 2017-10-17 | 常州大学 | 一种大脑功能网络的关键路径搜索方法 |
CN109859839A (zh) * | 2019-01-23 | 2019-06-07 | 桂林电子科技大学 | 一种基于分层集成学习的阿尔茨海默症分类方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110301431A1 (en) * | 2010-06-05 | 2011-12-08 | The Board Of Trustees Of The Leland Stanford Junior University | Methods of classifying cognitive states and traits and applications thereof |
US9510756B2 (en) * | 2012-03-05 | 2016-12-06 | Siemens Healthcare Gmbh | Method and system for diagnosis of attention deficit hyperactivity disorder from magnetic resonance images |
-
2019
- 2019-09-26 CN CN201910916563.8A patent/CN110718301B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104715261A (zh) * | 2015-03-23 | 2015-06-17 | 南京工业大学 | fMRI动态脑功能子网络构建及并联SVM加权识别方法 |
CN107256408A (zh) * | 2017-05-12 | 2017-10-17 | 常州大学 | 一种大脑功能网络的关键路径搜索方法 |
CN109859839A (zh) * | 2019-01-23 | 2019-06-07 | 桂林电子科技大学 | 一种基于分层集成学习的阿尔茨海默症分类方法 |
Non-Patent Citations (3)
Title |
---|
fMRI动态功能网络构建及其在脑部疾病识别中的应用;马士林;梅雪;李微微;周宇;;计算机科学(第10期);全文 * |
基于fMRI动态功能连接的抑郁症患者分类研究;皇甫浩然;杨剑;杨阳;;计算机应用研究(第03期);全文 * |
基于动态图论特征的阿尔茨海默病早期预测;董国昭;曾晓天;杨柳;张逸鹤;张勇;唐晓英;;北京生物医学工程(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110718301A (zh) | 2020-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110718301B (zh) | 基于动态脑功能网络的阿尔茨海默病辅助诊断装置及方法 | |
CN107133651B (zh) | 基于超网络判别子图的功能磁共振影像数据分类方法 | |
CN109528197B (zh) | 基于脑功能图谱进行精神疾病的个体化预测方法和系统 | |
CN113724880A (zh) | 一种异常脑连接预测系统、方法、装置及可读存储介质 | |
CN110840468B (zh) | 自闭症风险评估方法及装置、终端设备、存储介质 | |
CN109948647A (zh) | 一种基于深度残差网络的心电图分类方法及系统 | |
CN111783887B (zh) | 基于fMRI小世界脑网络计算机分类测谎识别方法 | |
CN103646183A (zh) | 一种基于人工神经网络和多模态mri的阿尔茨海默病智能判别分析方法 | |
CN108122221A (zh) | 弥散加权成像图像中脑缺血区域的分割方法及装置 | |
US8271414B2 (en) | Network characterization, feature extraction and application to classification | |
CN115272295A (zh) | 基于时域-空域联合状态的动态脑功能网络分析方法及系统 | |
CN113069117A (zh) | 一种基于时间卷积神经网络的脑电情绪识别方法及系统 | |
CN110502989A (zh) | 一种小样本高光谱人脸识别方法及系统 | |
CN114926396B (zh) | 一种精神障碍类磁共振图像初步筛查模型构建方法 | |
CN110889501B (zh) | 一种基于非平稳动态贝叶斯网的人脑效应连接网络构建方法 | |
CN108805181B (zh) | 一种基于多分类模型的图像分类装置及分类方法 | |
CN112336369B (zh) | 一种多通道心音信号的冠心病风险指数评估系统 | |
Mousavian et al. | Depression detection using atlas from fMRI images | |
CN112861881A (zh) | 一种基于改进MobileNet模型的蜂窝肺识别方法 | |
CN116543154A (zh) | 一种基于多层次语义特征的医学影像分割方法 | |
Duran et al. | Dual-HINet: dual hierarchical integration network of multigraphs for connectional brain template learning | |
CN108846407B (zh) | 基于独立成分高序不确定脑网络的磁共振影像分类方法 | |
CN107194918B (zh) | 数据分析方法及装置 | |
CN114224354B (zh) | 心律失常分类方法、装置及可读存储介质 | |
CN112052882B (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 |