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

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

THE EFFECT OF OUT-OF-PLANE PATIENT SHIELDING ON CT RADIATION EXPOSURE AND TUBE CURRENT MODULATIONS: A PHANTOM STUDY ACROSS THREE VENDORS

The aim of this study was to evaluate how out-of-plane patient shielding affects radiation exposure parameters and tube current modulation on different vendors' computed tomography (CT) scanners. Helical CT scans were performed using two homogenous phantoms to mimic patient attenuation. Four CT scanners from three vendors were investigated by varying the distance of the patient shield from the border of the imaging volume. Scans were performed with a shield placed before and after the localizer. Changes in volume computed tomography dose index (CTDIvol), dose-length product (DLP) and tube current-time products were studied. Out-of-field lead shield increased the CTDIvol and DLP values for each scanner at least for one scan setting when the shield was present in the localizer. The most notable changes were recorded with >1.3 pitch values when the shield was closest to the scanned volume (2.5 cm), and the scan direction was towards the shield. The usage of patient shields in the localizer CT scans can disturb TCM even when placed 7.5 cm away from the edge of the scan.

Virtual monochromatic imaging reduces beam hardening artefacts in cardiac interior photon counting computed tomography: a phantom study with cadaveric specimens

In interior cardiac computed tomography (CT) imaging, the x-ray beam is collimated to a limited field-of-view covering the heart volume, which decreases the radiation exposure to surrounding tissues. Spectral CT enables the creation of virtual monochromatic images (VMIs) through a computational material decomposition process. This study investigates the utility of VMIs for beam hardening (BH) reduction in interior cardiac CT, and further, the suitability of VMIs for coronary artery calcium (CAC) scoring and volume assessment is studied using spectral photon counting detector CT (PCD-CT).Ex vivocoronary artery samples (N = 18) were inserted in an epoxy rod phantom. The rod was scanned in the conventional CT geometry, and subsequently, the rod was positioned in a torso phantom and re-measured in the interior PCD-CT geometry. The total energy (TE) 10-100 keV reconstructions from PCD-CT were used as a reference. The low energy 10-60 keV and high energy 60-100 keV data were used to perform projection domain material decomposition to polymethyl methacrylate and calcium hydroxylapatite basis. The truncated basis-material sinograms were extended using the adaptive detruncation method. VMIs from 30-180 keV range were computed from the detruncated virtual monochromatic sinograms using filtered back projection. Detrending was applied as a post-processing method prior to CAC scoring. The results showed that BH artefacts from the exterior structures can be suppressed with high (≥100 keV) VMIs. With appropriate selection of the monoenergy (46 keV), the underestimation trend of CAC scores and volumes shown in Bland-Altman (BA) plots for TE interior PCD-CT was mitigated, as the BA slope values were -0.02 for the 46 keV VMI compared to -0.21 the conventional TE image. To conclude, spectral PCD-CT imaging using VMIs could be applied to reduce BH artefacts interior CT geometry, and further, optimal selection of VMI may improve the accuracy of CAC scoring assessment in interior PCD-CT.

Generative adversarial networks improve interior computed tomography angiography reconstruction

In interior computed tomography (CT), the x-ray beam is collimated to a limited field-of-view (FOV) (e.g. the volume of the heart) to decrease exposure to adjacent organs, but the resulting image has a severe truncation artifact when reconstructed with traditional filtered back-projection (FBP) type algorithms. In some examinations, such as cardiac or dentomaxillofacial imaging, interior CT could be used to achieve further dose reductions. In this work, we describe a deep learning (DL) method to obtain artifact-free images from interior CT angiography. Our method employs the Pix2Pix generative adversarial network (GAN) in a two-stage process: (1) An extended sinogram is computed from a truncated sinogram with one GAN model, and (2) the FBP reconstruction obtained from that extended sinogram is used as an input to another GAN model that improves the quality of the interior reconstruction. Our double GAN (DGAN) model was trained with 10 000 truncated sinograms simulated from real computed tomography angiography slice images. Truncated sinograms (input) were used with original slice images (target) in training to yield an improved reconstruction (output). DGAN performance was compared with the adaptive de-truncation method, total variation regularization, and two reference DL methods: FBPConvNet, and U-Net-based sinogram extension (ES-UNet). Our DGAN method and ES-UNet yielded the best root-mean-squared error (RMSE) (0.03 ± 0.01), and structural similarity index (SSIM) (0.92 ± 0.02) values, and reference DL methods also yielded good results. Furthermore, we performed an extended FOV analysis by increasing the reconstruction area by 10% and 20%. In both cases, the DGAN approach yielded best results at RMSE (0.03 ± 0.01 and 0.04 ± 0.01 for the 10% and 20% cases, respectively), peak signal-to-noise ratio (PSNR) (30.5 ± 2.6 dB and 28.6 ± 2.6 dB), and SSIM (0.90 ± 0.02 and 0.87 ± 0.02). In conclusion, our method was able to not only reconstruct the interior region with improved image quality, but also extend the reconstructed FOV by 20%.

