CN101968830A - 利用gpu加速确定电子束剂量的方法 - Google Patents

利用gpu加速确定电子束剂量的方法 Download PDF

Info

Publication number
CN101968830A
CN101968830A CN 201010296330 CN201010296330A CN101968830A CN 101968830 A CN101968830 A CN 101968830A CN 201010296330 CN201010296330 CN 201010296330 CN 201010296330 A CN201010296330 A CN 201010296330A CN 101968830 A CN101968830 A CN 101968830A
Authority
CN
China
Prior art keywords
gpu
electron beam
dose
determining
acceleration
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.)
Pending
Application number
CN 201010296330
Other languages
English (en)
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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN 201010296330 priority Critical patent/CN101968830A/zh
Publication of CN101968830A publication Critical patent/CN101968830A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radiation-Therapy Devices (AREA)

Abstract

本发明公开了一种利用GPU加速确定电子束剂量的方法,目的是提供一种加速方法,使基于CPU平台的确定电子束剂量方法在CPU+GPU异构平台上获得加速。包含以下步骤:1)在CPU平台上形成电子束照射野的强度矩阵和加速器电子束的参数;2)选定用于确定电子束剂量的GPU,为GPU分配存储器空间,向GPU传送计算数据;3)利用GPU并行计算三维剂量矩阵的每个单元的剂量。依照本发明确定电子束在均匀体模中的三维剂量分布,计算速度快、计算精度与基于CPU平台的确定电子束剂量分布的方法相当。

Description

