Implementation of FFT convolution and multigrid superposition models in the FOCUS RTP system.

In radiotherapy treatment planning, convolution/superposition algorithms currently represent the best practical approach for accurate photon dose calculation in heterogeneous tissues. In this work, the implementation, accuracy and performance of the FFT convolution (FFTC) and multigrid superposition (MGS) algorithms are presented. The FFTC and MGS models use the same 'TERMA' calculation and are commissioned using the same parameters. Both models use the same spectra, incorporate the same off-axis softening and base incident lateral fluence on the same measurements. In addition, corrections are explicitly applied to the polyenergetic and parallel kernel approximations, and electron contamination is modelled. Spectra generated by Monte Carlo (MC) modelling of treatment heads are used. Calculations using the MC spectra were in excellent agreement with measurements for many linear accelerator types. To speed up the calculations, a number of calculation techniques were implemented, including separate primary and scatter dose calculation, the FFT technique which assumes kernel invariance for the convolution calculation and a multigrid (MG) acceleration technique for the superposition calculation. Timing results show that the FFTC model is faster than MGS by a factor of 4 and 8 for small and large field sizes, respectively. Comparisons with measured data and BEAM MC results for a wide range of clinical beam setups show that (a) FFTC and MGS doses match measurements to better than 2% or 2 mm in homogeneous media; (b) MGS is more accurate than FFTC in lung phantoms where MGS doses are within 3% or 3 mm of BEAM results and (c) FFTC overestimates the dose in lung by a maximum of 9% compared to BEAM.

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

[2]  P Andreo,et al.  Determination of effective bremsstrahlung spectra and electron contamination for photon dose calculations. , 1989, Physics in medicine and biology.

[3]  C. Ma,et al.  BEAM: a Monte Carlo code to simulate radiotherapy treatment units. , 1995, Medical physics.

[4]  N Papanikolaou,et al.  Investigation of the convolution method for polyenergetic spectra. , 1993, Medical physics.

[5]  A Ahnesjö,et al.  Analytic modeling of photon scatter from flattening filters in photon therapy beams. , 1994, Medical physics.

[6]  A. Boyer,et al.  Fast Fourier transform convolution calculations of x-ray isodose distributions in homogeneous media. , 1989, Medical physics.

[7]  D A Jaffray,et al.  Extrafocal radiation: a unified approach to the prediction of beam penumbra and output factors for megavoltage x-ray beams. , 1995, Medical physics.

[8]  K L Lam,et al.  Flattening-filter-based empirical methods to parametrize the head scatter factor. , 1996, Medical physics.

[9]  E. Mok,et al.  A photon dose distribution model employing convolution calculations. , 1985, Medical physics.

[10]  J J Battista,et al.  Generation of photon energy deposition kernels using the EGS Monte Carlo code. , 1988, Physics in medicine and biology.

[11]  M. Karlsson,et al.  Electron contamination in clinical high energy photon beams. , 1996, Medical physics.

[12]  J. E. O'Connor.  The variation of scattered x-rays with density in an irradiated body. , 1957, Physics in medicine and biology.

[13]  R Mohan,et al.  Energy and angular distributions of photons from medical linear accelerators. , 1985, Medical physics.

[14]  P Keall,et al.  Accounting for primary electron scatter in x-ray beam convolution calculations. , 1995, Medical physics.

[15]  A Fenster,et al.  X-ray sources of medical linear accelerators: focal and extra-focal radiation. , 1993, Medical physics.

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

[17]  Leonard Eyges,et al.  Multiple Scattering with Energy Loss , 1948 .

[18]  E L Chaney,et al.  A Monte Carlo study of accelerator head scatter. , 1994, Medical physics.

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

[20]  Y. Zhu,et al.  Accuracy requirements of the primary x-ray spectrum in dose calculations using FFT convolution techniques. , 1995, Medical physics.

[21]  R Mohan,et al.  Investigation of buildup dose from electron contamination of clinical photon beams. , 1983, Medical physics.

[22]  A. Boyer,et al.  X-ray dose computations in heterogeneous media using 3-dimensional FFT convolution , 1990 .

[23]  J J Battista,et al.  Dose calculations using convolution and superposition principles: the orientation of dose spread kernels in divergent x-ray beams. , 1993, Medical physics.