CTA图像冠状动脉开口层与开口点的全自动定位方法

罗 洋，齐 林，徐礼胜，3

(1. 东北大学 医学与生物信息工程学院，辽宁 沈阳 110169；2. 鞍山师范学院 高等职业技术学院，辽宁 鞍山 114056；3. 东北大学 教育部医学影像计算重点实验室, 辽宁 沈阳 110169)

2016年，全世界1 790万人死于心血管疾病，占全球死亡人数的31%[1].早期诊断和治疗是降低冠心病发病率和死亡率的关键.计算机断层冠状动脉造影(computed tomography coronary angiography, CTCA)已越来越多地用于冠心病的无创诊断.冠状动脉分割是CTCA图像定量分析的重要前提条件.然而，人工分割繁琐费时，诊断精度完全依赖于医生的先验知识、图像上下文信息和医生的主观判断，观察者内部和观察者之间的变异性相对较高.因此，需要一种冠状动脉自动分割方法来加速和简化诊断、治疗和手术的过程.最近十几年产生了很多冠状动脉分割的方法，其中包括阈值方法[2]、区域生长方法[3]、Hessian矩阵方法[4]、匹配滤波方法[5]、图割方法[6]、动态轮廓方法[7]、水平集方法[8-9]、机器学习方法[10]和其他方法[11-14]等.然而这些方法对于冠状动脉开口层和开口点的定位主要采取人工标注的方式，这种方法需要花费放射科医生大量时间和精力.最近十几年提出的自动检测冠脉开口的方法包括[15-18]：Hennemuth等[15]首先对主动脉进行明确的分割，然后把开口作为主动脉的管状结构附着物，检测率较低；Wang等[16]采用基于Hough变换的二维圆检测算法对升主动脉进行检测和跟踪，并选取主动脉根部两侧最大的连通成分作为左、右冠状动脉，检测率低；Zheng等[17]针对非搭桥冠脉定义一个边界盒(立方体)，冠状动脉口位于立方体两个对面的中心.针对搭桥冠脉利用先验知识训练一个分类器来识别冠脉开口；Gong等[18]首先提取升主动脉，然后利用管状物增强和形态学操作求最大响应值来提取冠脉开口.本文提出一种CTA图像冠状动脉开口层与开口点的全自动定位方法.这是一种准确的方法，不需要任何标注、训练.

1 数据与方法

1.1 数据

本文的数据集采用从北部战区总医院随机选取的CTCA图像.数据集由57个对象组成，其中包含30个正常人和27个病人，每个病例的图像数量为290到428不等.CTCA图像的大小为512×512(像素).每个病例的左右冠状动脉开口的位置和开口点由北部战区总医院的放射科医生手工标注.

1.2 算法原理

本文介绍了在多条件区域增长基础上通过圆度的变化确定左冠脉开口层，并在指定范围内确定左冠脉开口点，然后计算升主动脉中心和右冠脉开口范围，通过求出距离曲线，从而确定右冠脉开口层和开口点的方法.算法框图见图1，图1a是求左冠脉开口层和开口点的算法框图，图1b是求右冠脉开口层和开口点的算法框图，图1a要先于图1b实现.

心脏冠状动脉树有两个分支:左冠状动脉(左冠脉)和右冠状动脉(右冠脉).其中每个分支只有一个根节点，位于整个冠状动脉偏上的位置、升主动脉的开口处.左冠脉(left coronary artery，LCA)的开口处位于升主动脉左后方，开口呈横位的椭圆形，位置略高于右冠脉开口.右冠脉(right coronary artery，RCA)的开口处位于右冠状动脉的外侧中上部，位于升主动脉根部的右前方[19]，具体结构如图2所示.

通常在患者的CTCA图像序列开始，升主动脉和降主动脉与周围区域灰度对比明显，且升主动脉和降主动脉都近似圆形，同时升主动脉面积明显大于降主动脉面积，如图3a所示.但有的患者CTCA图像序列开始显示的是主动脉弓，如图3b所示.对于这种情况从它的总数的第10%张图像开始能看到正常的升、降主动脉，如图3c所示.

本文算法如下:

1) 每名患者从第10%张图像开始，利用霍夫变换检测圆的方法找到升主动脉.

2) 把当前层的升主动脉中心(圆心)作为下一层的种子点，在种子点为中心周围1111区域内求平均灰度Vs，对当前升主动脉面积Ac、左边界Lc,上一层升主动脉面积Ap、左边界Lp，进行区域增长，执行如下代码:

fori=0∶10∶160

ifAc/Ap>1.2‖(Lc-Lp)<-5

%下层升主动脉面积/上层升主动脉面积小于等于1.2,下层升主动脉左边界相对于上层升主动脉左边界向左移动不超过5个像素(2.07 mm)

Regiongrow(Vs-(160-i)，Vs+(160-i))

%缩小阈值范围，重新区域增长

else

break；%找到升主动脉

end

end

3) 求升主动脉的圆度roundness:

(1)

其中:π是圆周率;area代表面积;perimeter代表周长;roundness介于0和1之间，值越大越圆.

4) 如果圆度≥0.82，重复步骤2),3),4);如果圆度<0.82，则找到左冠脉开口层.过中心点画斜率为1的直线l1，过圆心点画垂线l2，轮廓线上才可能出现左冠脉开口，如图4a所示.利用解剖结构和形态学运算，求出开口处左冠脉，如图4b所示.距离中心点最近的点(白色)即为左冠脉开口点，如图4c所示.

