Task-based image quality assessment of an intraoperative CBCT for spine surgery compared with conventional CT

PURPOSE

To analyze the image quality of a novel, state-of-the art platform for CBCT image-guided spine surgery, focusing particularly on the dose-effectiveness compared with conventional CT (the gold standard for postoperative assessment).

METHODS

The ClarifEye platform (Philips Healthcare) with integrated augmented-reality surgical navigation, has been compared with a GE Revolution CT (GE Healthcare). The 3D spatial resolution (TTF) and noise (NPS) were evaluated considering relevant feature contrasts (200-900 HU) and background noise for differently sized patients (200-300 mm water-equivalent diameter). These measures were used to determine the noise equivalent quanta (NEQ) and observer model detectability.

RESULTS

The CBCT system exhibited a linear response with 50% TTF at 5.7 cycles/cm (10% TTF at 9.2 cycles/cm), and the axial noise power peaking at about 3.6 cycles/cm (average frequency of 4.1 cycles/cm). The noise magnitude and texture differed markedly compared to iteratively reconstructed CT images (GE ASiR-V). The CBCT system had 26% lower detectability for a high-frequency task (related to edge detection) compared with CT images reconstructed using the Bone kernel combined with ASiR-V 50%. Likewise, it had 18% lower detectability for low- and mid-frequency tasks compared with CT images reconstructed using the Standard kernel. This difference translates to 50%-80% higher CBCT imaging doses required to match the CT image quality.

CONCLUSIONS

The ClarifEye platform demonstrates intraoperative CBCT-imaging capabilities that under certain circumstances are comparable with conventional CT. However, due to limited dose-effectiveness, a trade-off between timeliness and radiation exposure must be considered if end-of-procedure CBCT is to replace postoperative CT.

The impact of automatic tube current modulation related settings of a modern GE CT scanner on image quality and patient dose; details do matter

PURPOSE

To investigate the operation principles of the automatic tube current modulation (ATCM) of a modern GE healthcare CT scanner, and the impact of related settings on image quality and patient dose.

MATERIAL & METHODS

A dedicated phantom (Mercury 4.0) was scanned using two of the most frequently used clinical scanning protocols (chest and abdomen-pelvis). The preset protocol settings were used as starting points (reference conditions). Scan direction, scan mode (helical vs. axial), total beam width, tube potential (kVp), and ATCM settings were then modified individually to understand their impact on radiation dose and image quality. Regarding the ATCM settings, the SmartmA minimum and maximum mA limits, and the noise index (NI) values were varied. As surrogates of patient dose, the CTDIvol and DLP values of each scan were used. As surrogates of image quality were used the image noise and the detectability index (d') of five different materials (air, solid water, polystyrene, iodine, and bone) embedded in the Mercury phantom calculated with the ImQuest software.

RESULTS

The scanning direction did not have any effect on ATCM curves, unlike what has been observed in CT scanners from other manufacturers. Total beam width does matter, however, the SmartmA limit settings and kVp selection had the greatest impact on image quality and dose. It was seen that improper minimum mA limit settings practically invalidated the ATCM operation. In contrast, when full modulation was allowed without restrictions, noise standard deviation, and detectability index became much more consistent across the wide range of phantom diameters. For lower kVp settings an impressive dose reduction was observed that requires further investigation.

CONCLUSION

SmartmA is a tool that if not properly used may increase the patient doses considerably. Therefore, its settings should be carefully adjusted for each preset different clinical protocol.

Investigation of a potential upstream harmonization based on image appearance matching to improve radiomics features robustness: a phantom study

Objective. Radiomics is a promising valuable analysis tool consisting in extracting quantitative information from medical images. However, the extracted radiomics features are too sensitive to variations in used image acquisition and reconstruction parameters. This limited robustness hinders the generalizable validity of radiomics-assisted models. Our aim is to investigate a possible harmonization strategy based on matching image quality to improve feature robustness.Approach.We acquired CT scans of a phantom with two scanners across different dose levels and percentages of Iterative Reconstruction algorithms. The detectability index was used as a comprehensive task-based image quality metric. A statistical analysis based on the Intraclass Correlation Coefficient was performed to determine if matching image quality/appearance could enhance the robustness of radiomics features extracted from the phantom images. Additionally, an Artificial Neural Network was trained on these features to automatically classify the scanner used for image acquisition.Main results.We found that the ICC of the features across protocols providing a similar detectability index improves with respect to the ICC of the features across protocols providing a different detectability index. This improvement was particularly noticeable in features relevant for distinguishing between scanners.Significance.This preliminary study demonstrates that a harmonization based on image quality/appearance matching could improve radiomics features robustness and heterogeneous protocols can be used to obtain a similar image appearance in terms of the detectability index. Thus protocols with a lower dose level could be selected to reduce the amount of radiation dose delivered to the patient and simultaneously obtain a more robust quantitative analysis.

Technology Characterization Through Diverse Evaluation Methodologies: Application to Thoracic Imaging in Photon-Counting Computed Tomography

OBJECTIVE

Different methods can be used to condition imaging systems for clinical use. The purpose of this study was to assess how these methods complement one another in evaluating a system for clinical integration of an emerging technology, photon-counting computed tomography (PCCT), for thoracic imaging.

METHODS

Four methods were used to assess a clinical PCCT system (NAEOTOM Alpha; Siemens Healthineers, Forchheim, Germany) across 3 reconstruction kernels (Br40f, Br48f, and Br56f). First, a phantom evaluation was performed using a computed tomography quality control phantom to characterize noise magnitude, spatial resolution, and detectability. Second, clinical images acquired using conventional and PCCT systems were used for a multi-institutional reader study where readers from 2 institutions were asked to rank their preference of images. Third, the clinical images were assessed in terms of in vivo image quality characterization of global noise index and detectability. Fourth, a virtual imaging trial was conducted using a validated simulation platform (DukeSim) that models PCCT and a virtual patient model (XCAT) with embedded lung lesions imaged under differing conditions of respiratory phase and positional displacement. Using known ground truth of the patient model, images were evaluated for quantitative biomarkers of lung intensity histograms and lesion morphology metrics.

RESULTS

For the physical phantom study, the Br56f kernel was shown to have the highest resolution despite having the highest noise and lowest detectability. Readers across both institutions preferred the Br56f kernel (71% first rank) with a high interclass correlation (0.990). In vivo assessments found superior detectability for PCCT compared with conventional computed tomography but higher noise and reduced detectability with increased kernel sharpness. For the virtual imaging trial, Br40f was shown to have the best performance for histogram measures, whereas Br56f was shown to have the most precise and accurate morphology metrics.

CONCLUSION

The 4 evaluation methods each have their strengths and limitations and bring complementary insight to the evaluation of PCCT. Although no method offers a complete answer, concordant findings between methods offer affirmatory confidence in a decision, whereas discordant ones offer insight for added perspective. Aggregating our findings, we concluded the Br56f kernel best for high-resolution tasks and Br40f for contrast-dependent tasks.

The effect of a pre-reconstruction process in a filtered back projection reconstruction on an image quality of a low tube voltage computed tomography

This study aims to evaluate the effect of pre-reconstruction process for low tube voltage computed tomography (CT) on image quality of filtered back projection (FBP) reconstruction. Small and large quality assurance water phantoms (19 and 33 cm diameter) were scanned on a third-generation dual-source CT with 70 kVp and 120 kVp at various dose levels. Image quality was assessed in terms of the noise power spectrum (NPS) and task-based transfer function (TTF). NPSs and TTFs in the small phantom were comparable between 70 and 120 kVp protocols. In the large phantom, the curves of the NPS changed and the TTF decreased even at the high-dose levels for 70 kVp protocol compared to 120 kVp protocol. Our results indicated that the pre-reconstruction process is performed in low tube voltage CT for large objects even for the FBP reconstruction and has an effect on the image quality.

Usefulness of Ag Additional Filter on Image Quality and Radiation Dose for Low-Dose Chest Computed Tomography

OBJECTIVE

This study aimed to evaluate the effect of a silver (Ag) additional filter on dose characteristics and image quality in low-dose chest computed tomography (CT).

METHODS

A dose evaluation phantom, physical evaluation phantom, and chest phantom were scanned with and without an Ag additional filter. The doses were adjusted so that the displayed the volume CT dose indexes (CTDI vol ) were from 0.3 to 1.6 mGy. For dose characteristics, the spectrum of photon energies and the measured CTDI vol were calculated for each scanning condition. For task-based image quality analysis, task transfer function, noise power spectrum, and system performance were evaluated. Streak artifacts, image noise, and contrast-to-noise ratio were quantified using a chest phantom.

RESULTS

With the Ag additional filter, mean energy was 22% higher and the CTDI vol was approximately 30% lower than those without the Ag additional filter. The task transfer function and noise power spectrum with the Ag additional filter were lower than those without the Ag additional filter. The system performance with the Ag additional filter was similar to that without the Ag additional filter. The Ag additional filter reduced streak artifact near the lung apex and image noise in the lung fields. The contrast-to-noise ratio was slightly higher with the Ag additional filter than that without the Ag additional filter.

CONCLUSIONS

The output dose and spatial resolution with the Ag additional filter were lower than those without the Ag additional filter. However, this filter helped reduce the radiation dose, image noise, and streak artifacts, particularly when scanning at ultralow doses.

A systematic task-based image quality assessment of photon-counting and energy integrating CT as a function of reconstruction kernel and phantom size

BACKGROUND

Image quality of photon-counting and energy integrating CT scanners changes with object size, dose to the object, and kernel selection.

PURPOSE

To comprehensively compare task-generic image quality of photon-counting CT (PCCT) and energy integrating CT (EICT) systems as a function of phantom size, dose, and reconstruction kernel.

METHODS

A size-variant phantom (Mercury Phantom 3.0) was used to characterize the image quality of PCCT and EICT systems as a function of object size. The phantom contained five cylinders attached by slanted tapered sections. Each cylinder contained two sections: one uniform for noise, and the other with inserts for spatial resolution and contrast measurements. The phantom was scanned on Siemens' SOMATOM Force and NAEOTOM Alpha at 1.18 and 7.51 mGy without tube current modulation. CTDIvol was matched across two systems by setting the required tube currents, else, all other acquisition settings were fixed. CT sinograms were reconstructed using FBP and iterative (ADMIRE2 - Force; QIR2 - Alpha) algorithms with Body regular (Br) kernels. Noise Power Spectrum (NPS), Task Transfer Function (TTF), contrast-to-noise ratio (CNR), and detectability index (d') for a task of identifying 2-mm disk were evaluated based on AAPM TG-233 metrology using imQuest, an open-source software package. Averaged noise frequency (fav ) and 50% cut-off frequency for TTF (f50 ) were used as scalar metrics to quantify noise texture and spatial resolution, respectively. The difference between image quality metrics' measurements was calculated as IQPCCT - IQEICT .