利用GPU加速确定电子束剂量的方法 
技术领域
本发明涉及一种确定电子束剂量分布的方法,具体的说,是用于在电子束外照射放射治疗中,快速、准确地确定电子线在被照射的均匀体模中的三维剂量分布的方法。 
背景技术
目前比较成熟的肿瘤治疗和控制手段包括手术治疗、放射治疗和化疗。其中放射治疗占了相当大的比例。按照源的位置的不同,放射治疗分为外照射放射治疗和内照射放射治疗。按照源的种类,外照射放射治疗包括X线放射治疗、电子线放射治疗、质子放射治疗和重粒子放射治疗。其中电子线放射治疗对于浅表肿瘤具有一定的优势,受到临床应用的关注。 
在电子线放射治疗的临床应用上,快速、准确地确定电子束在人体或体模中的三维剂量分布是一个基本的问题。目前使用的确定电子束剂量的算法,包括Monte Carlo方法和解析算法。前者包括MMC(Macro Monte Carlo)方法、VMC(Voxel Monte Carlo)方法,后者包括PBRA(Electron Pencil-Beam Redefinition Algorithm)和HPBM(Hybrid Pencil Beam Model,发明专利:确定电子束剂量分布的方法,ZL200610022709.7)等。 
GPU(Graphic Processing Unit)最初应用于图形显示的加速。GPU的单指令多数据流(SIMD:Single Instruction Multiple Data)的处理方式可并行地对大规模的数据进行操作,可大大缩短计算的时间。 
本发明将GPU应用于HPBM的剂量计算。通过同时进行大量剂量矩阵单元的计算,其计算精度和使用CPU(Central Process Unit)的计算精度相当,其计算时间较使用CPU平台计算大大减少。 
发明内容
本发明的目的在于提供一种能够快速、准确地确定电子束在均匀体模中的三维剂量分布的方法。 
本发明的技术方案是: 
一种利用GPU加速确定电子束剂量分布的方法,其特征在于包含以下步骤: 
1、在CPU平台上形成电子束照射野的强度矩阵I=Ii,j和加速器电子束的参数。电子束参数包括:电子束能量(E)、SSD(电子束照射野的源皮距)、测量的深度剂量(DM(z))、不同深度电子束的散射参数A(z)。 
2、选定用于确定电子束剂量的GPU,为GPU分配存储器空间,向GPU传送计算数据。需要传送的数据包括强度矩阵、测量的深度剂量、散射参数; 
3、利用GPU并行计算三维剂量矩阵的每个单元的剂量分布。对于给定单元的剂量分布,是由所有有限电子笔束剂量的加权累积得到,权重为强度矩阵。 
能量为E的有限电子笔束(笔束的截面为有限大,截面大小为2a×2b)垂直入射在介质表面。以电子束入射方向为深度方向(z轴正方向),以入射点为坐标原点。该电子束在介质中产生的三维剂量分布为: 
D P ( x , y , z , E ) = D M ( z ) 4 * [ erf ( a + x A ( z ) ) + erf ( a - x A ( z ) ) ] * [ erf ( b + x A ( z ) ) + erf ( b - x A ( z ) ) ]
Dp(x,y,z,E)表示能量为E的有限电子笔束在空间点P(x,y,z)处沉积的能量(即介质在空间点P的吸收剂量),a和b分别是有限电子笔束的束截面的长度和宽度的一半。 
为了利用GPU的并行数据处理功能,加快计算速度,对三维剂量矩阵的每个单元分配一个线程进行计算。对于给定的单元V(i,j,k),其剂量D(i,j,k)是由所有有限电子笔束剂量的加权累积得到,权重为强度矩阵,其剂量为: 
D ( i , j , k ) = Σ m = LX HX Σ n = LY HY I ( m , n ) SSD 2 ( SSD + z k ) 2 D p ( x m , y n , z k )
式中,SSD(Source Surface Distance)是电子束照射野的源皮距,i,j,k分别是该单元在三维剂量矩阵的x,y,z方向的编号,I(m,n)是强度矩阵的(m,n)单元的值,即第(m,n)有限电子笔束的强度,m,n分别是强度矩阵单元在x,y方向上的编号,x方向上的强度矩阵的范围为[LX,HX],y方向强度矩阵的范围为[LY,HY],xm,yn分别是第(m,n)有限电子笔束的束截面的中心点与体元V(i,j,k)的中心点在x,y方向上的距离,zk是体元V(i,j,k)的深度, 
xm=(m-i)*δx
yn=(n-j)*δy
zk=k*δz
δx,δy分别是强度矩阵和三维剂量矩阵在x,y方向上的离散化步长,δz是三维剂量矩阵在z方向上的离散化步长。 
本发明的有益效果是:和基于CPU平台确定电子线三维剂量分布的方法相比,在保证计算精度的前提下,可以大大减少计算时间,对于10cm×10cm的电子束照射野,计算时间可以减少约28倍。 
附图说明
图1是利用GPU加速确定电子束剂量的流程图。 
图2是本发明在10cm×10cm的照射野,SSD=100cm,E=20MeV的照射条件下,在15cm×15cm×15cm的均匀水箱的中心平面(y=0)上,使用CPU平台和CPU+GPU异构平台计算得到的中心平面的二维剂量分布的比较。圆点为基于CPU平台计算得到的结果,实线为使用基于CPU+GPU的异构平台上得到的计算结果。 
具体实施方式
实例一:E=20MeV、10cm×10cm照射野、SSD=100cm的电子线照射 
本实例的测量数据(DM(z))引用ECWG(高能电子束治疗计划联合工作组:Collaborative Working Group contract on high energy electron beam treatment planning)的实验测量数据,散射参数A的计算使用发明专利“确定电子束剂量分布的方法”中使用的方法(其中定义为A2(z,E))。本实例计算了在SSD=100cm、10cm×10cm射野的照射条件下,20MeV电子束照射在均匀水箱中的剂量分布,如图2所示。图中给出了在中心平面y=0的等剂量曲线的比较,圆点为只使用CPU计算得到的结果,实线为使用GPU得到的计算结果。在CPU( 
Figure BSA00000289105000031
Dual-core 2.80GHz)平台上运行时,需要时间为6.157秒;在CPU+GPU(型号:GeForce GTX 470)的异构平台上运行时,需要时间是0.219秒。 
上述实施例说明使用本发明可以在和基于CPU平台的计算相当的精度下,减少确定剂量所需要的时间。 

Claims (1)

1.一种利用GPU加速确定电子束剂量的方法,其特征在于,包含以下步骤:
(1)对于电子束照射野,采用强度矩阵描述,强度矩阵的每个单元对应一个有限电子笔束;
(2)对于每个有限电子笔束在三维体模中的剂量分布,采用HPBM(混合笔束模型)计算,计算在GPU上完成;
(3)对于整个照射野在三维体模中的剂量分布,是由有限电子笔束剂量的加权累积得到,权重为强度矩阵,该加权累积过程在GPU上完成。
CN 201010296330 2010-09-29 2010-09-29 利用gpu加速确定电子束剂量的方法 Pending CN101968830A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010296330 CN101968830A (zh) 2010-09-29 2010-09-29 利用gpu加速确定电子束剂量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010296330 CN101968830A (zh) 2010-09-29 2010-09-29 利用gpu加速确定电子束剂量的方法

Publications (1)

Publication Number Publication Date
CN101968830A true CN101968830A (zh) 2011-02-09

Family

ID=43547986

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010296330 Pending CN101968830A (zh) 2010-09-29 2010-09-29 利用gpu加速确定电子束剂量的方法

Country Status (1)

Country Link
CN (1) CN101968830A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103425875A (zh) * 2013-07-29 2013-12-04 南京航空航天大学 一种光子束在目标体内剂量分布的快速获取方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100998497A (zh) * 2006-12-29 2007-07-18 四川大学 确定电子束剂量分布的方法
CN101120871A (zh) * 2006-12-29 2008-02-13 成都川大奇林科技有限责任公司 精确放射治疗计划系统
CN101664583A (zh) * 2009-09-09 2010-03-10 深圳市海博科技有限公司 基于cuda的剂量计算优化方法和系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100998497A (zh) * 2006-12-29 2007-07-18 四川大学 确定电子束剂量分布的方法
CN101120871A (zh) * 2006-12-29 2008-02-13 成都川大奇林科技有限责任公司 精确放射治疗计划系统
CN101664583A (zh) * 2009-09-09 2010-03-10 深圳市海博科技有限公司 基于cuda的剂量计算优化方法和系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103425875A (zh) * 2013-07-29 2013-12-04 南京航空航天大学 一种光子束在目标体内剂量分布的快速获取方法

Similar Documents

Publication Publication Date Title
Yepes et al. A GPU implementation of a track-repeating algorithm for proton radiotherapy dose calculations
CN100431485C (zh) 确定电子束剂量分布的方法
Jia et al. GPU-based fast Monte Carlo dose calculation for proton therapy
Kohno et al. Clinical implementation of a GPU-based simplified Monte Carlo method for a treatment planning system of proton beam therapy
CN107073286B (zh) 放射线治疗计划装置、放射线治疗计划方法以及放射线治疗系统
CN103405852A (zh) 立体定向放射治疗剂量场分布的快速计算方法
CN107708803B (zh) 剂量分布运算装置、粒子射线治疗装置及剂量分布运算方法
CN103127623B (zh) 一种放射治疗中在线验证加速器出束准确性的方法
Da Silva et al. Sub-second pencil beam dose calculation on GPU for adaptive proton therapy
Fujimoto et al. GPU-based fast pencil beam algorithm for proton therapy
CN101954148A (zh) 放射治疗中基于gpu的剂量计算加速方法
WO2017081768A1 (ja) 粒子線治療装置及び治療計画補正方法
Plante et al. Calculation of the energy deposition in nanovolumes by protons and HZE particles: geometric patterns of initial distributions of DNA repair foci
Debus et al. A point kernel algorithm for microbeam radiation therapy
CN110310720A (zh) 一种蒙特卡罗剂量计算方法、设备和储存介质
US9251302B2 (en) System and method for planning a radiation therapy treatment
Egashira et al. Application of the pencil-beam redefinition algorithm in heterogeneous media for proton beam therapy
CN102201036A (zh) Gpu加速剂量计算中微分卷积积分算法的实现
CN103425875A (zh) 一种光子束在目标体内剂量分布的快速获取方法
CN104298871B (zh) 一种辐射防护中基于gpu加速的体素人体模型剂量评估加速方法
CN102426377A (zh) 一种三维剂量反演方法
US9486644B2 (en) Method and system for dose determination of radiation therapy
CN107073288B (zh) 照射计划装置以及照射计划校正方法
CN101968830A (zh) 利用gpu加速确定电子束剂量的方法
CN107290774A (zh) 中子剂量快速确定方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20110209