Application of the 4-D XCAT Phantoms in Biomedical Imaging and Beyond.

The four-dimensional (4-D) eXtended CArdiac-Torso (XCAT) series of phantoms was developed to provide accurate computerized models of the human anatomy and physiology. The XCAT series encompasses a vast population of phantoms of varying ages from newborn to adult, each including parameterized models for the cardiac and respiratory motions. With great flexibility in the XCAT's design, any number of body sizes, different anatomies, cardiac or respiratory motions or patterns, patient positions and orientations, and spatial resolutions can be simulated. As such, the XCAT phantoms are gaining a wide use in biomedical imaging research. There they can provide a virtual patient base from which to quantitatively evaluate and improve imaging instrumentation, data acquisition, techniques, and image reconstruction and processing methods which can lead to improved image quality and more accurate clinical diagnoses. The phantoms have also found great use in radiation dosimetry, radiation therapy, medical device design, and even the security and defense industry. This review paper highlights some specific areas in which the XCAT phantoms have found use within biomedical imaging and other fields. From these examples, we illustrate the increasingly important role that computerized phantoms and computer simulation are playing in the research community.

A SPECT Camera for Combined MRI and SPECT for Small Animals.

We describe an MR-compatible SPECT camera for small animals. The SPECT camera system can be inserted into the bore of a state-of-the-art MRI system and allows researchers to acquire tomographic images from a mouse in-vivo with the MRI and the SPECT acquiring simultaneously. The SPECT system provides functional information, while MRI provides anatomical information. Until today it was impossible to operate conventional SPECT inside the MRI because of mutual interference. The new SPECT technology is based on semiconductor radiation sensors (CZT, ASICs), and it fits into conventional high field MRI systems with a minimum 12-cm bore size. The SPECT camera has an MR-compatible multi-pinhole collimator for mice with a ø25-mm field-of-view. For the work reported here we assembled a prototype SPECT camera system and acquired SPECT and MRI data from radioactive sources and resolution phantoms using the camera outside and inside the MRI.

MicroCT with energy-resolved photon-counting detectors.

The goal of this paper was to investigate the benefits that could be realistically achieved on a microCT imaging system with an energy-resolved photon-counting x-ray detector. To this end, we built and evaluated a prototype microCT system based on such a detector. The detector is based on cadmium telluride (CdTe) radiation sensors and application-specific integrated circuit (ASIC) readouts. Each detector pixel can simultaneously count x-ray photons above six energy thresholds, providing the capability for energy-selective x-ray imaging. We tested the spectroscopic performance of the system using polychromatic x-ray radiation and various filtering materials with K-absorption edges. Tomographic images were then acquired of a cylindrical PMMA phantom containing holes filled with various materials. Results were also compared with those acquired using an intensity-integrating x-ray detector and single-energy (i.e. non-energy-selective) CT. This paper describes the functionality and performance of the system, and presents preliminary spectroscopic and tomographic results. The spectroscopic experiments showed that the energy-resolved photon-counting detector was capable of measuring energy spectra from polychromatic sources like a standard x-ray tube, and resolving absorption edges present in the energy range used for imaging. However, the spectral quality was degraded by spectral distortions resulting from degrading factors, including finite energy resolution and charge sharing. We developed a simple charge-sharing model to reproduce these distortions. The tomographic experiments showed that the availability of multiple energy thresholds in the photon-counting detector allowed us to simultaneously measure target-to-background contrasts in different energy ranges. Compared with single-energy CT with an integrating detector, this feature was especially useful to improve differentiation of materials with different attenuation coefficient energy dependences.

4D XCAT phantom for multimodality imaging research.

