CN108175953B - 一种基于动态光栅的旋转调强优化方法及放射治疗设备 - Google Patents

一种基于动态光栅的旋转调强优化方法及放射治疗设备 Download PDF

Info

Publication number
CN108175953B
CN108175953B CN201711498167.5A CN201711498167A CN108175953B CN 108175953 B CN108175953 B CN 108175953B CN 201711498167 A CN201711498167 A CN 201711498167A CN 108175953 B CN108175953 B CN 108175953B
Authority
CN
China
Prior art keywords
intensity
subfield
field
oar
target
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711498167.5A
Other languages
English (en)
Other versions
CN108175953A (zh
Inventor
姚毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Linatech Medical Science And Technology Co ltd
Original Assignee
Suzhou Linatech Medical Science And Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Suzhou Linatech Medical Science And Technology Co ltd filed Critical Suzhou Linatech Medical Science And Technology Co ltd
Priority to CN201711498167.5A priority Critical patent/CN108175953B/zh
Publication of CN108175953A publication Critical patent/CN108175953A/zh
Application granted granted Critical
Publication of CN108175953B publication Critical patent/CN108175953B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103—Treatment planning systems
    • A61N5/1031—Treatment planning systems using a specific method of dose optimization
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1071—Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A61N2005/1054—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using a portal imaging system
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048—Monitoring, verifying, controlling systems and methods
    • A61N5/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A61N2005/1055—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using magnetic resonance imaging [MRI]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00—Radiation therapy
    • A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1092—Details

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

本发明公开了一种基于动态光栅的旋转调强优化方法及放射治疗设备，本发明与现有技术相比的优点在于：本发明结合了IMRT的Sliding Window的快速优化和肿瘤形态学特征算法的优点，又有效的克服了两者的缺点。针对旋转射野，采用基于射野强度分布和子野分割后的结果，转化成动态再使用共轭梯度法优化子野权重，扰动技术优化子野形状。这样既可以满足计算时间的要求，又能满足优化效果要求。

Description

一种基于动态光栅的旋转调强优化方法及放射治疗设备
技术领域
本发明涉及加速器治疗装置，具体涉及一种基于动态光栅的旋转调强优化方法及放射治疗设备。
背景技术
放射治疗作为肿瘤的一种局部治疗手段，一直在不断寻求解决的一个基本问题是，如何较好地处理肿瘤组织和周围正常组织的剂量关系，使肿瘤得到最大限度的局部控制而周围正常组织和器官的放射损伤最小。临床经验证明，肿瘤的局部控制与正常组织的放射损伤有一定的关系。多数情况下，肿瘤控制率与正常组织损伤成正比例，即提高肿瘤的局部控制率必然造成对正常组织的更多的损伤。临床经验也证明，通过改进照射技术、选择合理的时间剂量因子，在保持同等水平的肿瘤控制率的情况下，可减低正常组织的放射损伤。治疗方案的优化是实现上述目的的途径之一。
旋转调强是目前发展的最先进的调强技术，能够很好的避免传统IMRT治疗时间长的缺点，同时治疗效果优于IMRT，大大提高计划执行的效率。传统的IMRT治疗，在治疗过程中，加速器大机架是静止的，每个治疗野切换时，都要停止出束，机架到位之后再重新出束，导致治疗效率相对旋转治疗较为低下；最初的旋转治疗没有优化过程，将一个弧划分成若干个方向，每个方向上的子野形状根据BEV(Beam Eye View射野方向视图)中的形状确定。后来发展了旋转调强优化，采用一些优化方法，优化得到每个方向上的子野形状和权重。但是目前旋转调强优化的难点在于优化参数众多，优化时间长，很难得到最优的方案。
发明内容
为了克服现有技术的不足，解决现有技术中优化参数众多，优化时间长，很难得到最优的方案的缺陷，本发明提供了一种基于动态光栅的旋转调强优化方法和放射治疗设备。能够快速解决旋转调强参数优化问题，其优化结果能够完全满足放疗计划制定的要求，而又使得计算时间能够在本领域技术人员可以接受的范围内。
本文中的Arc即弧形野。
为达到上述目的，本发明的技术方案是：
一种基于动态光栅的旋转调强优化方法，包括以下步骤：
第1步，设置以下参数：弧形野数量、起始角度、角度间隔、剂量约束；
第2步，根据上述参数将Arc离散化，获得若干静态临时射野Beam；
第3步，基于第2步，对上述静态临时射野Beam进行强度优化，再对强度进行子野分割，将得到的静态临时射野Beam的子野按顺序均匀分布到原来的弧形野上；
第4步，使用gDPM计算每个子野的剂量分布，最后再使用共轭梯度法优化每个子野的权重，得到最优的剂量分布。
本发明与现有技术相比的优点在于：本发明结合了IMRT的Sliding Window的快速优化和肿瘤形态学特征算法的优点，又有效的克服了两者的缺点。针对旋转射野，采用基于射野强度分布和子野分割后的结果，转化成动态再使用共轭梯度法优化子野权重，扰动技术优化子野形状。这样既可以满足计算时间的要求，又能满足优化效果要求。
进一步的，所述第2步具体如下：
根据Arc的起始角度G_start和结束角度G_end，角度间隔G_step、每个射野子野数SegNumOneBeam，计算：
子野个数：
NSeg＝(G_end-G_start)/G_step；
射野个数：
NBeam＝NSeg/SegNumOneBeam；
射野角度间隔：
Step＝(G_end-G_start)/NBeam；
射野角度：
Ang(i)＝G_start+i*Step-Step/2进一步的，所述第3步中的强度优化为采用共轭梯度法对其进行强度优化，具体如下：
设定目标函数
Figure GDA0002447946740000031
其中，
Figure GDA0002447946740000032
Figure GDA0002447946740000033
Figure GDA0002447946740000034
其中，
Figure GDA0002447946740000035
为与所有OAR有关的目标函数项；
Figure GDA0002447946740000036
为与靶区有关的目标函数项；NOAR为OAR的个数；NTi为第i个OAR中的采样点数；NPTV为靶区中的采样点数；当器官中的点剂量违反约束时，δj＝1，否则δj＝0；Wj为第j个点的权重；dj为第j点的计算剂量；Pj为第j点的处方剂量；a、β为控制靶区和OAR的重要性的归一化因子；Nray为射束的总个数；ajm为单位强度的第m个射束对第j个点的剂量贡献；
Figure GDA0002447946740000041
为第m个射束的强度。
进一步的，所述第3步中的对强度进行子野分割，具体如下：
对输入的强度图，采用SlidingWindow的分割算法进行分割，得到连续运动的子野Seg1。
进一步的，所述第3步中的对强度进行子野分割时，每个临时射野Beam的子野数量应保持一致，记为m。
进一步的，对于上述临时射野Beam，将其分割得到的SlidingWindow的Seg1，均匀的分布到当前角度的(-step/2，step/2)；则整个Arc的m*n个子野的覆盖范围为：G_end-G_start。
进一步的，所述第4步中，使用gDPM计算每个子野的剂量分布，具体如下：
对每一个子野，使用基于GPU的蒙特卡洛算法gDPM，计算剂量分布Dseg(i)。
进一步的，所述第4步中，最后再使用共轭梯度法优化每个子野的权重，具体如下：
设定目标函数：
F′obj(Wi)＝F′target(Wi)+F′oar(Wi)；其中，
Figure GDA0002447946740000042
Figure GDA0002447946740000043
Figure GDA0002447946740000044
其中，
F′oar(Wi)为所有与OAR有关的目标函数项，
F′target(Wi)则为与靶区有关的目标函数项，
Noar为OAR的个数，Ntarget为靶区的个数。
本发明还提供一种放射治疗设备，其特征在于，包括：
参数设置模块，用以设置以下参数：弧形野数量、起始角度、剂量约束；
Arc离散化模块，用于根据上述参数将Arc离散化，获得若干静态临时射野Beam；
强度优化模块，用于对上述静态临时射野Beam进行强度优化，再对强度进行子野分割，将得到的静态临时射野Beam的子野按顺序均匀分布到原来的弧形野上；
剂量分布优化模块，使用gDPM计算每个子野的剂量分布，最后再使用共轭梯度法优化每个子野的权重，得到最优的剂量分布。
进一步的，所述Arc离散化模块具体工作原理如下：
根据Arc的起始角度G_start和结束角度G_end，角度间隔G_step、每个射野子野数SegNumOneBeam，计算：
子野个数：
NSeg＝(G_end-G_start)/G_step；
射野个数：
NBeam＝NSeg/SegNumOneBeam；
射野角度间隔：
Step＝(G_end-G_start)/NBeam；
射野角度：
Ang(i)＝G_start+i*Step-Step/2；
所述强度优化模块采用共轭梯度法对其进行强度优化，具体如下：
设定目标函数
Figure GDA0002447946740000061
其中，
Figure GDA0002447946740000062
Figure GDA0002447946740000063
Figure GDA0002447946740000064
其中，
Figure GDA0002447946740000065
为与所有OAR有关的目标函数项；
Figure GDA0002447946740000066
为与靶区有关的目标函数项；NOAR为OAR的个数；NTi为第i个OAR中的采样点数；NPTV为靶区中的采样点数；当器官中的点剂量违反约束时，δj＝1，否则δj＝0；Wj为第j个点的权重；dj为第j点的计算剂量；Pj为第j点的处方剂量；α、β为控制靶区和OAR的重要性的归一化因子；Nray为射束的总个数；ajm为单位强度的第m个射束对第j个点的剂量贡献；
Figure GDA0002447946740000067
为第m个射束的强度；
所述强度优化模块对强度进行子野分割，具体如下：
对输入的强度图，采用SlidingWindow的分割算法进行分割，得到连续运动的子野Seg1；并对强度进行子野分割时，每个临时射野Beam的子野数量应保持一致，记为m；并将其分割得到的SlidingWindow的Seg1，均匀的分布到当前角度的(-step/2，step/2)；则整个Arc的m*n个子野的覆盖范围为：G_end-G_start。所述剂量分布优化模块，使用gDPM计算每个子野的剂量分布，具体如下：
对每一个子野，使用基于GPU的蒙特卡洛算法gDPM，计算剂量分布Dseg(i)，最后再使用共轭梯度法优化每个子野的权重，具体如下：
设定目标函数：
F′obj(Wi)＝F′target(Wi)+F′oar(Wi)；其中，
Figure GDA0002447946740000071
Figure GDA0002447946740000072
Figure GDA0002447946740000073
其中，
F′oar(Wi)为所有与OAR有关的目标函数项，
F′target(Wi)则为与靶区有关的目标函数项，
Noar为OAR的个数，Ntarget为靶区的个数。
附图说明
图1是剂量点与射野笔形束之间的关系图；
图2是本发明的过程示意图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
在实际射疗的过程中，医生需要先根据患者病症的具体情况制定治疗计划，这个治疗计划包括射线从哪几个角度进行照射，每个角度的射野形状及照射时间(剂量控制)，照射次数等等。例如，医生为某一病患给出以下治疗计划：计划进行2个Arc射线，一个从90度开始至270度结束，一个从270度开始至90度结束；给定靶区最小剂量为6000cGy，最大65cGy，危及器官剂量不超过临床需求。最初，主机架旋转到某个角度后，静止下来，完成该角度的照射剂量，再旋转到下一个角度，静止并完成对应的剂量，以此类推，即静态调强。现有技术中，在每一个角度之间，开始谋求连续的调强治疗，即动态调强，在从一个角度到另一个角度之间，进行连续的调强，难度很大，现有技术还存在很多问题，如优化参数多，优化时间过长及优化的结果不是太理想等的。
本发明提出一种解决方案，基于动态光栅的旋转调强优化方法，包括以下步骤：
第1步，设置以下参数：弧形野数量、起始角度、角度间隔、剂量约束；
第2步，根据上述参数将Arc离散化，获得若干静态临时射野Beam；
第3步，基于第2步，对上述静态临时射野Beam进行强度优化，再对强度进行子野分割，将得到的静态临时射野Beam的子野按顺序均匀分布到原来的弧形野上；
第4步，使用gDPM计算每个子野的剂量分布，最后再使用共轭梯度法优化每个子野的权重，得到最优的剂量分布。
本发明与现有技术相比的优点在于：本发明结合了IMRT的Sliding Window的快速优化和肿瘤形态学特征算法的优点，又有效的克服了两者的缺点。针对旋转射野，采用基于射野强度分布和子野分割后的结果，转化成动态再使用共轭梯度法优化子野权重，扰动技术优化子野形状。这样既可以满足计算时间的要求，又能满足优化效果要求。
在实际应用中，所述第2步具体如下：
根据Arc的起始角度G_start和结束角度G_end，角度间隔G_step、每个射野子野数SegNumOneBeam，计算：
子野个数：
NSeg＝(G_end-G_start)/G_step；
射野个数：
NBeam＝NSeg/SegNumOneBeam；
射野角度间隔：
Step＝(G_end-G_start)/NBeam；
射野角度：
Ang(i)＝G_start+i*Step-Step/2。
在实际应用中，所述第3步中的强度优化为采用共轭梯度法对其进行强度优化，具体如下：
设定目标函数
Figure GDA0002447946740000091
其中，
Figure GDA0002447946740000092
Figure GDA0002447946740000093
Figure GDA0002447946740000094
其中，
Figure GDA0002447946740000095
为与所有OAR有关的目标函数项；
Figure GDA0002447946740000096
为与靶区有关的目标函数项；NOAR为OAR的个数；NTi为第i个OAR中的采样点数；NPTV为靶区中的采样点数；当器官中的点剂量违反约束时，δj＝1，否则δj＝0；Wj为第j个点的权重；dj为第j点的计算剂量；Pj为第j点的处方剂量；a、β为控制靶区和OAR的重要性的归一化因子；Nray为射束的总个数(ray又称为笔形束pencil beam)；ajm为单位强度的第m个射束对第j个点的剂量贡献，整个ajm构成的矩阵称为能量特征矩阵(或称为沉积矩阵deposition matrix)；
Figure GDA0002447946740000097
为第m个射束的强度。剂量点与射野笔形束的关系如图1所示。图1中，对应于j点(图1中标号为1的点)的为第m个强度为xm的笔形束(图1中标号为2)，对应于j点的点剂量值由以下公式决定：
Figure GDA0002447946740000101
使用CG方法优化出射野的强度分布，记为I(i)。
在实际应用中，所述第3步中的对强度进行子野分割，具体如下：
对输入的强度图，采用SlidingWindow的分割算法进行分割，得到连续运动的子野Seg1。
在实际应用中，所述第3步中的对强度进行子野分割时，每个临时射野Beam的子野数量应保持一致，记为m。
在实际应用中，对于上述临时射野Beam，将其分割得到的SlidingWindow的Seg1，均匀的分布到当前角度的(-step/2，step/2)；则整个Arc的m*n个子野的覆盖范围为：G_end-G_start。在实际应用中，所述第4步中，使用gDPM计算每个子野的剂量分布，具体如下：
对每一个子野，使用基于GPU的蒙特卡洛算法gDPM，计算剂量分布Dseg(i)。
在实际应用中，所述第4步中，最后再使用共轭梯度法优化每个子野的权重，具体如下：
设定目标函数：
F′obj(Wi)＝F′taregt(Wi)+F′oar(Wi)；其中，
Figure GDA0002447946740000111
Figure GDA0002447946740000112
Figure GDA0002447946740000113
其中，
F′oar(Wi)为所有与OAR有关的目标函数项，
F′target(Wi)则为与靶区有关的目标函数项，
Noar为OAR的个数，Ntarget为靶区的个数。
例如，如图2所示，开机，得到初始化的系统，医生根据治疗计划设置Arc数量、其实角度、各器官的剂量约束等参数。
然后装置将Arc离散化获得临时射野Beam(i)。
进而按照上述方法对其进行强度分布优化。
再按照上述方法进行子野分割，并保证每个临时射野Beam的子野数目一致。
根据实际需要，完成子野转换，可将每个临时射野Beam的子野在其角度范围内均匀分布。目的是保证每个子野覆盖的范围一致。
对子野剂量进行计算。
完成子野权重优化，保证剂量分布最优，获得最佳的执行治疗方案(旋转动态调强方案)。
本发明还提供一种放射治疗设备，其特征在于，包括：
参数设置模块，用以设置以下参数：弧形野数量、起始角度、剂量约束；
Arc离散化模块，用于根据上述参数将Arc离散化，获得若干静态临时射野Beam；
强度优化模块，用于对上述静态临时射野Beam进行强度优化，再对强度进行子野分割，将得到的静态临时射野Beam的子野按顺序均匀分布到原来的弧形野上；
剂量分布优化模块，使用gDPM计算每个子野的剂量分布，最后再使用共轭梯度法优化每个子野的权重，得到最优的剂量分布。
进一步的，所述Arc离散化模块具体工作原理如下：
根据Arc的起始角度G_start和结束角度G_end，角度间隔G_step、每个射野子野数SegNumOneBeam，计算：
子野个数：
NSeg＝(G_end-G_start)/G_step；
射野个数：
NBeam＝NSeg/SegNumOneBeam；
射野角度间隔：
Step＝(G_end-G_start)/NBeam；
射野角度：
Ang(i)＝G_start+i*Step-Step/2；
所述强度优化模块采用共轭梯度法对其进行强度优化，具体如下：
设定目标函数
Figure GDA0002447946740000121
其中，
Figure GDA0002447946740000122
Figure GDA0002447946740000131
Figure GDA0002447946740000132
其中，
Figure GDA0002447946740000133
为与所有OAR有关的目标函数项；
Figure GDA0002447946740000134
为与靶区有关的目标函数项；NOAR为OAR的个数；NTi为第i个OAR中的采样点数；NPTV为靶区中的采样点数；当器官中的点剂量违反约束时，δj＝1，否则δj＝0；Wj为第j个点的权重；dj为第j点的计算剂量；Pj为第j点的处方剂量；α、β为控制靶区和OAR的重要性的归一化因子；Nray为射束的总个数；ajm为单位强度的第m个射束对第j个点的剂量贡献；
Figure GDA0002447946740000135
为第m个射束的强度；
所述强度优化模块对强度进行子野分割，具体如下：
对输入的强度图，采用SlidingWindow的分割算法进行分割，得到连续运动的子野Segl；并对强度进行子野分割时，每个临时射野Beam的子野数量应保持一致，记为m；并将其分割得到的SlidingWindow的Seg1，均匀的分布到当前角度的(-step/2，step/2)；则整个Arc的m*n个子野的覆盖范围为：G_end-G_start。所述剂量分布优化模块，使用gDPM计算每个子野的剂量分布，具体如下：
对每一个子野，使用基于GPU的蒙特卡洛算法gDPM，计算剂量分布Dseg(i)，最后再使用共轭梯度法优化每个子野的权重，具体如下：
设定目标函数：
F′obj(Wi)＝F′taregt(Wi)+F′oar(Wi)；其中，
Figure GDA0002447946740000141
Figure GDA0002447946740000142
Figure GDA0002447946740000143
其中，
F′oar(Wi)为所有与OAR有关的目标函数项，
F′target(Wi)则为与靶区有关的目标函数项，
Noar为OAR的个数，Ntarget为靶区的个数。
本发明与现有技术相比的优点在于：本发明结合了IMRT的Sliding Window的快速优化和肿瘤形态学特征算法的优点，又有效的克服了两者的缺点。针对旋转射野，采用基于射野强度分布和子野分割后的结果，转化成动态再使用共轭梯度法优化子野权重，扰动技术优化子野形状。这样既可以满足计算时间的要求，又能满足优化效果要求。
以上所述的仅是本发明的优选实施方式，应当指出，对于本领域的普通技术人员来说，在不脱离本发明创造构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。

Claims (5)

1.一种基于动态光栅的旋转调强优化方法，其特征在于，包括以下步骤：
第1步，设置以下参数：弧形野数量、起始角度、角度间隔、剂量约束；
第2步，根据上述参数将Arc离散化，获得若干静态临时射野Beam；
第3步，基于第2步，对上述静态临时射野Beam进行强度优化，再对强度进行子野分割，将得到的静态临时射野Beam的子野按顺序均匀分布到原来的弧形野上，所述第3步中的强度优化为采用共轭梯度法对其进行强度优化，具体如下：
设定目标函数
Figure FDA0002447946730000011
其中
Figure FDA0002447946730000012
Figure FDA0002447946730000013
Figure FDA0002447946730000014
其中，
Figure FDA0002447946730000015
为与所有OAR有关的目标函数项；
Figure FDA0002447946730000016
为与靶区有关的目标函数项；NOAR为OAR的个数；NTi为第i个OAR中的采样点数；NPTV为靶区中的采样点数；当器官中的点剂量违反约束时，δj＝1，否则δj＝0；Wj为第j个点的权重；dj为第j点的计算剂量；Pj为第j点的处方剂量；α、β为控制靶区和OAR的重要性的归一化因子；Nray为射束的总个数；αjm为单位强度的第m个射束对第j个点的剂量贡献；
Figure FDA0002447946730000021
为第m个射束的强度；
所述第3步中的对强度进行子野分割，具体如下：对输入的强度图，采用SlidingWindow的分割算法进行分割，得到连续运动的子野Seg1；
第4步，使用gDPM计算每个子野的剂量分布，最后再使用共轭梯度法优化每个子野的权重，得到最优的剂量分布；
使用gDPM计算每个子野的剂量分布，具体如下：
对每一个子野，使用基于GPU的蒙特卡洛算法gDPM，计算剂量分布Dseg(i)；
所述第4步中，最后再使用共轭梯度法优化每个子野的权重，具体如下：
设定目标函数：
F′obj(Wi)＝F′target(Wi)+F′oar(Wi)；其中，
Figure FDA0002447946730000022
Figure FDA0002447946730000023
Figure FDA0002447946730000024
其中，
F′oar(Wi)为所有与OAR有关的目标函数项，
F′target(Wi)则为与靶区有关的目标函数项，
Noar为OAR的个数，Ntarget为靶区的个数。
2.根据权利要求1所述的基于动态光栅的旋转调强优化方法，其特征在于，所述第2步具体如下：
根据Arc的起始角度G_start和结束角度G_end，角度间隔G_step、每个射野子野数SegNumOneBeam，计算：
子野个数：
NSeg＝(G_end-G_start)/G_step；
射野个数：
NBeam＝NSeg/SegNumOneBeam；
射野角度间隔：
Step＝(G_end-G_start)/NBeam；
射野角度：
Ang(i)＝G_start+i*Step-Step/2。
3.根据权利要求1所述的基于动态光栅的旋转调强优化方法，其特征在于，所述第3步中的对强度进行子野分割时，每个临时射野Beam的子野数量应保持一致，记为m。
4.根据权利要求3所述的基于动态光栅的旋转调强优化方法，其特征在于，对于上述临时射野Beam，将其分割得到的SlidingWindow的Seg1，均匀的分布到当前角度的(-step/2，step/2)；则整个Arc的m*n个子野的覆盖范围为：G_end-G_start。
5.一种放射治疗设备，其特征在于，包括：
参数设置模块，用以设置以下参数：弧形野数量、起始角度、剂量约束；
Arc离散化模块，用于根据上述参数将Arc离散化，获得若干静态临时射野Beam；
强度优化模块，用于对上述静态临时射野Beam进行强度优化，再对强度进行子野分割，将得到的静态临时射野Beam的子野按顺序均匀分布到原来的弧形野上；
剂量分布优化模块，使用gDPM计算每个子野的剂量分布，最后再使用共轭梯度法优化每个子野的权重，得到最优的剂量分布；
所述Arc离散化模块具体工作原理如下：
根据Arc的起始角度G_start和结束角度G_end，角度间隔G_step、每个射野子野数SegNumOneBeam，计算：
子野个数：
NSeg＝(G_end-G_start)/G_step；
射野个数：
NBeam＝NSeg/SegNumOneBeam：
射野角度间隔：
Step＝(G_end-G_start)/NBeam；
射野角度：
Ang(i)＝G_start+i*Step-Step/2；
所述强度优化模块采用共轭梯度法对其进行强度优化，具体如下：
设定目标函数
Figure FDA0002447946730000041
其中，
Figure FDA0002447946730000042
Figure FDA0002447946730000043
Figure FDA0002447946730000044
其中，
Figure FDA0002447946730000045
为与所有OAR有关的目标函数项；
Figure FDA0002447946730000046
为与靶区有关的目标函数项；NOAR为OAR的个数；NTi为第i个OAR中的采样点数；NPTV为靶区中的采样点数；当器官中的点剂量违反约束时，δj＝1，否则δj＝0；Wj为第j个点的权重；dj为第j点的计算剂量；Pj为第j点的处方剂量；α、β为控制靶区和OAR的重要性的归一化因子；Nray为射束的总个数；αjm为单位强度的第m个射束对第j个点的剂量贡献；
Figure FDA0002447946730000051
为第m个射束的强度；
所述强度优化模块对强度进行子野分割，具体如下：
对输入的强度图，采用SlidingWindow的分割算法进行分割，得到连续运动的子野Seg1；并对强度进行子野分割时，每个临时射野Beam的子野数量应保持一致，记为m；并将其分割得到的SlidingWindow的Seg1，均匀的分布到当前角度的(-step/2，step/2)；则整个Arc的m*n个子野的覆盖范围为：G_end-G_start；
所述剂量分布优化模块，使用gDPM计算每个子野的剂量分布，具体如下：
对每一个子野，使用基于GPU的蒙特卡洛算法gDPM，计算剂量分布Dseg(i)，最后再使用共轭梯度法优化每个子野的权重，具体如下：
设定目标函数：
F′obj(Wi)＝F′target(Wi)+F′oar(Wi)；其中，
Figure FDA0002447946730000052
Figure FDA0002447946730000053
Figure FDA0002447946730000054
其中，
F′oar(Wi)为所有与OAR有关的目标函数项，
F′target(Wi)则为与靶区有关的目标函数项，
Noar为OAR的个数，Ntarget为靶区的个数。
CN201711498167.5A 2017-12-29 2017-12-29 一种基于动态光栅的旋转调强优化方法及放射治疗设备 Active CN108175953B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711498167.5A CN108175953B (zh) 2017-12-29 2017-12-29 一种基于动态光栅的旋转调强优化方法及放射治疗设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711498167.5A CN108175953B (zh) 2017-12-29 2017-12-29 一种基于动态光栅的旋转调强优化方法及放射治疗设备

Publications (2)

Publication Number Publication Date
CN108175953A CN108175953A (zh) 2018-06-19
CN108175953B true CN108175953B (zh) 2020-07-21

Family

ID=62549684

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711498167.5A Active CN108175953B (zh) 2017-12-29 2017-12-29 一种基于动态光栅的旋转调强优化方法及放射治疗设备

Country Status (1)

Country Link
CN (1) CN108175953B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109499011B (zh) * 2018-12-21 2020-11-10 苏州雷泰医疗科技有限公司 用于正交双层光栅装置的动态调强子野分割方法
CN110170109B (zh) * 2018-12-21 2020-12-18 苏州雷泰医疗科技有限公司 用于放疗设备的正交双层光栅装置及其子野分割控制方法
CN109499012B (zh) * 2018-12-27 2021-05-11 苏州雷泰医疗科技有限公司 一种优化剂量引导的tps自动迭代优化方法
CN110215623B (zh) * 2019-06-11 2020-11-10 苏州雷泰智能科技有限公司 基于子野优化的正交双层光栅旋转调强的实现方法及装置
CN110404185A (zh) * 2019-07-09 2019-11-05 苏州雷泰智能科技有限公司 一种正交双层光栅的旋转调强分割方法及装置
CN111681736B (zh) * 2020-06-08 2023-05-12 苏州雷泰医疗科技有限公司 一种基于分象限的正交双层光栅动态调强分割方法
CN113521561A (zh) * 2021-06-22 2021-10-22 苏州雷泰医疗科技有限公司 正交双层光栅逆向旋转调强的实现方法、存储介质及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101636199A (zh) * 2007-03-19 2010-01-27 皇家飞利浦电子股份有限公司 治疗的优化
CN101848745A (zh) * 2007-10-16 2010-09-29 皇家飞利浦电子股份有限公司 用于调强弧形治疗序列化与最优化的方法和装置
CN103083821A (zh) * 2013-01-27 2013-05-08 合肥超安医疗科技有限公司 一种直接优化静态调强子野形状及其机器跳数的系统和方法
CN104318122A (zh) * 2014-11-06 2015-01-28 中国科学院合肥物质科学研究院 一种基于扰动的旋转调强优化方法
CN104338240A (zh) * 2014-10-31 2015-02-11 章桦 一种在线自适应放疗计划自动优化方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101636199A (zh) * 2007-03-19 2010-01-27 皇家飞利浦电子股份有限公司 治疗的优化
CN101848745A (zh) * 2007-10-16 2010-09-29 皇家飞利浦电子股份有限公司 用于调强弧形治疗序列化与最优化的方法和装置
CN103083821A (zh) * 2013-01-27 2013-05-08 合肥超安医疗科技有限公司 一种直接优化静态调强子野形状及其机器跳数的系统和方法
CN104338240A (zh) * 2014-10-31 2015-02-11 章桦 一种在线自适应放疗计划自动优化方法及装置
CN104318122A (zh) * 2014-11-06 2015-01-28 中国科学院合肥物质科学研究院 一种基于扰动的旋转调强优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于GPU的蒙特卡洛模拟在辐射剂量计算中的研究进展及应用;黄莹;《中国医学物理杂志》;20171031;第34卷(第10期);全文 *
容积弧形调强放疗的技术分析;冯丽娜;《硕士学位论文》;20131231;全文 *
调强放射治疗中的优化技术研究;李永杰;《博士学位论文》;20041231;正文第3章 *

Also Published As

Publication number Publication date
CN108175953A (zh) 2018-06-19

Similar Documents

Publication Publication Date Title
CN108175953B (zh) 一种基于动态光栅的旋转调强优化方法及放射治疗设备
CN108211136B (zh) 基于动态光栅的旋转调强优化方法及放射治疗设备
Keller-Reichenbecher et al. Intensity modulation with the “step and shoot” technique using a commercial MLC: A planning study
Tseung et al. Clinically applicable Monte Carlo–based biological dose optimization for the treatment of head and neck cancers with spot-scanning proton therapy
RU2605998C2 (ru) Одновременная оптимизация мультимодальной инверсии для планирования радиотерапевтического лечения
JP2019532787A (ja) 放射線療法の治療計画を選択するシステムおよび方法
Cagni et al. Knowledge-based treatment planning: An inter-technique and inter-system feasibility study for prostate cancer
Longobardi et al. Comparing 3DCRT and inversely optimized IMRT planning for head and neck cancer: equivalence between step-and-shoot and sliding window techniques
CN107073284A (zh) 用于辐射逆向治疗计划的系统和计算机程序产品
US9495513B2 (en) GPU-based fast dose calculator for cancer therapy
Mišić et al. Neighborhood search approaches to non-coplanar beam orientation optimization for total marrow irradiation using IMRT
WO2020248489A1 (zh) 基于子野优化的正交双层光栅旋转调强的实现方法及装置
Onizuka et al. Monte Carlo dose verification of VMAT treatment plans using Elekta Agility 160-leaf MLC
Li et al. Linear energy transfer incorporated spot-scanning proton arc therapy optimization: a feasibility study
Balvert et al. A framework for inverse planning of beam-on times for 3D small animal radiotherapy using interactive multi-objective optimisation
Song et al. Fully automatic volumetric modulated arc therapy plan generation for rectal cancer
James et al. Considerations when treating lung cancer with passive scatter or active scanning proton therapy
Kan et al. The performance of the progressive resolution optimizer (PRO) for RapidArc planning in targets with low‐density media
US11590362B2 (en) Radiotherapy treatment planning based on treatment delivery efficiency
JP7423628B2 (ja) パッシブイオン放射線治療計画作成及び送達のためのシステム及び方法
CN107998519A (zh) 一种用于imrt的分子动力学再优化算法
Al-Yahya et al. Monte Carlo based modulated electron beam treatment planning using a few-leaf electron collimator—feasibility study
Kalet et al. Dosimetric comparison of single-beam multi-arc and 2-beam multi-arc VMAT optimization in the Monaco treatment planning system
CN112883516A (zh) 一种基于叶片速度时间曲线的旋转调强优化装置
CN113797447A (zh) 放射治疗系统及其治疗计划生成方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant