Fast Monte Carlo simulation for patient-specific CT/CBCT imaging dose calculation

Direct action of radiation on mummified cells: modeling of computed tomography by Monte Carlo algorithms

X-ray imaging is a nondestructive and preferred method in paleopathology to reconstruct the history of ancient diseases. Sophisticated imaging technologies such as computed tomography (CT) have become common for the investigation of skeletal disorders in human remains. Researchers have investigated the impact of ionizing radiation on living cells, but never on ancient cells in dry tissue. The effects of CT exposure on ancient cells have not been examined in the past and may be important for subsequent genetic analysis. To remedy this shortcoming, we developed different Monte Carlo models to simulate X-ray irradiation on ancient cells. Effects of mummification were considered by using two sizes of cells and three different phantom tissues, which enclosed the investigated cell cluster. This cluster was positioned at the isocenter of a CT scanner model, where the cell hit probabilities P(0,1,…, n) were calculated according to the Poisson distribution. To study the impact of the dominant physics process, CT scans for X-ray spectra of 80 and 120 kVp were simulated. Comparison between normal and dry tissue phantoms revealed that the probability of unaffected cells increased by 21 % following cell shrinkage for 80 kVp, while for 120 kVp, a further increase of unaffected cells of 23 % was observed. Consequently, cell shrinkage caused by dehydration decreased the impact of X-ray radiation on mummified cells significantly. Moreover, backscattered electrons in cortical bone protected deeper-lying ancient cells from radiation damage at 80 kVp X-rays.

Feasibility of an image planning system for kilovoltage image-guided radiation therapy

PURPOSE

Image guidance has become a standard of care for many treatment scenarios in radiation therapy. This is most typically accomplished by use of kV x-ray devices mounted onto the linear accelerator (Linac) gantry that yield planar, fluoroscopic, and cone-beam computed tomography (CBCT) images. Image acquisition parameters are chosen via preset techniques that rely on broad categorizations in patient anatomy and imaging goal. However, the optimal imaging technique results in detectability of the features of interest while exposing the patient to minimum dose. Herein, the authors present an investigation into the feasibility of developing an image planning system (IPS) for radiotherapy.

METHODS

In this first phase, the authors focused on developing an algorithm to predict tissue contrast produced by a common radiotherapy planar imaging chain. Input parameters include a CT dataset and simulated planar imaging technique settings that include kV and mAs. Energy-specific attenuation through each voxel of the CT dataset was calculated in the algorithm to derive a net transmitted intensity. The response of the flat panel detector was integrated into the image simulation algorithm. Verification was conducted by comparing simulated and measured images using four phantoms. Comparisons were made in both high and low contrast settings, as well as changes in the geometric appearance due to image saturation.

RESULTS

The authors studied a lung nodule test object to assess the planning system's ability to predict object contrast and detectability. Verification demonstrated that the slope of the pixel intensities is similar, the presence of the nodule is evident, and image saturation at high mAs values is evident in both images. The appearance of the lung nodule is a function of the image detector saturation. The authors assessed the dimensions of the lung nodule in measured and simulated images. Good quantitative agreement affirmed the algorithm's predictive capabilities. The invariance of contrast with kVp and mAs prior to saturation was predicted, as well as the gradual loss of object detectability as saturation was approached. Small changes in soft tissue density were studied using a mammography step wedge phantom. Data were acquired at beam qualities of 80 and 120 kVp and over exposure values ranging from 0.04 to 500 mAs. The data showed good agreement in terms of the absolute value of pixel intensities predicted, as well as small variations across the step wedge pattern. The saturation pixel intensity was consistent between the two beam qualities studied. Boney tissue contrast was assessed using two abdominal phantoms. Measured and calculated values agree in terms of predicting the mAs value at which detector saturation, and subsequent loss of contrast occurs. The lack of variation in contrast over mAs values lower than 10 suggests that there is wide latitude for minimizing patient dose.

CONCLUSIONS

The authors developed and tested an algorithm that can be used to assist in kV imaging technique selection during localization for radiotherapy. Phantom testing demonstrated the algorithm's predictive accuracy for both low and high contrast imaging scenarios. Detector saturation with subsequent loss of imaging detail, both in terms of object size and contrast were accurately predicted by the algorithm.

A GPU tool for efficient, accurate, and realistic simulation of cone beam CT projections

PURPOSE

Simulation of x-ray projection images plays an important role in cone beam CT (CBCT) related research projects, such as the design of reconstruction algorithms or scanners. A projection image contains primary signal, scatter signal, and noise. It is computationally demanding to perform accurate and realistic computations for all of these components. In this work, the authors develop a package on graphics processing unit (GPU), called gDRR, for the accurate and efficient computations of x-ray projection images in CBCT under clinically realistic conditions.

METHODS

The primary signal is computed by a trilinear ray-tracing algorithm. A Monte Carlo (MC) simulation is then performed, yielding the primary signal and the scatter signal, both with noise. A denoising process specifically designed for Poisson noise removal is applied to obtain a smooth scatter signal. The noise component is then obtained by combining the difference between the MC primary and the ray-tracing primary signals, and the difference between the MC simulated scatter and the denoised scatter signals. Finally, a calibration step converts the calculated noise signal into a realistic one by scaling its amplitude according to a specified mAs level. The computations of gDRR include a number of realistic features, e.g., a bowtie filter, a polyenergetic spectrum, and detector response. The implementation is fine-tuned for a GPU platform to yield high computational efficiency.