PURPOSE: The authors develop the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research. METHODS: Highly detailed whole-body anatomies for the adult male and female were defined in the XCAT using nonuniform rational B-spline (NURBS) and subdivision surfaces based on segmentation of the Visible Male and Female anatomical datasets from the National Library of Medicine as well as patient datasets. Using the flexibility of these surfaces, the Visible Human anatomies were transformed to match body measurements and organ volumes for a 50th percentile (height and weight) male and female. The desired body measurements for the models were obtained using the PEOPLESIZE program that contains anthropometric dimensions categorized from 1st to the 99th percentile for US adults. The desired organ volumes were determined from ICRP Publication 89 [ICRP, "Basic anatomical and physiological data for use in radiological protection: reference values," ICRP Publication 89 (International Commission on Radiological Protection, New York, NY, 2002)]. The male and female anatomies serve as standard templates upon which anatomical variations may be modeled in the XCAT through user-defined parameters. Parametrized models for the cardiac and respiratory motions were also incorporated into the XCAT based on high-resolution cardiac- and respiratory-gated multislice CT data. To demonstrate the usefulness of the phantom, the authors show example simulation studies in PET, SPECT, and CT using publicly available simulation packages. RESULTS: As demonstrated in the pilot studies, the 4D XCAT (which includes thousands of anatomical structures) can produce realistic imaging data when combined with accurate models of the imaging process. With the flexibility of the NURBS surface primitives, any number of different anatomies, cardiac or respiratory motions or patterns, and spatial resolutions can be simulated to perform imaging research. CONCLUSIONS: With the ability to produce realistic, predictive 3D and 4D imaging data from populations of normal and abnormal patients under various imaging parameters, the authors conclude that the XCAT provides an important tool in imaging research to evaluate and improve imaging devices and techniques. In the field of x-ray CT, the phantom may also provide the necessary foundation with which to optimize clinical CT applications in terms of image quality versus radiation dose, an area of research that is becoming more significant with the growing use of CT.

Implementation of angular response function modeling in SPECT simulations with GATE.

Among Monte Carlo simulation codes in medical imaging, the GATE simulation platform is widely used today given its flexibility and accuracy, despite long run times, which in SPECT simulations are mostly spent in tracking photons through the collimators. In this work, a tabulated model of the collimator/detector response was implemented within the GATE framework to significantly reduce the simulation times in SPECT. This implementation uses the angular response function (ARF) model. The performance of the implemented ARF approach has been compared to standard SPECT GATE simulations in terms of the ARF tables' accuracy, overall SPECT system performance and run times. Considering the simulation of the Siemens Symbia T SPECT system using high-energy collimators, differences of less than 1% were measured between the ARF-based and the standard GATE-based simulations, while considering the same noise level in the projections, acceleration factors of up to 180 were obtained when simulating a planar 364 keV source seen with the same SPECT system. The ARF-based and the standard GATE simulation results also agreed very well when considering a four-head SPECT simulation of a realistic Jaszczak phantom filled with iodine-131, with a resulting acceleration factor of 100. In conclusion, the implementation of an ARF-based model of collimator/detector response for SPECT simulations within GATE significantly reduces the simulation run times without compromising accuracy.

Analytic system matrix resolution modeling in PET: an application to Rb-82 cardiac imaging.

This work explores application of a novel resolution modeling technique based on analytic physical models which individually models the various resolution degrading effects in PET (positron range, photon non-collinearity, inter-crystal scattering and inter-crystal penetration) followed by their combination and incorporation within the image reconstruction task. In addition to phantom studies, the proposed technique was particularly applied to and studied in the task of clinical Rb-82 myocardial perfusion imaging, which presently suffers from poor statistics and resolution properties in the reconstructed images. Overall, the approach is able to produce considerable enhancements in image quality. The reconstructed FWHM for a Discovery RX PET/CT scanner was seen to improve from 5.1 mm to 7.7 mm across the field-of-view (FoV) to approximately 3.5 mm nearly uniformly across the FoV. Furthermore, extended-source phantom studies indicated clearly improved images in terms of contrast versus noise performance. Using Monte Carlo simulations of clinical Rb-82 imaging, the resolution modeling technique was seen to significantly outperform standard reconstructions qualitatively, and also quantitatively in terms of contrast versus noise (contrast between the myocardium and other organs, as well as between myocardial defects and the left ventricle).

Realistic CT simulation using the 4D XCAT phantom.