5) 右冠脉开口接近升主动脉根部，而在根部升主动脉可能与其他组织有连接，连接后形状无规则，如图5a和图5b所示.但是本文发现右冠脉开口会出现在升主动脉的左上方，如图5b所示.冠状动脉图像序列中，升主动脉横截面随序列递增向下向右缓慢移动，本文提出计算当前层升主动脉中心点坐标的公式：

(2)

其中:(Xl,Yl,Zl)表示左冠脉开口层的升主动脉中心点坐标;(Xc,Yc,Zc)代表当前层升主动脉的中心点坐标.

6) 过初始中心点(黑色“+”号)斜率为1和-1的两条直线l1和l2之间的区域，如图6a中l2与l1间加粗曲线轮廓所示.重新计算中心点(白色“+”号)，新的中心点在轮廓两个端点连线(白色)的垂直平分线(白色)上，到两个端点的距离等于首层升主动脉的半径，如图6a所示.求从新中心点到黑色轮廓上每一点的距离曲线，如图6b所示.如果某个波峰左右各10个值之内最小值都要比波峰值小3以上，表明波峰比较明显，则右冠脉开口层出现了.通过形态学运算，如图6c所示，距离新中心点最近的点(白色)即为右冠脉开口点，如图6d所示.

2 结 果

对57个病例做了定位测试.垂直于横截面的方向从上向下，左右冠脉的开口层通常包括2层以上.本文认为定位到刚发现开口的前两层都是准确的.如表1所示，在左冠脉开口层定位的结果中，有56个病例的定位准确，准确率达到98.25%.而对右冠脉开口层的定位中，有55个病例的定位准确，准确率达到96.49%.

表1 定位左右冠脉开口的准确率

Hennemuth等[15]在61个数据集上测试的冠脉开口的检测率只有57%.Wang等[16]提出的方法，在16个数据集上的检测率仅为80%左右.Zheng等[17]提出方法的几种检测率都低于96%.Gong等[18]的平均检测率为96.88%.与这几种方法相比，本文方法的平均检测率为97.37%，高于其他几种方法.

如表2所示，在左右冠脉开口点的定位中，左冠脉开口点的位置(定位准确的56个病例)与左冠脉金标准开口点的位置平均距离是(1.22±0.37) mm，右冠脉开口点的位置(定位准确的55个病例)与金标准开口点的位置平均距离是(1.76±0.49) mm.

表2 定位左右冠脉开口点的平均距离

3 讨 论

在数据集中，本文统计了左冠脉开口层号占总层数的比例关系，如图7所示.左冠脉开口层号占总层数百分比最小都在15%以上.选择每个病例开口层号为总层数的10%，这样既保证没错过左冠脉开口层，又减少了需要处理图像层数，同时保证开口层显示的是近似圆形的升主动脉和降主动脉.进行区域增长时，选取以中心点(种子点)为中心的11×11区域求平均灰度，这样既保证不会超出升主动脉区域，又保证区域增长灰度区间的合理性.

根据对57个病例图像的统计发现，在左冠脉开口出现前，每张图像的升主动脉都近似圆形，圆度都大于等于0.82，如图8a所示.左冠脉开口层的升主动脉圆度都小于0.82，如图8b所示.

左冠脉开口出现前，相邻图像升主动脉左边界变化很小，当前图像升主动脉左边界相比上一张图像升主动脉左边界向左移动小于5个像素(2.07 mm)，设置这么一个限制是防止区域增长时把升主动脉左侧不相干的组织增长进来.根据对数据集的统计发现，左冠脉开口层前所有相邻图像升主动脉面积之比(下层面积比上层面积)的均值非常接近1，如图9a所示.左冠脉开口层与前一层图像升主动脉面积之比的均值比1稍大，最小值也大于1，如图9b所示.图9说明了升主动脉在左冠脉开口层(含开口层)前面积、形状都变化很小，比较稳定.比值最大接近1.14，如果在区域增长时阈值不合适，通常会把周围某些组织增长进去，面积会有一定增加，那么当前图像与上层图像的升主动脉面积之比通常会明显大于1.14，所以在判断区域增长的阈值是否恰当时，把是否大于1.2作为一个条件是合理的.由于升主动脉每层的面积变化很小，所以在求右冠脉开口层和开口点的过程中，重新计算中心点(白色“+“号)时，为了新的中心点尽可能在升主动脉的中心，保证重新计算的中心点与加粗轮廓两个端点的距离等于首张图像中升主动脉的半径，如图6a所示.

4 结 论

1) 定位左右冠脉开口层的准确率分别是98.25%和96.49%，平均定位时间是41.25 s.定位左右冠脉开口点的平均距离误差分别是(1.22±0.37) mm和(1.76±0.49) mm.就准确率来说，本文方法较高，优于其他方法.

2) 对本文提出的方法所涉及到的参数做了讨论分析，说明了该方法的可行性.

3) 本文提出的方法对于左右冠脉开口位置异常的情况无法做出准确定位，如果定位左冠脉开口位置错误，则可能影响到定位右冠脉开口位置和开口点的准确率，本文方法还需要进一步改进.