Iterative reconstruction for photon-counting CT using prior image constrained total generalized variation

Polyenergetic reconstruction mitigates streak artifacts by dual source imaging in chest photon counting detector computed tomography

OBJECTIVE

This study aims to assess the efficacy of polyenergetic reconstruction methods in reducing streak artifacts caused by dual source imaging in Photon Counting Detector Computed Tomography (PCD-CT) imaging, thereby improving image quality and diagnostic accuracy.

METHODS

A retrospective cohort study was conducted, involving 50 patients who underwent chest Computed Tomography Angiography with PCD-CT, focusing on those with streak artifacts. Quantitative and qualitative analyses were performed on images reconstructed using monoenergetic and polyenergetic techniques. Quantitative evaluations measured the attenuation of tracheal air density in regions affected by streak artifacts, while qualitative assessments employed a modified Likert scale to rate image quality. Statistical analyses included Wilcoxon's signed-rank tests and Spearman's correlation, alongside assessments of inter-rater reliability.

RESULTS

There was significantly lower attenuation of tracheal air density on the polyenergetic reconstructions (Median - 1010 ± 62 HU vs -930 ± 110 HU; P < 0.001), and significantly decreased variation on the polyenergetic reconstructions (Median 65.2 ± 79.5 HU vs 38.8 ± 33.9 HU; P < 0.001). The median modified-Likert scale were significantly better for the polyenergetic reconstructions (median modified-Likert 4 ± 0.5 vs 2.5 ± 1; P < 0.001). The inter-rater agreement was substantial and not significantly different between reconstructions (Gwet's ACPolyenergetic = 0.78 vs Gwet's ACVMI = 0.775).

CONCLUSION

Polyenergetic reconstruction significantly mitigates streak artifacts in PCD-CT imaging, enhancing quantitative and qualitative image quality. This advancement addresses a known limitation of current PCD-CT reconstruction techniques, offering a promising approach to improving diagnostic reliability and accuracy in clinical practice. We demonstrate that future software implementations can resolve this artifact.

Multi-energy CT material decomposition using graph model improved CNN

In spectral CT imaging, the coefficient image of the basis material obtained by the material decomposition technique can estimate the tissue composition, and its accuracy directly affects the disease diagnosis. Although the precision of material decomposition is increased by employing convolutional neural networks (CNN), extracting the non-local features from the CT image is restricted using the traditional CNN convolution operator. A graph model built by multi-scale non-local self-similar patterns is introduced into multi-material decomposition (MMD). We proposed a novel MMD method based on graph edge-conditioned convolution U-net (GECCU-net) to enhance material image quality. The GECCU-net focuses on developing a multi-scale encoder. At the network coding stage, three paths are applied to capture comprehensive image features. The local and non-local feature aggregation (LNFA) blocks are designed to integrate the local and non-local features from different paths. The graph edge-conditioned convolution based on non-Euclidean space excavates the non-local features. A hybrid loss function is defined to accommodate multi-scale input images and avoid over-smoothing of results. The proposed network is compared quantitatively with base CNN models on the simulated and real datasets. The material images generated by GECCU-net have less noise and artifacts while retaining more information on tissue. The Structural SIMilarity (SSIM) of the obtained abdomen and chest water maps reaches 0.9976 and 0.9990, respectively, and the RMSE reduces to 0.1218 and 0.4903 g/cm3. The proposed method can improve MMD performance and has potential applications.

Adaptive prior image constrained total generalized variation for low-dose dynamic cerebral perfusion CT reconstruction

BACKGROUND

Dynamic cerebral perfusion CT (DCPCT) can provide valuable insight into cerebral hemodynamics by visualizing changes in blood within the brain. However, the associated high radiation dose of the standard DCPCT scanning protocol has been a great concern for the patient and radiation physics. Minimizing the x-ray exposure to patients has been a major effort in the DCPCT examination. A simple and cost-effective approach to achieve low-dose DCPCT imaging is to lower the x-ray tube current in data acquisition. However, the image quality of low-dose DCPCT will be degraded because of the excessive quantum noise.

OBJECTIVE

To obtain high-quality DCPCT images, we present a statistical iterative reconstruction (SIR) algorithm based on penalized weighted least squares (PWLS) using adaptive prior image constrained total generalized variation (APICTGV) regularization (PWLS-APICTGV).

METHODS

APICTGV regularization uses the precontrast scanned high-quality CT image as an adaptive structural prior for low-dose PWLS reconstruction. Thus, the image quality of low-dose DCPCT is improved while essential features of targe image are well preserved. An alternating optimization algorithm is developed to solve the cost function of the PWLS-APICTGV reconstruction.

RESULTS

PWLS-APICTGV algorithm was evaluated using a digital brain perfusion phantom and patient data. Compared to other competing algorithms, the PWLS-APICTGV algorithm shows better noise reduction and structural details preservation. Furthermore, the PWLS-APICTGV algorithm can generate more accurate cerebral blood flow (CBF) map than that of other reconstruction methods.

CONCLUSIONS

PWLS-APICTGV algorithm can significantly suppress noise while preserving the important features of the reconstructed DCPCT image, thus achieving a great improvement in low-dose DCPCT imaging.

