基于OPTA 细化算法与有关脑动脉瘤检测研究吴建, 张广明, 谢杰, 崔志明2009 年8 月21 日~ 23 日中国黄山20.docVIP

Research on Cerebral Aneurysm Detection Based on OPTA Algorithm Jian Wu, Guangming Zhang, Jie Xia, and Zhiming Cui Proceedings of the 2009 International Symposium on Information Processing Huangshan, P. R. China, August 21-23, 2009, pp. 037-040 基于OPTA细化算法的有关脑动脉瘤检测的研究 吴建，张广明，谢杰，崔志明 2009年8月21日~ 23日中国黄山2009年信息处理国际研讨会 论文集 037页~040页 摘 要 在脑动脉瘤识别系统中，这是能够准确、快速地在图像中定位出脑动脉瘤位置的关键步骤。经过了脑动脉瘤的形态特征的详细分析后，本文提出了一种新的脑动脉瘤检测方法，这种检测方法是基于改进了的细化算法。在这种新的检测方法中，首先利用OPTA细化算法提取出血管的骨架，然后根据血管的骨架信息来检测脑动脉瘤。经过大量实验之后发现，脑动脉瘤识别过程在以这种新的检测方法为前提的情况下，可以很好的检测到脑动脉瘤的位置。 关键词： OPTA， 细化算法， 脑动脉瘤检测， 模板匹配 Ⅰ.前言 脑血管疾病，尤其是脑动脉瘤，是导致成年人生病最后死亡的关键因素之一，它严重威胁着人们得生命安全。随着计算机技术的不断发展和成熟，信息技术和医学成像技术结合而产生的CAD(应用计算机辅助诊断)系统在脑血管疾病的检测与治疗中起到越来越重要的作用，它已经成为了医学成像上的一个研究重点。 脑血管瘤一般位于血管的交叉位置，尤其是在脑动脉周围。原因是血液的流动对血管交叉位置的影响很大。脑血管的影响类似于河流形成的网络，会出现许多分支的动脉血管。一般来说，血管是对称的，它的两侧的血管壁是相互平行的。而脑动脉瘤是由于血管壁损坏而导致的突出的部分。脑动脉瘤的示意图如图1所示。 图1 脑动脉瘤的示意图 如图1，图中用方格标志的地方就是脑动脉瘤。正常的血管出已经形成了突出的部分，大致平行的血管壁被破坏了。被破坏的位置明显是在交叉处。所以我们可以确定脑动脉瘤的位置应该在血管骨架结构的交叉处。 在基于DSA(数字减法血管造影术)的脑动脉瘤CAD系统中，它是特征提取及识别的前提和重要步骤，应用在检测脑动脉瘤位置的DSA中。本文分析了脑动脉瘤的形态特征，包括通过细化算法得到的血管骨架的拓扑结构，然后对骨架树进行深度优先遍历，最后定位出脑动脉瘤的位置。 Ⅱ.OPTA算法 OPTA(一次通过细化算法)是一种典型的基于模板的图像细化算法，其核心是通过应用消去和保留模板实现细化过程。OPTA是一个迭代的过程。如果当前的点满足消去模板同时不满足保留模板，那么这个点就被消去，否则，保留此点。不断遍历原始图片，直到再没有点满足上述要求。 对OPTA算法的主要改进的地方就在于根据原算法的效果和速度增加了消去和保留模板。其中最典型的参考文件是[1~3]。文献[1]的作者发现原始的OPTA算法细化的不完整，有一定的缺陷，细化出的图像有毛刺且不够平滑。改进的算法中已经将消去和保留模板进行了改善。图2和图3是改进了的消去和保留模板。其中有阴影的方格是当前点。 图2 消去模板[1] 图3 保留模板[1] 从文献[1]可以看出，在原始算法中，保留模板太松弛了，分析图像可以得知，会容易产生右斜线，且使迭代次数会增加，速度会变慢。因此，在文献[1]的基础上突出了新的保留模板，并增至9个。从图4可以看出，保留模板的条件变得更加严格，从而解决了右斜线的问题，速度也有明显的改善。 图4 保留模板[2] 通过深入的改进，由袁梅和其他作者发现的右斜线的影响已经被文献[2]中的算法很好的解决了，且速度有了很大的提升。由于增强了的保留条件，使得图像出现了断点。从而导致了拓扑结构和原始图像被破坏了，细化算法也被减弱了。在以上分析的基础上，袁梅和其他的作者提出了组合模板的概念，利用文献[1]的消去模板，放弃其中的保留模板，并且结合了对当前点的判断，判断当前点是否符合消去模板，符合的话保留它相邻点的位置。图5是当前点以及他的邻接点的示意图。 图5 点p的邻接点 基于组合模板的细化算法不需要匹配保留模板的操作。每个消去模板对应着一种情况，并且每个消去模板利用“&”和“||”操作来完成判断当前点是否要被保留。下面的8个条件式对应图2中a~h的8个消去模板。如果该点满足消去模板以及消去条件，则消去此点。否则暴力此点。 (q4=0 & q5=0 & q6=0 & q7=0) || (q5=1 & q6=1 & q12=0 & q13=0 & q14=0) (1) (q2=1 & q4=0 &