基于Snake 模型的血管树骨架三维重建技术

基于Snake模型的血管树骨架三维重建技术

曹治国，彭博*，桑农，张天序

华中科技大学图像识别与人工智能研究所，湖北武汉430074

摘要：从两幅不同角度的造影图像实现血管树的三维骨架重建，传统重建方法常常需要较多的人工干预，才能为每一点找到准确的对应点。文中研究了一种基于Snake模型的重建方法。它在采用多尺度Gabor滤波提取造影血管中轴线的基础上，选取血管树分叉点，优化几何变换矩阵，提高用以初始化snake的对应点的准确性；然后，针对血管的特性，文中采用GVF（Gradient Vector Flow）流量场作为三维Snake的外部能量场，并给出最小化血管snake模型能量函数的表达式，使Snake保持自身平滑连续的同时在空间中发生形变以逼近真实轮廓。实验结果表明，该方法与传统的重建方法相比，不但减少了人工干预，而且有效地提高了重建精度。

关键词：血管造影; 三维重建; snake模型; GVF（Gradient Vector Flow）

1.引言

血管造影图像（DSA）是临床诊治冠心病的主要依据。血管的三维重建技术不仅能为医生提供形象、直观的三维血管形状图像，而且可以辅助测量血管的有关参数(如直径大小、血管长度和截面积等)，从而有助于冠心病的诊断和治疗。因此根据血管造影图像重建血管的三维骨架具有很好的临床意义和很高的应用价值。

在传统的双平面造影图像重建方法中，大部分都是采用自底向上的重建方法，这种方法是基于图象的特征提取并通过找到每一对对应点进行重建[1][2]。这类方法最大的问题是，找到全部完全准确的对应点很困难，从而使得重建结果不精确。一些研究人员在寻找对应点方面做了大量的改进[3][4]，这些方法虽然提高了寻找的对应点的精确度，但同时也大大增加了计算的复杂性和计算量。

文献[5]提出了采用Snake模型[6]重建血管造影图像的自顶向下的方法，与传统的方法相比，它只需选用少量可以描述血管轮廓的点作为对应点即可，并且Snake的自动性较好，它通过同时结合两幅图像上的特征，使Snake在空间中不断发生形变，从而逼近真实轮廓，获得较高的重建精度。Radeva和Toledo等在此基础上重新定义了Snake的外部能量场，并给出了基于三维Snake模型重建任意角度的造影图象的一般方法[7]。上述这些工作都只给出单

*基金项目：国家重点研究发展计划973项目（2003CB716105）

作者联系方式：彭博Email: *******************.cn

通讯地址：华中科技大学图像识别与人工智能研究所，湖北武汉430074

根血管的重建方法，并且他们是将血管分割与重建同时进行的，由于成像质量的影响和分割技术的局限性，影响了寻找的对应点之间的准确性，并最终影响重建精度。

文献[8]则是采用Snake 模型实现血管树骨架三维重建，不过，它与上述文献均未考虑对成

像参数的校正，而未经校正的造影系统的参数（如造影角度，X 射线源到接受平面的距离

等）总是不可避免的存在一定的误差，从而造成在实际应用中所得结果也会产生一定误

差。

针对这些问题，本文研究了一种基于Snake 模型的血管树骨架三维重建方法。它在采用多尺度Gabor 滤波提取造影血管中轴线的基础[9]上，根据血管中轴的树状结构选取血管树分叉点，用以优化几何变换矩阵，提高用以初始化snake 的对应点的准确性；然后，本文针对血管的特性，提出采用GVF （Gradient Vector Flow ）流量场[10]作为三维Snake 模型的外部能量场，并给出血管的非闭合轮廓能量函数最小化的最佳表达式，使Snake 保持自身平滑连续的同时在空间中发生形变以逼近真实轮廓。由于采用多尺度Gabor 滤波器提取中轴的方法我们已在文献[9]中进行了详细的描述，因此，本文不再赘述，其它内容则按以下结构组织：第二节简单给出X 射线造影系统的几何描述，和基于几何关系重建的一般方法；第三节介绍初始化Snake ，并给出利用对应点对优化几何变换矩阵方法；第四节分析利用Snake 重建三维血管树的过程；最后，对实验结果进行讨论并出结论。

2造影图象几何关系

标准的造影设备由两个独立的，可以自由旋转的X 光系统组成[11]。两个旋转轴可以在水平和垂直方向上运动。两光轴交于光心O ，X 光源A,B 到两图象平面的距离SIDA,SIDB 已知。三维的坐标重建需要一个世界坐标系。我们假设两光轴互相垂直，世界坐标系的原点在光心，X,Y,Z 轴如图1所示，A 点发出的光线与X 轴重合，B 点发出的光线与Z 轴重合。A,B 各自的坐标系如图1，A A u v 为平面A 上的坐标系统，B B u v 为平面B 上的坐标系统。另外存在世界坐标系的同时，光源A,B 也存在各自的坐标系。在图1中，A 光源坐标系的x 轴与世界坐标系的X 轴重合，y,z 轴分别与Y,Z 轴平行。假设P 点在光源A 坐标系下的坐标为(,,)A A A x y z ，在光源B 坐标系下的坐标为(,,)B B B x y z 。则他们之间满足的关系可以写为：

B A B A B A x x y R y t z z ⎡⎤⎡⎤⎧

⎫

⎢⎥⎪⎪

⎢⎥=⋅-⎨

⎬⎢⎥⎢⎥⎪⎪⎢⎥⎢⎥⎩⎭

⎣⎦

⎣⎦

（1）

其中(,)R t 为几何变换矩阵，与造影角度和X 射线源到接受平面的距离等有关[11]。R 是旋转矩阵，t 是平移向量。

在此三维空间中，一个点P 在一幅图象上的投影可以定义为光源发出的光束经过该点与图象平面的交点A P ，B P ，如图即为直线A AP ，B BP 的交点。当两幅图象上的对应点已知，两条光线将在空间中交于一点。在这种情况下，如果成像系统的几何关系已知，并且一对对应点已找到，空间中的三维点是可以重建出来的。

然而实际的情况是，这两束光线可能并不相交，如图2，这主要是因为系统的误差和图像间几何关系的不精确度造成的。为了解决这个问题，一般是将到两条直线距离最近的点作为真实的点的存在位置[2]，如图2。

B

B y B

A

u

B

B

ImageA

ImageB

图1 双平面造影系统结构（两光源夹角为90度） 图2 由A P ，B P 重建三维点P

3 基于几何变换矩阵优化的Snake 模型初始化

对于传统的重建方法而言，找到准确的对应点是十分重要的工作。其准确与否直接影响到重建的结果。文献[3]和[4]在寻找对应点方面做了大量的改进。而一些研究人员[5][7][8]则把三维的Snake 模型应用于血管的重建上，避免了为血管上的每一点寻找对应点的问题，只需要寻找若干对对应点（如起止和分叉点）即可描述出血管的大致轮廓。再结合Snake 投影到两幅图像上的信息使Snake 不断发生形变，直到它的投影与两幅图像完全重合。