基于亮度和梯度的点状钙化区域检测方法及系统

摘要

本发明公开了一种基于亮度和梯度的点状钙化区域检测方法，包括：S1、对血管区域图像进行图像处理，获得钙化候选区域；S2、对钙化候选区域进行两次膨胀处理，获得一次膨胀区域与二次膨胀区域；S3、获得钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；S4、获得包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点的序列点集；S5、对序列点集进行亮度和梯度分析，获得符合筛选条件的序列点集；S6、求解标记为符合的序列点集与总序列点集数的比值，以将钙化候选区域确认为钙化区域。同时，本发明还公开了一种基于亮度和梯度的点状钙化区域检测系统。本发明对于点状钙化区域有准确的判定效果。

权利要求

1.基于亮度和梯度的点状钙化区域检测方法，其特征在于，包括：S1、对血管区域图像进行图像处理，获得钙化候选区域；

S2、对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；

S3、获得钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；

S4、序列点集的获得：

S41、按设定步长选取钙化候选区域轮廓线的轮廓点；

S42、对选取的各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对；

S43、以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，所述序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；

S5、对序列点集进行亮度和梯度分析，获得符合筛选条件的序列点集：S51、对各组序列点集中的轮廓点，分析其亮度是否高于设定阈值，若是，执行S52，若否，则标记为不符合；

S52、分析序列点集中的中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点是否满足梯度衰减，若是则标记为符合，若否，则标记为不符合；

S6、求解标记为符合的序列点集与总序列点集数的比值，当比值超过预设值，将该钙化候选区域确认为钙化区域。

2.基于亮度和梯度的点状钙化区域检测系统，其特征在于，包括：血管区域图像生成模块，其对输入图像进行处理，生成血管区域图像；

钙化候选区域生成模块，其对所述血管区域图像进行处理，获得钙化候选区域；

钙化候选区域预处理模块，其对所述钙化候选区域进行如下处理：对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；提取钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；计算钙化候选区域的面积；

序列点集生成模块，其基于钙化候选区域预处理模块的预处理结果，对钙化候选区域轮廓线上的各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对，以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，各序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；

亮度、梯度分析模块，其对各组序列点集中的轮廓点分析亮度是否大于预设阈值，若是，分析该序列点集中的中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点是否梯度衰减，将符合条件的点标记为符合，输出标记为符合的序列点集与总序列点集数的比值；

判定模块，其基于所述比值及所述钙化候选区域预处理模块输出的钙化候选区域的面积，判定当前钙化候选区域是否为钙化区域。

说明书

基于亮度和梯度的点状钙化区域检测方法及系统

技术领域

[0001] 本发明涉及冠脉医学图像处理领域，具体涉及一种基于亮度和梯度的点状钙化区域检测方法及系统。

背景技术

[0002] 自动化冠脉医学图像处理对医生具有重要的临床价值和实际意义，其能为医生反馈直观的检测结果，从而作为医生进行病情诊断的参考，将医生从解读医学图像的繁琐工作中解脱出来，降低了医生的诊断时间，提高诊断效率，缓解当前就医难的问题。

[0003] 钙化区域的识别是自动化冠脉医学图像处理中的重要一环，钙化一般在医学图像上的表现形式通常是其亮度值高于周围血管的亮度值，据此，大部分算法会设定一个阈值或动态阈值进行区分，进而识别出钙化区域。

[0004] 此种方法对于块状钙化区域，检测效果良好，但对于点状钙化区域，由于其周围的干扰因素多，常常出现误检。

发明内容

[0005] 本发明的目的在于提供一种基于亮度和梯度的点状钙化区域检测方法及系统。

[0006] 为实现上述目的，本发明采用以下技术方案：

[0007] 基于亮度和梯度的点状钙化区域检测方法，包括：

[0008] S1、对血管区域图像进行图像处理，获得钙化候选区域；

[0009] S2、对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；

[0010] S3、获得钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；

[0011] S4、序列点集的获得：

[0012] S41、按设定步长选取钙化候选区域轮廓线的轮廓点；

[0013] S42、对各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对；

[0014] S43、以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，所述序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；

[0015] S5、对序列点集进行亮度和梯度分析，获得符合筛选条件的序列点集：

[0016] S51、对各组序列点集中的轮廓点，分析其亮度是否高于设定阈值，若是，执行S52，若否，则标记为不符合；

[0017] S52、分析序列点集中的中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点是否满足梯度衰减，若是则标记为符合，若否，则标记为不符合；

[0018] S6、求解标记为符合的序列点集与总序列点集数的比值，当比值超过预设值，将该钙化候选区域确认为钙化区域。

[0019] 进一步地，所述预设值与所述钙化候选区域的面积相关，则有：

[0020]

[0021] 式中，A为钙化候选区域的面积，A1为面积阈值，R为预设值。

[0022] 本发明还公开了一种基于亮度和梯度的点状钙化区域检测系统，包括：

[0023] 血管区域图像生成模块，其对输入图像进行处理，生成血管区域图像；

[0024] 钙化候选区域生成模块，其对所述血管区域图像进行处理，获得钙化候选区域；

[0025] 钙化候选区域预处理模块，其对所述钙化候选区域进行如下处理：对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；提取钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；计算钙化候选区域的面积；

