CN111785378A - 放疗计划系统和存储介质 - Google Patents

放疗计划系统和存储介质 Download PDF

Info

Publication number
CN111785378A
CN111785378A CN202010568528.4A CN202010568528A CN111785378A CN 111785378 A CN111785378 A CN 111785378A CN 202010568528 A CN202010568528 A CN 202010568528A CN 111785378 A CN111785378 A CN 111785378A
Authority
CN
China
Prior art keywords
rotating
layer
target
target area
dose distribution
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.)
Granted
Application number
CN202010568528.4A
Other languages
English (en)
Other versions
CN111785378B (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.)
Xinlicheng Medical Accelerator Wuxi Co ltd
Original Assignee
Xinlicheng Medical Accelerator Wuxi 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 Xinlicheng Medical Accelerator Wuxi Co ltd filed Critical Xinlicheng Medical Accelerator Wuxi Co ltd
Priority to CN202010568528.4A priority Critical patent/CN111785378B/zh
Publication of CN111785378A publication Critical patent/CN111785378A/zh
Application granted granted Critical
Publication of CN111785378B publication Critical patent/CN111785378B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1092Details

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Epidemiology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Primary Health Care (AREA)
  • Business, Economics & Management (AREA)
  • General Business, Economics & Management (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

本发明实施例公开了一种放疗计划系统和存储介质,该系统包括处理器:所述处理器在工作时实现以下步骤:根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,所述目标靶区由至少两个旋转层同心层叠而成,所述旋转层的截面为圆环或部分圆环,所述圆环或部分圆环的圆心位于所述目标靶区的中心轴上;根据所述目标靶区的处方剂量分布确定所述目标靶区的每个旋转层的目标剂量分布;根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成所述目标靶区的放疗计划,其中,所述计划参数包括束流参数和旋转层的叶片对匹配信息。解决了现有放射治疗方法很难兼顾靶区适形度和剂量利用率的问题。

Description

放疗计划系统和存储介质
技术领域
本发明实施例涉及计算机软件领域,尤其涉及一种放疗计划系统和存储介质。
背景技术
肿瘤放射治疗一般分为多次实施(大约35次左右),为使被辐照的健康组织得到恢复,整个疗期可能持续数周。如果通过更精确的控制,使得健康组织的被辐照的剂量比率更少,则每次对肿瘤组织施加更多剂量就成为可能,这样总的治疗次数将会减少,甚至最少一次就完成整个疗期的剂量,这种方案被称为大剂量分割(Hypofractionation)方案。
目前,常用的放射治疗方法是IMRT(intensity-modulated radiation therapy,简称IMRT,强度调制放射治疗,即调强放疗),是由计算机控制的X光加速器向肿瘤内的特定区域发射精确的辐射剂量。其从诞生至今已经过多次革新,最开始的静态调强(Staticintensity-modulated radiation therapy,简称S-IMRT)技术,然后是动态调强(Dynamicintensity-modulated radiation therapy,简称D-IMRT)技术。近些年,又发展出两种更先进的动态调强技术:一种称为容积旋转调强放射治疗(Volumetric Modulated ArcTherapy,简称VMAT)技术,另一种称为螺旋断层放射治疗(Helical Tomotherapy)技术。
其中,螺旋断层放射治疗方法,从剂量适形度上来看无疑是一种很优秀的方法,但这种方法的剂量利用率非常低,大约为1%左右。容积旋转调强放射治疗相较于螺旋断层放射治疗,剂量利用效率高一些,相对治疗时会更短一些,但仍很难满足大剂量分割方案的要求。
综上,现有放射治疗方法存在很难兼顾靶区适形度和剂量利用率的问题。
发明内容
本发明实施例提供了一种放疗计划系统和存储介质,解决了现有放射治疗方法很难兼顾靶区适形度和剂量利用率的问题。
第一方面,本发明实施例提供了一种放疗计划系统,该系统包括处理器:所述处理器在工作时实现以下步骤:
根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,所述目标靶区由至少两个旋转层同心层叠而成,所述旋转层的截面为圆环或部分圆环,所述圆环或部分圆环的圆心位于所述目标靶区的中心轴上;
根据所述目标靶区的处方剂量分布确定所述目标靶区的每个旋转层的目标剂量分布;
根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成所述目标靶区的放疗计划,其中,所述计划参数包括束流参数和旋转层的叶片对匹配信息。
第二方面,本发明实施例还提供了一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行如下方法:
根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,所述目标靶区由至少两个旋转层同心层叠而成,所述旋转层的截面为圆环或部分圆环,所述圆环或部分圆环的圆心位于所述目标靶区的中心轴上;
根据所述目标靶区的处方剂量分布确定所述目标靶区的每个旋转层的目标剂量分布;
根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成所述目标靶区的放疗计划,其中,所述计划参数包括束流参数和旋转层的叶片对匹配信息。
本发明实施例提供的放射治疗系统的技术方案,相较于现有技术,将不规则的目标靶区分解成若干形状规则且同心层叠旋转层,且旋转层的截面为环形和/或部分圆环形,以及根据目标靶区的处方剂量分布确定每个旋转层的目标剂量分布,这样即可根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,从而确定目标靶区的放疗计划。由于任何形状的目标靶区的截面均可由同心层叠的环形和/或部分圆环形拼凑成,那么在环形厚度较小时,各个旋转层的剂量分布具有较高的均匀性,从而使得目标靶区的各个评价指标数据都比较好。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图做一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的由放疗计划系统执行的放疗方法的流程图;
图2是本发明实施例一提供的旋转层形成示意图;
图3A是本发明实施例一提供的单个旋转层的形成示意图;
图3B是本发明实施例一提供的又一单个旋转层的形成示意图;
图4是本发明实施例一提供的叶片对与旋转层的对应关系图;
图5是本发明实施例一提供的包含10个旋转层的目标靶区截面示意图;
图6是本发明实施例一提供的不规则目标靶区的旋转层分布示意图;
图7是本发明实施例一提供的哑铃状目标靶区的旋转层分布图;
图8A是本发明实施例一提供的极坐标系下的剂量分布示意图;
图8B是本发明实施例一提供的直角坐标系下的剂量分布示意图;
图9是本发明实施例一提供的同时向不同分区的同一旋转层输出射束的示意图;
图10是本发明实施例二提供的由放疗计划系统执行的放疗方法的流程图;
图11A是本发明实施例二提供的旋转层之间的剂量影响示意图;
图11B是本发明实施例二提供的又一旋转层之间的剂量影响示意图;
图12是本发明实施例三提供的由放疗计划系统执行的放疗方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下将参照本发明实施例中的附图,通过实施方式清楚、完整地描述本发明的技术方案,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
本发明实施例一提供了一种放疗计划系统,该系统包括处理器和存储有计算机程序的存储器:该处理器在执行该计算机程序时,实现图1所示的放疗方法的步骤,该放疗方法适用于基于多圈分层方法制定放疗计划,具体如下:
S101、根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,目标靶区由至少两个旋转层同心层叠而成,旋转层的截面为圆环或部分圆环,圆环或部分圆环的圆心位于目标靶区的中心轴上。
其中,目标靶区是由医生勾画的治疗靶区,即束流辐照的目标区域。束流是指放疗设备的束流输出装置向目标靶区输出的射线束,比如X射线束等。
其中,径向为目标靶区截面的半径方向。径向厚度可基于放疗设备的多叶光栅装置的叶片的厚度或等效厚度确定。本实施例各个旋转层的径向厚度优选均为一个叶片的等效厚度。其中,叶片的等效厚度为叶片能够用于约束束流形状的厚度。
多叶光栅装置通过每个叶片对中的两个叶片间的缝隙长度和叶片厚度来约束射束形状。如图2、图3A、图3B和图4所示,处于打开状态的叶片对输出的束流截面为矩形,该矩形的长边等于叶片对的缝隙宽度,对应于旋转层在目标靶区中心轴方向上的长度;该矩形的宽边为叶片对的厚度或等效厚度,对应叶片对输出的束流在目标靶区上形成的高剂量区的宽度,即旋转层的厚度。可以理解的是,如果不同叶片对在绕目标靶区中心轴的旋转过程中保持位置不变,并分别向目标靶区输出束流,那么不同叶片对所输出的束流会在目标靶区中形成半径不同的空心柱形高剂量区。图3A和图3B示出了一个叶片对绕目标靶区旋转形成旋转层的过程。其中,图3A中的环状半径比图3B中的环状半径大,但不同半径的旋转层的圆心均在目标靶区的中心轴上。其中,目标靶区的中心轴即为束流输出装置的旋转轴,其方向为垂直于束流输出装置的旋转面,并经过放疗设备的等中心。
由于束流输出装置的旋转轴与旋转层的中心轴是相同的,因此叶片对与旋转层的对应关系是固定不变的。图4和图5示出了叶片对与旋转层的对应关系,叶片对A1/B1始终对应旋转层1,叶片对A2/B2始终对应旋转层2,叶片对A3/B3始终对应旋转层3。图5示出了由10个叶片对形成的10个旋转层的分布。
可以理解的是,如果目标靶区的截面不是圆形,那么目标靶区的截面包括位于中心位置的多个同心圆环,以及位于最大圆环外周的多个部分圆环,参见图6所示。由此可见,任何形状的目标靶区的截面均可由圆环和/或部分圆环填满,即任何形状的目标靶区均可由截面为圆环和部分圆环的旋转层层叠而成。
由于肿瘤形状的千变万化,其旋转层在旋转轴方向可能出现断裂,参见图7所示,该目标靶区为哑铃状,两端粗,中间部分细,两粗端的部分相同半径的旋转层在中间部分出现了断裂缺失。虽然二者在空间上出现了断裂,但临床放疗中,通常要求二者具有相同的剂量分布。因此在一些实施例中,将具有相同剂量分布但位于不同分区的旋转层作为同一个旋转层。具体到图7中的哑铃状目标靶区,将两粗端中半径相同的旋转层作为一个旋转层。
S102、根据目标靶区的处方剂量分布确定目标靶区的每个旋转层的目标剂量分布。
在确定了目标靶区的旋转层分布之后,根据目标靶区的处方剂量确定该目标靶区的各个旋转层的剂量分布。其中,剂量分布可以是极坐标下的剂量分布,比如角度-y轴坐标值-剂量,参见图8A所示;也可以是直角坐标系下的剂量分布,比如,x轴表示旋转角度,Y轴表示距离,z轴表示剂量,参见图8B所示。图8A中的二维极坐标截面的中心位于Y轴(旋转轴)上,各个二维极坐标截面之间的距离优选相同,比如1mm。每个二维极坐标截面与剂量分布曲面相交得到此截面下的剂量分布曲线。剂量分布曲线与所有角度剂度线相交于一点,该点可用(θ,y,d)表示。
S103、根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成目标靶区的放疗计划,计划参数包括束流参数和旋转层的叶片对匹配信息。
其中,计划参数包括但不限于束流参数、旋转层的叶片对匹配信息、叶片状态信息等,束流参数包括但不限于束流的能量、剂量率等。
旋转层的叶片对匹配信息包括建立的叶片对与旋转层之间的对应关系,以使放疗计划被执行时,束流输出装置控制叶片对向对应的旋转层输出束流。叶片对状态信息用于限定每个叶片对的各种状态的改变时间和维持时间等。
由于每个叶片对仅能行成一个缝隙,因此只能向一个连续的旋转层,或者向断裂旋转层的一个分区输出束流。可以理解的是,在旋转层出现断裂时,其各个分区可以被设置为同时接收束流,也可以被设置为非同时接收束流,本实施例优选前者。前者的设置方法可选为:先确定断旋转层;然后分别为该旋转层的至少两个分区匹配叶片对,以使该旋转层的叶片对匹配信息被放疗设备执行时,相应的至少两个叶片对在旋转过程中同时向该旋转层的至少两个分区输出束流,参见图9。从图9还可以看出,对于目标靶区来说,某个断裂旋转层的两个分区分别在不同的角度接收对应叶片对输出的束流,但对于束流输出装置来说,两个叶片对在同一旋转角度同时向对应分区的相同旋转层输出射线束。
可以理解的是,在使多个叶片对向同一旋转层的多个分区同时输出束流时,需保证该同一旋转层的多个分区所需的束流参数与该多个叶片对输出的束流参数相同。
由于束流参数与旋转层的角度对应,而且每个叶片对可以独立控制,因此在目标靶区的旋转层数少于叶片对数目时,优选使得多个在同一角度且具有相同束流参数的旋转层在该角度同时接收束流。为了实现该目的,在一些实施例中,如果两个旋转层在同一角度的束流参数相同,则该两个旋转层的该同一角度的叶片对匹配信息被设置为,在该叶片对匹配信息被放疗设备执行时,至少两个相应的叶片对在该同一角度同时向该两个旋转层输出该相同束流参数的束流。可以理解的是,同时向多个旋转层输出束流,可以显著地减少束流输出装置向目标靶区输出束流的总时间,即减少患者接收放射治疗的时间。
本发明实施例提供的放射治疗系统的技术方案,相较于现有技术,将不规则的目标靶区分解成若干形状规则且同心层叠旋转层,且旋转层的截面为环形和/或部分圆环形,以及根据目标靶区的处方剂量分布确定每个旋转层的目标剂量分布,这样即可根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,从而确定目标靶区的放疗计划。由于任何形状的目标靶区的截面均可由同心层叠的环形和/或部分圆环形拼凑成,那么在环形厚度较小时,各个旋转层的剂量分布具有较高的均匀性,从而使得目标靶区的各个评价指标数据都比较好。
实施例二
图10是本发明实施例二提供的放疗计划制作方法的流程图。本发明实施例在上述实施例的基础上,对旋转层的剂量分布确定方法进行了详细说明。
相应地,本实施例的方法包括:
S201、根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,目标靶区由至少两个旋转层同心层叠而成,旋转层的截面为圆环或部分圆环,圆环或部分圆环的圆心位于所述目标靶区的中心轴上。
S2021、根据目标靶区的处方剂量分布确定目标靶区的每个旋转层的处方剂量分布。
整个目标靶区的处方剂量分布确定之后,即可根据目标靶区的处方剂量分布确定目标靶区的每个旋转层的处方剂量分布。
S2022、根据每个旋转层的处方剂量分布和已有剂量分布确定每个旋转层的目标剂量分布,其中,已有剂量分布由当前旋转层所包围的旋转层对应的束流产生。
参见图11A所示,束流必须经过当前旋转层的外部旋转层才能到达该当前旋转层,因此束流在向当前旋转层输入剂量分布时,会对当前旋转层的所有外部旋转层贡献剂量分布。而且相比更外层的旋转层而言,对其相邻的外层旋转层的剂量贡献最大。在图11A中,当前旋转层中的阴影区域表示当前旋转层吸收的剂量,相邻外层旋转层中的阴影区域表示该相邻外层旋转层吸收的剂量,两个阴影区的面积与两个旋转层的平均剂量成正比,与两个旋转层的面积比成反比。换言之,越是靠近中心区域,旋转层的面积比越大,两旋转层间的梯度分布也就越大;越是远离中心区域,两旋转层的面积比越小,两旋转层间的剂量分布梯度也就越小,参见表1所示。
表1层厚为1cm的旋转层之间的剂量分布梯度表
旋转层序号 半径 面积 外圈相邻层剂量衰减比例 外圈相邻层剂量吸收比例
0 0 0
1 10 314 -66.67% 33.33%
2 20 942 -50.00% 50.00%
3 30 1884 -40.00% 60.00%
4 40 3140 -33.33% 66.67%
5 50 4710 -28.57% 71.43%
6 60 6594 -25.00% 75.00%
7 70 8792 -22.22% 77.78%
8 80 11304 -20.00% 80.00%
9 90 14130 -18.18% 81.82%
10 100 17270 -16.67% 83.33%
11 110 20724 -15.38% 84.62%
12 120 24492 -14.29% 85.71%
13 130 28574 -13.33% 86.67%
14 140 32970 -12.50% 87.50%
15 150 37680 -11.76% 88.24%
16 160 42704
该表以层厚为1cm的旋转层为例来说明不同旋转层对应的剂量衰减情况。由该表中不同旋转层的相邻外部旋转层的剂量衰减比例可知,第1个旋转层到第2个“旋转层”的剂量衰减比例高达66.67%,而第15个旋转层到第16个旋转层的剂量衰减比例仅为11.76%。如果计划靶区在第15个旋转层之内,为了让计划靶区得到100%的剂量,其相邻层正常组织吸收的剂量达到了88.4%,这个值显然是临床不允许的。
如图11B所示,当束流的宽度变窄后,两个阴影区的面积比例会明显改变。束流的宽度越窄,当前旋转层的吸收剂量比例越大。因此可以通过束流宽度的改变,使得外部旋转层的剂量分布衰减梯度的增加,以及使正常组织的剂量分布大幅减少。可以理解的是,对于给定的放疗设备,其束流宽度最窄为叶片的厚度。
由束流传输特性可知,任一旋转层的束流对该外层旋转层所包围的内层旋转层的影响较小,为了便于技术方案的阐述,本实施例设定任一旋转层的束流对该外层旋转层所包围的内层旋转层的剂量分布没有影响。
基于此,除了最内层的旋转层,在计算其他任意旋转层的剂量分布时,需要考虑由内部旋转层的束流产生的已有剂量分布。在确定该已有剂量分布时,先确定当前旋转层所包围的各个内部旋转层分别对应的束流,以及所确定的束流在到达对应旋转层时分别对当前旋转层贡献的剂量分布;然后计算所有内部旋转层对当前旋转层贡献的剂量分布的累加结果,该累加结果即为已有剂量分布。
可以理解的是,半径最小的旋转层的已有剂量分布为空,在计算该旋转层的各个位置的计划信息时,无需考虑已有剂量分布的影响。
S203、根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成目标靶区的放疗计划,其中,计划参数包括束流参数和旋转层的叶片对匹配信息。
本发明实施例根据每个旋转层的已有剂量分布和处方剂量分布逐一确定每个旋转层的目标剂量分布,即在计算当前旋转层的剂量分布时,从当前旋转层的处方剂量减去由内部旋转层贡献的已有剂量分布,防止外部旋转层的剂量超限,提高了各个旋转层剂量分布计算的准确性。
实施例三
图12是本发明实施例三提供的放疗计划制作方法的流程图。本发明实施例在上述实施例的基础上,对旋转层的束流参数设计进行了详细说明。
相应地,本实施例的方法包括:
S301、根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,目标靶区由至少两个旋转层同心层叠而成,旋转层的截面为圆环或部分圆环,圆环或部分圆环的圆心位于目标靶区的中心轴上。
S302、根据目标靶区的处方剂量分布确定目标靶区的每个旋转层的目标剂量分布。
S3031、根据每个旋转层的目标剂量分布确定每个旋转层的各个角度对应的束流参数。
在确定旋转层的各个角度对应的束流参数时,优选先确定旋转层的各个角度是否可以采用已有的束流参数,即先确定当前旋转层的各个角度是否可以使用内部旋转层已使用的束流参数,或者当前旋转层的其他角度已使用的束流参数,如果可以,则为可以使用已有束流参数的角度匹配已有束流参数,如果不可以,则为不能使用已有束流参数的角度匹配新的束流参数。
S2032、将每个旋转层的各个角度对应的束流参数划分为至少两个扫描分组,以使束流输出装置在任一扫描分组的扫描过程中,在从一个旋转角度旋转至另一旋转角度时能够完成束流参数的调整。
可以理解的是,由于目标靶区各个方向距离人体表面的深度不一样,而且目标靶区通常不是规则的圆柱体,因此目标靶区的各个旋转层的各个角度对应束流参数会不一样。而且在一些位置上,各个角度的束流参数会存在较大幅度的跳变,这使得束流输出装置在从一个旋转角度旋转至另一旋转角度时,来不及完成束流参数的调整。为了解决该问题,本实施例将旋转层的各个角度对应的束流参数划分为至少两个扫描分组,并使每个扫描分组内的束流参数相同,或是符合预设变化趋势,比如逐渐增大、逐渐减小等,以使束流输出装置在绕目标靶区旋转一周的过程中,在从每个扫描分组的一个旋转角度旋转至相邻的另一旋转角度时,可以完成束流参数的调整。这样束流输出装置可以不停顿地依次完成各个扫描分组的旋转角度的束流的输出。
本实施例在为不能使用已有束流参数的旋转角度匹配新的束流参数时,优选匹配不会增加新的扫描分组的束流参数。
可以理解的是,束流输出装置绕目标靶区旋转一周,至少可以完成一个扫描分组内各角度的束流输出。为了降低输出装置的旋转次数,本实施例在为不能使用已有束流参数的角度匹配新的束流参数时,优选匹配不会增射线加输出装置的旋转圈数的束流参数。换言之,如果必须增加扫描分组,那么所增加的扫描分组不会增加束流输出装置的旋转圈数。即所增加的扫描分组与已有扫描分组中的至少一个扫描分组的束流,可由束流输出装置在旋转一圈的过程中都输出完。
S2033、根据所有旋转层的扫描分组结果确定剂量模拟结果。
所有旋转层的扫描分组完成之后,根据所有旋转层的扫描分组结果对目标靶区进行剂量模拟以生成剂量模拟结果。
需要说明的是,使用现有技术的剂量模拟方法对目标靶区进行剂量模拟即可,本实施例在此不作具体限定。
S2034、在剂量模拟结果的评价指标数据符合预设评价条件时,根据各个旋转层的各个角度对应的束流参数和扫描分组结果生成目标靶区的放疗计划。
其中,评价指标数据包括但不限于靶区适形指数、均一性、危及器官剂量等。靶区适形度指数的预设评价条件是VPTV95%的体积,即95%等剂量线所包围的体积至少应接受处方剂量或更高的剂量。均一性的预设评价条件为:PTV可以接受的最低剂量为95%的处方剂量,可以接受的最大剂量为110%的处方剂量。危及器官剂量的预设评价条件为,是否超出医生处方允许的最大耐受量。
目标靶区的剂量模拟结果得到之后,确定该剂量模拟结果的评价指标数据是否符合预设评价条件,即确定当前的扫描分组结果对应的剂量分布是否符合临床要求,如果符合,则根据各个旋转层的各个角度对应的束流参数、扫描分组结果和前述实施例中的叶片对匹配信息生成目标靶区的放疗计划;如果不符合,则对至少一个旋转层的扫描分组结果进行优化,然后确定重新分组后的扫描分组结果对应的剂量模拟结果,并判断该剂量模拟结果是否符合临床要求,直至所确定的扫描分组结果对应的剂量模拟结果符合临床要求。
对于扫描分组的优化,以射线束参数为剂量率为例进行说明。在一些实施例中,基于放疗计划均一性要求,扫描分组优化方法包括:对旋转层的至少一个角度的剂量率进行调整,以在剂量偏差允许的范围内降低扫描分组中的各个角度的剂量率的差异性,或者减少变化剂量率的使用,从而降低输出装置调整剂量率的幅度和次数,以及提高束流输出装置输出的束流剂量率的准确性。
本发明实施例提供的放疗计划的技术方案,从内容到外层,逐步确定每个旋转层的各个位置对应的束流参数,并将每个旋转层的角度-束流参数划分为至少两组,使得束流输出装置在任一扫描分组的扫描过程中,在从一个旋转角度旋转至另一旋转角度时能够完成束流参数的调整,从而使得束流输出装置在旋转过程中无需任何停歇即可完成任一旋转角度的束流的输出,有助于提高束流的输出效率和稳定性。
实施例四
本发明实施例四还提供了一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行以下方法:
根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,所述目标靶区由至少两个旋转层同心层叠而成,所述旋转层的截面为圆环或部分圆环,所述圆环或部分圆环的圆心位于所述目标靶区的中心轴上;
根据所述目标靶区的处方剂量分布确定所述目标靶区的每个旋转层的目标剂量分布;
根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成所述目标靶区的放疗计划,其中,所述计划参数包括束流参数和旋转层的叶片对匹配信息。
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作,还可以执行本发明任意实施例所提供的放疗方法中的相关操作。
通过以上关于实施方式的描述,所属领域的技术人员可以清楚地了解到,本发明可借助软件及必需的通用硬件来实现,当然也可以通过硬件实现,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如计算机的软盘、只读存储器(Read-Only Memory,简称ROM)、随机存取存储器(RandomAccess Memory,简称RAM)、闪存(FLASH)、硬盘或光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的放疗方法。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。

Claims (10)

1.一种放疗计划系统,其特征在于,该系统包括处理器:所述处理器在工作时实现以下步骤:
根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,所述目标靶区由至少两个旋转层同心层叠而成,所述旋转层的截面为圆环或部分圆环,所述圆环或部分圆环的圆心位于所述目标靶区的中心轴上;
根据所述目标靶区的处方剂量分布确定所述目标靶区的每个旋转层的目标剂量分布;
根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成所述目标靶区的放疗计划,其中,所述计划参数包括束流参数和旋转层的叶片对匹配信息。
2.根据权利要求1所述的系统,其特征在于,将半径相同的旋转层作为同一个旋转层。
3.根据权利要求1所述的系统,其特征在于,在放疗计划确定时,按照半径尺寸从小到大的顺序,依次确定每个旋转层的目标剂量分布和各个位置对应的束流参数。
4.根据权利要求1所述的系统,其特征在于,旋转层的目标剂量分布确定方法,包括:
根据所述目标靶区的处方剂量分布确定所述目标靶区的每个旋转层的处方剂量分布;
根据每个旋转层的处方剂量分布和已有剂量分布确定每个旋转层的目标剂量分布,其中,已有剂量分布由当前旋转层所包围的旋转层对应的束流产生。
5.根据权利要求1所述的系统,其特征在于,所述根据每个旋转层的剂量分布确定每个旋转层的各个位置对应的计划参数,以生成所述目标靶区的放疗计划,所述计划参数包括束流参数和旋转层的叶片对匹配信息,包括:
根据每个旋转层的目标剂量分布确定每个旋转层的各个角度对应的束流参数;
将每个旋转层的各个角度对应的束流参数划分为至少两个扫描分组,以使束流输出装置在任一扫描分组的扫描过程中,在从一个旋转角度旋转至另一旋转角度时能够完成束流参数的调整;
根据所有旋转层的扫描分组结果确定剂量模拟结果;
在所述剂量模拟结果的评价指标数据符合预设评价条件时,根据各个旋转层的各个角度对应的束流参数和扫描分组结果生成所述目标靶区的放疗计划。
6.根据权利要求5所述的系统,其特征在于,确定旋转层的各个角度对应的束流参数的方法,包括:
确定当前旋转层的各个角度是否可采用已有的束流参数;
若是,则为对应角度匹配已有的束流参数;
若否,则为对应角度匹配新的束流参数。
7.根据权利要求6所述的系统,其特征在于,同一扫描分组的束流参数相同,或者束流参数随着角度的变化符合预设变化趋势;
所述新的束流参数为不增加新的扫描分组或不增加束流输出装置的旋转圈数的束流参数。
8.根据权利要求5所述的系统,其特征在于,在所述根据所有旋转层的扫描分组结果确定剂量模拟结果之后,还包括:
如果该剂量模拟结果的评价指标数据不符合预设评价条件,则重新优化至少一个旋转层的扫描分组结果,直至所确定的扫描分组结果对应的剂量模拟结果的评价指标数据符合预设评价条件。
9.根据权利要求1-8任一所述的系统,其特征在于,所述旋转层的叶片对匹配信息的确定方法包括:
如果两个旋转层在同一角度的束流参数相同,则该两个旋转层的该同一角度的叶片对匹配信息被设置为,在该叶片对匹配信息被放疗设备执行时,至少两个相应的叶片对在该同一角度同时向该两个旋转层输出该相同束流参数的束流。
10.一种包含计算机可执行指令的存储介质,其特征在于,所述计算机可执行指令在由计算机处理器执行时用于执行如下方法:
根据目标靶区的径向半径和旋转层的径向厚度确定目标靶区的旋转层分布,其中,所述目标靶区由至少两个旋转层同心层叠而成,所述旋转层的截面为圆环或部分圆环,所述圆环或部分圆环的圆心位于所述目标靶区的中心轴上;
根据所述目标靶区的处方剂量分布确定所述目标靶区的每个旋转层的目标剂量分布;
根据每个旋转层的目标剂量分布确定每个旋转层的各个位置对应的计划参数,以生成所述目标靶区的放疗计划,其中,所述计划参数包括束流参数和旋转层的叶片对匹配信息。
CN202010568528.4A 2020-06-19 2020-06-19 放疗计划系统和存储介质 Active CN111785378B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010568528.4A CN111785378B (zh) 2020-06-19 2020-06-19 放疗计划系统和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010568528.4A CN111785378B (zh) 2020-06-19 2020-06-19 放疗计划系统和存储介质

Publications (2)

Publication Number Publication Date
CN111785378A true CN111785378A (zh) 2020-10-16
CN111785378B CN111785378B (zh) 2022-04-26

Family

ID=72757511

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010568528.4A Active CN111785378B (zh) 2020-06-19 2020-06-19 放疗计划系统和存储介质

Country Status (1)

Country Link
CN (1) CN111785378B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105457172A (zh) * 2015-12-25 2016-04-06 上海联影医疗科技有限公司 放射治疗中照射角度的设置方法及装置
CN105617536A (zh) * 2015-12-24 2016-06-01 上海联影医疗科技有限公司 旋转逆向调强放疗优化方法及装置
CN106659906A (zh) * 2015-11-16 2017-05-10 数码医疗集团 治疗计划制定方法、装置及放疗系统
CN107545137A (zh) * 2017-08-16 2018-01-05 强深智能医疗科技(昆山)有限公司 肿瘤放射治疗计划智能优化方法
US20180243584A1 (en) * 2017-02-24 2018-08-30 Varian Medical Systems International Ag Radiation treatment planning and delivery using collision free regions
US20180304098A1 (en) * 2017-04-21 2018-10-25 Varian Medical Systems, Inc. Dual-axis ring gantry radiotherapy systems
CN108744313A (zh) * 2018-06-25 2018-11-06 西安大医数码科技有限公司 放疗计划规划方法、放疗计划系统及放射治疗系统
CN109513119A (zh) * 2018-11-12 2019-03-26 武汉大学 一种容积旋转调强放疗计划准直器角度优化方法
CN110349665A (zh) * 2019-07-01 2019-10-18 复旦大学附属肿瘤医院 基于深度学习的直肠癌放疗计划自动化设计方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106659906A (zh) * 2015-11-16 2017-05-10 数码医疗集团 治疗计划制定方法、装置及放疗系统
CN105617536A (zh) * 2015-12-24 2016-06-01 上海联影医疗科技有限公司 旋转逆向调强放疗优化方法及装置
CN105457172A (zh) * 2015-12-25 2016-04-06 上海联影医疗科技有限公司 放射治疗中照射角度的设置方法及装置
US20180243584A1 (en) * 2017-02-24 2018-08-30 Varian Medical Systems International Ag Radiation treatment planning and delivery using collision free regions
US20180304098A1 (en) * 2017-04-21 2018-10-25 Varian Medical Systems, Inc. Dual-axis ring gantry radiotherapy systems
CN107545137A (zh) * 2017-08-16 2018-01-05 强深智能医疗科技(昆山)有限公司 肿瘤放射治疗计划智能优化方法
CN108744313A (zh) * 2018-06-25 2018-11-06 西安大医数码科技有限公司 放疗计划规划方法、放疗计划系统及放射治疗系统
CN109513119A (zh) * 2018-11-12 2019-03-26 武汉大学 一种容积旋转调强放疗计划准直器角度优化方法
CN110349665A (zh) * 2019-07-01 2019-10-18 复旦大学附属肿瘤医院 基于深度学习的直肠癌放疗计划自动化设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LI, MING-HSIEN等: "Variations in dosimetric distribution and plan complexity with collimator angles in hypofractionated volumetric arc radiotherapy for treating prostate cancer", 《JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS》 *
王彬冰: "放射治疗中角度最优化问题的研究", 《中国优秀硕士学位论文全文数据库(医药卫生科技辑)》 *
袁青青: "加速器运行误差对容积旋转调强放射治疗影响的研究", 《中国优秀硕士学位论文全文数据库(医药卫生科技辑)》 *

Also Published As

Publication number Publication date
CN111785378B (zh) 2022-04-26

Similar Documents

Publication Publication Date Title
Audet et al. Evaluation of volumetric modulated arc therapy for cranial radiosurgery using multiple noncoplanar arcs
Gevaert et al. Evaluation of a dedicated brain metastases treatment planning optimization for radiosurgery: a new treatment paradigm?
Cedric et al. Clinical implementation of intensity-modulated arc therapy
US7839974B2 (en) ARC-sequencing technique for intensity modulated ARC therapy
US7333591B2 (en) Method for the planning and delivery of radiation therapy
Rauschenbach et al. A dosimetric comparison of three‐dimensional conformal radiotherapy, volumetric‐modulated arc therapy, and dynamic conformal arc therapy in the treatment of non‐small cell lung cancer using stereotactic body radiotherapy
Lee et al. Comparative analysis of SmartArc‐based dual arc volumetric‐modulated arc radiotherapy (VMAT) versus intensity‐modulated radiotherapy (IMRT) for nasopharyngeal carcinoma
Kang et al. A method for optimizing LINAC treatment geometry for volumetric modulated arc therapy of multiple brain metastases
Pokhrel et al. FFF‐VMAT for SBRT of lung lesions: Improves dose coverage at tumor‐lung interface compared to flattened beams
Merrow et al. A dosimetric evaluation of VMAT for the treatment of non‐small cell lung cancer
Liu et al. Small‐spot intensity‐modulated proton therapy and volumetric‐modulated arc therapies for patients with locally advanced non‐small‐cell lung cancer: a dosimetric comparative study
Bedford et al. Constrained segment shapes in direct‐aperture optimization for step‐and‐shoot IMRT
Diot et al. Biological‐based optimization and volumetric modulated arc therapy delivery for stereotactic body radiation therapy
Skórska et al. Optimization of treatment planning parameters used in tomotherapy for prostate cancer patients
Oliver et al. Analysis of RapidArc optimization strategies using objective function values and dose‐volume histograms
MacDonald et al. Dynamic collimator trajectory algorithm for multiple metastases dynamic conformal arc treatment planning
Knapp et al. The role of volumetric modulated arc therapy (VMAT) in gynaecological radiation therapy: a dosimetric comparison of intensity modulated radiation therapy versus VMAT
US20140330064A1 (en) Station parameter optimized radiation therapy (SPORT): a novel scheme for treatment planning and delivery in radiation therapy
Smith et al. Optimization of spot and trimmer position during dynamically collimated proton therapy
Otto et al. Enhancement of IMRT delivery through MLC rotation
Milette et al. Maximizing the potential of direct aperture optimization through collimator rotation
CN111714790B (zh) 放疗计划系统和存储介质
Lim Optimization in radiation treatment planning
CN111785378B (zh) 放疗计划系统和存储介质
Qi et al. Relationship of segment area and monitor unit efficiency in aperture‐based IMRT optimization

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