[发明专利]一种基于径向基点插值法与质点弹簧法的软组织建模方法

说明书

本发明公开了一种基于径向基点插值法与质点弹簧法的软组织建模方法，其包括以下步骤：利用三维建模技术为采集的医学图像信息构建软组织模型，所述软组织模型按功能区划分为手术区和非手术区，并基于支持域的边缘提取策略建立过渡区；预处理模型参数以及建模方法中的参数；手术区模拟软组织模型的形变行为，力经过过渡区传递至非手术区；计算形变的单步迭代结束，进入下一轮循环。本发明将软组织模型分为手术区与非手术区，手术区采用高精度、便于施加边界条件的RPIM，引入了隐式积分法Newmark法对位移的变化进行动态分析并保证系统的稳定性；在非手术区使用实现简单，计算效率高的MSM。

技术领域

本发明属于软组织形变模型建模技术领域，具体涉及一种基于径向基点插值法与质点弹簧法的软组织建模方法，用于解决虚拟手术实验中软组织建模问题。

背景技术

虚拟手术提供了较理想的手术仿真效果和教学环境，替代传统的手术技能培训，极大地减少训练成本，有助于术前计划制定和手术操作练习，同时对传统的手术教学模式提出了新改变。计算机技术的迅速发展及现代医学的不断进步，促进了这一多学科交叉的先进领域在医学培训上的强大作用。因此虚拟手术的相关研究，例如软组织建模方法研究具有重要的现实意义。

软组织模型的建立是虚拟手术仿真系统中的基础，模拟形变和切割仿真是虚拟手术仿真系统的关键技术。在研究物体形变的早期，由于基于非物理的几何算法实现简单、易于操作等优点，使其占据相关研究的主流地位。几何模型从源数据CT断层图像获取表面模型的基本单元信息，没有大量复杂的物理计算，但是缺乏材料信息，难以对拓扑结构复杂的模型进行形变模拟，影响逼真度。

随后基于物理的建模方法被不断提出且日渐成熟，从最初的仅仅包含速度、角速度的运动学方程发展为包含质量力，转动惯量等物理特性的运动学方程。与非物理模型相比，物理模型能够更真实地反映虚拟手术过程中的形变过程。目前，有限元方法(FiniteElement Method，FEM)、质点弹簧模型(Mass-Spring Method，MSM)、无网格(Meshless)等物理模型成为了软组织建模方法研究的主流。

有限元法是一种数值解法，寻求难以获得解析结果的问题的场变量分布的近似解。在有限元法中，一个复杂形状的连续体被划分为不同单元，既有限单元。各个单元通过一个被称为网格的拓扑图连接在一起。各个有限单元的单元体形状各异，单元体之间的连接方式不尽相同，有限元法可以应用于复杂模型的求解域进行模块化处理，基于FEM的软组织模型具有较高的精确性。因此，有限元法在软组织建模中得到了广泛的应用。然而，有限元的局限性十分明显：形成FEM网络时的计算成本高；应力精度低；自适应分析信息困难，在求解特大形变的问题时，可能会产生扭曲网格，导致网格需要重构，严重影响解的精度。其根源是有限元法系统方程阶段利用了单元或网格信息。

质点弹簧模型自提出之初就收到研究者们的青睐。质点弹簧模型是一种基于胡克定律的弹性变形模型。在质点弹簧模型中，模型被离散为若干个被约束连接的质点，这种具有一定物理特性的约束称为弹簧。当单个质点或者多个质点受到外力时，质点的位移变化间相互影响并传递从而达到形变的目的，实现简单，计算量小。由于质点弹簧的动力学方程没有集成材料的物理参数，模型的结构方程中所需参数由经验得到。另外除了质点间的弹簧张量，质点所受的力矢量还与网格产生变化时对质点在不同方向上的力和位移变化补偿相关。例如模型中出现的弹簧扭转，约束变化及体积变化。

无网格法是建立在整个问题域的系统代数方程，不需利用预定义的网格信息进行离散的方法。无网格法利用一组散布的节点被称为场节点，他们并不构成网格，及不需要任何事先定义的节点连接信息用于构造变量未知函数的插值或者近似表达式。无网格方法基于点的近似，可以彻底或者部分消除网格，不需要网格的初始划分重构，不仅可以保证计算的精度，而且相较于有限元法，可以减小计算的难度。