A digitally reconstructed radiograph algorithm calculated from first principles.

PURPOSE To develop an algorithm for computing realistic digitally reconstructed radiographs (DRRs) that match real cone-beam CT (CBCT) projections with no artificial adjustments. METHODS The authors used measured attenuation data from cone-beam CT projection radiographs of different materials to obtain a function to convert CT number to linear attenuation coefficient (LAC). The effects of scatter, beam hardening, and veiling glare were first removed from the attenuation data. Using this conversion function the authors calculated the line integral of LAC through a CT along rays connecting the radiation source and detector pixels with a ray-tracing algorithm, producing raw DRRs. The effects of scatter, beam hardening, and veiling glare were then included in the DRRs through postprocessing. RESULTS The authors compared actual CBCT projections to DRRs produced with all corrections (scatter, beam hardening, and veiling glare) and to uncorrected DRRs. Algorithm accuracy was assessed through visual comparison of projections and DRRs, pixel intensity comparisons, intensity histogram comparisons, and correlation plots of DRR-to-projection pixel intensities. In general, the fully corrected algorithm provided a small but nontrivial improvement in accuracy over the uncorrected algorithm. The authors also investigated both measurement- and computation-based methods for determining the beam hardening correction, and found the computation-based method to be superior, as it accounted for nonuniform bowtie filter thickness. The authors benchmarked the algorithm for speed and found that it produced DRRs in about 0.35 s for full detector and CT resolution at a ray step-size of 0.5 mm. CONCLUSIONS The authors have demonstrated a DRR algorithm calculated from first principles that accounts for scatter, beam hardening, and veiling glare in order to produce accurate DRRs. The algorithm is computationally efficient, making it a good candidate for iterative CT reconstruction techniques that require a data fidelity term based on the matching of DRRs and projections.

[1]  P. Joseph,et al.  A Method for Correcting Bone Induced Artifacts in Computed Tomography Scanners , 1978, Journal of computer assisted tomography.

[2]  M. Levoy,et al.  Fast volume rendering using a shear-warp factorization of the viewing transformation , 1994, SIGGRAPH.

[3]  Martin J. Murphy,et al.  Reconstruction of 4D deformed CT for moving anatomy , 2008, International Journal of Computer Assisted Radiology and Surgery.

[4]  S. Leng,et al.  High temporal resolution and streak-free four-dimensional cone-beam computed tomography , 2008, Physics in medicine and biology.

[5]  M. Murphy,et al.  4D Cone-beam CT reconstruction using a motion model based on principal component analysis. , 2011, Medical physics.

[6]  Wenli Cai,et al.  DRR volume rendering using splatting in shear-warp context , 2000, 2000 IEEE Nuclear Science Symposium. Conference Record (Cat. No.00CH37149).

[7]  Zhong Su,et al.  Experimental validation of a Monte Carlo-based kV x-ray projection model for the Varian linac-mounted cone-beam CT imaging system , 2008, SPIE Medical Imaging.

[8]  S Giannouli,et al.  CT imaging based digitally reconstructed radiographs and their application in brachytherapy , 2000, Physics in medicine and biology.

[9]  Martin J Murphy,et al.  Assessing the intrinsic precision of 3D/3D rigid image registration results for patient setup in the absence of a ground truth. , 2010, Medical physics.

[10]  Chye Hwang Yan,et al.  Reconstruction algorithm for polychromatic CT imaging: application to beam hardening correction , 2000, IEEE Transactions on Medical Imaging.

[11]  O Nalcioglu,et al.  Post-reconstruction method for beam hardening in computerised tomography. , 1979, Physics in medicine and biology.

[12]  Shinichi Minohara,et al.  Development of a GPU-based multithreaded software application to calculate digitally reconstructed radiographs for radiotherapy , 2009, Radiological physics and technology.

[13]  Ann-Katherine Carton,et al.  The effect of scatter and glare on image quality in contrast-enhanced breast imaging using an a-Si/CsI(TI) full-field flat panel detector. , 2009, Medical physics.

[14]  Yuemin Zhu,et al.  Digitally reconstructed radiograph generation by an adaptive Monte Carlo method , 2006, Physics in medicine and biology.

[15]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[16]  Martin J. Murphy,et al.  Reconstruction of a cone-beam CT image via forward iterative projection matching , 2010 .

[17]  György Cserey,et al.  Fast DRR generation for 2D to 3D registration on GPUs. , 2012, Medical physics.

[18]  Ruola Ning,et al.  Cone beam CT image quality measurements: PSF de-convolution , 2007, SPIE Medical Imaging.

[19]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[20]  Kyehyun Kim,et al.  Fast 2D-3D Point-Based Registration Using GPU-Based Preprocessing for Image-Guided Surgery , 2006, CIARP.

[21]  Wolfgang Birkfellner,et al.  Fast DRR splat rendering using common consumer graphics hardware. , 2007, Medical physics.

[22]  L. Lucy.  An iterative technique for the rectification of observed distributions , 1974 .

[23]  S Webb,et al.  Removal and effects of scatter-glare in cone-beam CT with an amorphous-silicon flat-panel detector , 2011, Physics in medicine and biology.

[24]  Max A. Viergever,et al.  Mutual-information-based registration of medical images: a survey , 2003, IEEE Transactions on Medical Imaging.

[25]  G. Liney,et al.  A method to produce and validate a digitally reconstructed radiograph-based computer simulation for optimisation of chest radiographs acquired with a computed radiography imaging system. , 2011, The British journal of radiology.

[26]  Bostjan Likar,et al.  Evaluation of similarity measures for reconstruction-based registration in image-guided radiotherapy and surgery. , 2006, International journal of radiation oncology, biology, physics.

[27]  Klaus Klingenbeck,et al.  A general framework and review of scatter correction methods in x-ray cone-beam computerized tomography. Part 1: Scatter compensation approaches. , 2011, Medical physics.

[28]  Marc Levoy,et al.  Frequency domain volume rendering , 1993, SIGGRAPH.

[29]  Jeffrey F Williamson,et al.  Monte Carlo evaluation of scatter mitigation strategies in cone-beam CT. , 2010, Medical physics.

[30]  Daniel Rueckert,et al.  Fast generation of digitally reconstructed radiographs using attenuation fields with application to 2D-3D image registration , 2005, IEEE Transactions on Medical Imaging.

[31]  Bostjan Likar,et al.  Comparative evaluation of similarity measures for the rigid registration of multi-modal head images. , 2007, Physics in medicine and biology.

[32]  Wolfgang Birkfellner,et al.  Wobbled splatting—a fast perspective volume rendering method for simulation of x-ray images from CT , 2005, Physics in medicine and biology.

[33]  Jürgen Weese,et al.  A comparison of similarity measures for use in 2-D-3-D medical image registration , 1998, IEEE Transactions on Medical Imaging.

[34]  Been-Der Yang,et al.  FAST AND HIGH PRECISION DIGITALLY RECONSTRUCTED RADIOGRAPH GENERATION ON GRAPHIC PROCESSING UNIT , 2012 .