CN106504254B - 一种基于空间混合模型的颅脑磁共振图像分割方法 - Google Patents
一种基于空间混合模型的颅脑磁共振图像分割方法 Download PDFInfo
- Publication number
- CN106504254B CN106504254B CN201610903015.8A CN201610903015A CN106504254B CN 106504254 B CN106504254 B CN 106504254B CN 201610903015 A CN201610903015 A CN 201610903015A CN 106504254 B CN106504254 B CN 106504254B
- Authority
- CN
- China
- Prior art keywords
- vector
- model
- magnetic resonance
- image
- resonance image
- 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
Images
Classifications
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于空间混合模型的颅脑磁共振图像分割方法,包括:对输入的颅脑磁共振图像进行预处理,得到剔除非脑部组织后的图像;将预处理后的图像数据转换成行向量作为输入向量;对向量化后的图像数据利用基于空间关系的Inverted Dirichlet混合模型进行建模,并采用k‑means方法进行模型的初始化;采用贝叶斯变分推导技术估算混合模型的相关参数;通过计算后验概率从而确定各个像素点的所属类别,并得到新的标签向量;将标签向量作为输出向量并转换成灰度矩阵,从而得到最终的分割结果。本发明的技术方案对颅脑磁共振图像的分割效果出色、鲁棒性强,能够提高医学影像诊断的准确性。
Description
技术领域
本项发明属于计算机医学图像分析领域,具体涉及一种基于空间关系的InvertedDirichlet混合模型的颅脑磁共振医学图像的分割方法。
背景技术
近年来,医学影像技术发展迅速,特别是磁共振成像(Magnetic ResonanceImaging,MRI)技术,由于其具有无创伤性、检查范围覆盖人体各个系统、成像数据丰富等优点,应用最为广泛。医学图像分割是医学影像分析中的一个重要研究内容,其主要目的是根据医学图像区域内的相似性和区域间的不同特性将图像分割成若干区域。此外,由于颅脑是人体的一个重要器官,其拥有支配和控制人类思维活动的中枢神经系统,也是肿瘤、炎症和多发性硬化症等疾病的多发部位,因此颅脑磁共振图像分割在医学图像分割领域尤为重要,并能够在辅助医师临床诊断、计算机辅助手术等方面起到重要作用。
颅脑磁共振图像分割技术对脑肿瘤、脑炎性病变、脑白质病变、脑梗塞、脑先天性异常等疾病的诊断能起到关键性辅助作用。颅脑磁共振图像分割方法主要包括以下三类:1)基于区域的方法,此方法利用图像相邻像素间的局部关系来进行分割,因此适用于图像灰度值均匀条件下的目标分割,对于噪声高的医学图像分割效果显著下降;2)基于先验形状的方法,此方法的分割效果依赖于初始区域的选择,因此该方法存在一定的局限性,往往只适合特定形状的图像分割;3)基于像素的方法,此方法依据图像的灰度信息,多利用统计模型实现分割,其方法灵活,并能适应各种形状的医学图像分割。
混合模型是统计概率模型的一种,可将其看作为几个选定分布的叠加,每个分布都有各自的加权系数。因此,基于混合模型的图像分割其核心就是估算混合模型其成员分布的参数及它们的加权系数。最近几十年来,基于混合模型的颅脑磁共振图像分割的研究也取得了很多的成果,代表性的有:G.McLachlan等人(G.McLachlan and D.Peel,FiniteMixture Models.New York:Wiley,2000.)提出的利用期望最大化(EM)算法求解有限高斯混合模型(GMM)进行图像分割。但是该方法却没有考虑到一幅图像内相邻像素的空间关系,因此对噪声很敏感,分割效果不理想。D.Peel等人(D.Peel and G.McLachlan,“Robustmixture modeling using the t dis-tribution”,Stat.Comput.,vol.10,pp.339–348,2000.)利用EM算法对学生分布混合模型(SMM)求解。该方法相比于GMM有较强的鲁棒性,因此分割效果进一步提高,但是该算法复杂度较大。K.Blekas等人(K.Blekas,A.Likas,N.P.Galatsanos,and I.E.Lagaris,“A spatially constrained mixture model forimage segmentation”,IEEE Trans.Neural Netw.,vol.16,no.2,pp.494–498,2005)提出了一种空间变化有限混合模型(SVFMM),该方法利用Gibbs进行采样,并采用最大后验概率进行参数估计,但该方法需要额外的计算步骤来保证混合模型的混合系数满足数值为正且相加为一的条件,因此造成模型复杂度高、计算量大等缺点。
总体来说,在颅脑磁共振图像分割中,需要解决计算量大、噪声影响这两大难题,前者属于空间跟时间开销问题,后者属于准确度问题。因此,提出一种简单高效且抗噪音能力强的颅脑磁共振图像分割方法是当前的热点研究方向之一。
发明内容
本发明的目的在于克服现有技术的不足,提出一种基于空间混合模型的颅脑磁共振图像分割方法,通过基于空间关系的有限Inverted Dirichlet混合模型的建立和标签向量的获取,在计算开销低的情况下快速准确地得到颅脑磁共振图像的分割结果。
本发明解决其技术问题所采用的技术方案是:
一种基于空间混合模型的颅脑磁共振图像分割方法,包括:
对待分割的原颅脑磁共振图像进行预处理操作:利用分水岭算法去除非脑部组织,获得预处理图像;
将所述预处理图像转换成行向量;
采用K-means方法对所述行向量进行初始聚类,并设置聚类数K=3,分别表示灰质、白质和脑脊液;
基于空间关系的Inverted Dirichlet混合模型对行向量建立模型,并将所述行向量作为建立模型的输入向量;
利用贝叶斯变分推导方法计算所述建立模型的参数;
利用贝叶斯最大后验概率准则得到输入向量对应的标签向量的后验概率矩阵;
将后验概率矩阵转换成标签向量,然后转换成与原颅脑磁共振图像相同的矩阵形式,得到最终分割结果。
所述标签向量表示所述预处理图像中某一像素点属于白质、灰质或脑脊液的概率。
基于空间关系的Inverted Dirichlet混合模型对向量数据建立的模型用如下方式表示:
所述标签向量用如下方式表示:
本发明具有如下有益效果:
(1)本发明所提出的数据模型基于Inverted Dirichlet混合模型,InvertedDirichlet分布是一个多元泛化的inverted Beta分布,相比其它常用的分布(如高斯分布)它能够包含对称与非对称模态,因此对于图像分割更具有灵活性、准确性;
(2)本发明在Inverted Dirichlet混合模型的基础上加入了像素之间的空间关系,从而使得本分割方法具有更强的抗噪性;
(3)求解Inverted Dirichlet混合模型过程中,早期的方法多是采用期望最大化(EM)或复杂的采样法(如:Gibbs采样),这些方法往往使得模型求解复杂度高,甚至有可能陷入局部最优;针对这些问题,本发明采用了贝叶斯变分推导方法进行模型参数的估算,使得模型求解更加准确高效。
以下结合附图及实施例对本发明作进一步详细说明,但本发明的一种基于空间混合模型的颅脑磁共振图像分割方法不局限于实施例。
附图说明
图1为本发明方法的流程示意图;
图2是IBSR05数据集中,标签号为47的真实颅脑磁共振图像的分割效果对比图;
图3是加入噪声的模拟图的各个算法分割效果对比图。
具体实施方式
如图1所示,一种基于空间混合模型的颅脑磁共振图像分割方法,步骤如下:
步骤101,对待分割的原颅脑磁共振图像进行预处理操作:利用分水岭算法去除非脑部组织,获得N行M列的预处理图像;
步骤102,将所述预处理图像转换成行向量;
具体的,按照如下步骤将预处理图像转换成行向量:
先按照列的形式读取预处理图像矩阵,得到一个N*M维的列向量;
将N*M维的列向量转置成相应的行向量。
步骤103,采用K-means方法对所述行向量进行初始聚类,并设置聚类数K=3,分别表示灰质、白质和脑脊液;
步骤104,基于空间关系的Inverted Dirichlet混合模型对行向量建立模型,并将所述行向量作为建立模型的输入向量;
步骤105,利用贝叶斯变分推导方法计算所述建立模型的参数;
步骤106,利用贝叶斯最大后验概率准则得到输入向量对应的标签向量的后验概率矩阵;
步骤107,将后验概率矩阵转换成标签向量,然后转换成与原颅脑磁共振图像相同的矩阵形式,得到最终分割结果。
具体的,采用变分贝叶斯推导进行混合系数的相关参数估计,其中基于空间关系的有限Inverted Dirichlet混合模型如下:
采用Dirichlet compound multinomial分布作为混合系数的先验概率,对于每一个向量都有一个隐藏变量Zij∈{0,1},当属于第j个分布时候或者说是属于某一类(白质、灰质或脑脊液),Zij=1,否则Zij=0,其基于ξi的后验分布定义为:
其中,M为数据的个数,这里指像素点的个数,ξi服从Dirichlet的先验分布,定义如下:
此处的K值表示我们需要分的类别数目,即白质、灰质和脑脊液三类,其中Λij为空间关系模型,定义如下:
公式(5)的关系模型实际上是考虑周围像素的属于同一类的概率的影响,从而鲁棒性或者去噪能力增强,其中Ωi表示第i个像素的邻域大小,表示第t-1次循环后第i个像素点部分属于第j类(白质、灰质还是脑脊液)的概率,本发明采用领域大小为5*5的区域即周围的25个像素点进行考虑,为平滑参数,如下:
对于Zij参数估计采用如下方式:
公式(7)式是采用随机场理论和变分贝叶斯得出的用于估计(3)式的一个近似分布,它是潜在变量。其中rij的第t次迭代,即rij (t)采用如下式更新:
上述计算式中的数学期望值可以通过以下计算公式(采用变分贝叶斯方法得出的中间值)得出:
最终得到期望值后验概率矩阵,也就是所求的标签向量,它表示一幅图中的某一个像素点属于白质、灰质、脑脊液的概率分别是多少,然后我们取概率大的那个作为最终标签:
循环迭代以下步骤1~步骤2,直到收敛或者达到预设值。
步骤1,利用当前模型的预设值,计算等式(15)~(18);
本实施例中,发明的有益效果可以通过以下仿真实验得到验证:
1、仿真条件:
在CPU为Intel(R)Core(TM)i5-4590处理器:主频3.30GHZ,内存8G,操作系统:WINDOWS 10,仿真平台:MATLAB 2014Ra。
2、评价指标:
错分率(MCR)被用来对各种算法的分割效果进行量化评价,MCR的取值范围为[0,100],MCR越小表示分割结果越准确。MCR计算公式如下:
另外,为了分别衡量灰质、白质、脑脊液的分割正确率,引入了Dice metric方法,以Dice metric衡量灰质分割的准确性为例,TP表示被正确分类的灰质的数量,FP表示灰质被分成白质或脑脊液的数量,FN表示白质或脑脊液被分成灰质的数量。同理,Dice metric可以用来衡量白质和脑脊液部分的分割效果。Dice metric的取值范围为0到100,Dicemetric的值越大表示分割结果越准确。计算公式如下:
3、仿真对比方法
对比方法1:
G.McLachlan等人在文献“G.McLachlan and D.Peel,Finite MixtureModels.New York:Wiley,2000.”提出利用期望最大化(EM)算法求解有限高斯混合模型进行图像分割的方法,简称GMM。
对比方法2:
D.Peel等人在文献“Robust mixture modeling using the t dis-tribution,Stat.Comput.,vol.10,pp.339–348,2000.”提出利用EM算法求解学生分布混合模型的方法,简称SMM。
对比方法3:
K.Blekas等人在文献“A spatially constrained mixture model for imagesegmentation”,IEEE Trans.Neural Netw.,vol.16,no.2,pp.494–498,2005.”提出了一种利用Gibbs采样求解空间变化有限混合模型的方法,简称SVFMM。
4、仿真图像
仿真图像采用真实图2(a)(IBSR05数据集,标签号:47)和模拟图3(a)(Brainweb数据集,标签号:90,噪声:9%)为去除非脑部组织后的测试磁共振图像;其中图2(b)和图3(b)分别为对应的标准分割结果;图2(c)-图2(f)和图3(c)-图3(f)是分别采用GMM、SMM、SVFMM和本发明方法(IDMM-SC)分割后的效果图。
5、仿真内容
仿真1,利用现有的3种对比方法和本发明的方法对图2(a)进行图像分割,得到对应的分割效果图见图2(c)-图2(f),另外仿真得到的数据结果见表1。
图2(a)在四种分割方法下的分割结果对比,其中CSF、GM和WM分别对应脑脊液、灰质和白质:
表1
结合图2(c)-图2(f)以及表1,不难看出本发明的方法无论在CSF、GM、WM中的Dicemetric还是MCR,都表现出不错的分割效果,尤其是在CSF跟GM方面提升更大。错分率MCR提高近14%,由此可以充分的说明本发明方法较其它三种方法对颅脑磁共振图像分割更具有优势。
仿真2,利用现有的3种对比方法和本发明的方法对图3(a)进行图像分割,得到对应的分割效果图见图3(c)-图3(f),另外仿真得到的数据结果见表2。
表2
结合图3(c)-图3(f)跟表2,本发明的方法无论在CSF、GM、WM中的Dice metric还是MCR,都表现出理想的分割效果,每个指标的正确分割率都有90%以上,较其它三种方法有近3%正确分割率的提升,正确分割率达到94%左右(或者错分率MCR 6%)。
综上,本发明采用了加入空间关系的inverted Dirichlet混合模型作为图像分割方法,具有抗噪强、分割结果准确等优点,非常适合应用于颅脑磁共振图像的分割。
以上仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于空间混合模型的颅脑磁共振图像分割方法,其特征在于,包括:
对待分割的原颅脑磁共振图像进行预处理操作:利用分水岭算法去除非脑部组织,获得预处理图像;
将所述预处理图像转换成行向量;
采用K-means方法对所述行向量进行初始聚类,并设置聚类数K=3,分别表示灰质、白质和脑脊液;
基于空间关系的Inverted Dirichlet混合模型对行向量建立模型,并将所述行向量作为建立模型的输入向量;
利用贝叶斯变分推导方法计算所述建立模型的参数;
利用贝叶斯最大后验概率准则得到输入向量对应的标签向量的后验概率矩阵;
将后验概率矩阵转换成标签向量,然后转换成与原颅脑磁共振图像相同的矩阵形式,得到最终分割结果;
所述标签向量表示所述预处理图像中某一像素点属于白质、灰质或脑脊液的概率;
基于空间关系的Inverted Dirichlet混合模型对向量数据建立的模型用如下方式表示:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610903015.8A CN106504254B (zh) | 2016-10-17 | 2016-10-17 | 一种基于空间混合模型的颅脑磁共振图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610903015.8A CN106504254B (zh) | 2016-10-17 | 2016-10-17 | 一种基于空间混合模型的颅脑磁共振图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106504254A CN106504254A (zh) | 2017-03-15 |
CN106504254B true CN106504254B (zh) | 2020-03-10 |
Family
ID=58293804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610903015.8A Active CN106504254B (zh) | 2016-10-17 | 2016-10-17 | 一种基于空间混合模型的颅脑磁共振图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106504254B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110503655B (zh) * | 2019-08-21 | 2022-04-01 | 山东师范大学 | 一种基于稀疏贝叶斯模型和多图谱融合的mri神经胶质瘤分割方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015039054A1 (en) * | 2013-09-13 | 2015-03-19 | The Regents Of The University Of California | Method and system for analysis of volumetric data |
CN106023236A (zh) * | 2016-06-16 | 2016-10-12 | 华侨大学 | 基于截断Dirichlet过程无限Student’s t 混合模型的脑部核磁共振图像分割方法 |
-
2016
- 2016-10-17 CN CN201610903015.8A patent/CN106504254B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015039054A1 (en) * | 2013-09-13 | 2015-03-19 | The Regents Of The University Of California | Method and system for analysis of volumetric data |
CN106023236A (zh) * | 2016-06-16 | 2016-10-12 | 华侨大学 | 基于截断Dirichlet过程无限Student’s t 混合模型的脑部核磁共振图像分割方法 |
Non-Patent Citations (2)
Title |
---|
A non-local fuzzy segmentation method:Application to brain MRI;Benoit Caldairou等;《Pattern Recognition 44》;20110930;1916-1927 * |
Gaussian mixture model based segmentation methods for brain MRI images;M.A.Balafar;《Artificial Intelligence Review》;20140331;429-439 * |
Also Published As
Publication number | Publication date |
---|---|
CN106504254A (zh) | 2017-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Khalili et al. | Automatic brain tissue segmentation in fetal MRI using convolutional neural networks | |
US20220122250A1 (en) | Brain feature prediction using geometric deep learning on graph representations of medical image data | |
CN111488914B (zh) | 一种基于多任务学习的阿尔茨海默症分类及预测系统 | |
CN113616184B (zh) | 基于多模态磁共振图像的脑网络建模与个体预测方法 | |
Biffi et al. | Explainable anatomical shape analysis through deep hierarchical generative models | |
US20230301542A1 (en) | Brain atlas individualization method and system based on magnetic resonance and twin graph neural network | |
Zhu et al. | Multi-modal AD classification via self-paced latent correlation analysis | |
Wu et al. | Registration of longitudinal brain image sequences with implicit template and spatial–temporal heuristics | |
Zöllner et al. | Kidney segmentation in renal magnetic resonance imaging-current status and prospects | |
Goncalves et al. | Self-supervised MRI tissue segmentation by discriminative clustering | |
Galimzianova et al. | Stratified mixture modeling for segmentation of white-matter lesions in brain MR images | |
Yang et al. | DBAN: adversarial network with multi-scale features for cardiac MRI segmentation | |
Li et al. | Learning brain connectivity sub-networks by group-constrained sparse inverse covariance estimation for Alzheimer's disease classification | |
Qiu et al. | 3D MR ventricle segmentation in pre-term infants with post-hemorrhagic ventricle dilatation (PHVD) using multi-phase geodesic level-sets | |
Zhu et al. | Statistical learning methods for neuroimaging data analysis with applications | |
Liu et al. | AHU-MultiNet: adaptive loss balancing based on homoscedastic uncertainty in multi-task medical image segmentation network | |
Ossenberg-Engels et al. | Conditional generative adversarial networks for the prediction of cardiac contraction from individual frames | |
CN106504254B (zh) | 一种基于空间混合模型的颅脑磁共振图像分割方法 | |
Fan et al. | A survey of emerging applications of diffusion probabilistic models in mri | |
Xu et al. | 3D fetal pose estimation with adaptive variance and conditional generative adversarial network | |
Hu et al. | Accurate brain age prediction model for healthy children and adolescents using 3d-cnn and dimensional attention | |
CN106709921B (zh) | 一种基于空间Dirichlet混合模型的彩色图像分割方法 | |
Kasim et al. | Gaussian mixture model-expectation maximization algorithm for brain images | |
da Silva et al. | Enhanced pre-processing for deep learning in MRI whole brain segmentation using orthogonal moments | |
CN112950654A (zh) | 基于多核学习与超像素核低秩表示的脑肿瘤图像分割方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |