CN111529051B - 一种预测电脉冲消融区域的系统 - Google Patents
一种预测电脉冲消融区域的系统 Download PDFInfo
- Publication number
- CN111529051B CN111529051B CN202010301885.4A CN202010301885A CN111529051B CN 111529051 B CN111529051 B CN 111529051B CN 202010301885 A CN202010301885 A CN 202010301885A CN 111529051 B CN111529051 B CN 111529051B
- Authority
- CN
- China
- Prior art keywords
- electric field
- pulse
- field intensity
- ablation
- electrode needle
- 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 211
- 230000005684 electric field Effects 0.000 claims abstract description 265
- 238000011282 treatment Methods 0.000 claims abstract description 94
- 238000000034 method Methods 0.000 claims abstract description 29
- 238000010276 construction Methods 0.000 claims abstract description 16
- 238000009826 distribution Methods 0.000 claims description 37
- 230000014509 gene expression Effects 0.000 claims description 25
- 238000003062 neural network model Methods 0.000 claims description 14
- 230000036278 prepulse Effects 0.000 claims description 12
- 230000004044 response Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 10
- 238000012512 characterization method Methods 0.000 claims description 8
- 238000000926 separation method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 11
- 230000008569 process Effects 0.000 abstract description 7
- 230000006378 damage Effects 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000000875 corresponding effect Effects 0.000 description 92
- 230000006870 function Effects 0.000 description 41
- 238000010586 diagram Methods 0.000 description 13
- 238000004520 electroporation Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 10
- 210000004027 cell Anatomy 0.000 description 10
- 206010028980 Neoplasm Diseases 0.000 description 6
- 239000013598 vector Substances 0.000 description 6
- 210000000170 cell membrane Anatomy 0.000 description 5
- 230000002427 irreversible effect Effects 0.000 description 5
- 238000013528 artificial neural network Methods 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 4
- 230000008823 permeabilization Effects 0.000 description 4
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 4
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 3
- 230000002441 reversible effect Effects 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 2
- 230000030833 cell death Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 239000012528 membrane Substances 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 239000000232 Lipid Bilayer Substances 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000010317 ablation therapy Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000002421 cell wall Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
Images
Classifications
-
- 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/04—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating
- A61B18/12—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating by passing a current through the tissue to be heated, e.g. high-frequency current
-
- 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
- A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
- A61B2018/00636—Sensing and controlling the application of energy
- A61B2018/00773—Sensed parameters
- A61B2018/00827—Current
-
- 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/00636—Sensing and controlling the application of energy
- A61B2018/00773—Sensed parameters
- A61B2018/00892—Voltage
-
- 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/00636—Sensing and controlling the application of energy
- A61B2018/00898—Alarms or notifications created in response to an abnormal condition
Landscapes
- Health & Medical Sciences (AREA)
- Surgery (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Otolaryngology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Plasma & Fusion (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种预测电脉冲消融区域的系统，包括：绘制点选取模块，用于基于电极针裸露长度范围，沿电极针长度方向选取多个绘制点；二维区域获取模块，用于针对每个选取的绘制点，基于约束条件，从每个选取的绘制点对应的电场强度数据库中获取所述每个选取的绘制点对应的目标电场强度等值线；三维区域构建模块，用于基于所述每个选取的绘制点对应的目标电场强度等值线，构建沿电极针长度方向的三维模型。本发明提出的预测系统能够适应不同患者的个体差异，能够在治疗前快速有效地预测电脉冲消融区域，提高术前预测的准确性和完整性，有助于治疗过程中的合理有效布针，保证治疗计划的顺利实施，最大限度地确保治疗效果，减少正常组织消融与损伤。
Description
技术领域
本发明属于医疗器械技术领域，具体涉及一种预测电脉冲消融区域的系统。
背景技术
癌症是危害人类健康的主要疾病。传统的癌症治疗方法和新近发展起来的以微创消融为 特征的热消融疗法，由于受适应症、禁忌症、治疗副作用、热效应等因素的限制，使得其临 床应用受到一定的局限性。近年来，随着脉冲生物电学的不断发展，电场脉冲以其非热、微 创的生物医学效用受到广泛关注，并逐渐应用于肿瘤的临床治疗。电穿孔是一种利用电脉冲 来克服细胞膜这一屏障的过程，即在电场作用下，细胞壁组成和结构发生一定的变化，细胞 膜脂双层上形成瞬时孔洞，细胞膜的通透性和膜电导瞬时增大，在正常情况下不能通过细胞 膜的病毒颗粒、DNA、蛋白质以及染料颗粒等得以进出细胞的一个生物物理过程。电场造成细 胞的穿孔有两种形式：可逆电穿孔和不可逆电穿孔。如果脉冲条件不超过某一临界限度，细 胞膜的通透性是可逆的，属于可逆电穿孔；如果脉冲条件超过某一临界限度，细胞会遭到不 可逆的损伤甚至死亡，属于不可逆电穿孔。在不可逆电穿孔技术应用于肿瘤的临床治疗方法 中，利用电脉冲消融生物组织，是通过对靶区细胞发送电场脉冲，使细胞内外的离子运动并 聚集在外膜两侧，引起跨膜电位的急剧变化，细胞发生不可逆电穿孔，打破细胞内外平衡， 最终引起细胞死亡。
为了保证治疗方法的顺利实施，最大限度地确保治疗效果以及减少对正常组织的消融与 损伤，在治疗前对电脉冲消融生物组织区域进行预测是必要的。由于生物组织电学特性的差 异，使得对于不同的组织，需要电脉冲消融的肿瘤区域各不相同，可见，在治疗前对电脉冲 消融区域进行预测也一直是该治疗方法实施前的难点。目前电脉冲消融区域的预测主要集中 在二维消融区域的预测，依靠研究人员的经验，一方面通过动物实验、切片对比等方法来预 测消融区域的大小，另一方面通过对术后的核磁影响等方法来判断消融的效果。由于每位患 者病情各不相同，病灶各异，治疗方法也各不相同，单纯依靠经验来预测消融区域，而缺乏 适用于不同患者的通用预测手段，因此难以在治疗前针对患者的个体情况制定出合理有效的 治疗计划。
发明内容
为了解决上述缺乏适用于不同患者的通用预测手段的技术问题，本发明实施例提出了一 种预测电脉冲消融区域的系统，能够适用于不同的患者，在治疗前快速有效地预测三维电脉 冲消融区域，提高术前预测的准确性和完整性，并且对于治疗过程中的布针，保证治疗计划 的顺利实施，最大限度地确保治疗效果，以及减少正常组织消融与损伤至关重要。
本发明提出的预测电脉冲消融区域的系统，包括：
绘制点选取模块，用于基于电极针裸露长度范围，沿电极针长度方向选取多个绘制点；
二维区域获取模块，用于针对每个选取的绘制点，基于约束条件，从每个选取的绘制点 对应的电场强度数据库中获取所述每个选取的绘制点对应的目标电场强度等值线；以及，
三维区域构建模块，用于基于所述每个选取的绘制点对应的目标电场强度等值线，构建 沿电极针长度方向的三维模型，所述构建的三维模型内的区域作为预测的电脉冲消融区域。
在某些实施例中，所述约束条件包括第一约束条件和第二约束条件；
所述二维区域获取模块包括：集合获取模块，用于针对每个选取的绘制点，基于第一约 束条件，从所述每个选取的绘制点对应的电场强度数据库中获取所述每个选取的绘制点对应 的电场强度等值线集合；以及，目标获取模块，用于针对每个选取的绘制点，基于第二约束 条件，从所述每个选取的绘制点对应的电场强度等值线集合中获取所述每个选取的绘制点对 应的目标电场强度等值线。
在某些实施例中，所述第一约束条件包括电场强度消融阈值。
在某些实施例中，所述电场强度消融阈值通过拟合函数或神经网络模型确定；其中，所 述拟合函数包括用于表示所述拟合函数的表达式的拟合变量和用于选择拟合函数的表达式的 条件变量。
在某些实施例中，所述拟合函数的条件变量包括子脉冲宽度和电极针裸露长度；和/或，
所述拟合函数的拟合变量包括电脉冲的电场强度、电脉冲中包括的脉冲串个数和消融区 域的电导率比率，或者，所述拟合变量包括电脉冲的脉冲电压、电极针间距、电脉冲中包括 的脉冲串个数和消融区域的电导率比率。
在某些实施例中，响应于所述拟合变量包括电脉冲的电场强度、电脉冲中包括的脉冲串 个数和消融区域的电导率比率，所述拟合函数的表达式为：
Eth＝a1*E+b1*N+c1*R+d1*E*N+e1*E*R+f1*N*R+g1*E*N*R+h1，其中，Eth表示电场强度消融阈值，E表示电脉冲的电场强度，N表示电脉冲中包括的脉冲串个数，R表示消融区域的电导率比率， a1～h1表示消融区域预测参数；
或者，响应于所述拟合变量包括电脉冲的脉冲电压、电极针间距、电脉冲中包括的脉冲 串个数和消融区域的电导率比率，所述拟合函数的表达式为：
Eth＝a2*U+b2*N+c2*D+d2*R+e2*U*N+f2*U*D+g2*U*R+h2*N*D+i2*N*R+j2*L*R+k2*N*D*R+l2*U*N*D +m2*U*N*R+n2*U*D*R+o2*N*D*R+p2，其中，Eth表示电场强度消融阈值，U表示电脉冲的脉冲电 压，D表示电极针间距，N表示电脉冲中包括的脉冲串个数，R表示消融区域的电导率比率， a2～p2为消融区域预测参数。
在某些实施例中，所述消融区域预测参数基于历史数据采用拟合的方式确定，所述约束 条件与治疗计划相关。
在某些实施例中，所述选取多个绘制点是从预设的绘制点中选取多个绘制点，所述预设 的绘制点对应的电场强度数据库的构建包括：基于约束方程，以参数组合模型中的参数组合 为单位获取电场强度分布参数组；基于所述获取的电场强度分布参数组和每个预设的绘制点， 构建所述每个预设的绘制点对应的电场强度数据库。
在某些实施例中，所述约束方程包括电场方程、电导率方程和边界条件方程。
在某些实施例中，构建所述参数组合模型的构建参数包含：消融区域的第一电导率比率、 第一电脉冲的脉冲电压、第一电极针裸露长度和电极针间距。
在某些实施例中，所述电场强度等值线集合中的电场强度等值线采用卡西尼曲线或样条 拟合曲线表示。
在某些实施例中，所述二维区域获取模块还包括：坐标变换模块，用于响应于布针位置 不在同一水平线上，对所述目标电场强度等值线进行坐标变换。
在某些实施例中，所述二维区域获取模块还包括：区域叠加模块，用于响应于多次治疗， 叠加针对每次治疗所获取的目标电场强度等值线内的区域，作为所述多次治疗的预测电脉冲 消融区域。
在某些实施例中，所述第二约束条件包括第二电脉冲的脉冲电压、消融区域的第二电导 率比率、第二电极针裸露长度和电极针坐标。
在某些实施例中，所述消融区域的第二电导率比率通过释放预脉冲的方法获取。
在某些实施例中，所述释放预脉冲的方法，包括：对每对电极针施加预脉冲至电流稳定， 基于一表征参数在所述每对电极针的稳定状态和在所述每对电极针的初始状态的数值比，获 取所述消融区域的第二电导率比率；或者，包括：对每对电极针施加两次预脉冲至电流稳定， 基于一表征参数在所述每对电极针的两次稳定状态的数值比，获取所述消融区域的第二电导 率比率。
在某些实施例中，所述表征参数包括电流或电阻。
在某些实施例中，所述基于一表征参数在所述每对电极针的两次稳定状态的数值比，获 取所述消融区域的第二电导率比率，包括：基于一表征参数在所述每对电极针的两次稳定状 态的数值比，通过查询列表的方式获取所述消融区域的第二电导率比率；所述列表表征在不 同的电极针裸露长度时的电极针间距、所述表征参数的数值比和消融区域的电导率比率的关 系。
在某些实施例中，所述基于所述每个选取的绘制点对应的目标电场强度等值线，构建沿 电极针长度方向的三维模型，包括：基于所述每个选取的绘制点对应的目标电场强度等值线， 以网格的形式构建沿电极针长度方向的三维模型。
在某些实施例中，所述网格在每对相邻的绘制点对应的目标电场强度等值线之间构建。
在某些实施例中，选取[-L,L]作为所述电极针裸露长度范围，其中，L表示电极针裸露 长度。
本发明的有益效果：本发明实施例提出的预测电脉冲消融区域的系统，可以在治疗前， 通过约束条件对每个选取的绘制点对应的电场强度数据库进行约束，从而可以获得每个选取 的绘制点对应的特定的目标电场强度等值线，并由此构建出电脉冲消融区域的三维模型，从 而能够适应不同患者的个体差异，并且能够在治疗前快速有效地预测电脉冲消融区域，提高 术前预测的准确性和完整性，而且有助于治疗过程中的合理有效布针，保证治疗计划的顺利 实施，最大限度地确保治疗效果，以及减少正常组织消融与损伤。
附图说明
图1示出本发明实施例提出的预测电脉冲消融区域的系统的功能框图；
图2示出本发明实施例提出的预测电脉冲消融区域的系统所采用的三维坐标系；
图3示出本发明实施例提出的预测电脉冲消融区域的系统确定生物组织的电场强度消融 阈值采用的神经网络模型的示意图；
图4示出本发明实施例提出的预测电脉冲消融区域的系统中电场强度等值线的示意图；
图5a和图5b示出本发明实施例提出的预测电脉冲消融区域的系统中匹配电场强度等值 线的卡西尼曲线的示意图；
图6a和图6b示出本发明实施例提出的预测电脉冲消融区域的系统中电场强度等值线的 连通情形的示意图；
图7a示出本发明实施例提出的预测电脉冲消融区域的系统中样条拟合曲线在第一象限 的示意图；
图7b示出本发明实施例提出的预测电脉冲消融区域的系统中样条拟合曲线的示意图；
图8示出本发明实施例提出的预测电脉冲消融区域的系统中施加同一预设电场强度时的 初始电流和稳定电流示意图；
图9示出本发明实施例提出的预测电脉冲消融区域的系统中施加不同预设电场强度时的 稳定电流示意图；
图10示出本发明实施例提出的预测电脉冲消融区域的系统中坐标变换的示意图；
图11示出本发明实施例提出的预测电脉冲消融区域的系统中多次治疗的预测电脉冲消 融区域叠加的示意图；
图12a、图12b、图12c和图12d示出本发明实施例提出的预测电脉冲消融区域的系统中 构建三维区域的示意图；
图13示出本发明实施例提出的预测电脉冲消融区域的系统的一具体应用场景的结构示 意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图， 对本发明进一步详细说明。但本领域技术人员知晓，本发明并不局限于附图和以下实施例。
如前所述，目前单纯依靠研究人员的经验来预测消融区域，因此在治疗前针对患者的个 体情况制定出合理有效的治疗计划存在较大难度，并且对于治疗计划可能获得的治疗效果也 无法科学预见。基于此，本发明实施例提出了一种预测电脉冲消融区域的系统，通过在治疗 前，对电场强度数据库进行阶段性约束，从而可以从电场强度数据库中获取到不同约束条件 下的预测电脉冲消融区域，为合理有效制定治疗计划提供一定的依据，并且有助于在治疗前 科学地预见的治疗计划可能获得的治疗效果。
参考图1所示，本发明实施例提出的一种预测电脉冲消融区域的系统，包括：
绘制点选取模块，用于基于电极针裸露长度范围，沿电极针长度方向选取多个绘制点；
二维区域获取模块，用于针对每个选取的绘制点，基于约束条件，从每个选取的绘制点 对应的电场强度数据库中获取所述每个选取的绘制点对应的目标电场强度等值线；以及，
三维区域构建模块，用于基于所述每个选取的绘制点对应的目标电场强度等值线，构建 沿电极针长度方向的三维模型，所述构建的三维模型内的区域作为预测的电脉冲消融区域。
在本发明实施例中，采用电场强度等值线表征电脉冲消融区域的二维平面形状，在此基 础上增加电极针裸露长度作为建模维度，用以表征电脉冲消融的电能量于电极针裸露段向四 周传递而形成的三维消融区域，并且通过约束条件对基于电极针裸露长度范围选取的每个选 取的绘制点对应的电场强度数据库进行约束，以获得每个选取的绘制点对应的特定的目标电 场强度等值线，并由此构建出电脉冲消融区域的三维模型，从而能够适应不同患者的个体差 异，并且能够在治疗前快速有效地预测电脉冲消融区域，提高术前预测的准确性和完整性， 有助于治疗过程中的合理有效布针，保证治疗计划的顺利实施，最大限度地确保治疗效果， 以及减少正常组织消融与损伤。
在一实施例中，所述约束条件包括第一约束条件和第二约束条件。所述二维区域获取模 块包括：
集合获取模块，用于基于第一约束条件，从电场强度数据库中获取电场强度等值线集合； 以及，
目标获取模块，用于基于第二约束条件，从所述获取的电场强度等值线集合中获取目标 电场强度等值线，所述获取的目标电场强度等值线内的区域作为预测电脉冲消融区域。
通过所述第一约束条件和所述第二约束条件对所述电场强度数据库进行阶段性约束，从 而可以快速有效地获得特定的目标电场强度等值线。
在一实施例中，所述系统还包括绘制模块和显示模块，绘制模块用于绘制所述预测电脉 冲消融区域，显示模块用于显示所述预测电脉冲消融区域。
下面对本发明实施例提出的方案进行具体说明，其中的示例仅是为了更清楚地说明本发 明实施例，不应理解为对本发明实施例的限制。
关于电极针裸露长度范围
由于电脉冲消融的电能量于电极针裸露段向四周传递，所以电脉冲消融区域长度(消融体 积的长、宽或高的最大值)一般会大于电极针裸露长度，在一实施例中，为准确覆盖电脉冲消 融区域,可以选取[-L,L]作为所述电极针裸露长度范围，其中，L表示电极针裸露长度。所述 三维模型的三维坐标系可以采用xyz三维坐标系，所述电极针长度方向与所述三维坐标系的 z轴方向一致，所述电场强度等值线所在的平面为与Z轴垂直的平面，即为xy平面或与xy 平面平行的平面，从而可以针对每个选取的绘制点(可以理解为电极针裸露长度范围内的任 何一个Z轴上的点)获取其对应的目标电场强度等值线。在使用电穿孔技术进行肿瘤的临床 治疗中，一般是采用一对电极针同时作用于生物组织，因此在一可选实施例中，参考图2所 示，可以以这两个电极针的裸露段的中点的连线的中点作为三维模型的三维坐标系的坐标原 点，可以使得构建出的三维模型相对于坐标原点对称，所述电场强度等值线相对于坐标原点 对称，在某些应用场景中，可以有效减少计算量。
所述沿电极针长度方向选取多个绘制点，可以沿电极针长度方向有规律地选取多个绘制 点，例如采用等间隔的绘制点选取方式，或者间隔距离从中心向外逐渐增加或逐渐减小的绘 制点选取方式，或者遵循某种函数关系的绘制点选取方式；另外，也可以沿电极针长度方向 无规律地选取多个绘制点。
所述选取的绘制点的数量可以自行确定，可以理解的是，一般来说，选取的绘制点的数 量越多，预测的准确性越高，计算量越大。在实际选取时，可综合考虑计算量和准确性等因 素选取合理的绘制点的数量。在某些应用场景中，可以根据计划选取的绘制点的数量选择合 理的绘制点选取方式，以尽量覆盖三维模型的外部轮廓特征。绘制点的数量和绘制点选取方 式可以由用户决定，也可以由系统决定，例如采用系统默认值，或者由系统基于电极针裸露 长度范围等因素来决定。
在本发明实施例中，在用电场强度等值线表征电脉冲消融区域的二维平面形状的基础上， 增加电极针裸露长度作为建模维度，从而能够表征电脉冲消融的电能量于电极针裸露段向四 周传递而形成的三维消融区域。在预测三维电脉冲消融区域时，可以首先沿电极针长度方向 选取多个绘制点，然后针对每个选取的绘制点，获取约束条件下的目标电场强度等值线(在 某些应用场景中，每个选取的绘制点对应的目标电场强度等值线即为针对该绘制点预测的二 维电脉冲消融区域)，之后基于所有绘制点的目标电场强度等值线构建三维模型，从而可以得 到预测的三维电脉冲消融区域。
关于约束条件
所述约束条件包含的条件可以是单一条件，也可以是多个条件的组合。通过合理设置约 束条件，可以从电场强度数据库中快速筛选出电场强度等值线，获得电场强度等值线集合。 所述包含的条件可以为具体的数值，也可以为数值范围等。在某些应用场景中，所述约束条 件与治疗计划相关。一般来说，绘制点不同，约束条件会存在差异。
对于所述多个绘制点，响应于一次治疗，一般采用相同的约束条件。针对多次治疗，每 次治疗的约束条件可以相同，也可以不同，在约束条件包含的条件是多个条件的组合时，所 述多个条件可以相同或不同，或者其中部分条件相同。
在一实施例中，所述约束条件包括第一约束条件和第二约束条件。通过所述第一约束条 件和所述第二约束条件对所述电场强度数据库进行阶段性约束，从而可以快速有效地获得特 定的目标电场强度等值线。
关于第一约束条件
所述第一约束条件包含的条件可以是单一条件，也可以是多个条件的组合。通过合理设 置第一约束条件，可以从电场强度数据库中快速筛选出电场强度等值线，获得电场强度等值 线集合。所述包含的条件可以为具体的数值，也可以为数值范围等。在某些应用场景中，所 述第一约束条件与治疗计划相关。
在一实施例中，所述第一约束条件包括电场强度消融阈值，电场强度消融阈值可以用Eth表示。电场强度消融阈值表示施加到生物组织使生物组织细胞死亡的电场强度的阈值。每位 患者都是一个个体，每个个体都存在差异，生物组织的电场强度消融阈值不尽相同，相应的 治疗方案也不尽相同。
关于生物组织的电场强度消融阈值
电场强度消融阈值表示施加到生物组织使生物组织细胞死亡的电场强度的阈值。在一实 施例中，生物组织的电场强度消融阈值可以由多个变量共同决定，因此可以采用基于多个变 量建立的拟合函数关系来确定生物组织的电场强度消融阈值。所述多个变量包括来表示所述 拟合函数的表达式的拟合变量，拟合变量出现在拟合函数的表达式中。在一实施例中，拟合 变量可以包括电脉冲的电场强度、电脉冲中包括的脉冲串个数和消融区域的电导率比率。电 导率比率R等于完全电穿孔的电导率σmax与初始电导率σ0的比值，σ0一般等于1。由于电脉 冲的电场强度、电脉冲中包括的脉冲串个数和消融区域的电导率比率能够有效反映不同患者 的个体差异，因此，基于电脉冲的电场强度、电脉冲中包括的脉冲串个数和消融区域的电导 率比率建立的拟合函数确定的生物组织的电场强度消融阈值，能够简单且直接地建立生物组 织的电场强度消融阈值与不同患者的关联性，提高电脉冲消融区域预测的合理性和有效性。 电脉冲的电场强度可以用E表示，电脉冲中包括的脉冲串个数可以用N表示，消融区域的电 导率比率可以用R表示。电脉冲的电场强度可以由电脉冲的脉冲电压与电极针间距之比获取， 其中电脉冲的脉冲电压用U表示，电极针间距用D表示，即E＝U/D。在一优选实施例中，所 述拟合函数的表达式为：Eth(E,N,R)＝a1*E+b1*N+c1*R+d1*E*N+e1*E*R+f1*N*R+g1*E*N*R+h1，其 中，Eth表示电场强度消融阈值，E表示电脉冲的电场强度，N表示电脉冲中包括的脉冲串个数， R表示消融区域的电导率比率，a1～h1表示消融区域预测参数。在另一优选实施例中，所述拟 合函数的表达式为：Eth(U,D,N,R) ＝a2*U+b2*N+c2*D+d2*R+e2*U*N+f2*U*D+g2*U*R+h2*N*D+i2*N*R+j2*L*R+k2*N*D*R+l2*U*N*D+m2*U*N*R+n2*U*D*R+o2*N*D*R+p2，其中，Eth表示电场强度消融阈值，U表示电脉冲的脉冲电压(即电极针间的电压)，N表示电脉冲中包括的脉冲串个数，D表示电极针间距，R表示消融区域的电导率比率，a2～p2为消融区域预测参数。所述消融区域预测参数可以更新，在某些应用场 景中，所述消融区域预测参数可以实时更新。
所述消融区域预测参数可以基于历史数据采用拟合的方式确定。历史数据反映出不同数 值的电脉冲的电场强度E、电脉冲中包括的脉冲串个数N和消融区域的电导率比率R对应的 电场强度消融阈值Eth的数值，或者反映出不同数值的电脉冲的脉冲电压U、电极针间距D、 电脉冲中包括的脉冲串个数N和消融区域的电导率比率R对应的电场强度消融阈值Eth的数 值。在一实施例中，历史数据来自于实际治疗记录数据形成的治疗记录数据库，通过大量的 实际治疗记录数据对拟合函数进行拟合，确定消融区域预测参数，这种确定消融区域预测参 数的方法可以有效提高消融区域预测参数的合理性、有效性和准确性，为准确确定生物组织 的电场强度消融阈值提供可能，从而提高电脉冲消融区域预测的准确性。
为了更加准确地建立生物组织的电场强度消融阈值与不同患者的关联性，所述多个变量 还可以包括条件变量，条件变量用于选择对应的拟合函数的表达式，不出现在拟合函数的表 达式中，在采用拟合函数确定生物组织的电场强度消融阈值时，基于条件变量选择对应的拟 合函数的表达式，基于选择的对应的拟合函数的表达式和拟合变量来确定生物组织的电场强 度消融阈值。在一实施例中，条件变量包括因电极针确定而具有确定离散值的变量，条件变 量不出现在拟合函数的表达式中，可以有效降低拟合函数的表达式的复杂度，减少计算量。 在一可选实施例中，所述条件变量可以包括子脉冲宽度和电极针裸露长度，子脉冲宽度可以 用W表示，电极针裸露长度可以用L表示。其中子脉冲宽度与脉冲串内脉冲个数的乘积为固 定值，脉冲串内脉冲个数可以用n表示，固定值可以按照行业要求或行业规范设定，例如为 100，此时W*n＝100。子脉冲宽度W和电极针裸露长度L均为离散值，并且其数值随电极针确 定而确定，因此可以响应于不同的子脉冲宽度W和电极针裸露长度L，建立确定生物组织的 电场强度消融阈值的拟合函数的表达式，而子脉冲宽度W和电极针裸露长度L不出现在拟合 函数的表达式中，即子脉冲宽度W和电极针裸露长度L作为条件变量，电脉冲的电场强度E、 电脉冲中包括的脉冲串个数N和消融区域的电导率比率R作为拟合变量，基于条件变量选择 对应的拟合函数的表达式，基于拟合变量和选择的对应的拟合函数的表达式确定电场强度消 融阈值，以提高拟合函数的准确性。可以理解，对于不同的条件变量，拟合函数的表达式中 的消融区域预测参数一般来说不同。在一实施例中，历史数据反映出不同数值的子脉冲宽度 W、电极针裸露长度L、电脉冲的脉冲电压E、脉冲串个数N和消融区域的电导率比率R对应 的电场强度消融阈值Eth的数值，或者反映出不同数值的子脉冲宽度W、电极针裸露长度L、 电脉冲的脉冲电压U、电极针间距D、电脉冲中包括的脉冲串个数N和消融区域的电导率比率 R对应的电场强度消融阈值Eth的数值。可以基于历史数据采用拟合的方式确定拟合函数的表 达式中的消融区域预测参数。
例如对于特定的脉冲发生器，假设子脉冲宽度W可以为2us、5us、10us、20us、50us或 100us，电极针裸露长度L可以为1.0cm、1.5cm、2.0cm、2.5cm、3.0cm、3.5cm、4.0cm，子 脉冲宽度W与脉冲串内脉冲个数n的乘积为100。在基于多个变量建立拟合函数时，如果子 脉冲宽度W为2us和电极针裸露长度L为1.0cm，基于对应于子脉冲宽度W为2us和电极针 裸露长度L为1.0cm的实际治疗记录数据，确定拟合函数的消融区域预测参数；如果子脉冲 宽度W为2us和电极针裸露长度L为1.5cm，基于对应于子脉冲宽度W为2us和电极针裸露 长度L为1.5cm的实际治疗记录数据，确定拟合函数的消融区域预测参数；如果子脉冲宽度 W为2us和电极针裸露长度L为2.0cm，基于对应于子脉冲宽度W为2us和电极针裸露长度L 为2.0cm的实际治疗记录数据，确定拟合函数的消融区域预测参数；……；如果子脉冲宽度 W为50us和电极针裸露长度L为1.0cm，基于对应于子脉冲宽度W为50us和电极针裸露长度 L为1.0cm的实际治疗记录数据，确定拟合函数的消融区域预测参数；……。以此类推，确 定出不同的子脉冲宽度W和电极针裸露长度L对应的拟合函数的表达式的消融区域预测参数。 在采用拟合函数确定生物组织的电场强度消融阈值时，首先选择子脉冲宽度W和电极针裸露 长度L，基于所述选择的子脉冲宽度W和电极针裸露长度L，选择对应的拟合函数的表达式， 基于所述选择的拟合函数的表达式和拟合变量计算生物组织的电场强度消融阈值。可以理解，基于不同的拟合函数的表达式计算出的电场强度消融阈值的数值可以出现相同的情形。
所述历史数据可以用治疗记录表的列表形式进行记录，例如可以参考下表1所示。
表1治疗记录表(示例)
在另一实施例中，还可以采用神经网络模型来确定生物组织的电场强度消融阈值。采用 神经网络模型来确定生物组织的电场强度消融阈值一般适用于历史数据或样本数据较多的情 况。神经网络模型的输入节点可以包括决定生物组织的电场强度消融阈值的所有变量或者部 分变量，可以理解的是，一般来说，神经网络模型的输入节点越多，神经网络模型输出的预 测结果越准确。示例性地，可参考图3所示，所述神经网络模型可以包括三层神经网络，第 一层神经网络为输入层，输入层的输入节点包括决定生物组织的电场强度消融阈值的变量， 所述输入节点可以包括但不局限于子脉冲宽度W、电极针裸露长度L、电脉冲的脉冲电压U、 电极针间距D、电脉冲中包括的脉冲串个数N和消融区域的电导率比率R中的多个变量；第 二层神经网络为中间层，中间层也可以称为隐藏层，中间层可以设定为包括三个中间节点， 或者根据计算和处理能力等设定为包括四个或其他数量的中间节点；第三层神经网络为输出 层，输出层包括一个输出节点，输出节点表示神经网络模型预测的电场强度消融阈值。在一 可选实施例中，可以基于历史数据训练神经网络模型。在一实施例中，历史数据反映出不同 数值的子脉冲宽度W、电极针裸露长度L、电脉冲的脉冲电压E、脉冲串个数N和消融区域的 电导率比率R对应的电场强度消融阈值Eth的数值，或者反映出不同数值的子脉冲宽度W、电 极针裸露长度L、电脉冲的脉冲电压U、电极针间距D、电脉冲中包括的脉冲串个数N和消融 区域的电导率比率R对应的电场强度消融阈值Eth的数值。例如以若干组临床消融数据作为训 练神经网络模型的原始数据,每一组临床消融数据可以包括：治疗的电脉冲的电场强度(可以 通过两电极针间电位差与两电极针间距的比值来获取，即电脉冲的电场强度E＝电脉冲的脉冲 电压U/电极针间距D)、电极针间距D、电极针裸露长度L、电导率比率R(电导率比率R＝完 全电穿孔的电导率σmax/初始电导率σ0，σ0一般等于1)和电场强度消融阈值Eth。关于神经 网络模型的训练方法可以参考现有技术，在此不再赘述。
在采用拟合函数或神经网络模型确定生物组织的电场强度消融阈值时，基于治疗计划中 的相应变量的数值确定生物组织的电场强度消融阈值；基于所述确定的生物组织的电场强度 消融阈值，从电场强度数据库中获取电场强度等值线集合。
关于电场强度数据库
所述电场强度数据库存放电场强度。可以基于每个预设的绘制点单独管理电场强度数据 库。在某些应用场景中，所述绘制点选取模块是从预设的绘制点中选取所述多个绘制点。构 建电场强度数据库时，可以针对每个预设的绘制点构建其对应的电场强度数据 库，所有的绘制点对应的电场强度数据库组合成电场强度数据库。所述系统还可以包括电场强度数据库构 建模块，用于构建所述电场强度数据库。电场强度数据库可以更新。在一实施例中，所述电 场强度数据库的构建包括：针对每个预设的绘制点，基于约束方程，以参数组合模型中的参 数组合为单位获取电场强度分布参数组；基于所述获取的电场强度分布参数组，构建所述电 场强度数据库。参数组合模型表示由多个构建参数的预设值组合而成的参数组合所构成的模 型，即参数组合是从每个构建参数中各选取一个预设值组合而成，参数组合模型是由所有参 数组合构建而成的模型。不同参数组合中的各构建参数预设值不完全相同，参数组合的数量 由构建参数的个数以及各个构建参数的预设值个数决定。举例来说，假设构建参数的个数为 三个，三个构建参数的预设值个数分别为2、3和4，那么参数组合的数量不超过24个。针 对每一个参数组合，结合约束方程进行数值计算，得到对应该参数组合的场强分布参数组， 场强分布参数组中的每个场强分布参数代表一个场强等值线。针对每个预设的绘制点的参数 组合模型中的所有参数组合对应的场强分布参数组一起构成每个预设的绘制点对应的电场强 度数据库，所有预设的绘制点对应的电场强度数据库构成电场强度数据库。
在一可选实施例中，构建所述参数组合模型的构建参数包含：消融区域的第一电导率比 率、第一电脉冲的脉冲电压、第一电极针裸露长度和电极针间距。所述构建所述参数组合模 型的构建参数可以是预设值。例如对于特定的脉冲发生器，电极针间距的预设值可以包括 0.5cm、1.0cm、1.5cm、2.0cm、2.5cm、3.0cm、3.5cm、4.0cm，第一电极针裸露长度的预设 值可以包括1.0cm、1.5cm、2.0cm、2.5cm、3.0cm、3.5cm、4.0cm，第一电脉冲的脉冲电压的预设值可以包括500V、1000V、1500V、2000V、2500V、3000V，消融区域的第一电导率比 率的预设值可以包括1.0、1.5、2.0、2.5、3.0、3.5、4.0，在此示例中，参数组合的数量 不超过2352个。
在一可选实施例中，每个[消融区域的第一电导率比率、第一电脉冲的脉冲电压、第一电 极针裸露长度、电极针间距]参数组合对应一个场强分布，可以记为E＝f(xi,yj,zk)，其中，xi、 yj、zk为xyz三维坐标系中按照一定规律分布的离散点，zk对应各个预设的绘制点，E为离散 点(xi,yj,zk)对应的电场强度。两个电极针与z轴平行，且两电极针在xy平面内关于z轴 对称。坐标原点位于电极针裸露段中心连线的中点处。针对每个预设的绘制点zk，所有[消融 区域的第一电导率比率、第一电脉冲的脉冲电压、第一电极针裸露长度、电极针间距]参数组 合的场强分布的集合即为预设的绘制点zk对应的电场强度数据库，而所有预设的绘制点对应 的电场强度数据库的集合构成电场强度数据库。xyz三维坐标系可参考图2所示，每个选取 的绘制点对应的电场强度等值线即z＝zk的截面上的电场强度等值线可参考图4所示。
关于约束方程
约束方程可以用来获取每个参数组合模型对应的场强分布参数组。
在一优选实施例中，所述约束方程包括电场方程、电导率方程和边界条件方程。所述边 界条件方程为：
u1＝U u2＝0，其中，
表示参数组合模型的边 界上的每个点的法向量，参数组合模型的边界表面为一个曲面，边界表面上的每一个点的法 向量都不同，但都与该点所在的平面垂直，J表示电流密度，
表示电流密度矢量，u1和u2分别表示两个电极针上的电位，U表示电脉冲的脉冲电压，边界条件方程用来表示电流不会 流出参数组合模型的边界表面。电场方程为：
其中，J表示电流密度，
表示电流密度矢量，σ表示电导率，E表示电脉冲的电场强度，
表示电脉冲的电场强度矢量，
表示电流密度矢量的散度，u表示电脉冲的电位，▽u表 示u的梯度。电导率方程为：σ(E)＝σ0+(σmax-σ0)*exp[-A*exp(-B*E)]，其中，σ0为不发生 电击穿的组织电导率(也可以称为初始电导率)，σmax为生物组织细胞发生完全透化的组织电 导率(也可以称为完全电击穿、透化电导率)，A和B为决定组织电导率曲线的位置和增长率 的系数。
关于电场强度等值线
在利用电脉冲消融肿瘤细胞的治疗中，对生物组织施加的电脉冲的脉冲电压U、电极针 间距D、电场强度消融阈值Eth是判断电脉冲消融区域的大小和范围的关键因素。由于生物组 织的电导率σ是电脉冲的电场强度E的函数σ(E)，目前应用较多的电导率方程为：
σ(E)＝σ0+(σmax-σ0)*exp[-A*exp(-B*E)]
其中，σ0为不发生电击穿的组织电导率(也可以称为初始电导率)，σmax为生物组织细 胞发生完全透化的组织电导率(也可以称为完全电击穿、透化电导率)，A和B为决定组织电 导率曲线的位置和增长率的系数。
由于σ(E)不是常量，因此无法获得电场强度的解析解，需要通过数值计算需求的数值解， 图4示出了数值计算得出的生物组织内电极针垂直方向上的中间截面(即z＝0)上的场强等 值线，其中一条(可参考图5a的情形)或两条(可参考图5b的情形)封闭线代表一条电场 强度等值线，不同的封闭线代表电场强度不同。场强等值线的形状呈多种形态，可参考图4 所示，因此需要对场强等值线进行重建。在一实施例中，所述电场强度等值线集合中的电场 强度等值线采用卡西尼曲线或样条拟合曲线表示。
关于卡西尼(Cassini)曲线匹配
假定卡西尼曲线关于x轴和y轴对称，中心位于坐标原点，且两焦点位于x轴上，分别 为(-a,0)和(a,0)，则卡西尼曲线的方程为：
[(x-a)2+y2]*[(x+a)2+y2]＝b4
其中，(x,y)为卡西尼曲线上点的坐标，卡西尼曲线形状取决于形状参数a、b，卡西尼 曲线的两焦点对应两电极针的位置。
卡西尼曲线与x轴和y轴的交点有两种情形：
i)交点为(M,0)、(0,N)，即卡西尼曲线与x轴和y轴均有两个交点，可参考图5a所示；
ii)交点为(M,0)、(N,0)，即卡西尼曲线与x轴有四个交点，与y轴无交点，可参考图5b所示。
形状参数a、b的获取方式如下：
对于第i)种情形：a＝sqrt[(M2-N2)/2] b＝sqrt[(M2+N2)/2]；
对于第ii)种情形：a＝sqrt[(M2+N2)/2] b＝sqrt[(M2-N2)/2]。
在不引起歧义的情况下，所述场强分布参数可以用(M，N)来表示。
针对每个预设的绘制点，每个[消融区域的第一电导率比率、第一电脉冲的脉冲电压、第 一电极针裸露长度、电极针间距]参数组合都会有相应的场强分布参数(M、N)；那么针对每 个预设的绘制点，所有[消融区域的第一电导率比率、第一电脉冲的脉冲电压、第一电极针裸 露长度、电极针间距]参数组合相应的场强分布参数(M、N)就构成了每个预设的绘制点对应 的电场强度等值线，可参考图4所示。
关于样条拟合曲线
针对每个预设的绘制点，电场强度等值线关于x轴和y轴对称，因此只需画出第一象限 的曲线就可得到整个电场强度等值线。首先，以一定间距的x和y离散值提取出电场强度分 布，可以根据具体模型合理分配间距，由提取的x和y离散值即可复现对应截面上的电场强 度分布，这里的间距表示的是相邻点之间的距离，在实际操作时可以以规律间距的x和y离 散值提取出电场强度分布。合理分配间距的所述模型表示的是数值计算的几何模型。计算方 法包括：
判断电场强度等值线的连通情形，包括：检测y＝0的射线上的所求电场强度在第一象限 的点的数量，基于点的数量判断电场强度等值线的连通情形，若点的数量为1，如图6a所示， 则属于前述第i)种情形；若点的数量为2，如图6b所示，则属于前述第ii)种情形；
获取电场强度等值线的取值范围，包括：由线性插值，基于所述连通情形，计算所求电 场强度等值线在第一象限与x轴的交点(即y＝0的点)，则将x的取值范围表示为x∈[M,N]， 其中，M和N表示所求电场强度等值线在第一象限与x轴的交点的x坐标，如果交点为1个， 即点(N,0)，则M＝0，N≠0，可参考图6a所示；如果交点为两个，即点(M,0)和点(N,0)，则M ≠0，N≠0，可参考图6b所示；
求取电场强度等值线上的点，包括：在取值范围x∈[M,N]内，以步长间隔求取所求电场 强度等值线在第一象限对应的y坐标，则可获得一组电场强度等值线上的点集(xi,yi)，i代 表点集中点的序列号，例如i＝1、2、3、…，可参考图7a所示；步长间隔的大小可以视情况 选择，例如步长间隔可以为0.1mm；
获取样条拟合曲线，包括：根据获得的一组电场强度等值线上的点集(xi,yi)进行样条 拟合，得出第一象限的电场强度等值线，并对第一象限的电场强度等值线进行关于x轴和y 轴的对称操作，形成封闭的曲线，从而获得样条拟合曲线，可参考图7b所示。
在不引起歧义的情况下，所述场强分布参数可以用(M，N)来表示。
基于上述方法，可以获得针对每个预设的绘制点的电场强度等值线。
关于第二约束条件
所述第二约束条件包含的条件可以是单一条件，也可以是多个条件的组合。所述包含的 条件可以为具体的数值，也可以为数值范围等。在某些应用场景中，所述第二约束条件与治 疗计划相关。
在一实施例中，所述第二约束条件包括：第二脉冲电压、消融区域的第二电导率比率、 第二电极针裸露长度和电极针坐标。
根据第二脉冲电压、消融区域的第二电导率比率、第二电极针裸露长度和电极针坐标， 从所述获取的电场强度等值线集合中获取目标电场强度等值线，所述获取的目标电场强度等 值线内的区域作为预测电脉冲消融区域。
关于消融区域的第二电导率比率
第二电导率比率R等于完全电穿孔的电导率σmax与初始电导率σ0的比值，σ0一般等于1。
在一实施例中，第二电导率比率R可以通过释放预脉冲的方法获取，包括：对每对电极 针施加预脉冲至电流稳定，基于一表征参数在所述每对电极针的稳定状态和在所述每对电极 针的初始状态的数值比，获取所述消融区域的第二电导率比率；
或者，
所述消融区域的第二电导率比率通过释放预脉冲的方法获取，包括：对每对电极针施加 两次预脉冲至电流稳定，基于一表征参数在所述每对电极针的两次稳定状态的数值比，获取 所述消融区域的第二电导率比率。
其中，所述表征参数包括电流或电阻。
在一实施例中，所述获取所述消融区域的第二电导率比率包括通过查询列表的方式获取 所述消融区域的第二电导率比率。进一步地，在一实施例中，所述列表表征在不同的电极针 裸露长度时的电极针间距、所述表征参数的数值比和消融区域的电导率比率的关系。
下面结合具体示例对消融区域的第二电导率比率进行说明。
在治疗过程中，首先会施加预设电场强度(假定为500V/cm)的电脉冲，电流随治疗时 间缓慢上升，响应于预设电场强度的电脉冲，脉冲串与电流的关系可以参考图8所示，图中，I0为初始电流，Is500为预设电场强度为500V/cm时的稳定电流，不同患者产生的初始电流I0和稳定电流Is500均不同。由于患者个体的差异，难以确定初始电导率σ0，因此在一实施例中， 采用电流比S来抵消初始电导率σ0产生的影响，电流比S可以采用施加同一预设电场强度时 的稳定电流Is和初始电流I0的比值来计算。在该实施例中，第二电导率比率的获取方法包括：
构建关于电流的数据库，包括：根据预设参数组合建立计算模型，通过仿真软件(例如 Comsol仿真软件)再结合电场方程、电导率方程以及边界条件构建电流数据库；其中预设参 数包括：
预设电极针间距，取值例如可以为：0.5cm，1.0cm，1.5cm，2.0cm，2.5cm，3.0cm，3.5cm， 4.0cm；
预设电极针裸露长度，取值例如可以为：1.0cm、1.5cm、2.0cm、2.5cm、3.0cm、3.5cm、 4.0cm；
预设电场强度：500V/cm；
预设电导率比率：1.0、1.5、2.0、2.5、3.0、3.5、4.0；
电流数据库表征了在预设电场强度的条件下，预设电极针间距、预设电极针裸露长度和 预设电导率比率与电流的关系。
在构建的电流数据库的基础上，计算电流比的数值集合。所述电流比表示在预设电场强 度的条件下，预设电极针间距、预设电极针裸露长度和预设电导率比率组合下计算所得电流 与和预设电导率比率为1.0时计算所得电流的比值。
确定治疗数据，治疗数据包括：电极针间距D、电极针裸露长度L和电流比S。治疗数据 与治疗计划相关。其中电流比S的获取包括：对每对电极针施加预设电场强度的电脉冲至电 流稳定，假设预设电场强度(通过电极针间的电压U与电极针间距D的比值计算得到)为 500V/cm，得到初始电流I0和稳定电流Is500，计算稳定电流Is500和初始电流I0的比值，根据电 流比S＝Is500/I0，得到电流比S。
获取第二电导率比率，包括：根据治疗数据在电流比的数值集合中查找对应的电导率比 率。可以参考下表所示，表中给出了在某电极针裸露长度下，预设电极针间距和电流比与预 设电导率比率的关系。在使用该表查找第二电导率比率时，将治疗数据包括的电极针裸露长 度L作为预设电极针裸露长度找到电极针裸露长度L对应的关系表，将治疗数据包括的电极 针间距D作为预设电极针间距，并根据治疗数据包括的电极针间距D和电流比在所述对应的 关系表中查找对应的电导率比率，作为第二电导率比率。
| D<sub>1</sub> | D<sub>2</sub> | … | D<sub>n</sub> | |
| R<sub>sig1</sub> | S<sub>11</sub> | S<sub>12</sub> | … | S<sub>1n</sub> |
| R<sub>sig2</sub> | S<sub>21</sub> | S<sub>22</sub> | … | S<sub>2n</sub> |
| … | … | … | … | … |
| R<sub>sigm</sub> | S<sub>m1</sub> | S<sub>m2</sub> | … | S<sub>mn</sub> |
表2：某预设电极针裸露长度时预设电极针间距、预设电导率比率与电流比的关系
可替换地，电流比S可以采用施加不同预设电场强度时的稳定电流的比值来计算，且电 流比S大于1。例如，分别施加预设电场强度为500V/cm和1000V/cm的电脉冲，可参考图9 所示，各自的稳定电流分别记为Is500和Is1000，则电流比S可以用Is1000/Is500表示。
可替换地，电流比S可以采用施加同一预设电场强度时的稳定电阻Rs和初始电阻R0的比 值来计算。例如，施加预设电场强度为500V/cm的电脉冲时，初始电阻记为R0，稳定电阻记 为Rs500，电流比可以用Rs500/R0表示。
另可替换地，电流比S可以采用施加不同预设电场强度时的稳定电阻的比值来计算，且 电流比S大于1。例如，分别施加预设电场强度为500V/cm和1000V/cm的电脉冲，各自的稳 定电阻分别记为Rs500和Rs1000，则电流比S可以用Rs1000/Rs500表示。
关于坐标变换
所述系统还可以包括坐标变换模块，用于响应于布针位置不在同一水平线上，对所述目 标电场强度等值线进行坐标变换。可参考图10所示。
在一实施例中，响应于布针位置(即两电极针的位置)不在同一水平线上时，可通过平 移变换矩阵将目标电场强度等值线变换至布针位置。假设布针位置的坐标为(x1,y1)和(x2,y2)， 两电极针的中心记为(x0,y0)，目标电场强度等值线的旋转角度记为θ。
则，x0＝(x1+x2)/2，y0＝(y1+y2)/2，θ＝arctan[(y2-y1)/(x2-x1)；
卡西尼曲线方程中的x和y可由如下公式计算：
x＝cosθ*(x’-x0)+sinθ*(y’-y0)
y＝-sinθ*(x’-x0)+cosθ*(y’-y0)
在制定治疗计划的过程中，根据电场强度消融阈值Eth、电极针裸露长度L、电导率比率 R、电脉冲的脉冲电压U、电极针坐标等参数，即可计算出治疗计划对应的目标消融区域形状。
关于区域叠加
本发明实施例的预测电脉冲消融区域的系统还可以包括区域叠加模块，用于针对每个选 取的绘制点，响应于多次治疗，叠加针对每次治疗所获取的目标电场强度等值线内的区域， 作为所述多次治疗的预测电脉冲消融区域。可参考图11所示。
在一实施例中，可以增加电极针的布针、电脉冲的脉冲电压、治疗顺序等，将每次治疗 计算得出的目标消融区域叠加，得到最终的目标消融区域。治疗顺序例如对于包括多于两电 极针的情形，第一次施加电脉冲使用针号1、2，第二次施加电脉冲使用针号2、3，……，如 此按序进行治疗。
关于三维模型的构建
在一实施例中，所述基于所述每个选取的绘制点对应的目标电场强度等值线，构建沿电 极针长度方向的三维模型，包括：基于所述每个选取的绘制点对应的目标电场强度等值线， 以网格的形式构建沿电极针长度方向的三维模型。进一步地，所述网格在每对相邻的绘制点 对应的目标电场强度等值线之间构建。可以理解，在构建三维模型时，可以对三维模型的外 表面进行平滑处理，以使构建的三维模型更能准确地反映电脉冲在生物组织上的消融效果。
在一实施例中，所述基于所述每个选取的绘制点对应的目标电场强度等值线，以网格的 形式构建沿电极针长度方向的三维模型，可以包括：在每个选取的绘制点对应的目标电场强 度等值线分别选取m个点(m为自然数)，构建等值线点集，其中选取的m个点的示例可以参 考图12a所示(图12a中示出了m＝2n的情形，可参考后续的描述)，对于xyz三维坐标系， 图12a所示的平面为xy平面或与xy平面平行的平面；构建的等值线点集可以参考图12b所 示，图12b中的Qi和Qi+1表示所述相邻的绘制点对应的目标电场强度等值线，图12b所示的 立体图为xyz三维坐标系下的立体图，后续的图12c和图12d与图12b相同。基于构建的等 值线点集，在相邻等值线的点集间构建网格，可参考图12c所示；形成的最终网格内的区域 即为三维消融区域，可参考图12d所示。所述m个点可以等间隔选取，也可以任意选取。
在一可选实施例中，所述在每个选取的绘制点对应的目标电场强度等值线分别选取m个 点，可以采用在相邻的绘制点对应的目标电场强度等值线Qi和Qi+1分别取间隔相等的m个点， 在图12a所示的情形中，所述m个点为间隔相等的2n个点，n为自然数。更可选地，所述m 个点为从等值线中心(在xyz三维坐标系中即x＝0，y＝0)向外延伸的m条射线(在间隔相等 时，m条射线中相邻两条射线的夹角相等)分别与等值线Qi和Qi+1的交点，在图12a所示的情 形中，所述2n个点为经过等值线中心的n条直线(在间隔相等时，n条直线中相邻两条直线 的夹角相等)分别与等值线Qi和Qi+1的交点。这里取间隔相等的m个点或2n个点，主要是为 了计算简单，并不排除任意取m个点或者任意取2n个点或2n+1个点的情形。
为了提高预测的准确性，所述在每个选取的绘制点对应的目标电场强度等值线分别选取 m个点，包括：由m个从等值线中心线(在xyz三维坐标系中即为z轴)向外延伸的平面与 所述每个选取的绘制点对应的目标电场强度等值线的交点。在图12a所示的情形中，所述在 每个选取的绘制点对应的目标电场强度等值线分别选取2n个点，包括：由n个经过等值线中 心线(在xyz三维坐标系中即为z轴)的平面与所述每个选取的绘制点对应的目标电场强度 等值线的交点。相应地，所述网格在每对相邻的绘制点对应的目标电场强度等值线之间构建， 包括：基于经过等值线中心线的平面的位置顺序，在每对相邻的绘制点对应的目标电场强度 等值线之间构建网格。举例来说，假设第一平面和第二平面为经过所述等值线中心线的两个 相邻平面，等值线Qi与第一平面的交点为A1，与第二平面的交点为A2；等值线Qi+1与第一平 面的交点为B1，与第二平面的交点为B2，在等值线Qi和Qi+1之间构建的网格可以为A1、A2、B1和B2构建而成的网格。可以理解，网格的构建方式不局限于此，在某些应用场景中，还可以 在具有一定间隔的一对绘制点对应的目标电场强度等值线之间构建网格，所述的间隔可以是 间隔一个或多个绘制点。
在一可选实施例中，可以采用样条拟合的方法构建三维模型，具体来说，可以先构建一 个象限的三维模型，然后对所述一个象限的三维模型关于x轴、y轴和z轴的对称操作，形 成封闭的三维模型。
下面对本发明实施例的具体应用场景进行示例性说明，可参考图13所示。
在100，构建电场强度数据库。可以通过电场强度数据库构建模块实现。
例如，针对每个预设的绘制点，根据电极针间距D的预设值、电极针裸露长度L的预设 值、电脉冲的脉冲电压U的预设值以及电导率比率σmax/σ0(记为R)的预设值的不同组合，建立由D、L、U和R的各种参数预设值组合而成的参数组合模型；针对参数组合模型中的每个参数组合，结合电场方程、电导率方程以及边界条件方程进行数值计算，得到与每个参数组合对应的场强分布参数组，每一个场强分布参数代表一条电场强度等值线；参数组合模型 中的所有参数组合对应的场强分布参数组一起形成每个预设的绘制点对应的电场强度数据 库，从而可以获得所有预设的绘制点对应的电场强度数据库。
在200，选取绘制点。可以通过绘制点选取模块实现。
例如基于电极针裸露长度范围，沿电极针长度方向选取多个绘制点。所述选取的绘制点 一般来自所述预设的绘制点。
在300，针对每个选取的绘制点，获取与约束条件对应的目标电场强度等值线。可以通 过二维区域获取模块实现。
所述约束条件可以包括第一约束条件和第二约束条件。例如在集合获取模块中，以电场 强度消融阈值Eth作为第一约束条件，将电场强度消融阈值Eth的值固定。根据固定的电场强 度消融阈值Eth计算出在某电导率比率R(σmax/σ0)和电极针裸露长度L组合时，不同电脉冲 的电场强度E(其中，E＝U/D)和电极针间距D的场强分布参数，可以采用表3所示的表格形 式进行管理；同理，计算出各个不同的电导率比率R和电极针裸露长度L组合时，各个不同 电脉冲的电场强度E和电极针间距D的场强分布参数；所有对应该固定的电场强度消融阈值 Eth的场强分布参数，构成了针对每个选取的绘制点的与第一约束条件对应的电场强度等值线 集合。表格中M和N的含义可参考前述的描述。
| E<sub>1</sub> | E<sub>2</sub> | … | E<sub>n</sub> | |
| D<sub>1</sub> | M<sub>11</sub>,N<sub>11</sub> | M<sub>12</sub>,N<sub>12</sub> | … | M<sub>1n</sub>,N<sub>1n</sub> |
| D<sub>2</sub> | M<sub>21</sub>,N<sub>21</sub> | M<sub>22</sub>,N<sub>22</sub> | … | M<sub>2n</sub>,N<sub>2n</sub> |
| … | … | … | … | … |
| D<sub>m</sub> | M<sub>m1</sub>,N<sub>m1</sub> | M<sub>m2</sub>,N<sub>m2</sub> | … | M<sub>mn</sub>,N<sub>mn</sub> |
表3电场强度消融阈值Eth的值固定，某电导率比率R和电极针裸露长度L组合时对应的场强分布参数表
在目标获取模块中，基于第二约束条件，获取每个选取的绘制点对应的目标电场强度等 值线。例如，将电极针裸露长度L的值固定，对所述与第一约束条件对应的电场强度等值线 集合中与第一约束条件对应的电场强度等值线集合，应用拉格朗日插值公式，计算出所述固 定的电极针裸露长度L对应的电脉冲的电场强度E和电极针间距D的场强分布参数，此时不 同电导率比率R对应的场强分布参数表可以简化为表4；将电导率比率R固定，对所述固定 的电极针裸露长度L对应的电脉冲的电场强度E和电极针间距D的场强分布参数，应用拉格 朗日插值公式，得到目标场强分布参数，即可得到目标电场强度等值线。
| R<sub>1</sub> | R<sub>2</sub> | … | R<sub>n</sub> |
| M<sub>1</sub>,N<sub>1</sub> | M<sub>2</sub>,N<sub>2</sub> | … | M<sub>n</sub>,N<sub>n</sub> |
表4电场强度消融阈值Eth和电极针裸露长度L的值固定，不同电导率比率R对应的场强分布参数表
在400，针对每个选取的绘制点，绘制目标电场强度等值线。可以通过绘制模块实现。
由于两电极针可能不在同一水平线(同一水平线对于前述的坐标系而言可以理解为x轴 或与x轴平行的线)上，因此会出现以下两种情形。
第一种情况：两电极针在同一水平线上,此时可以根据目标场强分布参数绘制目标电场强 度等值线。
第二种情况：两电极针不在同一水平线上，此时可以通过平移变换矩阵将目标电场强度 等值线变换至布针位置。可以通过坐标变换模块实现。假设布针坐标为(x1,y1)和(x2,y2)，两 电极针的中心记为(x0,y0)，目标电场强度等值线(或者目标场强分布参数)的旋转角度记为 θ。
则，x0＝(x1+x2)/2，y0＝(y1+y2)/2，θ＝arctan[(y2-y1)/(x2-x1)；
卡西尼曲线方程中的x，y可由如下公式计算：
x＝cosθ*(x’-x0)+sinθ*(y’-y0)
y＝-sinθ*(x’-x0)+cosθ*(y’-y0)
在制定治疗计划的过程中，根据消融阈值Eth、电极针裸露长度L、电导率比率R、治疗 场强E、电极针坐标等参数，即可以计算某治疗方案下的目标电场强度等值线，在治疗前了 解二维目标消融区域的形状。
针对每个选取的绘制点，在治疗方案包括多次治疗的情况下，叠加针对每次治疗所获取 的目标电场强度等值线内的区域，作为针对每个选取的绘制点的目标电场强度等值线。可以 通过区域叠加模块实现。多次治疗可以是增加电极针的布针、电脉冲的脉冲电压、治疗顺序 (例如第一次施加电脉冲使用针号1、2，第二次施加电脉冲使用针号2、3……)等。
另外，为了方便查看，可以显示绘制出的目标电场强度等值线。可以通过显示模块实现。
在500，构建三维模型。可以通过三维区域构建模块实现。
例如，在电极针长度方向上绘制不同高度z(z∈[-L,L]，L为电极针裸露长度)时的二 维消融区域等值线组Q，之后分别在某等值线Qi和等值线Qi+1(Qi与Qi+1相邻)上取间隔相等 的2n个点，举例来说可以计算经过等值线中心(x＝0,y＝0)的n条直线(相邻两条直线的夹 角相等)与某等值线Qi的2n个交点，如图12a所示，计算经过等值线中心(x＝0,y＝0)的n条直线(相邻两条直线的夹角相等)与等值线Qi+1的2n个交点，如此类推，即可获取二维消融区域等值线组Q对应的等值线点集，如图12b所示；然后，在相邻等值线Qi和Qi+1的点集 间构建网格，如图12c所示，并形成最终的网格，即三维消融区域，如图12d所示。
本领域技术人员可以理解，在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例 如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机 可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或 其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执 行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存 储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置 或设备而使用的装置。
计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电 连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(RAM)，只读存储器(ROM)，可擦除可编辑只读存储器(EPROM或闪速存储器)，光纤装置，以及便携式光盘 只读存储器(CDROM)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其 他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必 要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施 方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件 来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术 中的任一项或它们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻 辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(PGA)，现场可编程门 阵列(FPGA)等。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、 或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包 含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指 的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个 或多个实施例或示例中以合适的方式结合。
以上，对本发明的实施方式进行了说明。但是，本发明不限定于上述实施方式。凡在本 发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范 围之内。
Claims (19)
1.一种预测电脉冲消融区域的系统，其特征在于，包括：
绘制点选取模块，用于基于电极针裸露长度范围，沿电极针长度方向选取多个绘制点；
二维区域获取模块，用于针对每个选取的绘制点，基于约束条件，从每个选取的绘制点对应的电场强度数据库中获取所述每个选取的绘制点对应的目标电场强度等值线；以及，
三维区域构建模块，用于基于所述每个选取的绘制点对应的目标电场强度等值线，构建沿电极针长度方向的三维模型，所述构建的三维模型内的区域作为预测的电脉冲消融区域；
其中，所述约束条件包括第一约束条件和第二约束条件；
所述二维区域获取模块包括：
数据库构建模块，用于构建电场强度数据库；
集合获取模块，用于针对每个选取的绘制点，基于第一约束条件，从所述每个选取的绘制点对应的电场强度数据库中获取所述每个选取的绘制点对应的电场强度等值线集合；所述第一约束条件包括电场强度消融阈值；以及，
目标获取模块，用于针对每个选取的绘制点，基于第二约束条件，从所述每个选取的绘制点对应的电场强度等值线集合中获取所述每个选取的绘制点对应的目标电场强度等值线。
2.根据权利要求1所述的系统，其特征在于，所述电场强度消融阈值通过拟合函数或神经网络模型确定；
其中，所述拟合函数包括用于表示所述拟合函数的表达式的拟合变量和用于选择拟合函数的表达式的条件变量。
3.根据权利要求2所述的系统，其特征在于，所述拟合函数的条件变量包括子脉冲宽度和电极针裸露长度；
和/或，
所述拟合函数的拟合变量包括电脉冲的电场强度、电脉冲中包括的脉冲串个数和消融区域的电导率比率，或者，所述拟合变量包括电脉冲的脉冲电压、电极针间距、电脉冲中包括的脉冲串个数和消融区域的电导率比率。
4.根据权利要求3所述的系统，其特征在于，响应于所述拟合变量包括电脉冲的电场强度、电脉冲中包括的脉冲串个数和消融区域的电导率比率，所述拟合函数的表达式为：
Eth＝a1*E+b1*N+c1*R+d1*E*N+e1*E*R+f1*N*R+g1*E*N*R+h1，
其中，Eth表示电场强度消融阈值，E表示电脉冲的电场强度，N表示电脉冲中包括的脉冲串个数，R表示消融区域的电导率比率，a1～h1表示消融区域预测参数；
或者，
响应于所述拟合变量包括电脉冲的脉冲电压、电极针间距、电脉冲中包括的脉冲串个数和消融区域的电导率比率，所述拟合函数的表达式为：
Eth＝a2*U+b2*N+c2*D+d2*R+e2*U*N+f2*U*D+g2*U*R+h2*N*D+i2*N*R+j2*D*R+k2*N*D*R+l2*U*N*D+m2*U*N*R+n2*U*D*R+o2*N*D*R*U+p2，
其中，Eth表示电场强度消融阈值，U表示电脉冲的脉冲电压，D表示电极针间距，N表示电脉冲中包括的脉冲串个数，R表示消融区域的电导率比率，a2～p2为消融区域预测参数。
5.根据权利要求4所述的系统，其特征在于，所述消融区域预测参数基于历史数据采用拟合的方式确定，所述约束条件与治疗计划相关。
6.根据权利要求1所述的系统，其特征在于，所述选取多个绘制点是从预设的绘制点中选取多个绘制点，所述预设的绘制点对应的电场强度数据库的构建包括：
基于约束方程，以参数组合模型中的参数组合为单位获取电场强度分布参数组；
基于所述获取的电场强度分布参数组和每个预设的绘制点，构建所述每个预设的绘制点对应的电场强度数据库。
7.根据权利要求6所述的系统，其特征在于，所述约束方程包括电场方程、电导率方程和边界条件方程。
8.根据权利要求6所述的系统，其特征在于，构建所述参数组合模型的构建参数包含：消融区域的第一电导率比率、第一电脉冲的脉冲电压、第一电极针裸露长度和电极针间距。
9.根据权利要求1所述的系统，其特征在于，所述电场强度等值线集合中的电场强度等值线采用卡西尼曲线或样条拟合曲线表示。
10.根据权利要求1所述的系统，其特征在于，所述二维区域获取模块还包括：
坐标变换模块，用于响应于布针位置不在同一水平线上，对所述目标电场强度等值线进行坐标变换。
11.根据权利要求1所述的系统，其特征在于，所述二维区域获取模块还包括：
区域叠加模块，用于响应于多次治疗，叠加针对每次治疗所获取的目标电场强度等值线内的区域，作为所述多次治疗的预测电脉冲消融区域。
12.根据权利要求1所述的系统，其特征在于，所述第二约束条件包括第二电脉冲的脉冲电压、消融区域的第二电导率比率、第二电极针裸露长度和电极针坐标。
13.根据权利要求12所述的系统，其特征在于，所述消融区域的第二电导率比率通过释放预脉冲的方法获取。
14.根据权利要求13所述的系统，其特征在于，所述释放预脉冲的方法，包括：对每对电极针施加预脉冲至电流稳定，基于一表征参数在所述每对电极针的稳定状态和在所述每对电极针的初始状态的数值比，获取所述消融区域的第二电导率比率；
或者，包括：对每对电极针施加两次预脉冲至电流稳定，基于一表征参数在所述每对电极针的两次稳定状态的数值比，获取所述消融区域的第二电导率比率。
15.根据权利要求14所述的系统，其特征在于，所述表征参数包括电流或电阻。
16.根据权利要求14所述的系统，其特征在于，所述基于一表征参数在所述每对电极针的两次稳定状态的数值比，获取所述消融区域的第二电导率比率，包括：基于一表征参数在所述每对电极针的两次稳定状态的数值比，通过查询列表的方式获取所述消融区域的第二电导率比率；
所述列表表征在不同的电极针裸露长度时的电极针间距、所述表征参数的数值比和消融区域的电导率比率的关系。
17.根据权利要求1至16中任一项所述的系统，其特征在于，所述基于所述每个选取的绘制点对应的目标电场强度等值线，构建沿电极针长度方向的三维模型，包括：基于所述每个选取的绘制点对应的目标电场强度等值线，以网格的形式构建沿电极针长度方向的三维模型。
18.根据权利要求17所述的系统，其特征在于，所述网格在每对相邻的绘制点对应的目标电场强度等值线之间构建。
19.根据权利要求1至16中任一项所述的系统，其特征在于，选取[-L,L]作为所述电极针裸露长度范围，其中，L表示电极针裸露长度。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010301885.4A CN111529051B (zh) | 2020-04-16 | 2020-04-16 | 一种预测电脉冲消融区域的系统 |
| PCT/CN2020/108971 WO2021208306A1 (zh) | 2020-04-16 | 2020-08-13 | 一种预测电脉冲消融区域的系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010301885.4A CN111529051B (zh) | 2020-04-16 | 2020-04-16 | 一种预测电脉冲消融区域的系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111529051A CN111529051A (zh) | 2020-08-14 |
| CN111529051B true CN111529051B (zh) | 2021-05-07 |
Family
ID=71952229
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010301885.4A Active CN111529051B (zh) | 2020-04-16 | 2020-04-16 | 一种预测电脉冲消融区域的系统 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN111529051B (zh) |
| WO (1) | WO2021208306A1 (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111529052B (zh) * | 2020-04-16 | 2021-05-07 | 上海睿刀医疗科技有限公司 | 一种预测电脉冲消融区域的系统 |
| CN112489741B (zh) * | 2020-11-13 | 2021-11-05 | 上海睿刀医疗科技有限公司 | 一种脉冲电场消融参数优化系统 |
| CN112315578B (zh) * | 2020-11-20 | 2021-06-15 | 上海睿刀医疗科技有限公司 | 确定电极针布针组合的装置、方法及电极针布针优化系统 |
| CN112315579B (zh) * | 2020-11-25 | 2021-07-02 | 上海睿刀医疗科技有限公司 | 一种基于病灶区域的电极针布针装置及方法 |
| CN112603536A (zh) * | 2020-12-29 | 2021-04-06 | 北京华科恒生医疗科技有限公司 | 三维模型中生成电极热凝参数的方法和系统 |
| CN112998844B (zh) * | 2021-02-05 | 2021-12-21 | 上海睿刀医疗科技有限公司 | 一种球囊电极导管 |
| CN112998845B (zh) * | 2021-02-08 | 2021-12-21 | 上海睿刀医疗科技有限公司 | 一种多电极消融针以及多电极消融针电极间距的确定方法 |
| CN113380394B (zh) * | 2021-06-18 | 2022-04-12 | 上海睿刀医疗科技有限公司 | 确定电极针消融边界的方法、装置、电子设备及存储介质 |
| CN113705807B (zh) * | 2021-08-26 | 2022-06-10 | 上海睿刀医疗科技有限公司 | 神经网络的训练装置及方法，消融布针规划装置及方法 |
| CN116983085A (zh) * | 2023-09-25 | 2023-11-03 | 上海微创电生理医疗科技股份有限公司 | 一种消融灶尺寸预测系统及消融设备 |
| CN117695005B (zh) * | 2024-02-05 | 2024-05-07 | 浙江伽奈维医疗科技有限公司 | 一种具有布针引导功能的陡脉冲治疗系统 |
| CN118356247B (zh) * | 2024-06-19 | 2024-09-13 | 杭州睿笛生物科技有限公司 | 一种基于单针多极型电极针的脉冲电场肿瘤消融治疗系统和装置 |
| CN118370600B (zh) * | 2024-06-25 | 2024-09-24 | 杭州睿笛生物科技有限公司 | 一种基于单针多极型电极针的术前方案生成系统和装置 |
| CN118542728B (zh) * | 2024-07-29 | 2024-10-15 | 天津市鹰泰利安康医疗科技有限责任公司 | 一种脉管内的不可逆电穿孔消融的系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103717167A (zh) * | 2011-07-28 | 2014-04-09 | 皇家飞利浦有限公司 | 消融规划系统 |
| CN108830017A (zh) * | 2018-07-05 | 2018-11-16 | 上海交通大学 | 基于个体阻抗的射频加热温度场预测方法及其系统 |
| CN110432977A (zh) * | 2019-08-07 | 2019-11-12 | 杭州睿笛生物科技有限公司 | 一种电脉冲消融设备以及适用其仿真方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9198733B2 (en) * | 2008-04-29 | 2015-12-01 | Virginia Tech Intellectual Properties, Inc. | Treatment planning for electroporation-based therapies |
| US10117707B2 (en) * | 2008-04-29 | 2018-11-06 | Virginia Tech Intellectual Properties, Inc. | System and method for estimating tissue heating of a target ablation zone for electrical-energy based therapies |
| US11272979B2 (en) * | 2008-04-29 | 2022-03-15 | Virginia Tech Intellectual Properties, Inc. | System and method for estimating tissue heating of a target ablation zone for electrical-energy based therapies |
| US9283051B2 (en) * | 2008-04-29 | 2016-03-15 | Virginia Tech Intellectual Properties, Inc. | System and method for estimating a treatment volume for administering electrical-energy based therapies |
| EP2413833A1 (en) * | 2009-03-31 | 2012-02-08 | AngioDynamics, Inc. | System and method for estimating a treatment region for a medical treatment device and for interactively planning a treatment of a patient |
| US9700368B2 (en) * | 2010-10-13 | 2017-07-11 | Angiodynamics, Inc. | System and method for electrically ablating tissue of a patient |
| CN105534593B (zh) * | 2014-10-29 | 2019-04-23 | 深圳迈瑞生物医疗电子股份有限公司 | 介入消融模拟系统及方法 |
| CN109157280A (zh) * | 2018-08-10 | 2019-01-08 | 重庆大学 | 不可逆电穿孔组织消融效果动态实时评估设备 |
-
2020
- 2020-04-16 CN CN202010301885.4A patent/CN111529051B/zh active Active
- 2020-08-13 WO PCT/CN2020/108971 patent/WO2021208306A1/zh active Application Filing
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103717167A (zh) * | 2011-07-28 | 2014-04-09 | 皇家飞利浦有限公司 | 消融规划系统 |
| CN108830017A (zh) * | 2018-07-05 | 2018-11-16 | 上海交通大学 | 基于个体阻抗的射频加热温度场预测方法及其系统 |
| CN110432977A (zh) * | 2019-08-07 | 2019-11-12 | 杭州睿笛生物科技有限公司 | 一种电脉冲消融设备以及适用其仿真方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111529051A (zh) | 2020-08-14 |
| WO2021208306A1 (zh) | 2021-10-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111529051B (zh) | 一种预测电脉冲消融区域的系统 | |
| CN111529052B (zh) | 一种预测电脉冲消融区域的系统 | |
| KR102607603B1 (ko) | 종양 치료장의 용량을 정량화하기 위한 전력 손실 밀도 및 관련 측정값의 사용 | |
| US11944821B2 (en) | System and method to estimate region of tissue activation | |
| US7467075B2 (en) | Three-dimensional finite-element code for electrosurgery and thermal ablation simulations | |
| Pavliha et al. | Patient-specific treatment planning of electrochemotherapy: Procedure design and possible pitfalls | |
| EP3300769B1 (en) | Radiation therapy planning optimization and visualization | |
| CA2955962A1 (en) | High definition coloring of heart chambers | |
| KR102280367B1 (ko) | 10-20 시스템 기반의 위치 정보 제공 방법, 장치 및 프로그램 | |
| WO2022100019A1 (zh) | 一种脉冲电场消融参数优化系统 | |
| CN107007279A (zh) | 一种基于堆栈式自编码器的无创心内异常激动点定位方法 | |
| Ding et al. | Treatment planning optimization in irreversible electroporation for complete ablation of variously sized cervical tumors: a numerical study | |
| CN103140875A (zh) | 利用实况反馈对内部组织进行多模态分割的系统和方法 | |
| Forssell et al. | Effect of skull thickness and conductivity on current propagation for noninvasively injected currents | |
| KR20180110857A (ko) | 고강도 집속 초음파의 스캐닝 장치 및 방법 | |
| CN106902480A (zh) | 一种并行量子退火靶点分布计算方法 | |
| Franceries et al. | Solution of Poisson’s equation in a volume conductor using resistor mesh models: application to event related potential imaging | |
| Berthommier et al. | Asymptotic behavior of neural networks and image processing | |
| CN117954050B (zh) | 电极片贴敷方案确定方法及装置、设备及存储介质 | |
| Morán | Classifying Human Atrial Electrograms and Generating Patient-specific Models of the Left Atrial Posterior Wall Using Point Cloud Data to Simulate Electrograms Arising from Different Tissue Substrates | |
| Labarbera | Anisotropic Conductivity and Uncertainty Quantification in Irreversible Electroporation Simulations | |
| Athawale et al. | A Statistical Framework for Visualization of Positional Uncertainty in Deep Brain Stimulation Electrodes | |
| Sannholm | Automated treatment planning in magnetic resonance guided high intensity focused ultrasound | |
| Ferree et al. | Primer on electroencephalography for functional connectivity | |
| JPH01256932A (ja) | 生体内等価電流双極子の表示装置 |
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 | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A system for predicting ablation region of electric pulse Effective date of registration: 20221009 Granted publication date: 20210507 Pledgee: Industrial Bank Co.,Ltd. Shanghai Pengpu Sub branch Pledgor: SHANGHAI REMEDICINE Co.,Ltd. Registration number: Y2022310000282 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 200135 3rd floor, 2123 Pudong Avenue, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai Patentee after: Shanghai RuiDao Medical Technology Co.,Ltd. Address before: 200135 3rd floor, 2123 Pudong Avenue, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai Patentee before: SHANGHAI REMEDICINE Co.,Ltd. |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20231114 Granted publication date: 20210507 Pledgee: Industrial Bank Co.,Ltd. Shanghai Pengpu Sub branch Pledgor: SHANGHAI REMEDICINE Co.,Ltd. Registration number: Y2022310000282 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A system for predicting the ablation area of electric pulses Effective date of registration: 20231123 Granted publication date: 20210507 Pledgee: Industrial Bank Co.,Ltd. Shanghai Wuning Branch Pledgor: Shanghai RuiDao Medical Technology Co.,Ltd. Registration number: Y2023310000770 |