Gregory M. Sturgeon

Patient-specific radiation dose and cancer risk estimation in CT: Part I. Development and validation of a Monte Carlo program

PURPOSE Radiation-dose awareness and optimization in CT can greatly benefit from a dose-reporting system that provides dose and risk estimates specific to each patient and each CT examination. As the first step toward patient-specific dose and risk estimation, this article aimed to develop a method for accurately assessing radiation dose from CT examinations. METHODS A Monte Carlo program was developed to model a CT system (LightSpeed VCT, GE Healthcare). The geometry of the system, the energy spectra of the x-ray source, the three-dimensional geometry of the bowtie filters, and the trajectories of source motions during axial and helical scans were explicitly modeled. To validate the accuracy of the program, a cylindrical phantom was built to enable dose measurements at seven different radial distances from its central axis. Simulated radial dose distributions in the cylindrical phantom were validated against ion chamber measurements for single axial scans at all combinations of tube potential and bowtie filter settings. The accuracy of the program was further validated using two anthropomorphic phantoms (a pediatric one-year-old phantom and an adult female phantom). Computer models of the two phantoms were created based on their CT data and were voxelized for input into the Monte Carlo program. Simulated dose at various organ locations was compared against measurements made with thermoluminescent dosimetry chips for both single axial and helical scans. RESULTS For the cylindrical phantom, simulations differed from measurements by -4.8% to 2.2%. For the two anthropomorphic phantoms, the discrepancies between simulations and measurements ranged between (-8.1%, 8.1%) and (-17.2%, 13.0%) for the single axial scans and the helical scans, respectively. CONCLUSIONS The authors developed an accurate Monte Carlo program for assessing radiation dose from CT examinations. When combined with computer models of actual patients, the program can provide accurate dose estimates for specific patients.

Patient-specific radiation dose and cancer risk for pediatric chest CT.

PURPOSE To estimate patient-specific radiation dose and cancer risk for pediatric chest computed tomography (CT) and to evaluate factors affecting dose and risk, including patient size, patient age, and scanning parameters. MATERIALS AND METHODS The institutional review board approved this study and waived informed consent. This study was HIPAA compliant. The study included 30 patients (0-16 years old), for whom full-body computer models were recently created from clinical CT data. A validated Monte Carlo program was used to estimate organ dose from eight chest protocols, representing clinically relevant combinations of bow tie filter, collimation, pitch, and tube potential. Organ dose was used to calculate effective dose and risk index (an index of total cancer incidence risk). The dose and risk estimates before and after normalization by volume-weighted CT dose index (CTDI(vol)) or dose-length product (DLP) were correlated with patient size and age. The effect of each scanning parameter was studied. RESULTS Organ dose normalized by tube current-time product or CTDI(vol) decreased exponentially with increasing average chest diameter. Effective dose normalized by tube current-time product or DLP decreased exponentially with increasing chest diameter. Chest diameter was a stronger predictor of dose than weight and total scan length. Risk index normalized by tube current-time product or DLP decreased exponentially with both chest diameter and age. When normalized by DLP, effective dose and risk index were independent of collimation, pitch, and tube potential (<10% variation). CONCLUSION The correlations of dose and risk with patient size and age can be used to estimate patient-specific dose and risk. They can further guide the design and optimization of pediatric chest CT protocols. SUPPLEMENTAL MATERIAL http://radiology.rsna.org/lookup/suppl/doi:10.1148/radiol.11101900/-/DC1.

Development of realistic physical breast phantoms matched to virtual breast phantoms based on human subject data

PURPOSE Physical phantoms are essential for the development, optimization, and evaluation of x-ray breast imaging systems. Recognizing the major effect of anatomy on image quality and clinical performance, such phantoms should ideally reflect the three-dimensional structure of the human breast. Currently, there is no commercially available three-dimensional physical breast phantom that is anthropomorphic. The authors present the development of a new suite of physical breast phantoms based on human data. METHODS The phantoms were designed to match the extended cardiac-torso virtual breast phantoms that were based on dedicated breast computed tomography images of human subjects. The phantoms were fabricated by high-resolution multimaterial additive manufacturing (3D printing) technology. The glandular equivalency of the photopolymer materials was measured relative to breast tissue-equivalent plastic materials. Based on the current state-of-the-art in the technology and available materials, two variations were fabricated. The first was a dual-material phantom, the Doublet. Fibroglandular tissue and skin were represented by the most radiographically dense material available; adipose tissue was represented by the least radiographically dense material. The second variation, the Singlet, was fabricated with a single material to represent fibroglandular tissue and skin. It was subsequently filled with adipose-equivalent materials including oil, beeswax, and permanent urethane-based polymer. Simulated microcalcification clusters were further included in the phantoms via crushed eggshells. The phantoms were imaged and characterized visually and quantitatively. RESULTS The mammographic projections and tomosynthesis reconstructed images of the fabricated phantoms yielded realistic breast background. The mammograms of the phantoms demonstrated close correlation with simulated mammographic projection images of the corresponding virtual phantoms. Furthermore, power-law descriptions of the phantom images were in general agreement with real human images. The Singlet approach offered more realistic contrast as compared to the Doublet approach, but at the expense of air bubbles and air pockets that formed during the filling process. CONCLUSIONS The presented physical breast phantoms and their matching virtual breast phantoms offer realistic breast anatomy, patient variability, and ease of use, making them a potential candidate for performing both system quality control testing and virtual clinical trials.

