Modeling "Textured" Bones in Virtual Human Phantoms

Development of physiologically-informed computational coronary artery plaques for use in virtual imaging trials

BACKGROUND

As a leading cause of death, worldwide, cardiovascular disease is of great clinical importance. Among cardiovascular diseases, coronary artery disease (CAD) is a key contributor, and it is the attributed cause of death for 10% of all deaths annually. The prevalence of CAD is commensurate with the rise in new medical imaging technologies intended to aid in its diagnosis and treatment. The necessary clinical trials required to validate and optimize these technologies require a large cohort of carefully controlled patients, considerable time to complete, and can be prohibitively expensive. A safer, faster, less expensive alternative is using virtual imaging trials (VITs), utilizing virtual patients or phantoms combined with accurate computer models of imaging devices.

PURPOSE

In this work, we develop realistic, physiologically-informed models for coronary plaques for application in cardiac imaging VITs.

METHODS

Histology images of plaques at micron-level resolution were used to train a deep convolutional generative adversarial network (DC-GAN) to create a library of anatomically variable plaque models with clinical anatomical realism. The stability of each plaque was evaluated by finite element analysis (FEA) in which plaque components and vessels were meshed as volumes, modeled as specialized tissues, and subjected to the range of normal coronary blood pressures. To demonstrate the utility of the plaque models, we combined them with the whole-body XCAT computational phantom to perform initial simulations comparing standard energy-integrating detector (EID) CT with photon-counting detector (PCD) CT.

RESULTS

Our results show the network is capable of generating realistic, anatomically variable plaques. Our simulation results provide an initial demonstration of the utility of the generated plaque models as targets to compare different imaging devices.

CONCLUSIONS

Vast, realistic, and variable CAD pathologies can be generated to incorporate into computational phantoms for VITs. There they can serve as a known truth from which to optimize and evaluate cardiac imaging technologies quantitatively.

Development and Application of a Virtual Imaging Trial Framework for Longitudinal Quantification of Emphysema in CT

Pulmonary emphysema is a progressive lung disease that requires accurate evaluation for optimal management. This task, possible using quantitative CT, is particularly challenging as scanner and patient attributes change over time, negatively impacting the CT-derived quantitative measures. Efforts to minimize such variations have been limited by the absence of ground truth in clinical data, thus necessitating reliance on clinical surrogates, which may not have one-to-one correspondence to CT-based findings. This study aimed to develop the first suite of human models with emphysema at multiple time points, enabling longitudinal assessment of disease progression with access to ground truth. A total of 14 virtual subjects were modeled across three time points. Each human model was virtually imaged using a validated imaging simulator (DukeSim), modeling an energy-integrating CT scanner. The models were scanned at two dose levels and reconstructed with two reconstruction kernels, slice thicknesses, and pixel sizes. The developed longitudinal models were further utilized to demonstrate utility in algorithm testing and development. Two previously developed image processing algorithms (CT-HARMONICA, EmphysemaSeg) were evaluated. The results demonstrated the efficacy of both algorithms in improving the accuracy and precision of longitudinal quantifications, from 6.1±6.3% to 1.1±1.1% and 1.6±2.2% across years 0-5. Further investigation in EmphysemaSeg identified that baseline emphysema severity, defined as >5% emphysema at year 0, contributed to its reduced performance. This finding highlights the value of virtual imaging trials in enhancing the explainability of algorithms. Overall, the developed longitudinal human models enabled ground-truth based assessment of image processing algorithms for lung quantifications.

Technical note: Towards more realistic 4DCT(MRI) numerical lung phantoms

BACKGROUND

Numerical 4D phantoms, together with associated ground truth motion, offer a flexible and comprehensive data set for realistic simulations in radiotherapy and radiology in target sites affected by respiratory motion.

PURPOSE

We present an openly available upgrade to previously reported methods for generating realistic 4DCT lung numerical phantoms, which now incorporate respiratory ribcage motion and improved lung density representation throughout the breathing cycle.

METHODS

Density information of reference CTs, toget her with motion from multiple breathing cycle 4DMRIs have been combined to generate synthetic 4DCTs (4DCT(MRI)s). Inter-subject correspondence between the CT and MRI anatomy was first established via deformable image registration (DIR) of binary masks of the lungs and ribcage. Ribcage and lung motions were extracted independently from the 4DMRIs using DIR and applied to the corresponding locations in the CT after post-processing to preserve sliding organ motion. In addition, based on the Jacobian determinant of the resulting deformation vector fields, lung densities were scaled on a voxel-wise basis to more accurately represent changes in local lung density. For validating this process, synthetic 4DCTs, referred to as 4DCT(CT)s, were compared to the originating 4DCTs using motion extracted from the latter, and the dosimetric impact of the new features of ribcage motion and density correction were analyzed using pencil beam scanned proton 4D dose calculations.

