一种消融模拟方法及系统与流程

1.本发明涉及医学数据处理领域，具体而言涉及一种消融模拟方法及系统。


背景技术：

2.消融是指通过物理或化学手段使肿瘤组织失活，利用射频消融技术治疗恶性肿瘤时通常是将不同形状的消融电极针插入到患者恶性肿瘤组织中，利用射频发生仪产生的交变电流使得组织内导电离子和极化分子沿射频电流方向做高速运动产生焦耳热，热能随时间增加逐渐向外传导给肿瘤细胞，利用肿瘤细胞对高温的承受能力较差的特点，完成对肿瘤细胞的原位灭活。
3.消融手术需要将消融针刺入患者的病变器官内，因此手术过程存在一定的危险性，为了能够提前进行手术模拟、教学等，目前已出现部分消融手术的模拟仿真技术，如王笑茹的“肝肿瘤射频消融温度场仿真技术研究综述”、 洪焦的“脑组织射频消融的有限元仿真与分析”等，但在上述现有技术中，消融过程均是通过消融病灶的温度场变化进行模拟的。一方面，上述方面的计算过程非常复杂，另一方面，上述方法完全基于理论，可能存在部分未考虑到的影响因素，模拟结果与真实的消融过程不符。


技术实现要素：

4.为了解决上述背景技术中的问题，本发明提供一种消融模拟方法及系统，采用计算机设备实时评估当前的消融效果，对消融的参数调节提供参考。
5.在本发明的一个方面，提供一种消融模拟方法，其特征在于所述方法包括如下步骤：步骤s01，获取历史消融手术病历中的消融半径、消融功率、消融时间；步骤s02，根据所述消融半径、消融功率、消融时间数据拟合成能量-消融半径曲线；步骤s03，获取待模拟病人目标消融位置三维模型；步骤s04，根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化。
6.进一步地，通过图像识别算法识别所述消融手术病历中的术后图像中的阴影深度达到一定阈值的范围来确定所述消融半径。
7.进一步地，所述根据所述消融半径、消融功率、消融时间数据拟合成能量-消融半径曲线包括：消融针输入的能量为其中为消融针热效率，为消融针功率，为消融时间当消融针功率各阶段的功率不同时输入的能量为其中为每阶段的功率，为每一阶段的时间；对于半径为的消融半径，其吸收的热量为
其中为消融位置的比热，为消融位置升高的温度对于半径为的消融半径，其放出的热量为其中为散热系数，为中心温度与半径r处的温度差；因此根据能量守恒原理可知其中为消融半径为r的病灶的吸热比例的值近似不变，近似不变，近似不变，则有其中为待定系数；根据所述步骤s01得到的消融半径、消融功率、消融时间的数据拟合出参数，得到的曲线为能量-消融半径曲线。
8.进一步地，所述根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化包括：在所述三维模型选择一个消融点，输入消融功率、消融时间，通过能量-消融半径曲线计算出消融半径的大小，根据所述消融半径在三维模型中绘制一个消融球。
9.进一步地，根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化包括: 修改消融功率、消融时间参数，三维模型中的消融球随着消融功率、消融时间参数的改变而改变。
10.本发明另一个方法还提供一种消融模拟系统，包括如下模块：获取模块，用于获取历史消融手术病历中的消融半径、消融功率、消融时间；拟合模块，用于根据所述消融半径、消融功率、消融时间数据拟合成能量-消融半径曲线；建模模块，用于获取待模拟病人目标消融位置三维模型；模拟模块，用于根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化。
11.进一步地，在所述获取模块中通过图像识别算法识别所述消融手术病历中的术后图像中的阴影深度达到一定阈值的范围来确定所述消融半径。
12.进一步地，根据所述消融半径、消融功率、消融时间数据拟合成能量-消融半径曲线包括：消融针输入的能量为其中为消融针热效率，为消融针功率，为消融时间
当消融针功率各阶段的功率不同时输入的能量为其中为每阶段的功率，为每一阶段的时间；对于半径为的消融半径，其吸收的热量为其中为消融位置的比热，为消融位置升高的温度对于半径为的消融半径，其放出的热量为其中为散热系数，为中心温度与半径r处的温度差；因此根据能量守恒原理可知其中为消融半径为r的病灶的吸热比例的值近似不变，近似不变，近似不变，则有其中为待定系数；根据所述获取模块得到的消融半径、消融功率、消融时间的数据拟合出参数，得到的曲线为能量-消融半径曲线。
13.进一步地，根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化包括：在所述三维模型选择一个消融点，输入消融功率、消融时间，通过能量-消融半径曲线计算出消融半径的大小，根据所述消融半径在三维模型中绘制一个消融球。
14.进一步地，根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化包括: 修改消融功率、消融时间参数，三维模型中的消融球随着消融功率、消融时间参数的改变而改变。
15.本发明通过上述技术方案，通过分析现有的消融病例获取真实的消融数据，通过真实的消融数据获取消融过程中消融目标随着消融功率、时间的变化，从而更接近真实地模拟出消融过程。
附图说明
16.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获
得其他的附图。
17.图1 本发明方法示意图；图2 能量-消融半径曲线示意图。
具体实施方式
18.下面，结合附图以及具体实施方式，对发明做出优选的描述。
19.本实施例通过如下步骤解决上术问题:在一个实施例中，参考图1，本发明提供一种消融模拟方法，包括如下步骤：步骤s01，获取历史消融手术病历中的消融半径、消融功率、消融时间。
20.消融手术时将消融针插入目标病灶，消融针开始工作之后，针尖的周围的温度开始上升，由于针尖处的温度上升导致针尖周围形成一个升温的圆球，在医学图像（超声、ct等）中呈现出为一个圆形，并且圆球内部温度高，越往外温度越低，在医学图像中通常呈现为一个由深至浅的阴影。
21.在现有医院中，由于存储空间的限制医院不会存储每一个消融过程（消融手术中的消融范围变化的实时影像视频）的全部过程数据，因此很难获得比较足量的消融全过程图像数据。但目前的医疗系统中有大量消融手术病历的记载，相关的病历记载中会记载消融手术使用的消融功率、消融时间，以及术后检查的医疗图像数据，医疗图像数据就包括了手术后产生的消融范围（消融手术的影响范围，通常是以温度升高至一预定值进行确定的）。
22.为了能够获得足量的数量进行模拟，并且使用足够多的病例以获得比较具有统计意义的数据；本发明首先获取历史消融手术病历中的消融半径、消融功率、消融时间。其中消融半径可以通过计算程序读取病例库中的术后病灶影像数据来获得；本发明一方面可以通过计算机系统进行自动识别消融半径，如通过图像识别算法识别阴影深度达到一定阈值的范围，则该范围内的半径为消融手术半径。另外，还可以采用人工的方式进行勾画，如通过人工进行识别勾画出消融的范围，通过人工测量消融半径。同时还通过计算机程序读取病例文件中医生采用的消融功率、消融时间数据；消融功率、消融时间通常由医生手写或输入计算机系统，可以采用ocr或字符匹配等方式实现具体数据的提取。
23.由于不同病人的器官含水量、脂肪量等不同，因此消融的效果可能会不同，获取的历史消融手术病历为不同病人的数据，可以获得更为广泛的数据，以使得模拟更具有统计意义。获取的历史消融手术可采用数据库、表格等方式进行存储。示例性地，分别获取了不同病历中消融半径为2cm、2.2cm、2.5cm
……
的一系列数据以及对应的消融功率、消融时间。
24.步骤s02，根据所述消融半径、消融功率、消融时间数据拟合成能量-消融半径曲线。
25.消融手术通过消融针在消融位置产生热量，以达到加热消融目标位置的目的，因此消融目标位置的温度为消融针产生的热量减去散热量的热量产生的温升。
26.消融针输入的能量为其中为消融针热效率，为消融针功率，为消融时间
很显然地，当消融针功率各阶段的功率不同时输入的能量为其中为每阶段的功率，为每一阶段的时间。
27.对于半径为的消融半径，其吸收的热量为其中为消融位置的比热，为消融位置升高的温度对于半径为的消融半径，其放出的热量为其中为散热系数，为中心温度与半径r处的温度差；因此根据能量守恒原理可知其中为消融半径为r的病灶的吸热比例的值近似不变，近似不变，近似不变，则有其中为待定系数。
28.其中为功率与时间的乘积，本发明称之为能量。
29.在前述步骤中获取了大量的消融半径r以及其对应的消融功率p、消融时间t的数据，因此通过可以拟合出参数，本发明称拟合出的曲线为能量-消融半径曲线，如图2所示。
30.进一步，由于数据可能存在错误或明显的误差，因此在进行拟合前去除离群数据。
31.在获取能量-消融半径曲线后，可以输入任意的消融功率、消融时间求得消融半径r。
32.步骤s03，获取待模拟病人目标消融位置三维模型。
33.进行消融模拟是为了模拟使用一定的消融参数在病人体内进行消融手术时对人体产生的影响，在手术前观察消融球的大小，以便提前预判消融效果以及消融风险（消融球过小可能达不到手术效果，消融球过大可能影响正常器官或血管）。
34.目标消融位置可以是人体的一个器官对应的部位，如目前消融手术应用比较多的器官有甲状腺、肝等，手术过程中需要通过颈部、腹部刺入消融针，因此目标位置可以是包括了目标器官的颈部、腹部等位置。
35.对手术目标位置进行三维建模可以使用现有技术中的任意方法，如ct扫描三维重建、mri三维重建等，通过现有的三维建模软件得到相应的三维数字图像，三维数字图像可以显示于显示器上以供医务人员查看。
36.步骤s04，根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化。
37.在前述步骤中求得了能量-消融半径曲线，因此可以可以输入不同消融功率、消融时间以求得消融球的大小。
38.在所述三维模型选择一个消融点（消融针的插入点，也就是消融球的球心），消融点的位置可以通医生手动选择，也可以通过计算机程序自动选择；输入消融功率、消融时间，通过能量-消融半径曲线计算出消融半径的大小，根据所述消融半径在三维模型中绘制一个消融球；医生可以通过消融球直观地观察消融范围是否完整地覆盖病灶，观察消融球是否会影响到周围的正常人体组织。
39.进一步地，医生可以修改消融功率、消融时间参数，三维模型中的消融球随着消融功率、消融时间参数的改变而改变，从而直观地查看消融范围的变化，以便医生选择最佳的手术参数。
40.在另一种实施方式中，本发明还公开了一种消融模拟系统，包括如下模块：获取模块，用于获取历史消融手术病历中的消融半径、消融功率、消融时间；拟合模块，用于根据所述消融半径、消融功率、消融时间数据拟合成能量-消融半径曲线；建模模块，用于获取待模拟病人目标消融位置三维模型；模拟模块，用于根据所述能量-消融半径曲线，在所述三维模型中模拟消融手术过程中使用不同消融功率、消融时间进行消融时消融范围的变化。
41.需要说明的是上述消融模拟系统的详细实现原理以及进一步的改进措施都与前述的消融模拟方法相同，本实施例中不再进行详细描述，本领域技术人员可以根据现有技术消融模拟方法在消融模拟系统中进行具体实现。
42.最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。
43.本发明未特别明确的部分模块结构，以现有技术记载的内容为准。本发明在前述背景技术部分以及具体实施例部分提及的现有技术可作为本发明的一部分，用于理解部分技术特征或者参数的含义。本发明的保护范围以权利要求实际记载的内容为准。