Cone-beam computed tomography (CT) for a megavoltage linear accelerator (LINAC) using an electronic portal imaging device (EPID) and the algebraic reconstruction technique (ART)

This study investigates the feasibility and applicability of cone beam CT for a megavoltage therapeutic LINAC. A Rando head phantom was irradiated using the 15 MV beam of a Varian 2100C LINAC at MSKCC for gantry angles from -102/spl deg/ to 102/spl deg/ with a 2/spl deg/ increment. The projection image for each gantry angle was obtained using a Varian Mark II EPID. Reference images without phantom were also collected at different angles. Pixel readings of each image were converted to dose rate using the EPID calibration curve. The ray sum (sum of linear attenuation coefficients along the ray from the source to a pixel) is calculated as the negative logarithm of the ratio of dose rate of that pixel to that of the corresponding pixel in the reference image. The ray sums were then used for volumetric reconstruction using ART. ART is an iterative method that solves a system of linear equations by iteratively updating the volume to reduce the errors between the measured and calculated ray sums. The authors' results indicate that reasonably good tomographic images can be obtained using projections at every 8/spl deg/, even after only one iteration (typically, 2-3 iterations are used.) The image quality depends on the number of projections and the number of iterations. The reconstruction can be achieved within a reasonable time (/spl sim/2 hours for pure software and /spl sim/5 min with the help of graphics hardware). The authors thus conclude that cone beam CT for megavoltage therapeutic LINAC has a potential to obtain useful tomographic images.

[1]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[2]  G. N. Ramachandran,et al.  Three-dimensional reconstruction from radiographs and electron micrographs: application of convolutions instead of Fourier transforms. , 1971, Proceedings of the National Academy of Sciences of the United States of America.

[3]  H. Tuy AN INVERSION FORMULA FOR CONE-BEAM RECONSTRUCTION* , 1983 .

[4]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[5]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[6]  A. Murat Tekalp,et al.  Image Recovery: Theory and Application (Henry Stark, ed.) , 1989, SIAM Rev..

[7]  M. van Herk,et al.  Physical aspects of a liquid-filled ionization chamber with pulsed polarizing voltage. , 1991 .

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

[9]  Joab R Winkler,et al.  Numerical recipes in C: The art of scientific computing, second edition , 1993 .

[10]  H Guan,et al.  Computed tomography using algebraic reconstruction techniques (ARTs) with different projection access schemes: a comparison study under practical situations , 1996, Physics in medicine and biology.

[11]  Klaus Mueller,et al.  The weighted-distance scheme: a globally optimizing projection ordering method for ART , 1997, IEEE Transactions on Medical Imaging.

[12]  T. Mackie,et al.  Megavoltage CT on a tomotherapy system. , 1999, Physics in medicine and biology.

[13]  Klaus Mueller,et al.  Anti-Aliased 3D Cone-Beam Reconstruction of Low-Contrast Objects with Algebraic Methods , 1999, IEEE Trans. Medical Imaging.

[14]  Klaus Mueller,et al.  Fast Implementation of Algebraic Methods for 3D Reconstruction from Cone-Beam Data , 1999, IEEE Trans. Medical Imaging.

[15]  C C Ling,et al.  Relative profile and dose verification of intensity-modulated radiation therapy. , 2000, International journal of radiation oncology, biology, physics.

[16]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[17]  C C Ling,et al.  An iterative EPID calibration procedure for dosimetric verification that considers the EPID scattering factor. , 2001, Medical physics.