RESULTS

Lung density scaling led to a reduction of maximum mean lung Hounsfield units (HU) differences from 45 to 12 HU when comparing simulated 4DCT(CT)s to their originating 4DCTs. Comparing 4D dose distributions calculated on the enhanced 4DCT(CT)s to those on the original 4DCTs yielded 2%/2 mm gamma pass rates above 97% with an average improvement of 1.4% compared to previously reported phantoms.

CONCLUSIONS

A previously reported 4DCT(MRI) workflow has been successfully improved and the resulting numerical phantoms exhibit more accurate lung density representations and realistic ribcage motion.

Surface-based anthropomorphic bone structures for use in high-resolution simulated medical imaging

Objective.Virtual imaging trials enable efficient assessment and optimization of medical image devices and techniques via simulation rather than physical studies. These studies require realistic, detailed ground-truth models or phantoms of the relevant anatomy or physiology. Anatomical structures within computational phantoms are typically based on medical imaging data; however, for small and intricate structures (e.g. trabecular bone), it is not reasonable to use existing clinical data as the spatial resolution of the scans is insufficient. In this study, we develop a mathematical method to generate arbitrary-resolution bone structures within virtual patient models (XCAT phantoms) to model the appearance of CT-imaged trabecular bone.Approach. Given surface definitions of a bone, an algorithm was implemented to generate stochastic bicontinuous microstructures to form a network to define the trabecular bone structure with geometric and topological properties indicative of the bone. For an example adult male XCAT phantom (50th percentile in height and weight), the method was used to generate the trabecular structure of 46 chest bones. The produced models were validated in comparison with published properties of bones. The utility of the method was demonstrated with pilot CT and photon-counting CT simulations performed using the accurate DukeSim CT simulator on the XCAT phantom containing the detailed bone models.Main results. The method successfully generated the inner trabecular structure for the different bones of the chest, having quantiative measures similar to published values. The pilot simulations showed the ability of photon-counting CT to better resolve the trabecular detail emphasizing the necessity for high-resolution bone models.Significance.As demonstrated, the developed tools have great potential to provide ground truth simulations to access the ability of existing and emerging CT imaging technology to provide quantitative information about bone structures.

Development and validation of a 3D anthropomorphic phantom for dental CBCT imaging research

BACKGROUND

Optimization of dental cone beam computed tomography (CBCT) imaging is still in a preliminary stage and should be addressed using task-based methods. Dedicated models containing relevant clinical tasks for image quality studies have yet to be developed.

PURPOSE

To present a methodology to develop and validate a virtual adult anthropomorphic voxel phantom for use in task-based image quality optimization studies in dental CBCT imaging research, focusing on root fracture (RF) detection tasks in the presence of metal artefacts.

METHODS

The phantom was developed from a CBCT scan with an isotropic voxel size of 0.2 mm, from which the main dental structures, mandible and maxilla were segmented. The missing large anatomical structures, including the spine, skull and remaining soft tissues, were segmented from a lower resolution full skull scan. Anatomical abnormalities were absent in the areas of interest. Fine detailed dental structures, that could not be segmented due to the limited resolution and noise in the clinical data, were modelled using a-priori anatomical knowledge. Model resolution of the teeth was therefore increased to 0.05 mm. Models of RFs as well as dental restorations to create the artefacts, were developed, and could be inserted in the phantom in any desired configuration. Simulated CBCT images of the models were generated using a newly developed multi-resolution simulation framework that incorporated the geometry, beam quality, noise and spatial resolution characteristics of a real dental CBCT scanner. Ray-tracing and Monte Carlo techniques were used to create the projection images, which were reconstructed using the classical FDK algorithm. Validation of the models was assessed by measurements of different tooth lengths, the pulp volume and the mandible, and comparison with reference values. Additionally, the simulated images were used in a reader study in which two oral radiologists had to score the realism level of the model's normal anatomy, as well as the modelled RFs and restorations.

RESULTS

A model of an adult head, as well as models of RFs and different types of dental restorations were created. Anatomical measurements were consistent with ranges reported in literature. For the tooth length measurements, the deviations from the mean reference values were less than 20%. In 77% of all the measurements, the deviations were within 10.1%. The pulp volumes, and mandible measurements were within one standard deviation of the reference values. Regarding the normal anatomy, both readers considered the realism level of the dental structures to be good. Background structures received a lower realism score due to the lack of detailed enough trabecular bone structure, which was expected but not the focus of this study. All modelled RFs were scored at least adequate by at least one of the readers, both in appearance and position. The realism level of the modelled restorations was considered to be good.

