Image Reconstruction In Spectral CT Research Articles

Anisotropic sparse transformation for spectral CT image reconstruction

AbstractPhoton‐counting detector‐based computed tomography (PCD‐CT) is an advanced realization of spectral CT, the multi‐energy projection data is captured from the same object, hence, CT images can provide additional spectral resolution, making it possible to perform material decomposition. However, multiple projections may have a low signal‐to‐noise ratio (SNR), such that CT images suffer from noise. To handle this problem, a spectral CT image reconstruction method based on anisotropic sparse transformation (AST) is proposed. To increase the quality of reconstruction, AST through an anisotropic guided filter (AGF) and quasi norm is proposed. Then, as a new regularization, AST is introduced into an iterative reconstruction process, generating an AST‐based method. Moreover, to utilize the correlation among projection data, the average image serves as the guidance image of AGF, it varies with the iterative index, resulting in a technique of dynamic average image (DAI). The AST‐based model involves quasi norm minimization, hence an effective strategy is employed to solve the corresponding problem. A series of experiments were performed. The experiment showed that, compared to other listed methods, the result of the AST‐based method can achieve better visualization and higher quantitative indexes, hence it has application potential in the medical imaging field.

Nonlocal Tensor Decomposition With Joint Low Rankness and Smoothness for Spectral CT Image Reconstruction

Nonlocal Tensor Decomposition With Joint Low Rankness and Smoothness for Spectral CT Image Reconstruction

Spectral CT image reconstruction using a constrained optimization approach-An algorithm for AAPM 2022 spectral CT grand challenge and beyond.

CT reconstruction is of essential importance in medical imaging. In 2022, the American Association of Physicists in Medicine (AAPM) sponsored a Grand Challenge to investigate the challenging inverse problem of spectral CT reconstruction, with the aim of achieving the most accurate reconstruction results. The authors of this paper participated in the challenge and won as a runner-upteam. This paper reports details of our PROSPECT algorithm (Prior-based Restricted-variable Optimization for SPEctral CT) and follow-up studies regarding the algorithm's accuracy and enhancement of its convergencespeed. We formulated the reconstruction task as an optimization problem. PROSPECT employed a one-step backward iterative scheme to solve this optimization problem by allowing estimation of and correction for the difference between the actual polychromatic projection model and the monochromatic model used in the optimization problem. PROSPECT incorporated various forms of prior information derived by analyzing training data provided by the Grand Challenge to reduce the number of unknown variables. We investigated the impact of projection data precision on the resulting solution accuracy and improved convergence speed of the PROSPECT algorithm by incorporating a beam-hardening correction (BHC) step in the iterative process. We also studied the algorithm's performance under noisy projectiondata. Prior knowledge allowed a reduction of the number of unknown variables by . PROSPECT algorithm achieved the average root of mean square error (RMSE) of in the test data set provided by the Grand Challenge. Performing the reconstruction with the same algorithm but using double-precision projection data reduced RMSE to . Including the BHC step in the PROSPECT algorithm accelerated the iteration process with a 40% reduction in computationtime. PROSPECT algorithm achieved a high degree of accuracy and computationalefficiency.

Dual energy CT reconstruction using the constrained one step spectral image reconstruction algorithm.

