基于特征向量的颅脑CT图像分割方法
- 格式:pdf
- 大小:447.20 KB
- 文档页数:6
医学图像配准与分割方法研究与实践医学图像配准与分割是医学影像处理中的两个重要任务,可以帮助医生更准确地诊断和治疗疾病。
图像配准旨在将多个医学图像对齐,使得它们在空间上或时间上完全或部分重叠。
而图像分割则是将医学图像中的组织、器官或异常区域进行分离和区分。
本文将介绍医学图像配准和分割的一些常见方法,并探讨其在实际应用中的研究与实践。
一、医学图像配准方法研究与实践医学图像配准是通过将多个医学图像进行准确对齐,实现多图像之间的一致性。
在医学影像处理中,常见的配准方法包括基于特征的方法和基于相似度度量的方法。
1. 基于特征的方法基于特征的配准方法通过提取图像中的特征点或特征区域来进行图像匹配。
其中,最常见的特征包括角点、边缘、纹理等。
常用的算法包括SIFT(尺度不变特征变换)、SURF(加速稳健特征)和ORB(旋转不变特征)等。
这些方法通过对特征点的检测、描述和匹配来实现图像配准。
在实践中,基于特征的方法往往具有较好的配准效果和鲁棒性。
2. 基于相似度度量的方法基于相似度度量的配准方法通过计算图像间的相似度来进行匹配。
常用的相似度度量指标包括互信息、相关系数和互相关等。
这些方法通过优化相似度度量指标来寻找最优的变换参数,从而实现图像的配准。
在实践应用中,基于相似度度量的方法通常能够得到较好的配准结果,并且具有较高的计算效率。
二、医学图像分割方法研究与实践医学图像分割是将医学图像中的组织、器官或异常区域进行分离和区分的过程。
常见的医学图像分割方法包括基于阈值的方法、基于边缘的方法和基于区域的方法。
1. 基于阈值的方法基于阈值的分割方法是最简单也是最直观的方法之一。
它通过选取适当的阈值来将图像进行二值化,将感兴趣的组织或区域与背景进行区分。
常见的阈值分割方法有全局阈值法、自适应阈值法和区域生长法等。
尽管基于阈值的方法简单易用,但对于复杂图像和噪声较多的情况下,效果不佳。
2. 基于边缘的方法基于边缘的分割方法通过检测图像中的边缘信息进行分割。
医学影像处理中的像分割技术医学影像是医生诊断和治疗疾病的重要工具之一。
在医学影像处理中,像分割技术被广泛应用于从医学图像中分离和提取感兴趣的解剖结构或病变区域。
本文将介绍像分割技术的基本原理、常用方法以及在医学影像处理中的应用。
一、像分割技术的基本原理像分割是指将医学图像中的像素划分为不同的区域或类别。
它是基于图像的特征和像素间的差异来实现的。
在医学影像处理中,常用的像分割技术包括阈值分割、区域生长、边缘检测和基于模型的分割等。
1. 阈值分割阈值分割是一种直观简单的像分割方法,根据像素的灰度值与预设的阈值进行比较来实现分割。
例如,在CT图像中,我们可以通过设置阈值将肿瘤区域与正常组织分开。
2. 区域生长区域生长是通过种子点开始,根据一定的准则将相邻像素合并成区域的过程。
该方法基于像素之间的相似性,将相似的像素聚集成为相同的区域。
这种方法常用于分割有连续性的结构,比如血管和肿瘤。
3. 边缘检测边缘检测是通过检测图像中灰度值的变化来找到物体的边缘。
常用的边缘检测算法有Canny算子和Sobel算子等。
边缘检测在医学影像处理中被广泛应用于器官和组织的分割。
4. 基于模型的分割基于模型的分割是一种基于统计模型或数学模型的像分割方法。
它通过构建图像的概率模型来实现分割,常用的方法有水平集方法和活动轮廓模型等。
这种方法适用于复杂的图像分割任务,如肿瘤的分割。
二、常用的像分割方法除了基本的像分割技术外,还有一些常用的高级分割方法被应用于医学影像处理。
1. 基于机器学习的分割方法机器学习方法在医学影像处理中得到了广泛的应用。
例如,支持向量机(SVM)和随机森林(Random Forest)等分类算法可以通过学习训练数据集,实现医学图像中不同结构的分割。
2. 基于深度学习的分割方法深度学习方法近年来在医学影像处理领域取得了显著的成果。
卷积神经网络(CNN)被广泛应用于医学图像的像分割任务中。
通过训练大量的医学图像数据,深度学习网络可以自动学习特征,并实现高效准确的分割。
医疗影像处理中的图像分割与特征提取技术随着医疗技术的不断发展,医疗影像成为疾病诊断和治疗中重要的辅助工具。
医疗影像处理中的图像分割与特征提取技术在辅助医生进行准确诊断和治疗方案制定方面发挥着重要作用。
本文将重点介绍医疗影像处理中的图像分割和特征提取技术的原理和应用。
图像分割是将影像中的区域分离出来,将其划分为不同的组织或结构。
这对于医生来说非常重要,因为在分析和理解影像中的病变时,准确地定位和分割病变区域非常关键。
图像分割技术可以使用多种方法实现,包括阈值分割、边缘检测和区域生长等。
其中,阈值分割是最常用的方法之一,它通过选取一个阈值来将图像中的像素分为两个不同的区域。
在医疗影像处理中,图像分割的一个重要应用是肿瘤的定位和量化。
通过准确分割出肿瘤区域,医生可以更好地评估肿瘤的大小、形状和位置,从而制定更有效的治疗方案。
在肿瘤检测和分割中常用的技术包括基于阈值的直方图分割、基于区域生长的分割以及基于边缘检测的分割等。
特征提取是从分割后的图像中提取有用的特征信息,用于进一步的分析和诊断。
这些特征可以是形状、纹理、灰度分布等。
特征提取的目的是通过提取有关病变的定量和定性特征来帮助医生进行准确的诊断和治疗方案制定。
常用的特征提取方法包括基于灰度统计的方法、基于纹理的方法和基于形状的方法等。
在医疗影像处理中,特征提取的一个重要应用是癌症筛查和诊断。
通过从肿瘤区域中提取形状、纹理和灰度等特征,医生可以更好地评估肿瘤的类型和恶性程度。
例如,通过提取肿瘤的形状特征可以判断其是否呈现出典型的良性或恶性形态。
通过提取肿瘤的纹理特征可以识别不同类型的肿瘤,从而指导进一步的治疗方案制定。
除了肿瘤的定位和量化、癌症筛查和诊断,图像分割和特征提取技术在其他医疗领域也发挥着重要作用。
例如,在脑部影像分析中,图像分割和特征提取可以帮助医生识别脑部病变,并评估病变的严重程度。
在心脏影像分析中,图像分割和特征提取可以用于自动定位和分割心脏的不同组织结构,如心房和心室,从而帮助医生评估心脏的功能和病变。
医学CT图像分割方法研究共3篇医学CT图像分割方法研究1医学CT图像分割方法研究随着医学科技的不断发展,现代医学已经离不开医学图像。
其中,医学CT图像是医生在诊断和治疗疾病时使用最多的一种医学图像。
CT图像能够提供非常有用的解剖结构信息,但是在医学实践中,医生往往需要将图像中感兴趣的区域或者器官进行分割,这对于医生来说具有十分重要的临床意义。
因此,在现代医学中,医学CT图像的分割技术也得到了极大的关注。
目前,医学CT图像分割技术中应用比较广泛的方法主要有以下几种:1. 基于阈值分割法阈值分割法是常见的分割方法之一,它是根据医学图像灰度值的特性来分割图像的。
该方法通过设定阈值,将图像中大于和小于阈值的像素点分成两个区域,以此来实现图像分割。
该方法实现较为简单,但受到图像噪声和灰度逐渐变化等因素影响较大,因此在临床上并不是最理想的方法。
2. 基于边缘检测算法边缘检测算法也是医学CT图像分割中一种经典的方法。
它利用图像中的边缘信息来进行分割。
常用的边缘检测算法有Sobel、Prewitt、Canny等。
该方法具有较高的精度,但对于图像中灰度分布均匀、边缘不明显的情况下,其分割效果会受到较大的影响。
3. 基于区域生长的方法区域生长法是一种基于像素相似性的图像分割方法,它不需要设置任何参数,完全基于图像本身的特征来进行分割。
该方法首先选择种子点,然后通过像素相似性来将其它像素靠近种子点的分到同一个分割区域中。
该方法对图像中存在明显分界线的情况下有效,但对于噪声和灰度分布均匀的图像分割效果较差。
4. 基于机器学习的方法机器学习算法是近年来医学CT图像分割方法中应用十分广泛的一种技术。
包括支持向量机、随机森林等机器学习算法。
该方法通过采用大量已知的医学图像来训练分类器,然后使用训练好的分类器来对新的医学图像进行分割。
该方法准确度较高,但在训练分类器的过程中需要大量的数据和时间。
以上是当前医学CT图像分割方法中应用比较广泛的几种方法,其中每种方法都有其优缺点,并没有一种方法适用于所有的医学图像。
医学影像处理中的图像分割算法使用技巧医学影像处理是一门涉及医学图像采集、存储、处理和分析的学科。
医学图像中通常包含大量的信息,因此图像分割是医学影像处理中必不可少的一环。
图像分割是将医学图像中感兴趣的区域从背景中分离出来的过程,它通常用于检测病变区域、提取感兴趣的解剖结构或组织等。
在医学影像处理的图像分割中,有许多算法可供选择,下面将介绍一些常用的图像分割算法以及它们的使用技巧。
1. 阈值分割阈值分割是最简单且常用的图像分割方法之一。
它基于像素的灰度值,将图像中大于或小于特定阈值的像素分离出来。
阈值分割适用于图像中目标和背景的灰度值存在明显差异的情况,例如CT扫描中的骨骼分割。
在使用阈值分割时,需要根据图像的特点选择适当的阈值,并进行阈值的优化和调整,以获得更好的分割效果。
2. 区域生长区域生长是一种逐像素地将图像分割为几个连通区域的方法。
它通过选择种子点和定义生长准则来实现图像的分割。
区域生长适用于图像中目标的灰度值相似的情况,例如MRI图像中的脑部分割。
在使用区域生长时,需要选择适当的种子点,并根据具体情况设置生长准则,以获得准确的分割结果。
3. 边缘检测边缘检测是通过寻找图像中不连续的灰度值变化来实现图像分割的方法。
它可以准确地检测出图像中的边缘信息,并将其作为分割结果。
边缘检测适用于图像中目标的边界清晰的情况,例如X射线图像中的器官分割。
在使用边缘检测进行图像分割时,需要选择适当的边缘检测算法,并进行参数调整以获得满意的分割效果。
4. 水平线剖分水平线剖分是一种基于灰度值水平变化的分割方法。
它通过对图像的水平方向进行剖分和分析,将图像中的区域分隔开。
水平线剖分适用于图像中存在明显的水平变化的情况,例如胸部X射线图像中的肺部分割。
在使用水平线剖分进行图像分割时,需要选择适当的剖分方法,并进行参数的调整以获得理想的分割效果。
5. 基于机器学习的分割基于机器学习的分割方法是近年来发展起来的一种分割方法。
中北大学学位论文长,算法相对困难,但对于较复杂的狄度图像二值化具有良好的整体效果。
动态闽值二值化是比较复杂的方法,因为它的一些属性与局部阂值法相似,所以在一般的方法分类中也可以把它划分为局部闽值法一类。
局部和动态阈值二值化能够处理品质较差的图像,但同时具有费时和某些失真的缺点。
3.3改进的阈值分割算法颅骨的分割处理中,闽值的选择时一个关键问题。
阈值选取的合理与否直接影响到颅骨分割的结果。
理想的阈值应该是最大程度的突出颅骨中有用的边界特征而抑制干扰。
由于颅骨灰度分布较集中,颅骨区域灰度高于背景区域灰度,其直方图一般呈现背景和目标双峰,如图3.4所示,所以本节采用全局自动阈值分割法。
这主要是由于这种阈值分割算法简单,时间开销小,易于实现,划分目标和背景的错误分割概率最小。
原图幽3.4人脑cT直方剀中北大学学位论文图3.5改进的白适应迭代分割结果分析结果可以得出,改进前的迭代阈值法分割得到的内部区域完整,边界也完全封闭,几乎没有局部断裂的现象。
但是部分区域存在着过分割现象,处理结果偏大。
这是由于自适应迭代阈值分割算法在改进之前是对整幅图像求灰度均值,这样切片图像中的黑色背景会对其产生很大的影响,导致最佳阈值偏小,所以一些软组织没有被分割出去。
改进的自适应迭代阈值分割精度要比之前高很多,能更有效的分割骨骼和软组织,但是在部分区域出现了断裂的现象,如图3.6(图3.6为图3.5(c)局部放大)所示,红色标注区域边界断裂,边界信息丢失。
幽3_6改进的白适应迭代闽值法存在的问题Otsu方法及改进的Otsu方法对于图3.7(a)分割结果如下。
图3.7(b)为Otsu方法分割结果,T=47:图3.7(c)为改进的Otsu方法分割结果,改进的Otsu方法在第二次分中北大学学位论文割后找出了最优闽值z=76。
图3.7改进的Otsu闽值分割结果对图3.7(a)运用Otsu方法进行分割根本无法区别软组织和骨骼,但是改进的Otsu法较为改进之前的图像分割质量有了很大的提高,能够从背景中完全提取目标区域。
医学图像处理技术——CT和MRI图像的3D重建与分割技术在现代医学诊断中,医学图像处理技术已经成为必不可少的一部分。
医学图像处理技术可以通过对成像设备(如CT和MRI)获取的大量图像数据进行处理和分析,获取患者疾病的详细信息,从而为诊断和治疗提供关键的支持。
其中,CT和MRI图像的3D重建与分割技术是医学图像处理技术中的两个关键环节。
下面,本文将从技术原理、应用场景以及未来发展方向等方面综述CT和MRI图像的3D重建与分割技术。
技术原理3D重建技术是指将一系列二维图像数据通过一定的算法处理,从而还原成完整的三维模型。
而CT和MRI图像的3D重建主要是通过体素(voxel)的形式来完成的。
体素是三维空间(x、y、z)中的一个像素点,在体素极度密集的情况下,所构成的形状就趋近于真实的物体,可以达到较为真实的3D重建效果。
而3D分割技术,从字面上就能看出它的意义:将三维图像数据进行分离,实现对不同组织、不同器官、不同病变区域的有针对性的处理和分析。
在医学诊断中,正确、精准的分割技术能够提高治疗的效果,减少治疗的负担。
目前,基于深度学习(Deep Learning)和卷积神经网络(CNN)的3D分割技术也逐渐成为热点研究领域。
应用场景那么,在实际的医学诊断中,CT和MRI图像的3D重建与分割技术究竟能够发挥哪些作用呢?俯视全图,观察整体结构。
在医学图像处理中,仅能识别单张图片只能了解一部分结构,而通过多张CT和MRI图像,可以将一个器官或组织等的完整结构进行重建。
其中,3D重建技术能够快速准确重建三维模型,并依照组织器官分割的方式清晰地展示出图像结构的全貌。
指引精细区域,精准定位病灶。
在医学诊断中,CT和MRI图像的3D分割技术能够将患病组织和健康组织分隔开,帮忙医生更准确地定位病灶,促进后续治疗方案的制定和落实。
再者,对于某些难以定位的病灶,3D分割技术能够将其清晰可见,并辅以医生对其周围环境的分析,达到如实、精准、科学的治疗效果。
医疗影像处理中的图像分割算法图像分割是一种将图像划分成多个具有独立语义信息的区域的方法,它在医疗影像处理中扮演着重要的角色。
医疗影像分割的目标是将影像中感兴趣的区域从背景中分离出来,以帮助医生进行病变分析、诊断和治疗。
在医疗影像处理中,图像分割算法的准确性和效率至关重要。
准确性确保分割结果与医生给定的标注一致,而效率则决定了算法能否在实际应用中处理大规模医疗影像数据。
常见的医疗影像分割算法包括阈值分割、区域生长、边缘检测和基于机器学习的方法等。
下面将分别介绍这些算法的原理和应用。
1. 阈值分割阈值分割是一种简单但常用的图像分割方法。
它基于图像灰度值的不同,在图像上设置一个或多个阈值来实现分割。
通过选择合适的阈值,可以将感兴趣的区域与背景区域分离开来。
但该方法对于光照变化和噪声敏感,因此在复杂的医学影像中效果有限。
2. 区域生长区域生长是一种基于相似性的图像分割方法。
它从种子点开始,不断生长并将与种子点相似的像素合并为同一区域。
该方法通常需要人工提供种子点,并根据图像特点调整生长准则。
区域生长方法适用于具有明显边界的图像,但对于灰度均匀且模糊边界的图像分割效果较差。
3. 边缘检测边缘检测是一种通过检测图像中的强度变化来实现分割的方法。
它通过计算图像中像素灰度值的一阶或二阶导数来检测边缘。
常用的边缘检测算法有Sobel、Canny和Laplacian等。
该方法对于具有明显边界的图像分割效果较好,但在存在噪声或纹理较强的图像中容易产生误检。
4. 基于机器学习的方法基于机器学习的方法在医疗影像分割中越来越受到关注。
这些方法利用训练数据来构建分类器或分割模型,以实现自动分割。
常用的机器学习算法包括支持向量机(SVM)、随机森林和卷积神经网络(CNN)等。
这些算法通过学习大量的标注数据,可在医疗影像分割任务中取得较高的准确性和泛化能力。
除了以上常见的图像分割方法外,还有一些其他的算法,如基于水平集的方法、活动轮廓模型等。
医学影像中脑肿瘤自动检测与分割方法研究在医学影像领域,脑肿瘤的自动检测和分割是一项关键的任务,可以辅助医生进行初步诊断和治疗过程中的定位。
然而,由于脑肿瘤的复杂形态和变异性,以及大量的临床数据和图像信息,使得传统的手工分割方法无法满足高效和准确的需求。
因此,研究人员致力于开发自动检测和分割方法来提高脑肿瘤的诊断效率和准确性。
在医学影像中,脑肿瘤的自动检测和分割方法研究可以分为两个主要的方向:基于机器学习的方法和基于深度学习的方法。
基于机器学习的方法主要依赖于特征提取和分类器的构建。
首先,研究人员需要利用图像处理技术对医学影像进行预处理,如去噪、增强和平滑化等操作。
然后,通过选择合适的特征提取方法来提取图像中的特征,例如形状、纹理和强度等。
最后,使用机器学习算法来训练分类器,将提取的特征与脑肿瘤进行分类。
常用的机器学习算法包括支持向量机(SVM)、随机森林(Random Forest)等。
基于深度学习的方法则是近年来医学影像领域的新热点。
深度学习模型可以自动地学习到特征和分类器,并且在大规模数据集上表现出了强大的性能。
一种常用的深度学习模型是卷积神经网络(Convolutional Neural Network,CNN)。
使用卷积层和池化层等结构,CNN可以有效地提取图像中的特征,并进行分类和分割。
针对脑肿瘤的检测和分割任务,研究人员可以构建不同的深度学习网络架构,如U-Net、FCN和SegNet等,来提高准确性和鲁棒性。
在脑肿瘤自动检测和分割方法的研究中,数据集的选择和标注是非常重要的环节。
通常,研究人员需要从医院获取大量的脑肿瘤影像数据,并进行手工标注。
这样的数据集可以帮助研究人员训练和评估他们的方法。
此外,数据增强和数据平衡等技术也被广泛应用于解决数据集不平衡和样本量不足的问题。
尽管基于机器学习和深度学习的方法在脑肿瘤自动检测和分割领域取得了令人瞩目的成就,但仍存在一些挑战和改进空间。
首先,医学影像中的噪声和伪影会对自动检测和分割的结果产生影响,因此需要进一步提高算法的鲁棒性。
第25卷第5期 20]0年9月 数 据 采 集 与 处 理
Journal of Data Acquisition&Processing Vo1.25 No.5
Sep.2010
引 文章编号:1004—9037(2010)05—0605—06 基于特征向量的颅脑CT图像分割方法
孙 涛 李传富h。 刘 伟 冯焕清 (1.中国科学技术大学电子科技系,合肥,230027;2.安徽中医学院第一附属医院影像中心,合肥,230031)
摘要:针对脑肿瘤和脑出血的检出应用,提出一种新的将颅脑CT图像分割为白质、灰质和脑脊液3个区域的方 法。首先用阈值法与数学形态学结合的方法提取出颅腔内脑部组织,然后利用灰度信息将其分割成脑脊液区域 和其他区域。计算其他区域像素的几何矩及几何矩的方差,并结合像素的边界特性为每个像素构建特征向量。用 MFCM算法对特征向量进行聚类将图像分割成白质和灰质,并用数学形态学对聚类结果进行滤波。实验结果验 证了本文方法的先进性。 关键词:图像分割;颅脑CT图像;特征向量;几何矩;模糊C均值聚类 中图分类号:TP391.41 文献标识码:A
Segmentation of Brain CT Images Based on Eigenvectors Sun Tao ,Li Chuanfu ~,Liu Wei ,Feng Huanqing (1.Department of Electronic Science and Technology,University of Science and Technology of China, Hefei,230027,China; 2.Medical Imaging Center,First Affiliated Hospital of Anhui TCM College,Hefei,230031,China)
Abstract:Aimed at the detection of brain tumors and the cerebral hemorrhage,a new approach for tissue segmentation of brain CT data is proposed to label the white matter(WM),the gray matter(GM)and the cerebrospinal fluid(CSF).By using thresholding and mathematical mot— phology,the cerebrum is firstly extracted.After the cerebrum is approximately classified into CSF and the other part according to the discrepancy of their intensity,the geometric moment, its standard deviation,and the edge strength of each voxel of the other part are calculated and constructed as the eigenvector of the voxe1.MFCM classifies the eigenvector into WM and GM,and the mathematical morphology modifies the clustering result.Experimental results demonstrate the good effect of the method. Key words:image segmentation;brain CT images;eigenvector;geometric moments;fuzzy C— means clustering
口 MR图像对软组织的分辨率远比CT图像高, 因而目前颅脑软组织的分割研究主要以MR图像 为对象。但在很多MR图像受限制的应用场合,CT 图像是重要的替代图像。此外,CT图像在钙化点、 脑出血和脑肿瘤等脑部疾病的检测方面具有MR 不可比拟的优势。因而研究基于CT图像的软组织 的分割具有一定的现实意义。 基金项目:国家自然科学基金(60771007)资助项目。 收稿日期:2008—09—28;修订日期:2009-01一l4 目前已有一些成功的颅脑MR图像软组织分 割的方法,但将它们直接用于颅脑CT图像则远达 不到同等的效果H]。也有作者提出过一些针对颅脑 CT的分割方法[2 ],比如阈值分割法、区域生长法 等,这些方法大多集中在研究分割颅骨组织和脑组 织或分割病灶区域等方面,专门针对CT颅脑数据 软组织的分割研究目前还没有见到。针对这个问 题,本文拟对像素计算特征向量,将特征向量作为 分类数据,以实现颅脑CT图像的分割 建立特征 向量的方法已成功用于颅脑MR图像的配准[5刮, 606 数 据 采 集 与 处 理 第2s卷 将它应用于颅脑CT数据的分割是一个尝试。首先 根据灰度信息用改进的模糊均值算法(Modified fuzzy c—mean method,MFCM)_7 将颅脑区域初分 成脑脊液(CSF)区域和其他区域。建立特征向量 后,对其他区域再次进行聚类划分,得到大致的白 质(WM)、灰质(GM),用数学形态学对聚类的结 果加入医学先验知识的约束,修正第二次聚类的结 果,最终实现所期望的分割,将颅腔内组织分割成 WM,GM,CSF3个区域。 1 算法原理 1.1 预处理 本文研究的对象是颅脑的软组织区域,所以应 先去除颅脑CT图像中的颅骨和支架,采用阈值法 与数学形态学结合的方法[8]剔除颅腔外的结构和 提取出颅腔内的部分。图1列出了去除颅骨前后的 CT图像。 圃■ (a)去除颅骨前 (b)去除颅骨后 图l去除颅骨前后的CT图像 1.2特征向量的构造 特征向量构造是最关键的一步,它们必须能反 映不同像素所代表的解剖结构特性。考虑到颅脑解 剖结构的复杂性,若仅以灰度为分割颅脑CT图像 的依据,并不能准确反映颅脑不同组织的差异,不 同组织的灰度值可能在某些区域极为相似,因而造 成分割的错误。图2是图l(b)的灰度分布直方图, 可见,颅腔内的灰质、自质、脑脊液3个部分的CT 值较均匀地分布在一2O~8o之间,且是连续的,无 窨 U CT值/Hu 图2颅腔内结构的灰度直方图 明显的波谷。特别是灰质和白质的CT值有重叠的 区域,只是根据CT值对白质和灰质进行分割容易 导致错误的分类,所以需要针对白质和灰质寻找合 适的分割方法。 克服此问题的~个方法是为每个像素计算其 小波系数[9],但计算代价很大,耗时长。为此,本文 采用从局部直方图提取每个像素的特征的方法。局 部直方图易于计算,且具有旋转不变性和平移不变 性,最重要的是它可以有效反映像素点周围的邻域 信息,从而对不同组织内具有相似灰度值的像素具 有良好的区别能力。 首先计算每个像素在一定邻域范围的局部直 方图h( ),v为像素点的坐标;然后针对h( )计算每 个点的局部直方图的三阶几何矩(Geometric mo— ments,GM)口0]及其在邻域内的方差。三阶几何矩 的计算公式为 Ⅳ ( )一∑ *h( , ) 一l,2,3 (1)
=1 式中: 为模板中心;h(v,i)是V为中心的模板内灰
度i出现的频率(局部直方图h(v)),i是v为中心的 模板内的不同灰度值;N为当前模板下所有不同 灰度的总数。 的取值范围为CT图像的颅脑区域 或者全部区域,这样就为每个像素点建立起各自的 三阶几何矩,即{M (v),M (v),M。(v)}。局部直方 图对于边缘信息的检测能力较弱,用canny边缘检 测算子计算每个像素的canny系数作为特征向量 的一维分量,从而在特征向量中引入像素的边界特 性。canny算子D1]可以根据图像边界强度的大小输 出相应数值的单像素宽的边界,并且能够通过滞后 阈值化处理来消除图像中噪声引起的对单个边缘 的虚假响应造成的“纹状”问题。 如前所述,CT图像对于软组织的成像分辨率 没有MR图像高,为更有效反映颅脑不同结构的差 异,本文在文献EG-I的特征向量形式基础上,增加三 维分量,即计算出每个像素的三阶几何矩后,再计 算其在邻域内各自的标准方差 , 1 、 I s (',)一{‘ ( . ( )一 。(',)) J‘
、‘ t=1
n一1,2,3 (2)
式中:Ⅳ是v为中心的模板下所有像素的个数; . (1,) 为当前模板下第 个像素的第 阶几何矩; (,,)为当 前模板下所有像素的第,z阶几何矩的均值。 构造的特征向量形式为 T一{ 1,M2, 3, 1,S2, 3,C} (3) 式中: ( 一1,2,3)为每个像素在模板邻域内的 608 数 据 采 集 与 处 理 第25卷 割算法。算法假设有限集X---{X , :,…, , }为n个 样本组成的集合, 为预定的类别数, 为每一个 聚类的中心,k=1,2,…,c; ( )为第i个样本关 于第k类的隶属度,则聚类目标函数由隶属函数定 义为
JFc 一∑∑( ) (置, ) (4) =1 1 式中: ( ,V )为 ,到 (是一1,2,…,f)的欧氏距
离;rr/为模糊加权幂指数,可以控制聚类结果的模 糊程度;U为 的模糊划分矩阵,同时满足 (蜀)∈
(0,1), :“ (蕾):1, 为X的聚类中心集合。 ^=1 FCM聚类算法的目的是使聚类目标达到最小。该
算法定义每个像素 对应的中心V 的隶属度 “ ( )为 1 : 一——— — (5)
f堂1\ 』
y 的更新公式为 N ∑u?kx
一 (6) ∑u?k l=1 在本算法中, 是每个点的特征向量 。根据
FCM的迭代算法,必须首先确定FCM初始聚类中 心 。通常聚类中心是随机产生的,如果它与最终 的迭代结果相差较大,那么迭代过程就会很长;若 选择的初始值和最终迭代结果很相近,就能显著减 少FCM的迭代次数。而K—means算法的收敛速度 比FCM算法要快,而且两者的聚类中心十分接 近[7],所以参考改进的MFCM算法,先用K—means 算法计算出待分类数据的硬聚类中心,然后将其作 为FCM算法的初始聚类中心,能减少FCM算法的 迭代次数,提高聚类速度。 文献[13-14]分别从算法的收敛性和聚类有效 性角度着手,提出MFCM算法中的参数m不作特 殊要求的情况下取值为 n一2。 因为要将数据分为wM,GM,CSF和背景4 类,所以c取4。实验中,本算法的主要时间消耗在 建立特征向量的步骤,为加速运算,可以利用脑脊 液的灰度值和颅脑其他区域的灰度差异较为显著