CONCLUSIONS

A methodology was proposed to develop and validate an anthropomorphic voxel phantom for image quality optimization studies in dental CBCT imaging, with a main focus on RF detection tasks. The methodology can be extended further to create more models representative of the clinical population.

MRXCAT2.0: Synthesis of realistic numerical phantoms by combining left-ventricular shape learning, biophysical simulations and tissue texture generation

BACKGROUND

Standardised performance assessment of image acquisition, reconstruction and processing methods is limited by the absence of images paired with ground truth reference values. To this end, we propose MRXCAT2.0 to generate synthetic data, covering healthy and pathological function, using a biophysical model. We exemplify the approach by generating cardiovascular magnetic resonance (CMR) images of healthy, infarcted, dilated and hypertrophic left-ventricular (LV) function.

METHOD

In MRXCAT2.0, the XCAT torso phantom is coupled with a statistical shape model, describing population (patho)physiological variability, and a biophysical model, providing known and detailed functional ground truth of LV morphology and function. CMR balanced steady-state free precession images are generated using MRXCAT2.0 while realistic image appearance is ensured by assigning texturized tissue properties to the phantom labels.

FINDING

Paired CMR image and ground truth data of LV function were generated with a range of LV masses (85-140 g), ejection fractions (34-51%) and peak radial and circumferential strains (0.45 to 0.95 and - 0.18 to - 0.13, respectively). These ranges cover healthy and pathological cases, including infarction, dilated and hypertrophic cardiomyopathy. The generation of the anatomy takes a few seconds and it improves on current state-of-the-art models where the pathological representation is not explicitly addressed. For the full simulation framework, the biophysical models require approximately two hours, while image generation requires a few minutes per slice.

CONCLUSION

MRXCAT2.0 offers synthesis of realistic images embedding population-based anatomical and functional variability and associated ground truth parameters to facilitate a standardized assessment of CMR acquisition, reconstruction and processing methods.

Harmonizing CT Images via Physics-based Deep Neural Networks

The rendition of medical images influences the accuracy and precision of quantifications. Image variations or biases make measuring imaging biomarkers challenging. The objective of this paper is to reduce the variability of computed tomography (CT) quantifications for radiomics and biomarkers using physics-based deep neural networks (DNNs). With the proposed framework, it is possible to harmonize the different renditions of a single CT scan (with variations in reconstruction kernel and dose) into an image that is in close agreement with the ground truth. To this end, a generative adversarial network (GAN) model was developed where the generator is informed by the scanner's modulation transfer function (MTF). To train the network, a virtual imaging trial (VIT) platform was used to acquire CT images, from a set of forty computational models (XCAT) serving as the patient model. Phantoms with varying levels of pulmonary disease, such as lung nodules and emphysema, were used. We scanned the patient models with a validated CT simulator (DukeSim) modeling a commercial CT scanner at 20 and 100 mAs dose levels and then reconstructed the images by twelve kernels representing smooth to sharp kernels. An evaluation of the harmonized virtual images was conducted in four different ways: 1) visual quality of the images, 2) bias and variation in density-based biomarkers, 3) bias and variation in morphological-based biomarkers, and 4) Noise Power Spectrum (NPS) and lung histogram. The trained model harmonized the test set images with a structural similarity index of 0.95±0.1, a normalized mean squared error of 10.2±1.5%, and a peak signal-to-noise ratio of 31.8±1.5 dB. Moreover, emphysema-based imaging biomarkers of LAA-950 (-1.5±1.8), Perc15 (13.65±9.3), and Lung mass (0.1±0.3) had more precise quantifications.

Development and Application of a Virtual Imaging Trial framework for Airway Quantifications via CT

Chronic obstructive pulmonary disease (COPD) is one of the top three causes of death worldwide, characterized by emphysema and bronchitis. Airway measurements reflect the severity of bronchitis and other airway-related diseases. Airway structures can be objectively evaluated with quantitative computed tomography (CT). The accuracy of such quantifications is limited by the spatial resolution and image noise characteristics of the imaging system and can be potentially improved with the emerging photon-counting CT (PCCT) technology. This study evaluated the quantitative performance of PCCT against energy-integrating CT (EICT) systems for airway measurements, and further identified optimum CT imaging parameters for such quantifications. The study was performed using a novel virtual imaging framework by developing the first library of virtual patients with bronchitis. These virtual patients were developed based on CT images of confirmed COPD patients with varied bronchitis severity. The human models were virtually imaged at 6.3 and 12.6 mGy dose levels using a scanner-specific simulator (DukeSim), synthesizing clinical PCCT and EICT scanners (NAEOTOM Alpha, FLASH, Siemens). The projections were reconstructed with two algorithms and kernels at different matrix sizes and slice thicknesses. The CT images were used to quantify clinically relevant airway measurements ("Pi10" and "WA%") and compared against their ground truth values. Compared to EICT, PCCT provided more accurate Pi10 and WA% measurements by 63.1% and 68.2%, respectively. For both technologies, sharper kernels and larger matrix sizes led to more reliable bronchitis quantifications. This study highlights the potential advantages of PCCT against EICT in characterizing bronchitis utilizing a virtual imaging platform.

Segars W, Samei E, M. Hertz H. Phase-contrast virtual chest radiography

Respiratory X-ray imaging enhanced by phase contrast has shown improved airway visualization in animal models. Limitations in current X-ray technology have nevertheless hindered clinical translation, leaving the potential clinical impact an open question. Here, we explore phase-contrast chest radiography in a realistic in silico framework. Specifically, we use preprocessed virtual patients to generate in silico chest radiographs by Fresnel-diffraction simulations of X-ray wave propagation. Following a reader study conducted with clinical radiologists, we predict that phase-contrast edge enhancement will have a negligible impact on improving solitary pulmonary nodule detection (6 to 20 mm). However, edge enhancement of bronchial walls visualizes small airways (< 2 mm), which are invisible in conventional radiography. Our results show that phase-contrast chest radiography could play a future role in observing small-airway obstruction (e.g., relevant for asthma or early-stage chronic obstructive pulmonary disease), which cannot be directly visualized using current clinical methods, thereby motivating the experimental development needed for clinical translation. Finally, we discuss quantitative requirements on distances and X-ray source/detector specifications for clinical implementation of phase-contrast chest radiography.

Synthetic 4DCT(MRI) lung phantom generation for 4D radiotherapy and image guidance investigations

Purpose

Respiratory motion is one of the major challenges in radiotherapy. In this work, a comprehensive and clinically plausible set of 4D numerical phantoms, together with their corresponding “ground truths,” have been developed and validated for 4D radiotherapy applications.

Methods

The phantoms are based on CTs providing density information and motion from multi‐breathing‐cycle 4D Magnetic Resonance imagings (MRIs). Deformable image registration (DIR) has been utilized to extract motion fields from 4DMRIs and to establish inter‐subject correspondence by registering binary lung masks between Computer Tomography (CT) and MRI. The established correspondence is then used to warp the CT according to the 4DMRI motion. The resulting synthetic 4DCTs are called 4DCT(MRI)s. Validation of the 4DCT(MRI) workflow was conducted by directly comparing conventional 4DCTs to derived synthetic 4D images using the motion of the 4DCTs themselves (referred to as 4DCT(CT)s). Digitally reconstructed radiographs (DRRs) as well as 4D pencil beam scanned (PBS) proton dose calculations were used for validation.

Results

Based on the CT image appearance of 13 lung cancer patients and deformable motion of five volunteer 4DMRIs, synthetic 4DCT(MRI)s with a total of 871 different breathing cycles have been generated. The 4DCT(MRI)s exhibit an average superior–inferior tumor motion amplitude of 7 ± 5 mm (min: 0.5 mm, max: 22.7 mm). The relative change of the DRR image intensities of the conventional 4DCTs and the corresponding synthetic 4DCT(CT)s inside the body is smaller than 5% for at least 81% of the pixels for all studied cases. Comparison of 4D dose distributions calculated on 4DCTs and the synthetic 4DCT(CT)s using the same motion achieved similar dose distributions with an average 2%/2 mm gamma pass rate of 90.8% (min: 77.8%, max: 97.2%).

Conclusion

We developed a series of numerical 4D lung phantoms based on real imaging and motion data, which give realistic representations of both anatomy and motion scenarios and the accessible “ground truth” deformation vector fields of each 4DCT(MRI). The open‐source code and motion data allow foreseen users to generate further 4D data by themselves. These numeric 4D phantoms can be used for the development of new 4D treatment strategies, 4D dose calculations, DIR algorithm validations, as well as simulations of motion mitigation and different online image guidance techniques for both proton and photon radiation therapy.

Anatomically- and physiologically-informed computational model of hepatic contrast perfusion for virtual imaging trials

PURPOSE

Virtual (in silico) imaging trials (VITs), involving computerized phantoms and models of the imaging process, provide a modern alternative to clinical imaging trials. VITs are faster, safer, and enable otherwise-impossible investigations. Current phantoms used in VITs are limited in their ability to model functional behavior such as contrast perfusion which is an important determinant of dose and image quality in CT imaging. In our prior work with the XCAT computational phantoms, we determined and modeled inter-organ (organ to organ) intravenous contrast concentration as a function of time from injection. However, intra-organ concentration, heterogeneous distribution within a given organ, was not pursued. We extend our methods in this work to model intra-organ concentration within the XCAT phantom with a specific focus on the liver.