Optimizing iterative reconstruction for quantification of calcium hydroxyapatite with photon counting flat-detector computed tomography: a cardiac phantom study

Purpose: Coronary artery calcium (CAC) scoring with computed tomography (CT) has been proposed as a screening tool for coronary artery disease, but concerns remain regarding the radiation dose of CT CAC scoring. Photon counting detectors and iterative reconstruction (IR) are promising approaches for patient dose reduction, yet the preservation of CAC scores with IR has been questioned. The purpose of this study was to investigate the applicability of IR for quantification of CAC using a photon counting flat-detector. Approach: We imaged a cardiac rod phantom with calcium hydroxyapatite (CaHA) inserts with different noise levels using an experimental photon counting flat-detector CT setup to simulate the clinical CAC scoring protocol. We applied filtered back projection (FBP) and two IR algorithms with different regularization strengths. We compared the air kerma values, image quality parameters [noise magnitude, noise power spectrum, modulation transfer function (MTF), and contrast-to-noise ratio], and CaHA quantification accuracy between FBP and IR. Results: IR regularization strength influenced CAC scores significantly ( p < 0.05 ). The CAC volumes and scores between FBP and IRs were the most similar when the IR regularization strength was chosen to match the MTF of the FBP reconstruction. Conclusion: When the regularization strength is selected to produce comparable spatial resolution with FBP, IR can yield comparable CAC scores and volumes with FBP. Nonetheless, at the lowest radiation dose setting, FBP produced more accurate CAC volumes and scores compared to IR, and no improved CAC scoring accuracy at low dose was demonstrated with the utilized IR methods.

Interior photon counting computed tomography for quantification of coronary artery calcium: pre-clinical phantom study

Computed tomography (CT) is the reference method for cardiac imaging, but concerns have been raised regarding the radiation dose of CT examinations. Recently, photon counting detectors (PCDs) and interior tomography, in which the radiation beam is limited to the organ-of-interest, have been suggested for patient dose reduction. In this study, we investigated interior PCD-CT (iPCD-CT) for non-enhanced quantification of coronary artery calcium (CAC) using an anthropomorphic torso phantom and ex vivo coronary artery samples. We reconstructed the iPCD-CT measurements with filtered back projection (FBP), iterative total variation (TV) regularization, padded FBP, and adaptively detruncated FBP and adaptively detruncated TV. We compared the organ doses between conventional CT and iPCD-CT geometries, assessed the truncation and cupping artifacts with iPCD-CT, and evaluated the CAC quantification performance of iPCD-CT. With approximately the same effective dose between conventional CT geometry (0.30 mSv) and interior PCD-CT with 10.2 cm field-of-view (0.27 mSv), the organ dose of the heart was increased by 52.3% with interior PCD-CT when compared to CT. Conversely, the organ doses to peripheral and radiosensitive organs, such as the stomach (55.0% reduction), were often reduced with interior PCD-CT. FBP and TV did not sufficiently reduce the truncation artifact, whereas padded FBP and adaptively detruncated FBP and TV yielded satisfactory truncation artifact reduction. Notably, the adaptive detruncation algorithm reduced truncation artifacts effectively when it was combined with reconstruction detrending. With this approach, the CAC quantification accuracy was good, and the coronary artery disease grade reclassification rate was particularly low (5.6%). Thus, our results confirm that CAC quantification can be performed with the interior CT geometry, that the artifacts are effectively reduced with suitable interior reconstruction methods, and that interior tomography provides efficient patient dose reduction.

Framework for Photon Counting Quantitative Material Decomposition

In this paper, the accuracy of material decomposition (MD) using an energy discriminating photon counting detector was studied. An MD framework was established and validated using calcium hydroxyapatite (CaHA) inserts of known densities (50 mg/cm³, 100 mg/cm³, 250 mg/cm³, 400 mg/cm³), and diameters (1.2, 3.0, and 5.0 mm). These inserts were placed in a cardiac rod phantom that mimics a tissue equivalent heart and measured using an experimental photon counting detector cone beam computed tomography (PCD-CBCT) setup. The quantitative coronary calcium scores (density, mass, and volume) obtained from the MD framework were compared with the nominal values. In addition, three different calibration techniques, signal-to-equivalent thickness calibration (STC), polynomial correction (PC), and projected equivalent thickness calibration (PETC) were compared to investigate the effect of the calibration method on the quantitative values. The obtained MD estimates agreed well with the nominal values for density (mass) with mean absolute percent errors (MAPEs) 8 ± 11% (9 ± 15%) and 4 ± 6% (9 ± 14%) for STC and PETC calibration methods, respectively. PC displayed large MAPEs for density (27 ± 9%), and mass (25 ± 12%). Volume estimation resulted in large deviations between true and measured values with notable MAPEs for STC (40 ± 90%), PC (40 ± 80%), and PETC (40 ± 90%). The framework demonstrated the feasibility of quantitative CaHA mass and density scoring using PCD-CBCT.