基于扩展收缩理论的放疗计划仿真设计方法及系统

技术领域

本发明涉及虚拟手术技术领域，特别是涉及一种基于扩展收缩理论的放疗计划仿真设计方法及系统。

背景技术

放射治疗在肿瘤治疗中的作用和地位日益突出，已成为治疗恶性肿瘤的主要手段之一。据统计，大约70％的癌症患者在治疗癌症的过程中需要用放射治疗，约有40％的癌症可以用放疗根治。肿瘤放射治疗是利用高能量放射线照射癌变的肿瘤，杀死或破坏癌细胞，抑制它们的生长、繁殖和扩散。

放射治疗根据放射源距患者的距离远近可以分为近距离放射治疗和外照射。近距离放疗手段是通过施源器或密封放射源将放射性粒子通过介入穿刺的方法植入肿瘤内，利用从放射性粒子持续发出的射线对病变区持续杀伤照射。放射性粒子植入使用尺寸为毫米级别的小型放射源(如碘125等)。为了合理放置粒子源并便于计算，将每根穿刺针离散化为多个驻留点，放射源可以在驻留位置作用一定的驻留时间。

传统的放疗手术计划仿真是一个不断试错的过程，物理师凭借经验根据放疗剂量分布构建放疗计划，对放疗计划进行评估，如果不满足手术要求，医师将调整放疗剂量分布。不断重复此过程，直至放疗计划满足手术要求的剂量分布。这种方法对物理师经验依赖程度较高，且获得的放疗计划难以达到最优的剂量分布。

鉴于此，本发明提出一种基于扩展收缩理论的放疗计划仿真设计方法及系统，以缓解现有技术的不足。

发明内容

第一方面，本发明提供了一种基于扩展收缩理论的放疗计划仿真设计方法，包括：在医学影像中标记肿瘤目标，建立包含肿瘤目标TV区域和OAR区域的三维影像，对三维影像均匀采样，获取目标点云数据；设置包含放射治疗入针点的模板，任一入针点与目标点云数据中任一点组合为针道，随机生成针道序列P，序列P中的任意两个针道所在线段均无交点；在序列P中的针道设置若干放射粒子，计算若干放射粒子的剂量分布；对放射粒子的驻留点进行扩展收缩操作，调整放射粒子驻留时间，以优化放射粒子的剂量分布；利用3D打印技术，打印三维影像的三维模型、针道序列P和放射粒子的驻留点。

进一步地，对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间的步骤，包括：在针道序列P中未设置放射粒子的任一驻留点进行扩展操作，扩展操作包括在驻留点设置放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，以优化放射粒子的剂量分布；在针道序列P中已设置放射粒子的任一驻留点进行收缩操作，收缩操作包括在驻留点撤销放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，若收缩操作不能优化放射粒子的剂量分布，则不进行收缩操作；重复进行扩展和收缩操作若干次，或通过扩展操作不能优化放射粒子的剂量分布时，停止扩展和收缩操作。

进一步地，在对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间的之后，还包括：记录剂量分布的目标值为第一优化值；随机替换序列P中的任一针道，替换后序列P中的作意两个针道所在线段均无交点；对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间，以获取针道替换后剂量分布的目标值，若针道替换后剂量分布的目标值优于第一优化值，将针道替换后剂量分布的目标值赋予第一优化值；继续随机替换序列P中的任一针道，直至第一优化值不再变化或针道替换次数超过第一阈值。

第二方面，本发明提供了一种基于扩展收缩理论的放疗计划仿真设计系统，包括：计算设备和3D打印设备，计算设备包括影像处理模块，计划生成模块和计划优化模块；影像处理模块，在医学影像中标记肿瘤目标，建立包含肿瘤目标TV区域和OAR区域的三维影像，对三维影像均匀采样，获取目标点云数据；计划生成模块，设置包含放射治疗入针点的模板，任一入针点与目标点云数据中任一点组合为针道，随机生成针道序列P，序列P中的任意两个针道所在线段均无交点；在序列P中的针道设置若干放射粒子，计算若干放射粒子的剂量分布；计划优化模块，对放射粒子的驻留点进行扩展收缩操作，调整放射粒子驻留时间，以优化放射粒子的剂量分布；3D打印设备，利用3D打印技术，打印三维影像的三维模型、针道序列P和放射粒子的驻留点。

进一步地，计划优化模块对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间的步骤，包括：在针道序列P中未设置放射粒子的任一驻留点进行扩展操作，扩展操作包括在驻留点设置放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，以优化放射粒子的剂量分布；在针道序列P中已设置放射粒子的任一驻留点进行收缩操作，收缩操作包括在驻留点撤销放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，若收缩操作不能优化放射粒子的剂量分布，则不进行收缩操作；重复进行扩展和收缩操作若干次，或通过扩展操作不能优化放射粒子的剂量分布时，停止扩展和收缩操作。

进一步地，计划优化模块被进一步配置为：记录剂量分布的目标值为第一优化值；随机替换序列P中的任一针道，替换后序列P中的作意两个针道所在线段均无交点；对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间，以获取针道替换后剂量分布的目标值，若针道替换后剂量分布的目标值优于第一优化值，将针道替换后剂量分布的目标值赋予第一优化值；继续随机替换序列P中的任一针道，直至第一优化值不再变化或针道替换次数超过第一阈值。

本发明的有益效果如下：

