一种混合滤波策略和模糊C均值的医学图像分割方法技术

本发明专利技术涉及一种混合滤波策略和模糊C均值的医学图像分割方法，包括以下步骤：读取医学图像，并对医学图像进行中值滤波；通过WFT将医学图像变换到频域；选择阈值函数进行噪声滤除；对滤波后的图像进行WFT逆变换；初始化聚类中心和样本点到聚类中心的距离，然后更新隶属度矩阵和聚类中心；计算目标函数值，若达到最大迭代次数或前后两次目标函数值的绝对值误差小于目标函数精度，则停止计算，否则继续迭代；输出聚类中心和各类隶属度值并完成模糊聚类的划分，输出分割后的图像。本发明专利技术利用直觉模糊集信息，达到了增强抗噪性、提高分割精确度的效果，解决了对医学图像进行高精度分割的实际问题。

A Hybrid Filtering Strategy and Fuzzy C-Mean Medical Image Segmentation Method

【技术实现步骤摘要】
一种混合滤波策略和模糊C均值的医学图像分割方法
本专利技术涉及数字图像处理
，特别是涉及一种混合滤波策略和模糊C均值的医学图像分割方法。
技术介绍
近年来，随着医学影像设备的迅猛发展和普及，全世界的医疗机构每天都要产生海量的医学影像资料，使得医学影像成为医学技术当中发展最快的领域之一。成像技术已经成为一种无侵害观察器官组织解剖结构的手段，例如计算机断层扫描(CT)，磁共振成像(MRI)等。新型成像技术以及设备的出现和计算机技术的发展，使得医学图像处理技术对医学科研和临床实践的作用、影响日益增大。医学图像处理的结果使科研人员和临床工作者对人体内部正常和病变位置的观察更加直观、更加清晰，提高了确诊率，为临床诊断和生物医学研究提供了强有力的保证。因此，国内外相关专家一直都非常重视医学图像处理技术，而医学图像分割是图像处理技术中的一个经典难题，在病理分析和手术计划方面，医学图像分割可以帮助实现医学图像的三维重建和可视化，这样有利于医生根据病人的三维信息，提前制定更好的手术方案，另外在影像信息处理方面，一个病人做一次MR或CT，可以生成大量的医学图像，会占据很大的存储空间，然而需要的仅仅是病人相关的病灶信息，所以医学图像分割有助于对病人关键信息的压缩，这在远程医疗中实现高效的医学图像信息的传输中有着很大的作用。医学图像分割在医学领域的各方面都起着重要的作用，有助于医学领域更进一步的发展，在医学研究方面，医学图像分割可以精确地分割目标并提取，这样有助有后续医生对这些医学图像进一步的分析和处理；在临床诊断方面，医学图像分割有助于医生对人体部分器官和病灶区域的大小的一个测量，这样有助于医生对这个病人的诊断和后续医疗方案的制定。以各种理论为基础，经过几十年的探索研究，人们提出了许多图像分割算法。图像分割虽然已经有了一定的研究成果，但提出的分割方法往往是针对某个特定的领域，适用面有限，没有一种通用的分割算法可以分割任意图像，通用的分割理论尚不能指定。
技术实现思路
本专利技术所要解决的技术问题是提供一种混合滤波策略和模糊C均值的医学图像分割方法，能够达到增强抗噪性、提高分割精确度的效果。本专利技术解决其技术问题所采用的技术方案是：提供一种混合滤波策略和模糊C均值的医学图像分割方法，包括以下步骤：(1)读取医学图像I(x,y)，并对医学图像I(x,y)进行中值滤波；(2)通过WFT将医学图像I(x,y)变换到频域；(3)选择阈值函数进行噪声滤除；(4)对滤除噪声后的图像进行WFT逆变换；(5)初始化聚类中心P和样本点到聚类中心的距离，然后根据更新隶属度矩阵U和聚类中心，其中，隶属度矩阵U为一个c×n的二维矩阵，隶属度矩阵U中的元素uij表示第j个数据点pj属于聚类i的隶属度，满足约束条件vi为聚类中心，||·||表示欧氏距离；m为模糊加权指数；(6)计算目标函数值，若达到最大迭代次数或前后两次目标函数值的绝对值误差小于目标函数精度，则停止计算，否则返回步骤(5)；(7)输出聚类中心和各类隶属度值并完成模糊聚类的划分，输出分割后的图像。所述步骤(1)中的中值滤波具体为：将模板中的各个像素点的灰度值进行排序，取中间值作为这个模板中心像素点的灰度值，使得医学图像中噪声呈非平稳高斯分布。所述步骤(2)中通过WFT将医学图像I(x,y)变换到频域，具体变换方式为：其中，(u,v)是时域坐标，(ξ,η)是频域坐标，SI(u,v；ξ,η)是WFT的谱函数，g(x,y)是窗函数，σx和σy分别是x轴和y轴的标准差。所述步骤(3)中的阈值函数为：其中，上横线表示滤波过程，SI(u,v；ξ,η)是WFT的谱函数，thr是阈值。所述步骤(4)中对滤除噪声后的图像进行WFT逆变换，具体的逆变换方式为：其中，为滤波后的医学图像。所述步骤(6)中通过计算目标函数值。有益效果由于采用了上述的技术方案，本专利技术与现有技术相比，具有以下的优点和积极效果：本专利技术采用组合滤波方法，首先采用中值滤波滤除孤立噪声点或脉冲知识点，然后进行窗口傅里叶滤波，最后结合利用了直觉模糊集理论，达到了增强抗噪性、提高分割精确度的效果，解决了对医学图像进行高精度分割的实际问题。本专利技术提出的改进算法具有较高的应用价值和广阔的市场前景。附图说明图1是本专利技术的流程图。具体实施方式下面结合具体实施例，进一步阐述本专利技术。应理解，这些实施例仅用于说明本专利技术而不用于限制本专利技术的范围。此外应理解，在阅读了本专利技术讲授的内容之后，本领域技术人员可以对本专利技术作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。本专利技术的实施方式涉及一种混合滤波策略和模糊C均值的医学图像分割方法，如图1所示，包括以下步骤：(1)给定聚类数目c，目标函数精度ε，模糊指数m，最大迭代次数maxTimes，以及人工需设定的参数；(2)读取医学图像I(x,y)，并对医学图像I(x,y)进行中值滤波；将模板中的各个像素点的灰度值进行排序，取中间值作为这个模板中心像素点的灰度值，中值滤波能够抑制图像中的颗粒噪声，保持了图像边缘及细节信息，提高了图像的质量。在处理过程中，模板窗口的大小很重要，一般选取3×3。(3)通过WFT将医学图像I(x,y)变换到频域；上式中各参数的含义如下：(u,v)是时域坐标，(ξ,η)是频域坐标，SI(u,v；ξ,η)是WFT的谱函数，g(x,y)是窗函数，σx和σy分别是x轴和y轴的标准差，高斯窗尺寸σ的选择时，大窗有利于分离低频分量和噪声的滤除，小窗则有着计算误差小的优点。(4)选择阈值函数用于噪声滤除；其中，上横线表示滤波过程，SI(u,v；ξ,η)是WFT的谱函数，thr是阈值。(5)对滤波后的图像进行WFT逆变换；即为滤波后的医学图像。(6)初始化聚类中心P和样本点到聚类中心的距离，然后根据以下的公式更新隶属度矩阵U和聚类中心；||·||表示欧氏距离；m为模糊加权指数，且m∈[1,∞)；定义一个c×n的二维隶属矩阵U，U中的元素uij表示第j个数据点pj属于聚类i的隶属度，并满足约束条件vi为聚类中心。(7)计算目标函数值若达到最大迭代次数或前后两次J的绝对值误差小于ε，则停止计算，否则转到步骤(6)；(8)输出聚类中心和各类隶属度值并完成模糊聚类的划分，输出分割后的图像。不难发现，相比于传统医学图像分割方法，本专利技术采用组合滤波方法，首先采用中值滤波滤除孤立噪声点或脉冲知识点，然后进行窗口傅里叶滤波，最后结合利用了直觉模糊集理论，提出了一种混合滤波策略和模糊聚类算法，达到了增强抗噪性、提高分割精确度的效果，解决了对医学图像进行高精度分割的实际问题。本专利技术提出的改进算法具有较高的应用价值和广阔的市场前景。本文档来自技高网...