A set of 4D pediatric XCAT reference phantoms for multimodality research.

PURPOSE The authors previously developed an adult population of 4D extended cardiac-torso (XCAT) phantoms for multimodality imaging research. In this work, the authors develop a reference set of 4D pediatric XCAT phantoms consisting of male and female anatomies at ages of newborn, 1, 5, 10, and 15 years. These models will serve as the foundation from which the authors will create a vast population of pediatric phantoms for optimizing pediatric CT imaging protocols. METHODS Each phantom was based on a unique set of CT data from a normal patient obtained from the Duke University database. The datasets were selected to best match the reference values for height and weight for the different ages and genders according to ICRP Publication 89. The major organs and structures were segmented from the CT data and used to create an initial pediatric model defined using nonuniform rational B-spline surfaces. The CT data covered the entire torso and part of the head. To complete the body, the authors manually added on the top of the head and the arms and legs using scaled versions of the XCAT adult models or additional models created from cadaver data. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from a template XCAT phantom (male or female 50th percentile adult) to the target pediatric model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. The masses of the organs in each phantom were matched to the reference values given in ICRP Publication 89. The new reference models were checked for anatomical accuracy via visual inspection. RESULTS The authors created a set of ten pediatric reference phantoms that have the same level of detail and functionality as the original XCAT phantom adults. Each consists of thousands of anatomical structures and includes parameterized models for the cardiac and respiratory motions. Based on patient data, the phantoms capture the anatomic variations of childhood, such as the development of bone in the skull, pelvis, and long bones, and the growth of the vertebrae and organs. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. CONCLUSIONS The development of patient-derived pediatric computational phantoms is useful in providing variable anatomies for simulation. Future work will expand this ten-phantom base to a host of pediatric phantoms representative of the public at large. This can provide a means to evaluate and improve pediatric imaging devices and to optimize CT protocols in terms of image quality and radiation dose.

Patient specific dosimetry phantoms using multichannel LDDMM of the whole body

This paper describes an automated procedure for creating detailed patient-specific pediatric dosimetry phantoms from a small set of segmented organs in a childs CT scan. The algorithm involves full body mappings from adult template to pediatric images using multichannel large deformation diffeomorphic metric mapping (MC-LDDMM). The parallel implementation and performance of MC-LDDMM for this application is studied here for a sample of 4 pediatric patients, and from 1 to 24 processors. 93.84% of computation time is parallelized, and the efficiency of parallelization remains high until more than 8 processors are used. The performance of the algorithm was validated on a set of 24 male and 18 female pediatric patients. It was found to be accurate typically to within 1-2 voxels (2–4 mm) and robust across this large and variable data set.

Patient Specific Computerized Phantoms to Estimate Dose in Pediatric CT

We create a series of detailed computerized phantoms to estimate patient organ and effective dose in pediatric CT and investigate techniques for efficiently creating patient-specific phantoms based on imaging data. The initial anatomy of each phantom was previously developed based on manual segmentation of pediatric CT data. Each phantom was extended to include a more detailed anatomy based on morphing an existing adult phantom in our laboratory to match the framework (based on segmentation) defined for the target pediatric model. By morphing a template anatomy to match the patient data in the LDDMM framework, it was possible to create a patient specific phantom with many anatomical structures, some not visible in the CT data. The adult models contain thousands of defined structures that were transformed to define them in each pediatric anatomy. The accuracy of this method, under different conditions, was tested using a known voxelized phantom as the target. Errors were measured in terms of a distance map between the predicted organ surfaces and the known ones. We also compared calculated dose measurements to see the effect of different magnitudes of errors in morphing. Despite some variations in organ geometry, dose measurements from morphing predictions were found to agree with those calculated from the voxelized phantom thus demonstrating the feasibility of our methods.

Population of 224 realistic human subject-based computational breast phantoms

PURPOSE To create a database of highly realistic and anatomically variable 3D virtual breast phantoms based on dedicated breast computed tomography (bCT) data. METHODS A tissue classification and segmentation algorithm was used to create realistic and detailed 3D computational breast phantoms based on 230 + dedicated bCT datasets from normal human subjects. The breast volume was identified using a coarse three-class fuzzy C-means segmentation algorithm which accounted for and removed motion blur at the breast periphery. Noise in the bCT data was reduced through application of a postreconstruction 3D bilateral filter. A 3D adipose nonuniformity (bias field) correction was then applied followed by glandular segmentation using a 3D bias-corrected fuzzy C-means algorithm. Multiple tissue classes were defined including skin, adipose, and several fractional glandular densities. Following segmentation, a skin mask was produced which preserved the interdigitated skin, adipose, and glandular boundaries of the skin interior. Finally, surface modeling was used to produce digital phantoms with methods complementary to the XCAT suite of digital human phantoms. RESULTS After rejecting some datasets due to artifacts, 224 virtual breast phantoms were created which emulate the complex breast parenchyma of actual human subjects. The volume breast density (with skin) ranged from 5.5% to 66.3% with a mean value of 25.3% ± 13.2%. Breast volumes ranged from 25.0 to 2099.6 ml with a mean value of 716.3 ± 386.5 ml. Three breast phantoms were selected for imaging with digital compression (using finite element modeling) and simple ray-tracing, and the results show promise in their potential to produce realistic simulated mammograms. CONCLUSIONS This work provides a new population of 224 breast phantoms based on in vivo bCT data for imaging research. Compared to previous studies based on only a few prototype cases, this dataset provides a rich source of new cases spanning a wide range of breast types, volumes, densities, and parenchymal patterns.

