Fluence-convolution broad-beam (FCBB) dose calculation

IMRT optimization requires a fast yet relatively accurate algorithm to calculate the iteration dose with small memory demand. In this paper, we present a dose calculation algorithm that approaches these goals. By decomposing the infinitesimal pencil beam (IPB) kernel into the central axis (CAX) component and lateral spread function (LSF) and taking the beam's eye view (BEV), we established a non-voxel and non-beamlet-based dose calculation formula. Both LSF and CAX are determined by a commissioning procedure using the collapsed-cone convolution/superposition (CCCS) method as the standard dose engine. The proposed dose calculation involves a 2D convolution of a fluence map with LSF followed by ray tracing based on the CAX lookup table with radiological distance and divergence correction, resulting in complexity of O(N(3)) both spatially and temporally. This simple algorithm is orders of magnitude faster than the CCCS method. Without pre-calculation of beamlets, its implementation is also orders of magnitude smaller than the conventional voxel-based beamlet-superposition (VBS) approach. We compared the presented algorithm with the CCCS method using simulated and clinical cases. The agreement was generally within 3% for a homogeneous phantom and 5% for heterogeneous and clinical cases. Combined with the 'adaptive full dose correction', the algorithm is well suitable for calculating the iteration dose during IMRT optimization.

[1]  N. Papanikolaou,et al.  Dose-calculation algorithms in the context of inhomogeneity corrections for high energy photon beams. , 2009, Medical physics.

[2]  A. Ahnesjö Collapsed cone convolution of radiant energy for photon dose calculation in heterogeneous media. , 1989, Medical physics.

[3]  B. Bjärngard,et al.  Equivalent fields and scatter integration for photon fields. , 1999, Physics in medicine and biology.

[4]  R. Mohan,et al.  Use of fast Fourier transforms in calculating dose distributions for irregularly shaped fields for three-dimensional treatment planning. , 1987, Medical physics.

[5]  Weiguo Lu,et al.  Accurate convolution/superposition for multi-resolution dose calculation using cumulative tabulated kernels , 2005, Physics in medicine and biology.

[6]  A. Ahnesjö,et al.  Dose calculations for external photon beams in radiotherapy. , 1999, Physics in medicine and biology.

[7]  J D Bourland,et al.  A finite-size pencil beam model for photon dose calculations in three dimensions. , 1992, Medical physics.

[8]  T. Bortfeld,et al.  Decomposition of pencil beam kernels for fast dose calculations in three-dimensional treatment planning. , 1993, Medical physics.

[9]  A. Ahnesjö,et al.  A pencil beam model for photon dose calculation. , 1992, Medical physics.

[10]  Weiguo Lu,et al.  Motion-encoded dose calculation through fluence/sinogram modification. , 2005, Medical physics.

[11]  H Hansson,et al.  Verification of a pencil beam based treatment planning system: output factors for open photon beams shaped with MLC or blocks. , 1999, Physics in medicine and biology.

[12]  J. Battista,et al.  A convolution method of calculating dose for 15-MV x rays. , 1985, Medical physics.

[13]  R. Mohan,et al.  Differential pencil beam dose computation model for photons. , 1986, Medical physics.

[14]  R. Mohan,et al.  Reducing dose calculation time for accurate iterative IMRT planning. , 2002, Medical physics.

[15]  Indrin J Chetty,et al.  A fluence convolution method to account for respiratory motion in three-dimensional dose calculations of the liver: a Monte Carlo study. , 2003, Medical physics.

[16]  T. Holmes,et al.  A filtered backprojection dose calculation method for inverse treatment planning. , 1994, Medical physics.

[17]  Weiguo Lu,et al.  A non-voxel-based broad-beam (NVBB) framework for IMRT treatment planning , 2010, Physics in medicine and biology.

[18]  Mingli Chen,et al.  Generalized equivalent field size for nonuniform fluence maps in IMRT dose calculation. , 2010, Medical physics.