METHODS

Intra-organ contrast perfusion depends on the organ's vessel network. We modeled the intricate vascular structures of the liver, informed by empirical and theoretical observations of anatomy and physiology. The developed vessel generation algorithm modeled a dual-input-single-output vascular network as a series of bifurcating surfaces to optimally deliver flow within the bounding surface of a given XCAT liver. Using this network, contrast perfusion was simulated within voxelized versions of the phantom by using knowledge of the blood velocities in each vascular structure, vessel diameters and length, and the time since the contrast entered the hepatic artery. The utility of the enhanced phantom was demonstrated through a simulation study with the phantom voxelized prior to CT simulation with the relevant liver vasculature prepared to represent blood and iodinated contrast media. The spatial extent of the blood-contrast mixture was compared to clinical data.

RESULTS

The vascular structures of the liver were generated with size and orientation which resulted in minimal energy expenditure required to maintain blood flow. Intravenous contrast was simulated as having known concentration and known total volume in the liver as calibrated from time-concentration curves (TCC). Measurements of simulated CT ROIs were found to agree with clinically-observed values of early arterial phase contrast enhancement of the parenchyma (∼5 HU). Similarly, early enhancement in the hepatic artery was found to agree with average clinical enhancement (180 HU).

CONCLUSIONS

The computational methods presented here furthered the development of the XCAT phantoms allowing for multi-timepoint contrast perfusion simulations, enabling more anthropomorphic virtual clinical trials intended for optimization of current clinical imaging technologies and applications. This article is protected by copyright. All rights reserved.

A truth-based primal-dual learning approach to reconstruct CT images utilizing the virtual imaging trial platform

Inherent to Computed tomography (CT) is image reconstruction, constructing 3D voxel values from noisy projection data. Modeling this inverse operation is not straightforward. Given the ill-posed nature of inverse problem in CT reconstruction, data-driven methods need regularization to enhance the accuracy of the reconstructed images. Besides, generalization of the results hinges upon the availability of large training datasets with access to ground truth. This paper offers a new strategy to reconstruct CT images with the advantage of ground truth accessible through a virtual imaging trial (VIT) platform. A learned primal-dual deep neural network (LPD-DNN) employed the forward model and its adjoint as a surrogate of the imaging's geometry and physics. VIT offered simulated CT projections paired with ground truth labels from anthropomorphic human models without image noise and resolution degradation. The models included a library of anthropomorphic, computational patient models (XCAT). The DukeSim simulator was utilized to form realistic projection data emulating the impact of the physics and geometry of a commercial-equivalent CT scanner. The resultant noisy sinogram data associated with each slice was thus generated for training. Corresponding linear attenuation coefficients of phantoms' materials at the effective energy of the x-ray spectrum were used as the ground truth labels. The LPD-DNN was deployed to learn the complex operators and hyper-parameters in the proximal primal-dual optimization. The obtained validation results showed a 12% normalized root mean square error with respect to the ground truth labels, a peak signal-to-noise ratio of 32 dB, a signal-to-noise ratio of 1.5, and a structural similarity index of 96%. These results were highly favorable compared to standard filtered-back projection reconstruction (65%, 17 dB, 1.0, 26%).

Photon-counting CT versus conventional CT for COPD quantifications: intra-scanner optimization and inter-scanner assessments using virtual imaging trials

Chronic obstructive pulmonary disease (COPD) is a chronic inflammatory lung disease and a major cause of death and disability worldwide. Quantitative CT is a powerful tool to better understand the heterogeneity and severity of this disease. Quantitative CT is being increasingly used in COPD research, and the recent advancements in CT technology have made it even more encouraging. One recent advancement has been the development of photon-counting detectors, offering higher spatial resolution, higher image contrast, and lower noise levels in the images. However, the quantification performance of this new technology compared to conventional scanners remains unknown. Additionally, different protocol settings (e.g., different dose levels, slice thicknesses, reconstruction kernels and algorithms) affect quantifications in an unsimilar fashion. This study investigates the potential advantages of photon-counting CT (PCCT) against conventional energy-integrating detector (EID) CT and explores the effects of protocol settings on lung density quantifications in COPD patients. This study was made possible using a virtual imaging platform, taking advantage of anthropomorphic phantoms with COPD (COPD-XCAT) and a scanner-specific CT simulator (DukeSim). Having the physical and geometrical properties of three Siemens commercial scanners (Flash, Force for EID and NAEOTOM Alpha for PCCT) modeled, we simulated CT images of ten COPD-XCAT phantoms at 0.63 and 3.17 mGy dose levels and reconstructed at three levels of kernel sharpness. The simulated CT images were quantified in terms of "Lung mean absolute error (MAE)," "LAA -950," "Perc 15," "Lung mass" imaging biomarkers and compared against the ground truth values of the phantoms. The intra-scanner assessment demonstrated the superior qualitative and quantitative performance of the PCCT scanner over the conventional scanners (21.01% and 22.74% mean lung MAE improvement, and 53.97% and 68.13% mean LAA -950 error improvement compared to Flash and Force). The results also showed that higher mAs, thinner slices, smoother kernels, and iterative reconstruction could lead to more accurate and precise quantification scores. This study underscored the qualitative and quantitative benefits of PCCT against conventional EID scanners as well as the importance of optimal protocol choice within scanners for more accurate quantifications.

