CN108597612A - 模拟出血模型的虚拟切割算法 - Google Patents

模拟出血模型的虚拟切割算法 Download PDF

Info

Publication number
CN108597612A
CN108597612A CN201810358659.2A CN201810358659A CN108597612A CN 108597612 A CN108597612 A CN 108597612A CN 201810358659 A CN201810358659 A CN 201810358659A CN 108597612 A CN108597612 A CN 108597612A
Authority
CN
China
Prior art keywords
particle
denotes
pressure
tension
acceleration
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.)
Granted
Application number
CN201810358659.2A
Other languages
English (en)
Other versions
CN108597612B (zh
Inventor
张小瑞
徐千雄
孙伟
宋爱国
赵沁平
牛建伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Information Science and Technology
Original Assignee
Nanjing University of Information Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN201810358659.2A priority Critical patent/CN108597612B/zh
Publication of CN108597612A publication Critical patent/CN108597612A/zh
Application granted granted Critical
Publication of CN108597612B publication Critical patent/CN108597612B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开模拟出血模型的虚拟切割算法,该方法能够在不降低变形精度的情况下,通过改进后的网格模型,同时考虑紧度的影响和血小板对出血的影响,加入了张力、粘滞力加速度的计算方法,使出血模拟更加真实。本发明还通过选取不同的核函数加强了出血的实际行为和实时性,提高了粒子相关性质的计算效率。

Description

模拟出血模型的虚拟切割算法
技术领域
本发明涉及模拟出血模型的虚拟切割算法。
背景技术
近年来,随着虚拟现实技术的不断发展,通过虚拟现实平台模拟手术训练得到可能。虚拟手术给予医护人员极大的便利,他们可以使用该技术进行反复的模拟练习,提升自己的技术。目前的虚拟手术通过建模、渲染、计算将模拟手术的过程呈现在平台上,使用了有限元模型、无网格模型等。然而,目前虚拟手术技术考虑的方面仍不全面,真实性仍旧不够。在虚拟手术中增加出血模拟可以大大提高手术环境的视觉真实感,也可以培养外科医生对手术出血的反应,以及使医师把握手术力度与进程。
众所周知,血液是一种粘稠的不可压缩的液体,由于血小板和其他细胞的存在而具有独特的机械性能。通常采用Navier-Stokes(N-S)方程来研究血液的物理性质,采用光滑粒子流体动力学(SPH)方法通过核近似得到N-S方程的数值解。然而,传统SPH的恒定芯半径将导致模拟出血区域的边缘扭曲。导致出血模拟效果不佳。
为了解决这些问题,本文提出模拟出血模型的虚拟切割算法,该方法能够在不降低变形精度的情况下,通过改进后的网格模型实现模拟出血并能使手术医生把控手术进程,所得到的模型更加具有真实性。
发明内容
发明目的:本发明的目的是为了解决现有技术的不足,提供基于出血模型的虚拟切割算法。
技术方案:为了实现以上目的,本发明所述的模拟出血模型的虚拟切割算法,具体包括如下步骤:
第一步:
虚拟手术刀与虚拟器官接触且满足产生切口的条件时,虚拟器官被切割,产生切口,粒子i在出血模型内;否则只发生形变,粒子i在出血模型外;
如果粒子i在出血模型外,则只计算粒子i的压力fi p产生的粒子i的压力加速度
如果粒子i在出血模型内,则计算粒子i的压力fi p产生的压力加速度和粒子i的张力fi t产生的粒子i的张力加速度
第三步:
如果粒子i在出血模型内,则计算粒子i的粘性力fi v产生的粒子i的粘性力加速度
第四步:
将上述粒子i的压力、粘性力产生的加速度叠加以获得粒子i的总加速度ai
第五步:
在出血模拟中,在计算粒子i本身的特性时,用核函数的核心半径h作为参数计算粒子i的密度,使用POLY6核函数计算;计算压力和张力时,使用SPIKY 内核函数计算;计算粘性力时,使用VISCOSITY核函数计算;从而得到每帧中粒子i的加速度,进而计算下一帧中粒子i的位置。
2、根据权利要求1所述的模拟出血模型的虚拟切割算法,其特征在于:计算压力fi p产生的压力加速度前,先计算粒子i的密度,其方程为:
其中ρi表示粒子i的密度,mi表示粒子i的质量,j表示在以粒子i的核心半径为半径的圆的范围内的粒子的标号,W表示核函数,ri与rj分别表示粒子i和粒子j的位置,ri-rj表示粒子i与粒子j之间的距离,h表示核函数的核心半径。
再计算每个粒子i上的压力fi p,其方程为:
fi p=k(ρi0)
其中fi p表示粒子i上的压力,k是与温度有关的常数,ρi表示粒子i的密度,ρ0是静态密度;
根据权利要求1或2所述的模拟出血模拟的虚拟切割算法,其特征在于:压力fi p产生加速度其方程为其中表示粒子i的压力产生的压力加速度,fi p表示粒子i上的压力,mi表示粒子i的质量,j表示粒子i的核心半径中的粒子,pi和pj分别表示粒子i与粒子j的压强,ρj表示粒子j的密度,为哈密顿算子,W表示核函数,ri-rj表示粒子i与粒子j之间的距离,h表示核函数的核心半径。
步骤二所述的张力fi t,定义一个Colorfield区域来计算流体表面的张力 fi t;当有粒子i时,其值为1,没有粒子i时,其值为0;
如果粒子i在出血模型内部,那么周围的Colorfield是1,并且张力相互抵消;
如果粒子i在出血模型的表面上,则外部的Colorfield为0,并且张力影响粒子i,SPH被用于计算张力的方程为:
其中C(ri)表示用于计算粒子i的张力的SPH函数,mj表示粒子j的质量,ρj表示粒子j的密度,Aj是表示粒子j的Colorfield的值的函数。
梯度是粒子i表面的法线对应的参数构成的向量,即其中C表示 SPH函数;表面曲率ω用计算,其中是单位法向量的模;
将出血表面粒子i的张力定义为:
其中,δ是用于规范计算的常数,是出血面的单位法向量,是单位法向量的模。
所述的张力fi t产生的粒子i的张力加速度如下:
其中,表示粒子i上的张力fi t产生的加速度。
所述粒子i的粘性系数定义为:
其中,μi表示粒子i的粘性系数,λ和γ是系数,Pi为表示粒子i上的血小板性质的函数,并且当Pi的值减小时,粘度系数μi增加;Pi随时间变化的值按照SPH方法计算:
其中P为表示血小板性质的函数,t表示时间,η是传递系数,Pi和Pj分别定义粒子i和粒子j的初始血小板性质。
所述粒子i的粘性力产生的粘性力加速度如下:
其中,表示粒子i上的粘性力产生的粘性力加速度,vi和vj分别定义粒子 i和粒子j的初始粘性力性质。
所述的粒子i总的加速度其中,fi p表示压力,fi t表示张力,fi d表示重力,fi v表示粘性力。
所述第五步中,使用POLY6核函数计算粒子i本身的特性:
采用核心半径h作为粒子密度的函数,其中,Wpoly6表示POLY6核函数,r表示粒子的位置,π是圆周率。
所述采用SPIKY核函数计算压力及张力其中Wspiky表示SPIKY核函数。梯度计算如下:
所述VISCOSITY核函数用于计算粘性力:
它的拉普拉斯算子是:
有益效果:本发明所述的模拟出血模型的虚拟切割算法,与现有技术相比,具有以下优点:
1、本文方法为现有网格纳入了出血模型,其中组织和器官出血是在手术中常见的而且不可避免的现象在虚拟手术中增加出血模拟可以大大提高虚拟手术的视觉真实感;
2、本发明考虑紧度的影响和血小板对出血的影响,加入了张力、粘滞力加速度的甲酸方法,使出血模拟更加真实;
3、选取不同的核函数加强了出血的实际行为和实时性,提高了粒子相关性质的计算效率。
附图说明
图1为虚拟手术步骤图;
图2为切割图;
图3为网格模型图。
具体实施方式
下面结合附图,进一步阐明本发明。
实施例1
如图1-3所示的模拟出血模型的虚拟切割算法,具体包括如下步骤:
第一步:
虚拟手术刀与虚拟器官接触且满足产生切口的条件时,虚拟器官被切割,产生切口,粒子i在出血模型内;否则只发生形变,粒子i在出血模型外;
第二步:
2.1如果粒子i在出血模型外,则只计算粒子i的压力fi p产生的粒子i的压力加速度
计算压力fi p产生的压力加速度前,先计算粒子i的密度,其方程为:
其中ρi表示粒子i的密度,mi表示粒子i的质量,j表示在以粒子i的核心半径为半径的圆的范围内的粒子的标号,W表示核函数,ri与rj分别表示粒子i和粒子j的位置,ri-rj表示粒子i与粒子j之间的距离,h表示核函数的核心半径。
再计算每个粒子i上的压力fi p,其方程为:
fi p=k(ρi0)
其中fi p表示粒子i上的压力,k是与温度有关的常数,ρi表示粒子i的密度,ρ0是静态密度;
根据权利要求1或2所述的模拟出血模拟的虚拟切割算法,其特征在于:压力fi p产生加速度其方程为其中表示粒子i的压力产生的压力加速度,fi p表示粒子i上的压力,mi表示粒子i的质量,j表示粒子i的核心半径中的粒子,pi和pj分别表示粒子i与粒子j的压强,ρj表示粒子j的密度,为哈密顿算子,W表示核函数,ri-rj表示粒子i与粒子j之间的距离,h表示核函数的核心半径。
2.2如果粒子i在出血模型内,则计算粒子i的压力fi p产生的压力加速度和粒子i的张力fi t产生的粒子i的张力加速度
压力fi p产生加速度其方程为其中表示粒子i的压力产生的压力加速度,fi p表示粒子i上的压力,mi表示粒子i的质量,j表示粒子i的核心半径中的粒子,pi和pj分别表示粒子i与粒子j 的压强,为哈密顿算子,W表示核函数,ri-rj表示粒子i与粒子j之间的距离, h表示核函数的核心半径;
所述的张力fi t,定义一个Colorfield区域来计算流体表面的张力fi t;当有粒子i时,其值为1,没有粒子i时,其值为0;
如果粒子i在出血模型内部,那么周围的Colorfield是1,并且张力相互抵消;
如果粒子i在出血模型的表面上,则外部的Colorfield为0,并且张力影响粒子i,SPH被用于计算张力的方程为:
其中C(ri)表示用于计算粒子i的张力的SPH函数,mj表示粒子j的质量,ρj表示粒子j的密度,Aj是表示粒子j的Colorfield的值的函数。
梯度是粒子i表面的法线对应的参数构成的向量,即其中C表示 SPH函数;表面曲率ω用计算,其中是单位法向量的模;
将出血表面粒子i的张力定义为:
其中,δ是用于规范计算的常数,是出血面的单位法向量,是单位法向量的模。
所述的张力fi t产生的粒子i的张力加速度如下:
其中,表示粒子i上的张力fi t产生的加速度。
第三步:
如果粒子i在出血模型内,则计算粒子i的粘性力fi v产生的粒子i的粘性力加速度
粒子i的粘性系数定义为:
其中,μi表示粒子i的粘性系数,λ和γ是系数,Pi为表示粒子i上的血小板性质的函数,并且当Pi的值减小时,粘度系数μi增加;Pi随时间变化的值按照SPH方法计算:
其中P为表示血小板性质的函数,t表示时间,η是传递系数,Pi和Pj分别定义粒子i和粒子j的初始血小板性质。
所述粒子i的粘性力产生的粘性力加速度如下:
其中,表示粒子i上的粘性力产生的粘性力加速度,vi和vj分别定义粒子 i和粒子j的初始粘性力性质。
第四步:
将上述粒子i的压力、张力、重力、粘性力产生的加速度叠加以获得粒子i的总加速度ai
粒子i总的加速度其中,fi p表示压力,fi t表示张力, fi d表示重力,fi v表示粘性力;
第五步:
在出血模拟中,在计算粒子i本身的特性时,使用POLY6核函数计算;计算压力和张力时,使用SPIKY内核函数计算;计算粘性力时,使用VISCOSITY核函数计算;从而得到每帧中粒子i的加速度,进而计算下一帧中粒子i的位置。
所述第五步中,使用POLY6核函数计算粒子i本身的特性:
用核函数的核心半径h作为参数计算粒子i的密度,使用POLY6核函数计算,其中,Wpoly6表示POLY6核函数,r表示粒子的位置,π是圆周率。
采用SPIKY核函数计算压力及张力其中 Wspiky表示SPIKY核函数。梯度计算如下:
所述VISCOSITY核函数用于计算粘性力:
它的拉普拉斯算子是:

Claims (10)

1.一种模拟出血模型的虚拟切割算法,其特征在于:具体包括如下步骤:
第一步:
虚拟手术刀与虚拟器官接触且满足产生切口的条件时,虚拟器官被切割,产生切口,粒子i在出血模型内;否则只发生形变,粒子i在出血模型外;
第二步:
如果粒子i在出血模型外,则只计算粒子i的压力fi p产生的粒子i的压力加速度
如果粒子i在出血模型内,则计算粒子i的压力fi p产生的压力加速度和粒子i的张力fi t产生的粒子i的张力加速度
第三步:
如果粒子i在出血模型内,则计算粒子i的粘性力fi v产生的粒子i的粘性力加速度
第四步:
将粒子i的压力、张力、粘性力以及重力产生的加速度叠加以获得粒子i的总加速度ai
第五步:
在出血模拟中,在计算粒子i本身的特性时,用核函数的核心半径h作为参数计算粒子i的密度,使用POLY6核函数计算;计算压力和张力时,使用SPIKY内核函数计算;计算粘性力时,使用VISCOSITY核函数计算;从而得到每帧中粒子i的总加速度ai,进而计算下一帧中粒子i的位置。
2.根据权利要求1所述的模拟出血模型的虚拟切割算法,其特征在于:计算压力fi p产生的压力加速度前,先计算粒子i的密度,其方程为:其中,ρi表示粒子i的密度,mi表示粒子i的质量,j表示在以粒子i的核心半径为半径的圆的范围内的粒子的标号,W表示核函数,ri与rj分别表示粒子i和粒子j的位置,ri-rj表示粒子i与粒子j之间的距离,h表示核函数的核心半径;
再计算每个粒子i上的压力fi p,其方程为:
fi p=k(ρi0)
其中fi p表示粒子i上的压力,k是与温度有关的常数,ρi表示粒子i的密度,ρ0是静态密度。
3.根据权利要求1或2所述的模拟出血模拟的虚拟切割算法,其特征在于:压力fi p产生压力加速度其方程为其中表示粒子i的压力产生的压力加速度,fi p表示粒子i上的压力,mi表示粒子i的质量,j表示粒子i的核心半径中的粒子,pi和pj分别表示粒子i与粒子j的压强,ρj表示粒子j的密度,为哈密顿算子,W表示核函数,ri-rj表示粒子i与粒子j之间的距离,h表示核函数的核心半径。
4.根据权利要求1所述的模拟出血模型的虚拟切割算法,其特征在于:步骤二所述的张力fi t,定义一个Colorfield区域来计算流体表面的张力当有粒子i时,其值为1,没有粒子i时,其值为0;
如果粒子i在出血模型内部,那么周围的Colorfield是1,并且张力相互抵消;
如果粒子i在出血模型的表面上,则外部的Colorfield为0,并且张力影响粒子i,SPH被用于计算张力的方程为:其中C(ri)表示用于计算粒子i的张力的SPH函数,mj表示粒子j的质量,ρj表示粒子j的密度,Aj是表示粒子j的Colorfield的值的函数;
梯度是粒子i表面的法线对应的参数构成的向量,即其中C表示SPH函数;表面曲率ω用计算,其中是单位法向量的模;
将出血表面粒子i的张力定义为:
其中,δ是用于规范计算的常数,是出血面的单位法向量,是单位法向量的模。
5.根据权利要求1或4所述的模拟出血模型的虚拟切割算法,其特征在于:所述的张力fi t产生的粒子i的张力加速度如下:
其中,表示粒子i上的张力fi t产生的加速度。
6.根据权利要求1所述的模拟出血模型的虚拟切割算法,其特征在于:所述粒子i的粘性系数μi定义为:
λ>0,γ>0,其中,μi表示粒子i的粘性系数,λ和γ是系数,Pi为表示粒子i上的血小板性质的函数,并且当Pi的值减小时,粘度系数μi增加;Pi随时间变化的值按照SPH方法计算:
其中P为表示血小板性质的函数,t表示时间,η是传递系数,Pi和Pj分别定义粒子i和粒子j的初始血小板性质。
7.根据权利要求1或6所述的模拟出血模型的虚拟切割算法,其特征在于:所述粒子i的粘性力产生的粘性力加速度如下:
其中,表示粒子i上的粘性力产生的粘性力加速度,vi和vj分别定义粒子i和粒子j的初始粘性力性质。
8.根据权利要求1所述的一种基于出血模型的虚拟切割算法,其特征在于:所述的粒子i总的加速度其中,fi p表示压力,fi t表示张力,fi d表示重力,fi v表示粘性力。
9.根据权利要求1所述的模拟出血模型的虚拟切割算法,其特征在于:
所述第五步中,使用POLY6核函数计算粒子i本身的特性:
采用核心半径h作为粒子密度的函数,其中,Wpoly6表示POLY6核函数,r表示粒子的位置,π是圆周率。
10.根据权利要求1所述的模拟出血模型的虚拟切割算法,其特征在于:所述采用SPIKY核函数计算压力及张力其中Wspiky表示SPIKY核函数。梯度计算如下:
所述VISCOSITY核函数用于计算粘性力:
它的拉普拉斯算子是:
CN201810358659.2A 2018-04-19 2018-04-19 模拟出血模型的虚拟切割方法 Active CN108597612B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810358659.2A CN108597612B (zh) 2018-04-19 2018-04-19 模拟出血模型的虚拟切割方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810358659.2A CN108597612B (zh) 2018-04-19 2018-04-19 模拟出血模型的虚拟切割方法

Publications (2)

Publication Number Publication Date
CN108597612A true CN108597612A (zh) 2018-09-28
CN108597612B CN108597612B (zh) 2022-09-23

Family

ID=63614180

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810358659.2A Active CN108597612B (zh) 2018-04-19 2018-04-19 模拟出血模型的虚拟切割方法

Country Status (1)

Country Link
CN (1) CN108597612B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111047707A (zh) * 2019-11-11 2020-04-21 南昌大学 一种基于sph的混合粒子血液模型的流血仿真方法
CN113343513A (zh) * 2021-05-11 2021-09-03 南京信息工程大学 一种用于模拟软组织形变和路径切割的方法及装置
CN114220547A (zh) * 2021-11-23 2022-03-22 南昌大学 一种虚拟手术中流血与抽吸交互模型的建模方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103729555A (zh) * 2013-12-20 2014-04-16 深圳先进技术研究院 一种模拟血流与血管壁作用的方法和装置
CN104462830A (zh) * 2014-12-12 2015-03-25 武汉大学 一种基于gpu加速的混合粒子的实时血流血管耦合方法
CN105095534A (zh) * 2014-04-23 2015-11-25 北京冠生云医疗技术有限公司 对血管中血流进行模拟的方法和系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103729555A (zh) * 2013-12-20 2014-04-16 深圳先进技术研究院 一种模拟血流与血管壁作用的方法和装置
CN105095534A (zh) * 2014-04-23 2015-11-25 北京冠生云医疗技术有限公司 对血管中血流进行模拟的方法和系统
CN104462830A (zh) * 2014-12-12 2015-03-25 武汉大学 一种基于gpu加速的混合粒子的实时血流血管耦合方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111047707A (zh) * 2019-11-11 2020-04-21 南昌大学 一种基于sph的混合粒子血液模型的流血仿真方法
CN111047707B (zh) * 2019-11-11 2021-09-28 南昌大学 一种基于sph的混合粒子血液模型的流血仿真方法
CN113343513A (zh) * 2021-05-11 2021-09-03 南京信息工程大学 一种用于模拟软组织形变和路径切割的方法及装置
CN114220547A (zh) * 2021-11-23 2022-03-22 南昌大学 一种虚拟手术中流血与抽吸交互模型的建模方法

Also Published As

Publication number Publication date
CN108597612B (zh) 2022-09-23

Similar Documents

Publication Publication Date Title
Basdogan et al. VR-based simulators for training in minimally invasive surgery
US9330502B2 (en) Mixed reality simulation methods and systems
CN103400023B (zh) 软组织形变仿真方法
Qin et al. Learning blood management in orthopedic surgery through gameplay
CN108597612B (zh) 模拟出血模型的虚拟切割方法
CN105205867B (zh) 一种微创虚拟腹主动脉血管手术中的碰撞检测方法
CN106875462A (zh) 一种基于元球模型和混合驱动方法的实时数字器官切割方法
Pan et al. Graphic and haptic simulation system for virtual laparoscopic rectum surgery
De et al. Physically realistic virtual surgery using the point-associated finite field (PAFF) approach
CN109192030A (zh) 基于虚拟现实的真实腔镜微创手术仿真训练系统及方法
CN111488670B (zh) 一种非线性的质点弹簧软组织形变仿真方法
CN110289104B (zh) 软组织按压和形变恢复的模拟方法
CN107992672A (zh) 一种基于虚拟弹簧的软组织形变建模方法
CN103729555A (zh) 一种模拟血流与血管壁作用的方法和装置
Wang et al. A review of virtual cutting methods and technology in deformable objects
KR100551201B1 (ko) 볼륨 모델 기반의 햅틱 인터페이스를 이용한 치과 치료훈련 및 평가 시스템
Banihani et al. POD for real‐time simulation of Hyperelastic soft biological tissue using the point collocation method of finite spheres
Radetzky et al. Elastodynamic shape modeler: A tool for defining the deformation behavior of virtual tissues
CN114267219B (zh) 虚拟手术仿真装置及仿真系统
CN114693892A (zh) 一种虚拟手术中切割流血模型的建模方法
Aloisio et al. Computer-based simulator for catheter insertion training
Furqan et al. Surface-based virtual dental surgical simulator using haptic display
De Paolis et al. Virtual model of the human brain for neurosurgical simulation
Ogata et al. A development of surgical simulator for training of operative skills using patient-specific data
CN109033641B (zh) 一种基于硅胶愈合模型的虚拟切割算法

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