专利名称：治疗计划逆向规划方法和治疗计划系统的制作方法
技术领域：
本发明涉及一种放射治疗规划技术，尤其涉及放射治疗计划逆向规划方法和治疗计划系统。
背景技术：
立体定向放射治疗手术或立体定向放射治疗是放射治疗中常见的两种放射治疗技术，常见的设备是基于钴-60放射源的伽玛刀和基于电子加速器的X刀。前者通常采用多个钴-60放射源聚焦照射的方式，使靶体接受高剂量的均勻照射而周围健康组织受量很低以达到控制或根除病变的目的。利用伽玛刀治疗设备实施放射治疗之前，通常需要制定出一个可接受的放射治疗计划。伽玛刀的治疗计划通常通过手工采用交互迭代的方式进行。这是一个正向规划过程，即医生或物理师根据靶体的体积和形状，采用试错方式，逐步增加靶点数目，交互调整各靶点的位置、准直器大小以及相对权重等参数，直到最终获得一个令人满意的治疗计划。由于伽玛刀可供选择的准直器大小有限，治疗计划通常需要采用多个靶点照射，这样需要调整的参数很多，尤其是当靶体的体积较大且形状不规则，或者靶体邻近有健康组织时，这是一个非常费时的过程，同时对计划设计人员的经验和技能要求很1 。为了解决这个问题提出了治疗计划的逆向规划，即由医生或物理师预先给出放射治疗的若干目标，然后通过数学优化技术，反求出满足这些放射治疗目标的最优的治疗计划。当前治疗计划的逆向规划方法通常是预设一个初始计划，然后根据计划设计人员设定的治疗目标，通过迭代优化方式，对初始计划进行优化以获得一个最优的治疗计划。因此对于逆向规划而言，其迭代过程是否高效成为一个关键问题。
发明内容
本发明要解决的技术问题是提供一种高效的治疗计划逆向规划方法；本发明要解决的另一技术问题是提供一种基于该规划方法的治疗计划系统。本发明要解决的技术问题通过以下技术方案加以解决—种治疗计划逆向规划方法，用于在放射治疗前对患者要进行放射治疗的区域进行剂量规划，包括步骤A 输入患者医学图像；步骤B 根据患者医学图像勾画患者体表、靶体、危及器官的组织轮廓；步骤C 设置治疗计划的逆向规划目标；步骤D 创建初始治疗计划；步骤E 设置迭代优化参数种群大小、种群进化次数；步骤F 对所述初始治疗计划进行随机抖动，产生种群中的个体治疗计划；步骤G 计算所述种群所有个体治疗计划对应的剂量场；步骤H 根据所述剂量场选择与所述逆向规划目标最接近的个体治疗计划作为当前最优计划；步骤I 若所述当前最优计划与所述逆向规划目标的差值小于预设的第一阈值， 转步骤M；否则进入下一步；步骤J 若进化次数达到所述种群进化次数，转到步骤M ；否则进入下一步；步骤K 将所述种群通过交配和/或变异进化到新一代种群，转到步骤G ；步骤M 停止种群进化并输出当前最优计划。其中所述步骤F包括步骤Fl 对所述初始治疗计划所包含的参数进行编码形成第一序列；步骤F2 对所述第一序列的位进行随机取反操作形成第二序列；步骤F3 对所述第二序列进行解码，获得个体治疗计划；步骤F4 重复F2至F3，直到形成种群大小数目的个体治疗计划。其中所述编码包括二进制编码、实数编码或格雷码编码中的任意一种。其中所述步骤K中所述的交配执行以下操作步骤Sl 在种群个体治疗计划中随机选择两个个体治疗计划；步骤S2 随机选择正整数M，M大于0且小于等于个体治疗计划的位长，将所述随机选择的两个个体治疗计划的第M位之后数位互换；步骤S3 重复步骤Sl至S2，直到形成种群大小数目的新一代个体治疗计划。其中所述步骤K中所述变异执行以下操作步骤Tl 在种群个体治疗计划中选择一个体治疗计划；步骤T2 随机选择正整数N，N大于0且小于等于个体治疗计划的位长，将所述选择的个体治疗计划的第N位数位取反；步骤T3 重复步骤Tl至T2，直到遍历种群中的所有个体治疗计划。一种治疗计划系统，用于在放射治疗前对患者要进行放射治疗的区域进行剂量规划，包括设置模块、优化模块、进化模块和迭代模块，所述设置模块用于，输入患者医学图像；根据患者医学图像勾画患者体表、靶体、危及器官的组织轮廓；设置治疗计划的逆向规划目标；创建初始治疗计划；设置迭代优化参数种群大小、种群进化次数；对所述初始治疗计划进行随机抖动，产生种群中的个体治疗计划；所述优化模块用于计算所述种群所有个体治疗计划对应的剂量场；根据所述剂量场选择与所述逆向规划目标最接近的个体治疗计划作为当前最优计划；所述进化模块用于将所述种群通过交配和/或变异进化到新一代种群；所述迭代模块用于若所述当前最优计划与所述逆向规划目标的差值小于于预设的第一阈值，或进化次数达到所述种群进化次数，输出当前最优计划并结束；否则调用进化模块产生新一代种群，再调用优化模块进行优化。其中所述设置模块包括抖动单元和第一迭代单元，所述抖动单元用于对所述初始治疗计划所包含的参数进行编码形成第一序列，对所述第一序列的位进行随机取反操作形成第二序列，对所述第二序列进行解码，获得个体治疗计划；所述第一迭代单元用于迭代调用所述抖动单元，直到形成种群大小数目的个体治疗计划。其中所述编码包括二进制编码、实数编码或格雷码编码中的任意一种。其中所述进化模块包括交配单元和第二迭代单元，所述交配单元用于在种群个体治疗计划中随机选择两个个体治疗计划，随机选择正整数M，M大于0且小于等于个体治疗计划的位长，将所述随机选择的两个个体治疗计划的第M位之后数位互换；所述第二迭代单元用于迭代调用所述交配单元，直到形成种群大小数目的新一代个体治疗计划。其中所述进化模块包括变异单元和第三迭代单元，所述变异单元用于在种群个体治疗计划中选择一个体治疗计划，随机选择正整数 N，N大于0且小于等于个体治疗计划的位长，将所述选择的个体治疗计划的第N位数位取反；所述第三迭代单元用于迭代调用所述变异单元，直到遍历种群中的所有个体治疗计划。由于采用了以上技术方案，使本发明具备的有益效果在于(1)本发明由于采用了基于基因算法的种群进化方法，并通过交配或变异的方法进行迭代，使得整个迭代过程简单、高效。(2)本发明对初始计划进行抖动形成种群中的个体计划，使得初始计划的敏感度降低，从而降低了对初始治疗计划的要求，提高了适应性。(3)本发明采用随机子代交配或随机数位取反变异来产生新一代种群，实现简单、高效。
图1示出根据本发明治疗计划逆向规划方法的一个实施例的流程图；图2示出根据本发明治疗计划逆向规划方法的另一个实施例的流程图；图3示出根据本发明治疗计划逆向规划方法的另一个实施例的靶体外轮廓插值示意图；图4示出根据本发明治疗计划逆向规划方法的另一个实施例的患者3D体素模型的示意图；图5示出根据本发明方法的一个实施例的交配前示意图；图6示出根据本发明方法的一个实施例的交配后示意图；图7示出根据本发明方法的一个实施例的变异前示意图；图8示出根据本发明方法的一个实施例的变异后示意图9示出根据本发明治疗计划系统的一个实施例的结构示意图。
具体实施例方式下面通过具体实施方式
结合附图对本发明作进一步详细说明。图1示出根据本发明治疗计划逆向规划方法的一个实施例的流程图，包括步骤102 输入患者医学图像；步骤104 根据患者医学图像勾画患者体表、靶体、危及器官的组织轮廓；步骤106 设置治疗计划的逆向规划目标；步骤108 创建初始治疗计划；步骤110 设置迭代优化参数种群大小、种群进化次数；步骤112 对初始治疗计划进行随机抖动，产生种群中的个体治疗计划；步骤114 计算种群所有个体治疗计划对应的剂量场；步骤116 根据剂量场选择与逆向规划目标最接近的个体治疗计划作为当前最优计划；步骤118 若当前最优计划与所述逆向规划目标的差值小于预设的第一阈值，转步骤124 ；否则进入下一步；步骤120 若进化次数达到种群进化次数，转到步骤124 ；否则进入下一步；步骤122 将种群通过交配和/或变异进化到新一代种群，转到步骤114 ；步骤124 停止种群进化并输出当前最优计划。剂量规划目标可包括靶体处方剂量、健康组织/危及器官的限制剂量等，也可以通过适形度等指标来表示。第一阈值是预先定义的可容许优化结果与剂量规划目标的差， 例如可定义为0.1等。种群大小是指种群中的个体个数，例如可以设为20个，种群进化次数是指从初始算起允许进化多少代，例如可以设为10。一种实施方式，步骤112包括步骤1120 对初始治疗计划所包含的参数进行编码形成第一序列；步骤1122 对第一序列的位进行随机取反操作形成第二序列；步骤IlM 对第二序列进行解码，获得个体治疗计划；步骤11 重复1120至1124，直到形成种群大小数目的个体治疗计划。上述编码可包括二进制编码、实数编码或格雷码编码中的任意一种。一种实施方式，交配可执行以下操作步骤Sl 在种群个体治疗计划中随机选择两个个体治疗计划；步骤S2 随机选择正整数M，M大于0且小于等于个体治疗计划的位长，将随机选择的两个个体治疗计划的第M位之后数位互换；步骤S3 重复步骤Sl至S2，直到形成种群大小数目的新一代个体治疗计划。本领域技术人员应该理解，交配操作中采用其他的规则选择两个个体治疗计划也是可以实现的。一种实施方式，变异可执行以下操作步骤Tl 在种群个体治疗计划中选择一个体治疗计划；
步骤T2 随机选择正整数N，N大于0且小于等于个体治疗计划的位长，将所述选择的个体治疗计划的第N位数位取反；步骤T3 重复步骤Tl至T2，直到遍历种群中的所有个体治疗计划。该遍历可顺序或采用其他规则来进行。图2示出根据本发明治疗计划逆向规划方法的另一个实施例的流程图，其使用 SGS-II型立体定向伽玛治疗系统进行放射治疗，包括步骤202 输入患者图像，可输入患者的CT或MRI图像序列；步骤204 勾画患者体表、靶体、危及器官等组织轮廓；步骤206 设置治疗计划逆向规划参数靶体PTV处方剂量Dp 一般选择50 %等剂量线健康组织/危及器官OARs限制剂量D。 (k)为危及器官/健康组织的剂量限制步骤208 建立患者3D体素模型；根据用户在定位序列图像上勾画的体表、靶体、危及器官等组织外轮廓，构造患者 3D体素模型。具体方法如下A 组织外轮廓插值通常定位扫描时，采用的层厚或层间距比定位图像的像素尺寸大很多，为了构造患者的3D体素模型，需要定位序列图像上勾画的体表、靶体、危及器官等所有外轮廓插值。 插值采用线性插值，图3所示为靶体外轮廓插值示意图。体表和危及器官等的外轮廓采用相同的方法进行插值。B 通过体素化构造患者3D体素模型对体表、靶体、危及器官等所有外轮廓进行插值后，将这些轮廓体素化即得到患者的3D体素模型。患者3D体素模型通常需要足够高的分辨率以确保后续治疗计划逆向规划获得好的结果。一种可选择的分辨率是采用患者定位图像的分辨率，这个分辨率一般为 0.5mm-lmm。另一种方法确定分辨率的方法是由用户定制分辨率大小。例如在SGS-II的治疗计划逆向规划中，3D体素模型的分辨率采用与剂量计算网格相同的分辨率。这样用户可以通过设置剂量计算矩阵网格的分辨率来调整3D体素模型的分辨率。图4示出一个患者3D体素模型的示意图。步骤210 创建初始“种子”治疗计划“种子”治疗计划用作创建一个种群的“种子”，即一个种群可以由该“种子”计划构
造出来。“种子”治疗计划可以通过手工交互的方式建立。另一种可选方法是通过自动靶点布置技术创建。在SGS的治疗计划逆向规划中，支持上述两种创建初始治疗计划方式。步骤212 治疗计划的逆向规划采用并行遗传算法，通过迭代优化技术进行治疗计划的逆向规划。具体流程如下1.设置迭代优化参数种群大小Np 种群中个体数目种群进化代数凡种群需要进化的代数交叉概率P。遗传优化中个体之间的交叉概率，一般由程序预设，不用用户设置。
变异概率Pm 遗传优化中个体变异概率，一般由程序预设，不用用户设置。2.种群初始化种群由若干个体组成，其中每一个个体对应一个候选的治疗计划。种群初始化就是创建一个包含若干候选治疗计划的初始种群。A 个体治疗计划编码一个治疗计划主要包括如下参数靶点数目Nf、靶点位置Pk(x，y，z)，靶点权重Wk， 靶点准直器规格Ck等。为了适应遗传优化算法，需要对上述参数进行编码。编码方式有多种，可以采用二进制编码、实数编码或者格雷码编码等。在SGS-II逆向治疗计划规划中，采用了标准的二进制编码方法。表1给出了一个个体治疗计划的二进制编码表1治疗计划二进制编码示例
权利要求
1.一种治疗计划逆向规划方法，用于在放射治疗前对患者要进行放射治疗的区域进行剂量规划，其特征在于，包括步骤A 输入患者医学图像；步骤B 根据患者医学图像勾画患者体表、靶体、危及器官的组织轮廓； 步骤C 设置治疗计划的逆向规划目标； 步骤D 创建初始治疗计划；步骤E 设置迭代优化参数种群大小、种群进化次数；步骤F 对所述初始治疗计划进行随机抖动，产生种群中的个体治疗计划；步骤G 计算所述种群所有个体治疗计划对应的剂量场；步骤H:根据所述剂量场选择与所述逆向规划目标最接近的个体治疗计划作为当前最优计划；步骤I 若所述当前最优计划与所述逆向规划目标的差值小于预设的第一阈值，转步骤M ；否则进入下一步；步骤J 若进化次数达到所述种群进化次数，转到步骤M ；否则进入下一步； 步骤K 将所述种群通过交配和/或变异进化到新一代种群，转到步骤G ； 步骤M 停止种群进化并输出当前最优计划。
2.如权利要求1所述的方法，其特征在于，其中所述步骤F包括 步骤Fl 对所述初始治疗计划所包含的参数进行编码形成第一序列； 步骤F2 对所述第一序列的位进行随机取反操作形成第二序列； 步骤F3 对所述第二序列进行解码，获得个体治疗计划；步骤F4 重复F2至F3，直到形成种群大小数目的个体治疗计划。
3.如权利要求2所述的方法，其特征在于，其中所述编码包括二进制编码、实数编码或格雷码编码中的任意一种。
4.如权利要求1所述的方法，其特征在于，其中所述步骤K中所述的交配执行以下操作步骤Sl 在种群个体治疗计划中随机选择两个个体治疗计划； 步骤S2 随机选择正整数M，M大于0且小于等于个体治疗计划的位长，将所述随机选择的两个个体治疗计划的第M位之后数位互换；步骤S3 重复步骤Sl至S2，直到形成种群大小数目的新一代个体治疗计划。
5.如权利要求1至4任一所述的方法，其特征在于，其中所述步骤K中所述变异执行以下操作步骤Tl 在种群个体治疗计划中选择一个体治疗计划；步骤T2 随机选择正整数N，N大于0且小于等于个体治疗计划的位长，将所述选择的个体治疗计划的第N位数位取反；步骤T3 重复步骤Tl至T2，直到遍历种群中的所有个体治疗计划。
6.一种治疗计划系统，用于在放射治疗前对患者要进行放射治疗的区域进行剂量规划，其特征在于，包括设置模块、优化模块、进化模块和迭代模块，所述设置模块用于， 输入患者医学图像；根据患者医学图像勾画患者体表、靶体、危及器官的组织轮廓； 设置治疗计划的逆向规划目标； 创建初始治疗计划；设置迭代优化参数种群大小、种群进化次数； 对所述初始治疗计划进行随机抖动，产生种群中的个体治疗计划； 所述优化模块用于计算所述种群所有个体治疗计划对应的剂量场；根据所述剂量场选择与所述逆向规划目标最接近的个体治疗计划作为当前最优计划；所述进化模块用于将所述种群通过交配和/或变异进化到新一代种群； 所述迭代模块用于若所述当前最优计划与所述逆向规划目标的差值小于于预设的第一阈值，或进化次数达到所述种群进化次数，输出当前最优计划并结束；否则调用进化模块产生新一代种群，再调用优化模块进行优化。
7.如权利要求6所述的系统，其特征在于，其中所述设置模块包括抖动单元和第一迭代单元，所述抖动单元用于对所述初始治疗计划所包含的参数进行编码形成第一序列，对所述第一序列的位进行随机取反操作形成第二序列，对所述第二序列进行解码，获得个体治疗计划；所述第一迭代单元用于迭代调用所述抖动单元，直到形成种群大小数目的个体治疗计划。
8.如权利要求7所述的系统，其特征在于，其中所述编码包括二进制编码、实数编码或格雷码编码中的任意一种。
9.如权利要求6所述的系统，其特征在于，其中所述进化模块包括交配单元和第二迭代单元，所述交配单元用于在种群个体治疗计划中随机选择两个个体治疗计划，随机选择正整数M，M大于0且小于等于个体治疗计划的位长，将所述随机选择的两个个体治疗计划的第 M位之后数位互换；所述第二迭代单元用于迭代调用所述交配单元，直到形成种群大小数目的新一代个体治疗计划。
10.如权利要求6至9任一所述的系统，其特征在于，其中所述进化模块包括变异单元和第三迭代单元，所述变异单元用于在种群个体治疗计划中选择一个体治疗计划，随机选择正整数N，N 大于0且小于等于个体治疗计划的位长，将所述选择的个体治疗计划的第N位数位取反； 所述第三迭代单元用于迭代调用所述变异单元，直到遍历种群中的所有个体治疗计划。
全文摘要
本发明公开了一种治疗计划逆向规划方法，包括A输入患者医学图像；B勾画患者体表、靶体、危及器官的组织轮廓；C设置逆向规划目标；D创建初始治疗计划；E设置迭代优化参数；F产生个体治疗计划；G计算个体治疗计划剂量场；H根据剂量场选择当前最优计划；I若当前最优计划与逆向规划目标的差值小于预设第一阈值，转L；否则进入下一步；J若进化次数达到种群进化次数，转到L；否则进入下一步；K将种群通过交配和/或变异进化到新一代种群，转到G；L停止种群进化并输出当前最优计划。本发明还公开了一种治疗计划系统。本发明采用基于基因算法的种群进化方法，并通过交配或变异的方法进行迭代，使得整个迭代过程简单、高效。
文档编号G06N3/12GK102184318SQ20111009767
公开日2011年9月14日 申请日期2011年4月18日 优先权日2011年4月18日
发明者刘启平, 卿侯, 崔智 申请人:深圳市海博科技有限公司