A scanner-specific framework for simulating CT images with tube current modulation

Although tube current modulation (TCM) is routinely implemented in modern computed tomography (CT) scans, no existing CT simulator is capable of generating realistic images with TCM. The goal of this study was to develop such a framework to (1) facilitate patient-specific optimization of TCM parameters and (2) enable future virtual imaging trials (VITs) with more clinically realistic image quality and x-ray flux distributions. The framework was created by developing a TCM module and integrating it with an existing CT simulator (DukeSim). The developed module utilizes scanner-calibrated TCM parameters and two localizer radiographs to compute the mAs for each simulated CT projection. This simulation pipeline was validated in two parts. First, DukeSim was validated in the context of a commercial scanner with TCM (SOMATOM Force, Siemens Healthineers) by imaging a physical CT phantom (Mercury, Sun Nuclear) and its computational analogue. Second, the TCM module was validated by imaging a computational anthropomorphic phantom (ATOM, CIRS) using DukeSim with real and module-generated TCM profiles. The validation demonstrated DukeSim's realism in terms of noise magnitude, noise texture, spatial resolution, and image contrast (with average differences of 0.38%, 6.31%, 0.43%, and -9 HU, respectively). It also demonstrated the TCM module's realism in terms of projection-level mAs and resulting noise magnitude (2.86% and -2.60%, respectively). Finally, the framework was applied to a pilot VIT simulating images of three computational anthropomorphic phantoms (XCAT, with body mass indices (BMIs) of 24.3, 28.2, and 33.0) under five different TCM settings. The optimal TCM for each phantom was characterized based on various criteria, such as minimizing mAs or maximizing image quality. 'Very Weak' TCM minimized noise for the 24.3 BMI phantom, while 'Very Strong' TCM minimized noise for the 33.0 BMI phantom. This illustrates the utility of the developed framework for future optimization studies of TCM parameters and, more broadly, large-scale VITs with scanner-specific TCM.

A GPU-accelerated framework for rapid estimation of scanner-specific scatter in CT for virtual imaging trials

Virtual imaging trials (VITs), defined as the process of conducting clinical imaging trials using computer simulations, offer a time- and cost-effective alternative to traditional imaging trials for CT. The clinical potential of VITs hinges on the realism of simulations modeling the image acquisition process, where the accurate scanner-specific simulation of scatter in a time-feasible manner poses a particular challenge. To meet this need, this study proposes, develops, and validates a rapid scatter estimation framework, based on GPU-accelerated Monte Carlo (MC) simulations and denoising methods, for estimating scatter in single source, dual-source, and photon-counting CT. A CT simulator was developed to incorporate parametric models for an anti-scatter grid and a curved energy integrating detector with an energy-dependent response. The scatter estimates from the simulator were validated using physical measurements acquired on a clinical CT system using the standard single-blocker method. The MC simulator was further extended to incorporate a pre-validated model for a PCD and an additional source-detector pair to model cross scatter in dual-source configurations. To estimate scatter with desirable levels of statistical noise using a manageable computational load, two denoising methods using a (1) convolutional neural network and an (2) optimized Gaussian filter were further deployed. The viability of this framework for clinical VITs was assessed by integrating it with a scanner-specific ray-tracer program to simulate images for an image quality (Mercury) and an anthropomorphic phantom (XCAT). The simulated scatter-to-primary ratios agreed with physical measurements within 4.4% ± 10.8% across all projection angles and kVs. The differences of ∼121 HU between images with and without scatter, signifying the importance of scatter for simulating clinical images. The denoising methods preserved the magnitudes and trends observed in the reference scatter distributions, with an averaged rRMSE value of 0.91 and 0.97 for the two methods, respectively. The execution time of ∼30 s for simulating scatter in a single projection with a desirable level of statistical noise indicates a major improvement in performance, making our tool an eligible candidate for conducting extensive VITs spanning multiple patients and scan protocols.

Virtual Imaging Trials for Coronavirus Disease (COVID-19)

