CN116561935B - 一种基于体弹簧和多尺度体积不变的软组织形变模拟方法 - Google Patents
一种基于体弹簧和多尺度体积不变的软组织形变模拟方法 Download PDFInfo
- Publication number
- CN116561935B CN116561935B CN202310841839.7A CN202310841839A CN116561935B CN 116561935 B CN116561935 B CN 116561935B CN 202310841839 A CN202310841839 A CN 202310841839A CN 116561935 B CN116561935 B CN 116561935B
- Authority
- CN
- China
- Prior art keywords
- soft tissue
- current
- deformation
- volume
- tissue model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 210000004872 soft tissue Anatomy 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000004088 simulation Methods 0.000 title claims abstract description 22
- 239000002245 particle Substances 0.000 claims abstract description 84
- 238000013334 tissue model Methods 0.000 claims abstract description 61
- 230000008569 process Effects 0.000 claims abstract description 30
- 238000013016 damping Methods 0.000 claims description 7
- 238000013459 approach Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000003993 interaction Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000012549 training Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于体弹簧和多尺度体积不变的软组织形变模拟方法,属于软组织模型技术领域,包括:S1、构建软组织模型;所述软组织模型由质子弹簧模型和虚拟体弹簧形成的若干四面体单元构成,每个所述四面体单元中包含3个虚拟体弹簧,用于获取软组织模型表面节点信息及内部体积信息;S2、计算软组织模型多尺度体积不变的体积;S3、基于软组织模型多尺度体积不变的体积,模拟软组织模型的形变过程,确定模型质点在形变过程中的位置坐标,实现软组织形变模拟。本发明方法为虚拟软组织模型快速回弹,提高计算速率,以及更好与力反馈设备交互,有效解决虚拟仿真模型真实性与实时性不高的问题。
Description
技术领域
本发明属于软组织模型技术领域,具体涉及一种基于体弹簧和多尺度体积不变的软组织形变模拟方法。
背景技术
随着虚拟现实技术的发展已经拓展到了医学领域,其中虚拟手术仿真技术就是其典型的应用之一,对于外科医生的培训,练习解决真实手术可能会面临的各种复杂情况起着至关重要的作用,针对现有的虚拟手术中存在的三维重建缺乏准确性,软组织力学参数和形变模型缺乏准确性,而力反馈设备的输入与软组织参数息息相关,软组织的形变效果和参数反馈直接决定了虚拟手术的真实性。虚拟手术仿真技术是虚拟现实与医学相融合的产物。
虚拟手术***由医学影像出发,重构出逼真的虚拟人体组织,搭建虚拟的医学环境,利用交互设备实现操作者与虚拟组织的交互。医生可以在虚拟的三维医学环境中对虚拟软组织实现触碰、切割、缝合等手术操作,实现逼真的外科手术。传统的手术训练主要以尸体标本解剖和临床手术观摩的方式进行,存在成本高、内容随机、训练周期长、效率低等诸多弊端。而虚拟手术***无需高成本的耗材,仅需一套计算机设备即可进行多次身临其境的、规范化的训练,可以作为常规教学以外的有力补充。此外,虚拟手术可以结合新兴的专家***进行远程指导,有助于提高训练质量与效率。
发明内容
针对现有技术中的上述不足,本发明提供的基于体弹簧和多尺度体积不变的软组织形变模拟方法解决了现有的软组织形变模拟过程中,虚拟仿真模型真实性与实时性不高的问题。
为了达到上述发明目的,本发明采用的技术方案为:一种基于体弹簧和多尺度体积不变的软组织形变模拟方法,包括以下步骤:
S1、构建软组织模型;
所述软组织模型由质子弹簧模型和虚拟体弹簧形成的若干四面体单元构成,每个所述四面体单元中包含3个虚拟体弹簧,用于获取软组织模型表面节点信息及内部体积信息;
S2、计算软组织模型多尺度体积不变的体积;
S3、基于软组织模型多尺度体积不变的体积,模拟软组织模型的形变过程,确定模型质点在形变过程中的位置坐标,实现软组织形变模拟。
进一步地,所述步骤S2具体为:
S21、使用AABB包围盒获取软组织模型的三维边界信息;
其中,三维边界信息包括X轴最大值、X轴最小值/>、Y轴最大值/>、Y轴最小值/>、Z轴最大值/>、Z轴最小值/>以及X轴、Y轴和Z轴的边界差/>,/>和/>;
S22、通过二分法对软组织模型进行网格划分,获得细分的立方体;
其中,进行网格划分时,每次取当前立方体三维边界的X轴长度、Y轴长度和Z轴长度的二分之一进行网格划分;
S23、判断当前立方体的八个顶点是否均在原立方体对应的凸多边形内;
其中,原立方体为当前立方体对应的网格划分前的立方体;
若是,则进入步骤S24;
若否,则进入步骤S25;
S24、设置当前软组织模型体积V=Vn-1+Vn,进入步骤S27;
其中,Vn为当前立方体的体积,Vn-1为前一次网格划分时原立方体的体积,下标n为当前网格划分次数的计数值,软组织模型体积V的初始值为0,n为大于等于1的自然数;
S25、判断当前立方体边界差是否小于等于软组织模型边界差的四十分之一;
若是,则进入步骤S26;
若否,则返回步骤S22;
S26、设置当前立方体的体积为,当前软组织模型体积V= Vn-1+/>;
S27、重复步骤S22~S26,直到网格划分次数达到预设值或划分后的所有立方体精度达到预设值,将当前软组织模型体积V作为软组织模型的多尺度体积不变的体积。
进一步地,所述步骤S3具体为:
S31、遍历软组织模型的所有质点,并计算当前质点的质量;
S32、遍历当前质点的所有相邻质点,计算当前质点与其相邻质点间在单步形变后的距离,并计算该距离下当前质点所受的弹簧力;
S33、根据当前质点所受的阻尼力和弹簧力,计算当前质点所受周围相邻质点的合力,并根据当前质点的运动方程计算其在单步形变后的位置;
S34、在有限形变次数的形变过程中内,重复步骤S32~S33,在软组织模型多尺度体积不变的基础上,计算当前质点在最近两次单步形变的位置之差,当位置之差逼近零时,则当前质点形变结束,进而确定当前质点形变后的位置坐标。
进一步地,所述步骤S31中,当前质点的质量为:
式中,为曲面的面密度,/>为当前质点i的一环邻域内的第j个三角面片的面积,j为质点i的一环邻域内的三角面片计数值,/>为三角片面的数量。
进一步地,所述步骤S32中,质点所受周围相邻质点的合力为:
式中,为与质点i相邻的质点数,/>为质点i受到相邻质点k的弹簧力,为质点i受到相邻质点k的阻尼力,/>为质点i受到的外力。
进一步地,所述步骤S33中,当前质点的运动方程为:
式中,为当前质点i的运动位移,/>为当前质点i的运动速度,/>为当前质点i所受周围相邻质点的合力,/>为当前质点i的质量。
进一步地,所述步骤S34中,在有限形变次数的形变过程中,当当前质点的实际形变次数超过设定的最大形变次数时,则当前质点形变结束,进而确定当前质点形变后的位置坐标。
本发明的有益效果为:
1)本发明方法为虚拟软组织模型快速回弹,提高计算速率,以及更好与力反馈设备交互,有效解决虚拟仿真模型真实性与实时性不高的问题。
2) 本发明将多尺度弹簧力学模型与质子弹簧模型进行比较,实验评价指标结果表明,所提出的方法在各向异性,模型仿真度,刷新率相对于原算法更加优异;
3) 本发明提高了软组织形变模拟过程的算法运算速度,降低了计算复杂度;
4)基于本发明方法使得软组织形变模型更贴近真实,更利于医生获取更加真实的力反馈效果;
5)本发明构建的软组织模型更贴近真实。
附图说明
图1为本发明提供的基于体弹簧和多尺度体积不变的软组织形变模拟方法流程图。
图2为本发明提供的加入了体弹簧的质子弹簧结构示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
本发明实施例提供了一种基于体弹簧和多尺度体积不变的软组织形变模拟方法,如图1所示,包括以下步骤:
S1、构建软组织模型;
所述软组织模型由质子弹簧模型和虚拟体弹簧形成的若干四面体单元构成,每个所述四面体单元中包含3个虚拟体弹簧,用于获取软组织模型表面节点信息及内部体积信息;
S2、计算软组织模型多尺度体积不变的体积;
S3、基于软组织模型多尺度体积不变的体积,模拟软组织模型的形变过程,确定模型质点在形变过程中的位置坐标,实现软组织形变模拟。
本发明实施例提供的上述形变模拟过程中,由质子弹簧模型(Mass-SpringModel, MSM)可以模拟真实的形变,因此使用加入体弹簧以及多尺度体积不变方法进行改进,通过加入体弹簧的形式优化软组织形变过程中的力学结构,该模型把对象用质点离散化,当一个质点在外力作用下发生运动时,产生的应力作用在其他相邻质点并带动相邻质点运动, 从而产生形变。传统质点弹簧模型中的弹簧是一种理想的弹簧,弹簧的伸长量与受力成正比,如果在软组织形变中采用传统的质点弹簧模型,很容易产生超弹性现象,影响形变的真实性。本实施例中采取限制体弹簧最大,最小阈值的方式解决弹簧超弹性问题,其中加入了体弹簧的质子弹簧结构如图2所示,其中的外框为最大阈值范围,内框为最小阈值范围。
本发明实施例的步骤S1中,基于待模拟的软组织三维数据和手术器械数据,构建上述软组织模型。
本发明实施例的步骤S2具体为:
S21、使用AABB包围盒获取软组织模型的三维边界信息;
其中,三维边界信息包括X轴最大值、X轴最小值/>、Y轴最大值/>、Y轴最小值/>、Z轴最大值/>、Z轴最小值/>以及X轴、Y轴和Z轴的边界差/>,/>和/>;
S22、通过二分法对软组织模型进行网格划分,获得细分的立方体;
其中,进行网格划分时,每次取当前立方体三维边界的X轴长度、Y轴长度和Z轴长度的二分之一进行网格划分;
S23、判断当前立方体的八个顶点是否均在原立方体对应的凸多边形内;
其中,原立方体为当前立方体对应的网格划分前的立方体;
若是,则进入步骤S24;
若否,则进入步骤S25;
S24、设置当前软组织模型体积V=Vn-1+Vn,进入步骤S27;
其中,Vn为当前立方体的体积,Vn-1为前一次网格划分时原立方体的体积,下标n为当前网格划分次数的计数值,软组织模型体积V的初始值为0,n为大于等于1的自然数;
S25、判断当前立方体边界差是否小于等于软组织模型边界差的四十分之一;
若是,则进入步骤S26;
若否,则返回步骤S22;
S26、设置当前立方体的体积为,当前软组织模型体积V= Vn-1+/>;
S27、重复步骤S22~S26,直到网格划分次数达到预设值或划分后的所有立方体精度达到预设值,将当前软组织模型体积V作为软组织模型的多尺度体积不变的体积。
本发明实施例的步骤S3具体为:
S31、遍历软组织模型的所有质点,并计算当前质点的质量;
S32、遍历当前质点的所有相邻质点,计算当前质点与其相邻质点间在单步形变后的距离,并计算该距离下当前质点所受的弹簧力;
S33、根据当前质点所受的阻尼力和弹簧力,计算当前质点所受周围相邻质点的合力,并根据当前质点的运动方程计算其在单步形变后的位置;
S34、在有限形变次数的形变过程中内,重复步骤S32~S33,在软组织模型多尺度体积不变的基础上,计算当前质点在最近两次单步形变的位置之差,当位置之差逼近零时,则当前质点形变结束,进而确定当前质点形变后的位置坐标。
在本实施例的步骤S31中,对于软组织模型的表面做均质处理,对于均质处理的米格三角面片,根据其面积把质量平均分配到三个顶点上,因此对于每个质点,其质量为:
式中,为曲面的面密度,通常为1,/>为当前质点i的一环邻域内的第j个三角面片的面积,j为质点i的一环邻域内的三角面片计数值,/>为三角片面的数量。
在本实施例的步骤S32中,假设一个质点弹簧***具有N个质点,每个质点同时受到内力和外力的作用,其中内力包括弹簧力和阻尼力,记一个弹簧S(i,k),该弹簧的两个端点质点记为i和k,在三维直角坐标系中的坐标分别为()和(/>),则质点所受周围相邻质点的合力/>为:
式中,为与质点i相邻的质点数,/>为质点i受到相邻质点k的弹簧力,为质点i受到相邻质点k的阻尼力,/>为质点i受到的外力。
本实施例的所述步骤S33中,当前质点i的运动方程为:
式中,为当前质点i的运动位移,/>为当前质点i的运动速度,/>为当前质点i所受周围相邻质点的合力,/>为当前质点i的质量。
在本实施例的步骤S34中,在有限形变次数的形变过程中,当当前质点的实际形变次数超过设定的最大形变次数时,则当前质点形变结束,进而确定当前质点形变后的位置坐标。
在本发明实施例中,在软组织形变模拟过程中,对于每个质点形变过程会实时发送到用户端的显示屏中展示软组织形变过程的实时渲染效果,并反馈质点所受力的信息。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (6)
1.一种基于体弹簧和多尺度体积不变的软组织形变模拟方法,其特征在于,包括以下步骤:
S1、构建软组织模型;
所述软组织模型由质子弹簧模型和虚拟体弹簧形成的若干四面体单元构成,每个所述四面体单元中包含3个虚拟体弹簧,用于获取软组织模型表面节点信息及内部体积信息;
S2、计算软组织模型多尺度体积不变的体积;
S3、基于软组织模型多尺度体积不变的体积,模拟软组织模型的形变过程,确定模型质点在形变过程中的位置坐标,实现软组织形变模拟;
所述步骤S2具体为:
S21、使用AABB包围盒获取软组织模型的三维边界信息;
其中,三维边界信息包括X轴最大值、X轴最小值/>、Y轴最大值/>、Y轴最小值/>、Z轴最大值/>、Z轴最小值/>以及X轴、Y轴和Z轴的边界差/>,/>和/>;
S22、通过二分法对软组织模型进行网格划分,获得细分的立方体;
其中,进行网格划分时,每次取当前立方体三维边界的X轴长度、Y轴长度和Z轴长度的二分之一进行网格划分;
S23、判断当前立方体的八个顶点是否均在原立方体对应的凸多边形内;
其中,原立方体为当前立方体对应的网格划分前的立方体;
若是,则进入步骤S24;
若否,则进入步骤S25;
S24、设置当前软组织模型体积V=Vn-1+Vn,进入步骤S27;
其中,Vn为当前立方体的体积,Vn-1为前一次网格划分时原立方体的体积,下标n为当前网格划分次数的计数值,软组织模型体积V的初始值为0,n为大于等于1的自然数;
S25、判断当前立方体边界差是否小于等于软组织模型边界差的四十分之一;
若是,则进入步骤S26;
若否,则返回步骤S22;
S26、设置当前立方体的体积为,当前软组织模型体积V= Vn-1+/>;
S27、重复步骤S22~S26,直到网格划分次数达到预设值或划分后的所有立方体精度达到预设值,将当前软组织模型体积V作为软组织模型的多尺度体积不变的体积。
2.根据权利要求1所述的基于体弹簧和多尺度体积不变的软组织形变模拟方法,其特征在于,所述步骤S3具体为:
S31、遍历软组织模型的所有质点,并计算当前质点的质量;
S32、遍历当前质点的所有相邻质点,计算当前质点与其相邻质点间在单步形变后的距离,并计算该距离下当前质点所受的弹簧力;
S33、根据当前质点所受的阻尼力和弹簧力,计算当前质点所受周围相邻质点的合力,并根据当前质点的运动方程计算其在单步形变后的位置;
S34、在有限形变次数的形变过程中内,重复步骤S32~S33,在软组织模型多尺度体积不变的基础上,计算当前质点在最近两次单步形变的位置之差,当位置之差逼近零时,则当前质点形变结束,进而确定当前质点形变后的位置坐标。
3.根据权利要求2所述的基于体弹簧和多尺度体积不变的软组织形变模拟方法,其特征在于,所述步骤S31中,当前质点的质量为:
式中,为曲面的面密度,/>为当前质点i的一环邻域内的第j个三角面片的面积,j为质点i的一环邻域内的三角面片计数值,/>为三角片面的数量。
4.根据权利要求2所述的基于体弹簧和多尺度体积不变的软组织形变模拟方法,其特征在于,所述步骤S32中,质点所受周围相邻质点的合力为:
式中,为与质点i相邻的质点数,/>为质点i受到相邻质点k的弹簧力,/>为质点i受到相邻质点k的阻尼力,/>为质点i受到的外力。
5.根据权利要求4所述的基于体弹簧和多尺度体积不变的软组织形变模拟方法,其特征在于,所述步骤S33中,当前质点的运动方程为:
式中,为当前质点i的运动位移,/>为当前质点i的运动速度,/>为当前质点i所受周围相邻质点的合力,/>为当前质点i的质量。
6.根据权利要求2所述的基于体弹簧和多尺度体积不变的软组织形变模拟方法,其特征在于,所述步骤S34中,在有限形变次数的形变过程中,当当前质点的实际形变次数超过设定的最大形变次数时,则当前质点形变结束,进而确定当前质点形变后的位置坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310841839.7A CN116561935B (zh) | 2023-07-11 | 2023-07-11 | 一种基于体弹簧和多尺度体积不变的软组织形变模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310841839.7A CN116561935B (zh) | 2023-07-11 | 2023-07-11 | 一种基于体弹簧和多尺度体积不变的软组织形变模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116561935A CN116561935A (zh) | 2023-08-08 |
CN116561935B true CN116561935B (zh) | 2023-10-03 |
Family
ID=87491898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310841839.7A Active CN116561935B (zh) | 2023-07-11 | 2023-07-11 | 一种基于体弹簧和多尺度体积不变的软组织形变模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116561935B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106156537A (zh) * | 2016-07-04 | 2016-11-23 | 南昌大学 | 基于麦夸特算法的径向基无网格软组织数据的力反馈模型建模方法 |
CN106504329A (zh) * | 2016-09-27 | 2017-03-15 | 西安科技大学 | 一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法 |
CN109389590A (zh) * | 2017-09-28 | 2019-02-26 | 上海联影医疗科技有限公司 | 结肠图像数据处理***及方法 |
CN110289103A (zh) * | 2019-06-19 | 2019-09-27 | 南京信息工程大学 | 软组织局部压缩形变模拟方法 |
CN111339654A (zh) * | 2020-02-24 | 2020-06-26 | 南京信息工程大学 | 一种虚拟手术***中软组织按压及形变恢复方法 |
CN112069933A (zh) * | 2020-08-21 | 2020-12-11 | 董秀园 | 基于体态识别和人体生物力学的骨骼肌肉受力估计方法 |
CN113409443A (zh) * | 2021-05-19 | 2021-09-17 | 南昌大学 | 一种基于位置约束和非线性弹簧的软组织建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2454384A (en) * | 2006-08-04 | 2009-05-06 | Auckland Uniservices Ltd | Biophysical virtual model database and applications |
-
2023
- 2023-07-11 CN CN202310841839.7A patent/CN116561935B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106156537A (zh) * | 2016-07-04 | 2016-11-23 | 南昌大学 | 基于麦夸特算法的径向基无网格软组织数据的力反馈模型建模方法 |
CN106504329A (zh) * | 2016-09-27 | 2017-03-15 | 西安科技大学 | 一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法 |
CN109389590A (zh) * | 2017-09-28 | 2019-02-26 | 上海联影医疗科技有限公司 | 结肠图像数据处理***及方法 |
CN110289103A (zh) * | 2019-06-19 | 2019-09-27 | 南京信息工程大学 | 软组织局部压缩形变模拟方法 |
CN111339654A (zh) * | 2020-02-24 | 2020-06-26 | 南京信息工程大学 | 一种虚拟手术***中软组织按压及形变恢复方法 |
CN112069933A (zh) * | 2020-08-21 | 2020-12-11 | 董秀园 | 基于体态识别和人体生物力学的骨骼肌肉受力估计方法 |
CN113409443A (zh) * | 2021-05-19 | 2021-09-17 | 南昌大学 | 一种基于位置约束和非线性弹簧的软组织建模方法 |
Non-Patent Citations (2)
Title |
---|
基于位置动力学的软组织建模研究;张禹;《北京生物医学工程》;第39卷(第4期);344-349 * |
面向虚拟手术的软组织建模及碰撞检测技术研究;赵瑞影;《中国优秀硕士学位论文全文数据库(医疗卫生科技辑)》(第8期);E080-41 * |
Also Published As
Publication number | Publication date |
---|---|
CN116561935A (zh) | 2023-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cotin et al. | Geometric and physical representations for a simulator of hepatic surgery | |
CN108648548A (zh) | 一种脑神经外科虚拟手术训练*** | |
Kim et al. | Haptic interaction and volume modeling techniques for realistic dental simulation | |
Zou et al. | A high-resolution model for soft tissue deformation based on point primitives | |
CN109192030A (zh) | 基于虚拟现实的真实腔镜微创手术仿真训练***及方法 | |
Kim et al. | Dental training system using multi-modal interface | |
Zhang et al. | A fast and accurate vascular tissue simulation model based on point primitive method | |
WO2006013813A1 (ja) | 情報処理装置およびプログラム | |
CN112587233B (zh) | 一种模拟气管内插管术的方法及*** | |
CN111026269A (zh) | 基于力反馈的生物组织结构的触觉反馈方法、装置和设备 | |
CN114333482B (zh) | 一种基于混合现实技术的虚拟解剖教学*** | |
CN116561935B (zh) | 一种基于体弹簧和多尺度体积不变的软组织形变模拟方法 | |
CN107342009B (zh) | 牙科备牙手术模拟方法及装置 | |
CN111276022A (zh) | 一种基于vr技术的胃镜模拟操作*** | |
CN113379929A (zh) | 一种基于物理仿真的骨组织修复虚拟现实解决方法 | |
CN117113452A (zh) | 一种肌肉模型及其构建方法和运动模拟方法 | |
Chen et al. | A displacement driven real-time deformable model for haptic surgery simulation | |
Wang et al. | Haptic rendering for dental training system | |
CN111105502A (zh) | 基于人工骨修复材料的生物肋骨鼻及下颌仿真整形技术 | |
CN114267219B (zh) | 虚拟手术仿真装置及仿真*** | |
KR20140050838A (ko) | 변형체 모델의 조각 시뮬레이션 방법 및 장치 | |
Guo et al. | A VR-based training system for vascular interventional surgery | |
Yan et al. | Real-time bone sawing interaction in orthopedic surgical simulation based on the volumetric object | |
Ye et al. | An improved collision detection and cutting algorithm of the soft tissue in virtual surgical simulation | |
Li et al. | A Method of Deformation Expression of Soft Tissue in Virtual Operation Based on Surface Mesh Model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |