Dose-volume histograms computation comparisons using conventional methods and optimized fast Fourier transforms algorithms for brachytherapy.

In anatomy based optimization procedures for large volume implants the calculation of dose-volume histograms (DVH) accounts for the major part of the time involved and can be as long as a few hours. This time is proportional to the number of seeds or source dwell positions required for the implant. A procedure for the calculation of brachytherapy seed dose distribution calculation employing fast Fourier transforms (FFT) and the convolution theorem has been described by others and was supposed to significantly improve the speed of the dose distribution computation. Using new significantly improved FFT algorithms and various other optimization techniques we have compared the calculated differential and integral DVHs in high dose rate (HDR) brachytherapy with a single stepping source using actual clinical implants. This is so that we could assess the efficiency and accuracy of the FFT method with that of conventional methods. Our results showed that the FFT based method of calculating DVHs in brachytherapy is comparable in speed with conventional dose calculation methods, but only for implants with more than 287 sources. It is therefore of limited practical use even for large implants. This result is in direct opposition to the claim by other authors.

[1]  R. Hartley.  A More Symmetrical Fourier Analysis Applied to Transmission Problems , 1942, Proceedings of the IRE.

[2]  D Baltas,et al.  Anatomy-based three-dimensional dose optimization in brachytherapy using multiobjective genetic algorithms. , 1999, Medical physics.

[3]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[4]  C. Sidney Burrus,et al.  Automatic generation of prime length FFT programs , 1996, IEEE Trans. Signal Process..

[5]  M H Loew,et al.  Use of the fast Hartley transform for three-dimensional dose calculation in radionuclide therapy. , 1998, Medical physics.

[6]  D A Low,et al.  A convolution-adapted ratio-TAR algorithm for 3D photon beam treatment planning. , 1995, Medical physics.

[7]  R Nath,et al.  Permanent prostate seed implant brachytherapy: report of the American Association of Physicists in Medicine Task Group No. 64. , 1999, Medical physics.

[8]  Steven G. Johnson,et al.  The Fastest Fourier Transform in the West , 1997 .

[9]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[10]  Clive Temperton,et al.  A Generalized Prime Factor FFT Algorithm for any N = 2p 3q 5r , 1992, SIAM J. Sci. Comput..

[11]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .

[12]  R. M. Sievert.  Die Intensitatsverteilung der primaren γ-Strahlung in der Nahe medizinischer Radiumpraparate , 1921 .

[13]  C. E. SHANNON,et al.  A mathematical theory of communication , 1948, MOCO.

[14]  L. Anderson,et al.  Dosimetry of interstitial brachytherapy sources: Recommendations of the AAPM Radiation Therapy Committee Task Group No. 43 , 1995 .

[15]  E. Mok,et al.  Brachytherapy seed dose distribution calculation employing the fast Fourier transform. , 1986, Medical physics.