一种伽玛刀放射治疗逆向计划设计方法技术

技术编号:10273426 阅读:127 留言:0更新日期:2014-07-31 15:50
本发明专利技术公开一种伽玛刀放射治疗逆向计划设计方法,首先,进行重建靶区和危及器官感兴趣区域三维模型;建立优化模型;确定优化参数与遗传算法参数;整个过程应用图形处理器(GPU)并行计算进行加速,随后根据改进的遗传算法进行求解,包括:基因编码初始化族群;计算评估个体适应度选择繁殖个体;个体繁殖变异更新族群;调整权重,最终得到可行的治疗计划。本发明专利技术的优点为:优化模型具有较好的适应性,可以通过调整目标函数中的权重适应不同的病症;针对伽玛刀治疗计划设计的特点,改进了遗传算法的基因编码,种群初始化,繁殖和变异规则,并根据优化模型提出了一种灵活有效的个体适应度评价函数。

【技术实现步骤摘要】

本专利技术涉及计算机
,具体来说,是。
技术介绍
放射治疗是治疗癌症最主要的手段之一。伽玛刀放射治疗使用伽玛射线照射癌变的肿瘤组织,在尽可能减少损害正常组织的前提下,抑制癌变细胞的生长、繁殖和扩散。为了达到治疗的目的,一般在对病人进行治疗之前,通过治疗计划系统(TPS)制定治疗计划,使革El区内的剂量具有较好的覆盖度(Coverage)和适形度(Conformity);危及器官(OAR)承受尽可能低的剂量(Avoidance);同时兼顾治疗计划的可行性。传统的治疗计划正向设计是一个不断试错修改(trial and error)的过程:物理师根据临床经验设计出初步的治疗计划,然后基于靶区和危及器官内的剂量分布和治疗计划的可行度进行调整,重复这个步骤直到得出一个可行的治疗计划。然而,计划设计通常需要综合权衡考虑多个目标,设计的过程繁重且耗时较长,而且在很大程度上依赖于物理师的临床经验。不同于正向计划设计,逆向计划设计根据预期的剂量分布定义优化目标和约束条件,建立优化模型,从而将计划设计转化为优化问题。借助计算机强大的运算能力,使用高效的优化算法,逆向计划设计可以快速的制定出治疗计划,还可以提高治疗计划的可行度,同时降低了对临床经验的依赖程度。现有的优化方法主要针对头部伽玛刀提出,主要包括非线性优化方法和基于几何特征的优化方法。其中, 针对头部伽玛刀的非线性优化方法,通过查表的方法计算剂量,进而建立优化模型并求解。对于体部伽玛刀,由于体部由于肋骨等骨骼的遮挡,不同位置靶点剂量分布相差较大,不能通过查表法计算剂量分布。此外,非线性优化方法求解时间较长,不适用于靶区较大的病例。由于头部伽玛刀使用直径较小的准直器,作用范围较小,不同靶点之间剂量分布相互影响较小。基于几何特征的优化方法近似的认为单个靶点的剂量分布为球形,不同直径的准直器对应大小不同的球体,在各个球体不相互交叠的前提下,利用尽可能少的球体填充形状体积已知的靶区,使靶区内的覆盖度尽可能的高。但是对于体部伽玛刀伽玛刀,准直器直径较大,距离较近的靶点之间剂量分布相互影响较大,因此,基于几何特征的优化方法也不适用于体部伽玛刀治疗计划设计。
技术实现思路
针对上述问题,本专利技术提出了一种适用于伽玛刀放射治疗逆向计划设计方法,建立了一种基于剂量分布的优化模型,并应用改进的遗传算法(Genetic Algorithm)进行求解。,具体通过下述步骤完成:步骤1:重建感兴趣区域的三维模型,并对靶区部分进行膨胀处理生成靶区边缘轮廓。步骤2:建立优化模型;优化目标为:1、保证靶区剂量的覆盖度;2、提高靶区剂量的适形度;3、降低危及器官承受的剂量。步骤3:确定步骤2建立的优化模型参数。步骤4:确定遗传算法参数。步骤5:基因编码及初始化族群。步骤6:对每个靶点序列进行剂量计算。步骤7:评价族群中各个个体,并通过计算适应度函数得到一个适应度数值。步骤8:选择繁殖个体;将族群中的个体按适应度数值由高到底排序,选取前m个个体两两任意配对,作为父代和母代个体繁殖下一代。m为族群中个体数与繁殖概率的乘积。步骤9:父代和母代个体繁殖;A、确定父代和母代个体的重心,得到重心位移方向向量d ;父代和母代个体的重心位置H为:本文档来自技高网
...

