Quantum noise properties of CT images with anatomical textured backgrounds across reconstruction algorithms: FBP and SAFIRE

Development, validation, and application of a generic image-based noise addition method for simulating reduced dose computed tomography images

BACKGROUND

Major efforts in computed tomography (CT) have been focused on reducing radiation dose to patients while maintaining adequate diagnostic quality. To that end, research tools have been developed to simulate reduced-dose images via either image-based or projection-based methods. The former is limited to fully capturing realistic texture, streak, and non-stationary characteristics of reduced dose, while the latter is impractical clinically.

PURPOSE

To develop and validate an image-based noise addition method that accounts for such attributes while being practical in clinical settings.

METHODS

A noise addition method was developed to add realistic noise in the image domain. The method first estimates the noise power spectrum (NPS) of CT images, which are also forward-projected to form synthetic projections. The projection data are supplemented with random white noise proportional to their attenuation values. The noise sinogram is then back-projected onto the image, filtered by the NPS, and scaled according to the desired dose reduction level. The tool was evaluated using both phantom images and patient data. The phantom images were acquired using a multi-sized image quality phantom (Mercury Phantom 3.0, Duke University), and a thorax anthropomorphic phantom (Lungman Phantom, Kyoto Kagaku) at different dose levels and reconstruction settings. The patient images consisted of two dose levels of various CT examinations and reconstruction settings. The simulated and real reduced-dose images were compared in terms of the noise magnitude and texture (i.e., NPS average frequency, NPS-fav). The utility of this methodology was also assessed for routine clinical use for CT protocol review.

RESULTS

For the phantom images, the percent errors in the noise magnitude between the simulated images and the actual images of the Mercury Phantom and anthropomorphic phantom images were 3.34% and 3.50%, respectively. The difference in fav was 0.07 mm-1 for the Mercury Phantom and 0.06 mm-1 for the anthropomorphic phantom between the simulated and actual images. The average noise magnitude percent error between the simulated and actual patient images was 4.61% with noise texture judged to be visually comparable with some kernel dependencies. When implemented clinically, the tool proved practical to simplify the process of estimating radiation dose reduction for CT protocols, resulting in a 50% dose reduction of our multiple myeloma protocol.

CONCLUSIONS

The method generated simulated CT images with realistic noise properties similar to images acquired at the same radiation exposure without needing access to raw projection data.

Detection of Pulmonary Nodules on Ultra-low Dose Chest Computed Tomography With Deep-learning Image Reconstruction Algorithm

PURPOSE

To evaluate the accuracy of ultra-low dose (ULD) chest computed tomography (CT), with a radiation exposure equivalent to a 2-view chest x-ray, for pulmonary nodule detection using deep learning image reconstruction (DLIR).

MATERIAL AND METHODS

This prospective cross-sectional study included 60 patients referred to our institution for assessment or follow-up of solid pulmonary nodules. All patients underwent low-dose (LD) and ULD chest CT within the same examination session. LD CT data were reconstructed using Adaptive Statistical Iterative Reconstruction-V (ASIR-V), whereas ULD CT data were reconstructed using DLIR and ASIR-V. ULD CT images were reviewed by 2 readers and LD CT images were reviewed by an experienced thoracic radiologist as the reference standard. Quantitative image quality analysis was performed, and the detectability of pulmonary nodules was assessed according to their size and location.

RESULTS

The effective radiation dose for ULD CT and LD CT were 0.13±0.01 and 1.16±0.6 mSv, respectively. Over the whole population, LD CT revealed 733 nodules. At ULD, DLIR images significantly exhibited better image quality than ASIR-V images. The overall sensitivity of DLIR reconstruction for the detection of solid pulmonary nodules from the ULD CT series was 93% and 82% for the 2 readers, with a good to excellent agreement with LD CT (ICC=0.82 and 0.66, respectively). The best sensitivities were observed in the middle lobe (97% and 85%, respectively).

CONCLUSIONS

At ULD, DLIR reconstructions, with minimal radiation exposure that could facilitate large-scale screening, allow the detection of pulmonary nodules with high sensitivity in an unrestricted BMI population.

Continuum topological derivative - a novel application tool for denoising CT and MRI medical images

BACKGROUND

CT and MRI modalities are important diagnostics tools for exploring the anatomical and tissue properties, respectively of the human beings. Several advancements like HRCT, FLAIR and Propeller have advantages in diagnosing the diseases very accurately, but still have enough space for improvements due to the presence of inherent and instrument noises. In the case of CT and MRI, the quantum mottle and the Gaussian and Rayleigh noises, respectively are still present in their advanced modalities of imaging. This paper addresses the denoising problem with continuum topological derivative technique and proved its trustworthiness based on the comparative study with other traditional filtration methods such as spatial, adaptive, frequency and transformation techniques using measures like visual inspection and performance metrics.

METHODS

This research study focuses on identifying a novel method for denoising by testing different filters on HRCT (High-Resolution Computed Tomography) and MR (Magnetic Resonance) images. The images were acquired from the Image Art Radiological Scan Centre using the SOMATOM CT and SIGNA Explorer (operating at 1.5 Tesla) machines. To compare the performance of the proposed CTD (Continuum Topological Derivative) method, various filters were tested on both HRCT and MR images. The filters tested for comparison were Gaussian (2D convolution operator), Wiener (deconvolution operator), Laplacian and Laplacian diagonal (2nd order partial differential operator), Average, Minimum, and Median (ordinary spatial operators), PMAD (Anisotropic diffusion operator), Kuan (statistical operator), Frost (exponential convolution operator), and HAAR Wavelet (time-frequency operator). The purpose of the study was to evaluate the effectiveness of the CTD method in removing noise compared to the other filters. The performance metrics were analyzed to assess the diligence of noise removal achieved by the CTD method. The primary outcome of the study was the removal of quantum mottle noise in HRCT images, while the secondary outcome focused on removing Gaussian (foreground) and Rayleigh (background) noise in MR images. The study aimed to observe the dynamics of noise removal by examining the values of the performance metrics. In summary, this study aimed to assess the denoising ability of various filters in HRCT and MR images, with the CTD method being the proposed approach. The study evaluated the performance of each filter using specific metrics and compared the results to determine the effectiveness of the CTD method in removing noise from the images.

RESULTS

Based on the calculated performance metric values, it has been observed that the CTD method successfully removed quantum mottle noise in HRCT images and Gaussian as well as Rayleigh noise in MRI. This can be evidenced by the PSNR (Peak Signal-to-Noise Ratio) metric, which consistently exhibited values ranging from 50 to 65 for all the tested images. Additionally, the CTD method demonstrated remarkably low residual values, typically on the order of e-09, which is a distinctive characteristic across all the images. Furthermore, the performance metrics of the CTD method consistently outperformed those of the other tested methods. Consequently, the results of this study have significant implications for the quality, structural similarity, and contrast of HRCT and MR images, enabling clinicians to obtain finer details for diagnostic purposes.

CONCLUSION

Continuum topological derivative algorithm is found to be constructive in removing prominent noises in both CT and MRI images and can serve as a potential tool for recognition of anatomical details in case of diseased and normal ones. The results obtained from this research work are highly inspiring and offer great promise in obtaining accurate diagnostic information for critical cases such as Thoracic Cavity Carina, Brain SPI Globe Lens 4th Ventricle, Brain-Middle Cerebral Artery, Brain-Middle Cerebral Artery and neoplastic lesions. These findings lay the foundation for implementing the proposed CTD technique in routine clinical diagnosis.

Optimization of Fast Non-Local Means Noise Reduction Algorithm Parameter in Computed Tomographic Phantom Images Using 3D Printing Technology

Noise in computed tomography (CT) is inevitably generated, which lowers the accuracy of disease diagnosis. The non-local means approach, a software technique for reducing noise, is widely used in medical imaging. In this study, we propose a noise reduction algorithm based on fast non-local means (FNLMs) and apply it to CT images of a phantom created using 3D printing technology. The self-produced phantom was manufactured using filaments with similar density to human brain tissues. To quantitatively evaluate image quality, the contrast-to-noise ratio (CNR), coefficient of variation (COV), and normalized noise power spectrum (NNPS) were calculated. The results demonstrate that the optimized smoothing factors of FNLMs are 0.08, 0.16, 0.22, 0.25, and 0.32 at 0.001, 0.005, 0.01, 0.05, and 0.1 of noise intensities, respectively. In addition, we compared the optimized FNLMs with noisy, local filters and total variation algorithms. As a result, FNLMs showed superior performance compared to various denoising techniques. Particularly, comparing the optimized FNLMs to the noisy images, the CNR improved by 6.53 to 16.34 times, COV improved by 6.55 to 18.28 times, and the NNPS improved by 10-2 mm2 on average. In conclusion, our approach shows significant potential in enhancing CT image quality with anthropomorphic phantoms, thus addressing the noise issue and improving diagnostic accuracy.

Noise power spectrum (NPS) in computed tomography: Enabling local NPS measurement without stationarity and ergodicity assumptions

BACKGROUND

Conventional methods for estimating the noise power spectrum (NPS) often necessitate multiple computed tomography (CT) data acquisitions and are required to satisfy stringent stationarity and ergodicity conditions, which prove challenging in CT imaging systems.

PURPOSE

The aim was to revisit the conventional NPS estimation method, leading to a new framework that estimates local NPS without relying on stationarity or ergodicity, thus facilitating experimental NPS estimations.

METHODS

The scientific foundation of the conventional CT NPS measurement method, based on the Wiener-Khintchine theorem, was reexamined, emphasizing the critical conditions of stationarity and ergodicity. This work proposes an alternative framework, characterized by its independence from stationarity and ergodicity, and its ability to facilitate local NPS estimations. A spatial average of local NPS over a Region of Interest (ROI) yields the conventional NPS for that ROI. The connections and differences between the proposed alternative method and the conventional method are discussed. Experimental studies were conducted to validate the new method.

RESULTS

(1) The NPS estimated using the conventional method was demonstrated to correspond to the spatial average of pointwise NPS from the proposed NPS estimation framework. (2) The NPS estimated over an ROI with the conventional method was shown to be the sum of the NPS estimated from the proposed method and a contribution from measurement uncertainty. (3) Local NPS estimations from the proposed method in this work elucidate the impact of surrounding image content on local NPS variations.

CONCLUSION

The NPS estimation method proposed in this work allows for the estimation of local NPS without relying on stationarity and ergodicity conditions, offering local NPS estimations with significantly improved precision.