RESULTS

From Br40 (soft) to Br64 (sharp), f50 for air insert increased from 0.35 mm-1  ± 0.04 (standard deviation) to 0.76 mm-1  ± 0.17, 0.34 mm-1  ± 0.04 to 0.77 mm-1  ± 0.17, 0.37 mm-1  ± 0.02 to 0.95 mm-1  ± 0.17 for PCCT-T3D-QIR2, PCCT-70keV-QIR2, and EICT-ADMIRE2, respectively, when averaged over all sizes and dose levels. Similarly, from Br40 to Br64, noise magnitude increased from 10.86 HU ± 7.12 to 38.61 HU ± 18.84, 10.94 HU ± 7.08 to 38.82 HU ± 18.70, 13.74 HU ± 11.02 to 52.11 HU ± 26.22 for PCCT-T3D-QIR2, PCCT-70keV-QIR2, and EICT-ADMIRE2, respectively. The difference in fav was consistent across all sizes and dose levels. PCCT-70keV-VMI showed better consistency in contrast measurements for iodine and bone inserts than PCCT-T3D and EICT; however, PCCT-T3D had higher contrast for both inserts. From Br40 to Br64, considering all sizes and dose levels, CNR for iodine insert decreased from 52.30 ± 46.44 to 12.18 ± 10.07, 40.42 ± 33.42 to 9.48 ± 7.16, 39.94 ± 37.60 to 7.84 ± 6.67 for PCCT-T3D-QIR2, PCCT-70keV-QIR2, and EICT-ADMIRE2, respectively.

CONCLUSIONS

Both PCCT image types, T3D and 70-keV-VMI exhibited similar or better noise, contrast, CNR than EICT when comparing kernels with similar names. For 512 × 512 matrix, PCCT's sharp kernels had lower resolution than EICT's sharp kernels. For all image quality metrics, except extreme low, every dose condition had a similar set of IQ-matching kernels. It suggests that considering patient size and dose level to determine IQ-matching kernel pairs across PCCT and EICT systems may not be imperative while translating protocols, except when the signal to the detector is extremely low.

Sparsier2Sparse: Self-supervised convolutional neural network-based streak artifacts reduction in sparse-view CT images

BACKGROUND

Sparse-view computed tomography (CT) has attracted a lot of attention for reducing both scanning time and radiation dose. However, sparsely-sampled projection data generate severe streak artifacts in the reconstructed images. In recent decades, many sparse-view CT reconstruction techniques based on fully-supervised learning have been proposed and have shown promising results. However, it is not feasible to acquire pairs of full-view and sparse-view CT images in real clinical practice.

PURPOSE

In this study, we propose a novel self-supervised convolutional neural network (CNN) method to reduce streak artifacts in sparse-view CT images.

METHODS

We generate the training dataset using only sparse-view CT data and train CNN based on self-supervised learning. Since the streak artifacts can be estimated using prior images under the same CT geometry system, we acquire prior images by iteratively applying the trained network to given sparse-view CT images. We then subtract the estimated steak artifacts from given sparse-view CT images to produce the final results.

RESULTS

We validated the imaging performance of the proposed method using extended cardiac-torso (XCAT) and the 2016 AAPM Low-Dose CT Grand Challenge dataset from Mayo Clinic. From the results of visual inspection and modulation transfer function (MTF), the proposed method preserved the anatomical structures effectively and showed higher image resolution compared to the various streak artifacts reduction methods for all projection views.

CONCLUSIONS

We propose a new framework for streak artifacts reduction when only the sparse-view CT data are given. Although we do not use any information of full-view CT data for CNN training, the proposed method achieved the highest performance in preserving fine details. By overcoming the limitation of dataset requirements on fully-supervised-based methods, we expect that our framework can be utilized in the medical imaging field.

Exploration of the pulse pileup effects in a clinical CdTe-based photon-counting computed tomography

BACKGROUND

High tube current generates a high flux of x-rays to photon counting detectors (PCDs) that can potentially result in the piling up of pulses formed by concurrent photons, which can cause count loss and energy resolution degradation.

PURPOSE

To evaluate the performance of clinical photon-counting CT (PCCT) systems in high flux, potentially influenced by pulse pileup effects, in terms of task-generic image quality metrics.

METHODS

A clinical phantom was scanned on a commercial PCCT scanner (NAEOTOM Alpha, Siemens) at 120 kV under fourteen different tube current levels (40-1000 mA) with a rotation time of 0.25 s and a pitch of 1. The dose levels corresponded to CTDIvol (32 cm phantom) of 0.79-19.8 mGy. CT sinograms were reconstructed using QIR-off mode (noniterative reconstruction algorithm), Br44 kernel, and a voxel size of0.4102 × 0.4102 × 3 mm 3 $0.4102 \times 0.4102 \times 3{\mathrm{\ mm}}^3$ . imQuest, an open-source MATLAB-based software package was used to calculate noise power spectrum (NPS), task transfer function (TTF), contrast-to-noise ratio (CNR), and CT number according to AAPM Task Group 233 metrology.

RESULTS

The 50% cut-off frequency of TTF (f50 ) remained mostly constant across all higher tube currents for all inserts, namely polyethylene, bone, air, and acrylic. Using the lowest two data points (40 and 80 mA), the expected relationship between noise magnitude and tube current was determined to be noise∝ $ \propto \ $ mA-0.47 . The measured noise magnitude were up to 11.1% higher than the expected value at the highest tube current. The average frequency of NPS (fav ) decreased from 0.32 to 0.29 mm-1 as tube current increased from 40 to 1000 mA. No considerable effects were observed in CT number measurement of any insert; however, CT numbers for air and bone changed almost monotonically as tube current increased. Absolute CNR increased monotonically for all inserts; however, the difference between measured and expected CNRs were approximately -6% to 12% across all tube currents.

CONCLUSIONS

Increasing tube currents did not affect the spatial resolution, but slightly affected the CT number and noise measurements of the clinical PCCT system. However, the effects were only considerable at clinically irrelevant tube currents used on a small 20-cm phantom. In general clinical practices, automatic exposure control techniques are used to decrease the variation of flux on the detector, which alleviates the chances of detector saturation due to high count rates. The observed effects could be due to pulse pileup, signal-dependent filtration of the system, or nonlinearities in the reconstruction algorithm. In conclusion, either the deadtime of the detector used in the photon-counting CT system is shorter such that count losses due to pulse pileup are negligible, or pulse pileup has inconsiderable effects on the image quality of clinical photon-counting CT systems in routine clinical practice due to possible corrections applied on the system.

Noise reduction performance of a deep learning-based reconstruction in brain computed tomography images acquired with organ-based tube current modulation

We aimed to evaluate the image quality of brain computed tomography (CT) images reconstructed using deep learning-based reconstruction (DLR) in organ-based tube current modulation (OB-TCM) acquisition. An anthropomorphic head phantom and a cylindrical low-contrast phantom were scanned at the standard dose level for adult brain CT in axial volume acquisition without OB-TCM. Moreover, image acquisition with OB-TCM was performed. The radiation dose on the eye lens was measured using a scintillation fibre-optic dosimeter placed on the anthropomorphic phantom's eye surface. The task transfer function (TTF), contrast-to-noise ratio (CNR), and low-contrast object specific CNR obtained from low-contrast phantom images reconstructed with filtered back projection (FBP), hybrid iterative reconstruction (HIR), and two types of DLR (DLRCTA and DLRLCD) were compared. In result, OB-TCM achieved a 32.5% dose reduction in the eye lens. Although HIR, DLRCTA, and DLRLCD showed lower TTF than FBP, the difference in TTF at the highest contributing spatial frequency corresponding to the contrast rod diameter was < 10%. Despite the OB-TCM acquisition, DLRCTA and DLRLCD achieved significantly lower noise and a higher CNR than FBP without OB-TCM (p < 0.05). However, low-contrast object specific CNR was equivalent among all reconstruction methods for the objective diameter of 5 mm and slightly improved in DLRLCD for the objective diameter of 7 mm. DLR with OB-TCM acquisition enabled dose reduction for the eye lens and high CNR image appearance, whereas the low contrast detectability evaluated by low-contrast object specific CNR did not always improve.

Evaluation of Image Quality and Detectability of Deep Learning Image Reconstruction (DLIR) Algorithm in Single- and Dual-energy CT

This study is aimed to evaluate effects of deep learning image reconstruction (DLIR) on image quality in single-energy CT (SECT) and dual-energy CT (DECT), in reference to adaptive statistical iterative reconstruction-V (ASIR-V). The Gammex 464 phantom was scanned in SECT and DECT modes at three dose levels (5, 10, and 20 mGy). Raw data were reconstructed using six algorithms: filtered back-projection (FBP), ASIR-V at 40% (AV-40) and 100% (AV-100) strength, and DLIR at low (DLIR-L), medium (DLIR-M), and high strength (DLIR-H), to generate SECT 120kVp images and DECT 120kVp-like images. Objective image quality metrics were computed, including noise power spectrum (NPS), task transfer function (TTF), and detectability index (d'). Subjective image quality evaluation, including image noise, texture, sharpness, overall quality, and low- and high-contrast detectability, was performed by six readers. DLIR-H reduced overall noise magnitudes from FBP by 55.2% in a more balanced way of low and high frequency ranges comparing to AV-40, and improved the TTF values at 50% for acrylic inserts by average percentages of 18.32%. Comparing to SECT 20 mGy AV-40 images, the DECT 10 mGy DLIR-H images showed 20.90% and 7.75% improvement in d' for the small-object high-contrast and large-object low-contrast tasks, respectively. Subjective evaluation showed higher image quality and better detectability. At 50% of the radiation dose level, DECT with DLIR-H yields a gain in objective detectability index compared to full-dose AV-40 SECT images used in daily practice.

A methodology for incorporating a photon-counting CT system into routine clinical use

Photon-counting computed tomography (PCCT) systems are increasingly available in the U.S. following Food and Drug Administration (FDA) approval of the first clinical PCCT system in Fall 2021. Consequently, there will be a need to incorporate PCCTs into existing fleets of traditional CT systems. The commissioning process of a PCCT was devised by evaluating the degree of agreement between the performance of the PCCT and that of established clinical CT systems. A PCCT system (Siemens NAEOTOM Alpha) was evaluated using the American College of Radiology(ACR) CT phantom (Gammex 464). The phantom was scanned on the system and on a 3rd Generation EID CT system (Siemens Force) at three clinical dose levels. Images were reconstructed across the range of available reconstruction kernels and Iterative Reconstruction (IR) strengths. Two image quality metrics-spatial resolution and noise texture-were calculated using AAPM TG233 software (imQuest), as well as a dose metric to achieve target image noise magnitude of 10 HU. For each pair of EID-PCCT kernel/IR strengths, the difference in metrics were calculated, weighted, and multiplied over all metrics to determine the concordance between systems. IR performance was characterized by comparing relative noise texture and reference dose as a function of IR strength for each system. In general, as kernel "sharpness" increased for each system, spatial resolution, noise spatial frequency, and reference dose increased. For a given kernel, EID reconstruction showed higher spatial resolution compared to PCCT in standard resolution mode. PCCT implementation of IR better preserved noise texture across all strengths compared to the EID, demonstrated by respective 20 and 7% shifts in noise texture from IR "Off" to IR "Max." Overall, the closest match for a given EID reconstruction kernel/IR strength was identified as a PCCT kernel with "sharpness" increased by 1 step and IR strength increased by 1-2 steps. Substantial dose reduction potential of up to 70% was found when targeting a constant noise magnitude.

High Spatial-Resolution Skull Base Imaging With Photon-Counting Computed Tomography and Energy-Integrating Computed Tomography: A Comparative Phantom Study

ABSTRACT

Photon-counting computed tomography (PCCT) offers better high-resolution and noise performance than energy integrating detector (EID) CT. In this work, we compared both technologies for imaging of the temporal bone and skull base. A clinical PCCT system and 3 clinical EID CT scanners were used to image the American College of Radiology image quality phantom using a clinical imaging protocol with matched CTDI vol (CT dose index-volume) of 25 mGy. Images were used to characterize the image quality of each system across a series of high-resolution reconstruction options. Noise was calculated from the noise power spectrum, whereas resolution was quantified with a bone insert by calculating a task transfer function. Images of an anthropomorphic skull phantom and 2 patient cases were examined for visualization of small anatomical structures. Across measured conditions, PCCT had a comparable or smaller average noise magnitude (120 Hounsfield units [HU]) to the EID systems (144-326 HU). Photon-counting CT also had comparable resolution (task transfer function f25 : 1.60 mm -1 ) to the EID systems (1.34-1.77 mm -1 ). Imaging results supported quantitative findings as PCCT more clearly showed the 12-lp/cm bars from the fourth section of the American College of Radiology phantom and better represented the vestibular aqueduct and oval and round windows when compared with the EID scanners. A clinical PCCT system was able to image the temporal bone and skull base with improved spatial resolution and lower noise than clinical EID CT systems at matched dose.

Noise power spectrum properties of deep learning-based reconstruction and iterative reconstruction algorithms: Phantom and clinical study

PURPOSE

To compare the noise power spectrum (NPS) properties and perform a qualitative analysis of hybrid iterative reconstruction (IR), model-based IR (MBIR), and deep learning-based reconstruction (DLR) at a similar noise level in clinical study and compare these outcomes with those in phantom study.

METHODS

A Catphan phantom with an external body ring was used in the phantom study. In the clinical study, computed tomography (CT) examination data of 34 patients were reviewed. NPS was calculated from DLR, hybrid IR, and MBIR images. The noise magnitude ratio (NMR) and the central frequency ratio (CFR) were calculated from DLR, hybrid IR, and MBIR images relative to filtered back-projection images using NPS. Clinical images were independently reviewed by two radiologists.

RESULTS

In the phantom study, DLR with a mild level had a similar noise level as hybrid IR and MBIR with strong levels. In the clinical study, DLR with a mild level had a similar noise level as hybrid IR with standard and MBIR with strong levels. The NMR and CFR were 0.40 and 0.76 for DLR, 0.42 and 0.55 for hybrid IR, and 0.48 and 0.62 for MBIR. The visual inspection of the clinical DLR image was superior to that of the hybrid IR and MBIR images.

CONCLUSION

Deep learning-based reconstruction improves overall image quality with substantial noise reduction while maintaining image noise texture compared with the CT reconstruction techniques.

Is it possible to use low-dose deep learning reconstruction for the detection of liver metastases on CT routinely?

OBJECTIVES

To compare the image quality and hepatic metastasis detection of low-dose deep learning image reconstruction (DLIR) with full-dose filtered back projection (FBP)/iterative reconstruction (IR).

METHODS

A contrast-detail phantom consisting of low-contrast objects was scanned at five CT dose index levels (10, 6, 3, 2, and 1 mGy). A total of 154 participants with 305 hepatic lesions who underwent abdominal CT were enrolled in a prospective non-inferiority trial with a three-arm design based on phantom results. Data sets with full dosage (13.6 mGy) and low dosages (9.5, 6.8, or 4.1 mGy) were acquired from two consecutive portal venous acquisitions, respectively. All images were reconstructed with FBP (reference), IR (control), and DLIR (test). Eleven readers evaluated phantom data sets for object detectability using a two-alternative forced-choice approach. Non-inferiority analyses were performed to interpret the differences in image quality and metastasis detection of low-dose DLIR relative to full-dose FBP/IR.

RESULTS

The phantom experiment showed the dose reduction potential from DLIR was up to 57% based on the reference FBP dose index. Radiation decreases of 30% and 50% resulted in non-inferior image quality and hepatic metastasis detection with DLIR compared to full-dose FBP/IR. Radiation reduction of 70% by DLIR performed inferiorly in detecting small metastases (< 1 cm) compared to full-dose FBP (difference: -0.112; 95% confidence interval [CI]: -0.178 to 0.047) and full-dose IR (difference: -0.123; 95% CI: -0.182 to 0.053) (p < 0.001).

CONCLUSION

DLIR enables a 50% dose reduction for detecting low-contrast hepatic metastases while maintaining comparable image quality to full-dose FBP and IR.

KEY POINTS

• Non-inferiority study showed that deep learning image reconstruction (DLIR) can reduce the dose to oncological patients with low-contrast lesions without compromising the diagnostic information. • Radiation dose levels for DLIR can be reduced to 50% of full-dose FBP and IR for detecting low-contrast hepatic metastases, while maintaining comparable image quality. • The reduction of radiation by 70% by DLIR is clinically acceptable but insufficient for detecting small low-contrast hepatic metastases (< 1 cm).

Spatial resolution, noise properties, and detectability index of a deep learning reconstruction algorithm for dual-energy CT of the abdomen

BACKGROUND

Iterative reconstruction (IR) has increasingly replaced traditional reconstruction methods in computed tomography (CT). The next paradigm shift in image reconstruction is likely to come from artificial intelligence, with deep learning reconstruction (DLR) solutions already entering the clinic. An enduring disadvantage to IR has been a change in noise texture, which can affect diagnostic confidence. DLR has demonstrated the potential to overcome this issue and has recently become available for dual-energy CT.

PURPOSE

To evaluate the spatial resolution, noise properties, and detectability index of a commercially available DLR algorithm for dual-energy CT of the abdomen and compare it to single-energy (SE) CT.

METHODS

An oval 25 cm x 35 cm custom-made phantom was scanned on a GE Revolution CT scanner (GE Healthcare, Waukesha, WI) at two dose levels (13 and 5 mGy) and two iodine concentrations (8 and 2 mg/mL), using three typical abdominal scan protocols: dual-energy (DE), SE 80 kV (SE-80 kV) and SE 120 kV (SE-120 kV). Reconstructions were performed with three strengths of IR (ASiR-V: AR0%, AR50%, AR100%) and three strengths of DLR (TrueFidelity: low, medium, high). The DE acquisitions were reconstructed as mono-energetic images between 40 and 80 keV. The noise power spectrum (NPS), task transfer function (TTF), and detectability index (d') were determined for the reconstructions following the recommendations of AAPM Task Group 233.

RESULTS

Noise magnitude reductions (relative to AR0%) for the SE protocols were on average (-29%, -21%) for (AR50%, TF-M), while for DE-70 keV were (-28%, -43%). There was less reduction in mean frequency (f_av ) for DLR than for IR, with similar results for SE and DE imaging. There was, however, a substantial change in the NPS shape when using DE with DLR, quantifiable by a marked reduction in the peak frequency (f_peak ) that was absent in SE mode. All protocols and reconstructions (including AR0%) exhibited slight to moderate shifts towards lower spatial frequencies at the lower dose (<12% in f_av ). Spatial resolution was consistently superior for DLR compared to IR for SE but not for DE. All protocols and reconstructions (including AR0%) showed decreased resolution with reduced dose and iodine concentration, with less decrease for DLR compared to IR. DLR displayed a higher d' than IR. The effect of energy was large: d' increased with lower keV, and SE-80 kV had higher d' than SE-120 kV. Using DE with DLR could provide higher d' than SE-80 kV at the higher dose but not at lower dose.

CONCLUSIONS

DE imaging with DLR maintained spatial resolution and reduced noise magnitude while displaying less change in noise texture than IR. The d' was also higher with DLR than IR, suggesting superiority in detectability of iodinated contrast. Despite these trends being consistent with those previously established for SE imaging, there were some noteworthy differences. For DE imaging there was no improvement in resolution compared to IR and a change in noise texture. DE imaging with low keV and DLR had superior detectability to SE DLR at the high dose but was not better than SE-80 kV at low dose.

Comparison of supervised-learning approaches for designing a channelized observer for image quality assessment in CT

BACKGROUND

The current paradigm for evaluating computed tomography (CT) system performance relies on a task-based approach. As the Hotelling observer (HO) provides an upper bound of observer performances in specific signal detection tasks, the literature advocates HO use for optimization purposes. However, computing the HO requires calculating the inverse of the image covariance matrix, which is often intractable in medical applications. As an alternative, dimensionality reduction has been extensively investigated to extract the task-relevant features from the raw images. This can be achieved by using channels, which yields the channelized-HO (CHO). The channels are only considered efficient when the channelized observer (CO) can approximate its unconstrained counterpart. Previous work has demonstrated that supervised learning-based methods can usually benefit CO design, either for generating efficient channels using partial least squares (PLS) or for replacing the Hotelling detector with machine-learning (ML) methods.

PURPOSE

Here we investigated the efficiency of a supervised ML-algorithm used to design a CO for predicting the performance of unconstrained HO. The ML-algorithm was applied either (1) in the estimator for dimensionality reduction, or (2) in the detector function.

METHODS

A channelized support vector machine (CSVM) was employed and compared against the CHO in terms of ability to predict HO performances. Both the CSVM and the CHO were estimated with channels derived from the singular value decomposition (SVD) of the system operator, principal component analysis (PCA), and PLS. The huge variety of regularization strategies proposed by CT system vendors for statistical image reconstruction (SIR) make the generalization capability of an observer a key point to consider upfront of implementation in clinical practice. To evaluate the generalization properties of the observers, we adopted a 2-step testing process: (1) achieved with the same regularization strategy (as in the training phase) and (2) performed using different reconstruction properties. We generated simulated- signal-known-exactly/background-known-exactly (SKE/BKE) tasks in which different noise structures were generated using Markov random field (MRF) regularizations using either a Green or a quadratic, function.

RESULTS

The CSVM outperformed the CHO for all types of channels and regularization strategies. Furthermore, even though both COs generalized well to images reconstructed with the same regularization strategy as the images considered in the training phase, the CHO failed to generalize to images reconstructed differently whereas the CSVM managed to successfully generalize. Lastly, the proposed CSVM observer used with PCA channels outperformed the CHO with PLS channels while using a smaller training data set.

CONCLUSION

These results argue for introducing the supervised-learning paradigm in the detector function rather than in the operator of the channels when designing a CO to provide an accurate estimate of HO performance. The CSVM with PCA channels proposed here could be used as a surrogate for HO in image quality assessment.

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.

A streak artifact reduction algorithm in sparse-view CT using a self-supervised neural representation

PURPOSE

Sparse-view computed tomography (CT) has been attracting attention for its reduced radiation dose and scanning time. However, analytical image reconstruction methods suffer from streak artifacts due to insufficient projection views. Recently, various deep learning-based methods have been developed to solve this ill-posed inverse problem. Despite their promising results, they are easily overfitted to the training data, showing limited generalizability to unseen systems and patients. In this work, we propose a novel streak artifact reduction algorithm that provides a system- and patient-specific solution.

METHODS

Motivated by the fact that streak artifacts are deterministic errors, we regenerate the same artifacts from a prior CT image under the same system geometry. This prior image need not be perfect but should contain patient-specific information and be consistent with full-view projection data for accurate regeneration of the artifacts. To this end, we use a coordinate-based neural representation that often causes image blur but can greatly suppress the streak artifacts while having multiview consistency. By employing techniques in neural radiance fields originally proposed for scene representations, the neural representation is optimized to the measured sparse-view projection data via self-supervised learning. Then, we subtract the regenerated artifacts from the analytically reconstructed original image to obtain the final corrected image.

RESULTS

To validate the proposed method, we used simulated data of extended cardiac-torso phantoms and the 2016 NIH-AAPM-Mayo Clinic Low-Dose CT Grand Challenge and experimental data of physical pediatric and head phantoms. The performance of the proposed method was compared with a total variation-based iterative reconstruction method, naive application of the neural representation, and a convolutional neural network-based method. In visual inspection, it was observed that the small anatomical features were best preserved by the proposed method. The proposed method also achieved the best scores in the visual information fidelity, modulation transfer function, and lung nodule segmentation.

CONCLUSIONS

The results on both simulated and experimental data suggest that the proposed method can effectively reduce the streak artifacts while preserving small anatomical structures that are easily blurred or replaced with misleading features by the existing methods. Since the proposed method does not require any additional training datasets, it would be useful in clinical practice where the large datasets cannot be collected.

[Evaluation of Noise Properties at Nonuniformity Area and Resolution Properties of CT Image Using Iterative Reconstruction Method]

PURPOSE

The purpose of this study was to investigate the resolution property and the noise properties at the nonuniformity area in the slice plane of iterative reconstruction (IR)-CT image.

METHODS

CT images of a phantom with nonuniformity areas including multiple contrast (medium, high, and ultra-high contrast) signals were acquired at various scan dose conditions and reconstructed with different iterative intensity levels. The noise properties at the nonuniformity and the uniformity areas were evaluated by measuring normalized noise power spectra (nNPSs) using subtracted images from sequential scanning data under the same scan conditions. We investigated the relationship between the noise properties and the resolution properties evaluated by measuring task transfer function (TTF) using multiple contrast signals before the subtraction.

RESULTS

There was a correlation between the nNPS at the nonuniformity area and the TTF because the nNPS values at high spatial frequency were increased with superior TTF (higher dose, mild iterative intensity, and higher contrast level). The differences in the high spatial frequency component of the nNPSs among each task were decreased with the inferior TTF tasks.

CONCLUSION

We conclude that the noise properties in the slice plane at the nonuniformity area and resolution properties of the IR-CT image were correlated similar to image quality properties of the linear imaging system due to the dependence of the nNPSs curve on the TTF value under each task.

[A Study of Longitudinal NPS Measurement in CT Images Based on the Central Cross-section Theorem]

PURPOSE

Various approaches in noise power spectrum (NPS) analysis are currently used for measuring a patient's longitudinal (z-direction) NPS from three-dimensional (3D) CT volume data. The purpose of this study was to clarify the relationship between those NPSs and 3D-NPS based on the central slice theorem.

METHODS

We defined the 3D-NPS(f_x, f_y, f_z) that was calculated by 3D Fourier transform (FT) from 3D noise data (3D-Noise(x, y, z), x-y scan plane). Here, f_x, f_y and f_z are spatial frequencies corresponding to the axes of x, y and z, respectively. Based on the central slice theorem, we described three relationships as follows. (1) The f_z-directional NPS calculated from the 3D-Noise(x=0, y=0, z) is equal to the profile obtained by projecting 3D-NPS(f_x, f_y, f_z) in f_x- and f_y-directions. (2) The f_z-directional NPS calculated from the profile obtained by projecting 3D-Noise(x=0, y, z) in the y-direction is equal to the profile at f_y=0 in the data obtained by projecting 3D-NPS(f_x, f_y, f_z) in the f_x-direction. (3) The f_z-directional NPS calculated from the profile obtained by projecting 3D-Noise(x, y, z) in x and y-directions is equal to the profile of 3D-NPS(f_x=0, f_y=0, f_z). To verify them, we compared the NPSs measured from actual 3D noise data that were obtained using a cylindrical water phantom.

RESULTS

In each relationship (1)-(3), the f_z-directional NPS matched the profile obtained from the 3D-NPS(f_x, f_y, f_z).

CONCLUSION

Based on the central slice theorem, we clarified the relationships between f_z-directional NPSs and 3D-NPS. We should understand them and then consider which method should be used for f_z-directional NPS measurement.

Optimization of a customized simultaneous algebraic reconstruction technique algorithm for phase-contrast breast computed tomography

Objective. To introduce the optimization of a customized GPU-based simultaneous algebraic reconstruction technique (cSART) in the field of phase-contrast breast computed tomography (bCT). The presented algorithm features a 3D bilateral regularization filter that can be tuned to yield optimal performance for clinical image visualization and tissues segmentation. Approach. Acquisitions of a dedicated test object and a breast specimen were performed at Elettra, the Italian synchrotron radiation (SR) facility (Trieste, Italy) using a large area CdTe single-photon counting detector. Tomographic images were obtained at 5 mGy of mean glandular dose, with a 32 keV monochromatic x-ray beam in the free-space propagation mode. Three independent algorithms parameters were optimized by using contrast-to-noise ratio (CNR), spatial resolution, and noise texture metrics. The results obtained with the cSART algorithm were compared with conventional SART and filtered back projection (FBP) reconstructions. Image segmentation was performed both with gray scale-based and supervised machine-learning approaches. Main results. Compared to conventional FBP reconstructions, results indicate that the proposed algorithm can yield images with a higher CNR (by 35% or more), retaining a high spatial resolution while preserving their textural properties. Alternatively, at the cost of an increased image ‘patchiness’, the cSART can be tuned to achieve a high-quality tissue segmentation, suggesting the possibility of performing an accurate glandularity estimation potentially of use in the realization of realistic 3D breast models starting from low radiation dose images. Significance. The study indicates that dedicated iterative reconstruction techniques could provide significant advantages in phase-contrast bCT imaging. The proposed algorithm offers great flexibility in terms of image reconstruction optimization, either toward diagnostic evaluation or image segmentation.

[A Study of 3D-NPS Analysis in CT Images Based on the Central Cross-section Theorem]

PURPOSE

The noise power spectrum (NPS) of a CT scanner is commonly measured from a single noise image. However, since CT images are three-dimensional (3D) volume data, they have 3D noise characteristics (3D-NPS). In this study, we clarify the relationship among NPSs measured by various approaches in NPS analysis based on the central slice theorem. Its validity is verified by the NPS measurements using actual 3D noise data.

METHODS

We defined the NPS_z-projection(f_x, f_y) that was calculated by the 2D Fourier transform (FT) from the 2D projection of 3D noise data in the patient longitudinal direction, the 3D-NPS(f_x, f_y, f_z) that was calculated by the 3D-FT from the 3D noise data, and the 2D-NPS(f_x, f_y) that was calculated by the 2D-FT from a single noise image; f_x, f_y, and f_z are spatial frequencies corresponding to the axes of x, y, and z in the reconstructed CT volume, respectively. Based on the central slice theorem, we described that the NPS_z-projection(f_x, f_y=0) was equal to the 3D-NPS(f_x, f_y=0, f_z=0), and the NPS(2D-NPS(f_x, f_y=0)) was different from the 3D-NPS(f_x, f_y=0, f_z=0). To verify them, we compared the NPSs calculated from actual 3D noise data that were obtained using a cylindrical water phantom.

RESULTS

The 3D-NPS(f_x, f_y=0, f_z=0) matched the NPS_z-projection(f_x, f_y=0) and was different from the 2D-NPS(f_x, f_y=0).

CONCLUSION

Based on the central slice theorem, we clarified the relationship among NPSs measured by various approaches in NPS analysis; it is important to understand this and then select an appropriate noise data handling and NPS measurement method.

Task-specific spatial resolution properties of iterative and deep learning-based reconstructions in computed tomography: Comparison using tasks assuming small and large enhanced vessels

PURPOSE

The present study aims to evaluate TTFs of deep-learning-based image reconstruction (DLIR) and iterative reconstruction (IR) in computed tomography (CT) using a conventional task with a rod object with a diameter of 30 mm and a newly-proposed task with a wire of 1 mm in diameter, simulating large and small enhanced vessels, respectively.

METHODS

The rod or wire phantom made of a material equivalent to diluted iodine that exhibits about 270 Hounsfield unit (HU) was placed inside a 30-cm water phantom. In-plane and z-directional TTFs were measured for the rod using the circular edge (CE) and plane edge (PE) methods, respectively. By using the wire (iodine wire: IW), in-plane and z-directional TTFs were measured using Fourier transform (IW method). TTFs of filtered back projection (FBP), IR, and DLIR of a 256-row CT system and FBP and IR of a 64-row CT system were evaluated with CT dose indices of 10 and 5 mGy.

RESULTS

For DLIR and IR, TTFs measured using the IW method were notably lower than those using the CE (or PE) method; moreover, they were also lower than those of corresponding FBP, indicating that the small enhanced vessels with a diameter of about 1 mm would be blurred with both DLIR and IR.

CONCLUSIONS

The proposed IW method has turned out to be effective to evaluate TTFs for small enhanced vessels, which have not been properly evaluated by the CE or PE method conventionally recommended.

Virtual monoenergetic images from dual-energy CT: systematic assessment of task-based image quality performance

Background

To compare task-based image quality (TB-IQ) among virtual monoenergetic images (VMI) and linear-blended images (LBI) from dual-energy CT as a function of contrast task, radiation dose, size, and lesion diameter.

Methods

A TB-IQ phantom (Mercury Phantom 4.0, Sun Nuclear Corporation) was imaged on a third-generation dual-source dual-energy CT with 100/Sn150 kVp at three volume CT dose levels (5, 10, 15 mGy). Three size sections (diameters 16, 26, 36 cm) with subsections for image noise and spatial resolution analysis were used. High-contrast tasks (e.g., calcium-containing stone and vascular lesion) were emulated using bone and iodine inserts. A low-contrast task (e.g., low-contrast lesion or hematoma) was emulated using a polystyrene insert. VMI at 40-190 keV and LBI were reconstructed. Noise power spectrum (NPS) determined the noise magnitude and texture. Spatial resolution was assessed using the task-transfer function (TTF) of the three inserts. The detectability index (d') served as TB-IQ metric.

Results

Noise magnitude increased with increasing phantom size, decreasing dose, and decreasing VMI-energy. Overall, noise magnitude was higher for VMI at 40-60 keV compared to LBI (range of noise increase, 3-124%). Blotchier noise texture was found for low and high VMIs (40-60 keV, 130-190 keV) compared to LBI. No difference in spatial resolution was observed for high contrast tasks. d' increased with increasing dose level or lesion diameter and decreasing size. For high-contrast tasks, d' was higher at 40-80 keV and lower at high VMIs. For the low-contrast task, d' was higher for VMI at 70-90 keV and lower at 40-60 keV.

Conclusions

Task-based image quality differed among VMI-energy and LBI dependent on the contrast task, dose level, phantom size, and lesion diameter. Image quality could be optimized by tailoring VMI-energy to the contrast task. Considering the clinical relevance of iodine, VMIs at 50-60 keV could be proposed as an alternative to LBI.

A novel method to assess the spatiotemporal image quality in fluoroscopy

Objectives. The planar formulation of the noise equivalent quanta (NEQ) and detective quantum efficiency (DQE) used to assess the image quality of projection images does not deal with the influence of temporal resolution on signal blurring and image noise. These metrics require correction factors based on temporal resolution when used for dynamic imaging systems such as fluoroscopy. Additionally, the standard NEQ and detector DQE are determined on pre-processed images in scatter-free conditions for effective energies produced by additional aluminium or copper filters that are not representative of clinical fluoroscopic procedures. In this work, we developed a method to measure 'frame NEQ' and 'frame system DQE' which include the temporal frequency bandwidth and consider the anti-scatter grid, the detector and the image processing procedures for beam qualities with scatter fractions representative of clinical use.Approach. We used a solid water phantom to simulate a patient and a thin copper disc to measure the spatial resolution. The copper disc, set in uniform rectilinear motion in the image plane, assessed the temporal resolution. These new metrics were tested on two fluoroscopy systems, a C-arm and a floor-mounted cardiology, for multiple parameters: phantom thicknesses from 5 to 20 cm, frame rates from 3 to 30 fps, spatial and temporal image processing of different weights.Main results.The frame NEQ correctly described the image quality for different scatter conditions, temporal resolutions and image processing techniques. The frame system DQE varied between 0.38 and 0.65 within the different beam and scatter conditions, and correctly mitigated the influence of spatial and temporal image processing.Significance.This study introduces and validates an unbiased formulation of in-plane NEQ and system DQE to assess the spatiotemporal image quality of fluoroscopy systems.

Comparison of Low Dose Performance of Photon-Counting and Energy Integrating CT

RATIONALE AND OBJECTIVES

The purpose of this study was to investigate the potential of photon-counting CT (PCCT) to improve quantitative image quality for low dose imaging compared to energy-integrating detector CT (EID CT).

MATERIALS AND METHODS

An investigational scanner (Siemens, Germany) with PCCT and EID CT subsystems was used to compare image quality performance at four dose levels: 1.7, 2, 4, 6 mGy CTDIvol, all at or below current dose values used for conventional abdominal CT. A CT quality control phantom with a homogeneous section for noise measurements and a section with cylindrical inserts of air (-910 HU), polystyrene (50 HU), acrylic (205 HU), and Teflon (1000 HU) was imaged and characterized in terms of noise, resolution, contrast-to-noise ratio (CNR), and detectability index. A second phantom with a 30 cm diameter was also imaged containing iodine solutions ranging from 0.125 to 8 mg I/mL. CNR of the iodine vials was computed as a function of CT dose and iodine concentration.

RESULTS

With resolution unaffected by dose in both PCCT and EID CT, PCCT images exhibited 22.1-24.0% improvement in noise across dose levels evaluated. This noise improvement translated into a 29-41% improvement in CNR and 20-36% improvement in detectability index. For iodine contrast, PCCT images had a higher CNR for all combinations of iodine contrast and dose evaluated.

CONCLUSION

For the conditions studied, PCCT exhibited superior image quality compared to EID CT. For iodine detection, PCCT offered a notable advantage with improved CNR at all doses and iodine concentration levels.

FREQUENCY RESPONSE AND DISTORTION PROPERTIES OF RECONSTRUCTION ALGORITHMS IN COMPUTED TOMOGRAPHY

Denoising reconstruction techniques can introduce nonlinear properties into computed tomography (CT) systems. These nonlinear algorithms introduce distortion which affects the assessment of the resolution of the system. The purpose of the present study was to decouple and investigate amplitude modulation and waveform distortion in reconstruction algorithms in CT. The methodology developed by Wells, J. R. and Dobbins, J. T. III [Frequency response and distortion properties of nonlinear image processing algorithms and the importance of imaging context. Med. Phys. 40, 091906 (2013)] was adapted to CT reconstruction algorithms. The CT simulating program ASTRA Toolbox© for MATLAB™ was used for the reconstruction of the sinusoidal wave functions. Filtered back projection and the simultaneous iterative reconstruction technique were investigated with simple nonlinear mechanisms: a median filter and a non-negative constraint, respectively. The native reconstruction algorithms were not free from nonlinear waveform distortion, however, none of the metrics showed any dependence on the contrast-to-noise ratio (CNR). Furthermore, the algorithms including nonlinear mechanisms showed a clear and specific CNR dependence, indicating the necessity for distortion analysis in nonlinear CT reconstruction.

Method for determining slice sensitivity profile of iterative reconstruction CT images using low-contrast sphere phantom

A novel method for measuring the slice sensitivity profile (SSP) of computed tomography (CT) images reconstructed using an iterative reconstruction (IR) algorithm is proposed herein. A phantom that included a low-contrast spherical object was scanned and consecutive cross-sectional images were reconstructed. The mean CT values in a region including the sphere were measured for all images and plotted as a function of slice position along the longitudinal [Formula: see text] direction to yield a mean CT value profile [Formula: see text]. Next, we numerically generated an object function corresponding to the sphere and obtained the mean CT value profile [Formula: see text]. Subsequently, the SSP was modeled as a product of the Gaussian and cosine functions. We convolved [Formula: see text] with the modeled SSP to obtain [Formula: see text]. The difference between [Formula: see text] and [Formula: see text] was evaluated using the root mean square error (RMSE), which was minimized via optimization of the SSP model parameters. To validate the methodology, we first used filtered back projection (FBP) images to compare the SSPs determined using the proposed and standard coin methods. Subsequently, the proposed method was applied to measure the SSPs of four types of IR algorithms in two scanners. The SSPs of the FBP images determined using the proposed and coin methods showed good agreement. Additionally, in the SSP measurements using the proposed method, [Formula: see text] agreed well with [Formula: see text] for every IR algorithm. The RMSEs for all measurements were less than 0.7 HU, indicating the accuracy of the SSPs. Thus, the proposed method is effective for obtaining valid SSPs.

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.

Photon-Counting Multienergy Computed Tomography With Spectrally Optimized Contrast Media for Plaque Removal and Stenosis Assessment

PURPOSE

The aim of this study was to systematically evaluate the potential to combine investigational contrast media with spectrally optimized energy-thresholding of photon-counting detector computed tomography (PCCT) for subtraction of calcified plaques in a coronary artery stenosis phantom.

METHODS

A small vessel phantom containing 3 fillable tubes (diameter, 3 mm each) with calcified plaques was placed into an anthropomorphic chest phantom. The plaques had incremental thicknesses ranging from 0.3 to 2.7 mm, simulating vessel stenoses ranging from 10% to 90% of the lumen diameter. The phantom was filled with 5 different investigational contrast media (iodine, bismuth, hafnium, holmium, and tungsten) at equal mass concentrations (15 mg/mL) and was imaged on a prototype PCCT at 140 kVp using optimized, contrast media-dependent energy thresholds. Contrast maps (CMs) were reconstructed for each contrast medium by applying a linear 2-material decomposition algorithm. Image noise magnitude and noise texture of CM were compared among the contrast media using the noise power spectrum. Two blinded readers independently rated the vessel lumen visualization on short-axis and the overall subjective image quality on long-axis CM relative to iodine as the reference standard. Four readers determined the highest degree of stenosis that could be assessed with high diagnostic confidence on long-axis CM.

RESULTS

Average image noise on CM was lower for tungsten (49 HU) and hafnium (62 HU) and higher for bismuth (81 HU) and holmium (165 HU) compared with iodine (78 HU). Noise texture of CM was similar among the contrast media. Interreader agreement for vessel lumen visualization on short-axis CM ranged from moderate to excellent (k = 0.567-0.814). Compared with iodine, lumen visualization of each reader was improved using tungsten (P < 0.001 for both readers), similar to improved using hafnium (P = 0.008, P = 0.29), similar using bismuth (P = 0.38, P = 0.69), and decreased using holmium (both, P < 0.001). Overall subjective image quality was similar for holmium and superior for tungsten, hafnium, and bismuth as compared with iodine. Higher-degree stenoses were evaluable with high confidence using tungsten (mean, 70%; interquartile range, 70%-70%), bismuth (70%; 60%-70%), and hafnium (75%; 70%-80%) compared with iodine (50%; 50%-60%) and holmium (50%; 50%-60%).

CONCLUSIONS

Spectral optimization in PCCT combined with investigational contrast media can improve calcium subtraction and stenosis assessment in small vessels. Contrast maps of tungsten and, to a lesser extent, hafnium as contrast media yielded superior image noise properties and improved vessel lumen visualization, along with a higher subjective image quality compared with the reference standard iodine.

Harmonization of in-plane resolution in CT using multiple reconstructions from single acquisitions

PURPOSE

To providea methodology that removes the spatial variability of in-plane resolution using different CT reconstructions. The methodology does not require any training, sinogram, or specific reconstruction method.

METHODS

The methodology is formulated as a reconstruction problem. The desired sharp image is modeled as an unobservable variable to be estimated from an arbitrary number of observations with spatially variant resolution. The methodology comprises three steps: (1) density harmonization, which removes the density variability across reconstructions; (2) point spread function (PSF) estimation, which estimates a spatially variant PSF with arbitrary shape; (3) deconvolution, which is formulated as a regularized least squares problem. The assessment was performed with CT scans of phantoms acquired with three different Siemens scanners (Definition AS, Definition AS+, Drive). Four low-dose acquisitions reconstructed with backprojection and iterative methods were used for the resolution harmonization. A sharp, high-dose (HD) reconstruction was used as a validation reference. The different factors affecting the in-plane resolution (radial, angular, and longitudinal) were studied with regression analysis of the edge decay (between 10% and 90% of the edge spread function (ESF) amplitude).

RESULTS

Results showed that the in-plane resolution improves remarkably and the spatial variability is substantially reduced without compromising the noise characteristics. The modulated transfer function (MTF) also confirmed a pronounced increase in resolution. The resolution improvement was also tested by measuring the wall thickness of tubes simulating airways. In all scanners, the resolution harmonization obtained better performance than the HD, sharp reconstruction used as a reference (up to 50 percentage points). The methodology was also evaluated in clinical scans achieving a noise reduction and a clear improvement in thin-layered structures. The estimated ESF and MTF confirmed the resolution improvement.

CONCLUSION

We propose a versatile methodology to reduce the spatial variability of in-plane resolution in CT scans by leveraging different reconstructions available in clinical studies. The methodology does not require any sinogram, training, or specific reconstruction, and it is not limited to a fixed number of input images. Therefore, it can be easily adopted in multicenter studies and clinical practice. The results obtained with our resolution harmonization methodology evidence its suitability to reduce the spatially variant in-plane resolution in clinical CT scans without compromising the reconstruction's noise characteristics. We believe that the resolution increase achieved by our methodology may contribute in more accurate and reliable measurements of small structures such as vasculature, airways, and wall thickness.

Dose reduction potential of vendor-agnostic deep learning model in comparison with deep learning-based image reconstruction algorithm on CT: a phantom study

Objectives

To compare the dose reduction potential (DRP) of a vendor-agnostic deep learning model (DLM, ClariCT.AI) with that of a vendor-specific deep learning–based image reconstruction algorithm (DLR, TrueFidelity™).

Methods

Computed tomography (CT) images of a multi-sized image quality phantom (Mercury v4.0) were acquired under six radiation dose levels (0.48/0.97/1.93/3.87/7.74/15.47 mGy) and were reconstructed using filtered back projection (FBP) and three strength levels of the DLR (low/medium/high). The FBP images were denoised using the DLM. For all DLM and DLR images, the detectability index (d′) (a task-based detection performance metric) was obtained, under various combinations of three target sizes (10/5/1 mm), five inlets (CT value difference with the background; −895/50/90/335/1000 HU), five phantom diameters (36/31/26/21/16 cm), and six radiation dose levels. Dose reduction potential (DRP) measures the dose reduction made by using DLM or DLR, while yielding d′ equivalent to that of FBP at full dose.

Results

The DRPs of the DLM, DLR-low, DLR-medium, and DLR-high were 86% (81–88%), 60% (46–67%), 76% (60–81%), and 87% (78–92%), respectively. For 10-mm targets, the DRP of the DLM (87%) was higher than that of all DLR algorithms (58–86%). However, for smaller targets (5 mm/1 mm), the DRPs of the DLR-high (89/88%) were greater than those of the DLM (87/84%).

Conclusion

The dose reduction potential of the vendor-agnostic DLM was shown to be comparable to that of the vendor-specific DLR at high strength and superior to those of the DLRs at medium and low strengths.

Key Points

• DRP of the vendor-agnostic model was comparable to that of high-strength vendor-specific model and superior to those of medium- and low-strength models.

• Under various radiation dose levels, the deep learning model shows higher detectability indexes compared to FBP.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-021-08199-9.

Determination of the optimal range for virtual monoenergetic images in dual-energy CT based on physical quality parameters

PURPOSE

Virtual monoenergetic images (VMI) obtained from Dual-Energy Computed Tomography (DECT) with iodinated contrast are used in radiotherapy of the Head and Neck to improve the delineation of target volumes and organs at-risk (OAR). The energies used to vary from 40 to 70 keV, but noise at low keV and the use of Single Energy CT (SECT) at low kV_p settings may shrink this interval. There is no guide about how to find out the optimal range where VMI has a significant improvement related to SECT images. Our study proposes a procedure to determine this optimal range, based on common image quality parameters, and establishes this range in a Siemens Somatom Confidence and a Head and Neck protocol.

METHODS

We compared the quality of the VMI series at 40-60 keV versus single X-ray tube voltage computed tomography (SECT) at 80 and 120 kV_p . Our reference was 120 kV_p . DECT images were sequentially acquired using the Siemens Somatom Confidence RT Pro CT according to the head and neck protocol in our department. VMI series were constructed using the Syngo Via software Monoenergetic+ algorithm. Quality parameters were: image uniformity, high- and low-contrast resolution, noise, and sensitivity to the iodinated contrast. We used the Catphan 604 phantom for quality control, except when assessing iodine sensitivity. To evaluate high contrast resolution, we calculated the modulation transfer function (MTF) using the point spread function estimation of a point bead and the slanted edge methods. For the low-contrast resolution, we used a statistical method for assessing differences between contrast structures and local noise. To measure the absolute value of noise and compare its texture, we used the standard deviation and the noise power spectrum. We measured iodine sensitivity by dissolving the Optiray Ultraject iodinated contrast in water in concentrations of 0 to 4500 mg/l and then compared the contrast to noise ratio (CNR) and analyzed the linear correlation between concentration and HU.

RESULTS

The entire series met the minimum quality requirements. However, the one at 40 keV presented uniformity at the limits of acceptability. The high- and low-contrast resolutions were similar between series. The noise of the VMI series decreased with increasing energy, while sensitivity to the contrast displayed the opposite behavior. All series showed linearity of HUs from very low iodine concentrations. Images at 60 keV presented lower iodine sensitivity than SECT at 80 kV_p , while those at 55 keV were similar to them.

CONCLUSIONS

Our method of image comparison based on standard quality parameters in phantom gave clear results about the optimal range and can be used as a guide to characterize any other DECT imaging protocols. The optimal range for using VMI images in iodinated contrasts in the Siemens system was 45-55 keV. Lower energies lacked noise and uniformity, while higher ones could be substituted by SECT images at low kilovoltage (80 kV_p ).

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.

Reduction of streak artifacts caused by low photon counts utilizing an image-based forward projection in computed tomography

BACKGROUND

The streak artifacts in computed tomography (CT) images caused by low photon counts are known to be effectively suppressed by raw-data-based techniques. This study aims to propose a technique to reduce the streak artifact without accessing the raw data.

METHODS

The proposed streak artifact reduction (SAR) technique consists of three steps: numerical forward projection to a CT image, adaptive filtering of the generated sinogram, and image reconstruction from the processed sinogram. The authors have expanded the two-dimensional method (2D-SAR) to three dimensions (3D-SAR) by using consecutive CT images. The modulation transfer function (MTF), the image noise (standard deviation), and the visibility of comb-shaped objects were evaluated at a low dose of 5 mGy. Using anthropomorphic abdominal and chest phantoms, CT images and the artifact index (AI) were compared between 3D-SAR and two types of iterative reconstruction (IR).

RESULTS

Sufficient artifact reductions associated with 54% and 61% reduction of noise for 2D- and 3D-SAR, respectively, were obtained in the phantom images, although the 50%MTF decreased by 28%. The visibility of the combs was improved with both the 2D- and 3D-SAR methods. The AI results of 3D-SAR were better than one type of IR and almost equal to the other type of IR, which was consistent with observed artifacts.

CONCLUSION

Both 2D-SAR and 3D-SAR have turned out to be effective in reducing streak artifacts. The proposed technique will be an effective tool since it needs no raw data, and thus can be applied to any CT images produced by a wide variety of CT systems.

The effect of tube focal spot size and acquisition mode on task-based image quality performance of a GE revolution HD dual energy CT scanner

PURPOSE

To assess the task-based performance of images obtained under different focal spot size and acquisition mode on a dual-energy CT scanner.

METHODS

Axial CT image series of the Catphan phantom were obtained using a tube focus at different sizes. Acquisitions were performed in standard single-energy, high resolution (HR) and dual-energy modes. Images were reconstructed using conventional and high definition (HD) kernels. Task-based transfer function at the 50% level (TTF_50%) for teflon, delrin, low density polyethylene (LDPE) and acrylic, as well as image noise and noise texture, were assessed across all focal spots and acquisition modes using Noise Power Spectrum (NPS) analysis. A non-prewhitening mathematical observer model was used to calculate detectability index (d_NPW^').

RESULTS

TTF_50% degraded with increasing focal spot size. TTF_50% ranged from 0.67 mm^-1 for teflon to 0.25 mm^-1 for acrylic. For standard kernel, image noise and NPS-determined average spatial frequency were 8.3 HU and 0.29 mm^-1, respectively in single-energy, 12.0 HU and 0.37 mm^-1 in HR, and 7.9 HU and 0.26 mm^-1 in dual-energy mode. For standard kernel, d_NPW^' was 61 in single-energy and HR mode and reduced to 56 in dual-energy mode.

CONCLUSIONS

The task-based image quality assessment metrics have shown that spatial resolution is higher for higher image contrast materials and detectability is higher in the standard single-energy mode compared to HR and dual-energy mode. The results of the current study provide CT operators the required knowledge to characterize their CT system towards the optimization of its clinical performance.

Toward a comparison and an optimization of CT protocols using new metrics of dose and image quality part I: prediction of human observers using a model observer for detection and discrimination tasks in low-dose CT images in various scanning conditions

In the context of reducing the patient dose coming from CT scanner examinations without penalizing the diagnosis, the assessment of both patient dose and image quality (IQ) with relevant metrics is crucial. The present study represents the first stage in a larger work, aiming to compare and optimize CT protocols using dose and IQ new metrics. We proposed here to evaluate the capacity of the Non-PreWhitening matched filter with an eye (NPWE) model observer to be a robust and accurate estimation of IQ. We focused our work on two types of clinical tasks: a low contrast detection task and a discrimination task. We designed a torso-shaped phantom, including Plastic Water^®slabs with cylindrical inserts of different diameters, sections and compositions. We led a human observer study with 13 human observers on images acquired in multiple irradiation and reconstruction scanning conditions (voltage, pitch, slice thickness, noise level of the reconstruction algorithm, energy level in dual-energy mode and dose), to evaluate the behavior of the model observer compared to the human responses faced to changing conditions. The model observer presented the same trends as the human observers with generally better results. We rescaled the NPWE model on the human responses by scanning conditions (kVp, pitch, slice thickness) to obtain the best agreement between both observer types, estimated using the Bland-Altman method. The impact of some scanning parameters was estimated using the correct answer rate given by the rescaled NPWE model, for both tasks and each insert size. In particular, the comparison between the dual-energy mode at 74 keV and the single-energy mode at 120 kVp showed that, if the 120 kVp voltage provided better results for the smallest insert at the lower doses for both tasks, their responses were equivalent in many cases.

Technical Note: A simple method for measuring the slice sensitivity profile of iteratively reconstructed CT images using a non-slanted edge plane

PURPOSE

A method for measuring the slice sensitivity profile (SSP) of computed tomography (CT) images reconstructed with iterative reconstruction (IR) algorithms was reported by the AAPM Task Group 233 (TG233). In this method, the phantom plane edge is slightly slanted with respect to the scan plane to obtain a composite oversampled edge-spread function (ESF). However, it is expected that a fine-sampled ESF can be obtained directly from images reconstructed with a small slice increment without slanting the edge plane. This study aimed to investigate the validity of using a non-slanted edge plane.

METHODS

In the proposed non-slanted edge method, the phantom was positioned so that the plane edge was perpendicular to the longitudinal z-axis, and images were reconstructed with a 1-mm slice thickness and 0.1-mm increment. The mean CT value was obtained in each slice and plotted as a function of slice position along the z-axis, thereby generating the ESF. The SSP was calculated from the ESF by differentiation. In the TG 233-recommended slanted edge method, the SSP was obtained by following the procedure described in the TG233 report. To validate the methodology, we first used filtered back projection (FBP) images to compare SSPs obtained using the non-slanted edge method, slanted edge method, and a standard method using a high-contrast thin object (coin). Next, for two types of IR algorithms, we compared the SSPs obtained using the non-slanted and slanted edge methods.

RESULTS

For the FBP images, the SSP measured using the non-slanted edge method agreed well with SSPs measured using the coin and slanted edge methods. For the IR images, the SSPs measured using the non-slanted and slanted edge methods showed good agreement.

CONCLUSIONS

The non-slanted edge method was demonstrated to be valid. The simplicity and practicality of the method allows routine and accurate determination of the SSP.

Assessment of the image quality of virtual monochromatic spectral computed tomography images: a phantom study considering object contrast, radiation dose, and frequency characteristics

Fast kilovoltage (kVp)-switching technology cannot obtain conventional 120 kVp images; thus, 70 keV virtual monochromatic spectral computed tomography (CT) images (VMSI) are generally used. The contrast-to-noise ratio (CNR) is used to evaluate the image quality of VMSI; however, CNR does not include frequency characteristics. The present study aimed to investigate the evaluation methods of VMSI considering frequency characteristics by comparing the image quality of 70 keV VMSI with that of conventional 120 kVp images. The evaluated object contrasts were 70 and 300 Hounsfield units (HU). Scans used two radiation dose levels: low (LD) and standard (SD). The volume CT dose index of LD and SD was 4.8- and 12 mGy, respectively. Images were reconstructed by filtered back projection, evaluating CNR, noise power spectrum (NPS), task transfer function (TTF), and system performance (SP) function calculated as TTF²/ NPS. The total NPS values (spatial frequency range: 0.2 ~ 0.4 mm^-1) of 70 keV VMSI were higher than those of 120 kVp images. The spatial frequency TTF values that reached 10% (f_10%) of the 70 keV VMSI changed based on object contrast. For the low-contrast condition, a lower f_10% was observed with 70 keV VMSI. The CNR of 70 keV VMSI was comparable to that of 120 kVp images in low- and high-contrast conditions. However, for 70 keV VMSI, SP of low-contrast was low, and SP of high-contrast was high, compared with those of 120 kVp images. This study suggested that only CNR was not sufficient to evaluate the image quality of VMSI; thus, evaluation methods considering frequency characteristics should be used.

Performance of clinically available deep learning image reconstruction in computed tomography: a phantom study

Purpose: To assess the physical performance of deep learning image reconstruction (DLIR) compared with those of filtered back projection (FBP) and iterative reconstruction (IR) and to estimate the dose reduction potential of the technique. Approach: A cylindrical water bath phantom with a diameter of 300 mm including two rods composed of acrylic and soft tissue-equivalent material was scanned using a clinical computed tomography (CT) scanner at four dose levels (CT dose index of 20, 15, 10, and 5 mGy). Phantom images were reconstructed using FBP, DLIR, and IR. The in-plane and z axis task transfer functions (TTFs) and in-plane noise power spectrum (NPS) were measured. The dose reduction potential was estimated by evaluating the system performance function calculated from TTF and NPS. The visibilities of a bar pattern phantom placed in the same water bath phantom were compared. Results: The use of DLIR resulted in a notable decrease in noise magnitude. The shift in peak NPS frequency was reduced compared with IR. Preservation of in-plane TTF was superior using DLIR than using IR. The estimated dose reduction potentials of DLIR and IR were 39% to 54% and 19% to 29%, respectively. However, the z axis resolution was decreased with DLIR by 6% to 21% compared with FBP. The bar pattern visibilities were approximately consistent with the TTF results in both planes. Conclusions: The in-plane edge-preserving noise reduction performance of DLIR is superior to that of IR. Moreover, DLIR enables approximately half-dose acquisitions with no deterioration in noise texture in cases that permit some z axis resolution reduction.

Are quantitative features of lung nodules reproducible at different CT acquisition and reconstruction parameters?

Consistency and duplicability in Computed Tomography (CT) output is essential to quantitative imaging for lung cancer detection and monitoring. This study of CT-detected lung nodules investigated the reproducibility of volume-, density-, and texture-based features (outcome variables) over routine ranges of radiation dose, reconstruction kernel, and slice thickness. CT raw data of 23 nodules were reconstructed using 320 acquisition/reconstruction conditions (combinations of 4 doses, 10 kernels, and 8 thicknesses). Scans at 12.5%, 25%, and 50% of protocol dose were simulated; reduced-dose and full-dose data were reconstructed using conventional filtered back-projection and iterative-reconstruction kernels at a range of thicknesses (0.6-5.0 mm). Full-dose/B50f kernel reconstructions underwent expert segmentation for reference Region-Of-Interest (ROI) and nodule volume per thickness; each ROI was applied to 40 corresponding images (combinations of 4 doses and 10 kernels). Typical texture analysis metrics (including 5 histogram features, 13 Gray Level Co-occurrence Matrix, 5 Run Length Matrix, 2 Neighboring Gray-Level Dependence Matrix, and 3 Neighborhood Gray-Tone Difference Matrix) were computed per ROI. Reconstruction conditions resulting in no significant change in volume, density, or texture metrics were identified as "compatible pairs" for a given outcome variable. Our results indicate that as thickness increases, volumetric reproducibility decreases, while reproducibility of histogram- and texture-based features across different acquisition and reconstruction parameters improves. To achieve concomitant reproducibility of volumetric and radiomic results across studies, balanced standardization of the imaging acquisition parameters is required.

A Clinically Driven Task-Based Comparison of Photon Counting and Conventional Energy Integrating CT for Soft Tissue, Vascular, and High-Resolution Tasks

Photon-counting CT detectors are the next step in advancing CT system development and will replace the current energy integrating detectors (EID) in CT systems in the near future. In this context, the performance of PCCT was compared to EID CT for three clinically relevant tasks: abdominal soft tissue imaging, where differentiating low contrast features is important; vascular imaging, where iodine detectability is critical; and, high-resolution skeletal and lung imaging. A multi-tiered phantom was imaged on an investigational clinical PCCT system (Siemens Healthineers) across different doses using three imaging modes: macro and ultra-high resolution (UHR) PCCT modes and EID CT. Images were reconstructed using filtered backprojection and soft tissue (B30f), vascular (B46f), or high-resolution (B70f; U70f for UHR) kernels. Noise power spectra, task transfer functions, and detectability index were evaluated. For a soft tissue task, PCCT modes showed comparable noise and resolution with improved contrast-to-noise ratio. For a vascular task, PCCT modes showed lower noise and improved iodine detectability. For a high resolution task, macro mode showed lower noise and comparable resolution while UHR mode showed higher noise but improved spatial resolution for both air and bone. PCCT offers competitive advantages to EID CT for clinical tasks.

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.

Image Quality Assessment of Abdominal CT by Use of New Deep Learning Image Reconstruction: Initial Experience

OBJECTIVE. The purpose of this study was to perform quantitative and qualitative evaluation of a deep learning image reconstruction (DLIR) algorithm in contrast-enhanced oncologic CT of the abdomen. MATERIALS AND METHODS. Retrospective review (April-May 2019) of the cases of adults undergoing oncologic staging with portal venous phase abdominal CT was conducted for evaluation of standard 30% adaptive statistical iterative reconstruction V (30% ASIR-V) reconstruction compared with DLIR at low, medium, and high strengths. Attenuation and noise measurements were performed. Two radiologists, blinded to examination details, scored six categories while comparing reconstructions for overall image quality, lesion diagnostic confidence, artifacts, image noise and texture, lesion conspicuity, and resolution. RESULTS. DLIR had a better contrast-to-noise ratio than 30% ASIR-V did; high-strength DLIR performed the best. High-strength DLIR was associated with 47% reduction in noise, resulting in a 92-94% increase in contrast-to-noise ratio compared with that of 30% ASIR-V. For overall image quality and image noise and texture, DLIR scored significantly higher than 30% ASIR-V with significantly higher scores as DLIR strength increased. A total of 193 lesions were identified. The lesion diagnostic confidence, conspicuity, and artifact scores were significantly higher for all DLIR levels than for 30% ASIR-V. There was no significant difference in perceived resolution between the reconstruction methods. CONCLUSION. Compared with 30% ASIR-V, DLIR improved CT evaluation of the abdomen in the portal venous phase. DLIR strength should be chosen to balance the degree of desired denoising for a clinical task relative to mild blurring, which increases with progressively higher DLIR strengths.

A pitfall of using the circular-edge technique with image averaging for spatial resolution measurement in iteratively reconstructed CT images

The circular‐edge technique using a low‐contrast cylindrical object is commonly used to measure the modulation transfer functions (MTFs) in computed tomography (CT) images reconstructed with iterative reconstruction (IR) algorithms. This method generally entails averaging multiple images of the cylinder to reduce the image noise. We suspected that the cylinder edge shape depicted in the IR images might exhibit slight deformation with respect to the true shape because of the intrinsic nonlinearity of IR algorithms. Image averaging can reduce the image noise, but does not effectively improve the deformation of the edge shape; thereby causing errors in the MTF measurements. We address this issue and propose a method to correct the MTF. We scanned a phantom including cylindrical objects with a CT scanner (Ingenuity Elite, Philips Healthcare). We obtained cylinder images with iterative model reconstruction (IMR) algorithms. The images suggested that the depicted edge shape deforms and fluctuates depending on slice positions. Because of this deformation, image averaging can potentially cause additional blurring. We define the deformation function D that describes the additional blurring, and obtain D by analyzing multiple images. The MTF measured by the circular‐edge method (referred to as MTF') can be thought of as the multiplication of the true MTF by the Fourier transformation (FT) of D. We thus obtain the corrected MTF (MTF_corrected) by dividing MTF' by the FT of D. We validate our correction method by comparing the calculated images based on the convolution theorem using MTF' and MTF_corrected with the actual images obtained with the scanner. The calculated image using MTF_corrected is more similar to the actual image compared with the image calculated using MTF', particularly in edge regions. We describe a pitfall in MTF measurement using the circular‐edge technique with image averaging, and suggest a method to correct it.