[0026] 序列点集生成模块，其基于钙化候选区域预处理模块的预处理结果，对钙化候选区域轮廓线上的各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对，以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，各序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；

[0027] 亮度、梯度分析模块，其对各组序列点集中的轮廓点分析亮度是否大于预设阈值，若是，分析该序列点集中的中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点是否梯度衰减，将符合条件的点标记为符合，输出标记为符合的序列点集与总序列点集数的比值；

[0028] 判定模块，其基于所述比值及所述钙化候选区域预处理模块输出的钙化候选区域的面积，判定当前钙化候选区域是否为钙化区域。

[0029] 采用上述技术方案后，本发明与背景技术相比，具有如下优点：本发明基于对钙化候选区域的膨胀操作，得出钙化候选区域及其膨胀区域的梯度变化，分析其是否符合衰减规律，以亮度和梯度两个指标综合判定其是否为真实钙化区域。其检测手段简单、图像计算量小，同时，对于点状钙化区域有准确的判定效果。

附图说明

[0030] 图1为血管与钙化候选区域的示意图；

[0031] 图2为钙化候选区域的中心线以及各区域的轮廓线的示意图；

[0032] 图3为序列点集中各点的示意图；

[0033] 图4为本发明系统组成框图。

具体实施方式

[0034] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

[0035] 在本发明中需要说明的是，术语“上”“下”“左”“右”“竖直”“水平”“内”“外”等均为基于附图所示的方位或位置关系，仅仅是为了便于描述本发明和简化描述，而不是指示或暗示本发明的装置或元件必须具有特定的方位，因此不能理解为对本发明的限制。

[0036] 实施例1

[0037] 本发明公开了一种基于亮度和梯度的点状钙化区域检测方法，核心包括5个步骤。

[0038] S1、对血管区域图像进行图像处理，获得钙化候选区域；

[0039] S1中，血管区域图像是以拉直图像与拉直分割图像作为输入而获取的。钙化候选区域可以通过基本阈值、对比度或极值等方法获取。具体示意图可参照图1所示。

[0040] S2、对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域。

[0041] S3、通过收形算法(如骨架提取方法)获得钙化候选区域的中心线，通过边缘算法获得钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线。其处理结果如图2所示，线1为中心线，线2为钙化候选区域轮廓线，线3为一次膨胀区域轮廓线，线4为二次膨胀区域轮廓线。

[0042] S4、序列点集的获得：

[0043] S41、按设定步长(可以为0，也可以为2-3个像素点)选取钙化候选区域轮廓线的轮廓点；

[0044] S42、对各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对；

[0045] S43、以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，所述序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；

[0046] 其处理结果如图3所示，点1为中心点、点2为轮廓点、点3为一次膨胀轮廓点、点4为二次膨胀轮廓点。

[0047] S5、对序列点集进行亮度和梯度分析，获得符合筛选条件的序列点集：

[0048] S51、对各组序列点集中的轮廓点，分析其亮度是否高于设定阈值，若是，执行S52，若否，则标记为不符合；

[0049] S52、分析序列点集中的中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点是否满足梯度衰减，若是则标记为符合，若否，则标记为不符合；

[0050] S6、求解标记为符合的序列点集与总序列点集数的比值，当比值超过预设值，将该钙化候选区域确认为钙化区域。

[0051] 其中，所述预设值与所述钙化候选区域的面积相关，则有：

[0052]

[0053] 式中，A为钙化候选区域的面积，A1为面积阈值，R为预设值。本实施例中，A1的取值范围为35～50个像素点。

[0054] 实施例2

[0055] 请参考图4所示，本发明还公开了一种基于亮度和梯度的点状钙化区域检测系统，包括血管区域图像生成模块、钙化候选区域生成模块、钙化候选区域预处理模块、序列点集生成模块、亮度、梯度分析模块及判定模块。

[0056] 血管区域图像生成模块，其对输入图像进行处理，生成血管区域图像；

[0057] 钙化候选区域生成模块，其对所述血管区域图像进行处理，获得钙化候选区域；

[0058] 钙化候选区域预处理模块，其对所述钙化候选区域进行如下处理：对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；提取钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；计算钙化候选区域的面积。

[0059] 序列点集生成模块，其基于钙化候选区域预处理模块的预处理结果，对钙化候选区域轮廓线上的各轮廓点(按设定步长选取)寻找距其最近的中心点，获得中心点-轮廓点对，以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，各序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点。

[0060] 亮度、梯度分析模块，其对各组序列点集中的轮廓点分析亮度是否大于预设阈值，若是，分析该序列点集中的中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点是否梯度衰减，将符合梯度衰减的点标记为符合，输出标记为符合的序列点集与总序列点集数的比值。

[0061] 判定模块，其基于所述比值及所述钙化候选区域预处理模块输出的钙化候选区域的面积，判定当前钙化候选区域是否为钙化区域。

[0062] 具体地，当比值超过预设值，将该钙化候选区域确认为钙化区域。所述预设值与所述钙化候选区域的面积相关，则有：

[0063]

[0064] 式中，A为钙化候选区域的面积，A1为面积阈值，R为预设值。

[0065] 以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。