Experimental measurement of local noise power spectrum (NPS) in photon counting detector-CT (PCD-CT) using a single data acquisition

BACKGROUND

Accurate noise power spectra (NPS) measurement in clinical X-ray CT exams is challenging due to the need for repeated scans, which expose patients to high radiation risks. A reliable method for single CT acquisition NPS estimation is thus highly desirable.

PURPOSE

To develop a method for estimating local NPS from a single photon counting detector-CT (PCD-CT) acquisition.

METHODS

A novel nearly statistical bias-free estimator was constructed from the raw counts data of PCD-CT scan to estimate the variance of sinogram projection data. An analytical algorithm is employed to reconstruct point-wise covariancecov ( x i , x j ) $\text{cov}({\bf x}_i,{\bf x}_j)$ between any two image pixel/voxel locationsx i ${\bf x}_i$ andx j ${\bf x_j}$ . A Fourier transform is applied to obtain the desired point-wise NPS for any chosen locationx i ${\bf x}_i$ . The method was validated using experimental data acquired from a benchtop PCD-CT system with various physical phantoms, and the results were compared with the conventional local NPS measurement method using repeated scans and statistical ensemble averaging.

RESULTS

The experimental results demonstrate that (1) the proposed method can achieve pointwise/local NPS measurement for a region of interest (ROI) located at any chosen position, accurately characterizing the NPS with spatial structures resulting from image content heterogeneity; (2) the local NPS measured using the proposed method show a higher precision in the measured NPS compared to the conventional measurement method; (3) spatial averaging of the local NPS yields the conventional NPS for a given local ROI.

CONCLUSION

A new method was developed to enable local NPS from a single PCD-CT acquisition.

3D printed phantom with 12 000 submillimeter lesions to improve efficiency in CT detectability assessment

BACKGROUND

The detectability performance of a CT scanner is difficult to precisely quantify when nonlinearities are present in reconstruction. An efficient detectability assessment method that is sensitive to small effects of dose and scanner settings is desirable. We previously proposed a method using a search challenge instrument: a phantom is embedded with hundreds of lesions at random locations, and a model observer is used to detect lesions. Preliminary tests in simulation and a prototype showed promising results.

PURPOSE

In this work, we fabricated a full-size search challenge phantom with design updates, including changes to lesion size, contrast, and number, and studied our implementation by comparing the lesion detectability from a nonprewhitening (NPW) model observer between different reconstructions at different exposure levels, and by estimating the instrument sensitivity to detect changes in dose.

METHODS

Designed to fit into QRM anthropomorphic phantoms, our search challenge phantom is a cylindrical insert 10 cm wide and 4 cm thick, embedded with 12 000 lesions (nominal width of 0.6 mm, height of 0.8 mm, and contrast of -350 HU), and was fabricated using PixelPrint, a 3D printing technique. The insert was scanned alone at a high dose to assess printing accuracy. To evaluate lesion detectability, the insert was placed in a QRM thorax phantom and scanned from 50 to 625 mAs with increments of 25 mAs, once per exposure level, and the average of all exposure levels was used as high-dose reference. Scans were reconstructed with three different settings: filtered-backprojection (FBP) with Br40 and Br59, and Sinogram Affirmed Iterative Reconstruction (SAFIRE) with strength level 5 and Br59 kernel. An NPW model observer was used to search for lesions, and detection performance of different settings were compared using area under the exponential transform of free response ROC curve (AUC). Using propagation of uncertainty, the sensitivity to changes in dose was estimated by the percent change in exposure due to one standard deviation of AUC, measured from 5 repeat scans at 100, 200, 300, and 400 mAs.

RESULTS

The printed insert lesions had an average position error of 0.20 mm compared to printing reference. As the exposure level increases from 50 mAs to 625 mAs, the lesion detectability AUCs increase from 0.38 to 0.92, 0.42 to 0.98, and 0.41 to 0.97 for FBP Br40, FBP Br59, and SAFIRE Br59, respectively, with a lower rate of increase at higher exposure level. FBP Br59 performed best with AUC 0.01 higher than SAFIRE Br59 on average and 0.07 higher than FBP Br40 (all P < 0.001). The standard deviation of AUC was less than 0.006, and the sensitivity to detect changes in mAs was within 2% for FBP Br59.

CONCLUSIONS

Our 3D-printed search challenge phantom with 12 000 submillimeter lesions, together with an NPW model observer, provide an efficient CT detectability assessment method that is sensitive to subtle effects in reconstruction and is sensitive to small changes in dose.

Tailoring the Mass Density of 3D Printing Materials for Accurate X-ray Imaging Simulation by Controlled Underfilling for Radiographic Phantoms

Additive manufacturing and 3D printing allow for the design and rapid production of radiographic phantoms for X-ray imaging, including CT. These are used for numerous purposes, such as patient simulation, optimization of imaging procedures and dose levels, system evaluation and quality assurance. However, standard 3D printing polymers do not mimic X-ray attenuation properties of tissues like soft, adipose, lung or bone tissue, and standard materials like liquid water. The mass density of printing polymers-especially important in CT-is often inappropriate, i.e., mostly too high. Different methods can be applied to reduce mass density. This work examines reducing density by controlled underfilling either realized by using 3D printing materials expanded through foaming during heating in the printing process, or reducing polymer flow to introduce microscopic air-filled voids. The achievable density reduction depends on the base polymer used. When using foaming materials, density is controlled by the extrusion temperature, and ranges from 33 to 47% of the base polymer used, corresponding to a range of -650 to -394 HU in CT with 120 kV. Standard filaments (Nylon, modified PLA and modified ABS) allowed density reductions by 20 to 25%, covering HU values in CT from -260 to 77 (Nylon), -230 to -20 (ABS) and -81 to 143 (PLA). A standard chalk-filled PLA filament allowed reproduction of bone tissue in a wide range of bone mineral content resulting in CT numbers from 57 to 460 HU. Controlled underfilling allowed the production of radiographic phantom materials with continuously adjustable attenuation in a limited but appropriate range, allowing for the reproduction of X-ray attenuation properties of water, adipose, soft, lung, and bone tissue in an accurate, predictable and reproducible manner.

COMBINING HI-RESOLUTION SCAN MODE WITH DEEP LEARNING RECONSTRUCTION ALGORITHMS IN CARDIAC CT

To investigate the impact of combining the high-resolution (Hi-res) scan mode with deep learning image reconstruction (DLIR) algorithm in CT. Two phantoms (Catphan600® and Lungman, small, medium, large size) were CT scanned using combinations of Hi-res/standard mode and high-definition (HD)/standard kernels. Images were reconstructed with ASiR-V and three levels of DLIR. Spatial resolution, noise and contrast-to-noise ratio (CNR) were assessed. The radiation dose was recorded. The spatial resolution increased using Hi-res & HD. Image noise in the Catphan600® (69%) and the Lungman (10-70%) significantly increased when Hi-res & HD was applied. DLIR reduced the mean noise (54%). The CNR was reduced (64%) for Hi-res & HD. The radiation dose increased for both small (+70%) and medium (+43%) Lungman phantoms but decreased slightly for the large ones (-3%) when Hi-res was applied. In conclusion, the Hi-res scan mode improved the spatial resolution. The HD kernel significantly increased the image noise. DLIR improved the image noise and CNR and did not affect the spatial resolution.

Pulmonary nodule volumetric accuracy of a deep learning-based reconstruction algorithm in low-dose computed tomography: A phantom study

PURPOSE

To compare the image properties and pulmonary nodule volumetric accuracies among deep learning-based reconstruction (DLR), filtered back projection (FBP), and hybrid iterative reconstruction (hybrid IR) in low-dose computed tomography (LDCT).

METHODS

A multipurpose chest phantom containing artificial spherical pulmonary nodules with 5-, 8-, 10-, and 12-mm diameters and Hounsfield units (HUs) of -630 and +100 HU was scanned 20 times at a standard dose, based on a low-dose screening CT trial, and at 1/2, 1/6, and 1/12 of the standard dose. To assess noise reduction performance and volumetric accuracy, the standard deviations (SDs) of the pixel values and volumetric percentage errors (PEs) were compared among FBP, hybrid IR, and DLR. The noise non-stationarity index (NNSI) was calculated from 20 image replicates and compared among FBP, hybrid IR, and DLR to evaluate noise stationarity.

RESULTS

The SD reduction rates for FBP in hybrid IR and DLR were 62 %-85 % and 79 %-90 %, respectively. For the four nodules with +100 HU, the PE of all reconstruction methods was <±25 % (not clinically relevant). For the four nodules with -630 HU, the PEs were equivalent or lower for hybrid IR and DLR than for FBP, and the PE difference between hybrid IR and DLR ranged from 0 % to 7%. The NNSI was significantly higher for DLR than for FBP and hybrid IR (p < 0.01).

CONCLUSIONS

Greater noise suppression was achieved with DLR than with hybrid IR without compromising nodule volumetric accuracy in LDCT despite the representative noise non-stationarity.

Comparative Assessment of Noise Properties for Two Deep Learning CT Image Reconstruction Techniques and Filtered Back Projection

BACKGROUND

Two deep-learning image reconstruction (DLIR) techniques from two different CT vendors have recently been introduced into clinical practice.

PURPOSE

To characterize the noise properties of two DLIR techniques with different training methods, using a phantom containing a simple uniform and a complex non-uniform region.

METHODS

A water-bath phantom with a diameter of 300 mm was used as a base phantom. A textured phantom with a diameter of 128 mm, which was made of two materials, one equivalent to water and the other being 12 mg/mL diluted iodine, irregularly mixed to create a complex texture (non-uniform region), was placed in the base phantom. Thirty repeated phantom scans were performed using two CT scanners (GE, Revolution CT with Apex Edition; Canon, Aquilion One PRISM Edition) at two dose levels (CTDI: 5 and 15 mGy). Images were reconstructed with each CT system's filtered back projection (FBP) and DLIR [GE, TrueFidelity (TF); Canon, Advanced intelligent Clear-IQ Engine Body Sharp (AC)] for three process strengths. For basic characteristics of noise, the standard deviation (SD) and noise power spectrum (NPS) were measured for the uniform (water) region. A noise magnitude map was generated by calculating the inter-image SD at each pixel position across the 30 images. Then, a noise reduction map (NRM), which visualizes the relative differences in noise magnitude between FBP and DLIR, was calculated. The NRM values ranged from 0.0 to 1.0. A low NRM value represents a less aggressive noise reduction. The histograms of the NRM value were analyzed for the uniform and non-uniform regions.

RESULTS