OBJECTIVE. The virtual imaging trial is a unique framework that can greatly facilitate the assessment and optimization of imaging methods by emulating the imaging experiment using representative computational models of patients and validated imaging simulators. The purpose of this study was to show how virtual imaging trials can be adapted for imaging studies of coronavirus disease (COVID-19), enabling effective assessment and optimization of CT and radiography acquisitions and analysis tools for reliable imaging and management of COVID-19. MATERIALS AND METHODS. We developed the first computational models of patients with COVID-19 and as a proof of principle showed how they can be combined with imaging simulators for COVID-19 imaging studies. For the body habitus of the models, we used the 4D extended cardiac-torso (XCAT) model that was developed at Duke University. The morphologic features of COVID-19 abnormalities were segmented from 20 CT images of patients who had been confirmed to have COVID-19 and incorporated into XCAT models. Within a given disease area, the texture and material of the lung parenchyma in the XCAT were modified to match the properties observed in the clinical images. To show the utility, three developed COVID-19 computational phantoms were virtually imaged using a scanner-specific CT and radiography simulator. RESULTS. Subjectively, the simulated abnormalities were realistic in terms of shape and texture. Results showed that the contrast-to-noise ratios in the abnormal regions were 1.6, 3.0, and 3.6 for 5-, 25-, and 50-mAs images, respectively. CONCLUSION. The developed toolsets in this study provide the foundation for use of virtual imaging trials in effective assessment and optimization of CT and radiography acquisitions and analysis tools to help manage the COVID-19 pandemic.

Virtual clinical trials in medical imaging: a review

The accelerating complexity and variety of medical imaging devices and methods have outpaced the ability to evaluate and optimize their design and clinical use. This is a significant and increasing challenge for both scientific investigations and clinical applications. Evaluations would ideally be done using clinical imaging trials. These experiments, however, are often not practical due to ethical limitations, expense, time requirements, or lack of ground truth. Virtual clinical trials (VCTs) (also known as in silico imaging trials or virtual imaging trials) offer an alternative means to efficiently evaluate medical imaging technologies virtually. They do so by simulating the patients, imaging systems, and interpreters. The field of VCTs has been constantly advanced over the past decades in multiple areas. We summarize the major developments and current status of the field of VCTs in medical imaging. We review the core components of a VCT: computational phantoms, simulators of different imaging modalities, and interpretation models. We also highlight some of the applications of VCTs across various imaging modalities.

Virtual clinical trial for quantifying the effects of beam collimation and pitch on image quality in computed tomography

Purpose: To utilize a virtual clinical trial (VCT) construct to investigate the effects of beam collimation and pitch on image quality (IQ) in computed tomography (CT) under different respiratory and cardiac motion rates. Approach: A computational human model [extended cardiac-torso (XCAT) phantom] with added lung lesions was used to simulate seven different rates of cardiac and respiratory motions. A validated CT simulator (DukeSim) was used in this study. A supplemental validation was done to ensure the accuracy of DukeSim across different pitches and beam collimations. Each XCAT phantom was imaged using the CT simulator at multiple pitches (0.5 to 1.5) and beam collimations (19.2 to 57.6 mm) at a constant dose level. The images were compared against the ground truth using three task-generic IQ metrics in the lungs. Additionally, the bias and variability in radiomics (morphological) feature measurements were quantified for task-specific lung lesion quantification across the studied imaging conditions. Results: All task-generic metrics degraded by 1.6% to 13.3% with increasing pitch. When imaged with motion, increasing pitch reduced motion artifacts. The IQ slightly degraded (1.3%) with changes in the studied beam collimations. Patient motion exhibited negative effects (within 7%) on the IQ. Among all features across all imaging conditions studies, compactness2 and elongation showed the largest ( - 26.5 % , 7.8%) and smallest ( - 0.8 % , 2.7%) relative bias and variability. The radiomics results were robust across the motion profiles studied. Conclusions: While high pitch and large beam collimations can negatively affect the quality of CT images, they are desirable for fast imaging. Further, our results showed no major adverse effects in morphology quantification of lung lesions with the increase in pitch or beam collimation. VCTs, such as the one demonstrated in this study, represent a viable methodology for experiments in CT.

Development of realistic multi-contrast textured XCAT (MT-XCAT) phantoms using a dual-discriminator conditional-generative adversarial network (D-CGAN)

