一种基于毗邻图像特征点集的dsa脑血管图像自动分割方法
【专利摘要】本发明公开了一种基于毗邻图像特征点集的DSA脑血管图像自动分割方法，包括1：导入若干对连续的DSA脑血管图像作为源图像数据；2：对每一对DSA脑血管图像均进行分区；3：对分区后的DSA脑血管设置图像阈值；4：基于sift算法提取特征点；5：对每一对提取特征点的DSA脑血管图像中的蒙片图像与活片图像均基于数字减影技术得到相应活片图像的特征点差值图像；6：提取所有所述特征点差值图像的图像特征点集，利用相邻图像关系对图像特征点集进行精确提取；7：对提取的图像特征点集进行区域生长，得到相应的脑血管图像。本发明采用图像分割技术，利用毗邻图像，将邻域间的像素信息进行综合，提高了特征点信息提取的准确性，有效减小了噪声。
【专利说明】一种基于毗邻图像特征点集的DSA脑血管图像自动分割方 法
【技术领域】
[0001] 本发明涉及医学图像处理领域，具体地说是涉及一种基于毗邻图像特征点集的 DSA脑血管图像自动分割方法。
【背景技术】
[0002] 在利用DSA机进行脑血管畸形和脑血管瘤的诊断与手术时，通常会遇到DSA图像 中伪影剧烈的困难，从而严重影响了脑血管诊断与治疗的可靠性。为了便于诊断病灶，得到 患者脑血管中的三维信息。其中通常需要减影操作，采集两个脑血管图像序列，即蒙片序列 图像和活片图像序列；两者对应相减，理想情况下应该得到只含有血管的减影图像。但由于 各种因素造成的误匹配经常使得所获取的造影图片的诊断大打折扣。在这样的条件下，就 需要一种快速有效配准算法来消除伪影的影响，从而得到高质量的脑血管减影图像，提高 脑血管疾病诊断和治疗的可靠性。
[0003] 目前出现的关于脑血管自动分割方法大部分是基于配准的方法，而医学图像配准 过程的运算量非常大，即使采用刚性变换，对于大规模高分辨率断层，图像数据集包含的庞 大数据量，也无法满足临床上实时处理的要求，因而必须采取优化措施；同时对于减影技术 的运动伪影的消除，目前尚没有能够大范围应用到临床医学中的方法。
[0004] 具体的，现有的常用脑血管自动分割方法：
[0005] 1 :人工进行配准，操作人员手动进行配准，其存在效率低下，图像数据集包含的数 据量十分庞大，无法满足临床上实时处理的要求等问题；
[0006] 2 :模块匹配自动配准技术，其中有框架的配准方法基于外部基准点特征，能够获 得较高的精度，可作为评估无框架配准算法的标准，但其植入式的特点会给患者带来很大 的痛苦，同时也不易对历史图像作回溯式的研宄。
[0007] 3 :采用基于规则网格模型的算法自动生成序列控制点，根据控制点空间排布的特 点，优化搜索相关点的过程，同时限制控制点在二维空间的搜索范围，然后采用基于逆拉伸 空间变换的像素映射填充算法来生成变形后的蒙片目标图像，实现减影，但是其存在算法 的鲁棒性、计算速度和正确性不够健全等缺陷。
【发明内容】
[0008] 鉴于已有技术存在的缺陷，本发明的目的是要提供一种基于毗邻图像特征点集的 DSA脑血管图像自动分割方法，该方法具有采用图像分割技术，利用毗邻图像，将邻域间的 像素信息进行综合，提高了特征点信息提取的准确性，有效减小了噪声。
[0009] 为了实现上述目的，本发明的技术方案：
[0010] -种基于毗邻图像特征点集的DSA脑血管图像自动分割方法，其特征在于：
[0011] 包括
[0012] 步骤1 :导入若干对连续的DSA脑血管图像作为源图像数据；
[0013] 步骤2 :对每一对DSA脑血管图像均进行分区，使得每一对DSA脑血管图像中的蒙 片图像与活片图像均分为上下两个区域；
[0014] 步骤3 :对步骤2中分区后的DSA脑血管设置图像阈值：对DSA脑血管图像的蒙片 图像与活片图像中的上区域与下区域分别设置不同图像阈值；
[0015] 步骤4 :基于sift算法提取特征点：对步骤3中设置图像阈值后的每一对DSA脑 血管图像中的蒙片图像与活片图像均基于sift算法，按照步骤3中所设置的上下区域不同 阈值提取特征点，得到对应的图像特征点集；
[0016] 步骤5 :对步骤4中每一对提取特征点的DSA脑血管图像中的蒙片图像与活片图 像均基于数字减影技术得到相应活片图像的特征点差值图像；
[0017] 步骤6 :提取步骤5中所有所述特征点差值图像的图像特征点集，利用相邻图像关 系对图像特征点集进行精确提取；
[0018] 步骤7 :对步骤6中提取的图像特征点集进行区域生长，得到相应的脑血管图像。
[0019] 进一步的，所述的基于数字减影技术得到相应活片图像的特征点差值图像是指 ①、提取蒙片图像的图像特征点集VI、活片图像的图像特征点集V2后；②、遍历蒙片图像的 图像特征点集Vl中的各个特征点，并在活片图像上找到相应特征点对应的位置坐标；③、 以所述位置坐标为圆心作圆后，遍历活片图像的图像特征点集V2,将处于所述圆的圆内且 包含圆上特征点去除，仅保留处于所述圆的圆外的特征点；④、将保留的所述特征点映射到 活片图像对应位置上，得到特征点集差值图像。
[0020] 进一步的所述利用相邻图像关系对图像特征点集进行精确提取是指：①、提取所 有映射到活片图像上的差值特征点集；②、依次将相邻三幅活片图像做以下步骤：i、所述 相邻三幅活片图像的差值特征点集分别设定为T1、T2、T3 ;ii、遍历Tl中的各个特征点；iii、 在T2中找到与ii中所述特征点对应的位置坐标，并以其为圆心作圆，保留所述圆圆内的特 征点，得到Tl、T2共同特征点集T4 ;iv、对T2、T3重复ii、iii步骤得到T2、T3共同特征点集 T5后，对共同特征点集T4、T5求并集，得到综合特征点集T6。
[0021] 与现有技术相比，本发明的有益效果：
[0022] 1.本发明不同于以往大多基于配准方法的研宄，采用图像分割技术，利用分区选 取不同的阈值的做法，使血管特征点的提取更加精确。
[0023] 2.本发明利用毗邻图像，将邻域间的像素信息进行综合，提高了特征点信息提取 的准确性，减小了噪声。
[0024] 3.本发明为自动分割，操作简单，不同操作人员得到的分割结果基本相同。
[0025] 4.本发明编程容易实现，算法复杂度较低。
【专利附图】
【附图说明】
[0026] 图1为本发明-基于毗邻图像特征点集的DSA脑血管图像自动分割方法步骤流程 示意图；
[0027] 图2为本发明-蒙片图像提取特征点示意图；
[0028] 图3为本发明-活片图像提取特征点示意图；
[0029] 图4为本发明-分区选取阈值后的提取特征点示意图；
[0030] 图5为本发明-特征点差值图像示意图；
[0031]图6为本发明-根据毗邻图像特征点集差值图像提取特征点示意图；
[0032]图7为本发明-步骤6中所述相邻三幅图像中的第n-1幅图像特征点集差值图像 示意图；
[0033] 图8为本发明-步骤6中所述相邻三幅图像中的第η幅图像特征点集差值图像示 意图；
[0034] 图9为本发明-步骤6中所述相邻三幅图像中的第η+1幅图像特征点集差值图像 示意图；
[0035] 图10为本发明-本地方窗与DSA图像上的移动方向示意图；
[0036] 图11为本发明-步骤7中的区域成长示意图；
[0037] 图12为本发明-区域成长完成后的完整脑血管图像示意图。
【具体实施方式】
[0038] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进 行进一步详细说明。
[0039] 如图1-图11所示，一种基于毗邻图像特征点集的DSA脑血管图像自动分割方法， 其特征在于：
[0040] 如图1所示，包括
[0041] 步骤1 :导入若干对连续的DSA脑血管图像作为源图像数据。
[0042] 步骤2 :对每一对DSA脑血管图像均进行分区，使得每一对DSA脑血管图像中的蒙 片图像与活片图像均分为上下两个区域。
[0043] 进一步的，所述分区是指将全部DSA脑血管图像（蒙片图像与活片图像）均按照 自上至下从所述图像的三分之一处水平分割成上下两个区域；之所以选取三分之一处是基 于实验观察发现脑血管蒙片活片图像在三分之一以上的脑血管较细且分布比较密集；下面 三分之二，血管较大且分布较稀疏的特点。
[0044] 步骤3 :对步骤2中分区后的DSA脑血管设置图像阈值：对DSA脑血管图像的蒙片 图像与活片图像中的上区域与下区域分别设置不同图像阈值，并使得全部的上区域以及下 区域统一采用相同的阈值尺度。即所述蒙片图像与活片图像中的上区域与下区域的图像阈 值应当尽可能选取相差较大的两个阈值，本例中所有蒙片图像与活片图像中的上区域的阈 值使用100,所有蒙片图像与活片图像中的下区域的阈值使用8。
[0045] 步骤4 :利用上述阈值策略，基于sift算法按照步骤3中所设置的上下区域不同 阈值提取特征点，得到特征点分布更为均匀的脑血管图像：对步骤3中设置图像阈值后的 每一对DSA脑血管图像中的蒙片图像与活片图像均基于sift算法提取特征点，得到对应的 图像特征点集。
[0046] 在图像处理中，特征点可以称兴趣点或者角点，三者经常相互使用，即图像的极值 点，线段的终点，曲线曲率最大的点或者水平或者竖直方向上属性最大的点等等，这些特征 点是图像很重要的特征，对图像图形的理解和分析有很重要的作用；因此本方法首先利用 常用的sift算法（包括尺度空间的生成；检测尺度空间极值点；精确定位极值点；为每个 关键点指定方向参数；关键点描述子的生成等，为领域内常用算法，这里不再赘述）对每一 对DSA脑血管图像中的蒙片图像与活片图像进行特征点提取，得到对应的图像特征点集， 即对应每一对的蒙片图像的图像特征点集VI、活片图像的图像特征点集V2。
[0047] 步骤5:对步骤4中每一对提取特征点的DSA脑血管图像中的蒙片图像与活片图 像均基于数字减影技术得到相应活片图像的特征点差值图像。
[0048] 所述的基于数字减影技术得到相应活片图像的特征点差值图像是指①、提取蒙片 图像的图像特征点集VI、活片图像的图像特征点集V2后；②、遍历蒙片图像的图像特征点 集Vl中的各个特征点，并在活片图像上找到相应特征点对应的位置坐标；③、以所述位置 坐标为圆心作圆后，遍历活片图像的图像特征点集V2,将处于所述圆的圆内且包含圆上特 征点去除，仅保留处于所述圆的圆外的特征点；④、将保留的特征点映射到活片图像对应位 置上，得到特征点集差值图像。
[0049] 如，提取蒙片图像的图像特征点集VI、活片图像的图像特征点集V2后；遍历蒙片 图像的图像特征点集Vl中的各个特征点，以特征点Pl为例，并在活片图像上找到特征点Pl 对应的位置坐标；以特征点Pl对应的位置坐标为圆心并选取合适半径rl作圆，其中半径 的选取则是通过不断的选取值实验，观察最后的图像中的特征点，若噪声特征点提取较少， 想要的血管上的特征点大部分都提取到即选择最好的情况下的值作为半径，选择无具体规 贝1J，以达到较好效果即可；遍历活片图像的图像特征点集V2,将V2中处于所述圆的圆内且 包含圆上特征点去除，仅保留处于圆外的特征点；将保留的特征点映射到活片图像对应位 置上，得到特征点集差值图像，上述过程可参看如图2-图6。
[0050] 步骤6 :提取步骤5中所有所述特征点差值图像的图像特征点集，利用相邻图像关 系对图像特征点集进行精确提取；
[0051] 所述利用相邻图像关系对图像特征点集进行精确提取是指：①、提取所有映射到 活片图像上的差值特征点集；②、依次分别将相邻三幅活片图像做以下步骤：i、所述相邻 三幅活片图像的差值特征点集分别设定为Tl、T2、T3 ;ii、遍历Tl中的各个特征点；iii、在 T2中找到与ii中所述特征点对应的位置坐标，并以其为圆心作圆，保留所述圆圆内的特征 点，得到Tl、T2共同特征点集T4 ;iv、对T2、T3重复ii、iii步骤即遍历T2中的各个特征点， 在T3中找到与ii中所述特征点对应的位置坐标，并以其为圆心作圆，保留所述圆圆内的特 征点，得到T2、T3共同特征点集T5 5后，对共同特征点集T4、T5求并集，得到综合特征点集 Τ6〇
[0052] 如，提取所有映射到活片图像上的差值特征点集，依次分别将相邻三幅图像做以 下步骤，以第η-1、η、η+1幅为例，如图7-图9所示，它们的特征点集分别为Tl、Τ2、Τ3 ;遍 历Tl中的特征点，以特征点nl为例，选取合适半径r2,作圆，半径选取原则与步骤5中的半 径选取原则一样；若T2中作圆后圆的内部有特征点则保留，否则除去；得到Tl、T2共同特 征点集T4 ;同样的将T2、T3做上述步骤，得到T2、T3共同特征点集T5 ;将两个共同特征点 集T4、T5求交得到综合特征点集T6。
[0053] 步骤7 :对步骤6中提取的综合图像特征点集进行区域生长，得到相应的脑血管图 像，如图12所示。
[0054] 如图11所示，区域生长为血管的详细过程包括三个步骤：
[0055] 第1步：初始化所有特征点作为种子像素集。标记所有这些种子像素点为"种子 像素未生长"的状态，标记图像中的其他像素是"普通像素未生长"的状态。
[0056] 第2步：搜索种子像素集中"种子像素未生长"状态的种子像素，标记状态为"种子 像素已生长"，转到步骤3。如果没有一个种子像素为"种子像素未生长"的状态时，整个遍 历处理结束。
[0057] 第3步：遍历种子像素的八邻域，如果它的灰度值比它的区域阈值g小，并且它的 状态为"普通像素未生长"，把它加入种子像素集，并标记状态为"种子像素未生长"，否则， 标记状态为"普通像素已生长"，然后遍历下一个像素，当八邻域遍历结束，转到第2步。
[0058] 为了得到所述区域阈值d，对每个所述种子像素都需要执行如下三个步骤：
[0059] 步骤1 :建立本地方形窗口宽度为w
[0060] 步骤2:计算阈值如下
【权利要求】
1. 一种基于毗邻图像特征点集的DSA脑血管图像自动分割方法，其特征在于： 包括 步骤1 :导入若干对连续的DSA脑血管图像作为源图像数据； 步骤2 :对每一对DSA脑血管图像均进行分区，使得每一对DSA脑血管图像中的蒙片图 像与活片图像均分为上下两个区域； 步骤3 :对步骤2中分区后的DSA脑血管设置图像阈值：对DSA脑血管图像的蒙片图像 与活片图像中的上区域与下区域分别设置不同图像阈值； 步骤4 :基于sift算法提取特征点：对步骤3中设置图像阈值后的每一对DSA脑血管 图像中的蒙片图像与活片图像均基于sift算法，按照步骤3中所设置的上下区域不同阈值 提取特征点，得到对应的图像特征点集； 步骤5 :对步骤4中每一对提取特征点的DSA脑血管图像中的蒙片图像与活片图像均 基于数字减影技术得到相应活片图像的特征点差值图像； 步骤6 :提取步骤5中所有所述特征点差值图像的图像特征点集，利用相邻图像关系对 图像特征点集进行精确提取； 步骤7 :对步骤6中提取的图像特征点集进行区域生长，得到相应的脑血管图像。
2. 根据权利要求1所述的DSA脑血管图像自动分割方法，其特征在于： 所述的基于数字减影技术得到相应活片图像的特征点差值图像是指①、提取蒙片图像 的图像特征点集VI、活片图像的图像特征点集V2后；②、遍历蒙片图像的图像特征点集VI 中的各个特征点，并在活片图像上找到相应特征点对应的位置坐标；③、以所述位置坐标为 圆心作圆后，遍历活片图像的图像特征点集V2,将处于所述圆的圆内且包含圆上特征点去 除，仅保留处于所述圆的圆外的特征点；④、将保留的特征点映射到活片图像对应位置上， 得到特征点集差值图像。
3. 根据权利要求1所述的DSA脑血管图像自动分割方法，其特征在于：所述利用相邻 图像关系对图像特征点集进行精确提取是指：①、提取所有映射到活片图像上的差值特征 点集；②、依次将相邻三幅活片图像做以下步骤：i、所述相邻三幅活片图像的差值特征点 集分别设定为Tl、T2、T3 ; ii、遍历T1中的各个特征点；iii、在T2中找到与ii中所述特征点 对应的位置坐标，并以其为圆心作圆，保留所述圆圆内的特征点，得到Tl、T2共同特征点集 T4 ;iv、对T2、T3重复ii、iii步骤得到T2、T3共同特征点集T5后，对共同特征点集T4、T5求 并集，得到综合特征点集T6。
【文档编号】G06T7/00GK104504708SQ201410830718
【公开日】2015年4月8日 申请日期:2014年12月26日 优先权日:2014年12月26日
【发明者】刘斌, 郝玲玲, 陈倩茹, 井晓彤, 朱琛, 江乾峰 申请人:大连理工大学