The reduction in noise magnitude compared with FBP tended to be greater with AC (45%-85%) than with TF (32%-65%). The average NPS frequencies of TF and AC were almost comparable to those of FBP, except for the low-dose condition and the high noise reduction strength for AC. The NRM values of TF and AC were higher in the uniform region than in the non-uniform region. In the non-uniform region, TF's average NRM values (0.21-0.48) tended to be lower than AC's (0.39-0.78). The histograms for TF showed a small overlap between the uniform and the non-uniform regions; in contrast, those for AC showed a greater overlap. This difference seems to indicate that TF processes the uniform and non-uniform regions more differently than AC does.

CONCLUSION

This study has revealed a distinct difference in characteristics between the two DLIR techniques: TF tends to offer less aggressive noise reduction in non-uniform regions and preserve the original signals, whereas AC tends to prioritize noise filtering over edge-preservation, especially at the low-dose condition and with the high noise reduction strength. This article is protected by copyright. All rights reserved.

Comparison of the Classification Results Accuracy for CT Soft Tissue and Bone Reconstructions in Detecting the Porosity of a Spongy Tissue

The aim of the study was to compare the accuracy of the classification pertaining to the results of two types of soft tissue and bone reconstructions of the spinal CT in detecting the porosity of L1 vertebral body spongy tissue. The dataset for each type of reconstruction (high-resolution bone reconstruction and soft tissue reconstruction) included 400 sponge tissue images from 50 healthy patients and 50 patients with osteoporosis. Texture feature descriptors were calculated based on the statistical analysis of the grey image histogram, autoregression model, and wavelet transform. The data dimensional reduction was applied by feature selection using nine methods representing various approaches (filter, wrapper, and embedded methods). Eleven methods were used to build the classifier models. In the learning process, hyperparametric optimization based on the grid search method was applied. On this basis, the most effective model and the optimal subset of features for each selection method used were determined. In the case of bone reconstruction images, four models achieved a maximum accuracy of 92%, one of which had the highest sensitivity of 95%, with a specificity of 89%. For soft tissue reconstruction images, five models achieved the highest testing accuracy of 95%, whereas the other quality indices (TPR and TNR) were also equal to 95%. The research showed that the images derived from soft tissue reconstruction allow for obtaining more accurate values of texture parameters, which increases the accuracy of the classification and offers better possibilities for diagnosing osteoporosis.

Comparison of a Deep Learning-Based Reconstruction Algorithm with Filtered Back Projection and Iterative Reconstruction Algorithms for Pediatric Abdominopelvic CT

Objective

To compare a deep learning-based reconstruction (DLR) algorithm for pediatric abdominopelvic computed tomography (CT) with filtered back projection (FBP) and iterative reconstruction (IR) algorithms.

Materials and Methods

Post-contrast abdominopelvic CT scans obtained from 120 pediatric patients (mean age ± standard deviation, 8.7 ± 5.2 years; 60 males) between May 2020 and October 2020 were evaluated in this retrospective study. Images were reconstructed using FBP, a hybrid IR algorithm (ASiR-V) with blending factors of 50% and 100% (AV50 and AV100, respectively), and a DLR algorithm (TrueFidelity) with three strength levels (low, medium, and high). Noise power spectrum (NPS) and edge rise distance (ERD) were used to evaluate noise characteristics and spatial resolution, respectively. Image noise, edge definition, overall image quality, lesion detectability and conspicuity, and artifacts were qualitatively scored by two pediatric radiologists, and the scores of the two reviewers were averaged. A repeated-measures analysis of variance followed by the Bonferroni post-hoc test was used to compare NPS and ERD among the six reconstruction methods. The Friedman rank sum test followed by the Nemenyi-Wilcoxon-Wilcox all-pairs test was used to compare the results of the qualitative visual analysis among the six reconstruction methods.

Results

The NPS noise magnitude of AV100 was significantly lower than that of the DLR, whereas the NPS peak of AV100 was significantly higher than that of the high- and medium-strength DLR (p < 0.001). The NPS average spatial frequencies were higher for DLR than for ASiR-V (p < 0.001). ERD was shorter with DLR than with ASiR-V and FBP (p < 0.001). Qualitative visual analysis revealed better overall image quality with high-strength DLR than with ASiR-V (p < 0.001).

Conclusion

For pediatric abdominopelvic CT, the DLR algorithm may provide improved noise characteristics and better spatial resolution than the hybrid IR algorithm.

Harmonization of technical image quality in computed tomography: comparison between different reconstruction algorithms and kernels from six scanners

Purpose. The radiology department faces a large number of reconstruction algorithms and kernels during their computed tomography (CT) optimization process. These reconstruction methods are proprietary and ensuring consistent image quality between scanners is becoming increasingly difficult. This study contributes to solving this challenge in CT image quality harmonization by modifying and evaluating a reconstruction algorithm and kernel matching scheme.Methods. The Catphan 600 phantom was scanned with six different CT scanners from four vendors. The phantom was scanned with volumetric CT dose indices (CTDIvols) of 10 mGy and 40 mGy, and the data were reconstructed using 1 mm and 5 mm slices with each combination of reconstruction algorithm, body region kernel, and iterative and deep learning reconstruction strength. A matching scheme developed in previous research, which utilizes the noise power spectrum (NPS) and modulation transfer function (MTF), was modified based on our organization's needs and used to identify the matching reconstruction algorithms and kernels between different scanners.Results. The matching paradigm produced good matching results, and the mean ± standard deviation (median) matching function values for the different acquisition settings were (a value of 1 indicates a perfect match): CTDIvol 10 mGy, 1 mm slice: 0.78 ± 0.31 (0.94); CTDIvol 10 mGy, 5 mm slice: 0.75 ± 0.33 (0.93); CTDIvol 40 mGy, 1 mm slice: 0.81 ± 0.28 (0.95); CTDIvol 40 mGy, 5 mm slice: 0.75 ± 0.33 (0.93). In general, soft reconstruction kernels, i.e., noise-reducing kernels that reduce sharpness, of one vendor were matched with the soft kernels of another vendor, and vice versa for sharper kernels. Conclusions. Combined quantitative assessment of NPS and MTF allows effective strategy for harmonization of technical image quality between different CT scanners. A software was also shared to support CT image quality harmonization in other institutions.

Evaluation of Apparent Noise on CT Images Using Moving Average Filters

This study aims to devise a simple method for evaluating the magnitude of texture noise (apparent noise) observed on computed tomography (CT) images scanned at a low radiation dose and reconstructed using iterative reconstruction (IR) and deep learning reconstruction (DLR) algorithms, and to evaluate the apparent noise in CT images reconstructed using the filtered back projection (FBP), IR, and two types of DLR (AiCE Body and AiCE Body Sharp) algorithms. We set a square region of interest (ROI) on CT images of standard- and obese-sized low-contrast phantoms, slid different-sized moving average filters in the ROI vertically and horizontally in steps of 1 pixel, and calculated the standard deviation (SD) of the mean CT values for each filter size. The SD of the mean CT values was fitted with a curve inversely proportional to the filter size, and an apparent noise index was determined from the curve-fitting formula. The apparent noise index of AiCE Body Sharp images for a given mAs value was approximately 58, 23, and 18% lower than that of the FBP, AIDR 3D, and AiCE Body images, respectively. The apparent noise index was considered to reflect noise power spectrum values at lower spatial frequency. Moreover, the apparent noise index was inversely proportional to the square roots of the mAs values. Thus, the apparent noise index could be a useful indicator to quantify and compare texture noise on CT images obtained with different scan parameters and reconstruction algorithms.

Method to estimate fan-beam CT noise power spectrum using two basis functions with a limited number of noise realizations

PURPOSE

The noise power spectrum (NPS) plays a key role in image quality (IQ) evaluation as it can be used for predicting detection performance or calculating detective quantum efficiency (DQE). Traditionally, the NPS is estimated by ensemble averaging multiple realizations of noise-only images. However, the estimation error increases when there are a limited number of images. Since the estimation error directly affects the IQ index, an accurate NPS estimation method is required. Recent works have proposed NPS estimation methods using the radial 1D NPS as the basis; however, when sharp kernels are used during image reconstruction, these methods cannot accurately estimate the amplitude of each angular spoke of the 2D NPS composed of different cutoff frequencies determined from the complementary projection magnification factors for different spatial regions. In this work, we propose a 2D NPS estimation method that reflects the accurate amplitude of each angular spoke for fan-beam CT images.

METHODS

An angular spoke of the 2D NPS is composed of two basis functions with different cutoff frequencies determined from the complementary projection magnification factors. The proposed method estimates these two weighting factors for each basis function by minimizing the mean-squared error (MSE) between the 2D NPS estimated from 10 noise realizations. Two noise profiles and two types of apodization filters (i.e., rectangular and Hanning) were used to reconstruct the noise-only images. To examine the nonstationary noise property of fan-beam CT images, the 2D NPS was estimated at three different local regions. The estimation accuracy of the proposed method was further improved by estimating the approximate weighting factors with sinusoidal functions, considering that the weighting factors vary slowly throughout the view angles. Regression orders of 1 to 4 were used during these estimations. The approximate weighting factors were then multiplied with each of the basis functions to estimate the 2D NPS. The normalized mean-squared error (NMSE) was used as an index to compare the performance of each NPS estimation method, with the analytical 2D NPS as the reference. Further validation was performed using XCAT phantom data.

RESULTS

We observed that the 2D NPS estimated using two basis functions reflected the accurate amplitude of each angular spoke, whereas the 2D NPS estimated using the radial 1D NPS as the basis could not. The 2D NPS estimated by applying the approximate weighting factors showed improved performance compared with that estimated using two basis functions. In addition, unlike the view-independent noise cases, where a lower regression order showed higher estimation performance, a higher regression order showed higher estimation performance in the view-dependent noise cases.

CONCLUSIONS

In this work, we propose a 2D NPS estimation method that reflects the accurate amplitude of each angular spoke for fan-beam CT images using two basis functions. We observed that the proposed 2D NPS estimation method using two basis functions achieved better estimation performance compared with the method using the radial 1D NPS as the basis. This article is protected by copyright. All rights reserved.

Deep-learning model observer for a low-contrast hepatic metastases localization task in computed tomography

PURPOSE

Conventional model observers (MO) in CT are often limited to a uniform background or varying background that is random and can be modeled in an analytical form. It is unclear if these conventional MOs can be readily generalized to predict human observer performance in clinical CT tasks that involve realistic anatomical background. Deep-learning-based model observers (DL-MO) have recently been developed, but have not been validated for challenging low contrast diagnostic tasks in abdominal CT. We consequently sought to validate a DL-MO for a low-contrast hepatic metastases localization task.

