本发明公开了一种钙化检测方法，包括以下步骤：获取血管区域图像；利用钙化检测算法获得钙化候选区域；基于梯度和亮度分析，检出点状钙化区域；基于亮度分析判断点状钙化区域是否为真钙化；钙化边界修正。本发明能够有效检测出血管图像上的点状钙化区域，通过形态分析剔除误报区域，有效避免了漏检、误报情况的发生。
S31、通过对钙化候选区域的膨胀操作，得出钙化候选区域及其膨胀区域的梯度变化，分析其是否符合衰减规律，结合亮度和梯度两个指标判定其是否为真实钙化区域；
S32、基于血管区域图像上像素点的亮度值排序，多次选取不同数量的像素点并进行区域分割，然后通过区域连续增长分析判断出是否为候选区域，最后分析候选区域与周围区域的亮度值差异，判断出是否为真钙化；
S311、对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；
S312、获得钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；
S313、获取序列点集，按设定步长选取钙化候选区域轮廓线的轮廓点，对选取的各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对，以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，所述序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；
S314、对序列点集进行亮度和梯度分析，获得符合筛选条件的序列点集；
S315、求解标记为符合的序列点集与总序列点集数的比值，当比值大于或等于预设值，将该钙化候选区域确认为钙化区域。
2.如权利要求1所述的一种钙化检测方法，其特征在于，其还包括步骤S4,所述步骤S4具体为：基于亮度分析判断点状钙化区域是否为真钙化。
3.如权利要求2所述的一种钙化检测方法，其特征在于，所述步骤S4通过以下分步骤实现：S41、针对点状钙化区域进行形态分析，探测是否存在小区域钙化发生在血管分叉处情况；
S42、对小区域钙化发生在血管分叉处情况，提升小区域钙化检测阈值，若仍被检测出钙化，则判定为真钙化。
4.如权利要求1所述的一种钙化检测方法，其特征在于，所述步骤S32包括以下子步骤：S321、统计血管区域图像中的亮度值，分别选取亮度值最高的M1个像素点、M2个像素点、M3个像素点，其中M1
S322、对所述候选区域与其前后左右区域进行对比分析，若存在弱对比则判定为假钙化，若存在强对比则判定为钙化区域；
计算钙化候选区域前后左右四个区域的亮度值均值并分别记为P1、P2、P3、P4，计算钙化候选区域的亮度值均值并记为V；分别计算P1、P2、P3、P4与V的差值，若符合至少三个差值大于预先设定的阈值，则存在强对比并判定为钙化区域，若少于三个差值大于预先设定的阈值，则存在弱对比并判定为假钙化。
5.如权利要求1所述的一种钙化检测方法，其特征在于，其还包括步骤S5,所述步骤S5具体为：通过边界进行多灰度级定义，引入中间灰度级进行边界重新计算，输出钙化分割区域。
6.如权利要求5所述的一种钙化检测方法，其特征在于，所述步骤S5通过以下分步骤实现：S51、对所述血管区域图像进行二值化分割，获得二值化图像；
S52、以3个灰度值对图像分割边界进行定义，分别为0、中间灰度级、255，中间灰度级为灰度值分别为0和255的两个像素点之间的调整值；
S53、针对二值化图像的边界，获取灰度值由0跃到255的两个像素点，灰度值为0的像素点记为P，灰度值为255的像素点记为Q，以像素点Q作为参考像素，将其变为灰度像素，参考像素的值即为所述中间灰度级，所述中间灰度级通过下公式求取：其中，V为中间灰度级，将边界内侧与像素点Q最相邻的灰度值为255的像素点记为R，L2为像素点R的原始灰度值，L1为像素点P的原始灰度值；
7.如权利要求1所述的一种钙化检测方法，其特征在于，所述步骤S2还包括：对所述钙化候选区域进行形态分析，剔除形态异常的钙化候选区域。
一种钙化检测方法
技术领域
[0001] 本发明涉及冠脉医学图像处理技术领域，特别涉及一种钙化检测方法。
背景技术
[0002] 自动化冠脉医学图像检测对医生而言具有重要的临床价值和实际意义，其能够能为医生反馈直观的检测结果，从而作为医生进行病情诊断的参考，将医生从解读医学图像的繁琐工作中解脱出来，从而降低了医生的诊断时间，提高诊断效率，缓解当前就医难的问题。
[0003] 钙化区域识别是自动化冠脉医学图像检测中的重要一环，钙化一般在医学图像上的表现形式通常是其亮度值高于周围血管的亮度值，据此，现有的算法大多通过设定一个固定阈值或动态阈值进行区分，进而识别出钙化区域。对于块状钙化区域，检测效果良好，但对于点状钙化区域，由于其周围的干扰因素多，容易出现漏检。另外对应一些复杂情况，也容易出现误报情况，例如在血管分叉处由于灰度值的变化，也容易将正常血管判断为小区域钙化。
发明内容
[0004] 为解决上述问题，本发明提供了一种钙化检测方法。
[0005] 本发明采用以下技术方案：
[0006] 一种钙化检测方法，包括以下步骤：
[0007] S1、获取血管区域图像；
[0008] S2、利用钙化检测算法获得钙化候选区域；
[0009] S3、基于梯度和亮度分析，检出点状钙化区域。
[0010] 优选地，其还包括步骤S4,所述步骤S4具体为：基于亮度分析判断点状钙化区域是否为真钙化。
[0011] 优选地，所述步骤S4通过以下分步骤实现：
[0012] S41、针对点状钙化区域进行形态分析，探测是否存在小区域钙化发生在血管分叉处情况；
[0013] S42、对小区域钙化发生在血管分叉处情况，提升小区域钙化检测阈值，若仍被检测出钙化，则判定为真钙化。
[0014] 优选地，所述步骤S3通过以下分步骤实现：
[0015] S31、通过对钙化候选区域的膨胀操作，得出钙化候选区域及其膨胀区域的梯度变化，分析其是否符合衰减规律，结合亮度和梯度两个指标判定其是否为真实钙化区域；
[0016] S32、基于血管区域图像上像素点的亮度值排序，多次选取不同数量的像素点并进行区域分割，然后通过区域连续增长分析判断出是否为候选区域，最后分析候选区域与周围区域的亮度值差异，判断出是否为真钙化。
[0017] 优选地，所述步骤S31包括以下子步骤：
[0018] S311、对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；
[0019] S312、获得钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；
[0020] S313、获取序列点集，按设定步长选取钙化候选区域轮廓线的轮廓点，对选取的各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对，以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，所述序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；
[0021] S314、对序列点集进行亮度和梯度分析，获得符合筛选条件的序列点集；
[0022] S315、求解标记为符合的序列点集与总序列点集数的比值，当比值超过预设值，将该钙化候选区域确认为钙化区域。
[0023] 优选地，所述预设值与所述钙化候选区域的面积相关，则有：
[0024]
[0025] 式中，A为钙化候选区域的面积，A1为面积阈值，R为预设值。
[0026] 优选地，所述步骤S32包括以下子步骤：
[0027] S321、统计血管区域图像中的亮度值，分别选取亮度值最高的M1个像素点、M2个像素点、M3个像素点，其中M1
[0028] S322、对所述候选区域与其前后左右区域进行对比分析，若存在弱对比则判定为假钙化，若存在强对比则判定为钙化区域。
[0029] 优选地，其还包括步骤S5,所述步骤S5具体为：通过边界进行多灰度级定义，引入中间灰度级进行边界重新计算，输出钙化分割区域。
[0030] 优选地，所述步骤S5通过以下分步骤实现：
[0031] S51、对所述血管区域图像进行二值化分割，获得二值化图像；
[0032] S52、以3个灰度值对图像分割边界进行定义，分别为0、中间灰度级、255，中间灰度级为灰度值分别为0和255的两个像素点之间的调整值；
[0033] S53、针对二值化图像的边界，获取灰度值由0跃到255的两个像素点，灰度值为0的像素点记为P，灰度值为255的像素点记为Q，以像素点Q作为参考像素，将其变为灰度像素，参考像素的值即为所述中间灰度级，所述中间灰度级通过下公式求取：
[0034]
[0035] 其中，V为中间灰度级，将边界内侧与像素点Q最相邻的灰度值为255的像素点记为R，L2为像素点R的原始灰度值，L1为像素点P的原始灰度值；
[0036] S54、基于中间灰度级计算边界宽度值；
[0037] S55、基于边界宽度值，输出钙化分割区域。
[0038] 优选地，所述步骤S2还包括：对所述钙化候选区域进行形态分析，剔除形态异常的钙化候选区域。
[0039] 采用上述技术方案后，本发明与背景技术相比，具有如下优点：
[0040] 本发明能够有效检测出血管图像上的点状钙化区域，通过形态分析剔除误报区域，有效避免了漏检、误报情况的发生。另外，通过对血管拉直图像上的边界进行多灰度级定义（以三值进行定义），引入中间灰度级来实现边界小数级计算，实现钙化边界修正，从而提升了冠脉狭窄度计算结果的准确性。
附图说明
[0041] 图1为本发明的流程示意图。
具体实施方式
[0042] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。实施例
[0043] 参考图1所示，本发明公开了一种钙化检测方法，包括以下步骤：
[0044] S1、获取血管区域图像。血管区域图像通过对原始血管提供进行拉直、分割获得。
[0045] S2、利用钙化检测算法获得钙化候选区域。钙化检测算法可采用基本阈值、对比度或极值算法。为便于后续处理，基于先验经验剔除形态上不符合钙化特征的钙化候选区域结果，通过对钙化候选区域域进行形态分析，剔除形态异常的钙化候选区域，这里提及的“形态异常”指的是基于先验经验在形态上不符合钙化特征，例如钙化候选区域在垂直血管方向上的宽度与沿血管方向上的宽度的比值较大，则认为是形态异常。
[0046] S3、基于梯度和亮度分析，检出点状钙化区域。
[0047] S31、通过对钙化候选区域的膨胀操作，得出钙化候选区域及其膨胀区域的梯度变化，分析其是否符合衰减规律，结合亮度和梯度两个指标判定其是否为真实钙化区域。
[0048] S311、对钙化候选区域进行一次膨胀与二次膨胀处理，获得一次膨胀区域与二次膨胀区域；
[0049] S312、获得钙化候选区域的中心线以及钙化候选区域轮廓线、一次膨胀区域轮廓线及二次膨胀区域轮廓线；
[0050] S313、获取序列点集，按设定步长（可以为0，也可以为2-3个像素点）选取钙化候选区域轮廓线的轮廓点，对选取的各轮廓点寻找距其最近的中心点，获得中心点-轮廓点对，以中心点-轮廓点对作穿过各膨胀区域的射线，获得多组序列点集，所述序列点集依次包括中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点；
[0051] S314、对序列点集进行亮度和梯度分析，获得符合筛选条件的序列点集。具体为：
[0052] a. 对各组序列点集中的轮廓点，分析其亮度是否高于设定阈值，若是，执行步骤b，若否，则标记为不符合；
[0053] b、分析序列点集中的中心点、轮廓点、一次膨胀轮廓点及二次膨胀轮廓点是否满足梯度衰减，若是则标记为符合，若否，则标记为不符合.
[0054] S315、求解标记为符合的序列点集与总序列点集数的比值，当比值超过预设值，将该钙化候选区域确认为钙化区域。
[0055] 在本实施例中，所述预设值与所述钙化候选区域的面积相关，则有：
[0056]
[0057] 式中，A为钙化候选区域的面积，A1为面积阈值，R为预设值。本实施例中，A1的取值范围为35 50个像素点。~
[0058] S32、基于血管区域图像上像素点的亮度值排序，多次选取不同数量的像素点并进行区域分割，然后通过区域连续增长分析判断出是否为候选区域，最后分析候选区域与周围区域的亮度值差异，判断出是否为真钙化。
[0059] S321、统计血管区域图像中的亮度值，分别选取亮度值最高的M1个像素点、M2个像素点、M3个像素点，其中M1
[0060] 统计血管区域图像中的亮度值，分别选取亮度值最高的M1个像素点、M2个像素点、M3个像素点，其中M1 M2区块的面积> M1区块的面积,则认定存在增长，M1区块作为钙化候选区域。
[0061] S322、对所述候选区域与其前后左右区域进行对比分析，若存在弱对比则判定为假钙化，若存在强对比则判定为钙化区域。该步骤具体为：
[0062] 计算钙化候选区域前后左右四个区域的亮度值均值并分别记为P1、P2、P3、P4，计算钙化候选区域的亮度值均值并记为V；分别计算P1、P2、P3、P4与V的差值，若符合至少三个差值大于预先设定的阈值，则存在强对比并判定为钙化区域，若少于三个差值大于预先设定的阈值，则存在弱对比并判定为假钙化。
[0063] 在本实施例中，为避免出现误报情况，在步骤S322之后还进行以下操作：计算钙化区域的尺寸，若存在大面积、纵向过长或横向过长则将其滤除。这样通过相应的形态分析，可以确保点状钙化区域检测过程中异常结果的误报。
[0064] S4、基于亮度分析判断点状钙化区域是否为真钙化。
[0065] S41、针对点状钙化区域进行形态分析，探测是否存在小区域钙化发生在血管分叉处情况。
[0066] S42、对小区域钙化发生在血管分叉处情况，提升小区域钙化检测阈值，若仍被检测出钙化，则判定为真钙化。
[0067] S5、通过边界进行多灰度级定义，引入中间灰度级进行边界重新计算，输出钙化分割区域。
[0068] S51、对所述血管区域图像进行二值化分割，获得二值化图像。
[0069] S52、以3个灰度值对图像分割边界进行定义，分别为0、中间灰度级、255，中间灰度级为灰度值分别为0和255的两个像素点之间的调整值。
[0070] S53、针对二值化图像的边界，获取灰度值由0跃到255的两个像素点，灰度值为0的像素点记为P，灰度值为255的像素点记为Q，以像素点Q作为参考像素，将其变为灰度像素，参考像素的值即为所述中间灰度级，所述中间灰度级通过下公式求取：
[0071]
[0072] 其中，V为中间灰度级，将边界内侧与像素点Q最相邻的灰度值为255的像素点记为R，L2为像素点R的原始灰度值，L1为像素点P的原始灰度值。
[0073] S54、基于中间灰度级计算边界宽度值。在本实施例中，边界宽度值的计算公式为：
[0074]
[0075] 其中，W为边界宽度值，V为中间灰度级。
[0076] S55、基于边界宽度值，输出钙化分割区域。这样，可以方便计算后续钙化狭窄程度。
[0077] 以上，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
