A multi-FPGA accelerator for radiation dose calculation in cancer treatment

Radiation therapy (RT) is a major modern modality for cancer treatment. However, the use of a large number of radiation beams in treatment plans and image-guided on-line replanning in newly emerging RT planning and delivery systems present a daunting challenge to radiation dose calculation even on the state-of-the-art quad-core computers as the dose calculation time becomes unacceptably long for clinical applications. Based on a popular dose calculation algorithm, the Collapsed-Cone Convolution/Superposition (CCCS) algorithm, this paper presents a multi-FPGA accelerator design for the dose calculation problem. Our performance-driven design strategy yields a fully pipelined architecture, which includes a resource-economic raytracing engine and high-performance energy deposition pipeline. The design is capable of processing 3 transporting lines simultaneously and accomplishing one dose deposition in every 2 cycles while raytracing along the test phantom. The evaluation based on a set of clinical treatment planning cases confirms that our FPGA design almost fully utilizes the available external memory bandwidth and achieves close to the best possible performance for the CCCS algorithm while using less resource. With the support of a memory-rich multi-FPGA platform, our floating-point based design working at 90Mhz obtains a speedup of 20X over a commercial multi-threaded software on a quad-core system and 15X performance improvement over the closely related results.

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

[2]  M. Aspradakis,et al.  Experimental verification of convolution/superposition photon dose calculations for radiotherapy treatment planning. , 2003, Physics in medicine and biology.

[3]  T. Krieger,et al.  Monte Carlo- versus pencil-beam-/collapsed-cone-dose calculation in a heterogeneous multi-layer phantom , 2005, Physics in medicine and biology.

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

[5]  An interactive, parallel, three-dimensional fast Fourier transform convolution dose calculation using a supercomputer. , 1997, Medical physics.

[6]  Cedric X. Yu,et al.  Arc-modulated radiation therapy (AMRT): a single-arc form of intensity-modulated arc therapy , 2008, Physics in medicine and biology.

[7]  B. Fraass,et al.  Summary and recommendations of a National Cancer Institute workshop on issues limiting the clinical use of Monte Carlo dose calculation algorithms for megavoltage external beam radiation therapy. , 2003, Medical physics.

[8]  Lothar Lilge,et al.  FPGA-based Monte Carlo Computation of Light Absorption for Photodynamic Cancer Therapy , 2009, 2009 17th IEEE Symposium on Field Programmable Custom Computing Machines.

[9]  Maria Grazia Pia,et al.  Monte Carlo dose calculation algorithm on a distributed system , 2003 .

[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]  Xiaobo Sharon Hu,et al.  Hardware Acceleration for 3-D Radiation Dose Calculation , 2007, 2007 IEEE International Conf. on Application-specific Systems, Architectures and Processors (ASAP).

[12]  Xiaobo Sharon Hu,et al.  An FPGA Solution for Radiation Dose Calculation , 2006, 2006 14th Annual IEEE Symposium on Field-Programmable Custom Computing Machines.

[13]  Peter Dunscombe,et al.  Experimental validation of the Eclipse AAA algorithm , 2007, Journal of applied clinical medical physics.