The constrained one-step spectral CT Image Reconstruction method (cOSSCIR) has been developed to estimate basis material maps directly from spectral CT data using a model of the polyenergetic x-ray transmissions and incorporating convex constraints into the inversion problem. This 'one-step' approach has been shown to stabilize the inversion in the case of photon-counting CT, and may provide similar benefits to dual-kV systems that utilize integrating detectors. Since the approach does not require the same rays be acquired for every spectral measurement, cOSSCIR can apply to dual energy protocols and systems used clinically, such as fast and slow kV switching systems and dual sourcescanning. The purpose of this study is to investigate the use of cOSSCIR applied to dual-kV data, using both registered and unregistered spectral acquisitions, specifically slow and fast kV switching imaging protocols. For this application, cOSSCIR is investigated using inverse crime simulations and dual-kV experiments. This study is the first demonstration of cOSSCIR on the dual-kV reconstructionproblem. An integrating detector model was developed for the purpose of reconstructing dual-kV data, and an inverse crime study was used to validate the detector model within the cOSSCIR framework using a simulated pelvic phantom. Experiments were also used to evaluate cOSSCIR on the dual energy problem. Dual-kV data was obtained from a physical phantom containing analogs of adipose, bone, and liver tissues, with the aim of recovering the material coefficients in the bone and adipose basis material maps. cOSSCIR was applied to acquisitions where all rays performed both spectral measurements (registered) and fast and slow kV switching acquisitions (unregistered). cOSSCIR was also compared to two image-domain decomposition approaches, where image-domain methods are the conventional approach for decomposing unregistered spectraldata. Simulations demonstrate the application of cOSSCIR to the dual-kV inversion problem by successfully recovering the material basis maps on ideal data, while further showing that unregistered data presents a more challenging inversion problem. In our experimental reconstructions, the recovered basis material coefficient errors were found to be less than 6.5% in the bone, adipose, and liver regions for both registered and unregistered protocols. Similarly, the errors were less than 4% in the 50 keV virtual mono-energetic images, and the recovered material decomposition vectors nearly overlap their corresponding ground-truth vectors. Additionally, a preliminary two material decomposition study of iodine quantification recovered an average concentration of 9.2 mg/mL from a 10 mg/mL experimental iodineanalog. Using our integrating detector and spectral models, cOSCCIR is capable of accurately recovering material basis maps from dual-kV data for both registered and unregistered data. The material decomposition quantification compare favorably to the image domain approaches, and our results were not affected by the imaging protocol. Our results also suggest the extension of cOSSCIR to iodine quantification using two materialdecomposition.

Constrained one-step material decomposition reconstruction of head CT data from a silicon photon-counting prototype.

Spectral CT material decomposition provides quantitative information but is challenged by the instability of the inversion into basis materials. We have previously proposed the constrained One-Step Spectral CT Image Reconstruction (cOSSCIR) algorithm to stabilize the material decomposition inversion by directly estimating basis material images from spectral CT data. cOSSCIR was previously investigated on phantomdata. This study investigates the performance of cOSSCIR using head CT datasets acquired on a clinical photon-counting CT (PCCT) prototype. This is the first investigation of cOSSCIR for large-scale, anatomically complex, clinical PCCT data. The cOSSCIR decomposition is preceded by a spectrum estimation and nonlinear counts correction calibration step to address nonideal detectoreffects. Head CT data were acquired on an early prototype clinical PCCT system using an edge-on silicon detector with eight energy bins. Calibration data of a step wedge phantom were also acquired and used to train a spectral model to account for the source spectrum and detector spectral response, and also to train a nonlinear counts correction model to account for pulse pileup effects. The cOSSCIR algorithm optimized the bone and adipose basis images directly from the photon counts data, while placing a grouped total variation (TV) constraint on the basis images. For comparison, basis images were also reconstructed by a two-step projection-domain approach of Maximum Likelihood Estimation (MLE) for decomposing basis sinograms, followed by filtered backprojection (MLE + FBP) or a TV minimization algorithm (MLE + TVmin ) to reconstruct basis images. We hypothesize that the cOSSCIR approach will provide a more stable inversion into basis images compared to two-step approaches. To investigate this hypothesis, the noise standard deviation in bone and soft-tissue regions of interest (ROIs) in the reconstructed images were compared between cOSSCIR and the two-step methods for a range of regularization constraintsettings. cOSSCIR reduced the noise standard deviation in the basis images by a factor of two to six compared to that of MLE + TVmin , when both algorithms were constrained to produce images with the same TV. The cOSSCIR images demonstrated qualitatively improved spatial resolution and depiction of fine anatomical detail. The MLE + TVmin algorithm resulted in lower noise standard deviation than cOSSCIR for the virtual monoenergetic images (VMIs) at higher energy levels and constraint settings, while the cOSSCIR VMIs resulted in lower noise standard deviation at lower energy levels and overall higher qualitative spatial resolution. There were no statistically significant differences in the mean values within the bone region of images reconstructed by the studied algorithms. There were statistically significant differences in the mean values within the soft-tissue region of the reconstructed images, with cOSSCIR producing mean values closer to the expectedvalues. The cOSSCIR algorithm, combined with our previously proposed spectral model estimation and nonlinear counts correction method, successfully estimated bone and adipose basis images from high resolution, large-scale patient data from a clinical PCCT prototype. The cOSSCIR basis images were able to depict fine anatomical details with a factor of two to six reduction in noise standard deviation compared to that of the MLE + TVmin two-stepapproach.