RESULTS

For a typical CBCT projection with a polyenergetic spectrum, the calculation time for the primary signal using the ray-tracing algorithms is 1.2-2.3 s, while the MC simulations take 28.1-95.3 s, depending on the voxel size. Computation time for all other steps is negligible. The ray-tracing primary signal matches well with the primary part of the MC simulation result. The MC simulated scatter signal using gDRR is in agreement with EGSnrc results with a relative difference of 3.8%. A noise calibration process is conducted to calibrate gDRR against a real CBCT scanner. The calculated projections are accurate and realistic, such that beam-hardening artifacts and scatter artifacts can be reproduced using the simulated projections. The noise amplitudes in the CBCT images reconstructed from the simulated projections also agree with those in the measured images at corresponding mAs levels.

CONCLUSIONS

A GPU computational tool, gDRR, has been developed for the accurate and efficient simulations of x-ray projections of CBCT with realistic configurations.

Four-dimensional cone beam CT reconstruction and enhancement using a temporal nonlocal means method

PURPOSE

Four-dimensional cone beam computed tomography (4D-CBCT) has been developed to provide respiratory phase-resolved volumetric imaging in image guided radiation therapy. Conventionally, it is reconstructed by first sorting the x-ray projections into multiple respiratory phase bins according to a breathing signal extracted either from the projection images or some external surrogates, and then reconstructing a 3D CBCT image in each phase bin independently using FDK algorithm. This method requires adequate number of projections for each phase, which can be achieved using a low gantry rotation or multiple gantry rotations. Inadequate number of projections in each phase bin results in low quality 4D-CBCT images with obvious streaking artifacts. 4D-CBCT images at different breathing phases share a lot of redundant information, because they represent the same anatomy captured at slightly different temporal points. Taking this redundancy along the temporal dimension into account can in principle facilitate the reconstruction in the situation of inadequate number of projection images. In this work, the authors propose two novel 4D-CBCT algorithms: an iterative reconstruction algorithm and an enhancement algorithm, utilizing a temporal nonlocal means (TNLM) method.

METHODS

The authors define a TNLM energy term for a given set of 4D-CBCT images. Minimization of this term favors those 4D-CBCT images such that any anatomical features at one spatial point at one phase can be found in a nearby spatial point at neighboring phases. 4D-CBCT reconstruction is achieved by minimizing a total energy containing a data fidelity term and the TNLM energy term. As for the image enhancement, 4D-CBCT images generated by the FDK algorithm are enhanced by minimizing the TNLM function while keeping the enhanced images close to the FDK results. A forward-backward splitting algorithm and a Gauss-Jacobi iteration method are employed to solve the problems. The algorithms implementation on GPU is designed to avoid redundant and uncoalesced memory access, in order to ensure a high computational efficiency. Our algorithms have been tested on a digital NURBS-based cardiac-torso phantom and a clinical patient case.

RESULTS

The reconstruction algorithm and the enhancement algorithm generate visually similar 4D-CBCT images, both better than the FDK results. Quantitative evaluations indicate that, compared with the FDK results, our reconstruction method improves contrast-to-noise-ratio (CNR) by a factor of 2.56-3.13 and our enhancement method increases the CNR by 2.75-3.33 times. The enhancement method also removes over 80% of the streak artifacts from the FDK results. The total computation time is 509-683 s for the reconstruction algorithm and 524-540 s for the enhancement algorithm on an NVIDIA Tesla C1060 GPU card.

CONCLUSIONS

By innovatively taking the temporal redundancy among 4D-CBCT images into consideration, the proposed algorithms can produce high quality 4D-CBCT images with much less streak artifacts than the FDK results, in the situation of inadequate number of projections.

GPU-based fast Monte Carlo dose calculation for proton therapy

Accurate radiation dose calculation is essential for successful proton radiotherapy. Monte Carlo (MC) simulation is considered to be the most accurate method. However, the long computation time limits it from routine clinical applications. Recently, graphics processing units (GPUs) have been widely used to accelerate computationally intensive tasks in radiotherapy. We have developed a fast MC dose calculation package, gPMC, for proton dose calculation on a GPU. In gPMC, proton transport is modeled by the class II condensed history simulation scheme with a continuous slowing down approximation. Ionization, elastic and inelastic proton nucleus interactions are considered. Energy straggling and multiple scattering are modeled. Secondary electrons are not transported and their energies are locally deposited. After an inelastic nuclear interaction event, a variety of products are generated using an empirical model. Among them, charged nuclear fragments are terminated with energy locally deposited. Secondary protons are stored in a stack and transported after finishing transport of the primary protons, while secondary neutral particles are neglected. gPMC is implemented on the GPU under the CUDA platform. We have validated gPMC using the TOPAS/Geant4 MC code as the gold standard. For various cases including homogeneous and inhomogeneous phantoms as well as a patient case, good agreements between gPMC and TOPAS/Geant4 are observed. The gamma passing rate for the 2%/2 mm criterion is over 98.7% in the region with dose greater than 10% maximum dose in all cases, excluding low-density air regions. With gPMC it takes only 6-22 s to simulate 10 million source protons to achieve ∼1% relative statistical uncertainty, depending on the phantoms and energy. This is an extremely high efficiency compared to the computational time of tens of CPU hours for TOPAS/Geant4. Our fast GPU-based code can thus facilitate the routine use of MC dose calculation in proton therapy.