METHODS

We adapted our recently developed DL-MO framework for the liver metastases localization task. Our previously-validated projection-domain lesion-/noise-insertion techniques were used to synthesize realistic positive and low-dose abdominal CT exams, using the archived patient projection data. Ten experimental conditions were generated, which involved different lesion sizes/contrasts, radiation dose levels, and image reconstruction types. Each condition included 100 trials generated from a patient cohort of 7 cases. Each trial was presented as liver image patches (160×160×5 voxels). The DL-MO performance was calculated for each condition and was compared with human observer performance, which was obtained by three sub-specialized radiologists in an observer study. The performance of DL-MO and radiologists was gauged by the area under localization receiver-operating-characteristic curves. The generalization performance of the DL-MO was estimated with the repeated twofold cross-validation method over the same set of trials used in the human observer study. A multi-slice Channelized Hoteling Observers (CHO) was compared with the DL-MO across the same experimental conditions.

RESULTS

The performance of DL-MO was highly correlated to that of radiologists (Pearson's correlation coefficient: 0.987; 95% CI: [0.942, 0.997]). The performance level of DL-MO was comparable to that of the grouped radiologists, that is, the mean performance difference was -3.3%. The CHO performance was poorer than the grouped radiologist performance, before internal noise could be added. The correlation between CHO and radiologists was weaker (Pearson's correlation coefficient: 0.812, and 95% CI: [0.378, 0.955]), and the corresponding performance bias (-29.5%) was statistically significant.

CONCLUSION

The presented study demonstrated the potential of using the DL-MO for image quality assessment in patient abdominal CT tasks.

The Discriminative Power and Stability of Radiomics Features With Computed Tomography Variations: Task-Based Analysis in an Anthropomorphic 3D-Printed CT Phantom

OBJECTIVES

The aims of this study were to determine the stability of radiomics features against computed tomography (CT) parameter variations and to study their discriminative power concerning tissue classification using a 3D-printed CT phantom based on real patient data.

MATERIALS AND METHODS

A radiopaque 3D phantom was developed using real patient data and a potassium iodide solution paper-printing technique. Normal liver tissue and 3 lesion types (benign cyst, hemangioma, and metastasis) were manually annotated in the phantom. The stability and discriminative power of 86 radiomics features were assessed in measurements taken from 240 CT series with 8 parameter variations of reconstruction algorithms, reconstruction kernels, slice thickness, and slice spacing. Pairwise parameter group and pairwise tissue class comparisons were performed using Wilcoxon signed rank tests.

RESULTS

In total, 19,264 feature stability tests and 8256 discriminative power tests were performed. The 8 CT parameter variation pairwise group comparisons had statistically significant differences on average in 78/86 radiomics features. On the other hand, 84% of the univariate radiomics feature tests had a successful and statistically significant differentiation of the 4 classes of liver tissue. The 86 radiomics features were ranked according to the cumulative sum of successful stability and discriminative power tests.

CONCLUSIONS

The differences in radiomics feature values obtained from different types of liver tissue are generally greater than the intraclass differences resulting from CT parameter variations.

Three-dimensional modeling of the mitral valve for surgical planning in a pediatric patient: A case-based discussion of the technical challenges of segmentation and printing from 3D transthoracic echocardiographic datasets

Abnormal atrioventricular valve present great challenges to the surgeon in achieving a successful repair, and thus present a great opportunity for enhanced 3D imaging to guide pre- and intra-operative management. Spatial and temporal resolution of 3D echocardiography enables 3D printing of valve morphology. However, non-linearity, angle dependence, speckle, blur, and resampling complicate segmentation compared to computed tomography (CT) and magnetic resonance imaging (MRI). A case of complex mitral valve disease in a pediatric patient is therefore presented to illustrate the technical challenges of segmentation and 3D printing from echocardiographic data.

Comparison of low-contrast detectability between uniform and anatomically realistic phantoms-influences on CT image quality assessment

Objectives

To evaluate the effects of anatomical phantom structure on task-based image quality assessment compared with a uniform phantom background.

Methods

Two neck phantom types of identical shape were investigated: a uniform type containing 10-mm lesions with 4, 9, 18, 30, and 38 HU contrast to the surrounding area and an anatomically realistic type containing lesions of the same size and location with 10, 18, 30, and 38 HU contrast. Phantom images were acquired at two dose levels (CTDIvol of 1.4 and 5.6 mGy) and reconstructed using filtered back projection (FBP) and adaptive iterative dose reduction 3D (AIDR 3D). Detection accuracy was evaluated by seven radiologists in a 4-alternative forced choice experiment.

Results

Anatomical phantom structure impaired lesion detection at all lesion contrasts (p < 0.01). Detectability in the anatomical phantom at 30 HU contrast was similar to 9 HU contrast in uniform images (91.1% vs. 89.5%). Detection accuracy decreased from 83.6% at 5.6 mGy to 55.4% at 1.4 mGy in uniform FBP images (p < 0.001), whereas AIDR 3D preserved detectability at 1.4 mGy (80.7% vs. 85% at 5.6 mGy, p = 0.375) and was superior to FBP (p < 0.001). In the assessment of anatomical images, superiority of AIDR 3D was not confirmed and dose reduction moderately affected detectability (74.6% vs. 68.2%, p = 0.027 for FBP and 81.1% vs. 73%, p = 0.018 for AIDR 3D).

Conclusions

A lesion contrast increase from 9 to 30 HU is necessary for similar detectability in anatomical and uniform neck phantom images. Anatomical phantom structure influences task-based assessment of iterative reconstruction and dose effects.

Key Points

• A lesion contrast increase from 9 to 30 HU is necessary for similar low-contrast detectability in anatomical and uniform neck phantom images.

• Phantom background structure influences task-based assessment of iterative reconstruction and dose effects.

• Transferability of CT assessment to clinical imaging can be expected to improve as the realism of the test environment increases.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-021-08248-3.

Development and validation of an automated methodology to assess perceptual in vivo noise texture in liver CT

Purpose: Developing, validating, and evaluating a method for measuring noise texture directly from patient liver CT images (i.e., in vivo). Approach: The method identifies target regions within patient scans that are least likely to have major contribution of patient anatomy, detrends them locally, and measures noise power spectrum (NPS) there using a previously phantom-validated technique targeting perceptual noise-non-anatomical fluctuations in the image that may interfere with the detection of focal lesions. Method development and validation used scanner-specific CT simulations of computational, anthropomorphic phantom (XCAT phantom, three phases of contrast-enhancement) with known ground truth of the NPS. Simulations were based on a clinical scanner (Definition Flash, Siemens) and clinically relevant settings (tube voltage of 120 kV at three dose levels). Images were reconstructed with filtered backprojection (kernel: B31, B41, and B50) and Sinogram Affirmed Iterative Reconstruction (kernel: I31, I41, and I50) using a manufacturer-specific reconstruction software (ReconCT, Siemens). All NPS measurements were made in the liver. Ground-truth NPS were taken as the sum of (1) a measurement in parenchymal regions of anatomy-subtracted (i.e., noise only) scans, and (2) a measurement in the same region of noise-free (pre-noise-insertion) images. To assess in vivo NPS performance, correlation of NPS average frequency (f avg ), was reported. Sensitivity of accuracy [root-mean-square-error (RMSE)] to number of pixels included in measurement was conducted via bootstrapped pixel-dropout. Sensitivity of NPS to dose and reconstruction kernel was assessed to confirm that ground truth NPS similarities were maintained in patient-specific measurements. Results: Pearson and Spearman correlation coefficients 0.97 and 0.96 forf avg indicated good correlation. Results suggested accurate NPS measurements (within 5% total RMSE) could be acquired with ∼ 10 6    pixels . Conclusions: Relationships of similar NPS due to reconstruction kernel and dose were preserved between gold standard and observed in vivo estimations. The NPS estimation method was further deployed on clinical cases to demonstrate the feasibility of clinical analysis.

Low dose cone beam CT for paediatric image-guided radiotherapy: Image quality and practical recommendations

PURPOSE

Cone beam CT (CBCT) is used in paediatric image-guided radiotherapy (IGRT) for patient setup and internal anatomy assessment. Adult CBCT protocols lead to excessive doses in children, increasing the risk of radiation-induced malignancies. Reducing imaging dose increases quantum noise, degrading image quality. Patient CBCTs also include 'anatomical noise' (e.g. motion artefacts), further degrading quality. We determine noise contributions in paediatric CBCT, recommending practical imaging protocols and thresholds above which increasing dose yields no improvement in image quality.

METHODS AND MATERIALS

Sixty CBCTs including the thorax or abdomen/pelvis from 7 paediatric patients (aged 6-13 years) were acquired at a range of doses and used to simulate lower dose scans, totalling 192 scans (0.5-12.8 mGy). Noise measured in corresponding regions of each patient and a 10-year-old phantom were compared, modelling total (including anatomical) noise, and quantum noise contributions as a function of dose. Contrast-to-noise ratio (CNR) was measured between fat/muscle. Soft tissue registration was performed on the kidneys, comparing accuracy to the highest dose scans.

RESULTS

Quantum noise contributed <20% to total noise in all cases, suggesting anatomical noise is the largest determinant of image quality in the abdominal/pelvic region. CNR exceeded 3 in over 90% of cases ≥ 1 mGy, and 57% of cases at 0.5 mGy. Soft tissue registration was accurate for doses > 1 mGy.

CONCLUSION

Anatomical noise dominates quantum noise in paediatric CBCT. Appropriate soft tissue contrast and registration accuracy can be achieved for doses as low as 1 mGy. Increasing dose above 1 mGy holds no benefit in improving image quality or registration accuracy due to the presence of anatomical noise.

Impact of iterative reconstruction algorithms on the applicability of Fourier-based detectability index for x-ray CT imaging

PURPOSE

The increasing application of iterative reconstruction algorithms in clinical computed tomography to improve image quality and reduce radiation dose, elicits strong interest, and needs model observers to optimize CT scanning protocols objectively and efficiently. The current paradigm for evaluating imaging system performance relies on Fourier methods, which presuppose a linear, wide-sense stationary system. Long-range correlations introduced by iterative reconstruction algorithms may narrow the applicability of Fourier techniques. Differences in the implementation of reconstruction algorithms between manufacturers add further complexity. The present work set out to quantify the errors entailed by the use of Fourier methods, which can lead to design decisions that do not correlate with detectability.