Non-convex primal-dual algorithm for image reconstruction in spectral CT

The work seeks to develop an algorithm for image reconstruction by directly inverting the non-linear data model in spectral CT. Using the non-linear data model, we formulate the image-reconstruction problem as a non-convex optimization program, and develop a non-convex primal-dual (NCPD) algorithm to solve the program. We devise multiple convergence conditions and perform verification studies numerically to demonstrate that the NCPD algorithm can solve the non-convex optimization program and under appropriate data condition, can invert the non-linear data model. Using the NCPD algorithm, we then reconstruct monochromatic images from simulated and real data of numerical and physical phantoms acquired with a standard, full-scan dual-energy configuration. The result of the reconstruction studies shows that the NCPD algorithm can correct accurately for the non-linear beam-hardening effect. Furthermore, we apply the NCPD algorithm to simulated and real data of the numerical and physical phantoms collected with non-standard, short-scan dual-energy configurations, and obtain monochromatic images comparable to those of the standard, full-scan study, thus revealing the potential of the NCPD algorithm for enabling non-standard scanning configurations in spectral CT, where the existing indirect methods are limited.

Iterative spectral CT reconstruction based on low rank and average-image-incorporated BM3D

In a photon counting detector integrated spectral CT scanner, the received photons are counted in several energy channels to generate the corresponding projections. Since the projection in each energy channel is generated using part of the received photons, the reconstructed channel image suffers from severe noise. Therefore, image reconstruction in spectral CT is considered to be a big challenge. Because the inter-channel images are all from the same object but in different energy bins, there exists a strong correlation among these images. Moreover, it is suggested that there are similarities among various patches of CT images in the spatial domain. In this work, we propose average-image-incorporated block-matching and 3D (aiiBM3D) filtering along with low rank regularization for iterative spectral CT reconstruction. The aiiBM3D method is based on filtered 3D data arrays formed by similar 2D blocks using the mapped version of the average image obtained from linear regression. The reconstruction procedure consists of two main steps. First, the alternating direction method of multipliers is employed to solve the problem with low rank regularization where the goal is to exploit the correlation in inter-channel images. Second, our proposed BM3D-based algorithm is applied to all the channel images to make use of the redundant information in the spatial domain and inter-channel. The two steps repeat until the stopping criteria are satisfied. The proposed method is validated on numerically simulated and preclinical datasets. Our results confirm its high performance in terms of signal to noise ratio and structural preservation.

MO‐FG‐204‐07: Nonlocal Total Variation Based Spectral CT Image Reconstruction

Purpose:When divided into narrow energy bins, the projection data for spectral CT contains spectral information with reduced signal‐to‐noise ratio (SNR). When used together, the total projection data has improved SNR without spectral information. The purpose of this work is to develop an image reconstruction method that maximizes the utility of spectral data in terms of synergized spectral information and SNR.Methods:The proposed method is based on nonlocal total variation (TV) (NLTV). That is, we first reconstruct a spectrally averaged image using total projection data without energy binning, and build the NLTV weights from this image that characterize the nonlocal image feature with good SNR, such as texture and fine structures. Then these NLTV weights are incorporated into a NLTV‐based iterative spectral image reconstruction scheme using the binned projection data with reduced SNR, so that the weights serve as the reference to regularize the image quality during spectral image reconstruction. In addition, the spectral reconstruction incorporates the material decomposition by directly reconstructing material compositions instead of the conventional approach, that separates image reconstruction and material decomposition. The spectral image reconstruction is formulated as an iterative reconstruction method with the NLTV regularization. The split Bergman method is developed for solving the above sparse optimization problem.Results:The proposed NLTV‐based method was validated using both simulated and experimental data. In comparison with FBP and TV‐ based method, NLTV had improved image quality, e.g., contrast‐to‐noise ratio (CNR) and modulation transfer function (MTF).Conclusion:NLTV‐based spectral CT image reconstruction is developed that maximizes the utility of spectral data in terms of synergized spectral information, and it offers significant improvement in CNR without scarifying spatial resolution of the image.Jiulong Liu and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000) and the Shanghai Pujiang Talent Program (#14PJ1404500).