TU‐A‐213CD‐01: GDRR: A GPU Tool for Cone‐Beam CT Projection Simulations

X. Jia, H. Yan, M. Folkerts, S. Jiang

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

Purpose: Realistic calculations of x‐ray projection images play an important role in many projects related to CBCT, such as the design of scanners and reconstruction algorithms. However, to yield a desired level of realism it usually requires a tremendously long computational time, which hinders the research process. The purpose of this project is to develop a realistic x‐ray projection image simulation package, gDRR, on a computer graphics processing unit (GPU) to achieve both high accuracy and efficiency. Methods: Primary signals in a projection is computed by GPU‐based ray‐tracing algorithms, with many features considered, e.g. source energy spectrum, fluence map due to a bowtie filter, and detector response. Scatter signals are obtained by Monte Carlo (MC) simulations on GPU under the aforementioned realistic setup followed by a denoising step. Noise signals are calculated by taking the difference between the MC simulated primary and the ray‐tracing primary signals, and the difference between the MC simulated scatter signals and the denoised scatter signals. The noise level is calibrated according to the mAs level in a scan. Results: The primary signals agree well with results from MC simulations. As for the scatter signals, our MC results for real patient cases are in good agreement with those from EGSnrc with less than 2% relative error using 10 million source photons. Various realistic artifacts can be observed in CBCT images reconstructed from the simulated projections, e.g. beam hardening, scatter, and noise. The computation time per projection is ∼20ms for primary signal per energy channel, ∼4sec for scatter signal, and 1∼2sec for all other steps. Conclusions: We have developed a complete simulation package on GPU to compute x‐ray projections in CBCT. Primary, scatter, and noise signals can be calculated with a high level of realism at a high efficiency. This work is supported in part by NIH (1R01CA154747‐01), Varian Medical Systems through a Master Research Agreement, and the Thrasher Research Fund.

Original languageEnglish (US)
Pages (from-to)3890
Number of pages1
JournalMedical Physics
Volume39
Issue number6
DOIs
StatePublished - 2012

Fingerprint

Noise
X-Rays
Research
Computer Graphics
Efficiency
Photons
Artifacts

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

TU‐A‐213CD‐01 : GDRR: A GPU Tool for Cone‐Beam CT Projection Simulations. / Jia, X.; Yan, H.; Folkerts, M.; Jiang, S.

In: Medical Physics, Vol. 39, No. 6, 2012, p. 3890.

Research output: Contribution to journalArticle

@article{ec961ae1d52a46fe80fdde0c1cb04102,
title = "TU‐A‐213CD‐01: GDRR: A GPU Tool for Cone‐Beam CT Projection Simulations",
abstract = "Purpose: Realistic calculations of x‐ray projection images play an important role in many projects related to CBCT, such as the design of scanners and reconstruction algorithms. However, to yield a desired level of realism it usually requires a tremendously long computational time, which hinders the research process. The purpose of this project is to develop a realistic x‐ray projection image simulation package, gDRR, on a computer graphics processing unit (GPU) to achieve both high accuracy and efficiency. Methods: Primary signals in a projection is computed by GPU‐based ray‐tracing algorithms, with many features considered, e.g. source energy spectrum, fluence map due to a bowtie filter, and detector response. Scatter signals are obtained by Monte Carlo (MC) simulations on GPU under the aforementioned realistic setup followed by a denoising step. Noise signals are calculated by taking the difference between the MC simulated primary and the ray‐tracing primary signals, and the difference between the MC simulated scatter signals and the denoised scatter signals. The noise level is calibrated according to the mAs level in a scan. Results: The primary signals agree well with results from MC simulations. As for the scatter signals, our MC results for real patient cases are in good agreement with those from EGSnrc with less than 2{\%} relative error using 10 million source photons. Various realistic artifacts can be observed in CBCT images reconstructed from the simulated projections, e.g. beam hardening, scatter, and noise. The computation time per projection is ∼20ms for primary signal per energy channel, ∼4sec for scatter signal, and 1∼2sec for all other steps. Conclusions: We have developed a complete simulation package on GPU to compute x‐ray projections in CBCT. Primary, scatter, and noise signals can be calculated with a high level of realism at a high efficiency. This work is supported in part by NIH (1R01CA154747‐01), Varian Medical Systems through a Master Research Agreement, and the Thrasher Research Fund.",
author = "X. Jia and H. Yan and M. Folkerts and S. Jiang",
year = "2012",
doi = "10.1118/1.4735880",
language = "English (US)",
volume = "39",
pages = "3890",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "6",

}