METHODS

To address this question, we evaluated the noise properties and the detectability index of the ideal linear observer using the spatial approach and the Fourier-based approach. For this purpose, a homogeneous phantom was imaged on two scanners: the Revolution CT (GE Healthcare) and the Somatom Definition AS+ (Siemens Healthcare) at different exposure levels. Images were reconstructed using different strength levels of IR algorithms available on the systems considered: Adaptative Statistical Iterative Reconstruction (ASIR-V) and Sinogram Affirmed Iterative Reconstruction (SAFIRE).

RESULTS

Our findings highlight that the spatial domain estimate of the detectability index is higher than the Fourier domain estimate. This trend is found to be dependent on the specific regularization used by IR algorithms as well as the signal to be detected. The eigenanalysis of the noise covariance matrix and of its circulant approximation yields explanation about the evoked trends. In particular, this analysis suggests that the predictive power of the Fourier-based ideal linear observer depends on the ability of each basis analyzed to be relevant to the signal to be detected.

CONCLUSION

The applicability of Fourier techniques is dependent on the specific regularization used by IR algorithms. These results argue for verifying the assumptions made when using Fourier methods since Fourier-task-based detectability index does not always correlate with signal detectability.

Quality control of radiomic features using 3D-printed CT phantoms

Purpose: The lack of standardization in quantitative radiomic measures of tumors seen on computed tomography (CT) scans is generally recognized as an unresolved issue. To develop reliable clinical applications, radiomics must be robust across different CT scan modes, protocols, software, and systems. We demonstrate how custom-designed phantoms, imprinted with human-derived patterns, can provide a straightforward approach to validating longitudinally stable radiomic signature values in a clinical setting. Approach: Described herein is a prototype process to design an anatomically informed 3D-printed radiomic phantom. We used a multimaterial, ultra-high-resolution 3D printer with voxel printing capabilities. Multiple tissue regions of interest (ROIs), from four pancreas tumors, one lung tumor, and a liver background, were extracted from digital imaging and communication in medicine (DICOM) CT exam files and were merged together to develop a multipurpose, circular radiomic phantom (18 cm diameter and 4 cm width). The phantom was scanned 30 times using standard clinical CT protocols to test repeatability. Features that have been found to be prognostic for various diseases were then investigated for their repeatability and reproducibility across different CT scan modes. Results: The structural similarity index between the segment used from the patients' DICOM image and the phantom CT scan was 0.71. The coefficient variation for all assessed radiomic features was < 1.0 % across 30 repeat scans of the phantom. The percent deviation (pDV) from the baseline value, which was the mean feature value determined from repeat scans, increased with the application of the lung convolution kernel, changes to the voxel size, and increases in the image noise. Gray level co-occurrence features, contrast, dissimilarity, and entropy were particularly affected by different scan modes, presenting with pDV > ± 15 % . Conclusions: Previously discovered prognostic and popular radiomic features are variable in practice and need to be interpreted with caution or excluded from clinical implementation. Voxel-based 3D printing can reproduce tissue morphology seen on CT exams. We believe that this is a flexible, yet practical, way to design custom phantoms to validate and compare radiomic metrics longitudinally, over time, and across systems.

Protocol Optimization Considerations for Implementing Deep Learning CT Reconstruction

OBJECTIVE. Previous advances over filtered back projection (FBP) have incorporated model-based iterative reconstruction. The purpose of this study was to characterize the latest advance in image reconstruction, that is, deep learning. The focus was on applying characterization results of a deep learning approach to decisions about clinical CT protocols. MATERIALS AND METHODS. A proprietary deep learning image reconstruction (DLIR) method was characterized against an existing advanced adaptive statistical iterative reconstruction method (ASIR-V) and FBP from the same vendor. The metrics used were contrast-to-noise ratio, spatial resolution as a function of contrast level, noise texture (i.e., noise power spectra [NPS]), noise scaling as a function of slice thickness, and CT number consistency. The American College of Radiology accreditation phantom and a uniform water phantom were used at a range of doses and slice thicknesses for both axial and helical acquisition modes. RESULTS. ASIR-V and DLIR were associated with improved contrast-to-noise ratio over FBP for all doses and slice thicknesses. No dose or contrast dependencies of spatial resolution were observed for ASIR-V or DLIR. NPS results showed DLIR maintained an FBP-like noise texture whereas ASIR-V shifted the NPS to lower frequencies. Noise changed with dose and slice thickness in the same manner for ASIR-V and FBP. DLIR slice thickness noise scaling differed from FBP, exhibiting less noise penalty with decreasing slice thickness. No clinically significant changes were observed in CT numbers for any measurement condition. CONCLUSION. In a phantom model, DLIR does not suffer from the concerns over reduction in spatial resolution and introduction of poor noise texture associated with previous methods.

A patient-specific approach for quantitative and automatic analysis of computed tomography images in lung disease: Application to COVID-19 patients

PURPOSE

Quantitative metrics in lung computed tomography (CT) images have been widely used, often without a clear connection with physiology. This work proposes a patient-independent model for the estimation of well-aerated volume of lungs in CT images (WAVE).

METHODS

A Gaussian fit, with mean (Mu.f) and width (Sigma.f) values, was applied to the lower CT histogram data points of the lung to provide the estimation of the well-aerated lung volume (WAVE.f). Independence from CT reconstruction parameters and respiratory cycle was analysed using healthy lung CT images and 4DCT acquisitions. The Gaussian metrics and first order radiomic features calculated for a third cohort of COVID-19 patients were compared with those relative to healthy lungs. Each lung was further segmented in 24 subregions and a new biomarker derived from Gaussian fit parameter Mu.f was proposed to represent the local density changes.

RESULTS

WAVE.f resulted independent from the respiratory motion in 80% of the cases. Differences of 1%, 2% and up to 14% resulted comparing a moderate iterative strength and FBP algorithm, 1 and 3 mm of slice thickness and different reconstruction kernel. Healthy subjects were significantly different from COVID-19 patients for all the metrics calculated. Graphical representation of the local biomarker provides spatial and quantitative information in a single 2D picture.

CONCLUSIONS

Unlike other metrics based on fixed histogram thresholds, this model is able to consider the inter- and intra-subject variability. In addition, it defines a local biomarker to quantify the severity of the disease, independently of the observer.

Effect of obesity on ability to lower exposure for detection of low-attenuation liver lesions

Purpose

The purpose of this study was to assess the effect of obesity and iterative reconstruction on the ability to reduce exposure by studying the accuracy for detection of low‐contrast low‐attenuation (LCLA) liver lesions on computed tomography (CT) using a phantom model.

Methods

A phantom with four unique LCLA liver lesions (5‐ to 15‐mm spheres, –24 to –6 HU relative to 90‐HU background) was scanned without (“thin” phantom) and with (“obese” phantom) a 5‐cm thick fat‐attenuation ring at 150 mAs (thin phantom) and 450 mAs (obese phantom) standard exposures and at 33% and 67% exposure reductions. Images were reconstructed using standard filtered back projection (FBP) and with iterative reconstruction (Adaptive Model‐Based Iterative Reconstruction strength 3, ADMIRE). A noninferiority analysis of lesion detection was performed.

Results

Mean area under the curve (AUC) values for lesion detection were significantly higher for the thin phantom than for the obese phantom regardless of exposure level (P < 0.05) for both FBP and ADMIRE. At 33% exposure reduction, AUC was noninferior for both FBP and ADMIRE strength 3 (P < 0.0001). At 67% exposure reduction, AUC remained noninferior for the thin phantom (P < 0.0035), but was no longer noninferior for the obese phantom (P ≥ 0.7353). There were no statistically significant differences in AUC between FBP and ADMIRE at any exposure level for either phantom.

Conclusions

Accuracy for lesion detection was not only significantly lower in the obese phantom at all relative exposures, but detection accuracy decreased sooner while reducing the exposure in the obese phantom. There was no significant difference in lesion detection between FBP and ADMIRE at equivalent exposure levels for either phantom.

Task-based assessment of neck CT protocols using patient-mimicking phantoms-effects of protocol parameters on dose and diagnostic performance

Objectives

To assess how modifying multiple protocol parameters affects the dose and diagnostic performance of a neck CT protocol using patient-mimicking phantoms and task-based methods.

Methods

Six patient-mimicking neck phantoms containing hypodense lesions of 1 cm diameter and 30 HU contrast and one non-lesion phantom were examined with 36 CT protocols. All possible combinations of the following parameters were investigated: 100- and 120-kVp tube voltage; tube current modulation (TCM) noise levels of SD 7.5, 10, and 14; pitches of 0.637, 0.813, and 1.388; filtered back projection (FBP); and iterative reconstruction (AIDR 3D). Dose-length products (DLPs) and lesion detectability (assessed by 14 radiologists) were compared with the clinical standard protocol (120 kVp, TCM SD 7.5, 0.813 pitch, AIDR 3D).

Results

The DLP of the standard protocol was 25 mGy•cm; the area under the curve (AUC) was 0.839 (95%CI: 0.790–0.888). Combined effects of tube voltage reduction to 100 kVp and TCM noise level increase to SD 10 optimized protocol performance by improving dose (7.3 mGy•cm) and detectability (AUC 0.884, 95%CI: 0.844–0.924). Diagnostic performance was significantly affected by the TCM noise level at 120 kVp (AUC 0.821 at TCM SD 7.5 vs. 0.776 at TCM SD 14, p = 0.003), but not at 100-kVp tube voltage (AUC 0.839 at TCM SD 7.5 vs. 0.819 at TCM SD 14, p = 0.354), the reconstruction method at 100 kVp (AUC 0.854 for AIDR 3D vs. 0.806 for FBP, p < 0.001), but not at 120-kVp tube voltage (AUC 0.795 for AIDR 3D vs. 0.793 for FBP, p = 0.822), and the tube voltage for AIDR 3D reconstruction (p < 0.001), but not for FBP (p = 0.226).

Conclusions

Combined effects of 100-kVp tube voltage, TCM noise level of SD 10, a pitch of 0.813, and AIDR 3D resulted in an optimal neck protocol in terms of dose and diagnostic performance. Protocol parameters were subject to complex interactions, which created opportunities for protocol improvement.

Key Points

• A task-based approach using patient-mimicking phantoms was employed to optimize a CT system for neck imaging through systematic testing of protocol parameters.

• Combined effects of 100-kVp tube voltage, TCM noise level of SD 10, a pitch of 0.813, and AIDR 3D reconstruction resulted in an optimal protocol in terms of dose and diagnostic performance.

