CN106447645B - 增强ct图像中冠脉钙化检测及量化装置和方法 - Google Patents
增强ct图像中冠脉钙化检测及量化装置和方法 Download PDFInfo
- Publication number
- CN106447645B CN106447645B CN201610209092.3A CN201610209092A CN106447645B CN 106447645 B CN106447645 B CN 106447645B CN 201610209092 A CN201610209092 A CN 201610209092A CN 106447645 B CN106447645 B CN 106447645B
- Authority
- CN
- China
- Prior art keywords
- image
- coronary
- region
- enhanced
- algorithm
- 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.)
- Expired - Fee Related
Links
Images
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
- 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/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及图像处理技术,为提供一种在CTA数据中基于模糊超像素聚类的全自动钙化斑块探测、分割及量化装置及方法,实现有效地抑制噪声,本发明采用的技术方案是,增强CT图像中冠脉钙化检测及量化方法,包含三个步骤:首先使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域;然后根据区域中像素点之间的欧式距离和灰度差异使用模糊C均值聚类算法将上述血管区域划分为有限个超像素集;最后使用一种基于灰度直方图的阈值选取办法对超像素集进行筛选,进而得到最终的钙化斑块探测量化结果,并根据分割结果完成对血管的钙化分计算。本发明主要应用于图像处理。
Description
技术领域
本发明涉及图像处理技术,特别涉及增强计算机断层扫描(ComputedTomography,CT)图像中的医学图像数据处理技术,尤其是一种用于在增强CT图像数据中对冠状动脉钙化斑块检测和量化的方法,即增强CT图像中冠脉钙化检测及量化方法。
背景技术
冠心病是指由于冠状动脉血管发生动脉粥样硬化、炎症、栓塞等病变引起血管腔狭窄或阻塞,进而造成心肌缺血、缺氧或坏死的现象,特别是当动脉硬化斑块突发破裂时可能引起心血管突发性堵塞,引起病人猝死,是世界上致死率最高的疾病之一。医学研究表明对冠脉钙化的探测和量化会为预测冠心病的发展提供重要的风险评估依据,从电子计算机断层扫描数据中获得的冠脉钙化积分已经被证明是冠心病的有效的预测参数,而冠脉钙化斑块的形状、大小与位置也是医生采取的介入式治疗手段时需要参考的重要信息。
作为一种非侵入式的而且能够提供三维信息的医学影像检查手段,CT被广泛应用于对冠心病的检测和评估。一般来说,针对冠心病做的标准CT扫描程序包含两个步骤:第一,使用低分辨率的非加强CT对病人进行扫描用以计算钙化分;第二,使用高分辨率的加强CT(Computed Tomography Angioplasty,CTA)扫描,所得数据用以观察冠脉的狭窄程度和评价血管软性斑块的性质。在非加强CT数据中,虽然较低的血管CT值(单位为:Hu)能够使得钙化斑块得到较好的分割,但是计算机却无法区分冠脉内钙化斑块与周围同样具有较高Hu值的软组织,故对于钙化斑块的分辨和标记需要放射医生的人工参与,耗费了大量的人力和时间。同时,为了减少患者受到的辐射剂量,非加强CT扫描时往往采取粗扫描的方式(CT切片的厚度一般为3mm),这在另一方面也使得对钙化斑块体积的量化变得不准确。
随着技术的不断发展,近年来CTA扫描所需辐射剂量不断减少,越来越多的研究团队开始研究在CTA数据中探测、量化冠状动脉钙化斑块并计算钙化积分的方法,这不仅能够减少患者接受扫描次数和受到辐射的剂量,而且由于CTA的分辨率一般高于非加强CT,所以使得对钙化的量化精度也会更高。一些研究人员通过研究CTA与CT中的钙化积分的关系,通过机器学习或其他方法获得两者之间的经验转化系数,进而得到在加强CT中量化钙化的方法。然而CTA数据往往受到CT协议、CT扫描仪生产厂商、造影剂用量等众多因素的影响,上述量化方法往往得不到具有较好的鲁棒性和准确性的结果。
为了克服上述问题,用于CTA图像的钙化探测及量化方法不断被提出。Teβmann等人提出了一种基于阈值的全自动分割和量化钙化病变的方法,该方法包含一种基于冠脉血管区域Hu值直方图的阈值选择方法,使用该方法得到的阈值被证实与人工选择的阈值相对接近。然而,该阈值获得方法依赖于相对光滑的直方图和经验参数,当CTA数据存在较多噪声或直方图波动强烈时,该方法的鲁棒性将不可避免地受到影响;而且由于造影剂的使用,靠近冠脉开口的冠状动脉血管的Hu值可能已经超过了局部钙化斑块,单纯依靠全局阈值对钙化进行探测可能会导致误分割或欠分割问题。
Eilot等提出一种基于模型的全自动的钙化分割及钙化分计算方法。该方法根据灰度分布特征对冠脉血管建立模型并使用该模型对血管钙化进行探测和分割,通过虚拟移除造影剂后的血管灰度分布来计算非增强CT图像中的钙化分。该方法对小型钙化斑块有着较好识别能力,但该方法执行过程相对复杂,基于模型的方法也加大了运算消耗。
Wolterink等提出了一种基于模型识别的全自动钙化识别和量化方法,首先利用神经网络算法提取潜在的钙化像素点,然后对这些候选点使用一系列特征描述,最后使用随机森林对候选点进行分类。该方法实现了在CTA数据中识别和量化钙化,但该方法本身需要大量的训练样本,且随着训练数据的增长,该方法对计算机硬件也有着较高的要求。
综上所述,这些钙化探测及量化方法普遍存在算法探测量化结果不精确、难度大、训练过程繁琐等问题。
发明内容
为克服现有技术的不足,提供一种在CTA数据中基于模糊超像素聚类的全自动钙化斑块探测、分割及量化装置及方法,实现有效地抑制噪声,本发明采用的技术方案是,增强CT图像中冠脉钙化检测及量化方法,包含三个步骤:首先使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域;然后根据区域中像素点之间的欧式距离和灰度差异使用模糊C均值聚类算法将上述血管区域划分为有限个超像素集;最后使用一种基于灰度直方图的阈值选取办法对超像素集进行筛选,进而得到最终的钙化斑块探测量化结果,并根据分割结果完成对血管的钙化分计算。
三个步骤细化为:
步骤1.对冠心病患者静脉注射造影剂之后,采用螺旋增强CT扫描仪对患者的胸腔进行扫描,进而获得患者冠脉的增强CT断层图像;
步骤2.对CT断层数据进行基于三次B样条曲线的层间插值;
步骤3.参考CT断层图像中提供的窗宽、窗位信息进行数据类型转化;
步骤4.检测并分割CT数据中升主动脉;
步骤5.检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割;
步骤6.对冠状动脉区域内的每个像素点构建特征向量;
步骤7.采用C均值模糊聚类算法对特征向量进行聚类,实现对冠状动脉的区域划分;
步骤8.根据灰度直方图计算自适应阈值,依靠阈值筛选步骤7中的冠脉血管区域;
步骤9.参考Agatston积分法对步骤8中获得的钙化区域进行钙化分计算,最终完成冠脉钙化量化。
步骤2中对断层数据进行基于三次B样条曲线的层间插值,以保证CT数据的横向分辨率与纵向分辨率保持统一,其中三次B样条曲线的基函数如下:
其中Gi,3表示三次B样条曲线的各项基函数,i=1,2,3;t为参变量,t∈[0,1];
步骤4中以二值化后的图像为掩模mask,以Hough圆检测的圆心点为标记,对CT断层图像进行形态学重构,从而分割出升主动脉截面,其中形态学重构的定义为:
若g是掩模,f为标记,将从f重构g记为R(f),则
1)f为初始化图像h1;
2)创建结构元素B;
3)重复迭代
hk+1=(hk⊕B)∩g
直至hk+1=hk;其中,⊕表示膨胀操作,标记f必须是g的子集,即f∈g。
步骤5采用动态轮廓—Snake模型拟合升主动脉横截面轮廓以实现对突出位置即冠脉开口的检测,其中,Snake模型下控制点的总能量函数为:
其中,Eint代表曲线弯曲造成的能量,Eimage代表图像力能量,Eext代表外部力能量,v(s)为控制点,v(s)=[x(s),y(s)],s∈[0,1];实际使用基于Balloon模型力改进的Snake模型:
步骤7中,采取模糊C均值聚类算法对对特征向量进行聚类,实现对冠状动脉的区域划分,其中,模糊C均值聚类算法的价值函数为:
其中,Zi为第i个的模糊类的聚类中心,uij∈[0 1],dij=||Zi-Xj||,表示第i个聚类中心与第j个样本之间的欧氏距离;α∈[1,∞),为算法的柔性参数,控制聚类的模糊程度,而算法的聚类中心与隶属度矩阵的迭代公式为:
算而算法聚类数目定义为:
其中,resolution为CT图像的分辨率,skelnum为对冠脉连通域进行骨架化后得到的血管骨架点数目。
增强CT图像中冠脉钙化检测及量化装置,包括:
螺旋增强CT扫描仪,包含扫描床、扫描架、高压发生器、操作台、图像数据处理平台等硬件设备,用于获得患者冠脉的增强CT断层图像;
计算机,用于处理增强CT断层图像;其中,计算上设置有如下模块:
冠状动脉区域生成模块,用于使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域;
聚类模块,用于根据区域中像素点之间的欧式距离和灰度差异使用模糊C均值聚类算法将上述血管区域划分为有限个超像素集;
筛选及输出模块,用于以基于灰度直方图的阈值选取办法对超像素集进行筛选,进而得到最终的钙化斑块探测量化结果,并根据分割结果完成对血管的钙化分计算。
使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域具体是,
对CT断层数据进行基于三次B样条曲线的层间插值;
参考CT断层图像中提供的窗宽、窗位信息进行数据类型转化;
检测并分割CT数据中升主动脉;
检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割。
检测并分割CT数据中升主动脉是以二值化后的图像为掩模mask,以Hough圆检测的圆心点为标记,对CT断层图像进行形态学重构,从而分割出升主动脉截面,其中形态学重构的定义为:
若g是掩模,f为标记,将从f重构g记为R(f),则
1)f为初始化图像h1;
2)创建结构元素B;
3)重复迭代
hk+1=(hk⊕B)∩g
直至hk+1=hk;其中,⊕表示膨胀操作,标记f必须是g的子集,即f∈g。
检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割,具体是采用动态轮廓—Snake模型拟合升主动脉横截面轮廓以实现对突出位置即冠脉开口的检测,其中,Snake模型下控制点的总能量函数为:
其中,Eint代表曲线弯曲造成的能量,Eimage代表图像力能量,Eext代表外部力能量,v(s)为控制点,v(s)=[x(s),y(s)],s∈[0,1];实际使用基于Balloon模型力改进的Snake模型:
聚类模块采取模糊C均值聚类算法对对特征向量进行聚类,实现对冠状动脉的区域划分,其中,模糊C均值聚类算法的价值函数为:
其中,Zi为第i个的模糊类的聚类中心,uij∈[0 1],dij=||Zi-Xj||,表示第i个聚类中心与第j个样本之间的欧氏距离;α∈[1,∞),为算法的柔性参数,控制聚类的模糊程度,而算法的聚类中心与隶属度矩阵的迭代公式为:
算而算法聚类数目定义为:
其中,resolution为CT图像的分辨率,skelnum为对冠脉连通域进行骨架化后得到的血管骨架点数目。
本发明的特点及有益效果是:
本发明提出了一种在CTA数据中基于模糊聚类与血管区域划分的全自动钙化斑块探测、分割及量化方法,该方法对小型钙化斑块有着较好的识别能力,分割结果相对精确,预定义参数少,具有较好的鲁棒性。并且由于采用了聚类的方法预先对冠状动脉区域内的像素点进行了划分,该方法对噪声有着良好的克服能力。研究成果可用于基于增强CT图像的冠状动脉钙化探测与量化工作,有助于减少患者所受到的辐射剂量,为医生提供有效的辅助诊断信息。
附图说明:
图1.CT图像采集系统硬件框架。图中,1扫描床,2扫描架,3高压发生器,4计算机系统,5操作台,6图像数据处理平台。
图2.全自动冠脉钙化探测及量化软件界面。
图3.钙化斑块探测及量化方法整体框架。
图4.增强CT数据中冠状动脉的全自动分割方法框架。
图5.基于模糊聚类算法的冠脉血管区域划分结果剖面示意图。
图6.基于高斯函数拟合的自适应分割阈值确定。
图7.钙化检测与量化结果三维示意图。
具体实施方式
本发明提出了一种CTA数据中基于模糊超像素聚类的全自动钙化斑块探测、分割及量化方法,首先使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域;然后根据区域中像素点之间的欧式距离,灰度差异使用模糊C均值聚类算法将上述血管区域划分为有限个超像素集;最后使用一种基于灰度直方图的阈值选取办法对超像素集进行筛选进而得到最终的钙化斑块探测量化结果,并根据分割结果完成对血管钙化分计算。
一种自动的CTA数据中钙化斑块探测及量化方法,该方法包括:
步骤1.对冠心病患者静脉注射造影剂之后,采用64排螺旋增强CT扫描仪对患者的胸腔进行扫描,进而获得患者冠脉的增强CT断层图像;
步骤2.对CT断层数据进行基于三次B样条曲线的层间插值;
步骤3.参考CT断层图像中提供的窗宽、窗位信息进行数据类型转化;
步骤4.检测并分割CT数据中升主动脉;
步骤5.检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割;
步骤6.对冠状动脉区域内的每个像素点构建特征向量;
步骤7.采用C均值模糊聚类算法对特征向量进行聚类,实现对冠状动脉的区域划分;
步骤8.根据灰度直方图计算自适应阈值,依靠阈值筛选步骤7中的冠脉血管区域;
步骤9.参考Agatston积分法对步骤8中获得的钙化区域进行钙化分计算,最终完成冠脉钙化量化。
步骤2中对断层数据进行基于三次B样条曲线的层间插值,以保证CT数据的横向分辨率与纵向分辨率保持统一,其中三次B样条曲线的基函数如下:
步骤4中以二值化后的图像为掩模(mask),以Hough圆检测的圆心点为标记,对CT断层图像进行形态学重构,从而分割出升主动脉截面,其中形态学重构的定义为:
若g是掩模,f为标记,将从f重构g记为R(f),则
1)f为初始化图像h1;
2)创建结构元素B;
3)重复迭代
hk+1=(hk⊕B)∩g
直至hk+1=hk。
其中,⊕表示膨胀操作,标记f必须是g的子集,即f∈g。
步骤5采用动态轮廓—Snake模型拟合升主动脉横截面轮廓以实现对突出位置即冠脉开口的检测。其中,Snake模型下控制点的总能量函数为:
实际使用基于Balloon模型力改进的Snake模型:
步骤7中,采取模糊C均值聚类算法对对特征向量进行聚类,实现对冠状动脉的区域划分,其中,模糊C均值聚类算法的价值函数为:
其中,Zi为第i个的模糊类的聚类中心,uij∈[0 1],dij=||Zi-Xj||,表示第i个聚类中心与第j个样本之间的欧氏距离;α∈[1,∞),为算法的柔性参数,控制聚类的模糊程度而算法的聚类中心与隶属度矩阵的迭代公式为:
算而算法聚类数目定义为:
其中,resolution为CT图像的分辨率,skelnum为对冠脉连通域进行骨架化后得到的血管骨架点数目。
本发明的具体实施方式如图1所示:
101.对冠心病患者静脉注射造影剂之后,采用64排螺旋增强CT扫描仪对患者的胸腔进行扫描,进而获得患者冠脉的增强CT断层图像;
102.使用B样条曲线拟合不同空间位置的断层图像对应像素点的CT值变化曲线,然后对该曲线进行插值,使用插值后的得到的新数据组成新的CT数据集,以达到增加断层数目,保持CT数据横向分辨率与纵向分辨率统一的目的。其中三次B样条曲线的基函数如下:
其中Gi,3表示三次B样条曲线的各项基函数,t为参变量。
103.参考CT断层图像中提供的窗宽、窗位信息信息,将DICOM格式图像中的int16型数据转化为八位二进制的灰度数据,以提高计算效率。
104.检测并分割CT数据中升主动脉;
采取Hough圆检测的方式实现对CT断层图像中升主动脉位置的检测,以Hough圆内像素平均灰度值的3/4作为分割阈值对图像进行二值化,以二值化后的图像为掩模(mask),以Hough圆检测的圆心点为标记,对CT断层图像进行形态学重构,从而分割出升主动脉截面,将每张CT图像中升主动脉截面叠加,可以得到升主动脉的三维分割结果,其中形态学重构的定义为:
若g是掩模,f为标记,将从f重构g记为R(f),则
1)f为初始化图像h1;
2)创建结构元素B;
3)重复迭代
hk+1=(hk⊕B)∩g
直至hk+1=hk。
其中,⊕表示膨胀操作,标记f必须是g的子集,即f∈g。
105.检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割;
采用动态轮廓—Snake模型拟合升主动脉横截面轮廓以实现对突出位置即冠脉开口的检测,同时也为后续的基于区域生长的冠状动脉分割工作提供了种子点位置;最后,本发明使用一种经验值确定的低阈值的区域生长方法获得冠状动脉区域。Snake模型由一组首尾相接的控制点v(s)=[x(s),y(s)],s∈[0,1]组成,定义Snake模型下控制点的总能量函数为:
其中,Eint代表曲线弯曲造成的能量,Eimage代表图像力能量,Eext代表外部力能量。在研究过程中本发明使用了基于Balloon模型力改进的Snake模型,增加了外部力:
其中,n(s)是Snake轮廓曲线在点v(s)处的单位法向量,k1是外力的大小,表示轮廓s处的梯度。该模型保证了活动的轮廓的膨胀特性。
106.对冠状动脉区域内的每个像素点构建特征向量;
对已经提取分割出的冠脉血管区域内的每个像素点构造特征向量Xi,定义:Xi=[r,c,h,v],其中r,c,h为该像素点的空间坐标,v为像素点CT值进行映射后的灰度值。
107.采用C均值模糊聚类算法对特征向量进行聚类,实现对冠状动脉的区域划分;
采取模糊C均值聚类算法对对特征向量进行聚类,实现对冠状动脉的区域划分,其中,模糊C均值聚类算法的价值函数为:
其中,Zi为第i个的模糊类的聚类中心,uij∈[0 1],dij=||Zi-Xj||,表示第i个聚类中心与第j个样本之间的欧氏距离;α∈[1,∞),为算法的柔性参数,控制聚类的模糊程度。而算法的聚类中心与隶属度矩阵的迭代公式为:
算法具体执行步骤为:
Step1:用取值在[0,1]区间的随机数初始化满足式1的隶属度矩阵U;
Step2:根据式(5)计算m个类的聚类中心Zi,i=1,…,m。
Step3:根据式(4)计算价值函数,如果它下雨某个确定的阈值ΔJ,则算法停止,否则进行步骤4;
Step4:使用式(6)计算新的隶属度矩阵,返回步骤2。
而算法聚类数目定义为:
其中,resolution为CT图像的分辨率,skelnum为对冠脉连通域进行骨架化后得到的血管骨架点数目。
108.根据灰度直方图计算自适应阈值,依靠阈值筛选步骤7中的冠脉血管区域;
使用一维高斯函数对直方图进行拟合,进而选取T=μ+1.5*σ为经验阈值,选取平均CT值大于阈值的血管区域为钙化区域。
109.参考Agatston积分法对步骤8中获得的钙化区域进行钙化分计算;
参考Agatston积分法对步骤8中获得的钙化区域进行钙化分计算,最终完成冠脉钙化量化。
Claims (8)
1.一种增强CT图像中冠脉钙化检测及量化方法,其特征是,包含三个步骤:首先使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域;然后根据区域中像素点之间的欧式距离和灰度差异使用模糊C均值聚类算法将血管区域划分为有限个超像素集;最后使用一种基于灰度直方图的阈值选取办法对超像素集进行筛选,进而得到最终的钙化斑块探测量化结果,并根据分割结果完成对血管的钙化分计算;
三个步骤细化为:
步骤1.对冠心病患者静脉注射造影剂之后,采用螺旋增强CT扫描仪对患者的胸腔进行扫描,进而获得患者冠脉的增强CT断层图像;
步骤2.对CT断层数据进行基于三次B样条曲线的层间插值;
步骤3.参考CT断层图像中提供的窗宽、窗位信息进行数据类型转化;
步骤4.检测并分割CT数据中升主动脉;
步骤5.检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割;
步骤6.对冠状动脉区域内的每个像素点构建特征向量;
步骤7.采用C均值模糊聚类算法对特征向量进行聚类,实现对冠状动脉的区域划分;
步骤8.根据灰度直方图计算自适应阈值,依靠阈值筛选步骤7中的冠状动脉区域,筛选出钙化区域;
步骤9.参考Agatston积分法对步骤8中获得的钙化区域进行钙化分计算,最终完成冠脉钙化量化。
2.如权利要求1所述的增强CT图像中冠脉钙化检测及量化方法,其特征是,步骤2中对断层数据进行基于三次B样条曲线的层间插值,以保证CT数据的横向分辨率与纵向分辨率保持统一,其中三次B样条曲线的基函数如下:
其中Gi,3表示三次B样条曲线的各项基函数,i=1,2,3;t为参变量,t∈[0,1];
步骤4中以二值化后的图像为掩模mask,以Hough圆检测的圆心点为标记,对CT断层图像进行形态学重构,从而分割出升主动脉截面,其中形态学重构的定义为:
若g是掩模,f为标记,将从f重构g记为R(f),则
1)h1为初始化图像;
2)创建结构元素B;
3)重复迭代
直至hk+1=hk;其中,⊕表示膨胀操作,标记f必须是g的子集,即f∈g。
3.如权利要求1所述的增强CT图像中冠脉钙化检测及量化方法,其特征是,步骤5采用动态轮廓—Snake模型拟合升主动脉横截面轮廓以实现对突出位置即冠脉开口的检测,其中,Snake模型下控制点的总能量函数为:
Etotal=∫0 1[Eint(v(s))+Eimage(v(s))+Eext(v(s))]ds
其中,Eint代表曲线弯曲造成的能量,Eimage代表图像力能量,Eext代表外部力能量,v(s)为控制点,v(s)=[x(s),y(s)],s∈[0,1];实际使用基于Balloon模型力改进的Snake模型:
5.一种增强CT图像中冠脉钙化检测及量化装置,其特征是,包括:
螺旋增强CT扫描仪,包含扫描床、扫描架、高压发生器、操作台、图像数据处理平台硬件设备,用于获得患者冠脉的增强CT断层图像;
计算机,用于处理增强CT断层图像;其中,计算机上设置有如下模块:
冠状动脉区域生成模块,用于使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域;
聚类模块,用于根据区域中像素点之间的欧式距离和灰度差异使用模糊C均值聚类算法将血管区域划分为有限个超像素集;
筛选及输出模块,用于以基于灰度直方图的阈值选取办法对超像素集进行筛选,进而得到最终的钙化斑块探测量化结果,并根据分割结果完成对血管的钙化分计算;使用种子点自动选择的低阈值区域生长获得包含钙化斑块的冠状动脉区域具体是,
对CT断层数据进行基于三次B样条曲线的层间插值;
参考CT断层图像中提供的窗宽、窗位信息进行数据类型转化;
检测并分割CT数据中升主动脉;
检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割。
7.如权利要求5所述的增强CT图像中冠脉钙化检测及量化装置,其特征是,检测升主动脉区域内冠状动脉开口位置,并使用区域生长方式实现对冠状动脉区域的分割,具体是:采用动态轮廓—Snake模型拟合升主动脉横截面轮廓以实现对突出位置即冠脉开口的检测,其中,Snake模型下控制点的总能量函数为:
Etotal=∫0 1[Eint(v(s))+Eimage(v(s))+Eext(v(s))]ds
其中,Eint代表曲线弯曲造成的能量,Eimage代表图像力能量,Eext代表外部力能量,v(s)为控制点,v(s)=[x(s),y(s)],s∈[0,1];实际使用基于Balloon模型力改进的Snake模型:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610209092.3A CN106447645B (zh) | 2016-04-05 | 2016-04-05 | 增强ct图像中冠脉钙化检测及量化装置和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610209092.3A CN106447645B (zh) | 2016-04-05 | 2016-04-05 | 增强ct图像中冠脉钙化检测及量化装置和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106447645A CN106447645A (zh) | 2017-02-22 |
CN106447645B true CN106447645B (zh) | 2019-10-15 |
Family
ID=58183869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610209092.3A Expired - Fee Related CN106447645B (zh) | 2016-04-05 | 2016-04-05 | 增强ct图像中冠脉钙化检测及量化装置和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106447645B (zh) |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107330888A (zh) * | 2017-07-11 | 2017-11-07 | 中国人民解放军第三军医大学 | 基于cta图像的动态心脏各腔室分割方法 |
CN107451995B (zh) * | 2017-07-26 | 2020-04-21 | 深圳先进技术研究院 | 在cta图像中提取心血管方法、装置、设备及存储介质 |
CN107527341B (zh) * | 2017-08-30 | 2020-05-19 | 上海联影医疗科技有限公司 | 血管造影图像的处理方法和系统 |
CN107590803B (zh) * | 2017-09-13 | 2020-11-06 | 东软医疗系统股份有限公司 | 一种在双能量ct图像中显示roi的装置 |
CN107610117A (zh) * | 2017-09-18 | 2018-01-19 | 广州慧扬健康科技有限公司 | 冠状动脉管腔图像自动分割系统 |
CN108198239B (zh) * | 2017-12-27 | 2021-07-23 | 中山大学 | 一种实现血管动态仿真的三维可视化方法 |
CN108416769B (zh) * | 2018-03-02 | 2021-06-04 | 成都斯斐德科技有限公司 | 基于预处理的ivoct图像易损斑块自动检测方法 |
CN108764531A (zh) * | 2018-04-28 | 2018-11-06 | 南京秀强信息科技有限公司 | 基于tls提取单木胸径的聚类圆拟合优化算法 |
CN108805923B (zh) * | 2018-06-01 | 2022-04-01 | 东软医疗系统股份有限公司 | 生理参数获取方法和装置 |
CN109389592B (zh) * | 2018-09-30 | 2021-01-26 | 数坤(北京)网络科技有限公司 | 计算冠脉钙化积分的方法、装置及系统 |
CN109544585A (zh) * | 2018-12-19 | 2019-03-29 | 中国石油大学(华东) | 一种基于轻量化卷积神经网络的胆石病ct医疗图像数据增强方法 |
CN111383231B (zh) * | 2018-12-28 | 2023-10-27 | 成都皓图智能科技有限责任公司 | 一种基于3d图像的图像分割方法、装置及系统 |
CN109979591B (zh) * | 2019-03-12 | 2021-01-01 | 众安信息技术服务有限公司 | 一种基于图神经网络分析斑块进展因子的方法及装置 |
CN110033466B (zh) * | 2019-04-01 | 2020-12-18 | 数坤(北京)网络科技有限公司 | 一种基于多灰度级的冠脉拉直图像分割边界确定方法 |
CN109846465B (zh) * | 2019-04-01 | 2021-09-28 | 数坤(北京)网络科技股份有限公司 | 一种基于亮度分析的血管钙化误报检测方法 |
CN110021016B (zh) * | 2019-04-01 | 2020-12-18 | 数坤(北京)网络科技有限公司 | 一种钙化检测方法 |
CN110189296B (zh) * | 2019-04-16 | 2022-05-10 | 上海鹰瞳医疗科技有限公司 | 眼底图像血管壁反光状态标记方法及设备 |
CN110264461B (zh) * | 2019-06-25 | 2020-10-27 | 南京工程学院 | 基于超声乳腺肿瘤图像的微小钙化点自动检测方法 |
CN113240684B (zh) * | 2019-07-05 | 2024-05-17 | 心医国际数字医疗系统(大连)有限公司 | 对提取出含冠脉中心线的掩膜数据转换为中心线点集方法 |
CN114126491B (zh) * | 2019-07-22 | 2023-04-04 | 西门子医疗有限公司 | 血管造影图像中的冠状动脉钙化的评估 |
CN110600126A (zh) * | 2019-09-19 | 2019-12-20 | 江苏大学附属医院 | 一种糖尿病足下肢动脉钙化斑块图像的辅助评价方法 |
CN110780080B (zh) * | 2019-11-08 | 2024-07-12 | 安邦(厦门)生物科技有限公司 | 一种血型分析仪及血型分析方法 |
CN111047607B (zh) * | 2019-12-09 | 2023-07-11 | 赛诺威盛科技(北京)股份有限公司 | 一种自动分割冠状动脉的方法 |
CN111353996B (zh) | 2020-04-08 | 2024-03-01 | 东软医疗系统股份有限公司 | 一种血管钙化检测方法及装置 |
CN111612743B (zh) * | 2020-04-24 | 2023-05-02 | 杭州电子科技大学 | 一种基于ct图像的冠状动脉中心线提取方法 |
CN111667467B (zh) * | 2020-05-28 | 2021-01-26 | 江苏大学附属医院 | 基于聚类算法的下肢血管钙化指数多参数累积计算方法 |
CN113744171B (zh) * | 2020-05-28 | 2023-11-14 | 上海微创卜算子医疗科技有限公司 | 血管钙化图像分割方法、系统和可读存储介质 |
CN111568459B (zh) * | 2020-06-03 | 2021-07-13 | 江苏省中医院 | 一种ct肺动脉成像扫描延迟时间的确定方法 |
CN112288752B (zh) * | 2020-10-29 | 2021-08-27 | 中国医学科学院北京协和医院 | 一种基于胸部平扫ct的冠脉钙化灶全自动分割方法 |
CN112686909B (zh) * | 2020-12-28 | 2024-04-19 | 淮阴工学院 | 基于异质多区域ct扫描数据处理的多相隐式曲面重建方法 |
CN113066064B (zh) * | 2021-03-29 | 2023-06-06 | 郑州铁路职业技术学院 | 基于人工智能的锥束ct图像生物结构识别及三维重建系统 |
CN113077432B (zh) * | 2021-03-30 | 2024-01-05 | 中国人民解放军空军军医大学 | 基于冠状动脉cta图像粥样硬化斑块综合特征的患者危险度分级系统 |
CN113096097A (zh) * | 2021-04-13 | 2021-07-09 | 上海商汤智能科技有限公司 | 血管图像检测方法及检测模型训练方法、相关装置、设备 |
CN113034491B (zh) * | 2021-04-16 | 2021-10-08 | 北京安德医智科技有限公司 | 一种冠脉钙化斑块检测方法及装置 |
CN114049282B (zh) * | 2022-01-07 | 2022-05-24 | 浙江大学 | 一种冠状动脉的构建方法、装置、终端及存储介质 |
CN114943699B (zh) * | 2022-05-16 | 2023-01-17 | 北京医准智能科技有限公司 | 分割模型训练方法、冠脉钙化斑块的分割方法、相关装置 |
CN116342588B (zh) * | 2023-05-22 | 2023-08-11 | 徕兄健康科技(威海)有限责任公司 | 一种脑血管图像增强方法 |
CN117765107B (zh) * | 2023-11-20 | 2024-07-12 | 河北港口集团有限公司秦皇岛中西医结合医院 | 冠状动脉疾病ct扫描的成像质量优化方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003046833A2 (en) * | 2001-11-24 | 2003-06-05 | Image Analysis, Inc. | Automatic detection and quantification of coronary and aortic calcium |
CN104091346A (zh) * | 2014-07-24 | 2014-10-08 | 东南大学 | 一种全自动ct图像冠状动脉钙化分数计算方法 |
CN104637044A (zh) * | 2013-11-07 | 2015-05-20 | 中国科学院深圳先进技术研究院 | 钙化斑块及其声影的超声图像提取系统 |
-
2016
- 2016-04-05 CN CN201610209092.3A patent/CN106447645B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003046833A2 (en) * | 2001-11-24 | 2003-06-05 | Image Analysis, Inc. | Automatic detection and quantification of coronary and aortic calcium |
CN104637044A (zh) * | 2013-11-07 | 2015-05-20 | 中国科学院深圳先进技术研究院 | 钙化斑块及其声影的超声图像提取系统 |
CN104091346A (zh) * | 2014-07-24 | 2014-10-08 | 东南大学 | 一种全自动ct图像冠状动脉钙化分数计算方法 |
Non-Patent Citations (3)
Title |
---|
Matthias Teßmann 等.Automatic detection and quantification of coronary calcium on 3D CT angiography data.《Computer Science-Research and Development》.2011,正文第3节. * |
基于 FCM 聚类算法的颅内出血 CT 图像分割;王海波 等;《CT理论与应用研究》;20090630;第18卷(第2期);正文第2、4节 * |
面向CTA图像的冠脉血管分割算法研究和血管狭窄度分析;杨栋;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20160115;正文第29-30页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106447645A (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106447645B (zh) | 增强ct图像中冠脉钙化检测及量化装置和方法 | |
CN107230206B (zh) | 一种基于多模态数据的超体素序列肺部图像的3d肺结节分割方法 | |
CN106600609B (zh) | 一种医学图像中脊柱的分割方法及系统 | |
CN103337096B (zh) | 一种冠状动脉ct造影图像钙化点检测方法 | |
CN101065776B (zh) | 分割解剖学结构的三维数字表示的方法和工作站 | |
Golemati et al. | Using the Hough transform to segment ultrasound images of longitudinal and transverse sections of the carotid artery | |
CN105279759B (zh) | 结合上下信息窄带约束的腹腔主动脉瘤外轮廓分割方法 | |
CN104794708B (zh) | 一种基于多特征学习的动脉粥样硬化斑块成分分割方法 | |
WO2018023917A1 (zh) | 一种下肢血管的提取方法及系统 | |
CN107545584A (zh) | 医学图像中定位感兴趣区域的方法、装置及其系统 | |
CN107590809A (zh) | 肺分割方法及医学成像系统 | |
CN111227822A (zh) | 基于心肌血流量和ct图像的冠状动脉血流储备分数计算方法 | |
US20100278405A1 (en) | Scoring Method for Imaging-Based Detection of Vulnerable Patients | |
Jeyavathana et al. | A survey: analysis on pre-processing and segmentation techniques for medical images | |
CN109215064A (zh) | 一种基于超像素向导的医学图像配准方法 | |
CN110969618B (zh) | 一种基于动态超声造影的肝肿瘤血管生成定量分析方法 | |
CN107909589A (zh) | 一种结合C‑V水平集和GrabCut算法的牙齿图像分割方法 | |
CN105139377A (zh) | 一种腹部ct序列图像肝脏的快速鲁棒自动分割方法 | |
CN109753997A (zh) | 一种ct图像中的肝脏肿瘤自动精确鲁棒分割方法 | |
CN108268897A (zh) | 一种乳腺肿瘤多模态超声多层次计算机辅助诊断方法 | |
Koundal et al. | Challenges and future directions in neutrosophic set-based medical image analysis | |
Liu et al. | Automatic segmentation algorithm of ultrasound heart image based on convolutional neural network and image saliency | |
CN112562058B (zh) | 一种基于迁移学习的颅内血管模拟三维模型快速建立方法 | |
CN110458859A (zh) | 一种基于多序列mri的多发性骨髓瘤病灶的分割系统 | |
Cheng et al. | Airway segmentation and measurement in CT images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191015 |