Finite‐element modeling of compression and gravity on a population of breast phantoms for multimodality imaging simulation

PURPOSE The authors are developing a series of computational breast phantoms based on breast CT data for imaging research. In this work, the authors develop a program that will allow a user to alter the phantoms to simulate the effect of gravity and compression of the breast (craniocaudal or mediolateral oblique) making the phantoms applicable to multimodality imaging. METHODS This application utilizes a template finite-element (FE) breast model that can be applied to their presegmented voxelized breast phantoms. The FE model is automatically fit to the geometry of a given breast phantom, and the material properties of each element are set based on the segmented voxels contained within the element. The loading and boundary conditions, which include gravity, are then assigned based on a user-defined position and compression. The effect of applying these loads to the breast is computed using a multistage contact analysis in FEBio, a freely available and well-validated FE software package specifically designed for biomedical applications. The resulting deformation of the breast is then applied to a boundary mesh representation of the phantom that can be used for simulating medical images. An efficient script performs the above actions seamlessly. The user only needs to specify which voxelized breast phantom to use, the compressed thickness, and orientation of the breast. RESULTS The authors utilized their FE application to simulate compressed states of the breast indicative of mammography and tomosynthesis. Gravity and compression were simulated on example phantoms and used to generate mammograms in the craniocaudal or mediolateral oblique views. The simulated mammograms show a high degree of realism illustrating the utility of the FE method in simulating imaging data of repositioned and compressed breasts. CONCLUSIONS The breast phantoms and the compression software can become a useful resource to the breast imaging research community. These phantoms can then be used to evaluate and compare imaging modalities that involve different positioning and compression of the breast.

Population of 100 Realistic, Patient-Based Computerized Breast Phantoms for Multi-modality Imaging Research

Breast imaging is an important area of research with many new techniques being investigated to further reduce the morbidity and mortality of breast cancer through early detection. Computerized phantoms can provide an essential tool to quantitatively compare new imaging systems and techniques. Current phantoms, however, lack sufficient realism in depicting the complex 3D anatomy of the breast. In this work, we created one-hundred realistic and detailed 3D computational breast phantoms based on high-resolution CT datasets from normal patients. We also developed a finiteelement application to simulate different compression states of the breast, making the phantoms applicable to multimodality imaging research. The breast phantoms and tools developed in this work were packaged into user-friendly software applications to distribute for breast imaging research.

The development of a population of 4D pediatric XCAT phantoms for imaging research and optimization

PURPOSE We previously developed a set of highly detailed 4D reference pediatric extended cardiac-torso (XCAT) phantoms at ages of newborn, 1, 5, 10, and 15 yr with organ and tissue masses matched to ICRP Publication 89 values. In this work, we extended this reference set to a series of 64 pediatric phantoms of varying age and height and body mass percentiles representative of the public at large. The models will provide a library of pediatric phantoms for optimizing pediatric imaging protocols. METHODS High resolution positron emission tomography-computed tomography data obtained from the Duke University database were reviewed by a practicing experienced radiologist for anatomic regularity. The CT portion of the data was then segmented with manual and semiautomatic methods to form a target model defined using nonuniform rational B-spline surfaces. A multichannel large deformation diffeomorphic metric mapping algorithm was used to calculate the transform from the best age matching pediatric XCAT reference phantom to the patient target. The transform was used to complete the target, filling in the nonsegmented structures and defining models for the cardiac and respiratory motions. The complete phantoms, consisting of thousands of structures, were then manually inspected for anatomical accuracy. The mass for each major tissue was calculated and compared to linearly interpolated ICRP values for different ages. RESULTS Sixty four new pediatric phantoms were created in this manner. Each model contains the same level of detail as the original XCAT reference phantoms and also includes parameterized models for the cardiac and respiratory motions. For the phantoms that were 10 yr old and younger, we included both sets of reproductive organs. This gave them the capability to simulate both male and female anatomy. With this, the population can be expanded to 92. Wide anatomical variation was clearly seen amongst the phantom models, both in organ shape and size, even for models of the same age and sex. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. CONCLUSIONS This work provides a large cohort of highly detailed pediatric phantoms with 4D capabilities of varying age, height, and body mass. The population of phantoms will provide a vital tool with which to optimize 3D and 4D pediatric imaging devices and techniques in terms of image quality and radiation-absorbed dose.