CN116721099B - 一种基于聚类的肝脏ct影像的图像分割方法 - Google Patents
一种基于聚类的肝脏ct影像的图像分割方法 Download PDFInfo
- Publication number
- CN116721099B CN116721099B CN202310993928.3A CN202310993928A CN116721099B CN 116721099 B CN116721099 B CN 116721099B CN 202310993928 A CN202310993928 A CN 202310993928A CN 116721099 B CN116721099 B CN 116721099B
- Authority
- CN
- China
- Prior art keywords
- gray value
- pixel points
- value interval
- gray
- clustering
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 210000004185 liver Anatomy 0.000 title claims abstract description 26
- 238000003709 image segmentation Methods 0.000 title claims abstract description 14
- 230000002776 aggregation Effects 0.000 claims abstract description 26
- 238000004220 aggregation Methods 0.000 claims abstract description 26
- 210000001015 abdomen Anatomy 0.000 claims abstract description 16
- 230000008859 change Effects 0.000 claims description 31
- 238000009825 accumulation Methods 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 claims description 3
- 238000011946 reduction process Methods 0.000 claims description 2
- 238000003064 k means clustering Methods 0.000 abstract description 15
- 230000011218 segmentation Effects 0.000 abstract description 15
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 4
- 230000009467 reduction Effects 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003703 image analysis method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/762—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using clustering, e.g. of similar faces in social networks
-
- 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/10081—Computed x-ray tomography [CT]
-
- 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/30056—Liver; Hepatic
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Artificial Intelligence (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于聚类的肝脏CT影像的图像分割方法,涉及医疗图像处理技术领域,包括获取腹部CT原始图像;获取CT原始图像的灰度直方图和灰度梯度直方图;将CT原始图像中的像素点灰度值分为多个灰度值区间;获取CT原始图像中的边缘特征像素点;获取每个灰度值区间的聚类匹配程度;获取每个灰度值区间的像素点聚集度;获取每个目标灰度值区间的聚类起始像素点;根据CT原始图像中每个目标灰度值区间中的聚类起始像素点对CT原始图像中的像素点进行聚类,完成对CT原始图像的分割;本发明解决了相关技术中,在使用K‑means聚类法对肝脏CT图像进行分割时分割效果差,以及计算冗余度大的技术问题。
Description
技术领域
本发明涉及医疗图像处理技术领域,具体涉及一种基于聚类的肝脏CT影像的图像分割方法。
背景技术
腹部CT图像则是目前医疗人员获取肝脏状态的最重要的手段之一,通过对腹部CT图像中肝脏区域的分割提取能够用于医学教学等多种领域。因此,从腹部CT图像中完整、精确的进行肝脏的分割,就显得尤为重要。但是,肝脏具有与周围毗邻器官的低对比度、弱边界以及病例间肝脏形态差异较大等特点,腹部CT图像中肝脏的精确分割仍面临较大困难。
现有的基于腹部CT序列图像的肝脏分割方法包括基于传统图像分析的方法、基于机器学习的方法和K-means聚类方法。基于传统图像分析的方法,分割流程复杂、耗时长、且需要人工干预,难以实现全自动分割,而且分割结果不理想。基于机器学习的方法,对判别性特征的设计要求非常高,而医学图像的复杂性高,判别特征设计难以得到理想结果;现有的K-means聚类方法在对肝脏CT图像进行分割时,由于类中心点处于一个均衡区域,会导致图像分割效果没有发生变化,同时在的K-means聚类迭代过程中,类中心点在灰度轴中的移动是无限制的,由此会增大K-means聚类过程中的计算冗余度。
发明内容
为了解决相关技术中,在使用K-means聚类法对肝脏CT图像进行分割时分割效果无法满足要求,以及计算冗余度大的技术问题,本发明提供了一种基于聚类的肝脏CT影像的图像分割方法,该方法通过获取CT图像的灰度直方图和灰度梯度直方图确定多个灰度值区间,并通过每个灰度值区间的聚类匹配程度和像素点聚集度确定多个目标灰度值区间,获取每个目标灰度值区间的初始聚类像素点,并根据每个目标灰度值区间的初始聚类像素点对CT图像中的像素点进行K-means聚类,完成对肝脏CT图像的分割;有鉴于此,本发明通过以下技术方案予以实现。
一种基于聚类的肝脏CT影像的图像分割方法,包括以下步骤:
获取腹部CT原始图像;
获取所述CT原始图像的灰度直方图和灰度梯度直方图;
在所述灰度梯度直方图中从最大灰度梯度值开始,向灰度梯度值减小的方向逐级进行像素点数量累加,每累加一次后获取累加像素点数量与CT原始图像中像素点数量的比值,根据当前累加后的比值与上一次累加后的比值获取比值的变化率,根据所述比值的变化率获取所述CT原始图像中的边缘特征像素点;
根据所述灰度直方图中像素点的灰度级分布将所述CT原始图像中的灰度值分为多个灰度值区间;根据每个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值,获取每个灰度值区间的聚类匹配程度;
根据每个灰度值区间中每个像素点与该灰度值区间中其余像素点在所述CT原始图像中的距离的平均值,获取每个灰度值区间中的像素点的聚集度;
根据灰度值区间的聚类匹配程度和聚集度从多个灰度值区间中选取目标灰度值区间;
获取目标灰度值区间的像素点形成的区域中每个像素点与该区域中其余像素点的灰度值的差值平方和;根据差值平方和获取该区域中的聚类起始像素点;依次获取每个目标灰度值区间的像素点形成区域的聚类起始像素点;
根据所述聚类起始像素点对每个目标灰度值区间的像素点形成区域进行聚类,获得多个聚类区域,得到分割后的CT影像图像。
进一步地,根据所述比值的变化率获取所述CT原始图像中的边缘特征像素点的过程中,还包括设定变化率阈值,当所述比值的变化率大于所述变化率阈值时,停止像素点数量的累加,并将已经累加的像素点确定为所述CT原始图像中的边缘特征像素点。
进一步地,所述变化率阈值为0.3,当所述比值的变化率大于0.3时,停止像素点数量的累加,并将已经累加的像素点确定为所述CT原始图像中的边缘特征像素点。
进一步地,所述灰度值区间的聚类匹配程度通过下式确定:
式中,为第/>个灰度值区间的聚类匹配程度;/>为第/>个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值;/>为第/>个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值;/>为灰度值区间个数;/>表示以自然常数e为底的指数函数。
进一步地,所述灰度值区间中的像素点的聚集度通过所述距离的平均值与该灰度值区间中全部像素点灰度值和的比值得到。
进一步地,所述获取每个目标灰度值区间的像素点形成区域的聚类起始像素点的过程为:
获取所述目标灰度值区间的像素点形成的区域中每个像素点与该区域中其余像素点灰度值的差值平方和;获取所述差值平方和中的最大值,并将所述差值平方和中的最大值对应的像素点作为该区域的聚类起始像素点;依次获取每个目标灰度值区间的像素点形成的每个区域的聚类起始像素点。
进一步地,所述根据灰度值区间的聚类匹配程度和聚集度从多个灰度值区间中选取目标灰度值区间的过程中,还包括:
根据每个灰度值区间的聚类匹配程度从多个灰度值区间中选取非必要灰度值区间;
将所述非必要灰度值区间相邻的灰度值区间的范围进行调整,获取调整后的灰度值区间的聚类匹配程度和聚集度;当调整后的灰度值区间的聚类匹配程度和聚集度均不再增大时,停止对灰度值区间的调整,并将剩余的全部灰度值区间确定为目标灰度值区间。
进一步地,所述灰度值区间的聚类匹配程度小于0.3时,该灰度值区间为非必要灰度值区间。
进一步地,所述获取腹部CT原始图像后,还包括对所述CT原始图像进行降噪处理。
进一步地,所述降噪处理采用高斯滤波降噪。
与现有技术相比,本发明的有益效果是:
本发明提供了一种基于聚类的肝脏CT影像的图像分割方法,该方法通过获取腹部CT原始图像;获取CT原始图像的灰度直方图和灰度梯度直方图;由此可根据灰度直方图中灰度值分布确定多个灰度值区间,灰度梯度直方图用于获取CT原始图像中的边缘特征像素点;在灰度直方图中从最大灰度梯度值开始,向灰度梯度值减小的方向逐级进行像素点数量累加,每累加一次后获取累加像素点数量与CT原始图像中像素点数量的比值,根据当前累加后的比值与上一次累加后的比值获取比值的变化率,根据比值的变化率获取CT原始图像中的边缘特征像素点;本发明之所以根据当前累加后的比值与上一次累加后的比值获取比值的变化率获取边缘特征像素点,是因为在腹部CT原始图像中边缘像素点灰度值与其他区域的灰度值特征较大,当比值的变化率较大时,说明边缘像素点与其余区域有明显的分界线;根据灰度直方图中像素点的灰度级分布将CT原始图像中的灰度值分为多个灰度值区间;根据每个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值,获取每个灰度值区间的聚类匹配程度;根据每个灰度值区间中每个像素点与该灰度值区间中其余像素点在CT原始图像中的距离的平均值,获取每个灰度值区间中的像素点的聚集度;获取每个灰度值区间的像素点聚集度是为了确定每个灰度值区间中的起始聚类像素点;在获得每个灰度值区间中的起始聚类像素点后可对CT原始图像中每个目标灰度值区间的像素点进行K-means聚类,可实现对CT原始图像中像素点的分类和对CT原始图像的分割;通过本发明的技术方案实现了对CT原始图像中每个目标灰度值区间的像素点进行K-means聚类,解决了相关技术中,在使用K-means聚类法对肝脏CT图像进行分割时分割效果差,以及计算冗余度大的技术问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为一种基于聚类的肝脏CT影像的图像分割方法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例提供了一种基于聚类的肝脏CT影像的图像分割方法,如图1所示,包括:
S101,获取腹部CT原始图像;获取CT原始图像的灰度直方图和灰度梯度直方图;
本实施例中,根据腹部CT原始图像获得了灰度直方图和灰度梯度直方图;在灰度直方图中,横轴是灰度值,纵轴是CT原始图像中每个灰度值对应的像素点数量;在灰度梯度直方图中,横轴是灰度梯度值,纵轴是每个灰度梯度值对应的像素点数量;获取腹部CT原始图像后,还包括对所述CT原始图像进行降噪处理;在本发明一个实施例中选择高斯滤波降噪。
S102,获取CT原始图像中的边缘特征像素点;
在所述灰度梯度直方图中从最大灰度梯度值开始,向灰度梯度值减小的方向逐级进行像素点数量累加,每累加一次后获取累加像素点数量与CT原始图像中像素点数量的比值,根据当前累加后的比值与上一次累加后的比值获取比值的变化率,根据所述比值的变化率获取所述CT原始图像中的边缘特征像素点;
需要说明的是,本实施例中根据灰度直方图中像素点灰度值分布将CT原始图像中的像素点灰度值分为多个灰度值区间;获取根据灰度直方图中像素点灰度值分布将CT原始图像中的像素点灰度值分为多个灰度值区间的过程为,将灰度直方图拟合成灰度曲线,获取灰度曲线中的极小值点,根据极小值点可将灰度曲线划分为多个灰度值区间;在灰度梯度直方图中从最大灰度级开始,向灰度级减小的方向逐级进行像素点数量累加,每累加一次获取累加值与CT原始图像中像素点数量的比值,比较相邻像素点数量比值的变化率,设定变化率阈值,当相邻像素点数量比值的变化率大于变化率阈值时,停止像素点数量累加,并将已经进行累加的像素点确定为边缘特征像素点;
需要说明的是,本实施例中根据具体实施条件设定变化率阈值为0.3,当相邻像素点数量比值的变化率大于0.3时,停止像素点数量累加,并将已经进行累加的像素点确定为边缘特征像素点,实施者可根据具体实施条件设定其它的数值作为变化率阈值;
S103,将CT原始图像中的灰度值分为多个灰度值区间,并获取目标灰度值区间,获取每个目标灰度值区间的聚类起始像素点;
本实施例中,根据所述灰度直方图中像素点的灰度级分布将所述CT原始图像中的灰度值分为多个灰度值区间;根据每个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值,获取每个灰度值区间的聚类匹配程度;根据每个灰度值区间中每个像素点与该灰度值区间中其余像素点在所述CT原始图像中的距离的平均值,获取每个灰度值区间中的像素点的聚集度;根据灰度值区间的聚类匹配程度和聚集度从多个灰度值区间中选取目标灰度值区间;获取目标灰度值区间的像素点形成的区域中每个像素点与该区域中其余像素点的灰度值的差值平方和;根据差值平方和获取该区域中的聚类起始像素点;依次获取每个目标灰度值区间的像素点形成区域的聚类起始像素点;
需要说明的是,本实施例在中获取每个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值;根据每个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值获取每个灰度值区间的聚类匹配程度;灰度值区间的聚类匹配程度通过下式确定:
式中,为第/>个灰度值区间的聚类匹配程度;/>为第/>个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值;/>为第/>个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值;/>为灰度值区间个数;/>表示以自然常数e为底的指数函数;
获取每个灰度值区间中单个像素点与其余像素点在CT原始图像中的平均距离;根据单个像素点与其余像素点在CT原始图像中的平均距离获取每个灰度值区间的像素点聚集度;灰度值区间的像素点聚集度为该灰度值区间中单个像素点与其余像素点在CT原始图像中的平均距离之和,与该灰度值区间中全部像素点灰度值之和的比值;灰度值区间的像素点聚集度通过下式确定:
式中,为第/>个灰度值区间的像素点聚集度;/>为CT原始图像中第/>个灰度值区间的像素点/>与第/>个像素点/>之间的距离;/>为第/>个灰度值区间中像素点个数;/>为第/>个灰度值区间中第/>个像素点的灰度值;
本实施例中,根据灰度值区间的聚类匹配程度和像素点聚集度从多个灰度值区间中获取多个目标灰度值区间;根据灰度值区间的聚类匹配程度和像素点聚集度从多个灰度值区间中获取多个目标灰度值区间的过程为,根据灰度值区间的聚类匹配程度获取非必要灰度值区间;将与非必要灰度值区间相邻的灰度值区间范围进行调整,获取调整后的灰度值区间的聚类匹配程度和聚集度均不再增大时,停止对灰度值区间的调整,获得多个目标灰度值区间;目标灰度值区间的像素点形成的区域中每个像素点与该区域中其余像素点的灰度值的差值平方和通过下式确定:
式中,为目标灰度值区间的像素点形成的区域中第/>个像素点与该区域中其余像素点灰度值的差值平方和;/>为目标灰度值区间的像素点形成的区域中第/>个像素点的灰度值;/>为目标灰度值区间的像素点形成的区域中第/>个像素点的灰度值;/>为目标灰度值区间的像素点形成的区域中像素点的个数;
S104,根据所述聚类起始像素点对每个目标灰度值区间的像素点形成区域进行聚类,获得多个聚类区域,得到分割后的CT原始图像;
本实施例中,根据CT原始图像中每个目标灰度值区间中的聚类起始像素点对CT原始图像中的像素点进行K-means聚类,得到多个聚类区域,完成对CT原始图像的分割;需要说明的是,在对CT原始图像中的像素点进行K-means聚类后,可对每个聚类区域用不同的颜色进行标记,以增强CT原始图像中每个聚簇区域的区分度;对于每个目标灰度值区间聚类起始像素点的选取也就是K-means聚类中k值的选取;对于目标灰度值区间的选取是通过对非必要灰度值区间的相邻的灰度值区间的调整来实现的;相较于传统的聚类中心点选取方法,如手肘法,轮廓系数法,本实施例的结合CT原始图像中像素点的灰度值和灰度梯度值分布获得的目标区间更适合于在单一灰度轴上找到最优的聚类中心点,本实施例的方法更符合CT原始图像中像素点的分布特点。
综上,本实施例提供了一种基于聚类的肝脏CT影像的图像分割方法,该方法通过获取腹部CT原始图像;获取CT原始图像的灰度直方图和灰度梯度直方图;由此可根据灰度直方图中灰度值分布确定多个灰度值区间,灰度梯度直方图用于获取CT原始图像中的边缘特征像素点;在灰度直方图中从最大灰度梯度值开始,向灰度梯度值减小的方向逐级进行像素点数量累加,每累加一次后获取累加像素点数量与CT原始图像中像素点数量的比值,根据当前累加后的比值与上一次累加后的比值获取比值的变化率,根据比值的变化率获取CT原始图像中的边缘特征像素点;本发明之所以根据当前累加后的比值与上一次累加后的比值获取比值的变化率获取边缘特征像素点,是因为在腹部CT原始图像中边缘像素点灰度值与其他区域的灰度值特征较大,当比值的变化率较大时,说明边缘像素点与其余区域有明显的分界线;根据灰度直方图中像素点的灰度级分布将CT原始图像中的灰度值分为多个灰度值区间;根据每个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值,获取每个灰度值区间的聚类匹配程度;根据每个灰度值区间中每个像素点与该灰度值区间中其余像素点在CT原始图像中的距离的平均值,获取每个灰度值区间中的像素点的聚集度;获取每个灰度值区间的像素点聚集度是为了确定每个灰度值区间中的起始聚类像素点;在获得每个灰度值区间中的起始聚类像素点后可对CT原始图像中每个目标灰度值区间的像素点进行K-means聚类,可实现对CT原始图像中像素点的分类和对CT原始图像的分割;通过本实施例的技术方案实现了对CT原始图像中每个目标灰度值区间的像素点进行K-means聚类,解决了相关技术中,在使用K-means聚类法对肝脏CT图像进行分割时分割效果差,以及计算冗余度大的技术问题。
以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于聚类的肝脏CT影像的图像分割方法,其特征在于,包括以下步骤:
获取腹部CT原始图像;
获取所述CT原始图像的灰度直方图和灰度梯度直方图;
在所述灰度梯度直方图中从最大灰度梯度值开始,向灰度梯度值减小的方向逐级进行像素点数量累加,每累加一次后获取累加像素点数量与CT原始图像中像素点数量的比值,根据当前累加后的比值与上一次累加后的比值获取比值的变化率,根据所述比值的变化率获取所述CT原始图像中的边缘特征像素点;
根据所述灰度直方图中像素点的灰度级分布将所述CT原始图像中的灰度值分为多个灰度值区间;根据每个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值,获取每个灰度值区间的聚类匹配程度;
根据每个灰度值区间中每个像素点与该灰度值区间中其余像素点在所述CT原始图像中的距离的平均值,获取每个灰度值区间中的像素点的聚集度;
根据灰度值区间的聚类匹配程度和聚集度从多个灰度值区间中选取目标灰度值区间;
获取目标灰度值区间的像素点形成的区域中每个像素点与该区域中其余像素点的灰度值的差值平方和;根据差值平方和获取该区域中的聚类起始像素点;依次获取每个目标灰度值区间的像素点形成区域的聚类起始像素点;
根据所述聚类起始像素点对每个目标灰度值区间的像素点形成区域进行聚类,获得多个聚类区域,得到分割后的CT原始图像;
所述灰度值区间的聚类匹配程度通过下式确定:
;式中,/>为第/>个灰度值区间的聚类匹配程度;/>为第/>个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值;/>为第/>个灰度值区间中像素点数量与该灰度值区间中边缘特征像素点数量的比值;/>为灰度值区间个数;表示以自然常数e为底的指数函数;
所述灰度值区间中的像素点的聚集度,公式表示为:
;式中,/>为第/>个灰度值区间的像素点聚集度;/>为CT原始图像中第/>个灰度值区间的像素点/>与第/>个像素点/>之间的距离;/>为第/>个灰度值区间中像素点个数;/>为第/>个灰度值区间中第/>个像素点的灰度值;
所述根据灰度值区间的聚类匹配程度和聚集度从多个灰度值区间中选取目标灰度值区间包括:
根据每个灰度值区间的聚类匹配程度从多个灰度值区间中选取非必要灰度值区间;
将所述非必要灰度值区间相邻的灰度值区间的范围进行调整,获取调整后的灰度值区间的聚类匹配程度和聚集度;当调整后的灰度值区间的聚类匹配程度和聚集度均不再增大时,停止对灰度值区间的调整,并将剩余的全部灰度值区间确定为目标灰度值区间。
2.根据权利要求1所述的基于聚类的肝脏CT影像的图像分割方法,其特征在于,根据所述比值的变化率获取所述CT原始图像中的边缘特征像素点的过程中,还包括设定变化率阈值,当所述比值的变化率大于所述变化率阈值时,停止像素点数量的累加,并将已经累加的像素点确定为所述CT原始图像中的边缘特征像素点。
3.根据权利要求2所述的基于聚类的肝脏CT影像的图像分割方法,其特征在于,所述变化率阈值为0.3,当所述比值的变化率大于0.3时,停止像素点数量的累加,并将已经累加的像素点确定为所述CT原始图像中的边缘特征像素点。
4.根据权利要求1所述的基于聚类的肝脏CT影像的图像分割方法,其特征在于,所述获取每个目标灰度值区间的像素点形成区域的聚类起始像素点的过程为:
获取所述目标灰度值区间的像素点形成的区域中每个像素点与该区域中其余像素点灰度值的差值平方和;获取所述差值平方和中的最大值,并将所述差值平方和中的最大值对应的像素点作为该区域的聚类起始像素点;依次获取每个目标灰度值区间的像素点形成的每个区域的聚类起始像素点。
5.根据权利要求1所述的基于聚类的肝脏CT影像的图像分割方法,其特征在于,所述灰度值区间的聚类匹配程度小于0.3时,该灰度值区间为非必要灰度值区间。
6.根据权利要求1所述的基于聚类的肝脏CT影像的图像分割方法,其特征在于,所述获取腹部CT原始图像后,还包括对所述CT原始图像进行降噪处理。
7.根据权利要求6所述的基于聚类的肝脏CT影像的图像分割方法,其特征在于,所述降噪处理采用高斯滤波降噪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310993928.3A CN116721099B (zh) | 2023-08-09 | 2023-08-09 | 一种基于聚类的肝脏ct影像的图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310993928.3A CN116721099B (zh) | 2023-08-09 | 2023-08-09 | 一种基于聚类的肝脏ct影像的图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116721099A CN116721099A (zh) | 2023-09-08 |
CN116721099B true CN116721099B (zh) | 2023-11-21 |
Family
ID=87870028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310993928.3A Active CN116721099B (zh) | 2023-08-09 | 2023-08-09 | 一种基于聚类的肝脏ct影像的图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116721099B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116993764B (zh) * | 2023-09-26 | 2023-12-08 | 江南大学附属医院 | 一种胃部ct智能分割提取方法 |
CN117057377B (zh) * | 2023-10-11 | 2024-01-12 | 青岛冠成软件有限公司 | 码的识别匹配方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101710425A (zh) * | 2009-12-25 | 2010-05-19 | 南京航空航天大学 | 基于图像灰度梯度和灰度统计直方图的自适应预分割方法 |
CN111539966A (zh) * | 2020-04-16 | 2020-08-14 | 华北电力大学 | 一种基于模糊c均值聚类的比色传感器阵列图像分割方法 |
CN112200800A (zh) * | 2020-10-30 | 2021-01-08 | 福州大学 | 一种基于灰度直方图梯度加权目标方差的电润湿显示器缺陷检测方法 |
WO2022089079A1 (zh) * | 2020-10-30 | 2022-05-05 | 湖北锐世数字医学影像科技有限公司 | 图像处理方法、装置、设备、系统及计算机可读存储介质 |
WO2022141178A1 (zh) * | 2020-12-30 | 2022-07-07 | 深圳市大疆创新科技有限公司 | 图像处理方法及装置 |
CN115311301A (zh) * | 2022-10-12 | 2022-11-08 | 江苏银生新能源科技有限公司 | 一种pcb板焊点缺陷检测方法 |
CN115359052A (zh) * | 2022-10-19 | 2022-11-18 | 南通鼎顺生物科技有限公司 | 基于聚类算法的医疗影像增强方法 |
WO2023134792A2 (zh) * | 2022-12-15 | 2023-07-20 | 苏州迈创信息技术有限公司 | 一种led灯芯缺陷检测方法 |
CN116563296A (zh) * | 2023-07-12 | 2023-08-08 | 吉林省禹语网络科技有限公司 | 一种用于腹部ct图像的识别方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050033139A1 (en) * | 2002-07-09 | 2005-02-10 | Deus Technologies, Llc | Adaptive segmentation of anatomic regions in medical images with fuzzy clustering |
CN109858570A (zh) * | 2019-03-08 | 2019-06-07 | 京东方科技集团股份有限公司 | 图像分类方法及系统、计算机设备及介质 |
-
2023
- 2023-08-09 CN CN202310993928.3A patent/CN116721099B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101710425A (zh) * | 2009-12-25 | 2010-05-19 | 南京航空航天大学 | 基于图像灰度梯度和灰度统计直方图的自适应预分割方法 |
CN111539966A (zh) * | 2020-04-16 | 2020-08-14 | 华北电力大学 | 一种基于模糊c均值聚类的比色传感器阵列图像分割方法 |
CN112200800A (zh) * | 2020-10-30 | 2021-01-08 | 福州大学 | 一种基于灰度直方图梯度加权目标方差的电润湿显示器缺陷检测方法 |
WO2022089079A1 (zh) * | 2020-10-30 | 2022-05-05 | 湖北锐世数字医学影像科技有限公司 | 图像处理方法、装置、设备、系统及计算机可读存储介质 |
WO2022141178A1 (zh) * | 2020-12-30 | 2022-07-07 | 深圳市大疆创新科技有限公司 | 图像处理方法及装置 |
CN115311301A (zh) * | 2022-10-12 | 2022-11-08 | 江苏银生新能源科技有限公司 | 一种pcb板焊点缺陷检测方法 |
CN115359052A (zh) * | 2022-10-19 | 2022-11-18 | 南通鼎顺生物科技有限公司 | 基于聚类算法的医疗影像增强方法 |
WO2023134792A2 (zh) * | 2022-12-15 | 2023-07-20 | 苏州迈创信息技术有限公司 | 一种led灯芯缺陷检测方法 |
CN116563296A (zh) * | 2023-07-12 | 2023-08-08 | 吉林省禹语网络科技有限公司 | 一种用于腹部ct图像的识别方法 |
Non-Patent Citations (2)
Title |
---|
Plant diseased leaf segmentation and recognition by fusion of superpixel, K-means and PHOG;Shanwen Zhang 等;《Optik》;第157卷;第866-872页 * |
基于k-medoids聚类的人脑DTI图像分割算法及其纤维追踪的研究;林涛;《知网》;第2017年卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116721099A (zh) | 2023-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116721099B (zh) | 一种基于聚类的肝脏ct影像的图像分割方法 | |
CN108053417B (zh) | 一种基于混合粗分割特征的3D U-Net网络的肺分割装置 | |
CN110889852B (zh) | 基于残差-注意力深度神经网络的肝脏分割方法 | |
CN101576997B (zh) | 基于二次三维区域生长的腹部器官分割方法 | |
CN109753997B (zh) | 一种ct图像中的肝脏肿瘤自动精确鲁棒分割方法 | |
CN115330800A (zh) | 基于图像处理的放疗靶区自动分割方法 | |
CN108615239B (zh) | 基于阈值技术和灰度投影的舌图像分割方法 | |
CN102831614B (zh) | 基于交互式字典迁移的序列医学图像快速分割方法 | |
CN110717900B (zh) | 基于改进Canny边缘检测算法的受电弓磨耗检测方法 | |
CN110276764A (zh) | 基于K值估计的K-Means水下图像背景分割改进算法 | |
CN113570628A (zh) | 一种基于活动轮廓模型的白细胞分割方法 | |
CN110473224B (zh) | 一种基于kl熵的rsf水平集图像自动分割方法 | |
CN117830308B (zh) | 一种介入手术前后血管造影智能对比分析方法 | |
CN117392155B (zh) | 基于图像处理的高通量基因测序数据处理方法 | |
CN110619633B (zh) | 一种基于多路滤波策略的肝脏图像分割方法 | |
CN113850792A (zh) | 一种基于计算机视觉的细胞分类计数方法及系统 | |
CN116503270A (zh) | 一种基于块匹配三维协同滤波的电子断层图像去噪方法 | |
CN106570837B (zh) | 基于高斯多尺度空间的脑组织mri图像偏移场校正方法 | |
CN112560808B (zh) | 一种基于灰度信息的体内静脉识别方法及装置 | |
CN109300138B (zh) | 一种基于两步纳什均衡改进c-v模型的医学图像分割方法 | |
CN118397023B (zh) | 一种基于mri的三叉神经区域分割方法 | |
CN110647870A (zh) | 一种基于滑动窗计算静息态fMRI数据近似熵的方法 | |
Zhang et al. | Adaptive segmentation of cervical smear image based on GVF Snake model | |
CN117974692B (zh) | 一种基于区域生长的眼科医学影像处理方法 | |
CN118333893B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |