CN106650251B - 一种针灸力反馈形变模型的建模方法 - Google Patents
一种针灸力反馈形变模型的建模方法 Download PDFInfo
- Publication number
- CN106650251B CN106650251B CN201611152098.8A CN201611152098A CN106650251B CN 106650251 B CN106650251 B CN 106650251B CN 201611152098 A CN201611152098 A CN 201611152098A CN 106650251 B CN106650251 B CN 106650251B
- Authority
- CN
- China
- Prior art keywords
- soft tissue
- coordinate
- model
- deformation
- particle
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT 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
Abstract
本发明公开了一种针灸力反馈形变模型的建模方法,首先重建软组织模型;然后计算用于记录软组织模型的全局形变信息的全局形变矩阵;接着计算针刺点形变前的初始坐标和针刺点坐标之间的距离,并据此判定针尖是否刺破表皮;如果针尖未刺破表皮,实时产生虚拟弹簧计算皮肤表层的力反馈;如果针尖已刺破表皮,则计算软组织模型的内部力反馈,同时对针刺后的皮肤进行复原控制。本发明通过虚拟弹簧和体元模型可以分别计算针在皮肤内部和外部的力反馈,在确保计算简单的同时确保高精度的变形模拟,在交互过程中,操作者可以感受到触觉信息的实时性和真实感。
Description
技术领域
本发明涉及针灸模拟,尤其涉及一种针灸力反馈形变模型的建模方法。
背景技术
针灸是传统中国养生治病的神奇医术,在现代养生医学中也有很广泛的应用。然而人体穴位比较难找,对于刚接触针灸的新手来说是一个很大的挑战,几乎无法实战练习,而扎错穴位也会对人体产生或大或小的影响,因此训练至关重要。传统的针灸培训成本高、风险高,虚拟针灸模拟系统可以帮助中医新手练习针灸,有效地解决这一问题。
柔性体表面形变是虚拟触觉交互过程的关键。常用的物理模型中,弹簧-质点模型具有建模简单、计算灵活、速度快的优点,但该模型精度有限,稳定性差。有限元模型精度高,但计算相对复杂。因此,具有良好的交互功能和形变效果的柔体模型对于虚拟手术系统来说至关重要。
发明内容
本发明所要解决的技术问题是针对背景技术的不足,提供了一种针灸力反馈形变模型的建模方法。
本发明为解决上述技术问题采用以下技术方案:
一种针灸力反馈形变模型的建模方法,包含以下步骤:
步骤1),重建软组织模型:采用体元模型重建软组织内部,采用网格模型重建软组织表皮;所述体元模型包含用于表示软组织内部的内部体元和用于表示模型针头的针体元;进行针灸模拟时,将模型针头在软组织内部所经过路径上的内部体元替换成针体元;
步骤2),根据以下公式计算用于记录软组织模型的全局形变信息的的全局形变矩阵A:
其中,n是软组织模型的网格质点数,mi是软组织模型中质点i的质量,i为大于等于1小于等于n的自然数,pi是软组织模型中质点i的坐标,pcm是形变后软组织模型的质心坐标,是质点i的初始坐标,是初始状态下软组织模型的质心坐标;
步骤3),计算针刺点形变前的初始坐标和针刺点坐标pc之间的距离
当且针尖坐标x与针刺点坐标pc重合时,判定皮肤表层为未被刺破的状态,执行步骤4),其中,dmax为预设的皮肤最大形变量;
当且针尖坐标x与针刺点坐标pc不重合时,判定皮肤表层为被刺破的状态,执行步骤5);
步骤4),根据以下公式实时产生虚拟弹簧计算皮肤表层的力反馈Fout:
其中,σ(i)为软组织模型中受形变影响的网格质点集合,k为预设的虚拟弹簧系数,为软组织模型中质点i的初始坐标和当前坐标pi之间的距离,Nout为以针尖为顶点沿针方向的单位向量;
步骤5),根据针刺深度计算软组织模型的内部力反馈,同时对针刺后的皮肤进行复原控制:
步骤5.1),根据针刺点坐标pc和针尖坐标x的距离din计算软组织模型的内部力反馈Fin:
Fin=ρ·din·Nin
其中,ρ是预设的的皮下组织密度,din是针刺点坐标pc和针尖坐标x的距离,Nin是从针尖坐标x到针刺点坐标pc的单位向量;
步骤5.2),根据以下公式计算软组织模型中各个质点目标位置的坐标:
其中,为质点i目标位置的坐标,R是旋转矩阵,S是对称系数,矩阵AT是矩阵A的转置矩阵;
步骤5.3),根据以下公式优化软组织模型中各质点目标位置的坐标:
其中,ηi为预设的恢复系数,ηi∈[0,1];
步骤5.4),对于软组织模型中的各个质点i,控制其在预设的速率Rr下以ηi=ηi-Rr的速度向其目标位置移动。
本发明可以在保证较好地表现针刺形变和恢复效果的前提下实现快速、逼真、实时的力反馈计算,与现有技术相比,具有以下技术效果:
1.使用体元和网格模型建模从而保证模型的稳定性,并且方便力反馈和形变的计算。
2.计算全局形变矩阵A以保证模型的全局稳定性,防止模型走样,并未之后的恢复算法做准备。
3.计算目标位置gi实现皮肤被针刺破后的恢复,使针灸模型更为逼真。
4.插入恢复系数ηi以控制表面的恢复速率,进一步模拟真实情景下的针灸形变,保证形变恢复的稳定性和真实性。
5.由针头x和针刺点pc之间的位置以及针刺点pc和针刺点初始位置的距离为根据进行刺破判定,可以快速找到针头刺破皮肤的时间点,并以此为分割线进行表面形变的模拟和力反馈的模拟,逼真地展现在两种不同状态下的针灸模型。
6.未刺破前使用虚拟弹簧计算外部力反馈,使用虚拟弹簧计算力反馈方便快捷,计算时通过形变位置各顶点与其初始位置间的距离的改变,保证了力反馈的刷新。
7.刺破后考虑皮下组织密度,进行内部力反馈的计算,让模型的力反馈更为逼真,力反馈随针刺入的深度实时刷新,满足了系统的实时性。
附图说明
图1是针灸力反馈形变模型的仿真流程图;
图2是改良的形状匹配算法示例图;
图3是刺破前手臂横截面示例图;
图4是刺破后手臂横截面示例图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
如图1所示,本发明公开了一种针灸力反馈形变模型的建模方法,包含以下步骤:
步骤1),采用体元模型重建软组织内部,采用网格模型重建软组织表皮;所述体元模型包含用于表示软组织内部的内部体元和用于表示模型针头的针体元;进行针灸模拟时,将模型针头在软组织内部所经过路径上的内部体元替换成针体元。
步骤2),计算全局形变矩阵A,全局形变矩阵A记录软组织的全局形变信息,也是形变恢复算法的基础,全局形变矩阵A需要获取软组织模型初始状态下的各质点坐标及质心坐标,在整个针灸模拟的过程中实时计算更新,根据以下公式计算软组织内部的全局形变矩阵A:
其中,n是软组织模型的网格质点数,mi是软组织模型中质点i的质量,i为大于等于1小于等于n的自然数,pi是软组织模型中质点i的坐标,pcm是形变后软组织模型的质心坐标,是质点i的初始坐标,是初始状态下软组织模型的质心坐标。
步骤3),进行刺破判定,计算针刺点形变前的初始坐标和针刺点坐标pc之间的距离其中,dmax为预设的皮肤最大形变量;
当且针尖坐标x与针刺点坐标pc重合时,视为皮肤表层未被刺破的状态,执行步骤4);
当且针尖坐标x与针刺点坐标pc不重合时,视为皮肤表层被刺破的状态,执行步骤5)。
步骤4),皮肤未刺破,根据皮肤形变实时产生虚拟弹簧计算皮肤表层的力反馈,如图3所示,皮肤表层的力反馈Fout计算如下:
其中,σ(i)为软组织模型受形变影响的网格质点集合,k为预设的虚拟弹簧系数,为软组织模型中质点i初始坐标和当前坐标pi之间的距离,Nout为以针尖为顶点沿针方向的单位向量。
步骤5),皮肤已刺破,根据针刺深度计算软组织模型的内部力反馈,同时对针刺后的皮肤形变进行优化,使皮肤在刺破后一定程度地恢复原状,如图4所示。
步骤5.1),计算软组织模型的内部接触力,通过针刺点坐标pc和针尖坐标x的距离din的改变生成持续且变化的力反馈,内部力反馈Fin计算如下:
Fin=ρ·din·Nin
其中,ρ是预设的皮下组织密度,din是针刺点坐标pc和针尖坐标x的距离,Nin是从针尖坐标x到针刺点坐标pc的单位向量。
计算刺破后的皮肤形变,当皮肤组织被刺破后,皮肤会一定程度地恢复形变,表现为各质点i向其对应的目标位置gi移动,从而恢复初始形状,如图2所示。
步骤5.2),目标位置gi是皮肤形变恢复的关键,通过全局形变矩阵A求解目标位置gi的坐标根据以下公式计算各质点的目标位置gi的坐标
其中R是旋转矩阵,S是对称系数,R和S都是求解目标位置gi的坐标所需要的重要参数,矩阵AT是矩阵A的转置矩阵;
在实际操作时各质点对应的目标位置gi位置固定,各质点i会快速回复形变,造成模拟的失真,为了防止软组织表面过快恢复,得到逼真的针灸形变效果,插入恢复系数ηi,ηi∈[0,1],如图4所示,通过控制目标位置gi随皮肤形变而移动,以控制各质点i恢复形变的速率,从而控制表面软组织的恢复速率。
步骤5.3),根据以下公式优化软组织模型中各质点目标位置的坐标:
步骤5.4),对于软组织模型中的各个质点i,控制其在预设的速率Rr下以ηi=ηi-Rr的速度向其目标位置移动。
当ηi=0时,质点i会直接移向目标位置gi;当ηi=1时,代表质点i不会立刻移动,i将在用户定义的速率Rr下逐渐以ηi=ηi-Rr的速度向目标位置gi移动。
本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种针灸力反馈形变模型的建模方法,其特征在于,包含以下步骤:
步骤1),重建软组织模型:采用体元模型重建软组织内部,采用网格模型重建软组织表皮;所述体元模型包含用于表示软组织内部的内部体元和用于表示模型针头的针体元;进行针灸模拟时,将模型针头在软组织内部所经过路径上的内部体元替换成针体元;
步骤2),根据以下公式计算用于记录软组织模型的全局形变信息的全局形变矩阵A:
其中,n是软组织模型的网格质点数,mi是软组织模型中质点i的质量,i为大于等于1小于等于n的自然数,pi是软组织模型中质点i的坐标,pcm是形变后软组织模型的质心坐标,是质点i的初始坐标,是初始状态下软组织模型的质心坐标;
步骤3),计算针刺点形变前的初始坐标和针刺点坐标pc之间的距离
当且针尖坐标x与针刺点坐标pc重合时,判定皮肤表层为未被刺破的状态,执行步骤4),其中,dmax为预设的皮肤最大形变量;
当且针尖坐标x与针刺点坐标pc不重合时,判定皮肤表层为被刺破的状态,执行步骤5);
步骤4),根据以下公式实时产生虚拟弹簧计算皮肤表层的力反馈Fout:
其中,σ(i)为软组织模型中受形变影响的网格质点集合,k为预设的虚拟弹簧系数,为软组织模型中质点i的初始坐标和当前坐标pi之间的距离,Nout为以针尖为顶点沿针方向的单位向量;
步骤5),根据针刺深度计算软组织模型的内部力反馈,同时对针刺后的皮肤进行复原控制;
步骤5.1),根据针刺点坐标pc和针尖坐标x的距离din计算软组织模型的内部力反馈Fin:
Fin=ρ·din·Nin
其中,ρ是预设的皮下组织密度,din是针刺点坐标pc和针尖坐标x的距离,Nin是从针尖坐标x到针刺点坐标pc的单位向量;
步骤5.2),根据以下公式计算软组织模型中各个质点目标位置的坐标:
其中,为质点i目标位置的坐标,R是旋转矩阵,S是对称系数,矩阵AT是矩阵A的转置矩阵;
步骤5.3),根据以下公式优化软组织模型中各质点目标位置的坐标:
其中,ηi为预设的恢复系数,ηi∈[0,1];
步骤5.4),对于软组织模型中的各个质点i,控制其在预设的速率Rr下以ηi=ηi-Rr的速度向其目标位置移动。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611152098.8A CN106650251B (zh) | 2016-12-14 | 2016-12-14 | 一种针灸力反馈形变模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611152098.8A CN106650251B (zh) | 2016-12-14 | 2016-12-14 | 一种针灸力反馈形变模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106650251A CN106650251A (zh) | 2017-05-10 |
CN106650251B true CN106650251B (zh) | 2018-08-14 |
Family
ID=58823416
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611152098.8A Expired - Fee Related CN106650251B (zh) | 2016-12-14 | 2016-12-14 | 一种针灸力反馈形变模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106650251B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107808698B (zh) * | 2017-11-08 | 2021-10-15 | 哈尔滨理工大学 | 一种柔性针穿刺软组织的力学及弯曲变形建模方法 |
CN109243575B (zh) * | 2018-09-17 | 2022-04-22 | 华南理工大学 | 一种基于移动交互和增强现实的虚拟针灸方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262699A (zh) * | 2011-07-27 | 2011-11-30 | 华北水利水电学院 | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 |
CN103793552A (zh) * | 2013-12-18 | 2014-05-14 | 哈尔滨工程大学 | 一种软组织形变的局部质点弹簧模型的实时动态生成方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102044086B (zh) * | 2010-11-30 | 2012-07-25 | 华北水利水电学院 | 一种软组织形变仿真方法 |
WO2013163762A1 (en) * | 2012-05-02 | 2013-11-07 | The Royal Institution For The Advancement Of Learning/Mcgill University | Device for soft tissue support and method for anchoring |
CN103400023B (zh) * | 2013-06-28 | 2016-11-02 | 华北水利水电大学 | 软组织形变仿真方法 |
CN105513130B (zh) * | 2016-02-01 | 2018-02-16 | 福建师范大学福清分校 | 一种基于网格与无网格混合的软组织形变方法 |
-
2016
- 2016-12-14 CN CN201611152098.8A patent/CN106650251B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102262699A (zh) * | 2011-07-27 | 2011-11-30 | 华北水利水电学院 | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 |
CN103793552A (zh) * | 2013-12-18 | 2014-05-14 | 哈尔滨工程大学 | 一种软组织形变的局部质点弹簧模型的实时动态生成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106650251A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Alterovitz et al. | Needle insertion and radioactive seed implantation in human tissues: Simulation and sensitivity analysis | |
CN105913718A (zh) | 一种埋线提升整形手术仿真系统 | |
CN106650251B (zh) | 一种针灸力反馈形变模型的建模方法 | |
CN106200982A (zh) | 一种基于虚拟现实技术的医学手术训练系统 | |
CN102968811B (zh) | 微创血管介入手术中导丝的实时运动仿真方法 | |
CN105069301B (zh) | 支持力觉交互的腰椎穿刺虚拟仿真与训练系统 | |
CN105825752A (zh) | 基于力反馈设备的虚拟眼角膜手术培训系统 | |
CN106408035B (zh) | 基于人力触觉感知特性的力触觉再现真实感客观评价方法 | |
CN107433021A (zh) | 一种基于镜像神经元的vr康复系统 | |
CN204029245U (zh) | 一种手术模拟系统 | |
Goksel et al. | 3D simulation of needle-tissue interaction with application to prostate brachytherapy | |
CN105138223A (zh) | 一种基于指套式力触觉交互装置的物体柔软度再现方法 | |
Fortmeier et al. | Image-based palpation simulation with soft tissue deformations using chainmail on the GPU | |
CN105719526A (zh) | 基于力反馈的埋线提眉整形手术模拟系统 | |
CN110211661A (zh) | 基于混合现实的手功能训练系统及数据处理方法 | |
CN106528993A (zh) | 基于由碟形弹簧片构成的组合弹簧虚拟模型建模方法 | |
CN104063057A (zh) | 一种支持增强力触觉反馈的涡卷弹簧仿真系统 | |
CN106570341A (zh) | 一种基于无网格径向基数据拟合的软组织受力形变模型建模方法 | |
CN107221032A (zh) | 虚拟腹腔镜肝囊肿切除开窗手术教学方法和系统 | |
CN103413485B (zh) | 一种用于骨髓穿刺术教学模具及其使用方法 | |
CN207886596U (zh) | 一种基于镜像神经元的vr康复系统 | |
Cenydd et al. | VCath: a tablet-based neurosurgery training tool | |
CN107038756B (zh) | 一种三维人体模型的旋转系统及方法 | |
Makhlouf et al. | A new TLM algorithm to solve the Pennes's equation for dosimetry applications | |
CN206021740U (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180814 Termination date: 20211214 |
|
CF01 | Termination of patent right due to non-payment of annual fee |