【技术保护点】
一种伽玛刀放射治疗逆向计划设计方法,其特征在于:具体通过下述步骤完成:步骤1:重建感兴趣区域的三维模型,并对靶区部分进行膨胀处理生成靶区边缘轮廓;步骤2:建立优化模型;优化目标为:1、保证靶区剂量的覆盖度;2、提高靶区剂量的适形度;3、降低危及器官承受的剂量;步骤3:确定步骤2建立的优化模型参数;步骤4:确定遗传算法参数;步骤5:基因编码及初始化族群;步骤6:对每个靶点序列进行剂量计算;步骤7:评价族群中各个个体,并通过计算适应度函数得到一个适应度数值;步骤8:选择繁殖个体;将族群中的个体按适应度数值由高到底排序,选取前m个个体两两任意配对,作为父代和母代个体繁殖下一代;m为族群中个体数与繁殖概率的乘积;步骤9:父代和母代个体繁殖;A、确定父代和母代个体的重心,得到重心位移方向向量d;父代和母代个体的重心位置H为:H=Σi=1nφi×Si×Pi---(1)]]>其中,φi为靶点序列中第i个靶点对应的靶点准直器的半径;Si为靶点序列中第i个靶点对应的靶点的使能状态;Pi为靶点序列中第i个靶点的位置坐标;n为靶点序列中靶点数量;B、以d为方向向量,随机的取正向或者反向,整体移动父代和母代中所有靶点,得到子代个体。步骤10:子代个体变异;令子代个体中的靶点i变异的平移方向向量Imi,则:Imi=Σj=1n1(Pi-Pj)×Sj×φj×Rj,j≠i---(2)]]>其中,R随机取±1,j为子代个体中靶点i外的其余靶点;步骤11:将族群中进行繁殖的父代母代个体由其繁殖的子代个体代替,而未进行繁殖的个体保持不变;步骤12:重复步骤5~11进行迭代运算,直至满足终止条件后,进行步骤9;步骤13:调整权重对更新后的族群中最优个体的每个靶点,增加权重,随后计算个体的适应度;若适应度减小,继续增加该靶点权重直至适应度值不再减小;若适应度增大,则需减小该靶点的权重以提高个体适应度;因此将靶点权重重置,并减小权重,直至适应度不再减小。...

【技术特征摘要】
1.一种伽玛刀放射治疗逆向计划设计方法,其特征在于:具体通过下述步骤完成: 步骤1:重建感兴趣区域的三维模型,并对靶区部分进行膨胀处理生成靶区边缘轮廓; 步骤2:建立优化模型; 优化目标为:1、保证靶区剂量的覆盖度;2、提高靶区剂量的适形度;3、降低危及器官承受的剂量; 步骤3:确定步骤2建立的优化模型参数; 步骤4:确定遗传算法参数; 步骤5:基因编码及初始化族群; 步骤6:对每个靶点序列进行剂量计算; 步骤7:评价族群中各个个体,并通过计算适应度函数得到一个适应度数值; 步骤8:选择繁殖个体; 将族群中的个体按适应度数值由高到底排序,选取前m个个体两两任意配对,作为父代和母代个体繁殖下一代为族群中个体数与繁殖概率的乘积; 步骤9:父代和母代个体繁殖; A、确定父代和母代个体的重心,得到重心位移方向向量d; 父代和母代个体的重心位置H为: 2.如权利要求1所述一种伽玛刀放射治疗逆向计划设计方法,其特征在于:对感兴趣区域的三维模型进行采样处理。3.如权利要求1所述一种伽玛刀放射治疗逆向计划设计...

【专利技术属性】
技术研发人员:梁斌周付根刘博郭斌许轩昂
申请(专利权)人:北京航空航天大学
类型:发明
国别省市:北京;11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1