The authors develop a unique CT simulation tool based on the 4D extended cardiac-torso (XCAT) phantom, a whole-body computer model of the human anatomy and physiology based on NURBS surfaces. Unlike current phantoms in CT based on simple mathematical primitives, the 4D XCAT provides an accurate representation of the complex human anatomy and has the advantage, due to its design, that its organ shapes can be changed to realistically model anatomical variations and patient motion. A disadvantage to the NURBS basis of the XCAT, however, is that the mathematical complexity of the surfaces makes the calculation of line integrals through the phantom difficult. They have to be calculated using iterative procedures; therefore, the calculation of CT projections is much slower than for simpler mathematical phantoms. To overcome this limitation, the authors used efficient ray tracing techniques from computer graphics, to develop a fast analytic projection algorithm to accurately calculate CT projections directly from the surface definition of the XCAT phantom given parameters defining the CT scanner and geometry. Using this tool, realistic high-resolution 3D and 4D projection images can be simulated and reconstructed from the XCAT within a reasonable amount of time. In comparison with other simulators with geometrically defined organs, the XCAT-based algorithm was found to be only three times slower in generating a projection data set of the same anatomical structures using a single 3.2 GHz processor. To overcome this decrease in speed would, therefore, only require running the projection algorithm in parallel over three processors. With the ever decreasing cost of computers and the rise of faster processors and multi-processor systems and clusters, this slowdown is basically inconsequential, especially given the vast improvement the XCAT offers in terms of realism and the ability to generate 3D and 4D data from anatomically diverse patients. As such, the authors conclude that the efficient XCAT-based CT simulator developed in this work will have applications in a broad range of CT imaging research.

Impact of respiratory motion on the detection of small pulmonary nodules in SPECT imaging.

The objective of this investigation is to determine the impact of respiratory motion on the detection of small solitary pulmonary nodules (SPN) in single photon emission computed tomographic (SPECT) imaging. We have previously modeled the respiratory motion of SPN based on the change of location of anatomic structures within the lungs identified on breath-held CT images of volunteers acquired at two different stages of respiration. This information on respiratory motion within the lungs was combined with the end-expiration and time-averaged NCAT phantoms to allow the creation of source and attenuation maps for the normal background distribution of Tc-99m NeoTect. With the source and attenuation distribution thus defined, the SIMIND Monte Carlo program was used to produce SPECT projection data for the normal background and separately for each of 150 end-expiration and time-averaged simulated 1.0 cm tumors. Normal and tumor SPECT projection sets each containing one lesion were combined with a clinically realistic noise level and counts. These were reconstructed with RBI-EM using 1) no correction (NC), 2) attenuation correction (AC), 3) detector response correction (RC), and 4) attenuation correction, detector response correction, and scatter correction (AC_RC_SC). The post-reconstruction parameters of number of iterations and 3-D Gaussian filtering were optimized by human-observer studies. Comparison of lesion detection by human-observer LROC studies reveals that respiratory motion degrades tumor detection for all four reconstruction strategies, and that the magnitude of this effect is greatest for NC and RC, and least for AC_RC_SC. Additionally, the AC_RC_SC strategy results in the best detection of lesions.

Fast modelling of the collimator-detector response in Monte Carlo simulation of SPECT imaging using the angular response function.

Interactions of incident photons with the collimator and detector, including septal penetration, scatter and x-ray fluorescence, are significant sources of image degradation in applications of SPECT including dual isotope imaging and imaging using radioisotopes that emit high- or medium-energy photons. Modelling these interactions using full Monte Carlo (MC) simulations is computationally very demanding. We present a new method based on the use of angular response functions (ARFs). The ARF is a function of the incident photon's direction and energy and represents the probability that a photon will either interact with or pass through the collimator, and be detected at the intersection of the photon's direction vector and the detection plane in an energy window of interest. The ARFs were pre-computed using full MC simulations of point sources that include propagation through the collimator-detector system. We have implemented the ARF method for use in conjunction with the SimSET/PHG MC code to provide fast modelling of both interactions in the patient and in the collimator-detector system. Validation results in the three cases studied show that there was good agreement between the projections generated using the ARF method and those from previously validated full MC simulations, but with hundred to thousand fold reductions in simulation time.