CN109558906A - 基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法 - Google Patents
基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法 Download PDFInfo
- Publication number
- CN109558906A CN109558906A CN201811412834.8A CN201811412834A CN109558906A CN 109558906 A CN109558906 A CN 109558906A CN 201811412834 A CN201811412834 A CN 201811412834A CN 109558906 A CN109558906 A CN 109558906A
- Authority
- CN
- China
- Prior art keywords
- cluster centre
- time series
- voxel
- value
- dimensionality reduction
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,本发明利用高斯过程隐变量分析,时间序列间的绝对值距离,功能磁共振实验设计模式相关分析和距离指标建立了一种功能磁共振大数据聚类成像模型。利用高斯过程隐变量分析对数据协方差矩阵进行降维,利用降维后的数据向量作为聚类中心,然后利用与聚类中心的最小绝对值距离来对磁共振数据集对应体素进行归类,同时结合实验设计模式相关分析识别任务相关激活区并构建激活指标。本发明通过对协方差矩阵降维后较小的数据集确定聚类中心,避免了对磁共振大数据集直接聚类的常规做法,提高了成像效率、节省了计算资源,利用此模型架构实现功能磁共振成像目前是一种新的技术尝试。
Description
技术领域
本发明涉及一种医学成像方法,尤其涉及一种基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法。
背景技术
使用血氧水平依赖对比的功能磁共振成像提供了一种大脑对任务或刺激响应的氧合血流量的测量方法,并且在解读大脑功能上,功能磁共振成像是一种重要的、非侵入性的、高时间和空间分辨率的技术。在典型的功能磁共振成像实验中,几秒间隔出现的外部刺激将引起体素信号强度的变化、血液动力学响应的延迟。通常,这些功能磁共振数据集能够使用两类基于体素的分析方法进行处理:基于估计的血液动力学响应函数模型的统计方法和数据驱动的方法。
常用的几种数据驱动聚类分析技术中,包括k中心分析,模糊聚类分析和层次聚类分析。这些聚类方法是重要的探索性数据分析工具,并且已经被直接用于功能磁共振数据集去直接识别具有相似脑激活模式的区域。它们的主要缺点是它们的可靠性只能通过重复运行来确定。这就是所谓的不平衡数据问题,即激活区域是大脑的一小部分,它们嵌入大量非活动体素中。通过聚类分析本身处理这样的问题并不容易。此外,聚类族的生理解释并不明确。fMRI数据集是四维的,可能有64×64×25个以上的体素,具有可能大于100个以上的时间点,若在功能磁共振数据集上采用上述聚类方法直接进行聚类,这意味着它们需要巨大的存储器资源并且需要大量的计算机负载。它无法在个人计算机上实现大型数据集聚类,并且会出现运算缓慢、内存不足。
因此,为避免上述问题本发明并不直接在功能磁共振数据上使用聚类算法进行聚类,而是利用高斯过程隐变量分析对数据的协方差矩阵进行降维,利用降维后的数据向量作为聚类中心,然后利用与这些聚类中心的欧式距离最小原则对原输入功能磁共振大数据集进行归类,然后利用相关分析甄别筛选任务相关脑激活体素类代表点、识别激活区和构建激活指标定量激活强度。
发明内容
本发明的目的就在于提供一种解决上述问题,利用高斯过程隐变量分析对数据的协方差矩阵进行降维,利用降维后的数据向量作为聚类中心,然后利用与这些聚类中心对数据集X中的体素进行归类,然后利用相关分析甄别筛选任务相关脑激活体素类代表点、识别激活区和构建激活指标定量激活强度的基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法。
为了实现上述目的,本发明采用的技术方案是这样的:一种基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,包括以下步骤:
一种基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,包括以下步骤:
一种基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,其特征在于:包括以下步骤:
(1)采集脑功能磁共振数据,并对脑功能磁共振的功能像数据进行预处理,得到预处理后的数据集X,其中,I为数据集X中体素个数,xi为第i个体素的时间序列,每个体素的时间序列长度为N,I>>N;
(2)对数据集X的Z分数的协方差矩阵的主成分特征值进行本质维数估计,确定最优维数M,利用最优维数对协方差矩阵进行高斯过程隐变量分析降维,得到降维矩阵Y,yk为长度为N的时间序列,M<<N;
(3)分别以降维矩阵Y的M个行向量为聚类中心,对数据集X进行聚类;
(31)以M个行向量为聚类中心,第k个聚类中心yk对应第k个分类,k=1,2,…,M;
(32)计算xi与矩阵Y中M个聚类中心时间序列的绝对值距离,找出最小绝对值距离对应的聚类中心,将xi归入该聚类中心对应的分类中,得到聚类结果,其中,i=1,…,I;
(4)甄别筛选任务相关脑激活体素的聚类中心;
(41)构造一实验设计模式的时间序列,分别计算该时间序列与M个聚类中心时间序列的相关系数值;
(42)将对应相关系数值最大的聚类中心作为与实验任务相关的脑激活体素聚类中心;
(5)识别激活区,构建激活指标量化激活强度;
(51)将归入脑激活体素聚类中心的体素作为激活点,激活点所在的区域为激活区;
(52)利用激活区各体素时间序列分别与它们的聚类中心时间序列的绝对值距离值计算激活区各体素的反向标准化值,并将反向标准化值做为激活强度值;
(6)将步骤(5)识别出的脑激活区叠加投射在结构像模板上显示成像。
作为优选:所述步骤(1)中,预处理为:将脑功能磁共振的功能像数据,先进行头动校正、准化到EPI模板、空间平滑、再滤除信号的低频噪声。
作为优选:步骤(32)中,数据集X中各体素对应的时间序列xi与聚类中心时间序列yk之间的绝对值距离d采用下式计算:
其中,n表示时间点,N表示时间序列长度;xi(n)表示第i个体素的时间序列xi在第n个时间点的值,yk(n)表示第k个聚类中心的时间序列yk在第n个时间点的值。
作为优选:步骤(41)中具体为:
构造实验设计模式的时间序列,所述该时间序列包括任务态和静息态,在序列的每个时间点上,若为任务态对应的值设置为1,若为静息态设置值为0;
设实验设计模式时间序列表示为s(n),n=1,2,…,N,本发明采用下式计算实验设计模式时间序列与各聚类中心时间序列的相关系数值ρ:
n和l表示时间点。
作为优选:步骤(52)采用下式计算各体素的反向标准化值
D为激活区各体素分别与它们的聚类中心时间序列的绝对值距离的集合,max(D)表示集合D中的最大值,min(D)表示集合D中的最小值,di表示激活区第i个体素点对应的绝对值距离,Vi为激活区第i个体素点的的反向标准化值。
本发明的整体思路为:
将预处理后的脑功能磁共振数据先利用高斯过程隐变量分析对数据的协方差矩阵进行降维,降维前确定的要降维到的最优维数是通过对协方差矩阵的主成分特征值进行本质维数估计得到的;
再利用降维后得到的矩阵的行向量作为聚类中心,然后利用与聚类中心的最小绝对值距离来对磁共振数据集对应体素进行归类;
再构建实验设计模式时间序列并利用相关分析甄别筛选任务相关脑激活体素类代表点、识别激活区和构建激活指标定量激活强度。本发明的整体成像的技术架构如图1所示。
与现有技术相比,本发明的优点在于:利用高斯过程隐变量分析,时间序列间的绝对值距离,功能磁共振实验设计模式相关分析和距离指标建立了一种功能磁共振大数据聚类成像模型。通过对协方差矩阵降维后较小的数据集确定聚类中心,避免了对磁共振大数据集直接聚类的常规做法,提高了成像效率、节省了计算资源,利用此模型架构实现功能磁共振成像目前是一种新的技术尝试。具体说来,还体现在以下几点:
(1)该技术用于确定聚类中心的数据集与用于聚类成像的数据集采用分离式不同数据矩阵,不同于常规使用的聚类技术,如k中心聚类分析,模糊聚类分析、仿射聚类和层次聚类等等,这些常规技术直接在功能磁共振成像数据上使用聚类算法确定聚类中心进行聚类。发明技术直接减少了处理数据量,避免在个人计算机上出现内存溢出。
(2)本发明针对协方差矩阵高斯过程隐变量分析降维,不直接对输入磁共振数据直接进行降维操作,由于协方差矩阵数据量远小于输入磁共振数据量,提高了降维处理速度、进一步减少了降维操作对计算资源的占用,该方法适合功能磁共振大数据成像数据处理。
(3)本发明构造先验任务实验模式时间序列,并利用来对降维后的数据确定脑激活代表点。该技术较直接聚类确定活动代表点计算量更小,更简单。
因此,基于高斯过程隐变量降维类中心的脑功能磁共振聚类成像技术,提供了一种功能磁共振成像优化模型,在功能磁共振大数据成像领域具有很大的潜在应用前景。
附图说明
图1为本发明流程图;
图2为采用本发明技术探测到视觉与双手指运动同时刺激的脑功能活动区成像范例图。
具体实施方式
下面将结合附图对本发明作进一步说明。
实施例1:参见图1,一种基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,包括以下步骤:
(1)采集脑功能磁共振数据,并对脑功能磁共振的功能像数据进行预处理,得到预处理后的数据集X,其中,I为数据集X中体素个数,xi为第i个体素的时间序列,每个体素的时间序列长度为N,I>>N;
其中预处理为:将脑功能磁共振的功能像数据,先进行头动校正、准化到EPI模板、空间平滑、再滤除信号的低频噪声;
(2)对数据集X的Z分数的协方差矩阵的主成分特征值进行本质维数估计,确定最优维数M,利用最优维数对协方差矩阵进行高斯过程隐变量分析降维,得到降维矩阵Y,yk为长度为N的时间序列,M<<N;
(3)分别以降维矩阵Y的M个行向量为聚类中心,对数据集X进行聚类;
(31)以M个行向量为聚类中心,第k个聚类中心yk对应第k个分类,k=1,2,…,M;
(32)计算xi与矩阵Y中M个聚类中心时间序列的绝对值距离,找出最小绝对值距离对应的聚类中心,将xi归入该聚类中心对应的分类中,得到聚类结果,其中,i=1,…,I;
步骤(32)中,数据集X中各体素对应的时间序列xi与聚类中心时间序列yk之间的绝对值距离d采用下式计算:
其中,n表示时间点,N表示时间序列长度;xi(n)表示第i个体素的时间序列xi在第n个时间点的值,yk(n)表示第k个聚类中心的时间序列yk在第n个时间点的值;
(4)甄别筛选任务相关脑激活体素的聚类中心;
(41)构造一实验设计模式的时间序列,分别计算该时间序列与M个聚类中心时间序列的相关系数值,具体为:
构造实验设计模式的时间序列,所述该时间序列包括任务态和静息态,在序列的每个时间点上,若为任务态对应的值设置为1,若为静息态设置值为0;
设实验设计模式时间序列表示为s(n),n=1,2,…,N,本发明采用下式计算实验设计模式时间序列与各聚类中心时间序列的相关系数值ρ:
n和l表示时间点;
(42)将对应相关系数值最大的聚类中心作为与实验任务相关的脑激活体素聚类中心;
(5)识别激活区,构建激活指标量化激活强度;
(51)将归入脑激活体素聚类中心的体素作为激活点,激活点所在的区域为激活区;
(52)利用激活区各体素时间序列分别与它们的聚类中心时间序列的绝对值距离值计算激活区各体素的反向标准化值,并将反向标准化值做为激活强度值,具体是:采用下式计算各体素的反向标准化值
D为激活区各体素分别与它们的聚类中心时间序列的绝对值距离的集合,max(D)表示集合D中的最大值,min(D)表示集合D中的最小值,di表示激活区第i个体素点对应的绝对值距离,Vi为激活区第i个体素点的的反向标准化值。
(6)将步骤(5)识别出的脑激活区叠加投射在结构像模板上显示成像。
实施例2:
(1)采样得到脑功能磁共振数据,将脑功能磁共振的功能项数据先进行头动校正、准化到EPI模板、空间平滑、再滤除信号的低频噪声,得到预处理后的数据集X,此时,数据集X为I行N列矩阵;
(2)然后进行以下几步:
(21)将X标准化为Z分数,并对Z分数矩阵计算协方差矩阵C,此时,协方差矩阵C为N行N列,数据量比矩阵X(I行N列)大大减小;
(22)对协方差矩阵C的主成分特征值进行本质维数估计,确定最优维数M,利用最优维数对协方差矩阵进行高斯过程隐变量分析降维,得到降维矩阵Y,降维矩阵Y为M行N列的矩阵,N依然表示时间序列长度,该步骤进一步降低数据量;
(3)分别以降维矩阵Y的M个行向量为聚类中心,对数据集X进行聚类,主要包括步骤;
(31)以M个行向量为聚类中心,分别对应M个类别,分别为第一类、第二类、第三类、……、第M类;
(32)计算xi与M个聚类中心的绝对值距离,找出最小绝对值距离对应的聚类中心,假设最小绝对值距离对应的是第二个行向量y2,则将其归入第二类中,依次类推,把所有xi都进行归类,分别归入第一类到第M类中;
(4)甄别筛选任务相关脑激活体素的聚类中心;
(41)构造一实验设计模式的时间序列,分别计算该时间序列与M个聚类中心时间序列的相关系数值ρ;
(42)将对应最大相关系数值的聚类中心点作为与实验任务相关的脑激活体素聚类中心。例如:假设第二个聚类中心计算出来的ρ值最大,则将第二个聚类中心作为与实验任务相关的脑激活体素聚类中心;
(5)识别激活区,构建激活指标量化激活强度;
(51)将归入脑激活聚类中心的体素作为激活点;
(52)将(51)中作为激活点的体素时间序列分别与它们的聚类中心时间序列的绝对值距离值利用来计算激活区各体素的反向标准化值,并做为探测活动强度值。
(6)将步骤(5)识别出的脑激活区叠加投射在结构像模板上显示成像。例如图2展示了视觉与双手指运动同时刺激的实验设计的功能磁共振成像数据利用本发明方法探测脑激活区成像结果,图中右边的色条代表激活强度值大小,左边的脑区域越偏白色它的激活强度越大、值越接近1。图片的第一行展示了视觉活动刺激对相关脑区域的刺激激活,第二行代表双手指运动对相关运动皮层等相关脑区域的激活。
Claims (5)
1.一种基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,其特征在于:包括以下步骤:
(1)采集脑功能磁共振数据,并对脑功能磁共振的功能像数据进行预处理,得到预处理后的数据集X,其中,I为数据集X中体素个数,xi为第i个体素的时间序列,每个体素的时间序列长度为N,I>>N;
(2)对数据集X的Z分数的协方差矩阵的主成分特征值进行本质维数估计,确定最优维数M,利用最优维数对协方差矩阵进行高斯过程隐变量分析降维,得到降维矩阵Y,yk为长度为N的时间序列,M<<N;
(3)分别以降维矩阵Y的M个行向量为聚类中心,对数据集X进行聚类;
(31)以M个行向量为聚类中心,第k个聚类中心yk对应第k个分类,k=1,2,…,M;
(32)计算xi与矩阵Y中M个聚类中心时间序列的绝对值距离,找出最小绝对值距离对应的聚类中心,将xi归入该聚类中心对应的分类中,得到聚类结果,其中,i=1,…,I;
(4)甄别筛选任务相关脑激活体素的聚类中心;
(41)构造一实验设计模式的时间序列,分别计算该时间序列与M个聚类中心时间序列的相关系数值;
(42)将对应相关系数值最大的聚类中心作为与实验任务相关的脑激活体素聚类中心;
(5)识别激活区,构建激活指标量化激活强度;
(51)将归入脑激活体素聚类中心的体素作为激活点,激活点所在的区域为激活区;
(52)利用激活区各体素时间序列分别与它们的聚类中心时间序列的绝对值距离值计算激活区各体素的反向标准化值,并将反向标准化值做为激活强度值;
(6)将步骤(5)识别出的脑激活区叠加投射在结构像模板上显示成像。
2.根据权利要求1所述的基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,其特征在于:所述步骤(1)中,预处理为:将脑功能磁共振的功能像数据,先进行头动校正、准化到EPI模板、空间平滑、再滤除信号的低频噪声。
3.根据权利要求1所述的基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,其特征在于:步骤(32)中,数据集X中各体素对应的时间序列xi与聚类中心时间序列yk之间的绝对值距离d采用下式计算:
其中,n表示时间点,N表示时间序列长度;xi(n)表示第i个体素的时间序列xi在第n个时间点的值,yk(n)表示第k个聚类中心的时间序列yk在第n个时间点的值。
4.根据权利要求1所述的基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,其特征在于:步骤(41)中具体为:
构造实验设计模式的时间序列,所述该时间序列包括任务态和静息态,在序列的每个时间点上,若为任务态对应的值设置为1,若为静息态设置值为0;
设实验设计模式时间序列表示为s(n),n=1,2,…,N,本发明采用下式计算实验设计模式时间序列与各聚类中心时间序列的相关系数值ρ:
n和l表示时间点。
5.根据权利要求1所述的基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法,其特征在于:步骤(52)采用下式计算各体素的反向标准化值
D为激活区各体素分别与它们的聚类中心时间序列的绝对值距离的集合,max(D)表示集合D中的最大值,min(D)表示集合D中的最小值,di表示激活区第i个体素点对应的绝对值距离,Vi为激活区第i个体素点的的反向标准化值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811412834.8A CN109558906B (zh) | 2018-11-23 | 2018-11-23 | 基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811412834.8A CN109558906B (zh) | 2018-11-23 | 2018-11-23 | 基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109558906A true CN109558906A (zh) | 2019-04-02 |
CN109558906B CN109558906B (zh) | 2022-05-03 |
Family
ID=65867175
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811412834.8A Active CN109558906B (zh) | 2018-11-23 | 2018-11-23 | 基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109558906B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113657527A (zh) * | 2021-08-23 | 2021-11-16 | 广东工业大学 | 一种基于贝叶斯非参模型的脑功能区划分方法 |
CN114854920A (zh) * | 2022-05-06 | 2022-08-05 | 浙江大学 | 高斯混合模型嵌入的gru自编码器高炉异常监测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5887074A (en) * | 1996-12-13 | 1999-03-23 | Siemens Corporate Research, Inc. | Local principal component based method for detecting activation signals in functional MR images |
WO2009058915A1 (en) * | 2007-10-29 | 2009-05-07 | The Trustees Of The University Of Pennsylvania | Computer assisted diagnosis (cad) of cancer using multi-functional, multi-modal in-vivo magnetic resonance spectroscopy (mrs) and imaging (mri) |
CN101706561A (zh) * | 2009-11-20 | 2010-05-12 | 电子科技大学 | 功能性磁共振图像的聚类方法 |
CN104700120A (zh) * | 2015-03-23 | 2015-06-10 | 南京工业大学 | 一种基于自适应熵投影聚类算法的fMRI特征提取及分类方法 |
CN106204545A (zh) * | 2016-06-28 | 2016-12-07 | 南方医科大学 | 一种基于区域划分和Fisher向量的医学病灶图像特征表达方法 |
CN107392907A (zh) * | 2017-09-01 | 2017-11-24 | 上海理工大学 | 基于静息态功能性磁共振成像的海马旁回功能分区方法 |
CN108492873A (zh) * | 2018-03-13 | 2018-09-04 | 山东大学 | 一种用于辅助诊断阿兹海默症的知识迁移学习方法 |
CN108596879A (zh) * | 2018-04-02 | 2018-09-28 | 北京工业大学 | 一种基于希尔伯特黄变换的fMRI时频域动态网络构建方法 |
-
2018
- 2018-11-23 CN CN201811412834.8A patent/CN109558906B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5887074A (en) * | 1996-12-13 | 1999-03-23 | Siemens Corporate Research, Inc. | Local principal component based method for detecting activation signals in functional MR images |
WO2009058915A1 (en) * | 2007-10-29 | 2009-05-07 | The Trustees Of The University Of Pennsylvania | Computer assisted diagnosis (cad) of cancer using multi-functional, multi-modal in-vivo magnetic resonance spectroscopy (mrs) and imaging (mri) |
CN101706561A (zh) * | 2009-11-20 | 2010-05-12 | 电子科技大学 | 功能性磁共振图像的聚类方法 |
CN104700120A (zh) * | 2015-03-23 | 2015-06-10 | 南京工业大学 | 一种基于自适应熵投影聚类算法的fMRI特征提取及分类方法 |
CN106204545A (zh) * | 2016-06-28 | 2016-12-07 | 南方医科大学 | 一种基于区域划分和Fisher向量的医学病灶图像特征表达方法 |
CN107392907A (zh) * | 2017-09-01 | 2017-11-24 | 上海理工大学 | 基于静息态功能性磁共振成像的海马旁回功能分区方法 |
CN108492873A (zh) * | 2018-03-13 | 2018-09-04 | 山东大学 | 一种用于辅助诊断阿兹海默症的知识迁移学习方法 |
CN108596879A (zh) * | 2018-04-02 | 2018-09-28 | 北京工业大学 | 一种基于希尔伯特黄变换的fMRI时频域动态网络构建方法 |
Non-Patent Citations (3)
Title |
---|
F. SEGOVIA AT AL.: "fMRI data analysis using a novel clustering technique", 《 2009 IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD (NSS/MIC)》 * |
JIANG ZHANG ET AL.: "Analysis of fMRI Data Using an Integrated Principal Component Analysis and Supervised Affinity Propagation Clustering Approach", 《 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
王秀美等: "一种基于高斯隐变量模型的分类算法", 《计算机学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113657527A (zh) * | 2021-08-23 | 2021-11-16 | 广东工业大学 | 一种基于贝叶斯非参模型的脑功能区划分方法 |
CN113657527B (zh) * | 2021-08-23 | 2023-07-07 | 广东工业大学 | 一种基于贝叶斯非参模型的脑功能区划分方法 |
CN114854920A (zh) * | 2022-05-06 | 2022-08-05 | 浙江大学 | 高斯混合模型嵌入的gru自编码器高炉异常监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109558906B (zh) | 2022-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109697718A (zh) | 一种基于图论的自闭症检测方法及装置 | |
Sui et al. | An ICA-based method for the identification of optimal FMRI features and components using combined group-discriminative techniques | |
Carlson et al. | Patterns of activity in the categorical representations of objects | |
Thirion et al. | Detection of signal synchronizations in resting-state fMRI datasets | |
CN102509123B (zh) | 一种基于复杂网络的脑功能磁共振图像分类方法 | |
CN101706561B (zh) | 功能性磁共振图像的聚类方法 | |
Reddy et al. | An efficient multilevel thresholding scheme for heart image segmentation using a hybrid generalized adversarial network | |
Paul et al. | RETRACTED ARTICLE: Computer aided diagnosis of brain tumor using novel classification techniques | |
CN109558906A (zh) | 基于高斯隐变量降维聚类中心的功能磁共振聚类成像方法 | |
CN114469120A (zh) | 一种基于相似度阈值迁移的多尺度Dtw-BiLstm-Gan心电信号生成方法 | |
Aurangzeb et al. | An efficient and light weight deep learning model for accurate retinal vessels segmentation | |
Amarapur | An automated approach for brain tumor identification using ANN classifier | |
Kushol et al. | Retinal blood vessel segmentation from fundus image using an efficient multiscale directional representation technique Bendlets | |
Sujit et al. | Automated image quality evaluation of structural brain magnetic resonance images using deep convolutional neural networks | |
Ahmmed et al. | An advanced algorithm combining SVM and ANN classifiers to categorize tumors with position from brain MRI images | |
Björnsdotter et al. | A memetic algorithm for selection of 3d clustered features with applications in neuroscience | |
Anand Kumar et al. | 3D deep learning for automatic brain MR tumor segmentation with T-spline intensity inhomogeneity correction | |
Mehena et al. | Medical image segmentation and detection of MR images based on spatial multiple-kernel fuzzy C-means algorithm | |
Arora et al. | Noise adaptive FCM algorithm for segmentation of MRI brain images using local and non-local spatial information | |
Anitha et al. | WML detection of brain images using fuzzy and possibilistic approach in feature space | |
Abdulmunem | Brain MR Images Classification for Alzheimer’s Disease | |
Cherradi et al. | Brain extraction and fuzzy tissue segmentation in cerebral 2D T1-weigthed magnetic resonance images | |
Balasubramanian et al. | Automatic diagnosis and classification of glaucoma using hybrid features and k-nearest neighbor | |
Zafar et al. | Classification of fMRI data using support vector machine and convolutional neural network | |
CN107194934A (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 |