[发明专利]微创血管介入手术中导丝的实时运动仿真方法

说明书

技术领域

本发明涉及一种虚拟手术的模拟方法，尤其是涉及一种微创血管介入手术中导丝的实时运动仿真方法。

背景技术

随着科学技术的高速发展，虚拟手术的应用为医疗的现代化提供了越来越多的帮助。虚拟的微创血管介入手术系统能够帮助实习医生快速提高手术技能，让更多的心脑血管疾病患者能够进行微创介入治疗。导丝的实时运动仿真作为微创血管介入手术系统中不可缺少的一部分，成为许多学者研究的课题。在导丝实时运动仿真的研究过程中，早期的一部分学者采用几何模型对导丝进行建模，这种方法完全忽略了导丝的物理性质，然而如果完全忽略了导丝的物理性质，与其他模型之间的交互将难以处理。也有一部分学者采用质点弹簧这样简单的物理模型来对导丝进行建模，这种方法虽然计算效率很高，但是由于模型过于简单，根本无法真实地模拟导丝在血管中的运动。本文实现的微创血管介入手术中导丝的实时运动仿真对导丝的模拟有很高的要求，不仅要求导丝要有高的真实性，还对导丝与血管之间交互的模拟要有可靠的物理模型，并且由于要应用于虚拟手术，这就要求模型系统达到高的更新速度，并且能达到实时。但是微创血管介入手术中的导丝的实时运动仿真对血管的模拟要求较低。一个完整的微创血管介入手术中的导丝的实时运动仿真包括了血管模拟、导丝模拟、碰撞检测及响应、导丝的不可伸缩和力反馈渲染。

对于血管的模拟，由于血管在仿真的过程中几乎不产生形变，只需要采用三角形表面网格这样简单的表面模型对血管建模就足够了。不需要利用四面体网格这样复杂的立体模型来模拟血管。这样有效地降低了计算成本和存储空间。

导丝的模拟，是微创血管介入手术中导丝的实时运动仿真中的重点部分，算法设计的好坏关系到整个仿真系统的真实性、准确性、高效性。目前，大部分的方法是多刚体运动模型，弹簧模型，有限元模型，B-Spline等理论，这其中弹簧模型和B-Spline被常用于导丝的模拟，由于其效率高能达到实时。然而弹簧模型过于简单，忽略了导丝很多重要的物理性质而导致真实感很差。B-Spline的缺点是缺乏物理支持，只是从几何的角度来仿真导丝在血管中的形态，无法知道导丝的受力情况。因此无法为后面的力反馈设备提供真实的力反馈。有限元模型能够提供很强的真实性。但是有限元方法有一个致命的弱点：实时性差。不能满足在仿真过程中与用户进行实时交互的要求。为了同时达到真实性和实时性的目标，本文采用了基尔霍夫弹性杆对导丝进行物理建模。

主流的碰撞检测方法之一是，采用层次结构(hierarchy)为基础核心，例如Octree，k-d树，BSP树等，利用包围体的hierarchy结构也被广泛地应用于各种碰撞检测算法(同时这种方法也被应用于其他领域，如：光线追踪)，包围体的选择，学者们也提出了许多方法，包括，球包围体，有轴向的长方体(AABBs)另外一些碰撞检测方法，包括采用space-time bounds以及采用四维几何体来包围物体，近些年，学者们也提出基于Voronoi图的碰撞检测方法。

导丝与血管之间的碰撞响应，目前的两种主流算法是惩罚法和约束法。惩罚法的真实性和稳定性差，然而约束法需要进行多次迭代才能确保导丝不刺穿血管壁，计算成本很高。结合惩罚法与约束法的优点，本文采用一种非迭代的约束方法来计算导丝与血管之间的接触力并进行碰撞响应。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种稳定可靠、实现方便、真实高效的微创血管介入手术中导丝的实时运动仿真方法。

本发明的目的可以通过以下技术方案来实现：

一种微创血管介入手术中导丝的实时运动仿真方法，其特征在于，该方法包括以下步骤：

1)由于血管在仿真的过程中不产生形变，采用三角形表面网格对血管进行建模；

2)采用基尔霍夫弹性杆模型，对导丝进行建模，并利用拉格朗日方程计算导丝在仿真过程产生的形变力；

3)进行导丝的渲染，相邻的两导丝质点之间绘制一个圆柱体，同时在相邻的两圆柱体之间绘制一个圆球进行连接；

4)采用K-DOPS树进行导丝与血管之间碰撞检测；

5)采用非迭代的约束方法计算导丝与血管之间的接触力并进行碰撞响应；

6)使用Verlet积分公式来迭代更新导丝的位置与速度；

7)利用拉格朗日乘数和快速投影方法来实现导丝的不可伸缩特性；

8)使用力反馈设备Phantom Omni来渲染导丝给用户的作用力。

所述的采用基尔霍夫弹性杆模型，对导丝进行建模，并利用拉格朗日方程计算导丝在仿真过程产生的形变力的具体实现过程为：