A Deep Learning Approach for Rapid and Generalizable Denoising of Photon-Counting Micro-CT Images

Photon-counting CT (PCCT) is powerful for spectral imaging and material decomposition but produces noisy weighted filtered backprojection (wFBP) reconstructions. Although iterative reconstruction effectively denoises these images, it requires extensive computation time. To overcome this limitation, we propose a deep learning (DL) model, UnetU, which quickly estimates iterative reconstruction from wFBP. Utilizing a 2D U-net convolutional neural network (CNN) with a custom loss function and transformation of wFBP, UnetU promotes accurate material decomposition across various photon-counting detector (PCD) energy threshold settings. UnetU outperformed multi-energy non-local means (ME NLM) and a conventional denoising CNN called UnetwFBP in terms of root mean square error (RMSE) in test set reconstructions and their respective matrix inversion material decompositions. Qualitative results in reconstruction and material decomposition domains revealed that UnetU is the best approximation of iterative reconstruction. In reconstructions with varying undersampling factors from a high dose ex vivo scan, UnetU consistently gave higher structural similarity (SSIM) and peak signal-to-noise ratio (PSNR) to the fully sampled iterative reconstruction than ME NLM and UnetwFBP. This research demonstrates UnetU's potential as a fast (i.e., 15 times faster than iterative reconstruction) and generalizable approach for PCCT denoising, holding promise for advancing preclinical PCCT research.

Low-dose spectral reconstruction with global, local, and nonlocal priors based on subspace decomposition

Background

Multienergy computed tomography (MECT) is a promising imaging modality for material decomposition, lesion detection, and other clinical applications. However, there is an urgent need to design efficient and accurate algorithms to solve the inverse problems related to spectral reconstruction and improve image quality, especially under low-dose and incomplete datasets. The key issue for MECT reconstruction is how to efficiently describe the interchannel and intrachannel priors in multichannel images.

Methods

In this model, in order to correlate the similarities of interchannel images and regularize the multichannel images, the global, local, and nonlocal priors are jointly integrated into the low-dose MECT reconstruction model. First, the subspace decomposition method first employs the global low-rankness to map the original MECT images to the low-dimensional eigenimages. Then, nonlocal self-similarity of the eigenimages is cascaded into the optimization model. Additionally, the L0 quasi-norm on gradient images is incorporated into the proposed method to further enhance the local sparsity of intrachannel images. The alternating direction method is applied to solve the optimization model in an iterative scheme.

Results

Simulation, preclinical, and real datasets were applied to validate the effectiveness of the proposed method. From the simulation dataset, the new method was found to reduce the root-mean-square error (RMSE) by 42.31% compared with the latest research fourth-order nonlocal tensor decomposition MECT reconstruction (FONT-SIR) method under 160 projection views. The calculation time of an iteration for the proposed method was 23.07% of the FONT-SIR method. The results of material decomposition in real mouse data further confirmed the accuracy of the proposed method for different materials.

Conclusions

We developed a method in which the global, local, and nonlocal priors are jointly used to develop the reconstruction model for low-dose MECT, where the global low-rankness and nonlocal prior are cascaded by subspace decomposition and block-matching, and the L0 sparsity is applied to express the local prior. The results of the experiments demonstrate that the proposed method based on subspace improves computational efficiency and has advantages in noise suppression and structure preservation over competing algorithms.

A multienergy computed tomography method without image segmentation or prior knowledge of X-ray spectra or materials

Many methods have been proposed for multienergy computed tomography (CT) imaging based on traditional CT systems. Usually, either prior knowledge of the X-ray spectra distribution or materials or the segmentation of the projection or reconstructed image is needed. To avoid these requirements, a multienergy CT method is proposed in this paper. A CT image can be seen as a linear combination of energy-dependent components and spatially dependent components. The latter components are the base images, while the former components are the coefficients. A blind decomposition model is constructed to decompose the multivoltage projections to obtain the base images and the energies. Multienergy CT images are computationally synthesized with the base images and the energies. Multivoltage projections can be acquired based on one scan with stepped voltages. X-ray scattering is considered an important factor in imaging errors and appears as a low-frequency signal. The variance is used to describe the low-frequency features and is minimized as the optimized objective function of the decomposition model. The solution of the model uses Karush-Kuhn-Tucker (KKT) conditions. In the experiments, the images reconstructed with the proposed method exhibit weak beam-hardening artifacts. Additionally, the X-ray energies of the different materials represented have small relative errors. Therefore, the reconstructed images have narrow energy intervals. This shows the effectiveness of the proposed method.

Performance of Spectral Photon-Counting Coronary CT Angiography and Comparison with Energy-Integrating-Detector CT: Objective Assessment with Model Observer

AIMS

To evaluate spectral photon-counting CT's (SPCCT) objective image quality characteristics in vitro, compared with standard-of-care energy-integrating-detector (EID) CT.

METHODS

