一种基于血管内超声图像的血管内腔分割方法
【专利摘要】本发明公开了一种基于血管内超声图像的血管内腔分割方法，包括一个确定血管内腔（即血管内膜所包含区域）种子点的过程；包括一个利用模糊连通(FuzzyConnectedness)算法计算图像中从每个象素点与血管内腔种子点模糊连通强度，获得模糊连通强度图像的过程；包括一个利用超声图像梯度信息确定模糊连通强度阈值，并根据模糊连通强度阈值和模糊连通强度图像确定血管内腔边界的过程。本发明通过自动确定种子点和模糊连通阈值，保证了分割过程的自动性。基于模糊连通算法的处理方法，不仅保证了分割方法的简单、有效性，而且避免了现有算法模型的复杂性和对成像条件的依赖性。
【专利说明】-种基于血管内超声图像的血管内腔分割方法
[0001]
【技术领域】
[0002] 本发明涉及医学图像处理领域，特别是涉及一种模糊连通（Fuzzy Connectedness)算法、应用于血管内超声（IVUS :Intravascular ultrasound)图像的血管 内腔分割方法。
[0003]
【背景技术】
[0004] 血管内超声（IVUS: Intravascular Ultrasound)图像，不仅能显不血管内腔形 态，还能显示血管壁分层结构，对动脉粥样硬化等心血管疾病的诊断和治疗具有非常重要 的价值。基于IVUS图像诊断动脉粥样硬化需要粥样硬化图像的特征如血管内腔面积、斑块 面积等量化指标，这些量化指标的准确提取依赖于有效的图像分割。人工分割即由医生手 动勾画血管内腔、中外膜边界等，不仅费时费力，而且受医生经验等主观性的限制。因此， 用计算机算法准确、快速、自动地分割IVUS图像就显得很有必要。目前，基于IVUS图像的 血管内腔计算机自动分割算法主要有三类：第一类为统计学方法（G. Mendizabal-Ruiz, M. Rivera, et al. , "A probabilistic segmentation method for the identification of luminal borders in intravascular ultrasound images"， IEEE Conference on Computer Vision and Pattern Recognition, pp. 1-8，2008.)，对图像的灰度分布进行统 计学建模实现IVUS图像分割，但IVUS图像中的伪影、钙化等复杂的图像特征将大大降低 统计建模的准确性；第二类是机器学习的方法（l.E. G. Bovenkamp, J. Dijkstra, J. G. Bosch, et al. , ^Multi-agent segmentation of IVUS images"，Patten Recognition, Vol. 37, No. 4, pp. 647-663, 2004; 2. G. Unal, S. Bucher, S. Carlier, et al., "Shape-driven segmentation of the arterial wall in intravascular ultrasound images"，IEEE Trans. On information technology in biomedicine, Vol. 12, No. 3, pp. 335-346，2008.)，该类方法模型复杂，实际应用时受到诸多限制；第三类是基于活动 轮廓线模型的方法（1.张麒，汪源源等，"活动轮廓模型和Contourlet多分辨率分析分割 血管内超声图像"，光学精密工程，Vol. 16，No. 11，pp. 2301-311，2008; 2. X. Zhu, P. Zhang, J. Shao, et al. , "A snake-based method for segmentation of intravascular ultrasound images and its in vivo validation"，Ultrasonics, Vol. 51, pp. 181-189, 2011.)，该类方法简单易行，但是往往需要给定初始轮廓线，而且，分割结果易受噪声等复 杂图像特征的影响。上述几类基于IVUS图像的血管内腔分割方法都使用了一定的预先 设定的模型，如图像灰度分布模型、形状模型等，而这些模型往往依赖于具体的成像条件。 为更好地实现对IVUS图像中血管内腔的有效分割，本发明提出一种基于模糊连通（Fuzzy Connectedness)算法的IVUS图像血管内腔自动分割算法。
[0005]
【发明内容】
[0006] 本发明的目的在于针对现有技术的不足，提供一种更简单、更有效、更通用的IVUS 图像血管内腔自动分割方法。
[0007] 本发明的技术解决方案如下： 一种基于血管内超声图像的血管内腔分割方法，包括一个确定血管内腔（即血管内膜 所包含区域）种子点的过程；包括一个利用模糊连通（Fuzzy Connectedness)算法计算 图像中从每个象素点与血管内腔种子点模糊连通强度，获得相应的模糊连通强度图像的过 程；包括一个利用超声图像梯度信息确定模糊连通强度阈值，并根据模糊连通强度阈值和 模糊连通强度图像确定血管内腔边界的过程。
[0008] 进一步的，在一个利用血管内超声图像确定血管内腔（即血管内膜所包含区域）种 子点的过程中，首先根据距离图像中心点同一距离的像素平均灰度值确定血管内腔种子 点，即在大于导管半径的距离中，以最小像素平均灰度值对应的距离为半径，图像中心为圆 心的圆周所在的像素点。
[0009] 进一步的，在一个利用模糊连通算法计算图像中从每个像素点与血管内腔种子点 模糊连通强度，获得相应的模糊连通强度图像，并根据模糊连通强度图像和超声图像梯度 信息确定模糊连通强度阈值的过程中，首先基于模糊连通算法计算图像中从每个象素点与 血管内腔种子点模糊连通强度，获得相应的模糊连通强度图像；然后计算待处理血管内超 声图像的梯度图像，最后根据血管内超声图像的梯度信息和模糊连通强度图像确定模糊 连通强度阈值。
[0010] 进一步的，在一个根据血管内超声图像梯度信息和模糊连通强度图像确定模糊 连通强度阈值的过程中，首先利用整幅梯度图像的平均值&_和标准差G td确定一梯度阈 值^s(如's4_n+0. 5*Gstd)，并将梯度高于梯度阈值^s的像素点看作高梯度象素点，最 后连续变化的模糊连通强度阈值(如从0.95逐步降低)对与模糊连通强度图像作阈值处理， 并考察高于连通强度阈值且与种子点相连的连通区域，若其外侧边界对应的高梯度象素点 的比例高于某一给定的设定的值(如0. 65)，则以此连通强度阈值作为最终的模糊连通强度 阈值。
[0011] 进一步的，在一个根据模糊连通强度阈值和模糊连通强度图像确定血管内腔边界 的过程中，首先以模糊连通强度阈值对模糊连通强度图像作阈值处理，然后将模糊连通强 度图像中模糊连通强度高于模糊连通强度阈值且与种子点相连的连通区域作为血管内腔 区域。
[0012] 与现有技术相比，本发明的有益效果如下： 本发明和已有技术相比，其效果是积极的和明显的。通过自动确定种子点和模糊连通 阈值，本发明保证了分割过程的自动性。基于模糊连通算法的处理方法，不仅保证了分割方 法的简单、有效性，而且避免了现有算法模型的复杂性和对成像条件的依赖性。
【专利附图】
【附图说明】
[0013] 图1是本发明基于血管内超声图像的血管内腔分割方法流程图； 图2是本发明血管内超声图像示意图； 图3是本发明极坐标变换后的血管内超声图像示意图； 图4是本发明不同半径像素平均灰度示意图； 图5是本发明血管内腔种子点示意图； 图6是本发明模糊连通强度示意图； 图7是本发明图像梯度示意图； 图8是本发明高梯度像素示意图； 图9是本发明血管内腔分割结果示意图(极坐标）； 图10是本发明血管内腔分割结果示意图(原始直角坐标)。
[0014]
【具体实施方式】
[0015] 下面结合实施案例和附图对本发明作进一步说明，但不应以此限制本发明的保护 范围。
[0016] 图1是本发明一种基于血管内超声图像的血管内腔分割方法的流程图。如图所 示，一种基于血管内超声（IVUS :Intravascular Ultrasound)图像的血管内腔分割方法， 包括一个确定血管内腔（即血管内膜所包含区域）种子点的过程；包括一个利用模糊连通 (Fuzzy Connectedness)算法计算图像中从每个象素点与血管内腔种子点模糊连通强度， 获得模糊连通强度图像的过程；包括一个利用超声图像梯度信息确定模糊连通强度阈值， 并根据模糊连通强度阈值和模糊连通强度图像确定血管内腔边界的过程。
[0017] 对于IVUS图像，在实际处理过程中，为方便计算，首先将原始图像（图2)转化成极 坐标图像（图3)，然后进行各类计算和处理。在极坐标图像中，原点为原始IVUS图像的中 心。纵坐标为采样半径，横坐标为扫描角度。
[0018] 进一步的，在一个利用IVUS图像确定血管内腔（即血管内膜所包含区域）种子点 的过程中，首先根据距离图像中心点同一距离的像素平均灰度值（图4)确定血管内腔种子 点，即在大于导管半径的距离中，以最小像素平均灰度值对应的距离为半径，图像中心为圆 心的圆周所在的像素为血管内腔种子点。因为半径相同，所以在极坐标图像中，内腔种子点 为一条直线（图5中白线所示)。
[0019] 进一步的，在一个利用模糊连通算法计算图像中从每个像素点与血管内腔种子点 模糊连通强度，获得相应的模糊连通强度图像，并根据模糊连通强度图像和超声图像梯度 信息确定模糊连通强度阈值的过程中，首先基于模糊连通算法计算图像中从每个象素点与 血管内腔种子点模糊连通强度，获得相应的模糊连通强度图像（图6);然后用式（1)作算子 计算极坐标图像梯度（图7)(考虑到IVUS图像自身的特征，梯度图像中，小于零的梯度置 零）；最后根据血管内超声图像的梯度信息和模糊连通强度图像确定模糊连通强度阈值。
[0020]
【权利要求】
1. 一种基于血管内超声图像的血管内腔分割方法，其特征在于，该方法包括一个确定 血管内腔(即血管内膜所包含区域）种子点的过程；包括一个利用模糊连通算法计算图像中 从每个象素点与血管内腔种子点模糊连通强度，获得相应的模糊连通强度图像的过程；包 括一个利用超声图像梯度信息确定模糊连通强度阈值，并根据模糊连通强度阈值和模糊连 通强度图像确定血管内腔边界的过程。
2. 根据权利要求1所述的一种基于血管内超声图像的血管内腔分割方法，其特征在 于：在一个利用血管内超声图像确定血管内腔种子点的过程中，首先根据距离图像中心点 同一距离的像素平均灰度值确定血管内腔种子点，即在大于导管半径的距离中，以最小像 素平均灰度值对应的距离为半径，图像中心为圆心的圆周所在的像素点。
3. 根据权利要求1所述的一种基于血管内超声图像的血管内腔分割方法，其特征在 于：在一个利用模糊连通算法计算图像中从每个像素点与血管内腔种子点模糊连通强度， 获得相应的模糊连通强度图像，并根据模糊连通强度图像和超声图像梯度信息确定模糊连 通强度阈值的过程中，首先基于模糊连通算法计算图像中从每个象素点与血管内腔种子点 模糊连通强度，获得相应的模糊连通强度图像；然后计算待处理血管内超声图像的梯度图 像，最后根据血管内超声图像的梯度信息和模糊连通强度图像确定模糊连通强度阈值。
4. 根据权利要求3所述的一种基于血管内超声图像的血管内腔分割方法，其特征在 于：在一个根据血管内超声图像梯度信息和模糊连通强度图像确定模糊连通强度阈值的过 程中，首先根据整幅图像梯度的平均值和标准差Gtd确定一梯度阈值Gthes，并将梯度高 于梯度阈值Gthes的像素点看作高梯度象素点，最后连续变化的模糊连通强度阈值对与模糊 连通强度图像作阈值处理，并考察高于连通强度阈值且与种子点相连的连通区域，若其外 侧边界对应的高梯度象素点的比例高于给定的设定的值，则以此连通强度阈值作为最终的 模糊连通强度阈值。
5. 根据权利要求4所述的一种基于血管内超声图像的血管内腔分割方法，其特征在 于：首先利用整幅梯度图像的平均值&_和标准差Gtd确定一梯度阈值，确定梯度阈值的公 式为 Gthes^mean+O. 5*6^。
6. 根据权利要求4所述的一种基于血管内超声图像的血管内腔分割方法，其特征在 于：将梯度高于梯度阈值的像素点看作高梯度象素点，最后连续变化的模糊连通强度 阈值对与模糊连通强度图像作阈值处理中，模糊连通强度阈值从〇. 95逐步降低。
7. 根据权利要求4所述的一种基于血管内超声图像的血管内腔分割方法，其特征 在于：若其外侧边界对应的高梯度象素点的比例高于某一给定的设定的值，设定的值为 0. 65〇
8. 根据权利要求1所述基于血管内超声图像的血管内腔分割方法，其特征在于：在一 个根据模糊连通强度阈值和模糊连通强度图像确定血管内腔边界的过程中，首先以模糊连 通强度阈值对模糊连通强度图像作阈值处理，然后将模糊连通强度图像中模糊连通强度高 于模糊连通强度阈值且与种子点相连的连通区域作为血管内腔区域。
【文档编号】G06T7/00GK104331881SQ201410566443
【公开日】2015年2月4日 申请日期:2014年10月23日 优先权日:2014年10月23日
【发明者】严加勇, 崔崤峣, 向永嘉, 韩志乐, 简小华 申请人:中国科学院苏州生物医学工程技术研究所