CN114557770B - 微波消融手术的交互式路径规划方法及系统 - Google Patents
微波消融手术的交互式路径规划方法及系统 Download PDFInfo
- Publication number
- CN114557770B CN114557770B CN202210169199.5A CN202210169199A CN114557770B CN 114557770 B CN114557770 B CN 114557770B CN 202210169199 A CN202210169199 A CN 202210169199A CN 114557770 B CN114557770 B CN 114557770B
- Authority
- CN
- China
- Prior art keywords
- ablation
- area
- ablated
- feasible
- point
- 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
- 238000002679 ablation Methods 0.000 title claims abstract description 181
- 230000002452 interceptive effect Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004088 simulation Methods 0.000 claims abstract description 32
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000003628 erosive effect Effects 0.000 claims abstract description 7
- 230000011218 segmentation Effects 0.000 claims description 24
- 238000011156 evaluation Methods 0.000 claims description 20
- 238000001356 surgical procedure Methods 0.000 claims description 13
- 230000005672 electromagnetic field Effects 0.000 claims description 12
- 238000012986 modification Methods 0.000 claims description 7
- 230000004048 modification Effects 0.000 claims description 7
- 238000011158 quantitative evaluation Methods 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 210000000056 organ Anatomy 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000006907 apoptotic process Effects 0.000 claims description 2
- 238000005485 electric heating Methods 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 2
- 230000001105 regulatory effect Effects 0.000 abstract description 2
- 206010028980 Neoplasm Diseases 0.000 description 3
- 208000014018 liver neoplasm Diseases 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 206010019695 Hepatic neoplasm Diseases 0.000 description 1
- 238000011298 ablation treatment Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000001640 apoptogenic effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 201000007270 liver cancer Diseases 0.000 description 1
- 210000005228 liver tissue Anatomy 0.000 description 1
- 238000000554 physical therapy Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000007674 radiofrequency ablation Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
- A61B18/18—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by applying electromagnetic radiation, e.g. microwaves
- A61B18/1815—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by applying electromagnetic radiation, e.g. microwaves using microwaves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
- A61B2018/00571—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body for achieving a particular surgical effect
- A61B2018/00577—Ablation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/101—Computer-aided simulation of surgical operations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/101—Computer-aided simulation of surgical operations
- A61B2034/105—Modelling of the patient, e.g. for ligaments or bones
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/107—Visualisation of planned trajectories or target regions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2065—Tracking using image or pattern recognition
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- General Physics & Mathematics (AREA)
- Robotics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Electromagnetism (AREA)
- Otolaryngology (AREA)
- Surgical Instruments (AREA)
Abstract
本发明提供了一种微波消融手术的交互式路径规划方法及系统,包括输入CT/MR影像数据,自动生成粗分割的待消融区域,交互确认待消融区域,获得精细待消融区域;根据获得的精细待消融区域,使用形态学中的侵蚀算法,获得消融区域内部的可行点区域,根据设定条件推算可行消融点位置;获得暂定消融点;根据暂定消融点,对最终产生的消融区域进行仿真计算,动态预测消融区域;评估并规划消融路径。本发明通过影像信息的引入,半自动地确定可行消融点,减少了医生规划路径的工作量。通过仿真模块的使用,提供了对消融结果更精细的预估,提高了消融路径规划的可靠性。通过交互式设计,使得每个关键流程都有调整确认,提高了路径规划的安全可靠性。
Description
技术领域
本发明涉及医疗领域,具体地,涉及微波消融手术的交互式路径规划方法及系统。
背景技术
微波消融手术是目前肿瘤治疗方法中最先进的物理治疗方法之一,传统微波消融手术中,术前医生需要根据病人的CT图像和经验确定消融点和入针路径规划,缺乏对规划结果的定量评估。另一方面,传统路径规划算法往往将消融范围假设为固定的形状,缺乏对于不同消融位置产生的消融范围变化的动态预测,易对规划产生不利影响。
专利文献为CN113456219A的发明专利公开了一种基于CT图像的肝癌射频消融手术路径规划方法及装置,方法包括:(1)数据预处理:对规划相关的器官模型提取分离并进行适当抽样,计算规划所需的相关点集,其中靶向点点集与皮肤点集之间的连接线构成原始的全部解空间;(2)路径筛选:对这组解空间,根据临床约束、是否能够完全覆盖肿瘤、是否是帕累托最优前沿部分等条件依次进行筛选;(3)计算消融区位置:应用整数规划找出剩余路径中最佳消融区位置,并通过设置的权重评分得出最优解并输出。但是上述方案使用预设的消融区域形状假设,没有考虑到不同位置消融区域形状的差异。鲍楠公开的论文:CT引导的肺穿刺路径规划与手术导航关键技术研究中,除了影像信息外,还依赖病人体表标记物提供位置信息,需要额外的CT扫描。张睿公开的论文:CT引导肝肿瘤热消融治疗手术路径规划方法与系统研究中,同样使用预设的消融区域形状假设,没有考虑到不同位置消融区域形状的差异。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种微波消融手术的交互式路径规划方法及系统。
根据本发明提供的一种微波消融手术的交互式路径规划方法,包括如下步骤:
待消融区域设置步骤:输入CT/MR影像数据,自动生成粗分割的待消融区域,交互确认待消融区域,获得精细待消融区域;
可行消融点规划步骤:根据获得的精细待消融区域,使用形态学中的侵蚀算法,获得消融区域内部的可行点区域,根据设定条件推算可行消融点位置;获得暂定消融点;
消融结果动态评估步骤:根据暂定消融点,对最终产生的消融区域进行仿真计算,动态预测消融区域;
交互式消融点调整步骤:根据可行消融点规划步骤和消融结果动态评估步骤的结果,迭代式地给予反馈,每次修改消融点之后将重新执行消融结果动态评估步骤,获得消融效果仿真结果和定量评估;
消融路径规划步骤:使用可行消融点和皮肤上的可行入针点的连线构成消融路径的原始解空间;将消融结果动态评估步骤中获得的预测消融区域旋转,获得解空间内不同可行解的仿真区域粗预测,排除掉无法覆盖待消融区域的路径,依照临床硬约束筛选出所有的可行路径,确定最终消融路径规划。
优选地,所述设定条件包括临床硬约束条件、是否能大致覆盖消融区域条件。
优选地,待消融区域设置步骤中,交互式分割使用MIDeepSeg分割算法,具体步骤包括:
步骤S1:输入影像数据,手动交互点击待消融区域的边界点,设定带分割区域的整体范围,算法中的预测网络计算粗分割结果并可视化;
步骤S2:在粗分割结果周围,对假阳性区域点击负面点,假阴性区域点击正面点,算法获得交互输入和粗分割结果,重新计算修正后的分割区域。
步骤S3:如果满意分割结果则结束,不满意则重复步骤S2直至满意。
优选地,消融结果动态评估步骤中,基于仿真模块,对最终产生的消融区域进行仿真计算,动态预测消融区域。
优选地,采用时域有限差分进行微波的仿真,使用有限差分法进行传热的仿真,交替计算电磁场和热场以求得电热耦合场的分布,具体步骤如下:
步骤F1:将待消融区域和所在脏器建模,根据不同组织类型设置不同的电/热参数,根据人体温度输入初始温度值;
步骤F2:根据激励源的设置开始计算,基于时域有限差分迭代求取当前时间步长下的电磁场分布情况,获得电磁场分布后计算焦耳热,焦耳热功率为其中J为电流体密度,E为电场强度;将焦耳热输入Pennes方程中以仿真有内热源场景下的生物热传导:
其中Qj为电磁场以焦耳热的形式产生的热量;
步骤F3:根据多次迭代后生成的结果,使用Arrhenius方程计算细胞累积收到的热损伤,当损伤超过阈值之后,细胞凋亡,修改所在位置电/热参数;
步骤F4:迭代计算完成后,计算预测消融区域和待消融区域的重叠程度,量化计算为Dice分数和覆盖率。
优选地,Dice分数计算如下:
优选地,覆盖率计算如下:
优选地,所述损伤阈值默认为1。
优选地,电/热参数包括热导率、比热容、相对介电常数、电导率、相对磁导率。
根据本发明提供的一种微波消融手术的交互式路径规划系统,包括如下模块:
待消融区域设置模块:输入CT/MR影像数据,自动生成粗分割的待消融区域,交互确认待消融区域,获得精细待消融区域;
可行消融点规划模块:根据获得的精细待消融区域,使用形态学中的侵蚀算法,获得消融区域内部的可行点区域,根据设定条件推算可行消融点位置;获得暂定消融点;
消融结果动态评估模块:根据暂定消融点,对最终产生的消融区域进行仿真计算,动态预测消融区域;
交互式消融点调整模块:根据可行消融点规划和消融结果动态评估的结果,迭代式地给予反馈,每次修改消融点之后将重新执行消融结果动态评估,获得消融效果仿真结果和定量评估;
消融路径规划模块:使用可行消融点和皮肤上的可行入针点的连线构成消融路径的原始解空间;将消融结果动态评估中获得的预测消融区域旋转,获得解空间内不同可行解的仿真区域粗预测,排除掉无法覆盖待消融区域的路径,依照临床硬约束筛选出所有的可行路径,确定最终消融路径规划。
与现有技术相比,本发明具有如下的有益效果:
1、本发明通过影像信息的引入,半自动地确定可行消融点,减少了医生规划路径的工作量。
2、本发明通过仿真模块的使用,提供了对消融结果更精细的预估,提高了消融路径规划的可靠性。
3、本发明通过交互式设计,使得每个关键流程都有医生调整和确认,提高了路径规划的安全性和可靠性。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为微波消融手术的交互式路径规划方法步骤流程图。
图2为微波消融手术的交互式路径规划方法交互界面示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1和图2所示,根据本发明提供的一种微波消融手术的交互式路径规划方法及系统,本发明提供的方法具体包括:根据影像信息,半自动规划可行消融点;基于仿真模块,动态预测消融区域,提供可视化的结果估计;交互式地接受医生反馈,迭代执行规划和仿真,使规划结果更好地服务于医生的经验和知识。
进一步具体说明,根据本发明提供的一种微波消融手术的交互式路径规划方法,包括以下4步:
S1:待消融区域设置。输入CT/MR影像数据,自动生成粗分割的待消融区域,提供交互窗口供医生进行待消融区域的修正和确认,最终获得影像数据中的精细待消融区域。
其中,交互式地进行分割使用的是MIDeepSeg,一种深度学习的交互式分割算法,具体使用步骤如下:
S1-D1:输入影像数据,手动交互点击待消融区域的边界点,设定带分割区域的整体范围,算法中的预测网络计算粗分割结果并可视化。
S1-D2:在粗分割结果周围,对假阳性区域点击负面点,假阴性区域点击正面点,算法获得交互输入和粗分割结果,重新计算修正后的分割区域。
S1-D3:如果满意分割结果则结束,不满意则重复S1-D2直至满意。
S2:可行消融点规划。根据S1中获得的精细待消融区域,使用形态学中的侵蚀算法,获得消融区域内部的可行点区域,依次根据临床硬约束、是否能大致覆盖消融区域等条件,推算可行消融点位置。将可行消融点位置可视化,供医生交互式地选择和微调消融点设置,最终获得经过医生确认的暂定消融点。
S3:消融结果动态评估。根据S2中获得的暂定消融点,基于仿真模块,对最终产生的消融区域进行仿真计算,动态预测消融区域。根据预测消融区域和待消融区域的重叠情况输出可视化结果和定量指标。
其中,为了兼顾预测的精度和速度,我们使用了时域有限差分进行微波的仿真,使用有限差分法进行传热的仿真,交替计算电磁场和热场以求得电热耦合场的分布,具体步骤如下:
S3-F1:将待消融区域和所在脏器建模,根据不同组织类型设置不同的电/热参数(热导率、比热容、相对介电常数、电导率、相对磁导率等),然后根据人体温度输入初始温度值(默认为37℃)。
S3-F2:根据激励源的设置开始计算,基于时域有限差分迭代求取当前时间步长下的电磁场分布情况,获得电磁场分布后计算焦耳热,焦耳热功率为其中J为电流体密度,E为电场强度。将焦耳热输入Pennes方程中以仿真有内热源场景下的生物热传导:
其中Qj为电磁场以焦耳热的形式产生的热量。
S3-F3:根据多次迭代后生成的结果,使用Arrhenius方程计算细胞累积收到的热损伤,当损伤超过阈值(默认为1)之后,细胞凋亡,修改所在位置电/热参数。以肝脏组织为例,具体电/热参数值如下:
S3-F4:迭代计算完成后,计算预测消融区域和待消融区域的重叠程度,量化计算为Dice分数和覆盖率,计算公式如下:
S4:交互式消融点调整。医生根据S2和S3的结果,迭代式地给予反馈,每次修改消融点之后将重新执行S3步骤,获得消融效果仿真结果和定量评估。医生可不断根据结果和指标的指示进行调整和修改,直至最终得到医生确认的可行消融点。
S5:消融路径规划。使用可行消融点和皮肤上的可行入针点的连线构成消融路径的原始解空间;将S3中获得的预测消融区域旋转,获得解空间内不同可行解的仿真区域粗预测,排除掉无法覆盖待消融区域的路径,然后依照临床硬约束筛选出所有的可行路径。经过医生的交互确定最终消融路径规划。
本发明还提供一种微波消融手术的交互式路径规划系统,包括待消融区域设置模块:输入CT/MR影像数据,自动生成粗分割的待消融区域,交互确认待消融区域,获得精细待消融区域;可行消融点规划模块:根据获得的精细待消融区域,使用形态学中的侵蚀算法,获得消融区域内部的可行点区域,根据设定条件推算可行消融点位置;获得暂定消融点;消融结果动态评估模块:根据暂定消融点,对最终产生的消融区域进行仿真计算,动态预测消融区域;交互式消融点调整模块:根据可行消融点规划和消融结果动态评估的结果,迭代式地给予反馈,每次修改消融点之后将重新执行消融结果动态评估,获得消融效果仿真结果和定量评估;消融路径规划模块:使用可行消融点和皮肤上的可行入针点的连线构成消融路径的原始解空间;将消融结果动态评估中获得的预测消融区域旋转,获得解空间内不同可行解的仿真区域粗预测,排除掉无法覆盖待消融区域的路径,依照临床硬约束筛选出所有的可行路径,确定最终消融路径规划。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以,本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构;也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (10)
1.一种微波消融手术的交互式路径规划方法,其特征在于,包括如下步骤:
待消融区域设置步骤:输入CT/MR影像数据,自动生成粗分割的待消融区域,交互确认待消融区域,获得精细待消融区域;
可行消融点规划步骤:根据获得的精细待消融区域,使用形态学中的侵蚀算法,获得消融区域内部的可行点区域,根据设定条件推算可行消融点位置;获得暂定消融点;
消融结果动态评估步骤:根据暂定消融点,对最终产生的消融区域进行仿真计算,动态预测消融区域;
交互式消融点调整步骤:根据可行消融点规划步骤和消融结果动态评估步骤的结果,迭代式地给予反馈,每次修改消融点之后将重新执行消融结果动态评估步骤,获得消融效果仿真结果和定量评估;
消融路径规划步骤:使用可行消融点和皮肤上的可行入针点的连线构成消融路径的原始解空间;将消融结果动态评估步骤中获得的预测消融区域旋转,获得解空间内不同可行解的仿真区域粗预测,排除掉无法覆盖待消融区域的路径,依照临床硬约束筛选出所有的可行路径,确定最终消融路径规划。
2.根据权利要求1所述的微波消融手术的交互式路径规划方法,其特征在于,所述设定条件包括临床硬约束条件、是否能大致覆盖消融区域条件。
3.根据权利要求1所述的微波消融手术的交互式路径规划方法,其特征在于,待消融区域设置步骤中,交互式分割使用MIDeepSeg分割算法,具体步骤包括:
步骤S1:输入影像数据,手动交互点击待消融区域的边界点,设定带分割区域的整体范围,算法中的预测网络计算粗分割结果并可视化;
步骤S2:在粗分割结果周围,对假阳性区域点击负面点,假阴性区域点击正面点,算法获得交互输入和粗分割结果,重新计算修正后的分割区域;
步骤S3:如果满意分割结果则结束,不满意则重复步骤S2直至满意。
4.根据权利要求1所述的微波消融手术的交互式路径规划方法,其特征在于,消融结果动态评估步骤中,基于仿真模块,对最终产生的消融区域进行仿真计算,动态预测消融区域。
5.根据权利要求4所述的微波消融手术的交互式路径规划方法,其特征在于,采用时域有限差分进行微波的仿真,使用有限差分法进行传热的仿真,交替计算电磁场和热场以求得电热耦合场的分布,具体步骤如下:
步骤F1: 将待消融区域和所在脏器建模,根据不同组织类型设置不同的电/热参数,根据人体温度输入初始温度值;
步骤F2:根据激励源的设置开始计算,基于时域有限差分迭代求取当前时间步长下的电磁场分布情况,获得电磁场分布后计算焦耳热,焦耳热功率为,其中/>为电流体密度,/>为电场强度;将焦耳热输入Pennes方程中以仿真有内热源场景下的生物热传导:
其中为电磁场以焦耳热的形式产生的热量;
步骤F3:根据多次迭代后生成的结果,使用Arrhenius方程计算细胞累积收到的热损伤,当损伤超过阈值之后,细胞凋亡,修改所在位置电/热参数;
步骤F4:迭代计算完成后,计算预测消融区域和待消融区域的重叠程度,量化计算为Dice分数和覆盖率。
6.根据权利要求5所述的微波消融手术的交互式路径规划方法,其特征在于,Dice分数计算如下:
。
7.根据权利要求5所述的微波消融手术的交互式路径规划方法,其特征在于,覆盖率计算如下:
。
8.根据权利要求5所述的微波消融手术的交互式路径规划方法,其特征在于,所述损伤阈值默认为1。
9.根据权利要求5所述的微波消融手术的交互式路径规划方法,其特征在于,电/热参数包括热导率、比热容、相对介电常数、电导率、相对磁导率。
10.一种微波消融手术的交互式路径规划系统,其特征在于,包括如下模块:
待消融区域设置模块:输入CT/MR影像数据,自动生成粗分割的待消融区域,交互确认待消融区域,获得精细待消融区域;
可行消融点规划模块:根据获得的精细待消融区域,使用形态学中的侵蚀算法,获得消融区域内部的可行点区域,根据设定条件推算可行消融点位置;获得暂定消融点;
消融结果动态评估模块:根据暂定消融点,对最终产生的消融区域进行仿真计算,动态预测消融区域;
交互式消融点调整模块:根据可行消融点规划和消融结果动态评估的结果,迭代式地给予反馈,每次修改消融点之后将重新执行消融结果动态评估,获得消融效果仿真结果和定量评估;
消融路径规划模块:使用可行消融点和皮肤上的可行入针点的连线构成消融路径的原始解空间;将消融结果动态评估中获得的预测消融区域旋转,获得解空间内不同可行解的仿真区域粗预测,排除掉无法覆盖待消融区域的路径,依照临床硬约束筛选出所有的可行路径,确定最终消融路径规划。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210169199.5A CN114557770B (zh) | 2022-02-23 | 2022-02-23 | 微波消融手术的交互式路径规划方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210169199.5A CN114557770B (zh) | 2022-02-23 | 2022-02-23 | 微波消融手术的交互式路径规划方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114557770A CN114557770A (zh) | 2022-05-31 |
CN114557770B true CN114557770B (zh) | 2023-10-03 |
Family
ID=81714157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210169199.5A Active CN114557770B (zh) | 2022-02-23 | 2022-02-23 | 微波消融手术的交互式路径规划方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114557770B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024032570A1 (zh) * | 2022-08-08 | 2024-02-15 | 武汉联影智融医疗科技有限公司 | 一种介入规划系统、方法、装置和存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103717167A (zh) * | 2011-07-28 | 2014-04-09 | 皇家飞利浦有限公司 | 消融规划系统 |
CN109077804A (zh) * | 2018-08-19 | 2018-12-25 | 天津大学 | 一种基于ct图像的微波消融治疗计划方法 |
CN110151309A (zh) * | 2018-02-14 | 2019-08-23 | 上海交通大学 | 多模态消融治疗术前规划方法及其设备 |
CN111012484A (zh) * | 2020-01-06 | 2020-04-17 | 南京康友医疗科技有限公司 | 一种实时消融区域成像系统 |
CN113456219A (zh) * | 2021-06-15 | 2021-10-01 | 艾瑞迈迪医疗科技(北京)有限公司 | 基于ct图像的肝癌射频消融手术路径规划方法及装置 |
-
2022
- 2022-02-23 CN CN202210169199.5A patent/CN114557770B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103717167A (zh) * | 2011-07-28 | 2014-04-09 | 皇家飞利浦有限公司 | 消融规划系统 |
CN110151309A (zh) * | 2018-02-14 | 2019-08-23 | 上海交通大学 | 多模态消融治疗术前规划方法及其设备 |
CN109077804A (zh) * | 2018-08-19 | 2018-12-25 | 天津大学 | 一种基于ct图像的微波消融治疗计划方法 |
CN111012484A (zh) * | 2020-01-06 | 2020-04-17 | 南京康友医疗科技有限公司 | 一种实时消融区域成像系统 |
CN113456219A (zh) * | 2021-06-15 | 2021-10-01 | 艾瑞迈迪医疗科技(北京)有限公司 | 基于ct图像的肝癌射频消融手术路径规划方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN114557770A (zh) | 2022-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103209653B (zh) | 用于概率消融规划的系统 | |
US10607738B2 (en) | System and method for minimally invasive thermal ablation treatment planning | |
US11908584B2 (en) | Methods and systems for modeling a necrotized tissue volume in an ablation procedure | |
CN103800075B (zh) | 用于肝脏肿瘤消融术的患者特定建模的系统和方法 | |
CN114007538B (zh) | 基于深度学习的规划组织消融的方法 | |
Payne et al. | Image-based multi-scale modelling and validation of radio-frequency ablation in liver tumours | |
CN107427280A (zh) | 用于基于弹性成像监测的自适应消融和治疗的系统和方法 | |
WO2014133924A1 (en) | System and method for interactive patient specific simulation of radiofrequency ablation therapy | |
Pavliha et al. | Planning of electroporation-based treatments using web-based treatment-planning software | |
Memon et al. | Segmentation of lungs from CT scan images for early diagnosis of lung cancer | |
Chen et al. | Semiautomatic radiofrequency ablation planning based on constrained clustering process for hepatic tumors | |
CN114557770B (zh) | 微波消融手术的交互式路径规划方法及系统 | |
Audigier et al. | Lattice Boltzmann method for fast patient-specific simulation of liver tumor ablation from CT images | |
CN110049715B (zh) | 用于可视化患者的内脏的设备及相关联的可视化方法 | |
US20240058062A1 (en) | System and method for ablation treatment of tissue with interactive guidance | |
McCreedy et al. | Radio frequency ablation registration, segmentation, and fusion tool | |
CN114299009A (zh) | 基于医学图像的消融区域确定方法、设备及存储介质 | |
CN113196411A (zh) | 用于辅助提供用于消融处置的模板处置参数的系统 | |
CN114699166A (zh) | 路径规划方法和可读存储介质 | |
WO2021029973A1 (en) | Methods and systems for modeling a necrotized tissue volume in an ablation procedure | |
Preusser et al. | Patient-specific planning for radio-frequency ablation of tumors in the presence of uncertainty | |
Deshazer et al. | Developing an open platform for evidence-based microwave ablation treatment planning and validation | |
CN116936033B (zh) | 一种经皮消融手术规划方法、装置、介质及设备 | |
Frackowiak et al. | Model-based hepatic percutaneous microwave ablation planning. First validation on a clinical dataset | |
Li et al. | Image processing and modeling for active needle steering in liver surgery |
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 | ||
PP01 | Preservation of patent right |
Effective date of registration: 20240621 Granted publication date: 20231003 |