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

《基于.snake模型的血管树骨架三维重建技术》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、.........................基于Snake模型的血管树骨架三维重建技术曹治国，彭博基金项目：国家重点研究发展计划973项目（2003CB716105）作者联系方式：彭博Email:pbo820722@yahoo.com.cn通讯地址：华中科技大学图像识别与人工智能研究所，湖北武汉430074，桑农，张天序华中科技大学图像识别与人工智能研究所，湖北武汉430074摘要：从两幅不同角度的造影图像实现血管树的三维骨架重建，传统重建方法常常需要较多的人工干预，才能为每一点找到准确的对应点。文中研究了一种基于Snake模型的重建方法。它在采用多尺度Ga

2、bor滤波提取造影血管中轴线的基础上，选取血管树分叉点，优化几何变换矩阵，提高用以初始化snake的对应点的准确性；然后，针对血管的特性，文中采用GVF（GradientVectorFlow）流量场作为三维Snake的外部能量场，并给出最小化血管snake模型能量函数的表达式，使Snake保持自身平滑连续的同时在空间中发生形变以逼近真实轮廓。实验结果表明，该方法与传统的重建方法相比，不但减少了人工干预，而且有效地提高了重建精度。关键词：血管造影;三维重建;snake模型;GVF（GradientVectorFlow）1.引言血管造影图像（DSA）是临床诊治冠心病的

3、主要依据。血管的三维重建技术不仅能为医生提供形象、直观的三维血管形状图像，而且可以辅助测量血管的有关参数(如直径大小、血管长度和截面积等)，从而有助于冠心病的诊断和治疗。因此根据血管造影图像重建血管的三维骨架具有很好的临床意义和很高的应用价值。在传统的双平面造影图像重建方法中，大部分都是采用自底向上的重建方法，这种方法是基于图象的特征提取并通过找到每一对对应点进行重建[1][2]。这类方法最大的问题是，找到全部完全准确的对应点很困难，从而使得重建结果不精确。一些研究人员在寻找对应点方面做了大量的改进[3][4]，这些方法虽然提高了寻找的对应点的精确度，但同时也大大

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

5、建同时进行的，由于成像质量的影响和分割技术的局限性，影响了寻找的对应点之间的准确性，并最终影响重建精度。文献[8]则是采用Snake模型实现血管树骨架三维重建，不过，它与上述文献均未考虑对成像参数的校正，而未经校正的造影系统的参数（如造影角度，X射线源到接受平面的距离等）总是不可避免的存在一定的误差，从而造成在实际应用中所得结果也会产生一定误差。针对这些问题，本文研究了一种基于Snake模型的血管树骨架三维重建方法。它在采用多尺度Gabor滤波提取造影血管中轴线的基础[9]上，根据血管中轴的树状结构选取血管树分叉点，用以优化几何变换矩阵，提高用以初始化snake的

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

7、程；最后，对实验结果进行讨论并出结论。2造影图象几何关系标准的造影设备由两个独立的，可以自由旋转的X光系统组成[11]。两个旋转轴可以在水平和垂直方向上运动。两光轴交于光心O，X光源A,B到两图象平面的距离SIDA,SIDB已知。三维的坐标重建需要一个世界坐标系。我们假设两光轴互相垂直，世界坐标系的原点在光心，X,Y,Z轴如图1所示，A点发出的光线与X轴重合，B点发出的光线与Z轴重合。A,B各自的坐标系如图1，为平面A上的坐标系统，为平面B上的坐标系统。另外存在世界坐标系的同时，光源A,B也存在各自的坐标系。在图1中，A光源坐标系的x轴与世界坐标系的X轴重合，y,

8、z轴分别与