• Interactions of protocol parameters affect diagnostic performance and should be considered when optimizing CT techniques.

Electronic supplementary material

The online version of this article (10.1007/s00330-020-07374-8) contains supplementary material, which is available to authorized users.

Spatial Distribution of Noise Reduction in Four Iterative Reconstruction Algorithms in CT—A Technical Evaluation

Iterative reconstruction (IR) is a computed tomgraphy (CT) reconstruction algorithm aiming at improving image quality by reducing noise in the image. During this process, IR also changes the noise properties in the images. To assess how IR algorithms from four vendors affect the noise properties in CT images, an anthropomorphic phantom was scanned and images reconstructed with filtered back projection (FBP), and a medium and high level of IR. Each image acquisition was performed 30 times at the same slice position, to create noise maps showing the inter-image pixel standard deviation through the 30 images. We observed that IR changed the noise properties in the CT images by reducing noise more in homogeneous areas than at anatomical edges between structures of different densities. This difference increased with increasing IR level, and with increasing difference in density between two adjacent structures. Each vendor’s IR algorithm showed slightly different noise reduction properties in how much noise was reduced at different positions in the phantom. Users need to be aware of these differences when working with optimization of protocols using IR across scanners from different vendors.

Noise and spatial resolution properties of a commercially available deep learning-based CT reconstruction algorithm

PURPOSE

To characterize the noise and spatial resolution properties of a commercially available deep learning-based computed tomography (CT) reconstruction algorithm.

METHODS

Two phantom experiments were performed. The first used a multisized image quality phantom (Mercury v3.0, Duke University) imaged at five radiation dose levels (CTDI_vol : 0.9, 1.2, 3.6, 7.0, and 22.3 mGy) with a fixed tube current technique on a commercial CT scanner (GE Revolution CT). Images were reconstructed with conventional (FBP), iterative (GE ASiR-V), and deep learning-based (GE True Fidelity) reconstruction algorithms. Noise power spectrum (NPS), high-contrast (air-polyethylene interface), and intermediate-contrast (water-polyethylene interface) task transfer functions (TTF) were measured for each dose level and phantom size and summarized in terms of average noise frequency (f_av ) and frequency at which the TTF was reduced to 50% (f_50% ), respectively. The second experiment used a custom phantom with low-contrast rods and lung texture sections for the assessment of low-contrast TTF and noise spatial distribution. The phantom was imaged at five dose levels (CTDI_vol : 1.0, 2.1, 3.0, 6.0, and 10.0 mGy) with 20 repeated scans at each dose, and images reconstructed with the same reconstruction algorithms. The local noise stationarity was assessed by generating spatial noise maps from the ensemble of repeated images and computing a noise inhomogeneity index, η , following AAPM TG233 methods. All measurements were compared among the algorithms.

RESULTS

Compared to FBP, noise magnitude was reduced on average (± one standard deviation) by 74 ± 6% and 68 ± 4% for ASiR-V (at "100%" setting) and True Fidelity (at "High" setting), respectively. The noise texture from ASiR-V had substantially lower noise frequency content with 55 ± 4% lower NPS f_av compared to FBP while True Fidelity had only marginally different noise frequency content with 9 ± 5% lower NPS f_av compared to FBP. Both ASiR-V and True Fidelity demonstrated locally nonstationary noise in a lung texture background at all radiation dose levels, with higher noise near high-contrast edges of vessels and lower noise in uniform regions. At the 1.0 mGy dose level η values were 314% and 271% higher in ASiR-V and True Fidelity compared to FBP, respectively. High-contrast spatial resolution was similar between all algorithms for all dose levels and phantom sizes (<3% difference in TTF f_50% ). Compared to FBP, low-contrast spatial resolution was lower for ASiR-V and True Fidelity with a reduction of TTF f_50% of up to 42% and 36%, respectively.

CONCLUSIONS

The deep learning-based CT reconstruction demonstrated a strong noise magnitude reduction compared to FBP while maintaining similar noise texture and high-contrast spatial resolution. However, the algorithm resulted in images with a locally nonstationary noise in lung textured backgrounds and had somewhat degraded low-contrast spatial resolution similar to what has been observed in currently available iterative reconstruction techniques.

Creation of an anthropomorphic CT head phantom for verification of image segmentation

Purpose

Many methods are available to segment structural magnetic resonance (MR) images of the brain into different tissue types. These have generally been developed for research purposes but there is some clinical use in the diagnosis of neurodegenerative diseases such as dementia. The potential exists for computed tomography (CT) segmentation to be used in place of MRI segmentation, but this will require a method to verify the accuracy of CT processing, particularly if algorithms developed for MR are used, as MR has notably greater tissue contrast.

Methods

To investigate these issues we have created a three‐dimensional (3D) printed brain with realistic Hounsfield unit (HU) values based on tissue maps segmented directly from an individual T1 MRI scan of a normal subject. Several T1 MRI scans of normal subjects from the ADNI database were segmented using SPM12 and used to create stereolithography files of different tissues for 3D printing. The attenuation properties of several material blends were investigated, and three suitable formulations were used to print an object expected to have realistic geometry and attenuation properties. A skull was simulated by coating the object with plaster of Paris impregnated bandages. Using two CT scanners, the realism of the phantom was assessed by the measurement of HU values, SPM12 segmentation and comparison with the source data used to create the phantom.

Results

Realistic relative HU values were measured although a subtraction of 60 was required to obtain equivalence with the expected values (gray matter 32.9–35.8 phantom, 29.9–34.2 literature). Segmentation of images acquired at different kVps/mAs showed excellent agreement with the source data (Dice Similarity Coefficient 0.79 for gray matter). The performance of two scanners with two segmentation methods was compared, with the scanners found to have similar performance and with one segmentation method clearly superior to the other.

Conclusion

The ability to use 3D printing to create a realistic (in terms of geometry and attenuation properties) head phantom has been demonstrated and used in an initial assessment of CT segmentation accuracy using freely available software developed for MRI.

Evaluation of low-contrast detectability for iterative reconstruction in pediatric abdominal computed tomography: a phantom study

BACKGROUND

Iterative reconstruction is offered by all vendors to achieve similar or better CT image quality at lower doses than images reconstructed with filtered back-projection.

OBJECTIVE

The purpose of this study was to investigate the dose-reduction potential for pediatric abdominal CT imaging when using either a commercially available hybrid or a commercially available model-based iterative reconstruction algorithm from a single manufacturer.

MATERIALS AND METHODS

A phantom containing four low-contrast inserts and a uniform background with total attenuation equivalent to the abdomen of an average 8-year-old child was imaged on a CT scanner (IQon; Philips Healthcare, Cleveland, OH). We reconstructed images using both hybrid iterative reconstruction (iDose⁴) and model-based iterative reconstruction (Iterative Model Reconstruction). The four low-contrast inserts had circular cross-section with diameters of 3 mm, 5 mm, 7 mm and 10 mm and contrasts of 14 Hounsfield units (HU), 7 HU, 5 HU and 3 HU, respectively. Helical scans with identical kilovoltage (kV), pitch, rotation time, and collimation were repeated on the phantom at volume CT dose index (CTDI_vol) of 2.0 milligrays (mGy), 3.0 mGy, 4.5 mGy and 6.0 mGy. We measured the contrast-to-noise ratio (CNR) in each rod across scans. Additionally, we collected sub-images containing each rod and sub-images containing the background and used them in two-alternative forced choice observer experiments with four observers (two radiologists and two physicists). We calculated the dose-reduction potential of both iterative reconstruction algorithms relative to a scan performed at 6 mGy and reconstructed with filtered back-projection.

RESULTS

We calculated dose-reduction potential by either matching average equal observer performance in the two-alternative forced choice experiments or matching CNR. When matching CNR, the dose-reduction potential was 34% to 45% for hybrid iterative reconstruction and 89% to 95% for model-based iterative reconstruction. When matching average observer performance, the dose-reduction potential was 9% to 30% for hybrid iterative reconstruction and 57% to 74% for model-based iterative reconstruction. The range in dose-reduction potential depended on the rod size and contrast level.

CONCLUSION

Observer performance in this phantom study indicates that the dose-reduction potential indicated by an observer study is less than that of CNR; extrapolating the results to clinical studies suggests that the dose-reduction potential would also be less.

Initial clinical experience with high-pitch dual-source CT as a rapid technique for thoraco-abdominal evaluation in awake infants and young children

AIM

To evaluate dual-source high-pitch computed tomography (HPCT) imaging of the chest and abdomen as a rapid scanning technique to obtain diagnostic-quality imaging evaluation of infants and young children without sedation.

MATERIALS AND METHODS

Fifty-three paediatric patients (age 24.1±2 months) who underwent chest or abdomen HPCT (≥1.5) and standard pitch CT (SPCT, <1.5) on a dual-source 128-row multidetector CT system were included in the study. Image quality assessment was performed by two paediatric radiologists for diagnostic confidence, image artefacts, and image noise. Objective image noise was measured.

RESULTS

Most of the CT examinations were performed in children who were >1 year old (n=15 and n=20) followed by ≤1 year old (n=8 and n=10) in SPCT and HPCT, respectively. The mean radiation dose (SSDE) from HPCT was 1.96±1 mGy compared to 2.2±1 mGy for SPCT (p=0.3). No major artefacts were reported and overall image quality of all HPCT examinations was acceptable diagnostically. In addition, objective image noise values were not significantly different between HPCT compared with SPCT (11±3 versus 11±5, p=0.7).

CONCLUSION

Ultra-fast, HPCT can be performed without the need for sedation as a potential alternative to anaesthetised magnetic resonance imaging in infants and young children.

Automatic quantitative analysis of pulmonary vascular morphology in CT images

PURPOSE

Vascular remodeling is a significant pathological feature of various pulmonary diseases, which may be assessed by quantitative computed tomography (CT) imaging. The purpose of this study was therefore to develop and validate an automatic method for quantifying pulmonary vascular morphology in CT images.

METHODS

The proposed method consists of pulmonary vessel extraction and quantification. For extracting pulmonary vessels, a graph-cuts-based method is proposed which considers appearance (CT intensity) and shape (vesselness from a Hessian-based filter) features, and incorporates distance to the airways into the cost function to prevent false detection of airway walls. For quantifying the extracted pulmonary vessels, a radius histogram is generated by counting the occurrence of vessel radii, calculated from a distance transform-based method. Subsequently, two biomarkers, slope α and intercept β, are calculated by linear regression on the radius histogram. A public data set from the VESSEL12 challenge was used to independently evaluate the vessel extraction. The quantitative analysis method was validated using images of a three-dimensional (3D) printed vessel phantom, scanned by a clinical CT scanner and a micro-CT scanner (to obtain a gold standard). To confirm the association between imaging biomarkers and pulmonary function, 77 scleroderma patients were investigated with the proposed method.

RESULTS

In the independent evaluation with the public data set, our vessel segmentation method obtained an area under the receiver operating characteristic (ROC) curve of 0.976. The median radius difference between clinical and micro-CT scans of a 3D printed vessel phantom was 0.062 ± 0.020 mm, with interquartile range of 0.199 ± 0.050 mm. In the studied patient group, a significant correlation between diffusion capacity for carbon monoxide and the biomarkers, α (R = -0.27, P = 0.018) and β (R = 0.321, P = 0.004), was obtained.

CONCLUSION

In conclusion, the proposed method was validated independently using a public data set resulting in an area under the ROC curve of 0.976 and using a 3D printed vessel phantom data set, showing a vessel sizing error of 0.062 mm (0.16 in-plane pixel units). The correlation between imaging biomarkers and diffusion capacity in a clinical data set confirmed an association between lung structure and function. This quantification of pulmonary vascular morphology may be helpful in understanding the pathophysiology of pulmonary vascular diseases.

Design and fabrication of heterogeneous lung nodule phantoms for assessing the accuracy and variability of measured texture radiomics features in CT

We aimed to design and fabricate synthetic lung nodules with patient-informed internal heterogeneity to assess the variability and accuracy of measured texture features in CT. To that end, 190 lung nodules from a publicly available database of chest CT images (Lung Image Database Consortium) were selected based on size ( > 3    mm ) and malignancy. The texture features of the nodules were used to train a statistical texture synthesis model based on clustered lumpy background. The model parameters were ascertained based on a genetic optimization of a Mahalanobis distance objective function. The resulting texture model defined internal heterogeneity within 24 anthropomorphic lesion models which were subsequently fabricated into physical phantoms using a multimaterial three-dimensional (3-D) printer. The 3-D-printed lesions were imbedded in an anthropomorphic chest phantom and imaged with a clinical scanner using different acquisition parameters including slice thickness, dose level, and reconstruction kernel. The imaged lesions were analyzed in terms of texture features to ascertain the impact of CT imaging on lesion texture quantification. The texture modeling method produced lesion models with low and stable Mahalanobis distance between real and synthetic textures. The virtual lesions were successfully printed into 3-D phantoms. The accuracy and variability of the measured features extracted from the CT images of the phantoms showed notable influence from the imaging acquisition parameters with contrast, energy, and texture entropy exhibiting most sensitivity in terms of accuracy, and contrast, dissimilarity, and texture entropy most variability. Thinner slice thicknesses yielded more accurate and edge reconstruction kernels more stable results. We conclude that printed textured models of lesions can be developed using a method that can target and minimize the mathematical distance between real and synthetic lesions. The synthetic lesions can be used as the basis to investigate how CT imaging conditions might affect radiomics features derived from CT images.

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.

Time stability of delta-radiomics features and the impact on patient analysis in longitudinal CT images

PURPOSE

This study first aims to show that the values of texture features extracted from phantoms are stable over clinical timescales. Second, that changes in patients' feature values over the course of radiation therapy (RT) are treatment induced and statistically significant.

METHODS

The CT datasets of a 3D printed anatomically informed texture phantom containing liver and low-contrast modules, and the homogeneous module of the Catphan 500-Series phantom, were acquired once per week over the course of a 6-week period, to simulate the timescale of conventional RT duration. A Definition AS Open CT scanner on rails (Siemens) and our institution's standard abdominal protocol were used. In each phantom module, 8 regions of interest (20 cm 3 ) were selected and 50 texture features were extracted from each module over the longitudinal dataset. The time stability of each feature was evaluated. The expected variation over the treatment timescale was quantified for each texture (module). Subsequently, the pancreas heads of 10 patients who underwent RT for adenocarcinoma of the pancreas head with a pathologic response of at least "moderate" (grade 2), were contoured on the daily CTs acquired using the same scanner. The pancreas heads were contoured on one image per week. Mean CT number, skewness, kurtosis, and coarseness were extracted from these data. The phantom modules were shown to be accurate representations of these features in the pancreas data. The change in the feature value between fractions 2 and 26 was compared with the phantom data in order to identify significant changes in feature value.

RESULTS

Of the 50 features examined in all 3 phantom modules, 47 were found to have zero time-trend when a fit assuming homogeneous variance was used. When a fit allowing for heterogeneous variance was used, 49 features were found to have zero time-trend. Features were stable and repeatable within a feature-specific confidence interval over the 6-week period of acquisition in all three phantom modules. Changes in feature value between fractions 2 and 26 were highly patient specific. Mean CT number was found to decrease significantly in 7 of 10 patients and increase significantly in one patient. Skewness increased significantly in one patient and decreased significantly in one patient. Kurtosis decreased significantly in four patients and increased significantly in one patient. Coarseness increased significantly in seven patients and decreased significantly in one patient. Only one patient experienced no significant changes in feature value.

CONCLUSION

The CT texture feature measurements of phantoms are stable and repeatable within a feature-specific confidence interval in all three phantom modules. This suggests that the changes observed in features extracted from longitudinal patient CT data may be treatment induced, and demonstrates their potentiality for early assessment of treatment response.

Predicting image properties in penalized-likelihood reconstructions of flat-panel CBCT

PURPOSE

Model-based iterative reconstruction (MBIR) algorithms such as penalized-likelihood (PL) methods exhibit data-dependent and shift-variant properties. Image quality predictors have been derived to prospectively estimate local noise and spatial resolution, facilitating both system hardware design and tuning of reconstruction methods. However, current MBIR image quality predictors rely on idealized system models, ignoring physical blurring effects and noise correlations found in real systems. In this work, we develop and validate a new set of predictors using a physical system model specific to flat-panel cone-beam CT (FP-CBCT).

METHODS

Physical models appropriate for integration with MBIR analysis are developed and parameterized to represent nonidealities in FP projection data including focal spot blur, scintillator blur, detector aperture effect, and noise correlations. Flat-panel-specific predictors for local spatial resolution and local noise properties in PL reconstructions are developed based on these realistic physical models. Estimation accuracy of conventional (idealized) and FP-specific predictors is investigated and validated against experimental CBCT measurements using specialized phantoms.

RESULTS

Validation studies show that flat-panel-specific predictors can accurately estimate the local spatial resolution and noise properties, while conventional predictors show significant deviations in the magnitude and scale of the spatial resolution and local noise. The proposed predictors show accurate estimations over a range of imaging conditions including varying x-ray technique and regularization strength. The conventional spatial resolution prediction is sharper than ground truth. Using conventional spatial resolution predictor, the full width at half maximum (FWHM) of local point spread function (PSF) is underestimated by 0.2 mm. This mismatch is mostly eliminated in FP-specific prediction. The general shape and amplitude of local noise power spectrum (NPS) FP-specific predictions are consistent with measurement, while the conventional predictions underestimated the noise level by 70%.

CONCLUSION

The proposed image quality predictors permit accurate estimation of local spatial resolution and noise properties for PL reconstruction, accounting for dependencies on the system geometry, x-ray technique, and patient-specific anatomy in real FP-CBCT. Such tools enable prospective analysis of image quality for a range of goals including novel system and acquisition design, adaptive and task-driven imaging, and tuning of MBIR for robust and reliable behavior.

The Use of Model-based Iterative Reconstruction to Optimize Chest CT Examinations for Diagnosing Lung Metastases in Patients with Sarcoma: A Phantom Study

RATIONALE AND OBJECTIVES

This phantom study aimed to evaluate low-dose (LD) chest computed tomography (CT) protocols using model-based iterative reconstruction (MBIR) for diagnosing lung metastases in patients with sarcoma.

MATERIALS AND METHODS

An adult female anthropomorphic phantom was scanned with a 64-slice CT using four LD protocols and a standard-dose protocol. Absorbed organ doses were measured with 10 metal-oxide-semiconductor field-effect transistor dosimeters. Furthermore, Monte Carlo simulations were performed to estimate organ and effective doses. Image quality in terms of image noise, contrast, and resolution was measured from the CT images reconstructed with conventional filtered back projection, adaptive statistical iterative reconstruction, and MBIR algorithms. All the results were compared to the performance of the standard-dose protocol.

RESULTS

Mean absorbed organ and effective doses were reduced by approximately 95% with the LD protocol (100-kVp tube voltage and a fixed 10-mA tube current) compared to the standard-dose protocol (120-kVp tube voltage and tube current modulation) while yielding an acceptable image quality for diagnosing round-shaped lung metastases. The effective doses ranged from 0.16 to 2.83 mSv in the studied protocols. The image noise, contrast, and resolution were maintained or improved when comparing the image quality of LD protocols using MBIR to the performance of the standard-dose chest CT protocol using filtered back projection. The small round-shaped lung metastases were delineated at levels comparable to the used protocols.

CONCLUSIONS

Radiation exposure in patients can be reduced significantly by using LD chest CT protocols and MBIR algorithm while maintaining image quality for detecting round-shaped lung metastases.

Spatial-temporal variability of radiomic features and its effect on the classification of lung cancer histology

The purpose of this research was to study the sensitivity of Computed Tomography (CT) radiomic features to motion blurring and signal-to-noise ratios (SNR), and investigate its downstream effect regarding the classification of non-small cell lung cancer (NSCLC) histology. Forty-three radiomic features were considered and classified into one of four categories: Morphological, Intensity, Fine Texture, and Coarse Texture. First, a series of simulations were used to study feature-sensitivity to changes in spatial-temporal resolution. A dynamic digital phantom was used to generate images with different breathing amplitudes and SNR, from which features were extracted and characterized relative to initial simulation conditions. Stage I NSCLC patients were then retrospectively identified, from which three different acquisition-specific feature-spaces were generated based on free-breathing (FB), average-intensity-projection (AIP), and end-of-exhalation (EOE) CT images. These feature-spaces were derived to cover a wide range of spatial-temporal tradeoff. Normalized percent differences and concordance correlation coefficients (CCC) were used to assess the variability between the 3D and 4D acquisition techniques. Subsequently, three corresponding acquisition-specific logistic regression models were developed to classify lung tumor histology. Classification performance was compared between the different data-dependent models. Simulation results demonstrated strong linear dependences (p  >  0.95) between respiratory motion and morphological features, as well as between SNR and texture features. The feature Short Run Emphasis was found to be particularly stable to both motion blurring and changes in SNR. When comparing FB-to-EOE, 37% of features demonstrated high CCC agreement (CCC  >  0.8), compared to only 30% for FB-to-AIP. In classifying tumor histology, EoE images achieved an average AUC, Accuracy, Sensitivity, and Specificity of [Formula: see text], respectively. FB images achieved respective values of [Formula: see text], and AIP images achieved respective values of [Formula: see text]. Several radiomic features have been identified as being relatively robust to spatial-temporal variations based on both simulation data and patient data. In general, features that were sensitive to motion blurring were not necessarily the same features that were sensitive to changes in SNR. Our modeling results suggest that the EoE phase of a 4DCT acquisition may provide useful radiomic information, particularly for features that are highly sensitive to respiratory motion.