【技术保护点】
1.一种混合滤波策略和模糊C均值的医学图像分割方法，其特征在于，包括以下步骤：(1)读取医学图像I(x,y)，并对医学图像I(x,y)进行中值滤波；(2)通过WFT将医学图像I(x,y)变换到频域；(3)选择阈值函数进行噪声滤除；(4)对滤除噪声后的图像进行WFT逆变换；(5)初始化聚类中心和样本点到聚类中心的距离，然后根据

【技术特征摘要】
1.一种混合滤波策略和模糊C均值的医学图像分割方法，其特征在于，包括以下步骤：(1)读取医学图像I(x,y)，并对医学图像I(x,y)进行中值滤波；(2)通过WFT将医学图像I(x,y)变换到频域；(3)选择阈值函数进行噪声滤除；(4)对滤除噪声后的图像进行WFT逆变换；(5)初始化聚类中心和样本点到聚类中心的距离，然后根据更新隶属度矩阵U和聚类中心，其中，隶属度矩阵U为一个c×n的二维矩阵，隶属度矩阵U中的元素uij表示第j个数据点pj属于聚类i的隶属度，满足约束条件vi为聚类中心，||·||表示欧氏距离；m为模糊加权指数；(6)计算目标函数值，若达到最大迭代次数或前后两次目标函数值的绝对值误差小于目标函数精度，则停止计算，否则返回步骤(5)；(7)输出聚类中心和各类隶属度值并完成模糊聚类的划分，输出分割后的图像。2.根据权利要求1所述的混合滤波策略和模糊C均值的医学图像分割方法，其特征在于，所述步骤(1)中的中值滤波具体为：将模板中的各个像素点的灰度值进...

【专利技术属性】
技术研发人员：李孝斐，赵曙光，周玲，
申请(专利权)人：东华大学，
类型：发明
国别省市：上海,31