Browsing by Author "Clark, Darin P"
Now showing 1 - 18 of 18
Results Per Page
Sort Options
Item Open Access A Deep Learning Approach for Rapid and Generalizable Denoising of Photon-Counting Micro-CT Images.(Tomography (Ann Arbor, Mich.), 2023-07) Nadkarni, Rohan; Clark, Darin P; Allphin, Alex J; Badea, Cristian TPhoton-counting CT (PCCT) is powerful for spectral imaging and material decomposition but produces noisy weighted filtered backprojection (wFBP) reconstructions. Although iterative reconstruction effectively denoises these images, it requires extensive computation time. To overcome this limitation, we propose a deep learning (DL) model, UnetU, which quickly estimates iterative reconstruction from wFBP. Utilizing a 2D U-net convolutional neural network (CNN) with a custom loss function and transformation of wFBP, UnetU promotes accurate material decomposition across various photon-counting detector (PCD) energy threshold settings. UnetU outperformed multi-energy non-local means (ME NLM) and a conventional denoising CNN called UnetwFBP in terms of root mean square error (RMSE) in test set reconstructions and their respective matrix inversion material decompositions. Qualitative results in reconstruction and material decomposition domains revealed that UnetU is the best approximation of iterative reconstruction. In reconstructions with varying undersampling factors from a high dose ex vivo scan, UnetU consistently gave higher structural similarity (SSIM) and peak signal-to-noise ratio (PSNR) to the fully sampled iterative reconstruction than ME NLM and UnetwFBP. This research demonstrates UnetU's potential as a fast (i.e., 15 times faster than iterative reconstruction) and generalizable approach for PCCT denoising, holding promise for advancing preclinical PCCT research.Item Open Access A neural network-based method for spectral distortion correction in photon counting x-ray CT.(Physics in medicine and biology, 2016-08) Touch, Mengheng; Clark, Darin P; Barber, William; Badea, Cristian TSpectral CT using a photon counting x-ray detector (PCXD) shows great potential for measuring material composition based on energy dependent x-ray attenuation. Spectral CT is especially suited for imaging with K-edge contrast agents to address the otherwise limited contrast in soft tissues. We have developed a micro-CT system based on a PCXD. This system enables both 4 energy bins acquisition, as well as full-spectrum mode in which the energy thresholds of the PCXD are swept to sample the full energy spectrum for each detector element and projection angle. Measurements provided by the PCXD, however, are distorted due to undesirable physical effects in the detector and can be very noisy due to photon starvation in narrow energy bins. To address spectral distortions, we propose and demonstrate a novel artificial neural network (ANN)-based spectral distortion correction mechanism, which learns to undo the distortion in spectral CT, resulting in improved material decomposition accuracy. To address noise, post-reconstruction denoising based on bilateral filtration, which jointly enforces intensity gradient sparsity between spectral samples, is used to further improve the robustness of ANN training and material decomposition accuracy. Our ANN-based distortion correction method is calibrated using 3D-printed phantoms and a model of our spectral CT system. To enable realistic simulations and validation of our method, we first modeled the spectral distortions using experimental data acquired from (109)Cd and (133)Ba radioactive sources measured with our PCXD. Next, we trained an ANN to learn the relationship between the distorted spectral CT projections and the ideal, distortion-free projections in a calibration step. This required knowledge of the ground truth, distortion-free spectral CT projections, which were obtained by simulating a spectral CT scan of the digital version of a 3D-printed phantom. Once the training was completed, the trained ANN was used to perform distortion correction on any subsequent scans of the same system with the same parameters. We used joint bilateral filtration to perform noise reduction by jointly enforcing intensity gradient sparsity between the reconstructed images for each energy bin. Following reconstruction and denoising, the CT data was spectrally decomposed using the photoelectric effect, Compton scattering, and a K-edge material (i.e. iodine). The ANN-based distortion correction approach was tested using both simulations and experimental data acquired in phantoms and a mouse with our PCXD-based micro-CT system for 4 bins and full-spectrum acquisition modes. The iodine detectability and decomposition accuracy were assessed using the contrast-to-noise ratio and relative error in iodine concentration estimation metrics in images with and without distortion correction. In simulation, the material decomposition accuracy in the reconstructed data was vastly improved following distortion correction and denoising, with 50% and 20% reductions in material concentration measurement error in full-spectrum and 4 energy bins cases, respectively. Overall, experimental data confirms that full-spectrum mode provides superior results to 4-energy mode when the distortion corrections are applied. The material decomposition accuracy in the reconstructed data was vastly improved following distortion correction and denoising, with as much as a 41% reduction in material concentration measurement error for full-spectrum mode, while also bringing the iodine detectability to 4-6 mg ml(-1). Distortion correction also improved the 4 bins mode data, but to a lesser extent. The results demonstrate the experimental feasibility and potential advantages of ANN-based distortion correction and joint bilateral filtration-based denoising for accurate K-edge imaging with a PCXD. Given the computational efficiency with which the ANN can be applied to projection data, the proposed scheme can be readily integrated into existing CT reconstruction pipelines.Item Open Access An Integrated Framework for Spectral and Temporal X-ray CT(2015) Clark, Darin PX-ray CT is widely used, both clinically and preclinically, for fast, high-resolution, anatomic imaging; however, compelling opportunities exist to expand its use in functional imaging applications. For instance, spectral information combined with nanoparticle contrast agents enables quantification of tissue perfusion levels, while temporal information details cardiac and respiratory dynamics. Common implementations of spectral and temporal (spectro-temporal) CT discretely sample the time points and energies to be reconstructed, proportionally increasing acquisition time and ionizing radiation dose with data dimensionality. Here, we propose and develop an integrated framework for spectro-temporal CT data acquisition, reconstruction, and analysis which drastically reduces the sampling time and radiation dose associated with spectro-temporal CT imaging. Specifically, we exploit the latent, gradient sparse and low rank structure of spectro-temporal CT data sets to recover their full dimensionality from highly undersampled projection measurements. We achieve reliable, high fidelity results through a novel combination of hierarchical projection sampling, the split Bregman optimization method, and piecewise-constant kernel regression. The integrated framework generalizes to arbitrary spectral and temporal CT reconstruction problems, while maintaining or even improving upon the sampling time and radiation dose associated with anatomic imaging protocols. We believe that this integrated framework will serve as the basis for a new generation of routine, functional CT imaging protocols.
Item Open Access Deep learning based spectral extrapolation for dual-source, dual-energy x-ray computed tomography.(Medical physics, 2020-09) Clark, Darin P; Schwartz, Fides R; Marin, Daniele; Ramirez-Giraldo, Juan C; Badea, Cristian TPurpose
Data completion is commonly employed in dual-source, dual-energy computed tomography (CT) when physical or hardware constraints limit the field of view (FoV) covered by one of two imaging chains. Practically, dual-energy data completion is accomplished by estimating missing projection data based on the imaging chain with the full FoV and then by appropriately truncating the analytical reconstruction of the data with the smaller FoV. While this approach works well in many clinical applications, there are applications which would benefit from spectral contrast estimates over the larger FoV (spectral extrapolation)-e.g. model-based iterative reconstruction, contrast-enhanced abdominal imaging of large patients, interior tomography, and combined temporal and spectral imaging.Methods
To document the fidelity of spectral extrapolation and to prototype a deep learning algorithm to perform it, we assembled a data set of 50 dual-source, dual-energy abdominal x-ray CT scans (acquired at Duke University Medical Center with 5 Siemens Flash scanners; chain A: 50 cm FoV, 100 kV; chain B: 33 cm FoV, 140 kV + Sn; helical pitch: 0.8). Data sets were reconstructed using ReconCT (v14.1, Siemens Healthineers): 768 × 768 pixels per slice, 50 cm FoV, 0.75 mm slice thickness, "Dual-Energy - WFBP" reconstruction mode with dual-source data completion. A hybrid architecture consisting of a learned piecewise linear transfer function (PLTF) and a convolutional neural network (CNN) was trained using 40 scans (five scans reserved for validation, five for testing). The PLTF learned to map chain A spectral contrast to chain B spectral contrast voxel-wise, performing an image domain analog of dual-source data completion with approximate spectral reweighting. The CNN with its U-net structure then learned to improve the accuracy of chain B contrast estimates by copying chain A structural information, by encoding prior chain A, chain B contrast relationships, and by generalizing feature-contrast associations. Training was supervised, using data from within the 33-cm chain B FoV to optimize and assess network performance.Results
Extrapolation performance on the testing data confirmed our network's robustness and ability to generalize to unseen data from different patients, yielding maximum extrapolation errors of 26 HU following the PLTF and 7.5 HU following the CNN (averaged per target organ). Degradation of network performance when applied to a geometrically simple phantom confirmed our method's reliance on feature-contrast relationships in correctly inferring spectral contrast. Integrating our image domain spectral extrapolation network into a standard dual-source, dual-energy processing pipeline for Siemens Flash scanner data yielded spectral CT data with adequate fidelity for the generation of both 50 keV monochromatic images and material decomposition images over a 30-cm FoV for chain B when only 20 cm of chain B data were available for spectral extrapolation.Conclusions
Even with a moderate amount of training data, deep learning methods are capable of robustly inferring spectral contrast from feature-contrast relationships in spectral CT data, leading to spectral extrapolation performance well beyond what may be expected at face value. Future work reconciling spectral extrapolation results with original projection data is expected to further improve results in outlying and pathological cases.Item Open Access Detection of Lung Nodules in Micro-CT Imaging Using Deep Learning(Tomography) Holbrook, Matthew D; Clark, Darin P; Patel, Rutulkumar; Qi, Yi; Bassil, Alex M; Mowery, Yvonne M; Badea, Cristian TWe are developing imaging methods for a co-clinical trial investigating synergy between immunotherapy and radiotherapy. We perform longitudinal micro-computed tomography (micro-CT) of mice to detect lung metastasis after treatment. This work explores deep learning (DL) as a fast approach for automated lung nodule detection. We used data from control mice both with and without primary lung tumors. To augment the number of training sets, we have simulated data using real augmented tumors inserted into micro-CT scans. We employed a convolutional neural network (CNN), trained with four competing types of training data: (1) simulated only, (2) real only, (3) simulated and real, and (4) pretraining on simulated followed with real data. We evaluated our model performance using precision and recall curves, as well as receiver operating curves (ROC) and their area under the curve (AUC). The AUC appears to be almost identical (0.76–0.77) for all four cases. However, the combination of real and synthetic data was shown to improve precision by 8%. Smaller tumors have lower rates of detection than larger ones, with networks trained on real data showing better performance. Our work suggests that DL is a promising approach for fast and relatively accurate detection of lung tumors in mice.Item Open Access Dual-energy micro-CT functional imaging of primary lung cancer in mice using gold and iodine nanoparticle contrast agents: a validation study.(PLoS One, 2014) Ashton, Jeffrey R; Clark, Darin P; Moding, Everett J; Ghaghada, Ketan; Kirsch, David G; West, Jennifer L; Badea, Cristian TPURPOSE: To provide additional functional information for tumor characterization, we investigated the use of dual-energy computed tomography for imaging murine lung tumors. Tumor blood volume and vascular permeability were quantified using gold and iodine nanoparticles. This approach was compared with a single contrast agent/single-energy CT method. Ex vivo validation studies were performed to demonstrate the accuracy of in vivo contrast agent quantification by CT. METHODS: Primary lung tumors were generated in LSL-Kras(G12D); p53(FL/FL) mice. Gold nanoparticles were injected, followed by iodine nanoparticles two days later. The gold accumulated in tumors, while the iodine provided intravascular contrast. Three dual-energy CT scans were performed-two for the single contrast agent method and one for the dual contrast agent method. Gold and iodine concentrations in each scan were calculated using a dual-energy decomposition. For each method, the tumor fractional blood volume was calculated based on iodine concentration, and tumor vascular permeability was estimated based on accumulated gold concentration. For validation, the CT-derived measurements were compared with histology and inductively-coupled plasma optical emission spectroscopy measurements of gold concentrations in tissues. RESULTS: Dual-energy CT enabled in vivo separation of gold and iodine contrast agents and showed uptake of gold nanoparticles in the spleen, liver, and tumors. The tumor fractional blood volume measurements determined from the two imaging methods were in agreement, and a high correlation (R(2) = 0.81) was found between measured fractional blood volume and histology-derived microvascular density. Vascular permeability measurements obtained from the two imaging methods agreed well with ex vivo measurements. CONCLUSIONS: Dual-energy CT using two types of nanoparticles is equivalent to the single nanoparticle method, but allows for measurement of fractional blood volume and permeability with a single scan. As confirmed by ex vivo methods, CT-derived nanoparticle concentrations are accurate. This method could play an important role in lung tumor characterization by CT.Item Open Access Enhancing in vivo preclinical studies with VivoVist™ and photon-counting micro-CT imaging(Medical Imaging 2024: Clinical and Biomedical Imaging, 2024-04-02) Badea, Cristian T; Rickard, Ashlyn; Allphin, Alex; Clark, Darin P; Ghaghada, Ketan B; Ridwan, S; Smilowitz, Henry M; Hainfeld, James; Mowery, Yvonne MItem Open Access Evaluating renal lesions using deep-learning based extension of dual-energy FoV in dual-source CT-A retrospective pilot study.(European journal of radiology, 2021-06) Schwartz, Fides R; Clark, Darin P; Ding, Yuqin; Ramirez-Giraldo, Juan Carlos; Badea, Cristian T; Marin, DanielePurpose
Dual-source (DS) CT, dual-energy (DE) field of view (FoV) is limited to the size of the smaller detector array. The purpose was to establish a deep learning-based approach to DE extrapolation by estimating missing image data using data from both tubes to evaluate renal lesions.Method
A DE extrapolation deep-learning (DEEDL) algorithm had been trained on DECT data of 50 patients using a DSCT with DE-FoV = 33 cm (Somatom Flash). Data from 128 patients with known renal lesions falling within DE-FoV was retrospectively collected (100/140 kVp; reference dataset 1). A smaller DE-FoV = 20 cm was simulated excluding the renal lesion of interest (dataset 2) and the DEEDL was applied to this dataset. Output from the DEEDL algorithm was evaluated using ReconCT v14.1 and Syngo.via. Mean attenuation values in lesions on mixed images (HU) were compared calculating the root-mean-squared-error (RMSE) between the datasets using MATLAB R2019a.Results
The DEEDL algorithm performed well reproducing the image data of the kidney lesions (Bosniak 1 and 2: 125, Bosniak 2F: 6, Bosniak 3: 1 and Bosniak 4/(partially) solid: 32) with RSME values of 10.59 HU, 15.7 HU for attenuation, virtual non-contrast, respectively. The measurements performed in dataset 1 and 2 showed strong correlation with linear regression (r2: attenuation = 0.89, VNC = 0.63, iodine = 0.75), lesions were classified as enhancing with an accuracy of 0.91.Conclusion
This DEEDL algorithm can be used to reconstruct a full dual-energy FoV from restricted data, enabling reliable HU value measurements in areas not covered by the smaller FoV and evaluation of renal lesions.Item Open Access Evaluation of the impact of a novel denoising algorithm on image quality in dual-energy abdominal CT of obese patients.(European radiology, 2023-04) Schwartz, Fides R; Clark, Darin P; Rigiroli, Francesca; Kalisz, Kevin; Wildman-Tobriner, Benjamin; Thomas, Sarah; Wilson, Joshua; Badea, Cristian T; Marin, DanieleObjectives
Evaluate a novel algorithm for noise reduction in obese patients using dual-source dual-energy (DE) CT imaging.Methods
Seventy-nine patients with contrast-enhanced abdominal imaging (54 women; age: 58 ± 14 years; BMI: 39 ± 5 kg/m2, range: 35-62 kg/m2) from seven DECT (SOMATOM Flash or Force) were retrospectively included (01/2019-12/2020). Image domain data were reconstructed with the standard clinical algorithm (ADMIRE/SAFIRE 2), and denoised with a comparison (ME-NLM) and a test algorithm (rank-sparse kernel regression). Contrast-to-noise ratio (CNR) was calculated. Four blinded readers evaluated the same original and denoised images (0 (worst)-100 (best)) in randomized order for perceived image noise, quality, and their comfort making a diagnosis from a table of 80 options. Comparisons between algorithms were performed using paired t-tests and mixed-effects linear modeling.Results
Average CNR was 5.0 ± 1.9 (original), 31.1 ± 10.3 (comparison; p < 0.001), and 8.9 ± 2.9 (test; p < 0.001). Readers were in good to moderate agreement over perceived image noise (ICC: 0.83), image quality (ICC: 0.71), and diagnostic comfort (ICC: 0.6). Diagnostic accuracy was low across algorithms (accuracy: 66, 63, and 67% (original, comparison, test)). The noise received a mean score of 54, 84, and 66 (p < 0.05); image quality 59, 61, and 65; and the diagnostic comfort 63, 68, and 68, respectively. Quality and comfort scores were not statistically significantly different between algorithms.Conclusions
The test algorithm produces quantitatively higher image quality than current standard and existing denoising algorithms in obese patients imaged with DECT and readers show a preference for it.Clinical relevance statement
Accurate diagnosis on CT imaging of obese patients is challenging and denoising algorithms can increase the diagnostic comfort and quantitative image quality. This could lead to better clinical reads.Key points
• Improving image quality in DECT imaging of obese patients is important for accurate and confident clinical reads, which may be aided by novel denoising algorithms using image domain data. • Accurate diagnosis on CT imaging of obese patients is especially challenging and denoising algorithms can increase quantitative and qualitative image quality. • Image domain algorithms can generalize well and can be implemented at other institutions.Item Open Access High-resolution hybrid micro-CT imaging pipeline for mouse brain region segmentation and volumetric morphometry.(PloS one, 2024-01) Nadkarni, Rohan; Han, Zay Yar; Anderson, Robert J; Allphin, Alex J; Clark, Darin P; Badea, Alexandra; Badea, Cristian TBackground
Brain region segmentation and morphometry in humanized apolipoprotein E (APOE) mouse models with a human NOS2 background (HN) contribute to Alzheimer's disease (AD) research by demonstrating how various risk factors affect the brain. Photon-counting detector (PCD) micro-CT provides faster scan times than MRI, with superior contrast and spatial resolution to energy-integrating detector (EID) micro-CT. This paper presents a pipeline for mouse brain imaging, segmentation, and morphometry from PCD micro-CT.Methods
We used brains of 26 mice from 3 genotypes (APOE22HN, APOE33HN, APOE44HN). The pipeline included PCD and EID micro-CT scanning, hybrid (PCD and EID) iterative reconstruction, and brain region segmentation using the Small Animal Multivariate Brain Analysis (SAMBA) tool. We applied SAMBA to transfer brain region labels from our new PCD CT atlas to individual PCD brains via diffeomorphic registration. Region-based and voxel-based analyses were used for comparisons by genotype and sex.Results
Together, PCD and EID scanning take ~5 hours to produce images with a voxel size of 22 μm, which is faster than MRI protocols for mouse brain morphometry with voxel size above 40 μm. Hybrid iterative reconstruction generates PCD images with minimal artifacts and higher spatial resolution and contrast than EID images. Our PCD atlas is qualitatively and quantitatively similar to the prior MRI atlas and successfully transfers labels to PCD brains in SAMBA. Male and female mice had significant volume differences in 26 regions, including parts of the entorhinal cortex and cingulate cortex. APOE22HN brains were larger than APOE44HN brains in clusters from the hippocampus, a region where atrophy is associated with AD.Conclusions
This work establishes a pipeline for mouse brain analysis using PCD CT, from staining to imaging and labeling brain images. Our results validate the effectiveness of the approach, setting a foundation for research on AD mouse models while reducing scanning durations.Item Open Access Hybrid spectral CT reconstruction.(PLoS One, 2017) Clark, Darin P; Badea, Cristian TCurrent photon counting x-ray detector (PCD) technology faces limitations associated with spectral fidelity and photon starvation. One strategy for addressing these limitations is to supplement PCD data with high-resolution, low-noise data acquired with an energy-integrating detector (EID). In this work, we propose an iterative, hybrid reconstruction technique which combines the spectral properties of PCD data with the resolution and signal-to-noise characteristics of EID data. Our hybrid reconstruction technique is based on an algebraic model of data fidelity which substitutes the EID data into the data fidelity term associated with the PCD reconstruction, resulting in a joint reconstruction problem. Within the split Bregman framework, these data fidelity constraints are minimized subject to additional constraints on spectral rank and on joint intensity-gradient sparsity measured between the reconstructions of the EID and PCD data. Following a derivation of the proposed technique, we apply it to the reconstruction of a digital phantom which contains realistic concentrations of iodine, barium, and calcium encountered in small-animal micro-CT. The results of this experiment suggest reliable separation and detection of iodine at concentrations ≥ 5 mg/ml and barium at concentrations ≥ 10 mg/ml in 2-mm features for EID and PCD data reconstructed with inherent spatial resolutions of 176 μm and 254 μm, respectively (point spread function, FWHM). Furthermore, hybrid reconstruction is demonstrated to enhance spatial resolution within material decomposition results and to improve low-contrast detectability by as much as 2.6 times relative to reconstruction with PCD data only. The parameters of the simulation experiment are based on an in vivo micro-CT experiment conducted in a mouse model of soft-tissue sarcoma. Material decomposition results produced from this in vivo data demonstrate the feasibility of distinguishing two K-edge contrast agents with a spectral separation on the order of the energy resolution of the PCD hardware.Item Open Access Material decomposition from photon-counting CT using a convolutional neural network and energy-integrating CT training labels.(Physics in medicine and biology, 2022-06-29) Nadkarni, Rohan; Allphin, Alex; Clark, Darin P; Badea, Cristian TObjective
Photon-counting CT (PCCT) has better dose efficiency and spectral resolution than energy-integrating CT, which is advantageous for material decomposition. Unfortunately, the accuracy of PCCT-based material decomposition is limited due to spectral distortions in the photon-counting detector (PCD).Approach
In this work, we demonstrate a deep learning (DL) approach that compensates for spectral distortions in the PCD and improves accuracy in material decomposition by using decomposition maps provided by high-dose multi-energy-integrating detector (EID) data as training labels. We use a 3D U-net architecture and compare networks with PCD filtered backprojection (FBP) reconstruction (FBP2Decomp), PCD iterative reconstruction (Iter2Decomp), and PCD decomposition (Decomp2Decomp) as the input.Main results
We found that our Iter2Decomp approach performs best, but DL outperforms matrix inversion decomposition regardless of the input. Compared to PCD matrix inversion decomposition, Iter2Decomp gives 27.50% lower root mean squared error (RMSE) in the iodine (I) map and 59.87% lower RMSE in the photoelectric effect (PE) map. In addition, it increases the structural similarity (SSIM) by 1.92%, 6.05%, and 9.33% in the I, Compton scattering (CS), and PE maps, respectively. When taking measurements from iodine and calcium vials, Iter2Decomp provides excellent agreement with multi-EID decomposition. One limitation is some blurring caused by our DL approach, with a decrease from 1.98 line pairs/mm at 50% modulation transfer function (MTF) with PCD matrix inversion decomposition to 1.75 line pairs/mm at 50% MTF when using Iter2Decomp.Significance
Overall, this work demonstrates that our DL approach with high-dose multi-EID derived decomposition labels is effective at generating more accurate material maps from PCD data. More accurate preclinical spectral PCCT imaging such as this could serve for developing nanoparticles that show promise in the field of theranostics (therapy and diagnostics).Item Open Access MCR toolkit: A GPU-based toolkit for multi-channel reconstruction of preclinical and clinical x-ray CT data.(Medical physics, 2023-06) Clark, Darin P; Badea, Cristian TBackground
The advancement of x-ray CT into the domains of photon counting spectral imaging and dynamic cardiac and perfusion imaging has created many new challenges and opportunities for clinicians and researchers. To address challenges such as dose constraints and scanning times while capitalizing on opportunities such as multi-contrast imaging and low-dose coronary angiography, these multi-channel imaging applications require a new generation of CT reconstruction tools. These new tools should exploit the relationships between imaging channels during reconstruction to set new image quality standards while serving as a platform for direct translation between the preclinical and clinical domains.Purpose
We outline and demonstrate a new Multi-Channel Reconstruction (MCR) Toolkit for GPU-based analytical and iterative reconstruction of preclinical and clinical multi-energy and dynamic x-ray CT data. To promote open science, open-source distribution of the Toolkit will coincide with the release of this publication (GPL v3; gitlab.oit.duke.edu/dpc18/mcr-toolkit-public).Methods
The MCR Toolkit source code is implemented in C/C++ and NVIDIA's CUDA GPU programming interface, with scripting support from MATLAB and Python. The Toolkit implements matched, separable footprint CT reconstruction operators for projection and backprojection in two geometries: planar, cone-beam CT (CBCT) and 3rd generation, cylindrical multi-detector row CT (MDCT). Analytical reconstruction is performed using filtered backprojection (FBP) for circular CBCT, weighted FBP (WFBP) for helical CBCT, and cone-parallel projection rebinning followed by WFBP for MDCT. Arbitrary combinations of energy and temporal channels are iteratively reconstructed under a generalized multi-channel signal model for joint reconstruction. We solve this generalized model algebraically using the split Bregman optimization method and the BiCGSTAB(l) linear solver interchangeably for both CBCT and MDCT data. Rank-sparse kernel regression (RSKR) and patch-based singular value thresholding (pSVT) are used to regularize the energy and time dimensions, respectively. Under a Gaussian noise model, regularization parameters are estimated automatically from the input data, dramatically reducing algorithm complexity for end users. Multi-GPU parallelization of the reconstruction operators is supported to manage reconstruction times.Results
Denoising with RSKR and pSVT and post-reconstruction material decomposition are illustrated with preclinical and clinical cardiac photon-counting (PC)CT data. A digital MOBY mouse phantom with cardiac motion is used to illustrate single energy (SE), multi-energy (ME), time resolved (TR), and combined multi-energy and time-resolved (METR) helical, CBCT reconstruction. A fixed set of projection data is used across all reconstruction cases to demonstrate the Toolkit's robustness to increasing data dimensionality. Identical reconstruction code is applied to in vivo cardiac PCCT data acquired in a mouse model of atherosclerosis (METR). Clinical cardiac CT reconstruction is illustrated using the XCAT phantom and the DukeSim CT simulator, while dual-source, dual-energy CT reconstruction is illustrated for data acquired with a Siemens Flash scanner. Benchmarking results with NVIDIA RTX 8000 GPU hardware demonstrate 61%-99% efficiency in scaling computation from one to four GPUs for these reconstruction problems.Conclusions
The MCR Toolkit provides a robust solution for temporal and spectral x-ray CT reconstruction problems and was built from the ground up to facilitate translation of CT research and development between preclinical and clinical applications.Item Open Access Photon Counting CT and Radiomic Analysis Enables Differentiation of Tumors Based on Lymphocyte Burden(Tomography) Allphin, Alex J; Mowery, Yvonne M; Lafata, Kyle J; Clark, Darin P; Bassil, Alex M; Castillo, Rico; Odhiambo, Diana; Holbrook, Matthew D; Ghaghada, Ketan B; Badea, Cristian TThe purpose of this study was to investigate if radiomic analysis based on spectral micro-CT with nanoparticle contrast-enhancement can differentiate tumors based on lymphocyte burden. High mutational load transplant soft tissue sarcomas were initiated in Rag2+/− and Rag2−/− mice to model varying lymphocyte burden. Mice received radiation therapy (20 Gy) to the tumor-bearing hind limb and were injected with a liposomal iodinated contrast agent. Five days later, animals underwent conventional micro-CT imaging using an energy integrating detector (EID) and spectral micro-CT imaging using a photon-counting detector (PCD). Tumor volumes and iodine uptakes were measured. The radiomic features (RF) were grouped into feature-spaces corresponding to EID, PCD, and spectral decomposition images. The RFs were ranked to reduce redundancy and increase relevance based on TL burden. A stratified repeated cross validation strategy was used to assess separation using a logistic regression classifier. Tumor iodine concentration was the only significantly different conventional tumor metric between Rag2+/− (TLs present) and Rag2−/− (TL-deficient) tumors. The RFs further enabled differentiation between Rag2+/− and Rag2−/− tumors. The PCD-derived RFs provided the highest accuracy (0.68) followed by decomposition-derived RFs (0.60) and the EID-derived RFs (0.58). Such non-invasive approaches could aid in tumor stratification for cancer therapy studies.Item Open Access Photon counting micro-CT imaging of Bi2WO6 nanoparticles(Medical Imaging 2024: Clinical and Biomedical Imaging, 2024-04-02) Badea, Cristian T; Bhavane, Rohan; Allphin, Alex; Nadkarni, Rohan; Clark, Darin P; Annapragada, Ananth; Ghaghada, KetanItem Open Access Photon-counting cine-cardiac CT in the mouse.(PloS one, 2019-01) Clark, Darin P; Holbrook, Matthew; Lee, Chang-Lung; Badea, Cristian TThe maturation of photon-counting detector (PCD) technology promises to enhance routine CT imaging applications with high-fidelity spectral information. In this paper, we demonstrate the power of this synergy and our complementary reconstruction techniques, performing 4D, cardiac PCD-CT data acquisition and reconstruction in a mouse model of atherosclerosis, including calcified plaque. Specifically, in vivo cardiac micro-CT scans were performed in four ApoE knockout mice, following their development of calcified plaques. The scans were performed with a prototype PCD (DECTRIS, Ltd.) with 4 energy thresholds. Projections were sampled every 10 ms with a 10 ms exposure, allowing the reconstruction of 10 cardiac phases at each of 4 energies (40 total 3D volumes per mouse scan). Reconstruction was performed iteratively using the split Bregman method with constraints on spectral rank and spatio-temporal gradient sparsity. The reconstructed images represent the first in vivo, 4D PCD-CT data in a mouse model of atherosclerosis. Robust regularization during iterative reconstruction yields high-fidelity results: an 8-fold reduction in noise standard deviation for the highest energy threshold (relative to unregularized algebraic reconstruction), while absolute spectral bias measurements remain below 13 Hounsfield units across all energy thresholds and scans. Qualitatively, image domain material decomposition results show clear separation of iodinated contrast and soft tissue from calcified plaque in the in vivo data. Quantitatively, spatial, spectral, and temporal fidelity are verified through a water phantom scan and a realistic MOBY phantom simulation experiment: spatial resolution is robustly preserved by iterative reconstruction (10% MTF: 2.8-3.0 lp/mm), left-ventricle, cardiac functional metrics can be measured from iodine map segmentations with ~1% error, and small calcifications (615 μm) can be detected during slow moving phases of the cardiac cycle. Given these preliminary results, we believe that PCD technology will enhance dynamic CT imaging applications with high-fidelity spectral and material information.Item Open Access Spectral diffusion: an algorithm for robust material decomposition of spectral CT data.(Phys Med Biol, 2014-11-07) Clark, Darin P; Badea, Cristian TClinical successes with dual energy CT, aggressive development of energy discriminating x-ray detectors, and novel, target-specific, nanoparticle contrast agents promise to establish spectral CT as a powerful functional imaging modality. Common to all of these applications is the need for a material decomposition algorithm which is robust in the presence of noise. Here, we develop such an algorithm which uses spectrally joint, piecewise constant kernel regression and the split Bregman method to iteratively solve for a material decomposition which is gradient sparse, quantitatively accurate, and minimally biased. We call this algorithm spectral diffusion because it integrates structural information from multiple spectral channels and their corresponding material decompositions within the framework of diffusion-like denoising algorithms (e.g. anisotropic diffusion, total variation, bilateral filtration). Using a 3D, digital bar phantom and a material sensitivity matrix calibrated for use with a polychromatic x-ray source, we quantify the limits of detectability (CNR = 5) afforded by spectral diffusion in the triple-energy material decomposition of iodine (3.1 mg mL(-1)), gold (0.9 mg mL(-1)), and gadolinium (2.9 mg mL(-1)) concentrations. We then apply spectral diffusion to the in vivo separation of these three materials in the mouse kidneys, liver, and spleen.Item Open Access Spectrotemporal CT data acquisition and reconstruction at low dose.(Med Phys, 2015-11) Clark, Darin P; Lee, Chang-Lung; Kirsch, David G; Badea, Cristian TPURPOSE: X-ray computed tomography (CT) is widely used, both clinically and preclinically, for fast, high-resolution anatomic imaging; however, compelling opportunities exist to expand its use in functional imaging applications. For instance, spectral information combined with nanoparticle contrast agents enables quantification of tissue perfusion levels, while temporal information details cardiac and respiratory dynamics. The authors propose and demonstrate a projection acquisition and reconstruction strategy for 5D CT (3D+dual energy+time) which recovers spectral and temporal information without substantially increasing radiation dose or sampling time relative to anatomic imaging protocols. METHODS: The authors approach the 5D reconstruction problem within the framework of low-rank and sparse matrix decomposition. Unlike previous work on rank-sparsity constrained CT reconstruction, the authors establish an explicit rank-sparse signal model to describe the spectral and temporal dimensions. The spectral dimension is represented as a well-sampled time and energy averaged image plus regularly undersampled principal components describing the spectral contrast. The temporal dimension is represented as the same time and energy averaged reconstruction plus contiguous, spatially sparse, and irregularly sampled temporal contrast images. Using a nonlinear, image domain filtration approach, the authors refer to as rank-sparse kernel regression, the authors transfer image structure from the well-sampled time and energy averaged reconstruction to the spectral and temporal contrast images. This regularization strategy strictly constrains the reconstruction problem while approximately separating the temporal and spectral dimensions. Separability results in a highly compressed representation for the 5D data in which projections are shared between the temporal and spectral reconstruction subproblems, enabling substantial undersampling. The authors solved the 5D reconstruction problem using the split Bregman method and GPU-based implementations of backprojection, reprojection, and kernel regression. Using a preclinical mouse model, the authors apply the proposed algorithm to study myocardial injury following radiation treatment of breast cancer. RESULTS: Quantitative 5D simulations are performed using the MOBY mouse phantom. Twenty data sets (ten cardiac phases, two energies) are reconstructed with 88 μm, isotropic voxels from 450 total projections acquired over a single 360° rotation. In vivo 5D myocardial injury data sets acquired in two mice injected with gold and iodine nanoparticles are also reconstructed with 20 data sets per mouse using the same acquisition parameters (dose: ∼60 mGy). For both the simulations and the in vivo data, the reconstruction quality is sufficient to perform material decomposition into gold and iodine maps to localize the extent of myocardial injury (gold accumulation) and to measure cardiac functional metrics (vascular iodine). Their 5D CT imaging protocol represents a 95% reduction in radiation dose per cardiac phase and energy and a 40-fold decrease in projection sampling time relative to their standard imaging protocol. CONCLUSIONS: Their 5D CT data acquisition and reconstruction protocol efficiently exploits the rank-sparse nature of spectral and temporal CT data to provide high-fidelity reconstruction results without increased radiation dose or sampling time.