CN104066479A - 使用异质性补偿迭加进行用于放射疗法的剂量计算 - Google Patents
使用异质性补偿迭加进行用于放射疗法的剂量计算 Download PDFInfo
- Publication number
- CN104066479A CN104066479A CN201280058918.0A CN201280058918A CN104066479A CN 104066479 A CN104066479 A CN 104066479A CN 201280058918 A CN201280058918 A CN 201280058918A CN 104066479 A CN104066479 A CN 104066479A
- Authority
- CN
- China
- Prior art keywords
- radiotherapy
- calculating
- health
- per unit
- positions
- 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
Links
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
- 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/103—Treatment planning systems
-
- 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
-
- 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/1064—Monitoring, verifying, controlling systems and methods for adjusting radiation treatment in response to monitoring
- A61N5/1065—Beam adjustment
- A61N5/1067—Beam adjustment in real time, i.e. during treatment
-
- 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/1075—Monitoring, verifying, controlling systems and methods for testing, calibrating, or quality assurance of the radiation treatment apparatus
-
- 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
- A61N2005/1032—Genetic optimization methods
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
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
一种用于放射疗法的***,包括放射计划***。所述放射计划***包括数据处理器,所述数据处理器适于接收与身体的预期放射治疗区域有关的信息,接收计算出的身体内多个位置每单位质量释放的初始能量,基于计算出的每单位物质释放的初始能量并包括由于身体内的其他位置的分散带来的放射剂量贡献而计算放射治疗区域内多个位置的放射剂量,以及基于在放射治疗区域内多个位置计算的放射剂量而确定放射疗法参数用于对预期放射治疗区域提供放射治疗。包括由于身体内的其他位置的分散带来的放射剂量贡献考虑到了身体内的密度不连续。
Description
相关申请的交叉引用
本申请要求于2011年9月29日提交的美国临时申请第61/540,773号的优先权,其全部公开内容通过引用的方式并入本文。
背景技术
1.发明领域
本发明涉及放射治疗***,更具体而言,涉及放射物质非匀质分布的放射治疗***。
2.相关技术的讨论
放射治疗是医学使用放射来治疗恶性的细胞,例如癌细胞。这种放射可以具有电磁的形式,例如高能光子,或者具有微粒的形式,例如电子、质子、中子或者阿尔法粒子。
迄今为止,在现今实践中使用得最普遍的放射形式是高能光子。光子在人体组织中的吸收是由放射的能量及相关组织的原子结构所确定的。在放射肿瘤学中采用的能量的基本单位是电子伏特(eV);103eV=1keV,106eV=1MeV。在治疗性能量的光子和组织之间的三个主要的相互作用是:光电效应、康普顿效应(Compton effect)和电子偶的产生(pair production)。
在光电效应中,传入的光子将能量传递给被紧密束缚着的电子。光子几乎将它所有的能量传递给电子,然后不复存在。电子利用来自光子的大部分能量而脱离,并且开始使周围的分子离子化。这种相互作用取决于传入的光子的能量以及组织的原子数量;能量越低且原子数量越高,则光电效应越有可能发生。使光电效应在组织中占优势的能量范围约为10keV到25keV。
康普顿效应是用于癌症治疗的最重要的光子与组织的相互作用。在这种情况下,光子与“自由电子”(即没有被紧密束缚到原子的电子)碰撞。与光电效应不同,在康普顿相互作用中,光子与电子都被散射。虽然能量较低,但光子之后仍可以继续承受另外的相互作用。电子利用光子给它的能量开始离子化。康普顿相互作用的概率与传入的光子的能量成反比,并且独立于材料的原子数量。康普顿效应在25keV到25MeV的范围内占优势,并且因为大多数放射治疗在大约6MeV到20MeV的能量下执行,所以康普顿效应是临床上最经常发生的相互作用。
在电子偶的产生中,光子与原子核相互作用。光子把能量给原子核,并且在此过程中创造粒子的正电子-电子偶。正的电子(正电子)离子化直到其在正电子-电子湮没中与自由电子结合。这种正电子-电子湮没产生在相反方向上运动的两个光子。电子偶产生的概率与传入的光子的能量的对数成正比,并且取决于材料的原子数量。使电子偶占优势的能量范围为大于等于25MeV。这种相互作用在某种程度上在使用高能光子射束的常规放射治疗中发生。
随着高能线性加速器的出现,在治疗深度达到大约5cm的浅表肿瘤方面电子变为可行的选择。电子深度剂量特性的独特之处在于其产生高皮肤剂量,但是仅在几厘米后就表现出衰减。
电子在人体组织中的吸收极大的受到存在气孔(air cavities)和骨骼的影响。电子射束最通常的临床使用包括对皮肤病变(比如基底细胞癌,以及先前接收过光子放射的面积(比如乳癌病人在***肿瘤切除术或***切除术后的伤疤)的增大,以及在头部和颈部选择的结节区域)的治疗。
快速、精确的剂量计算算法对于放射疗法计划而言十分重要,因为该放射疗法计划是确保对特定病人给予所希望的剂量的唯一可用的方法。剂量计算包括两部分:源模型和传输模型。源模型提供入射流量。传输模型计算由入射流量产生的剂量,并且目前是性能瓶颈。三个主要的传输算法以增加的精确性/减少的性能的次序依次是笔状射束、迭加/卷积和蒙特卡洛(Monte Carlo)。迭加/卷积是计算外部射束放射疗法的放射剂量的当前临床标准方法。
近年来,通过使用强度调节提高了治疗质量。这种技术使用多叶准直器限定根据单一射束方向的多孔径,以提供改变穿过射束的放射强度的能力。这种技术允许使放射治疗符合目标的形状,并且在射束参数的数量急剧增长时避免临界的结构。为了确定多叶准直器最佳的一组设置,治疗计划***必须通过剂量计算的多次迭代使射束参数的数量急剧增长的目标函数最优化。在实践中,为了对病人达到可能的最好结果,治疗计划者多次重复该最优化。因此,当对一组五个射束的单一最优化可能需要五分钟时,整个过程可能需要几小时来产生临床可接受的计划。这就限制了在临床工作流程中强度调节计划的数量和质量。
这种临床工作流程的局限性延伸到更复杂的技术,例如容积调节的弧形疗法(Otto,K.,Med.Phys.35,310-317,2008)、强度调节的弧形疗法(Yu,C.X.,Phys.Med.Biol.40,1435-1449,1995)和自适应放射疗法(Yan,D.,Vicini,F.,Wong,J.,Martinez,A,Phys.Med.Biol.42,123-132,1997)。此外,这种临床工作流程的局限性禁止实时的放射疗法;每天扫描、重新计划和治疗每个病人的能力。从Ahnesjo等人(Ahnesjo,A.,Aspradakis,M,Phys.Med.Biol.44,R99-Rl551999)可得到对放射疗法中的剂量计算的全面的讨论。
因此,剂量计算的计算性能是放射疗法的治疗计划的质量的限制因素。传统地,通过更快的硬件已经实现了治疗质量方面的改进。但是,摩尔定律(Moore′slaw)已改变。计算机不是每18个月速度加倍,而是处理内核(kernel)的数量加倍。并且随着处理器变为多核的,图形处理单元(GPU)的多核架构获得了运行一般用途的算法的灵活性。为了根据计算机硬件的近来趋势实现许诺的性能增益,在放射剂量计算中使用的传统的串行算法应该被并行算法取代。最近,Nucletron公司宣布了在其治疗计划***中的GPU加速,不过尚未获得其公布的细节。然而,将现有的串行算法直接分隔来产生用于多处理核的多线程行不通。这是因为线程被外包(farm out)来基于相同输入数据计算相同放射剂量,从而容易出现读/写冲突。当覆写(write on write,WOW)冲突出现时,例如可能发生仅存储了最后的写入内容,从而导致剂量计算不精确。此外,执行实时放射治疗计划的传统方法,对被辐照物质内的不连续难以提供好结果。因此,仍需改进放射剂量计划的方法和***,以将被辐照体的密度不均匀分布考虑在内。
发明内容
根据本发明的一些实施方式的用于放射疗法的***包括放射计划***。所述放射计划***包括数据处理器,所述数据处理器适于接收与身体的预期放射治疗区域有关的信息,接收计算出的身体内多个位置每单位质量(per unit mass)释放的初始能量,基于计算出的每单位物质释放的初始能量及包括由于身体内的其他位置的分散带来的放射剂量贡献而计算放射治疗区域内多个位置的放射剂量,以及基于在放射治疗区域内多个位置计算的放射剂量而确定放射疗法参数用于对预期放射治疗区域提供放射治疗。包括由于身体内的其他位置的分散带来的放射剂量贡献考虑到身体内的密度不连续。
根据本发明的一些实施方式的用于确定用于放射治疗的放射疗法参数的方法,包括接收与身体的预期放射治疗区域有关的信息,接收计算出的身体内多个位置每单位质量(per unit mass)释放的初始能量,基于计算出的单位物质释放的初始能量及包括由于身体内的其他位置的分散带来的放射剂量贡献而计算放射治疗区域内多个位置的放射剂量,以及基于在放射治疗区域内多个位置计算的放射剂量而确定放射疗法参数用于对预期放射治疗区域提供放射治疗。包括由于身体内的其他位置的分散带来的放射剂量贡献考虑到身体内的密度变化。
根据本发明的一些实施方式的计算机可读介质包括非临时性可执行代码。当由计算机执行时,所述非临时性可执行代码促使计算机接收与身体的预期放射治疗区域有关的信息,接收计算出的身体内多个位置每单位质量(per unitmass)释放的初始能量,基于计算出的单位物质释放的初始能量及包括由于身体内的其他位置的分散带来的放射剂量贡献而计算放射治疗区域内多个位置的放射剂量,以及基于在放射治疗区域内多个位置计算的放射剂量而确定放射疗法参数用于对预期放射治疗区域提供放射治疗。包括由于身体内的其他位置的分散带来的放射剂量贡献考虑到身体内的密度不连续。
附图说明
从对说明书、附图和示例的讨论中,更多目标和优点将会变得显而易见。
图1示出本发明一种实施方式的示意图。
图2示出传统的正向射线跟踪。
图3示出根据本发明一种实施方式的后投影射线跟踪。
图4示出使用市售的Pinnacle***(Philips Radiation Oncology SystemsMadison WI)和本发明一种实施方式以分析法计算出的沿中心轴的每单位质量释放的总能量(Total Energy Released per unit Mass,TERMA),其中每个TERMA以10cm的深度归一化。
图5A示出使用正向跟踪计算出的有离散化伪影的TERMA的半影片段。
图5B示出使用本发明一种实施方式计算出的没有离散化伪影的同一TERMA的半影片段。
图6A和6B分别示出根据本发明一种实施方式使用固定步长和精确放射线距离计算出的TERMA的片段。
图7A和7B示出在迭加操作中倾斜的内核和不倾斜的内核。
图8A和8B分别示出根据传统的统一采样的存储器访问模式和所计算出的小野(5mm)剂量沉积片段的图。
图9A和9B分别示出根据本发明一种实施方式的存储器访问模式和所计算出的小野(5mm)剂量沉积片段的图。
图10示出由市售的Pinnacle***(Philips Radiation Oncology SystemsMadison WI)和本发明一种实施方式计算出的沿中心轴的吸收剂量图,其中每个吸收剂量图以10cm的深度归一化。
图11示出由市售的Pinnacle***(Philips Radiation Oncology SystemsMadison WI)和本发明一种实施方式计算出的在10cm的深度的吸收剂量分布图,其中每个吸收剂量分布图在中点处归一化。
图12示出根据本发明另一实施方式的一种方法。
图13-42提供根据本发明一种实施方式的包括异质补偿迭加(HeterogeneityCompensated Superposition,HCS)结果的实例。
具体实施方式
下面将详细讨论本发明的一些实施方式。在描述实施方式时,为了清楚起见而采用了特定的术语。然而,本发明并非意在受限于由此而选择的特定术语。本领域的技术人员可认识到在不偏离本发明的广义构思的情况下可以采用其它等同的部件和开发出其它方法。通过引用将所举出的所有参考文献合并于此,如同每个参考文献是单独合并的一样。
图1是根据本发明一种实施方式的放射疗法***100的示意图。放射疗法***100包括:放射计划***101,放射计划***101进一步包括数据处理器102。数据处理器102可以是单一处理器、多元处理器、在网络上分布的分布式处理器和/或并行处理器,例如但不限于,一个或多个图形处理单元(GPUs)。数据处理器102适于接收与具有预期放射治疗区域的身体105有关的输入信息。数据处理器102也适于产生用于为身体105的预期放射治疗区域提供放射治疗的输出信息。数据处理器102可适于在确定用于提供放射治疗的输出信息时基于所接收到的输入信息执行多个反向射线跟踪计算。每个反向射线跟踪包括:计算与身体的预期放射治疗区域的第一子区域对应的第一物理性质,该第一子区域被从源位置穿过预期放射治疗区域行进的射线穿过;以及在上述首先提及的计算之后,计算与预期放射治疗区域的第二子区域对应的第二物理性质,该第二子区域与上述射线相交于比第一子区域更接近源位置的位置处。放射计划***101还可包括:存储器107、显示器108和输入/输出装置109。存储器107可以例如是硬盘驱动器、CD-ROM驱动器、DVD驱动器、闪存驱动器等。显示器108可以例如是液晶显示器(LCD)、阴极射线管(CRT)监视器、等离子显示器等。输入/输出装置109可包括例如鼠标、键盘、用于在网络或者数据总线上传输数据的接口等。
放射疗法***100还可包括:放射治疗***103,其与放射计划***101进行通信。放射治疗***103进一步包括放射源106。放射源106是直接发出放射束到待治疗的身体105上的源。放射源的实例可以包括:X射线源、伽马射线源、电子束源等。放射源106还可包括对射束进行校准的多叶准直器(MLC)。通过调节MLC的叶片的位置,剂量测定员可以使放射野与身体105的治疗区域的形状相匹配。在一些实施方式中可以包括其它的射束形状和/或轮廓。放射源106可以具有对应的源模型。放射***103可以由放射治疗计划***101控制,例如,以传递强度加以调节的放射能量,从而使得放射治疗符合身体105的预定放射治疗区域的形状。
放射疗法***100还可包括诊断***,该诊断***与放射计划***101进行通信,产生身体105的经验数据。该经验数据可被用作放射计划***101和数据处理器102的输入信息,并且可被用于反向射线跟踪计算。诊断***104包括传感器,以获得身体105的经验数据。诊断***的实例可以是计算机断层扫描(CT)扫描仪、磁共振成像(MRI)扫描仪、正电子发射型断层显像(PET)扫描仪等。身体105例如可以是人类或者动物。
根据本发明一种实施方式,用于放射疗法的***100包括放射计划***,所述放射计划***顺次包括数据处理器102,所述数据处理器102适于接收与身体的预期放射治疗区域有关的信息,接收计算出的身体内多个位置每单位质量(per unit mass)释放的初始能量,基于计算出的单位物质释放的初始能量及包括由于身体内的其他位置的分散带来的放射剂量贡献而计算放射治疗区域内多个位置的放射剂量,以及基于在放射治疗区域内多个位置计算的放射剂量而确定放射疗法参数用于对预期放射治疗区域提供放射治疗。包括由于身体内的其他位置的分散带来的放射剂量贡献考虑到身体内的密度不连续。在计算中被考虑在内的身体内的密度不连续可以包括身体内不同类型物质之间的界面。身体内不同类型物质之间的界面可以包括例如至少两个骨骼、组织,一个内部器官、空气、水和一个人工植入物间的界面。然而,不同密度的物质间的任何界面均意在被包括在本发明的一般性概念中。
在一些实施方式中,在放射治疗区域内多个位置中每一个的放射剂量的计算可以包括身体内多个位置上每单位质量释放的初始能量的综合,其由内核函数(kernel function)加权得到。这可以是例如卷积综合。在一些实施方式中,所述内核函数可以包括经验性确定的有效密度函数,其将身体内的密度变化考虑在内。
在一些实施方式中,所述综合可用于第一放射光束能量,而且在一些实施方式中,在放射治疗区域内多个位置中每一个的放射剂量计算可以进一步包括身体内多个位置上每单位质量释放的初始能量的第二综合,其由内核函数加权得到,用于第二放射光束能量,并且在放射治疗区域内多个位置中每一个的放射剂量计算进一步包括将第一和第二综合加和。在一些实施方式中,这可以延伸至两个以上的能量。例如,在本发明的一些实施方式中,可以执行对四个或更多个不同能量的四个或更多个综合。
在一些实施方式中,用于放射疗法的***可以进一步包括,在接收计算出的由身体内多个位置每单位质量释放的初始能量之前,计算每单位质量释放的初始能量。该计算可以是但不局限于如在下面更详细描述的TERMA计算。然而,在其它实施方式中,其可以是或可以包括每单位质量释放的动能(KERMA)和/或每单位质量释放的散射能(SCERMA)计算。在一些实施方案中,每单位质量释放的初始能量计算可以包括背投影射束追踪计算。
在一些实施方式中,确定放射疗法参数对于向预期放射治疗区域提供放射治疗可足够快以在放射治疗程序期间地实时执行。
示出的迭加/卷积算法用来产生对剂量分配的精确计算(Mackie,T.R.,Scrimger,J.W.,Battista,J.J.,Med.Phys.12,188-196,1985;Mackie,T.R.,Ahnesjo,A.,Dickof,P.,Snider,A,Use of Comp.In Rad.Ther.,107-1101987;Mackie,T.R.,Reckwerdt,P.J.,McNutt,T.R.,Gehring,M.,Sanders,C.,Proceedings of the1996AAPM Summer School,1996)。其包括两个阶段。第一,通过病人的密度表示传输入射流量来计算每个位置处的每单位质量释放的总能量(TERMA)。点r′处具有特定能量E的TERMA(即TE(r′))被定义为将能量E的流量ψE(r′)用相对于水的密度ρmd(r′)以及在点r′处的线性衰减μE(r′)加权,如下述方程1所示。
线性衰减系数μE(r′)也取决于原子材料。尽管康普顿散射在与放射疗法相关的百万伏能量范围内占优势,并且康普顿散射取决于电子密度而不是材料,但是,在临床上,对于正常人体组织,在标准CT数量和密度之间存在分段线性关系。通常而言,在此能量范围内非康普顿的相互作用被视为可以忽略。根据下述方程2,能量E在点r′处的流量ψE(r′)是由源焦点s和在方向r′上的能量E的入射流量ψE,0(r′)所决定的。
然后,迭加/卷积算法通过剂量沉积内核传播上述能量,以确定每个位置处的最终剂量。为了允许剂量沉积内核根据组织不均匀性来实际地进行调整,基于将迭加与卷积区别开来的下述方程3使用点之间的放射线(又名加权的电子密度ρed)距离dρ(r,r′)。
根据下述方程4,点r处的剂量D(r)是由在TERMA体积上的积分并用取决于能量的剂量沉积内核KE加权而计算出来的。标准压缩(collapsed)椎体内核是由在点和内核轴之间的放射线距离和相对角ω来标示的,并且缺少几何距离平方效应。
通常而言,不计算如方程4中的每个体素(voxel)的单能贡献。而是如基于方程5和6的方程7中一样,选择射线角度ω和方向v的离散组合,并且沿着使用的单个多能内核对其进行积分。由于距离平方效应,这由近似的指数内核衰减和任意单个远距离的体素急剧减少的贡献得到证明。射线表现出的立体角量上的增加抵消了距离平方效应。射线方向被选择来平衡几何形状和内核的能量因子。
传统地,根据方程1和2,TERMA是由投射一组将入射流量沉积至体积的射线来计算的。为了数字的精确性,大约四条射线应该穿过每个TERMA体素,应该在病人的表面处开始衰减,并且应该通过贡献给该体素的射线的总长度来对每个体素的流量进行归一化。这种归一化去除了由于发散源造成的正常反向平方流量的下降。因此必须再次应用该归一化,且通常将其应用于TERMA网格(grid)。临床上可接受的速度增强以不使剂量沉积内核倾斜来与每个体素处的射线轴对准。在这种情况下,将发散校正应用于剂量网格会更精确(Papanikolaou,N.,Mackie,T.R.,Meger-Wells,C.,Gehring,M.,Reckwerdt,P.,Med.Phys.20,1327-1336,1993)。
TERMA强烈取决于射束光谱。该光谱关于射束轴旋转对称,并且随着材料的深度而变得稳定。稳定的光谱指的是:由于较低的能量成分被优先散射出射束外并被吸收,因此光谱中较高的能量成分占优势。传统地,使用查找表来将衰减模型化,当该查找表与线性衰减查找表结合时,该查找表具有深度轴、密度轴和偏轴角。该表还需要使用固定步骤射线投射算法,该算法避免了在每一步骤都要进行传统的、成本昂贵的对指数的评估。尽管出于性能原因在临床上接受了该查找表,但该查找表假定介质是均匀的。由于不均匀的组织优先使不同的光谱衰减,因此这是不正确的。此外,固定步长和离散射线导致数字的伪象。
剂量沉积内核还依赖于每个体素处的能量谱。然而,单个多能内核已被展示对临床使用而言足够精确且该内核是当前医疗的标准。多能内核是通过对单能内核谱进行组合来创建的,单能内核是使用蒙特卡洛模拟产生的,该模拟强制单能光子在水球体中心处互相作用并且对沉积在该球体每处的剂量进行计数(Ahnesjo,A.,Andreo,P.,Brahme,A.,Acta.Oncol.,26,49-56,1987;Mackie,T.R.,Bielajew,A.F.,Rogers,D.W.O.,Battista,J.J.,Phys.Med.Biol.33,1-20,1988)。
图形处理算法明显地适于并行或者图形处理单元(GPU)的架构。存在适应于传统图形***的多个可视化算法,例如射线投射、射线跟踪和体积的可视化算法。然而,当可视化算法关注线路性质(比如其整数或最大化)时,用于放射疗法的剂量计算根本上与线路和体积的相互作用密切相关。这种区别致使多个为可视化开发的并行射线跟踪算法不适用于放射疗法。尽管如在(Kriiger,J.,Westermann,R.,2003)中所做的轻易改变先前GPU的执行是不切实际的,并且将不会得到精确的结果,但传统的TERMA算法的一部分可以与体积射线投射粗略地相似。因为剂量沉积主要处理电子的相互作用,因此其与可视化算法根本不同。
根据本发明的一种实施方式,将NVIDA的统一计算设备架构(CUDA)软件开发环境与Digital Mars的D编程语言结合起来用于执行本发明,可使迭加/卷积算法适应于GPU。以下描述该适应的细节。
如同所讨论的,迭加/卷积是两阶段的算法:首先计算TERMA,然后迭加剂量沉积内核。
在TERMA计算中,可以执行标准的正向TERMA算法。该方法需要将分散的积累进行重叠,这在并行运行时将导致读写冲突。当线程同时尝试读、处理然后写入相同的存储器时,发生读写冲突,导致仅有最后的更新被记录。这极大地减少了对每个体素有作用的射线的有效数量,使其低于数字精确性所需的极限。然而,可以使用源的扩散性来创造射线组,保证一直有一个体素是分离的。当单个的组大到足够在GPU上有效运行时,这种串行化引起过多的GPU呼叫开销。传统的3D谱衰减查找表超过结构缓存值,从而降低了性能。通过单个地使分开的谱箱衰减可以提升性能,使用共享存储器来减少寄存器的使用。这减少了对小的2D结构的线性衰减系数的查找,这些线性衰减系数由能量和密度使其参数化。这可以由在GPU上的硬件加速指数来实现。当在不均匀的介质中谱正确地衰减时,也可以提高精确度。然而,数字精确性所需的射线的数量与分辨率的关系并不清楚,离散化的效果会很明显,如同将要讨论的那样,并且计算会包括总的剂量计算时间的很大一部分。
根据本发明的一些实施方式,通过反向TERMA算法可以避免这些问题。通过对方程(1)、(2)和(5)进行整理,将入射流量ψE,0(r′)分为谱权wE和净流量因子ψ0(r′),并且定义衰减因子,获得下述方程8。
图2示出了传统的正向射线跟踪。图3示出了根据本发明一种实施方式的背投影射线跟踪。通过从每个TERMA体素背向源投射射线,沿着路线聚集净衰减,当到达病人的边界时,可以获得早期的射线终结,从而提高性能和精确度。在计算上述定义的衰减因子时,预先编译的查找表可以用于数量μE/ρmd。例如,基于放射治疗区域的以Hounsfield为单位的CT图像,可以从查找表查找这个数量。谱权wE例如可以包含在用于放射治疗的放射源的模型中,并且每个谱箱可以针对给定的射线独立地衰减。
然后可以根据下述方程9通过将针对每个体素的净衰减乘以从源朝向体素的净入射流量来计算TERMA。
T(r′)=Ψ0(r′)A(r′). (9)
图4示出使用市售的Pinnacle3***(Philips Radiation Oncology SystemsMadison WI)和本发明一种实施方式以分析法计算出的沿中心轴的TERMA能量,其中每个TERMA能量以10cm的深度归一化。
尽管方程9是O(n4)的算法,但由于与用于标准正向方法的O(n~3)不同,每个线程仅向其自身的体素写入,因此避免了读写冲突。该算法允许对结构存储器联合进行读访问,从而大幅度地提高了存储性能。当仅改变流量野时,只要在相互使用或者强度调节最优化期间,就可以在多种临床分辨率下使用标准正向和背投影射线跟踪算法,并且可以计算该计算方法的高额衰减量。根据本发明的一些实施方式,这种预先计算可以提供多个量级的性能改进。
大小 | 正向 | 背投影 | 近似的衰减| | 强度调节 |
323 | 23ms | 2ms | 1ms | 1ms |
643 | 44ms | 18ms | 2ms | 1ms |
1283 | 150ms | 220ms | 26ms | 3ms |
2563 | 1869ms | 3317ms | 454ms | 13ms |
表1.标准正向、背投影和近似的TERMA计算方法的性能。也包括在强度调节期间有用的仅更新入射流量野的性能。
表1中比较了多种分辨率下基于衰减量的TERMA算法相对于传统射线投射算法的性能。两种方法均物理结合了对多谱衰减的校正并且针对GPU进行了优化。在正向射线跟踪计算中使用的射线投射可能有读/写冲突。射线投射的可量测性表现出最佳位置和所需的参数调整,以维持穿过每个体素的射线的必要数量。这种调整费力且很不直观。通过背投影射线跟踪计算的基于衰减量的TERMA计算表现出经验性能O(n3.76),该性能在其理论的O(n4)上有微小的改进。根据本发明一些实施方式的背投影射线跟踪计算不会经受读/写冲突。基于快速的、近似的放射线深度的衰减算法使得对标准均匀的材料近似提供约为8倍的加速。在强度调节下介绍了如通常在治疗计划中的强度调节中使用的那样仅更新入射流量的性能。
背投影射线跟踪计算的一个优点是消除由于射线离散化导致的伪像。这是因为,与将每个体素存储器访问的射线减少到一个从而消除相关的离散化伪像的光栅化方法(Arnanatides,J.,Woo.,A.,Eurographics'87,ConferenceProceedings,1987)相似,在背投影射线跟踪方法中可以使用确切的放射线路径。图5A和5B分别示出了根据本发明一种实施方式使用固定步长和精确的放射线距离所计算出的TERMA片段。
此外,物理校正的多谱衰减(Liu,H.H.,Mackie,T.R.,McCullough,E.C.,Med.Phys.24,1729-1741,1997)可被应用于背投影射线跟踪计算。对于相互射束角度变化的情况,可以实现如通常用于当前临床***那样的无需物理校正多谱衰减的更简单的变型,从而可以得到显著的性能改进。在TERMA计算中,物理校正的多谱衰减(Liu,H.H.,Mackie,T.R.,McCullough,E.C.,Med.Phys.24,1729-1741,1997)可被用于背投影射线计算。当前的临床标准是固定步骤的、快速的、近似的正向方法,该方法没有物理校正的多谱衰减。正如根据本发明的一些实施方式所指出的,缓存衰减量的使用可被用来加速TERMA的计算。图6A示出了使用没有多谱衰减的标准正向射线跟踪计算出的有离散化伪像的TERMA的半影片段。图6B示出了有多谱衰减的根据本发明一种实施方式计算出的没有离散化伪像的同一TERMA的半影片段。
在计算出TERMA之后,可以使用剂量沉积内核的迭加来计算被吸收的放射能量的量。该迭加具有两个标准的方程式。正向方程将剂量从TERMA体素传播到周围的剂量体素。该正向方程需要计算给每个病人的体素的剂量并且会经受读写冲突,因为每个剂量体素都有多个TERMA体素做贡献。反向内核方程聚集从周围的TERMA体素给剂量体素的贡献。由于仅计算了感兴趣的体积的剂量,因此这种计算是高效的。因为剂量沉积内核是可逆的,所以这是可行的。
图7A和7B示出了在迭加操作期间倾斜的内核和非倾斜的内核。严格来说,在标准迭加中使用倾斜的内核破坏了内核可逆的假定。然而,在已知源的距离远和内核快速下降的情况下,可逆性仍是一个合理的假定从而在临床中得到使用。虽然如此,内核的快速下降也产生了在通常的临床分辨率下数字采样的问题。有两种标准的可选方案。根据下述方程10,累积内核(cumulative kernel,CK)(Ahnesjo,A.,Med.Phys.16,577-92,1989)表示从射线段到点的剂量沉积。
累积-累积内核(CCK)(Lu,W.,Olivera,G.H.,Chen,M.,Reckwerdt,P.J.,Mackie,T.R.,Phys.Med.Biol.50,655-680,2005)表示从射线段到射线段的剂量沉积。
对于特定的放射线深度dr和角度ω,这两者都源于标准的点到点内核K(dr,ω)的积分:
尽管更精确(尤其在粗糙的分辨率下),但CCK方程传统上比CK方程慢50%。然而,缓存存储器访问和提供专用的线性内插硬件的GPU的结构单元使得能够使用CCK方程且其性能降低可忽略。因此,使用累积-累积内核(Lu,W.,Olivera,G.H.,Chen,M.,Reckwerdt,P.J.,Mackie,T.R.,Phys.Med.Biol.50,655-680,2005)能够增强迭加计算。
通过移动所有索引越过一个体素,串行CPU的执行允许重新使用射线投射索引计算。但是,这阻止了内核倾斜,从而导致大的偏轴角误差,如将在表2中讨论的那样。相反,在GPU上执行的本发明的一些实施方式允许待计算的内核倾斜和非倾斜。内核倾斜传统上导致300%的性能损失(Liu,H.H.,Mackie,T.R.,McCullough,E.C.,Med.Phys.24,1729-1741,1997),但是在GPU上执行的本发明的一些实施方式的性能损失仅为19%。
由于迭加操作在粗糙的分辨率下维持精确性并且内核在边缘处表现出快速的下降,根据本发明的一些实施方式,可以采用将每个射线近似为真实的立体角的多分辨率迭加算法。与射线不同,立体角的宽度随着几何距离而增加。在离散量中,射线的宽度与体素的宽度是成正比的。因此,通过随着几何距离增加体素的宽度,射线可以近似为立体角。这种近似也以对数方式增加了步长,这通过将计算复杂度从O(ωDT1/3)减小到O(ωDlog(T1/3))而提高了性能;其中ω是角的数量,D是剂量体素数量,而T是TERMA体素的数量。
与标准方法相比,根据本发明一种实施方式的多分辨率迭加表现出令人感兴趣的精确性折衷。表2将根据本发明一些实施方式的多分辨率迭加算法的精确性在多个野尺寸和内核射线采样的情况下与标准方法相对比。由于很少有TERMA被射线几何地错过,因此根据本发明一些实施方式的多分辨率迭加算法通常在半影和用于小野尺寸的低剂量区域中能够执行得更好。然而,由于大步长可以引起射束界限被模糊,因此高剂量区域的精确性可以轻微地折衷。使用相同的步长、但不使用多分辨率数据结构的作为本发明的实施方式的一种变型的多分辨率迭加算法的一种变型可表现出缓存性能下降以及可平均增加60%的平均误差。
表2.平均沉积剂量相对于最大沉积剂量Dmax的误差,其野尺寸从1cm到23cm。通过具有被限定为具有比0.3Dmax更大的剂量梯度的半影区域和比0.2Dmax更低的低剂量区域的区域来减小误差。使用4608射线的倾斜内核采样来计算野是方形的、大小被限定在深度为10cm的参考剂量沉积。笔状射束的精确性是通过将迭加内核在半径上截短3cm来近似的。2%到5%的绝对放射量测定误差在临床上是可接受的(Ahnesjo,A.,Aspradakis,M.,Phys.Med.Biol.44,R99-R155,1999)。
通过使用根据本发明一些实施方式的多分辨率网格可以减少当由于射线采样稀疏而导致漏掉整个射束时产生的小野伪像。图8A和8B分别示出根据传统的统一采样的存储器访问模式和所计算出的小野(5mm)剂量沉积片段的图。图9A和9B分别示出使用多分辨率网格的根据本发明一种实施方式的存储器访问模式和所计算出的小野(5mm)剂量沉积片段的图。
然而,更大的步长可能减小剂量沉积内核的精确性。当邻近的体素横穿不同的粗糙分辨率的体素时产生伪像。应注意本实施具有固有的同向性,这将减少不均匀角度采样的益处。
使用测定体积的最大强度投影(MIP)映射(Williams,L.,SIGGRAPHComput.Graph.17,3,1-11,1983),可以实施多分辨率算法,因为其计算高效且具有良好的缓存性能。分辨率的改变可被限于之前每个步骤的最大值并且可仅允许其在分辨率较粗糙的体素边界处发生。这样可以防止TERMA体素多次给同一剂量体素做贡献。
可以使用多个策略来优化统一计算设备架构(CUDA)的性能。CUDA的执行模型是线程的3D块中的2D网格,该线程执行一项函数(称为内核)。每个块代表一个并行工作单元,并且因此在尺寸上有限制。使用NVIDIA的CUDA占用计算器优化了块线程计算。为了体积处理,可以使用在x和y方向上的线程到体素的1∶1的映射。通过以增加的z索引使体素函数循环可以处理z方向。步幅可以是z块大小,其维持了线程空间的凝聚力。CUDA现在也支持3D块中的3D网格。通过立方体状块尺寸增强的提高了的线程空间的凝聚力可以减小缓存遗漏从而提升性能。通常在结构中可以缓存所有的输入阵列数据,这在迭加的情况下可以提供约2倍的性能提升。可以使用共享的存储器来缓存阵列多分辨率体积结构。共享的存储器还可以提供用于多谱TERMA计算的寄存器方面的性能改进。最多可以选择21个能量箱,因为它们不仅足够用于高能量射束而且没有存储冲突。
在根据本发明实施方式的TERMA计算和迭加这两个阶段中,与线性光栅化相似(Amanatides,J.,Woo.,A.,Eurographics′87,Conference Proceedings,1987),可以用精确地放射线路径方法代替标准的固定步长射线投射算法。
图10示出由市售的Pinnacle***(Philips Radiation Oncology SystemsMadison WI)和本发明一种实施方式计算出的沿中心轴的吸收剂量图,其中每个吸收剂量图以10cm的深度归一化。图11示出由市售的Pinnacle***和本发明一种实施方式计算出的在10cm的深度的吸收剂量分布图,其中每个吸收剂量分布图在中点处归一化。
对传输算法(比如迭加/卷积)的定量分析可以由对来自用于放射治疗的放射源的源模型处的入射流量的强烈依赖而复杂化(Mohan,R.,Chui,C.,Lidofsky,L.,Med.Phys.12,592-597,1985;McNutt,T.R.:DoseCalculations,Philips white paper)。源模型转而需要被优化,以使测量剂量和计算剂量之间的误差最小。因此,在传输算法中,源模型经常可以掩盖误差。此外,市售***可以包括单独的电子污染模型。在使用了本发明一种实施方式引导进行的初步实验中,仅考虑简单的源模型。初步实验产生了与那些使用市售治疗计划***的模型参数相似的结果,同时提供了如表3中所示的数量级的速度上的改进。
表3.在立方水假体上的剂量引擎性能的比较。速率是每秒钟的体积(VPS)。倾斜的32射线内核比不倾斜的80内核更精确(参见表2)。自适应的多网格方法(未实施)在临床数据上具有2倍的加速因子。
在一些实验中,使用基于来自Varian6EX线性加速器的临床数据的Pinnacle3(Philips-Madison,WI)治疗计划***产生了临床参考数据,提供剂量沉积、TERMA、传输阵列、谱、质量衰减表和单一能量剂量沉积内核。所有实验在AMD处理器254(双核,2.8GHz)上运行。使用标准迭加/卷积引擎,在没有其它程序运行的情况下至少重复十次计时实验。
使用CPU的高性能硬件计数器,将针对本发明实施方式的计时结果重复多次。实验在有单个NVIDIA GeForce GTX280的3GHz的奔腾4上执行。80-射线内核使用了10顶点、用8方位角。32-射线内核使用了4个顶点乘以8个方位角方向。倾斜的72-射线内核使用了6个顶点乘以12个方位角方向,这提供了更高的精确性(参见表3)。总的剂量引擎性能是由源模型、TERMA和迭加性能的总和所确定的,不包括可再次使用的计算,比如修平滤波器结构和TERMA衰减量。
所有试验都在边长为25.6cm的立方水假体上执行。试验在代表标准临床工作量的具有643个体素的体积上、和在具有1283个体素的额外的高分辨率体积上执行。作为本发明一种实施方式的多分辨率迭加方法利用测定体积的mip-map(一种纹理映射技术)并且可以比传统的迭加执行起来快2倍到3倍以及如示出的那样的具有对高分辨率调节得更好的性能。
图12示出根据本发明另一实施方式的一种方法。本方法包括:获得关于身体的预期放射治疗区域的输入信息1201;在方框1202中,在从所述身体到源位置的方向上,基于输入信息1201计算沿射线的第一和第二子区域的物理性质,射线从源位置出发穿过身体并且与身体交叉于预期放射治疗区域,沿射线的第二子区域比沿射线的第一子区域更接近于源位置;以及在方框1203中,确定用于为预期放射治疗区域提供放射治疗的输出信息1204。
在方框1202中,基于在第一子区域处计算的物理性质计算第二子区域处的物理性质。物理性质是指表示从入射的放射束提取的能量的相对数量的衰减因子。如上讨论的衰减因子是一个实例。此外,输入信息1201可以包括来自身体的经验数据。测量数据可以包括身体的预期治疗区域的至少一个电子计算机X射线断层扫描技术(CT)图像。此外,在方框1202中可以使用至少一个预编译的查找表。
方框1203还包括计算被预期放射治疗区域吸收的放射能量的量的步骤1204。使用统一采样、可变采样、多分辨率网格或者其变型中的至少一个以覆盖预期放射治疗区域,可以计算出被吸收的放射能量。此外,多分辨率网格被指定使得来自源位置穿过身体行进的射线可以近似于立体角。
然而,本发明的另一实施方式可以包括计算机可读介质,其包括软件,当该软件被计算机执行时使计算机执行上述方法。
迭加
卷积剂量传输算法基于一个简单的近似:剂量DC,可通过对能量点响应函数(内核)与患者特定能量(TERMA)分布的卷积来计算。
然而,这种近似仅在有非现实平行光束的均匀体积中有效。根据卷积,通过使用放射性距离而不是几何距离索引到内核中,迭加算法得以改进。通过在极坐标中使用径向对称性来计算卷积,此改进成为可能。极坐标变换允许内核倾斜,这更好地模型化投影几何形状,并忽略对距离平方项的需要。通过采样射线的数量ν控制坐标变换的分辨率。根据本发明的一些实施方式,我们实施逆向内核构想,它是工作的并高速缓存有效的,倾斜和非倾斜内核均是,集成了模拟从射线段到点的剂量的累积内核(CK)和集成了模拟从射线段到射线段的剂量的累积-累积内核(CCK)均是。
多能迭加
公式(15)和(16)都含有一个多能近似值。传统上,计算单个多能TERMA量并应用单个多能迭加,而不是计算多重能量特异性TERMA量和迭加各个单能内核来实施。此近似值与人工规定的光谱很好地工作,所述人工规定的光谱比它们的蒙特卡罗(Monte-Carlo)衍生的对应物更难实现。用Monte-Carlo衍生的光谱执行迭加,导致在浅表深度剂量的过高估计。这由于被较高能量内核迭加的低能光子。此外,由于随着深度和楔如光束稳定,内核和TERMA光谱自然分叉。称为内核稳定的技术,通过预稳定多能内核间的线性插值,试图矫正深度依赖的光谱变化。
为解决这些问题,我们向量化了迭加算法。我们表示每一个频谱为累积函数,所述累积函数离散成NT和NS/C,使用单调样条插值均匀地加权能量箱。我们对TERMA、NT和NS/C使用单独的离散,因为TERMA能够使用单个射束跟踪有效计算NT能量箱,同时我们根据本发明的一种实施方式,限定每个射束跟踪迭加到四个能量箱,由于那是GPU内置结构的极限。对于NS/C大于4,我们执行多个通量传输。
我们用NT={16,32,64}且Ns/c={1,2,4,NT}作为例子实验。
异质性补偿迭加
迭加的精度在很大程度上取决于TERMA分布的电子平衡。所述电子平衡受密度和TERMA同步梯度的干扰。这导致接口附近,小型或高度调制域的显著错误。根据经验,这些错误表现为高密度到低密度转变之后的剂量高估和低密度到高密度转变之后的剂量低估。从理论角度来看,发生这些错误,是由于迭加模拟初级电子散射事件,而不是所有的二次散射碰撞。
根据本发明一种实施方式的一种新算法,异质性补偿迭加向迭加算法引入了二次散射碰撞的模型。为了计算放射性距离,我们通过随机游动的近似值过滤,而不是直接取样给定体素的密度。
我们选择使用一个指数函数来近似随机游动。该指数可以从一个定向光源(在我们的例子中,初级散射事件)提供电子剂量沉积的近似值。
也可以使用线性递归滤波器很有效地执行所述指数。传统上,设计线性递归滤波器用于已知固定步长的数字***。当精确的放射学路径射束跟踪产生一个样品间可变的步长,我们将递归指数滤波器的固有衰减常数修饰成dt的函数:
τdt=max(0,1-αdt) (20)
ρeff,t+dt(r′)=(1-τ)ρ(r′)+τdtρeff,t (21)
由于dt是就几何距离而言的,并且电子的随机游动依赖于放射距离,我们采用沉积体素的密度参数化我们的过滤器。由于来自直接围绕沉积点的体素的剂量很少依赖于电子平衡,更多依赖于原始密度,我们增加一个术语β,加强对周边体素密度的影响。此外,为了正确地将周围邻域的影响考虑在内,使用第二递归滤波器和τα值在整个射束跟踪中更新α。我们增加一个角参数F0,以将滤波器的任意可能方向灵敏度考虑在内。由于不平衡量依赖于密度差异,我们增加密度差异术语。
αt+dt=(1-τα,dt)α(r′)+τα,dt,tαt (22)
βt+dt=max(0,βt-B1dt(B0(B2+αt)+B3βt))
α0=α(r)
β0=B0(B2+α0(r))
为了研究这个参数空间,我们使用自动优化和人工图解检查的组合。我们将除了F0以外的所有参数限定到一个小的非负数,其限制在±1/π。我们采用了多种优化技术和目标。对于大部分优化,我们总结了成平方的平均平方误,作为一组示范性的深度剂量曲线。示范实例是ICCR2000、24MV水-肺-水(Water-Lung-Water)和24MV水-低密度-水(Water-Low Density-Water)虚像研究。选择这些虚像研究是因为它们最难同时实现。使用的误差不同于传统的伽玛度量,其使用百分比而不是绝对剂量测定误差。我们使用百分比误差,由于它不减少在深处或低密度材料中的误差权重。我们使用平方误,由于它增加较大误差的权重。我们将误差平均限制到目标区域,3cm至17cm,其在示范实例中包括棘手的材料转换。为了优化,我们最终设置选定多分辨率的贪婪-游动者(greedy-walker)重启,因为它非常有效地工作并提供局部极小值列表,其用于指导人工参数空间探测。然而,本发明的广义精神并不局限于此例。
τα,dt,t=max(0,1-αtdt)
α(r″)=0.5+2ρ(r″)
αt+dt=(1-τα,dt)α(r′)+τα,dt,tαt
βt+dt=max(0,βt(1-1.5dt))
α0=α(r)
β0=0.15α0(r)
在描述本发明的实施方式时,为了清楚起见而采用了特定的术语。然而,本发明并非意在受限于由此而选择的特定术语。本领域的技术人员可想到在不偏离本发明的情况下可根据上述构思对本发明的上述实施方式做出修改或变型。因此可以理解,可以在所附权利要求和与其等同方案的范围内以与所明确描述的方式不同的方式来实施本发明。
Claims (36)
1.一种用于放射疗法的***,包括放射计划***,其中所述放射计划***包括数据处理器,所述数据处理器适于:
接收与身体的预期放射治疗区域有关的信息;
接收计算出的所述身体内多个位置每单位质量释放的初始能量;
基于所述计算出的每单位物质释放的初始能量及包括由于所述身体内的其他位置的分散带来的放射剂量贡献而计算所述放射治疗区域内多个位置的放射剂量;和
基于在所述放射治疗区域内所述多个位置计算的所述放射剂量而确定放射疗法参数用于对所述预期放射治疗区域提供放射治疗,
其中所述包括由于所述身体内的其他位置的分散带来的放射剂量贡献考虑到所述身体内的密度不连续。
2.根据权利要求1所述的用于放射疗法的***,其中在所述计算中考虑到的所述身体内的密度不连续包括在所述身体内不同类型物质之间的界面。
3.根据权利要求2所述的用于放射疗法的***,其中所述身体内不同类型物质之间的所述界面包括骨骼、组织、内部器官、空气、水和人工植入物中的至少两个之间的界面。
4.根据权利要求1所述的用于放射疗法的***,其中所述计算在所述放射治疗区域内所述多个位置中每一个的所述放射剂量包括经内核函数加权的所述身体内的所述多个位置上每单位质量释放的所述初始能量的综合(intergration),
其中所述内核函数包括经验性确定的有效密度函数,目的是将所述身体内的所述密度变化考虑在内。
5.根据权利要求4所述的用于放射疗法的***,其中所述综合针对第一放射光束能量,
其中在所述放射治疗区域内所述多个位置中每一个处的所述计算所述放射剂量还包括针对第二放射光束能量的经内核函数加权的所述身体内的所述多个位置上每单位质量释放的所述初始能量的第二综合,并且
所述计算在所述放射治疗区域内所述多个位置中每一个处的所述放射剂量还包括将所述第一和第二综合加和。
6.根据权利要求1所述的用于放射疗法的***,还包括,在所述接收所述计算出的由所述身体内的多个位置每单位质量释放的初始能量之前,计算要在所述接收中提供的每单位质量释放的所述初始能量。
7.根据权利要求6所述的用于放射疗法的***,其中所述计算每单位质量释放的所述初始能量包括背投影射束追踪计算。
8.根据权利要求5所述的用于放射疗法的***,还包括,在所述接收所述计算出的由所述身体内的多个位置每单位质量释放的初始能量之前,计算要在所述接收中提供的每单位质量释放的所述初始能量。
9.根据权利要求8所述的用于放射疗法的***,其中所述计算每单位质量释放的所述初始能量包括背投影射束追踪计算。
10.根据权利要求1至9任一项所述的用于放射疗法的***,其中每单位质量释放的所述初始能量是TERMA、KERMA或SCERMA中的至少一个。
11.根据权利要求10所述的用于放射疗法的***,其中所述确定所述放射疗法参数对于向所述预期放射治疗区域提供放射治疗足够快以在放射治疗程序期间实时执行。
12.根据权利要求1至9任一项所述的用于放射疗法的***,其中所述确定所述放射疗法参数对于向所述预期放射治疗区域提供放射治疗足够快以在放射治疗程序期间实时执行。
13.一种用于放射治疗的用于确定放射疗法参数的方法,包括:
接收与身体的预期放射治疗区域有关的信息;
接收计算出的所述身体内的多个位置每单位质量释放的初始能量;
基于所述计算出的每单位物质释放的初始能量及包括由于所述身体内的其他位置的分散带来的放射剂量贡献而计算所述放射治疗区域内多个位置的放射剂量;和
基于在所述放射治疗区域内的所述多个位置计算的所述放射剂量而确定放射疗法参数以用于对所述预期放射治疗区域提供放射治疗,
其中所述包括由于所述身体内的其他位置的分散带来的放射剂量贡献考虑到了所述身体内的密度不连续。
14.根据权利要求13所述的用于确定放射疗法参数的方法,其中在所述计算中考虑到的在所述身体内的密度不连续包括在所述身体内不同类型物质之间的界面。
15.根据权利要求14所述的用于确定放射疗法参数的方法,其中所述身体内不同类型物质之间的所述界面包括骨骼、组织、、内部器官、空气、水和、人工植入物中的至少两个之间的界面。
16.根据权利要求13所述的用于确定放射疗法参数的方法,其中所述计算在所述放射治疗区域内的所述多个位置中每一个的所述放射剂量包括经内核函数加权的所述身体内的所述多个位置上每单位质量释放的所述初始能量的综合,
其中所述内核函数包括经验性确定的有效密度函数,目的是将所述身体内的所述密度变化考虑在内。
17.根据权利要求16所述的用于确定放射疗法参数的方法,其中所述综合针对第一放射光束能量,
其中所述计算在所述放射治疗区域内的所述多个位置中每一个的所述放射剂量还包括针对第二放射光束能量的经内核函数加权的所述身体内所述多个位置上每单位质量释放的所述初始能量的第二综合,并且
其中所述计算在所述放射治疗区域内的所述多个位置中每一个的所述放射剂量还包括将所述第一和第二综合加和。
18.根据权利要求13所述的用于确定放射疗法参数的方法,还包括,在所述接收所述计算出的所述身体内的多个位置每单位质量释放的初始能量之前,计算要在所述接收中提供的每单位质量释放的所述初始能量。
19.根据权利要求18所述的用于确定放射疗法参数的方法,其中所述计算每单位质量释放的所述初始能量包括背投影射束追踪计算。
20.根据权利要求17所述的用于确定放射疗法参数的方法,还包括,在所述接收所述计算出的所述身体内多个位置每单位质量释放的初始能量之前,计算要在所述接收中提供的每单位质量释放的所述初始能量。
21.根据权利要求20所述的用于确定放射疗法参数的方法,其中所述计算每单位质量释放的所述初始能量包括背投影射束追踪计算。
22.根据权利要求13至21任一项所述的用于确定放射疗法参数的方法,其中每单位质量释放的所述初始能量是TERMA、KERMA或SCERMA中的至少一个。
23.根据权利要求22所述的用于确定放射疗法参数的方法,其中所述确定所述放射疗法参数对于向所述预期放射治疗区域提供放射治疗足够快以在放射治疗程序期间实时执行。
24.根据权利要求13至21任一项所述的用于确定放射疗法参数的方法,其中所述确定所述放射疗法参数对于向所述预期放射治疗区域提供放射治疗足够快以在放射治疗程序期间实时执行。
25.一种计算机可读介质,包括非临时性可执行代码,所述非临时性可执行代码当由计算机执行时,使得计算机:
接收与身体的预期放射治疗区域有关的信息;
接收计算出的所述身体内多个位置每单位质量释放的初始能量;
基于所述计算出的每单位物质释放的初始能量并包括由于所述身体内的其他位置的分散带来的放射剂量贡献而计算所述放射治疗区域内的多个位置的放射剂量;和
基于在所述放射治疗区域内所述多个位置计算的所述放射剂量而确定放射疗法参数用于对所述预期放射治疗区域提供放射治疗,
其中所述包括由于所述身体内的其他位置的分散带来的放射剂量贡献考虑到所述身体内的密度不连续。
26.根据权利要求25所述的计算机可读介质,其中在所述计算中考虑到的在所述身体内的密度不连续包括在所述身体内不同类型物质之间的界面。
27.根据权利要求26所述的计算机可读介质,其中所述身体内不同类型物质之间的所述界面包括骨骼、组织、内部器官、空气、水和人工植入物中的至少两个之间的界面。
28.根据权利要求25所述的计算机可读介质,其中所述计算在所述放射治疗区域内的所述多个位置中每一个的所述放射剂量包括经内核函数加权的所述身体内的所述多个位置上每单位质量释放的所述初始能量的综合,
其中所述内核函数包括经验性确定的有效密度函数,目的是将所述身体内的所述密度变化考虑在内。
29.根据权利要求28所述的计算机可读介质,其中所述综合针对第一放射光束能量,
其中所述计算在所述放射治疗区域内的所述多个位置中每一个的所述放射剂量还包括针对第二放射光束能量的经内核函数加权的所述身体内所述多个位置上每单位质量释放的所述初始能量的第二综合,并且
所述计算在所述放射治疗区域内所述多个位置中每一个的所述放射剂量还包括将所述第一和第二综合加和。
30.根据权利要求25所述的计算机可读介质,还包括,在所述接收所述计算出的由所述身体内的多个位置每单位质量释放的初始能量之前,计算要在所述接收中提供的每单位质量释放的所述初始能量。
31.根据权利要求30所述的计算机可读介质,其中所述计算每单位质量释放的所述初始能量包括背投影射束追踪计算。
32.根据权利要求29所述的计算机可读介质,还包括,在所述接收所述计算出的由所述身体内的多个位置每单位质量释放的初始能量之前,计算要在所述接收中提供的每单位质量释放的所述初始能量。
33.根据权利要求32所述的计算机可读介质,其中所述计算每单位质量释放的所述初始能量包括背投影射束追踪计算。
34.根据权利要求25至33任一项所述的计算机可读介质,其中每单位质量释放的所述初始能量是TERMA、KERMA或SCERMA中的至少一个。
35.根据权利要求34所述的计算机可读介质,其中所述确定所述放射疗法参数对于向所述预期放射治疗区域提供放射治疗足够快以在放射治疗程序期间实时执行。
36.根据权利要求25至33任一项所述的计算机可读介质,其中所述确定所述放射疗法参数对于向所述预期放射治疗区域提供放射治疗足够快以在放射治疗程序期间实时执行。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161540773P | 2011-09-29 | 2011-09-29 | |
US61/540,773 | 2011-09-29 | ||
PCT/US2012/058345 WO2013049839A2 (en) | 2011-09-29 | 2012-10-01 | Dose computation for radiation therapy using heterogeneity compensated superposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104066479A true CN104066479A (zh) | 2014-09-24 |
CN104066479B CN104066479B (zh) | 2017-06-23 |
Family
ID=47996792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280058918.0A Active CN104066479B (zh) | 2011-09-29 | 2012-10-01 | 使用异质性补偿迭加进行用于放射疗法的剂量计算 |
Country Status (6)
Country | Link |
---|---|
US (1) | US9750955B2 (zh) |
EP (1) | EP2760543B1 (zh) |
JP (1) | JP6138804B2 (zh) |
CN (1) | CN104066479B (zh) |
IN (1) | IN2014KN00902A (zh) |
WO (1) | WO2013049839A2 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105445777A (zh) * | 2015-12-01 | 2016-03-30 | 南京航空航天大学 | 一种硼中子俘获治疗剂量测量方法 |
WO2017054316A1 (zh) * | 2015-09-29 | 2017-04-06 | 广州瑞多思医疗科技有限公司 | 放射治疗中在线剂量监测和验证的方法 |
CN106943679A (zh) * | 2017-04-24 | 2017-07-14 | 安徽慧软科技有限公司 | 基于gpu蒙特卡洛算法的磁场下光子和电子剂量计算方法 |
CN107544086A (zh) * | 2016-06-29 | 2018-01-05 | 长庚医疗财团法人林口长庚纪念医院 | 加马光子侦测成像装置与方法 |
CN111494815A (zh) * | 2020-05-14 | 2020-08-07 | 安徽慧软科技有限公司 | 基于混合变尺度模型的三维剂量计算方法、装置及介质 |
CN111789614A (zh) * | 2020-08-10 | 2020-10-20 | 上海联影医疗科技有限公司 | 一种成像***和方法 |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3460529A1 (en) | 2012-05-29 | 2019-03-27 | Sun Nuclear Corporation | Method and system for calorimetry probe |
US9842424B2 (en) * | 2014-02-10 | 2017-12-12 | Pixar | Volume rendering using adaptive buckets |
US11495355B2 (en) * | 2014-05-15 | 2022-11-08 | The Johns Hopkins University | Method, system and computer-readable media for treatment plan risk analysis |
JP6421194B2 (ja) * | 2014-09-24 | 2018-11-07 | 株式会社日立製作所 | 放射線治療計画装置、放射線治療計画方法および放射線治療システム |
EP3233186B1 (en) | 2014-12-19 | 2018-11-21 | Sun Nuclear Corporation | Radiation therapy dose calculation |
CN104750992B (zh) | 2015-04-01 | 2017-11-28 | 上海联影医疗科技有限公司 | 模拟粒子输运和确定放疗中人体剂量的方法、装置及*** |
US10617891B2 (en) | 2015-04-23 | 2020-04-14 | Sun Nuclear Corporation | Radiation detector calibration |
US10702708B2 (en) * | 2015-09-25 | 2020-07-07 | Varian Medical Systems, Inc. | Accounting for imaging-based radiation doses |
WO2018023049A1 (en) | 2016-07-28 | 2018-02-01 | Sun Nuclear Corporation | Beam angle direction determination |
WO2018115022A1 (en) * | 2016-12-23 | 2018-06-28 | Koninklijke Philips N.V. | Ray tracing for the detection and avoidance of collisions between radiotherapy devices and patient |
WO2018152332A1 (en) * | 2017-02-15 | 2018-08-23 | Henry Ford Health System | Applications of pair-production for improved radiotherapy |
US10918888B2 (en) | 2017-02-28 | 2021-02-16 | Sun Nuclear Corporation | Radiation therapy treatment verification with electronic portal imaging device transit images |
US10493299B2 (en) * | 2017-12-08 | 2019-12-03 | Elekta, Inc. | Determining parameters for a beam model of a radiation machine using deep convolutional neural networks |
CN108815721B (zh) * | 2018-05-18 | 2021-06-25 | 山东省肿瘤防治研究院(山东省肿瘤医院) | 一种照射剂量确定方法及*** |
US11278744B2 (en) | 2018-09-28 | 2022-03-22 | Sun Nuclear Corporation | Systems and methods to account for tilt of a radiation measurement system |
CN110215623B (zh) * | 2019-06-11 | 2020-11-10 | 苏州雷泰智能科技有限公司 | 基于子野优化的正交双层光栅旋转调强的实现方法及装置 |
WO2021007459A1 (en) | 2019-07-10 | 2021-01-14 | Sun Nuclear Corporation | Scintillator-based radiation therapy quality assurance |
US12011616B2 (en) | 2019-07-10 | 2024-06-18 | Sun Nuclear Corporation | Image-based radiation therapy quality assurance |
US11600004B2 (en) | 2019-07-10 | 2023-03-07 | Sun Nuclear Corporation | Image-based radiation therapy quality assurance |
WO2021155342A1 (en) * | 2020-01-30 | 2021-08-05 | Xie George Xu | Patient-specific organ dose quantification and inverse optimization for ct |
US11559701B2 (en) | 2021-06-17 | 2023-01-24 | Varian Medical Systems International Ag | Method and apparatus to facilitate administering therapeutic radiation to a heterogeneous body |
JP7345768B1 (ja) | 2022-03-30 | 2023-09-19 | 国立大学法人大阪大学 | 予測支援システム、予測支援方法及び予測支援プログラム |
US20240100360A1 (en) * | 2022-09-28 | 2024-03-28 | Siemens Healthineers International Ag | Radiation treatment plan optimization apparatus and method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005052721A2 (en) * | 2003-09-24 | 2005-06-09 | Radion Technologies | Deterministic computation of radiation doses delivered to tissues and organs of a living organism |
WO2008013956A1 (en) * | 2006-07-26 | 2008-01-31 | Nomos Corporation | System for enhancing intensity modulated radiation therapy, program product, and related methods |
WO2009137794A2 (en) * | 2008-05-08 | 2009-11-12 | The Johns Hopkins University | Real-time dose computation for radiation therapy using graphics processing unit acceleration of the convolution/superposition dose computation method |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5602892A (en) | 1996-03-21 | 1997-02-11 | Llacer; Jorge | Method for optimization of radiation therapy planning |
US6636622B2 (en) * | 1997-10-15 | 2003-10-21 | Wisconsin Alumni Research Foundation | Method and apparatus for calibration of radiation therapy equipment and verification of radiation treatment |
EP1055181B1 (en) | 1998-02-09 | 2002-05-22 | University Of Southampton | Treatment planning method and apparatus for radiation therapy |
US7289599B2 (en) * | 2002-10-04 | 2007-10-30 | Varian Medical Systems Technologies, Inc. | Radiation process and apparatus |
US20080091388A1 (en) * | 2003-03-14 | 2008-04-17 | Failla Gregory A | Method for calculation radiation doses from acquired image data |
US7202486B2 (en) * | 2004-08-04 | 2007-04-10 | Standard Imaging, Inc. | Treatment planning tool for multi-energy electron beam radiotherapy |
WO2006138513A1 (en) | 2005-06-16 | 2006-12-28 | Nomos Corporation | Variance reduction simulation system, program product, and related methods |
WO2007062178A2 (en) * | 2005-11-21 | 2007-05-31 | The Regents Of The University Of California | Method for computing patient radiation dose in computed tomoraphy |
WO2007106815A2 (en) * | 2006-03-14 | 2007-09-20 | University Of Notre Dame Du Lac | Methods and apparatus for hardware based radiation dose calculation |
US8483458B2 (en) * | 2009-09-10 | 2013-07-09 | General Electric Company | Method and system for measuring visceral fat mass using dual energy x-ray absorptiometry |
WO2013024534A1 (ja) * | 2011-08-17 | 2013-02-21 | 三菱電機株式会社 | 皮膚線量評価支援装置及び治療計画装置 |
-
2012
- 2012-10-01 CN CN201280058918.0A patent/CN104066479B/zh active Active
- 2012-10-01 JP JP2014533482A patent/JP6138804B2/ja active Active
- 2012-10-01 US US14/346,949 patent/US9750955B2/en active Active
- 2012-10-01 WO PCT/US2012/058345 patent/WO2013049839A2/en active Application Filing
- 2012-10-01 EP EP12834789.5A patent/EP2760543B1/en active Active
- 2012-10-01 IN IN902KON2014 patent/IN2014KN00902A/en unknown
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005052721A2 (en) * | 2003-09-24 | 2005-06-09 | Radion Technologies | Deterministic computation of radiation doses delivered to tissues and organs of a living organism |
WO2008013956A1 (en) * | 2006-07-26 | 2008-01-31 | Nomos Corporation | System for enhancing intensity modulated radiation therapy, program product, and related methods |
WO2009137794A2 (en) * | 2008-05-08 | 2009-11-12 | The Johns Hopkins University | Real-time dose computation for radiation therapy using graphics processing unit acceleration of the convolution/superposition dose computation method |
CN102036712A (zh) * | 2008-05-08 | 2011-04-27 | 约翰霍普金斯大学 | 卷积/迭加剂量计算方法的、使用图形处理单元加速的用于放射治疗的实时剂量计算 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017054316A1 (zh) * | 2015-09-29 | 2017-04-06 | 广州瑞多思医疗科技有限公司 | 放射治疗中在线剂量监测和验证的方法 |
CN105445777A (zh) * | 2015-12-01 | 2016-03-30 | 南京航空航天大学 | 一种硼中子俘获治疗剂量测量方法 |
CN105445777B (zh) * | 2015-12-01 | 2018-06-05 | 南京航空航天大学 | 一种硼中子俘获治疗剂量测量方法 |
CN107544086A (zh) * | 2016-06-29 | 2018-01-05 | 长庚医疗财团法人林口长庚纪念医院 | 加马光子侦测成像装置与方法 |
CN107544086B (zh) * | 2016-06-29 | 2020-03-10 | 长庚医疗财团法人林口长庚纪念医院 | 加马光子侦测成像装置与方法 |
CN106943679A (zh) * | 2017-04-24 | 2017-07-14 | 安徽慧软科技有限公司 | 基于gpu蒙特卡洛算法的磁场下光子和电子剂量计算方法 |
CN111494815A (zh) * | 2020-05-14 | 2020-08-07 | 安徽慧软科技有限公司 | 基于混合变尺度模型的三维剂量计算方法、装置及介质 |
CN111494815B (zh) * | 2020-05-14 | 2022-04-29 | 安徽慧软科技有限公司 | 基于混合变尺度模型的三维剂量计算方法、装置及介质 |
CN111789614A (zh) * | 2020-08-10 | 2020-10-20 | 上海联影医疗科技有限公司 | 一种成像***和方法 |
Also Published As
Publication number | Publication date |
---|---|
US9750955B2 (en) | 2017-09-05 |
IN2014KN00902A (zh) | 2015-10-09 |
CN104066479B (zh) | 2017-06-23 |
EP2760543A2 (en) | 2014-08-06 |
JP2014528797A (ja) | 2014-10-30 |
WO2013049839A3 (en) | 2013-05-30 |
EP2760543A4 (en) | 2015-04-08 |
JP6138804B2 (ja) | 2017-05-31 |
US20140235923A1 (en) | 2014-08-21 |
EP2760543B1 (en) | 2021-05-19 |
WO2013049839A2 (en) | 2013-04-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102036712B (zh) | 用于放射治疗的***和确定放射治疗参数的方法 | |
CN104066479A (zh) | 使用异质性补偿迭加进行用于放射疗法的剂量计算 | |
US7856082B2 (en) | System and method for optimization of a radiation therapy plan in the presence of motion | |
Reynaert et al. | Monte Carlo treatment planning for photon and electron beams | |
Jacques et al. | Towards real-time radiation therapy: GPU accelerated superposition/convolution | |
US20080091388A1 (en) | Method for calculation radiation doses from acquired image data | |
US20110122997A1 (en) | Non-voxel-based broad-beam (nvbb) algorithm for intensity modulated radiation therapy dose calculation and plan optimization | |
US9495513B2 (en) | GPU-based fast dose calculator for cancer therapy | |
Wang et al. | Development of a fast Monte Carlo dose calculation system for online adaptive radiation therapy quality assurance | |
Papanikolaou et al. | Dose‐calculation algorithms in the context of inhomogeneity corrections for high energy photon beams | |
US9251302B2 (en) | System and method for planning a radiation therapy treatment | |
Moskvin et al. | PENELOPE Monte Carlo engine for treatment planning in radiation therapy with very high energy electrons (VHEE) of 150–250 MeV | |
Dickmann et al. | Fluence-modulated proton CT optimized with patient-specific dose and variance objectives for proton dose calculation | |
CN110675932A (zh) | 一种基于蒙特卡罗的点剂量计算方法、设备和存储介质 | |
Muhammad et al. | Tracking organ doses for patient safety in radiation therapy | |
Jacques et al. | GPU accelerated real time kV/MV dose computation | |
Scholz | Development and evaluation of advanced dose calculations for modern radiation therapy techniques | |
Chatzipapas et al. | Standardization and Validation of Brachytherapy Seeds’ Modelling Using GATE and GGEMS Monte Carlo Toolkits. Cancers 2021, 13, 5315 | |
Muhammad et al. | Tracking 9 therapy organ doses for patient safety in radiation | |
Papanikolaou | Clinical photon beam treatment planning using convolution and superposition | |
Muhammad et al. | 9 Tracking doses for safety in radiation therapy | |
Rai | Image quality measures in proton computed tomography | |
Tamimi | Development of a Forward/Adjoint Hybrid Monte Carlo Absorbed Dose Calculational Method for Use in Radiation Therapy | |
CN116350959A (zh) | 放射治疗剂量分析方法和*** | |
Wang | Fast Monte Carlo Simulations for Quality Assurance in Radiation Therapy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |