技术领域
[0001] 本
发明涉及
放射治疗技术领域，具体涉及一种预测剂量引导的TPS优化方法及装置。
背景技术
[0002] 调强
放射治疗(Intensity modulated radiation therapy，IMRT)作为放疗的一种，针对靶区三维形状和危及器官与靶区的空间
位置关系，对射束强度进行调节，使放射剂量的分布与靶区形状一致，靶区内接受分布均匀的高的剂量，保证对
肿瘤细胞的杀伤，同时使得正常组织接受较低的剂量照射，达到保护正常器官、减少不良反应的目的。
[0003] 传统的两步法IMRT计划制作过程一般包括两个步骤，按先后顺序依次为：步骤1.优化射野强度；步骤2.子野分割。步骤1中，计划设计者根据经验设置射束方向，设置器官约束条件，使用逆向优化
算法求出射野强度图及对应的优化剂量，当优化剂量不满足治疗要求时，需调整约束条件，再次优化，如此反复，直至优化剂量满足治疗要求。在步骤2中，将步骤1得到的优化强度图作为输入，执行子野分割算法，由于光栅的特性和物理限制，譬如光栅运动最大速度、光栅半影参数等，分割后的剂量与优化剂量的分布存在一定的差异，当分割剂量不满足治疗要求时，需返回步骤1，调整约束条件，执行优化，然后再进行子野分割，如此反复，直至分割剂量满足治疗要求。整体而言，传统的两步法IMRT是个不断试错的过程，费时费
力，尤其是对经验较少的物理师更是耗时。
发明内容
[0004] 本发明的目的在于，针对上述
现有技术中的不足，提供一种预测剂量引导的TPS(放射治疗计划系统)优化方法及装置，以解决物理师在放射治疗设计优化阶段的耗时问题。
[0005] 为实现上述目的，本发明采用的技术方案如下：
[0006] 第一方面，本发明提供了一种预测剂量引导的TPS优化方法，用于优化计算放射治疗计划阶段的子野剂量分割，所述方法包括：
[0007] 针对预定肿瘤病人，采用预定剂量
预测模型计算放射治疗中所需要的剂量作为预测剂量；
[0008] 根据所述预测剂量，计算放射治疗的射野强度图贡献作为预测强度图；
[0009] 根据所述预测强度图，
迭代计算子野剂量分割，以获得优化的子野剂量分割结果。
[0010] 可选地，所述针对预定肿瘤病人，采用预定剂量预测模型计算放射治疗中所需要的剂量作为预测剂量，包括：
[0011] 采用训练好的
卷积神经网络模型作为预定剂量预测模型；
[0012] 获取预定肿瘤病人的器官勾画数据，并且输入所述预定剂量预测模型，以得到模型预测结果；
[0013] 将所述模型预测结果进行同步推量，使得放射治疗靶区的95％体积的接收剂量作为预测剂量。
[0014] 可选地，所述根据所述预测剂量，计算放射治疗的射野强度图贡献作为预测强度图，包括：
[0015] 设置目标函数：
[0016]
[0017]
[0018] 其中，Fobj为目标函数，k为迭代次数，Npt为病人体数据中
像素点个数，ωj为像素点的权重，dj为第j个点的计算剂量，pj为第j个点的预测剂量，Nray为射束的总个数，ajm为单位k强度的第m个射束对第j个点的剂量贡献， 为第m个射束的强度，x为所有射束的总强度；
[0019] 使用共轭梯度法，使所述目标函数最小化，从而得到预测剂量的射野强度图贡献作为预测强度图。
[0020] 可选地，所述根据所述预测强度图，迭代计算子野剂量分割，以获得优化的子野剂量分割结果，包括：
[0021] a、设置最大迭代次数，并且设置子野分割强度图、迭代次数、迭代目标函数和优化分割结果的初始值，其中所述子野分割强度图的初始值为所述预测强度图；
[0022] b、针对所述子野分割强度图进行子野剂量分割，得到分割剂量；
[0023] c、计算所述分割剂量与所述预测剂量之间的加权平方误差；
[0024] d、若所述加权平方误差小于所述迭代目标函数的值，则将所述加权平方误差作为所述迭代目标函数的值，并且更新优化分割结果，所述迭代次数的值增加1；若所述迭代目标函数小于所述预测剂量与所述预测强度图对应的沉积剂量的加权平方误差或者所述迭代次数的值达到最大迭代次数，则迭代停止；否则计算所述预测剂量与所述分割剂量的剂量差，并且根据所述剂量差更新所述子野分割强度图，并转到步骤b。
[0025] 可选地，所述根据所述剂量差更新所述子野分割强度图，具体包括：
[0026] 计算所述剂量差的强度图贡献；
[0027] 将子野分割强度图减去预定系数与所述剂量差的强度图贡献乘积的差作为新的子野分割强度图。
[0028] 可选地，所述预定系数为0.382。
[0029] 可选地，所述迭代次数的初始值为零、所述迭代目标函数的初始值为无穷大、所述优化分割结果的初始值为空。
[0030] 第二方面，本发明提供了一种预测剂量引导的TPS优化装置，用于优化计算放射治疗计划阶段的子野剂量分割，所述装置包括：
[0031] 预测
剂量计算模
块，用于针对预定肿瘤病人，采用预定剂量预测模型计算放射治疗中所需要的剂量作为预测剂量；
[0032] 预测强度图计算模块，用于根据预测剂量，计算放射治疗的射野强度图贡献作为预测强度图；
[0033] 子野剂量分割计算模块，用于根据预测强度图，迭代计算子野剂量分割，以获得优化的子野剂量分割结果。
[0034] 可选地，预测剂量计算模块，具体用于：
[0035] 采用训练好的卷积神经网络模型作为预定剂量预测模型；
[0036] 获取预定肿瘤病人的器官勾画数据，并且输入预定剂量预测模型，以得到模型预测结果；
[0037] 将模型预测结果进行同步推量，使得放射治疗靶区的95％体积的接收剂量作为预测剂量。
[0038] 可选地，预测强度图计算模块，具体用于：
[0039] 设置目标函数：
[0040]
[0041]
[0042] 其中，Fobj为目标函数，k为迭代次数，Npt为病人体数据中像素点个数，ωj为像素点的权重，dj为第j个点的计算剂量，pj为第j个点的预测剂量，Nray为射束的总个数，ajm为单位k强度的第m个射束对第j个点的剂量贡献， 为第m个射束的强度，x为所有射束的总强度；
[0043] 使用共轭梯度法，使所述目标函数最小化，从而得到预测剂量的射野强度图贡献作为预测强度图。
[0044] 可选地，子野剂量分割计算模块，具体用于：
[0045] a、设置最大迭代次数，并且设置子野分割强度图、迭代次数、迭代目标函数和优化分割结果的初始值，其中所述子野分割强度图的初始值为所述预测强度图；
[0046] b、针对所述子野分割强度图进行子野剂量分割，得到分割剂量；
[0047] c、计算所述分割剂量与所述预测剂量之间的加权平方误差；
[0048] d、若所述加权平方误差小于所述迭代目标函数的值，则将所述加权平方误差作为所述迭代目标函数的值，并且更新优化分割结果，所述迭代次数的值增加1；若所述迭代目标函数小于所述预测剂量与所述预测强度图对应的沉积剂量的加权平方误差或者所述迭代次数的值达到最大迭代次数，则迭代停止；否则计算所述预测剂量与所述分割剂量的剂量差，并且根据所述剂量差更新所述子野分割强度图，并转到步骤b。
[0049] 可选地，子野剂量分割计算模块，还具体用于：
[0050] 计算所述剂量差的强度图贡献；
[0051] 将子野分割强度图减去预定系数与所述剂量差的强度图贡献乘积的差作为新的子野分割强度图。
[0052] 可选地，预定系数为0.382。
[0053] 可选地，迭代次数的初始值为零，迭代目标函数的初始值为无穷大、优化分割结果的初始值为空。
[0054] 第三方面，本发明提供了一种放射治疗装置，所述装置采用根据第一方面所述的预测剂量引导的TPS优化方法来优化计算放射治疗计划阶段的子野剂量分割，或者包括根据第二方面所述的预测剂量引导的TPS优化装置。
[0055] 本发明的有益效果包括：
[0056] 本发明提供的预测剂量引导的TPS优化方法，用于优化计算放射治疗计划阶段的子野剂量分割，该方法包括：针对预定肿瘤病人，采用预定剂量预测模型计算放射治疗中所需要的剂量作为预测剂量；根据预测剂量，计算放射治疗的射野强度图贡献作为预测强度图；根据预测强度图，迭代计算子野剂量分割，以获得优化的子野剂量分割结果。通过根据预测强度图，迭代计算子野剂量分割结果，直至分割剂量与预测剂量尽可能小，从而可以自动获得优化的分割结果，避免了放射治疗物理师进行子野分割的试错过程，提高了物理师的工作效率。
附图说明
[0057] 为了更清楚地说明本发明
实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0058] 图1示出了本发明实施例提供的剂量预测模型的结构示意图；
[0059] 图2示出了本发明实施例提供的剂量预测的流程示意图；
[0060] 图3示出了本发明实施例提供的子野分割迭代优化方法的流程示意图；
[0061] 图4示出了本发明实施例提供的剂量预测模型在前四个测试病例上的表现结果。
具体实施方式
[0062] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0063] 本发明提供了一种预测剂量引导的TPS(放射治疗计划系统)优化方法，用于优化计算放射治疗计划阶段的子野剂量分割，从而减少物理师在优化阶段的时间消耗。
[0064] 该方法包括：针对预定肿瘤病人，采用预定剂量预测模型计算放射治疗中所需要的剂量(例如，经优化的最优剂量)作为预测剂量；根据预测剂量，计算放射治疗的射野强度图贡献作为预测强度图；根据预测强度图，迭代计算子野剂量分割，以获得优化的子野剂量分割结果。
[0065] 为了获得预测剂量，首先采用训练好的卷积神经网络模型作为预定剂量预测模型；获取预定肿瘤病人的器官勾画数据，并且输入预定剂量预测模型，以得到模型预测结果；将模型预测结果进行同步推量，使得放射治疗靶区的95％体积的接收剂量作为预测剂量。下面将对获得预测剂量的过程进行详细描述。
[0066] 可以选择卷积神经网络模型(模型的结构参见图1)作为肿瘤病人的剂量预测模型，下面以
前列腺癌病人为例叙述模型的训练以及预测流程，预测流程如图2所示。
[0067] 模型的输入是CT(
电子计算机
断层扫描)切片上的勾画数据，大小为256*256*6,256*256为切片大小，6为通道数，对应器官个数。通道图像是一个像素值只取0、1的二值图像，1代表CT切片上当前像素点属于该通道对应的器官，0则反之。模型的输出是一个大小为
256*256的浮点型图像，像素值为CT切片上对应像素点处的剂量预测值。
[0068] 选择88例临床前列腺IMRT病人，要求计划有相同的
加速器类型及
能量，以及相同的射野
角度。另外要求每个病人有6个器官的勾画：靶区(PTV)、膀胱、
皮肤、左股骨头、右股骨头、直肠。从88例病人中选择80例作为训练集，余下的8例作为测试集。
[0069] 使用剪裁与插值方法，使病人的体数据大小为256*256*64，空间
分辨率大小为2*2*2.5。对所有病人的体剂量进行归一化，使得所有病人PTV的平均剂量归一到1。
[0070] 将
训练数据(80例)分为训练集、验证集两个部分，训练集72例病人，验证集8例病人，训练集中的数据用于训练，验证集中的数据用于筛选权重，防止模型过拟合。训练时采用的loss函数(损失函数)为预测剂量与真实剂量的均方误差。为了防止模型过拟合，使用交叉验证的方法训练10个模型(训练集病人72例，验证集中的病人8例)，选择一个在验证集上loss值最小的模型。另外训练时，为了突出中间切片的重要性，可以采用高斯抽样。通过上述过程可以获得训练好的剂量预测模型。
[0071] 然后获取病人的勾画数据，将其输入到训练好的剂量预测模型以获得初始预测结果，将初始预测结果进行同步推量，使得PTV的95％体积的接收剂量为
处方剂量，从而获得预测剂量。具体地，可以通过将初始预测结果乘以如下的因子来得到预测剂量：
[0072]
[0073] 其中处方剂量为由专业治疗人员(例如，医师)对病人开具的剂量。
[0074] 然后，根据预测剂量，计算放射治疗的射野强度图贡献作为预测强度图，具体包括：
[0075] 设置目标函数：
[0076]
[0077]
[0078] 其中，Fobj为目标函数，k为迭代次数，Npt为病人体数据中像素点个数，ωj为像素点的权重，dj为第j个点的计算剂量，pj为第j个点的预测剂量，Nray为射束的总个数，ajm为单位强度的第m个射束对第j个点的剂量贡献， 为第m个射束的强度，xk为所有射束的总强度；
[0079] 使用共轭梯度法，使目标函数最小化，从而得到预测剂量的射野强度图贡献作为预测强度图。
[0080] 根据预测强度图迭代计算子野剂量分割以获得优化的子野剂量分割结果，具体包括：a、设置最大迭代次数，并且设置子野分割强度图、迭代次数、迭代目标函数和优化分割结果的初始值，其中所述子野分割强度图的初始值为所述预测强度图；b、针对子野分割强度图进行子野剂量分割，得到分割剂量；c、计算分割剂量与预测剂量之间的加权平方误差；d、若加权平方误差小于迭代目标函数的值，则将加权平方误差作为迭代目标函数的值，并且更新优化分割结果，迭代次数的值增加1；若迭代目标函数小于预测剂量与预测强度图对应的沉积剂量的加权平方误差或者迭代次数的值达到最大迭代次数，则迭代停止；否则计算预测剂量与分割剂量的剂量差，并且根据剂量差更新所述子野分割强度图，并转到步骤b继续进行。其中，迭代次数的初始值为零、迭代目标函数的初始值为无穷大、优化分割结果的初始值为空。
[0081] 根据剂量差更新子野分割强度图具体包括：计算剂量差的强度图贡献；将子野分割强度图减去预定系数与所述剂量差的强度图贡献乘积的差作为新的子野分割强度图。其中预定系数可以为0.382。
[0082] 具体地，参见图3，设Dosepred为预测剂量，Mappred为预测强度图，Fobjpred为预测剂量与预测强度图对应的沉积剂量的加权平方误差，即上述式(1)中的强度图取预测强度图时目标函数的取值，迭代计算分割结果的步骤如下：
[0083] 步骤1：设置迭代次数为iter，最大迭代次数itermax，初始化强度图Map＝Mappred，iter＝0，迭代目标函数Fobj＝inf(无穷大)，最优分割结果置为空；
[0084] 步骤2：对强度图Map进行子野分割，计算分割剂量Dosedel、分割强度图Mapdel；
[0085] 步骤3：计算预测剂量与分割剂量的加权平方误差Fobjdel，若Fobjdel＜Fobj，则令Fobj＝Fobjdel，更新最优分割结果，iter自增1；
[0086] 步骤4：若Fobj＜Fobjpred，或者iter＞＝itermax，迭代停止；否则，计算预测剂量Dosepred与分割剂量Dosedel的剂量差Dosediff，可以使用前述计算强度图贡献的方法来计算Dosediff的强度图贡献Mapdiff，并根据如下式来更新Map，
[0087] Map＝Mappred-0.382*Mapdiff
[0088] 并且转到步骤2继续进行。
[0089] 综上所述，本发明实施例提供的优化方法通过根据预测强度图，迭代计算子野剂量分割结果，直至分割剂量与预测剂量尽可能小，从而可以自动获得优化的分割结果，避免了放射治疗物理师进行子野分割的试错过程，节省了物理师的操作时间，提高了物理师的工作效率。
[0090] 下面以前列腺癌病人为例，具体描述实施过程。
[0091] 首先选择88例临床前列腺IMRT病人。计划的射野角度一致，都是7野共面，射野角度依次为150、70、40、0、320、290、210；计划的加速器类型一致，能量都是15MV。每个病人有6个器官的勾画：靶区(PTV)、膀胱、皮肤、左股骨头、右股骨头、直肠。从88例病人中选择80例作为训练集，余下的8例作为测试集。
[0092] 使用收集的数据训练剂量预测模型，本实施例中选取的模型为图1所示的卷积神经网络模型。模型在前4个测试病人上的表现参见图4(每个病人只
抽取了模型在PTV中间切片上表现)，从图4可以看出，模型预测结果与实际结果非常接近。
[0093] 模型训练好之后，为计划设置7个共面野，角度依次为150、70、40、0、320、290、210。使用剂量预测模型，预测病人的最优剂量，记为Dosepred；计算预测剂量的强度图贡献(称为预测强度图)，记为Mappred，计算预测剂量与预测强度图对应的沉积剂量的加权平方误差Fobjpred，最后迭代计算分割结果，迭代计算分割结果的步骤如上述步骤1至步骤4中所述，在此不再赘述。
[0094] 本发明还提供了一种预测剂量引导的TPS优化装置，用于优化计算放射治疗计划阶段的子野剂量分割，具体地，该装置用于执行本发明上述实施例中提供的方法。
[0095] 该装置包括：预测剂量计算模块，用于针对预定肿瘤病人，采用预定剂量预测模型计算放射治疗中所需要的剂量作为预测剂量；预测强度图计算模块，用于根据预测剂量，计算放射治疗的射野强度图贡献作为预测强度图；子野剂量分割计算模块，用于根据预测强度图，迭代计算子野剂量分割，以获得优化的子野剂量分割结果。
[0096] 可选地，预测剂量计算模块，具体用于：采用训练好的卷积神经网络模型作为预定剂量预测模型；获取预定肿瘤病人的器官勾画数据，并且输入预定剂量预测模型，以得到模型预测结果；将模型预测结果进行同步推量，使得放射治疗靶区的95％体积的接收剂量作为预测剂量。
[0097] 可选地，预测强度图计算模块，具体用于：
[0098] 设置目标函数：
[0099]
[0100]
[0101] 其中，Fobj为目标函数，k为迭代次数，Npt为病人体数据中像素点个数，ωj为像素点的权重，dj为第j个点的计算剂量，pj为第j个点的预测剂量，Nray为射束的总个数，ajm为单位强度的第m个射束对第j个点的剂量贡献， 为第m个射束的强度，xk为所有射束的总强度；
[0102] 使用共轭梯度法，使所述目标函数最小化，从而得到预测剂量的射野强度图贡献作为预测强度图。
[0103] 可选地，子野剂量分割计算模块，具体用于：a、设置最大迭代次数，并且设置子野分割强度图、迭代次数、迭代目标函数和优化分割结果的初始值，其中所述子野分割强度图的初始值为所述预测强度图；b、针对所述子野分割强度图进行子野剂量分割，得到分割剂量；c、计算所述分割剂量与所述预测剂量之间的加权平方误差；d、若所述加权平方误差小于所述迭代目标函数的值，则将所述加权平方误差作为所述迭代目标函数的值，并且更新优化分割结果，所述迭代次数的值增加1；若所述迭代目标函数小于所述预测剂量与所述预测强度图对应的沉积剂量的加权平方误差或者所述迭代次数的值达到最大迭代次数，则迭代停止；否则计算所述预测剂量与所述分割剂量的剂量差，并且根据所述剂量差更新所述子野分割强度图，并转到步骤b。
[0104] 可选地，子野剂量分割计算模块，还具体用于：计算所述剂量差的强度图贡献；将子野分割强度图减去预定系数与所述剂量差的强度图贡献乘积的差作为新的子野分割强度图。可选地，预定系数为0.382。可选地，迭代次数的初始值为零，迭代目标函数的初始值为无穷大、优化分割结果的初始值为空。
[0105] 另外，本发明还提供了一种放射治疗装置，所述装置采用本发明上述提供的预测剂量引导的TPS优化方法来优化计算放射治疗计划阶段的子野剂量分割，或者包括本发明上述提供的预测剂量引导的TPS优化装置。
[0106] 上述实施例只为说明本发明的技术构思及特点，其目的在于让本领域普通技术人员能够了解本发明的内容并加以实施，并不能以此限制本发明的保护范围，凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围内。