Develop a machine learning-based method to generate multi-contrast anatomical textures in the 4D extended cardiac-torso (XCAT) phantom for more realistic imaging simulations. As a pilot study, we synthesize CT and CBCT textures in the chest region. For training purposes, major organs and gross tumor volumes (GTVs) in chest region were segmented from real patient images and assigned to different HU values to generate organ maps, which resemble the XCAT images. A dual-discriminator conditional-generative adversarial network (D-CGAN) was developed to synthesize anatomical textures in the corresponding organ maps. The D-CGAN was uniquely designed with two discriminators, one trained for the body and the other for the tumor. Various XCAT phantoms were input to the D-CGAN to generate textured XCAT phantoms. The D-CGAN model was trained separately using 62 CT and 63 CBCT images from lung SBRT patients to generate multi-contrast textured XCAT (MT-XCAT). The MT-XCAT phantoms were evaluated by comparing the intensity histograms and radiomic features with those from real patient images using Wilcoxon rank-sum test. The visual examination demonstrated that the MT-XCAT phantoms presented similar general contrast and anatomical textures as CT and CBCT images. The mean HU of the MT-XCAT-CT and MT-XCAT-CBCT were [Formula: see text] and [Formula: see text], compared with that of real CT ([Formula: see text]) and CBCT ([Formula: see text]). The majority of radiomic features from the MT-XCAT phantoms followed the same distribution as the real images according to the Wilcoxon rank-sum test, except for limited second-order features. The study demonstrated the feasibility of generating realistic MT-XCAT phantoms using D-CGAN. The MT-XCAT phantoms can be further expanded to include other modalities (MRI, PET, ultrasound, etc) under the same scheme. This crucial development greatly enhances the value of the phantom for various clinical applications, including testing and optimizing novel imaging techniques, validation of radiomics analysis methods, and virtual clinical trials.

Development of a scanner-specific simulation framework for photon-counting computed tomography

The aim of this study was to develop and validate a simulation platform that generates photon-counting CT images of voxelized phantoms with detailed modeling of manufacturer-specific components including the geometry and physics of the x-ray source, source filtrations, anti-scatter grids, and photon-counting detectors. The simulator generates projection images accounting for both primary and scattered photons using a computational phantom, scanner configuration, and imaging settings. Beam hardening artifacts are corrected using a spectrum and threshold dependent water correction algorithm. Physical and computational versions of a clinical phantom (ACR) were used for validation purposes. The physical phantom was imaged using a research prototype photon-counting CT (Siemens Healthcare) with standard (macro) mode, at four dose levels and with two energy thresholds. The computational phantom was imaged with the developed simulator with the same parameters and settings used in the actual acquisition. Images from both the real and simulated acquisitions were reconstructed using a reconstruction software (FreeCT). Primary image quality metrics such as noise magnitude, noise ratio, noise correlation coefficients, noise power spectrum, CT number, in-plane modulation transfer function, and slice sensitivity profiles were extracted from both real and simulated data and compared. The simulator was further evaluated for imaging contrast materials (bismuth, iodine, and gadolinium) at three concentration levels and six energy thresholds. Qualitatively, the simulated images showed similar appearance to the real ones. Quantitatively, the average relative error in image quality measurements were all less than 4% across all the measurements. The developed simulator will enable systematic optimization and evaluation of the emerging photon-counting computed tomography technology.

DukeSim: A Realistic, Rapid, and Scanner-Specific Simulation Framework in Computed Tomography

The purpose of this study was to develop a CT simulation platform that is: 1) compatible with voxel-based computational phantoms; 2) capable of modeling the geometry and physics of commercial CT scanners; and 3) computationally efficient. Such a simulation platform is designed to enable the virtual evaluation and optimization of CT protocols and parameters for achieving a targeted image quality while reducing radiation dose. Given a voxelized computational phantom and a parameter file describing the desired scanner and protocol, the developed platform DukeSim calculates projection images using a combination of ray-tracing and Monte Carlo techniques. DukeSim includes detailed models for the detector quantum efficiency, quantum and electronic noise, detector crosstalk, subsampling of the detector and focal spot areas, focal spot wobbling, and the bowtie filter. DukeSim was accelerated using GPU computing. The platform was validated using physical and computational versions of a phantom (Mercury phantom). Clinical and simulated CT scans of the phantom were acquired at multiple dose levels using a commercial CT scanner (Somatom Definition Flash; Siemens Healthcare). The real and simulated images were compared in terms of image contrast, noise magnitude, noise texture, and spatial resolution. The relative error between the clinical and simulated images was less than 1.4%, 0.5%, 2.6%, and 3%, for image contrast, noise magnitude, noise texture, and spatial resolution, respectively, demonstrating the high realism of DukeSim. The runtime, dependent on the imaging task and the hardware, was approximately 2-3 minutes per rotation in our study using a computer with 4 GPUs. DukeSim, when combined with realistic human phantoms, provides the necessary toolset with which to perform large-scale and realistic virtual clinical trials in a patient and scanner-specific manner.