TY - JOUR

T1 - TU‐A‐213CD‐01

T2 - GDRR: A GPU Tool for Cone‐Beam CT Projection Simulations

AU - Jia, X.

AU - Yan, H.

AU - Folkerts, M.

AU - Jiang, S.

PY - 2012

Y1 - 2012

N2 - Purpose: Realistic calculations of x‐ray projection images play an important role in many projects related to CBCT, such as the design of scanners and reconstruction algorithms. However, to yield a desired level of realism it usually requires a tremendously long computational time, which hinders the research process. The purpose of this project is to develop a realistic x‐ray projection image simulation package, gDRR, on a computer graphics processing unit (GPU) to achieve both high accuracy and efficiency. Methods: Primary signals in a projection is computed by GPU‐based ray‐tracing algorithms, with many features considered, e.g. source energy spectrum, fluence map due to a bowtie filter, and detector response. Scatter signals are obtained by Monte Carlo (MC) simulations on GPU under the aforementioned realistic setup followed by a denoising step. Noise signals are calculated by taking the difference between the MC simulated primary and the ray‐tracing primary signals, and the difference between the MC simulated scatter signals and the denoised scatter signals. The noise level is calibrated according to the mAs level in a scan. Results: The primary signals agree well with results from MC simulations. As for the scatter signals, our MC results for real patient cases are in good agreement with those from EGSnrc with less than 2% relative error using 10 million source photons. Various realistic artifacts can be observed in CBCT images reconstructed from the simulated projections, e.g. beam hardening, scatter, and noise. The computation time per projection is ∼20ms for primary signal per energy channel, ∼4sec for scatter signal, and 1∼2sec for all other steps. Conclusions: We have developed a complete simulation package on GPU to compute x‐ray projections in CBCT. Primary, scatter, and noise signals can be calculated with a high level of realism at a high efficiency. This work is supported in part by NIH (1R01CA154747‐01), Varian Medical Systems through a Master Research Agreement, and the Thrasher Research Fund.

AB - Purpose: Realistic calculations of x‐ray projection images play an important role in many projects related to CBCT, such as the design of scanners and reconstruction algorithms. However, to yield a desired level of realism it usually requires a tremendously long computational time, which hinders the research process. The purpose of this project is to develop a realistic x‐ray projection image simulation package, gDRR, on a computer graphics processing unit (GPU) to achieve both high accuracy and efficiency. Methods: Primary signals in a projection is computed by GPU‐based ray‐tracing algorithms, with many features considered, e.g. source energy spectrum, fluence map due to a bowtie filter, and detector response. Scatter signals are obtained by Monte Carlo (MC) simulations on GPU under the aforementioned realistic setup followed by a denoising step. Noise signals are calculated by taking the difference between the MC simulated primary and the ray‐tracing primary signals, and the difference between the MC simulated scatter signals and the denoised scatter signals. The noise level is calibrated according to the mAs level in a scan. Results: The primary signals agree well with results from MC simulations. As for the scatter signals, our MC results for real patient cases are in good agreement with those from EGSnrc with less than 2% relative error using 10 million source photons. Various realistic artifacts can be observed in CBCT images reconstructed from the simulated projections, e.g. beam hardening, scatter, and noise. The computation time per projection is ∼20ms for primary signal per energy channel, ∼4sec for scatter signal, and 1∼2sec for all other steps. Conclusions: We have developed a complete simulation package on GPU to compute x‐ray projections in CBCT. Primary, scatter, and noise signals can be calculated with a high level of realism at a high efficiency. This work is supported in part by NIH (1R01CA154747‐01), Varian Medical Systems through a Master Research Agreement, and the Thrasher Research Fund.

UR - http://www.scopus.com/inward/record.url?scp=84878938143&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84878938143&partnerID=8YFLogxK

U2 - 10.1118/1.4735880

DO - 10.1118/1.4735880

M3 - Article

AN - SCOPUS:84878938143

VL - 39

SP - 3890

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 6

ER -