Development of an organ-specific insert phantom generated using a 3D printer for investigations of cardiac computed tomography protocols

INTRODUCTION

An ideal organ-specific insert phantom should be able to simulate the anatomical features with appropriate appearances in the resultant computed tomography (CT) images. This study investigated a 3D printing technology to develop a novel and cost-effective cardiac insert phantom derived from volumetric CT image datasets of anthropomorphic chest phantom.

METHODS

Cardiac insert volumes were segmented from CT image datasets, derived from an anthropomorphic chest phantom of Lungman N-01 (Kyoto Kagaku, Japan). These segmented datasets were converted to a virtual 3D-isosurface of heart-shaped shell, while two other removable inserts were included using computer-aided design (CAD) software program. This newly designed cardiac insert phantom was later printed by using a fused deposition modelling (FDM) process via a Creatbot DM Plus 3D printer. Then, several selected filling materials, such as contrast media, oil, water and jelly, were loaded into designated spaces in the 3D-printed phantom. The 3D-printed cardiac insert phantom was positioned within the anthropomorphic chest phantom and 30 repeated CT acquisitions performed using a multi-detector scanner at 120-kVp tube potential. Attenuation (Hounsfield Unit, HU) values were measured and compared to the image datasets of real-patient and Catphan® 500 phantom.

RESULTS

The output of the 3D-printed cardiac insert phantom was a solid acrylic plastic material, which was strong, light in weight and cost-effective. HU values of the filling materials were comparable to the image datasets of real-patient and Catphan® 500 phantom.

CONCLUSIONS

A novel and cost-effective cardiac insert phantom for anthropomorphic chest phantom was developed using volumetric CT image datasets with a 3D printer. Hence, this suggested the printing methodology could be applied to generate other phantoms for CT imaging studies.

Paper-based 3D printing of anthropomorphic CT phantoms: Feasibility of two construction techniques

OBJECTIVES

To develop and evaluate methods for assembling radiopaque printed paper sheets to realistic patient phantoms for CT dose and image quality testing.

METHODS

CT images of two patients were radiopaque printed with aqueous potassium iodide solution (0.6 g/ml) on paper. Two methods were developed for assembling the paper sheets to head and neck phantoms. (1) Printed sheets were fed to a paper-based 3D printer along with corresponding 3D printable STL files. (2) Paper stacks of 5-mm thickness were glued with toner, cut to the patient shape and assembled to a phantom. In a sample application study, both phantoms were examined with five different tube current settings. Images were reconstructed using filtered-back projection (FBP) and iterative reconstruction (AIDR 3D) with three strength levels. Dose length product (DLP), signal-to-noise ratios (SNR) and contrast-to-noise ratios (CNRs) were analysed. Data were analysed using 2-way analysis of variance (ANOVA).

RESULTS

Both methods achieved anthropomorphic phantoms with detailed patient anatomy. The 3D printer yielded a precise reproduction of the external patient shape, but caused visible glue artefacts. Gluing with toner avoided these artefacts and yielded more flexibility with regard to phantom size. In the sample application study, non-inferior SNR and CNR and up to 83.7% lower DLP were achieved on the phantoms with AIDR 3D compared with FBP.

CONCLUSIONS

Two methods for assembling radiopaque printed paper sheets to phantoms of individual patients are presented. The sample application demonstrates potential for simulation of patient imaging and systematic CT dose and image quality assessment.

KEY POINTS

• Two methods were developed to create realistic CT phantoms of individual patients from radiopaque printed paper sheets. • Analysis of five tube current and four reconstruction settings on two radiopaque 3D printed patient phantoms yielded non-inferior SNR and CNR and up to 83.7% lower dose with iterative reconstruction in comparison with filtered back projection. • Radiopaque 3D printed phantoms can simulate patients and allow systematic analysis of CT dose and image quality parameters.

Modeling Lung Architecture in the XCAT Series of Phantoms: Physiologically Based Airways, Arteries and Veins

The purpose of this paper was to extend the extended cardiac-torso (XCAT) series of computational phantoms to include a detailed lung architecture including airways and pulmonary vasculature. Eleven XCAT phantoms of varying anatomy were used in this paper. The lung lobes and initial branches of the airways, pulmonary arteries, and veins were previously defined in each XCAT model. These models were extended from the initial branches of the airways and vessels to the level of terminal branches using an anatomically-based volume-filling branching algorithm. This algorithm grew the airway and vasculature branches separately and iteratively without intersecting each other using cylindrical models with diameters estimated by order-based anatomical measurements. Geometrical features of the extended branches were compared with the literature anatomy values to quantitatively evaluate the models. These features include branching angle, length to diameter ratio, daughter to parent diameter ratio, asymmetrical branching pattern, diameter, and length ratios. The XCAT phantoms were then used to simulate CT images to qualitatively compare them with the original phantom images. The proposed growth model produced 46369 ± 12521 airways, 44737 ± 11773 arteries, and 39819 ± 9988 veins to the XCAT phantoms. Furthermore, the growth model was shown to produce asymmetrical airway, artery, and vein networks with geometrical attributes close to morphometry and model based studies. The simulated CT images of the phantoms were judged to be more realistic, including more airways and pulmonary vessels compared with the original phantoms. Future work will seek to add a heterogeneous parenchymal background into the XCAT lungs to make the phantoms even more representative of human anatomy, paving the way towards the use of XCAT models as a tool to virtually evaluate the current and emerging medical imaging technologies.

Characterization and validation of the thorax phantom Lungman for dose assessment in chest radiography optimization studies

This work concerns the validation of the Kyoto-Kagaku thorax anthropomorphic phantom Lungman for use in chest radiography optimization. The equivalence in terms of polymethyl methacrylate (PMMA) was established for the lung and mediastinum regions of the phantom. Patient chest examination data acquired under automatic exposure control were collated over a 2-year period for a standard x-ray room. Parameters surveyed included exposure index, air kerma area product, and exposure time, which were compared with Lungman values. Finally, a voxel model was developed by segmenting computed tomography images of the phantom and implemented in PENELOPE/penEasy Monte Carlo code to compare phantom tissue-equivalent materials with materials from ICRP Publication 89 in terms of organ dose. PMMA equivalence varied depending on tube voltage, from 9.5 to 10.0 cm and from 13.5 to 13.7 cm, for the lungs and mediastinum regions, respectively. For the survey, close agreement was found between the phantom and the patients' median values (deviations lay between 8% and 14%). Differences in lung doses, an important organ for optimization in chest radiography, were below 13% when comparing the use of phantom tissue-equivalent materials versus ICRP materials. The study confirms the value of the Lungman for chest optimization studies.

A novel method for quantification of beam's-eye-view tumor tracking performance

PURPOSE

In-treatment imaging using an electronic portal imaging device (EPID) can be used to confirm patient and tumor positioning. Real-time tumor tracking performance using current digital megavolt (MV) imagers is hindered by poor image quality. Novel EPID designs may help to improve quantum noise response, while also preserving the high spatial resolution of the current clinical detector. Recently investigated EPID design improvements include but are not limited to multi-layer imager (MLI) architecture, thick crystalline and amorphous scintillators, and phosphor pixilation and focusing. The goal of the present study was to provide a method of quantitating improvement in tracking performance as well as to reveal the physical underpinnings of detector design that impact tracking quality. The study employs a generalizable ideal observer methodology for the quantification of tumor tracking performance. The analysis is applied to study both the effect of increasing scintillator thickness on a standard, single-layer imager (SLI) design as well as the effect of MLI architecture on tracking performance.

METHODS

The present study uses the ideal observer signal-to-noise ratio (d') as a surrogate for tracking performance. We employ functions which model clinically relevant tasks and generalized frequency-domain imaging metrics to connect image quality with tumor tracking. A detection task for relevant Cartesian shapes (i.e., spheres and cylinders) was used to quantitate trackability of cases employing fiducial markers. Automated lung tumor tracking algorithms often leverage the differences in benign and malignant lung tissue textures. These types of algorithms (e.g., soft-tissue localization - STiL) were simulated by designing a discrimination task, which quantifies the differentiation of tissue textures, measured experimentally and fit as a power-law in trend (with exponent β) using a cohort of MV images of patient lungs. The modeled MTF and NPS were used to investigate the effect of scintillator thickness and MLI architecture on tumor tracking performance.

RESULTS

Quantification of MV images of lung tissue as an inverse power-law with respect to frequency yields exponent values of β = 3.11 and 3.29 for benign and malignant tissues, respectively. Tracking performance with and without fiducials was found to be generally limited by quantum noise, a factor dominated by quantum detective efficiency (QDE). For generic SLI construction, increasing the scintillator thickness (gadolinium oxysulfide - GOS) from a standard 290 μm to 1720 μm reduces noise to about 10%. However, 81% of this reduction is appreciated between 290 and 1000 μm. In comparing MLI and SLI detectors of equivalent individual GOS layer thickness, the improvement in noise is equal to the number of layers in the detector (i.e., 4) with almost no difference in MTF. Further, improvement in tracking performance was slightly less than the square-root of the reduction in noise, approximately 84-90%. In comparing an MLI detector with an SLI with a GOS scintillator of equivalent total thickness, improvement in object detectability is approximately 34-39%.

CONCLUSIONS

We have presented a novel method for quantification of tumor tracking quality and have applied this model to evaluate the performance of SLI and MLI EPID designs. We showed that improved tracking quality is primarily limited by improvements in NPS. When compared to very thick scintillator SLI, employing MLI architecture exhibits the same gains in QDE, but by mitigating the effect of optical Swank noise, results in more dramatic improvements in tracking performance.