CN111709957A - 一种基于二维最大熵阈值c-v模型的医学图像分割法 - Google Patents
一种基于二维最大熵阈值c-v模型的医学图像分割法 Download PDFInfo
- Publication number
- CN111709957A CN111709957A CN202010574395.1A CN202010574395A CN111709957A CN 111709957 A CN111709957 A CN 111709957A CN 202010574395 A CN202010574395 A CN 202010574395A CN 111709957 A CN111709957 A CN 111709957A
- Authority
- CN
- China
- Prior art keywords
- target
- background
- image
- noise
- quasi
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于二维最大熵阈值C‑V模型的医学图像分割方法,包括:首先,基于输入的原始图像的灰度图,计算出其点灰度‑区域灰度均值的二维直方图;其次,使用二维最大熵法确定最佳阈值,将图像预分割为目标、背景、噪声和边缘四部分;然后,根据预分割的结果定义一个初始的水平集函数,定义该水平集函数为只有两个函数值1和‑1的分段常值函数;最后,通过能量泛函来判断噪声、边缘区域的像素点属于目标或是背景。本发明提供的基于二维最大熵阈值C‑V模型的医学图像分割方法可有效的解决噪声和初始轮廓的位置、形状和大小对曲线演化的效果和速率影响,并通过逐点代入测试的办法来极小化能量泛函,提高了图像分割速度。
Description
技术领域
本发明涉及数字图像处理技术领域,特别是涉及一种基于二维最大熵阈值C-V模型的医学图像分割方法。
背景技术
医学图像的分析与处理是数字图像处理的一个重要应用,它通过计算机软件程序对医疗器械采集得到的各种影像进行分析处理,为医生提供相应辅助而对病人进行精准的诊疗。医学影像分析涉及多种技术,包括图像分割,三维可视化,计算机辅助远程诊疗等等。其中,图像分割是指将图像中具有特殊含义的不同区域区分开来,这些区域是互相不交叉的,每一个区域都满足特定区域的一致性。医学影像图像分割的研究具有重要的意义,结构分析、运动分析、三维可视化等后续操作,以及图像引导手术、肿瘤仿射治疗、治疗评估等应用研究都假设已对图像做了准确分割,或者都是以图像分割为基础的。医学图像分割是正常组织和病变组织的三维重建、定量分析等后续操作的基础,也是临床医学应用的瓶颈。由于医学影像的极其复杂的多样性和复杂性,及目前医学影像设备成像技术上的特点等,都使得医学影像图像存在一定的噪声,图像中目标边缘可能局部不清晰,这就使得医学图像的分割更加困难。
图像分割可用的特征包括图像的灰度、颜色、纹理、局部统计特征和频谱特征等,利用这些特征可以区分图像中不同的目标物体。但目前,在医学影像图像分割方面尚没有通用的理论和方法。因此,有必要针对医学应用这个领域,对图像分割方法进行研究。
当前,C-V模型被广泛用于图像分割的领域。然而,其必须周期性地重新初始化水平集函数,演化速度明显依赖于演化曲线初始位置,因此,还需对其进行进一步改进。
发明内容
本发明的目的是提供一种基于二维最大熵阈值C-V模型的医学图像分割方法,可以有效地解决噪声处理问题和初始位置影响曲线演化速度问题,并通过逐点代入测试的办法来极小化能量泛函,提高了图像分割速度。
为实现上述目的,本发明提供了如下方案:
一种基于二维最大熵阈值C-V模型的医学图像分割方法,该方法包括以下步骤:
S1,对待分割的原始图像进行预分割,采用二维最大熵阈值分割法确定最佳阈值,将原始图像分割为目标、背景、噪声和边缘四区域;
S2,根据预分割结果,定义一个水平集函数,利用水平集函数将图像划分为准目标和准背景两区域;
S3,通过能量泛函F(c1,c2,Φ)和C-V模型来判断噪声和边缘区域中各像素点属于背景还是目标,同时利用同质区域的全局信息对预分割结果进行微调,得到更好的分割结果。
优选的,所述步骤S1中,对待分割的原始图像进行预分割,采用二维最大熵阈值分割法确定最佳阈值,将原始图像分割为目标、背景、噪声和边缘四区域,具体包括:
S101,以原始图像中各像素及其邻域4个像素为一个区域,计算区域灰度均值图像,则原始图像中的每一个像素都对应于一个点灰度-区域灰度均值对,该数据对存在L×L种可能的取值;
S103,建立点灰度-区域灰度均值的二维直方图,在二维直方图的平面图中,用点灰度-区域灰度均值二维最大熵法确定最佳阈值,将原始图像分成目标、背景、噪声和边缘四区域。
优选的,所述步骤S2,根据预分割结果,定义一个水平集函数,利用水平集函数将图像划分为准目标和准背景两区域,具体包括:
水平集函数Φ的表达式如下:
其中,Ω1、Ω2、Ω3、Ω4分别对应目标、噪声、背景和边缘;该水平集函数只有两个函数值1和-1的分段常值函数,其中函数值1对应目标和边缘,函数值-1对应背景和噪声;利用水平集函数将图像划分为准目标和准背景两区域,准目标区域为要确定的目标和边缘,准背景区域为噪声和背景。
优选的,所述步骤S3,通过能量泛函F(c1,c2,Φ)和C-V模型来判断噪声和边缘区域中各像素点属于背景还是目标,同时利用同质区域的全局信息对预分割结果进行微调,得到更好的分割结果,具体包括:
S301,根据水平集函数划分的准目标和准背景两区域,分别计算两区域各自的加权像素平均值,并计算对应的能量泛函;
S302,通过水平集函数Φ寻找集合
S303,对集合K∪Ω2∪Ω4中的像素点,当Φ(x)的符号改变时,如果对应的能量泛函减小,则令Φ(x)=Φ(-x),如果没有减小,则令Φ(x)不变;
S304,不断重复S302至S303步骤,直到能量泛函的值不再变化。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供的基于二维最大熵阈值C-V模型的医学图像分割方法,建立二维直方图,使用二维最大熵阈值分割和C-V模型,不仅充分利用了图像像素点的信息,而且考虑到了像素点与邻域的空间相关信息;加上能量泛函F(c1,c2,Φ)的引入,边缘和噪声区域中各像素点能很好地被判断出来属于背景还是目标,可以有效地解决噪声处理问题和初始位置影响曲线演化速度问题,并通过逐点代入测试的办法来极小化能量泛函,提高了图像分割速度,从而实现快速自动分割医学影像,获得纹理清晰、细节分明的医学影像图。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明基于二维最大熵阈值C-V模型的医学图像分割方法流程图;
图2是二维灰度直方图的平面图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一区域实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于二维最大熵阈值C-V模型的医学图像分割方法,可以有效地解决噪声处理问题和初始位置影响曲线演化速度问题,并通过逐点代入测试的办法来极小化能量泛函,提高了图像分割速度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明提供的基于二维最大熵阈值C-V模型的医学图像分割方法,包括以下步骤:
S1,对待分割的原始图像进行预分割,采用二维最大熵阈值分割法确定最佳阈值,将原始图像分割为目标、背景、噪声和边缘四区域;
S2,根据预分割结果,定义一个水平集函数,利用水平集函数将图像划分为准目标和准背景两区域;
S3,通过能量泛函F(c1,c2,Φ)和c-v模型来判断噪声和边缘区域中各像素点属于背景还是目标,同时利用同质区域的全局信息对预分割结果进行微调,得到更好的分割结果。
其中,所述步骤S1中,对待分割的原始图像进行预分割,采用二维最大熵阈值分割法确定最佳阈值,将原始图像分割为目标、背景、噪声和边缘四区域,具体包括:
S101,以原始图像中各像素及其邻域4个像素为一个区域,计算区域灰度均值图像,则原始图像中的每一个像素都对应于一个点灰度-区域灰度均值对,该数据对存在L×L种可能的取值;
S103,建立点灰度-区域灰度均值的二维直方图,在二维直方图的平面图中,用点灰度-区域灰度均值二维最大熵法确定最佳阈值,将原始图像分成目标、背景、噪声和边缘四区域。
如图2所示,为二维直方图的xoy平面图,沿对角线分布的Ω1区和Ω2区分别代表目标和背景,远离对角线的Ω2区和Ω4区代表噪声和边缘。
设Ω1区和Ω2区各自具有不同的概率分布,用Ω1区和Ω3区的后验概率对各区域的Pi,j进行归一化处理。设阈值在(s,t),则:
Pi,j=fi,j/N
式中,fi,j为图像中灰度为i及其区域灰度均值为j的像素点数,N为图像的总像素数,Pi,j为灰度均值对(i,j)发生的概率,PΩ1为Ω1区的后验概率,PΩ3为Ω3区的后验概率。
Ω1区和Ω3区的二维熵分别为
式中,
熵的判别函数定义为:
通过最佳阈值,划分目标、背景、噪声、边缘,使真正代表目标和背景的信息量最大。
所述步骤S2,根据预分割结果,定义一个水平集函数,利用水平集函数将图像划分为准目标和准背景两区域,具体包括:
水平集函数Φ的表达式如下:
其中,Ω1、Ω2、Ω3、Ω4分别对应目标、噪声、背景和边缘;该水平集函数只有两个函数值1和-1的分段常值函数,其中,函数值1对应目标和边缘,函数值-1对应背景和噪声;利用水平集函数将图像划分为准目标和准背景两区域,准目标区域为要确定的目标和边缘,准背景区域为噪声和背景。通过该水平集函数解决了噪声和初始轮廓的位置、形状和大小对曲线演化效果和速率的影响。
所述步骤S3,通过能量泛函F(c1,c2,Φ)和C-V模型来判断噪声和边缘区域中各像素点属于背景还是目标,同时利用同质区域的全局信息对预分割结果进行微调,得到更好的分割结果,具体包括:
S301,根据水平集函数划分的准目标和准背景两区域,分别计算两区域各自的加权像素平均值,并计算对应的能量泛函;
S302,通过水平集函数Φ寻找集合
S303,对集合K∪Ω2∪Ω4中的像素点,当Φ(x)的符号改变时,如果对应的能量泛函减小,则令Φ(x)=Φ(-x),如果没有减小,则令Φ(x)不变;
S304,不断重复S302至S303步骤,直到能量泛函的值不再变化。
其中,步骤3中的c-v模型,使用水平集函数表示目标轮廓,并一般通过如下能量泛函求极小值对图像分割:
能量泛函F(c1,c2,Φ)为
式中,δ为单位冲击函数,Φ为水平集函数,H为Heaviside函数;μ、λ1、λ2为加权系数,要求μ≥0,λ1>0,λ2>0,c1、c2分别为演化曲线内部和外部区域灰度均值,u0(x,y)为图像的水平集函数形式,其函数表示为:
u0(x,y)=c1H(φ(x,y))+c2(1-H(φ(x,y)))
Heaviside函数表示为:
利用Heaviside函数,令Φ不变,最小化能量函数F(C1,C2,Φ),对当前零水平集内部和外部的图像灰度平均值C1和C2计算如下:
本发明提供的基于二维最大熵阈值C-V模型的医学图像分割方法,建立二维直方图,使用二维最大熵阈值分割和C-V模型,不仅充分利用了图像像素点的信息,而且考虑到了像素点与邻域的空间相关信息;加上能量泛函F(c1,c2,Φ)的引入,边缘和噪声区域中各像素点能很好地被判断出来属于背景还是目标,可以有效地解决噪声处理问题和初始位置影响曲线演化速度问题,并通过逐点代入测试的办法来极小化能量泛函,提高了图像分割速度,从而实现快速自动分割医学影像,获得纹理清晰、细节分明的医学影像图。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种基于二维最大熵阈值C-V模型的医学图像分割方法,其特征在于,包括以下步骤:
S1,对待分割的原始图像进行预分割,采用二维最大熵阈值分割法确定最佳阈值,将原始图像分割为目标、背景、噪声和边缘四区域;
S2,根据预分割结果,定义一个水平集函数,利用水平集函数将图像划分为准目标和准背景两区域;
S3,通过能量泛函F(c1,c2,Φ)和C-V模型来判断噪声和边缘区域中各像素点属于背景还是目标,同时利用同质区域的全局信息对预分割结果进行微调,得到更好的分割结果。
2.根据权利要求1所述的基于二维最大熵阈值C-V模型的医学图像分割方法,其特征在于,所述步骤S1中,对待分割的原始图像进行预分割,采用二维最大熵阈值分割法确定最佳阈值,将原始图像分割为目标、背景、噪声和边缘四区域,具体包括:
S101,以原始图像中各像素及其邻域4个像素为一个区域,计算区域灰度均值图像,则原始图像中的每一个像素都对应于一个点灰度-区域灰度均值对,该数据对存在L×L种可能的取值;
S103,建立点灰度-区域灰度均值的二维直方图,在二维直方图的平面图中,用点灰度-区域灰度均值二维最大熵法确定最佳阈值,将原始图像分成目标、背景、噪声和边缘四区域。
4.根据权利要求3所述的基于二维最大熵阈值C-V模型的医学图像分割方法,其特征在于,所述步骤S3,通过能量泛函F(c1,c2,Φ)和C-V模型来判断噪声和边缘区域中各像素点属于背景还是目标,同时利用同质区域的全局信息对预分割结果进行微调,得到更好的分割结果,具体包括:
S301,根据水平集函数划分的准目标和准背景两区域,分别计算两区域各自的加权像素平均值,并计算对应的能量泛函;
S303,对集合K∪Ω2∪Ω4中的像素点,当Φ(x)的符号改变时,如果对应的能量泛函减小,则令Φ(x)=Φ(-x),如果没有减小,则令Φ(x)不变;
S304,不断重复S302至S303步骤,直到能量泛函的值不再变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010574395.1A CN111709957A (zh) | 2020-06-22 | 2020-06-22 | 一种基于二维最大熵阈值c-v模型的医学图像分割法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010574395.1A CN111709957A (zh) | 2020-06-22 | 2020-06-22 | 一种基于二维最大熵阈值c-v模型的医学图像分割法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111709957A true CN111709957A (zh) | 2020-09-25 |
Family
ID=72541544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010574395.1A Pending CN111709957A (zh) | 2020-06-22 | 2020-06-22 | 一种基于二维最大熵阈值c-v模型的医学图像分割法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111709957A (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006126347A1 (ja) * | 2005-05-27 | 2006-11-30 | Konica Minolta Photo Imaging, Inc. | 撮像装置、画像処理方法及び画像処理プログラム |
CN102147922A (zh) * | 2011-05-05 | 2011-08-10 | 河南工业大学 | 灰度图像二维Otsu折线阈值分割法 |
US20120075440A1 (en) * | 2010-09-28 | 2012-03-29 | Qualcomm Incorporated | Entropy based image separation |
CN102831604A (zh) * | 2012-07-30 | 2012-12-19 | 常州大学 | 一种灰度图像的二维Renyi熵阈值分割方法 |
WO2016101279A1 (zh) * | 2014-12-26 | 2016-06-30 | 中国海洋大学 | 一种合成孔径雷达图像舰船目标快速检测方法 |
CN106023205A (zh) * | 2016-05-23 | 2016-10-12 | 杭州健培科技有限公司 | 基于简化粒子群优化二维最大熵阈值的医学影像分割方法 |
CN107292890A (zh) * | 2017-06-19 | 2017-10-24 | 北京理工大学 | 一种医学图像分割方法和装置 |
CN108830864A (zh) * | 2018-05-03 | 2018-11-16 | 上海电机学院 | 图像分割方法 |
CN110378924A (zh) * | 2019-07-19 | 2019-10-25 | 合肥学院 | 基于局部熵的水平集图像分割方法 |
-
2020
- 2020-06-22 CN CN202010574395.1A patent/CN111709957A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006126347A1 (ja) * | 2005-05-27 | 2006-11-30 | Konica Minolta Photo Imaging, Inc. | 撮像装置、画像処理方法及び画像処理プログラム |
US20120075440A1 (en) * | 2010-09-28 | 2012-03-29 | Qualcomm Incorporated | Entropy based image separation |
CN102147922A (zh) * | 2011-05-05 | 2011-08-10 | 河南工业大学 | 灰度图像二维Otsu折线阈值分割法 |
CN102831604A (zh) * | 2012-07-30 | 2012-12-19 | 常州大学 | 一种灰度图像的二维Renyi熵阈值分割方法 |
WO2016101279A1 (zh) * | 2014-12-26 | 2016-06-30 | 中国海洋大学 | 一种合成孔径雷达图像舰船目标快速检测方法 |
CN106023205A (zh) * | 2016-05-23 | 2016-10-12 | 杭州健培科技有限公司 | 基于简化粒子群优化二维最大熵阈值的医学影像分割方法 |
CN107292890A (zh) * | 2017-06-19 | 2017-10-24 | 北京理工大学 | 一种医学图像分割方法和装置 |
CN108830864A (zh) * | 2018-05-03 | 2018-11-16 | 上海电机学院 | 图像分割方法 |
CN110378924A (zh) * | 2019-07-19 | 2019-10-25 | 合肥学院 | 基于局部熵的水平集图像分割方法 |
Non-Patent Citations (3)
Title |
---|
刘正光等: "基于二维灰度直方图的模糊熵分割方法", 《天津大学学报》 * |
杨帆: "《数字图像处理与分析 第4版》", 北京航空航天大学出版社 * |
职占江等: "基于Chan-Vese模型的医学图像分割算法", 《计算机工程与应用》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210192758A1 (en) | Image processing method and apparatus, electronic device, and computer readable storage medium | |
CN110930416B (zh) | 一种基于u型网络的mri图像前列腺分割方法 | |
CN110321920A (zh) | 图像分类方法、装置、计算机可读存储介质和计算机设备 | |
TW202044198A (zh) | 一種圖像處理方法、裝置、電子設備及電腦可讀儲存媒介 | |
US10405834B2 (en) | Surface modeling of a segmented echogenic structure for detection and measurement of anatomical anomalies | |
Li et al. | Automated measurement network for accurate segmentation and parameter modification in fetal head ultrasound images | |
WO2021136368A1 (zh) | 钼靶图像中胸大肌区域自动检测方法及装置 | |
Ballester et al. | Segmentation and measurement of brain structures in MRI including confidence bounds | |
Song et al. | Kidney segmentation in CT sequences using SKFCM and improved GrowCut algorithm | |
Heyde et al. | Anatomical image registration using volume conservation to assess cardiac deformation from 3D ultrasound recordings | |
Wei et al. | Learning-based 3D surface optimization from medical image reconstruction | |
Göçeri et al. | A comparative performance evaluation of various approaches for liver segmentation from SPIR images | |
CN115830016B (zh) | 医学图像配准模型训练方法及设备 | |
CN111383759A (zh) | 一种肺炎自动诊断系统 | |
Zhao et al. | Segmentation of dermoscopy images based on deformable 3D convolution and ResU-NeXt++ | |
Goyal et al. | Multi-modality image fusion for medical assistive technology management based on hybrid domain filtering | |
Zhang et al. | Segmentation of brain tumor MRI image based on improved attention module Unet network | |
Zhu et al. | MRI enhancement based on visual-attention by adaptive contrast adjustment and image fusion | |
CN117911432A (zh) | 图像分割方法、装置及存储介质 | |
CN112750110A (zh) | 基于神经网络对肺部病灶区进行评估的评估系统和相关产品 | |
CN117115437A (zh) | 一种基于区域的多指标多器官医学图像分割模型评估系统 | |
Xu et al. | Deep convolutional neural networks for bias field correction of brain magnetic resonance images | |
Dong et al. | Hole-filling based on content loss indexed 3D partial convolution network for freehand ultrasound reconstruction | |
CN115439423B (zh) | 一种基于ct图像的识别方法、装置、设备及存储介质 | |
CN108447066B (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 |