本发明提供的技术方案可以包括以下有益效果：提出了一种基于扩展收缩理论的放疗计划仿真设计方法及系统，自动生成放射治疗手术实施计划，对放疗计划的剂量分布进行仿真优化，缓解了传统的放射治疗计划对物理师的经验高度依赖，提高了手术计划的制定效率。利用基于扩展收缩理论的优化方法对针道中放射粒子的驻留位置进行迭代优化，提升了治疗计划制定的效能。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本公开。

附图说明

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一种实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明第一实施例一种基于扩展收缩理论的放疗计划仿真设计方法流程示意图；

图2为本发明第二实施例一种基于扩展收缩理论的放疗计划仿真设计系统结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，所描述的实施例是本发明一部分实施例，而不是全部的实施例。

第一实施例：

图1为本发明第一实施例一种基于扩展收缩理论的放疗计划仿真设计方法流程示意图，如图1所示，该方法包括如下四个步骤。

步骤S11：在医学影像中标记肿瘤目标。具体地，在医学影像中标记肿瘤目标，建立包含肿瘤目标TV区域和OAR区域的三维影像，对三维影像均匀采样，获取目标点云数据。

需要进行说明的是，TV(Target Volume)治疗靶区为CTV(Clinical TargetVolume)临床靶区，包括已确定存在的肿瘤和潜在的受侵组织，GTV(Gross Tumor Volume)肿瘤区和周围的亚临床病灶构成CTV，放射治疗的目的在于杀灭TV区域的肿瘤细胞。OAR(Organ At Risk)危及器官区域是指放疗区域周边的正常器官，通常在放射治疗中会受到影响。

步骤S12：生成放射治疗手术计划。具体地，设置包含放射治疗入针点的模板，任一入针点与目标点云数据中任一点组合为针道，随机生成针道序列P，序列P中的任意两个针道所在线段均无交点；在序列P中的针道设置若干放射粒子，计算若干放射粒子的剂量分布。

步骤S13：优化放射治疗手术计划。具体地，对放射粒子的驻留点进行扩展收缩操作，调整放射粒子驻留时间，以优化放射粒子的剂量分布。

在一个具体的实施例中，对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间的步骤，包括：在针道序列P中未设置放射粒子的任一驻留点进行扩展操作，扩展操作包括在驻留点设置放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，以优化放射粒子的剂量分布；在针道序列P中已设置放射粒子的任一驻留点进行收缩操作，收缩操作包括在驻留点撤销放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，若收缩操作不能优化放射粒子的剂量分布，则不进行收缩操作；重复进行扩展和收缩操作若干次，或通过扩展操作不能优化放射粒子的剂量分布时，停止扩展和收缩操作。

在一个可选的实施例中，在对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间的之后，还包括：记录剂量分布的目标值为第一优化值；随机替换序列P中的任一针道，替换后序列P中的作意两个针道所在线段均无交点；对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间，以获取针道替换后剂量分布的目标值，若针道替换后剂量分布的目标值优于第一优化值，将针道替换后剂量分布的目标值赋予第一优化值；继续随机替换序列P中的任一针道，直至第一优化值不再变化或针道替换次数超过第一阈值。

步骤S14：打印三维模型和针道数据。具体地，利用3D打印技术，打印三维影像的三维模型、针道序列P和放射粒子的驻留点。

需要进行说明的是，打印实体的三维模型和针道数据，可以更直观地复现手术计划。

第二实施例：

图2是本发明实施例一种基于扩展收缩理论的放疗计划仿真设计系统结构示意图，如图2所示，系统包括：计算设备100和3D打印设备200，计算设备100包括影像处理模块101，计划生成模块102和计划优化模块103。

影像处理模块101，在医学影像中标记肿瘤目标，建立包含肿瘤目标TV区域和OAR区域的三维影像，对三维影像均匀采样，获取目标点云数据。

计划生成模块102，设置包含放射治疗入针点的模板，任一入针点与目标点云数据中任一点组合为针道，随机生成针道序列P，序列P中的任意两个针道所在线段均无交点；在序列P中的针道设置若干放射粒子，计算若干放射粒子的剂量分布。

计划优化模块103，对放射粒子的驻留点进行扩展收缩操作，调整放射粒子驻留时间，以优化放射粒子的剂量分布。

在一个具体的实施例中，对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间的步骤，包括：在针道序列P中未设置放射粒子的任一驻留点进行扩展操作，扩展操作包括在驻留点设置放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，以优化放射粒子的剂量分布；在针道序列P中已设置放射粒子的任一驻留点进行收缩操作，收缩操作包括在驻留点撤销放射粒子，调整针道序列P中已放置的放射粒子的驻留时间，若收缩操作不能优化放射粒子的剂量分布，则不进行收缩操作；重复进行扩展和收缩操作若干次，或通过扩展操作不能优化放射粒子的剂量分布时，停止扩展和收缩操作。

在一个可选的实施例中，计划优化模块103被进一步配置为：记录剂量分布的目标值为第一优化值；随机替换序列P中的任一针道，替换后序列P中的作意两个针道所在线段均无交点；对放射粒子驻留点进行扩展收缩操作，调整放射粒子驻留时间，以获取针道替换后剂量分布的目标值，若针道替换后剂量分布的目标值优于第一优化值，将针道替换后剂量分布的目标值赋予第一优化值；继续随机替换序列P中的任一针道，直至第一优化值不再变化或针道替换次数超过第一阈值。

3D打印设备200，利用3D打印技术，打印三维影像的三维模型、针道序列P和放射粒子的驻留点。

最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。