《隐马尔可夫场的脊柱CT图像分割算法范文》由会员分享,可在线阅读,更多相关《隐马尔可夫场的脊柱CT图像分割算法范文(7页珍藏版)》请在金锄头文库上搜索。
1、【Word版本下载可任意编辑】 隐马尔可夫场的脊柱CT图像分割算法范文 0 引 言 近年来不断发展的医学成像条件令成像数量不断增多,这使得医生逐片解读图像来判定疾病的过程越来越繁琐,又因为现代人脊柱畸形的患病率正逐年上升,所以脊柱图像的分割将会是医学图像处理的一个研究热点和方向。考虑到分割效果将直接影响到医生对患者病情的判断,所以寻找一种具有更高分割精度的算法迫在眉睫。 现有的自动分割技术主要包括边缘检测法、区域生长法、自适应阈值法、模糊聚类法、基于数学形态学方法、基于人工智能的方法、基于神经网络法等。由于传统的分割算法往往不注重像素点内在联系,即仅仅利用像素的灰度信息,不对像素的空间信息开展
2、考虑,所以导致对脊柱CT图像分割结果不理想。近年来为了解决这一缺陷隐马尔可夫模型(hidden Markov model, HMM)成为图像分割的研究热点之一。原因是这一方法将目标像素与其邻域像素间的相互信息作为先验知识,同时结合最大后验准则(maximum a posteriori, MAP),把图像分割问题转化成了能量函数最优问题。鉴于此,基于隐马尔可夫随机场的图像分割方法在医学图像分割方面更具有实用性。 文在图像与统计力学系统之间做了一个类比,提出了著名的 Hammersley-Clifford 定理,从理论上证明MRF与Gibbs 分布是等价的,使得 MRF理论可以被用到图像分割等众多
3、领域来解决实际的难题。在此之后众多研究者开始不断地改良MRF模型,文将其应用到院舜殴舱裢枷竦姆指钪腥Dagan等在隐马尔可夫随机场模型中引入了期望最大化算法(expectation-maximization,EM),然而这一方法对参数初始值的选定非常敏感。文应用马尔可夫随机场的方法来分割脊柱的CT图像从而自动推理出了脊柱关节,但该算法存在一定程度的过度分割问题。 本文中通过邻域来编码空间信息,在建立HMM的同时,将模拟退火算法与期望最大化算法融合并结合最大后验概率的思想来分割CT图像中胸椎的轮廓,为医生后续的诊断工作提供一定理论依据。 1 隐马尔可夫随机场模型 2 基于MSA-EM模型确定算法
4、 用EM算法来获得模型参数的优点是代码便捷,但因为隐马尔可夫模型的最大似然方程一般会有多个根,这导致EM算法容易陷入局部极值。模拟退火算法(simulated annealing,SA)在全局优化解决局部极值问题上有一定优势,但传统SA算法不能保证最终结果一定比搜索中曾放弃的值小,因此本文提出一种MSA-EM(EM algorithm based simulated annealing with memory, MSA-EM)算法,在执行过程中参加记忆器,用于记录执行过程中的最优解。基于模拟退火算法的EM算法(MSA-EM)具体步骤如下: 3 基于MAP的图像分割 最大后验准则的基本思想是用最
5、大后验概率估计来近似图像的真实分类情况。具体步骤如算法2所示: 算法第一步是初始化,把算法MSA-EM的结果作为MAP估计的初始值。因为后验概率与后验能量函数存在一种反比例关系,即: 所以第二步中用迭代的方法求解最小后验能量函数。 4 实验结果分析 本文实验的硬件环境Intel(R)Core(TM)i7-4790 CPU3.60GHz,4GB内存。算法在VS20*与MatlabR20*b混编环境下实现。实验数据采用合作医院的双源64排CT机开展扫描从而获取患者患病部位的数据,然后通过图像处理技术转化为512*512的二维图像,本章将利用多种评价指标,在95幅胸椎、20幅腰椎和5幅颈椎CT图像上
6、展开实验,对分割框架开展性能评估 采用Jaccard、Dice、Correct classification ration(CCR)三种评估指标开展实验评估。三种定量评估指标定义如下: 其中,SEG和GT分别表示算法自动分割的结果和专家分割的参考结果。三种指标系数的取值范围均为0到1,越接近1表示分割效果越好。 对图像用HMRF-EM、eHMRF算法、本文的HMRF- MSA-EM算法开展分割实验,如图2所示。 沿用Tohka等的处理方法对图2(a)中的左上图添加不同噪声层次的噪声后,分别用三种分割算法开展处理,对处理后的图像的Jaccard、Dice和CCR系数开展了比照。 从表1、表2、表
7、3可以得出,本文算法在不同噪声层次均取得更高的Jaccard、Dice、CCR值。随着噪声层次的不断增加, eHMRF算法的Dice系数值下降了0.0871,而本文算法下降了0.0869,说明本文算法不仅具有更高的准确度还有更好的稳定性。 对50幅胸椎CT图像采用eHMRF和本文的HMRF-MSA-EM算法开展分割,并计算Jaccard、Dice和CCR三种定量评估指标,结果如图3所示。 由图3可以看出本文算法对胸椎图像的分割效果具有普遍的提高。Jaccard系数值的中位数依次为0.8706(HMRF-EM),0.8778(eHMRF),0.8828(HMRF-MSA-EM),Dice系数的中
8、位数为0.9316(HMRF-EM),0.9361(eHMRF),0.9383(HMRF-MSA-EM),CCR系数值的中位数0.9347(HMRF-EM),0.9414(eHMRF),0.9449(HMRF-MSA-EM)。 经过对实验数据的分析和统计,得到三种分割方法的执行时间比照图如图4所示。 其中,y轴表示图像分割所用的时长,x轴以顺时针顺序分别对应图2(a)中的四幅图像。第一列代表HMRF-EM算法所用时长,第二列是eHMRF算法所用时长,第三列是HMRF-MSA-EM算法所用时长。在对四幅脊柱CT图像的处理中本文算法比HMRF-EM算法耗时少得多,且没有出现空洞和不连续区域。本文算
9、法在分割效果上与eHMRF算法相似,但在执行时间上有优势。 5 结 论 本文对基于隐马尔可夫随机场的医学影像分割算法开展了改良,首先构造了一个隐马尔可夫随机场模型,随后引入带有记忆器的模拟退火算法来优化求解模式,以此改良期望最大化算法,最后建立了与之对应的MAP准则。最终通过实验比照证明了本文算法对脊柱分割的准确性有一定程度的提高,实现了重现性,具有很好的可行性。虽然目前CT成像应用极其广泛,但因为成像原理等原因其仍在局部环节存在一些无法克服的劣势,因此怎样融合多种类型图像的优势,探寻多种模态图像的分割将是下一步的研究重点。 参 考 文 献: CHEN Hejun, DING Haiqiang
10、, HE Xiongxiong, et al. Color Image Segmentation Based on Seeded Region Growing with Canny Edge Detection. 20*:683-686. NARAYANASWAMY A, DWARAKAPURAM S, BJORNSSON C, et al. Robust Adaptive 3-D Segmentation of Vessel Laminae from Fluoresce-nce Confocal Microscope Images and Parallel GPU Implementatio
11、n. Microscopy & Microanalysis, 20*, 29(3):583-97. FOOLADIVANDA A, SHOKOUHI S, AHMADINEJAD N, et al. Automatic Segmenta-tion of Breast and Fibroglandular Tissue in Breast MRI Using Local Adaptive Thresholding/ Biomedical Engineering. IEEE, 20*. 马文萍, 黄媛媛, 李豪,等. 基于粗糙集与差分免疫模糊聚类算法的图像分割. 软件学报, 20*, 25(11)
12、:2675-2689. KUNDN A, CHATTERJEE R.K.Morphological Scale-Space Based Vessel Segmentation of Retinal Image/ India Conference. IEEE, 20*:986-990. JAISWAL V. A. Survey of Lmage Segmentation Based on Artificial Intelligence and Evolutionary Approach. Ai Magazine, 20*, 15(3):71-78. 郑欣, 彭真明. 基于活跃度的脉冲耦合神经网络
13、图像分割. 光学精细工程, 20*, 21(3):821-827. HELD K, KOPS E.R, KRAUSE B. J, et al. Markov Random Field Segmentation of Brain MR Images. IEEE Transactions on Medical Imaging, 1997, 16(6):878-86. GEMAN S. GEMAN D. Stochastic Relaxation Gibbs Distribution and the Bayesian Restoration of Image.IEEE Trans Pattern A
14、nal Machine Intell,1991;13(5):401-412. ZHANG Yongyue, M. Brady, and S. Smith. Segmentation of Brain MR Images Through a Hidden Markov Random Field Model and the Expectation-maximization Algorithm. IEEE Transaction on Medical Imaging, 20*,20(1):45-57. ZHANG Tong, XIA Yong. An Evolutionary HM-RFapproa
15、ch to Brain MR Image Segmentation Using Clonal Selection Algorithm.Biological and Medical Systems,20*, 8(1): 6-11. KADOURY S., LABELLE H., PARAGIOS N. Automatic Inference of Articulated Spine Models in CT Images Using High-order Markov Random Fields. Medical Image Analysis, 20*, 15(4):426-437. KADOURY S., LABELLE H., PARAGIOS N. Spine Segmentation in Medical Images Using Manifold Embeddings and Higher-order MRFs. IEEE Transactions on Medical Imaging, 20*, 32(7):1227-1238. BROGNAUX S.,DRUGMAN T. HMM-Based Spee