We scanned a thorax phantom with a coronary artery module at 10 mGy on a prototype SPCCT and a clinical dual-layer EID-CT under various conditions of simulated patient size (small, medium, and large). We used filtered back-projection with a soft-tissue kernel. We assessed noise and contrast-dependent spatial resolution with noise power spectra (NPS) and target transfer functions (TTF), respectively. Detectability indices (d') of simulated non-calcified and lipid-rich atherosclerotic plaques were computed using the non-pre-whitening with eye filter model observer.

RESULTS

SPCCT provided lower noise magnitude (9-38% lower NPS amplitude) and higher noise frequency peaks (sharper noise texture). Furthermore, SPCCT provided consistently higher spatial resolution (30-33% better TTF10). In the detectability analysis, SPCCT outperformed EID-CT in all investigated conditions, providing superior d'. SPCCT reached almost perfect detectability (AUC ≈ 95%) for simulated 0.5-mm-thick non-calcified plaques (for large-sized patients), whereas EID-CT had lower d' (AUC ≈ 75%). For lipid-rich atherosclerotic plaques, SPCCT achieved 85% AUC vs. 77.5% with EID-CT.

CONCLUSIONS

SPCCT outperformed EID-CT in detecting simulated coronary atherosclerosis and might enhance diagnostic accuracy by providing lower noise magnitude, markedly improved spatial resolution, and superior lipid core detectability.

Computed tomography with a full FOV photon-counting detector in a clinical setting, the first experience

OBJECTIVES

to assess the feasibility of CT with an integrated photon-counting-detector system (PC-CT) in the body imaging of clinical patients.

METHODS

120 examinations using photon counting detector CT were evaluated in six groups: 1/ a standard-dose lung, 2/ low-dose lung, 3/ ultra-high resolution (UHR) lung, 4/ standard-dose abdominal, 5/ lower-dose abdominal, 6/ UHR abdominal CTA. All CT examinations were performed on a single-source prototype device equipped with a photon counting detector covering a 50 cm scan field of view. Standard dose examinations were performed with the use of detector element size of 0.4 mm, ultra-high-resolution examinations with the detector element size of 0.2 mm, respectively. The stability of the system during imaging was tested. The diagnostic quality of the acquired images was assessed based on the imaging of key structures and the noise level in five-point scale, the effective dose equivalent, dose length product and noise level, and also relation to body mass index and body surface area were compared with three similar groups of CT images made with energy integrating high end scanner. The parameters were evaluated using Wilcoxon test for independent samples, the independence was tested using Kruskal-Wallis test.

RESULTS

When PC-CT images radiation dose is compared with the similar imaging using energy integrating CT, the PC-CT shows lower dose in ultra-high resolution mode, the dose is significantly lower (p < 0.0001), the standard dose examinations were performed with the comparable radiation doses. PC-CT exhibited the significantly higher ratio between parenchyma signal and background noise both in lung and in abdominal imaging (p < 0.0001).

CONCLUSIONS

PC-CT showed imaging stability and excellent diagnostic quality at dose values that are comparable or better to the quality of energy integrating CT, the better signal and improved resolution is most important advantage of photon counting detector CT over energy integrating detector CT.

Photon-counting x-ray detectors for CT

The introduction of photon-counting detectors is expected to be the next major breakthrough in clinical x-ray computed tomography (CT). During the last decade, there has been considerable research activity in the field of photon-counting CT, in terms of both hardware development and theoretical understanding of the factors affecting image quality. In this article, we review the recent progress in this field with the intent of highlighting the relationship between detector design considerations and the resulting image quality. We discuss detector design choices such as converter material, pixel size, and readout electronics design, and then elucidate their impact on detector performance in terms of dose efficiency, spatial resolution, and energy resolution. Furthermore, we give an overview of data processing, reconstruction methods and metrics of imaging performance; outline clinical applications; and discuss potential future developments.

Statistical image-based material decomposition for triple-energy computed tomography using total variation regularization

BACKGROUND

Triple-energy computed tomography (TECT) can obtain x-ray attenuation measurements at three energy spectra, thereby allowing identification of different material compositions with same or very similar attenuation coefficients. This ability is known as material decomposition, which can decompose TECT images into different basis material image. However, the basis material image would be severely degraded when material decomposition is directly performed on the noisy TECT measurements using a matrix inversion method.

OBJECTIVE

To achieve high quality basis material image, we present a statistical image-based material decomposition method for TECT, which uses the penalized weighted least-squares (PWLS) criteria with total variation (TV) regularization (PWLS-TV).

METHODS

The weighted least-squares term involves the noise statistical properties of the material decomposition process, and the TV regularization penalizes differences between local neighboring pixels in a decomposed image, thereby contributing to improving the quality of the basis material image. Subsequently, an alternating optimization method is used to minimize the objective function.

RESULTS

The performance of PWLS-TV is quantitatively evaluated using digital and mouse thorax phantoms. The experimental results show that PWLS-TV material decomposition method can greatly improve the quality of decomposed basis material image compared to the quality of images obtained using the competing methods in terms of suppressing noise and preserving edge and fine structure details.

CONCLUSIONS

The PWLS-TV method can simultaneously perform noise reduction and material decomposition in one iterative step, and it results in a considerable improvement of basis material image quality.