Accurate Iterative FBP Reconstruction Method for Material Decomposition of Dual Energy CT

Multi-energy blended CBCT spectral imaging and scatter-decoupled material decomposition using a spectral modulator with flying focal spot (SMFFS)

BACKGROUND

Cone-beam CT (CBCT) has been extensively employed in industrial and medical applications, such as image-guided radiotherapy and diagnostic imaging, with a growing demand for quantitative imaging using CBCT. However, conventional CBCT can be easily compromised by scatter and beam hardening artifacts, and the entanglement of scatter and spectral effects introduces additional complexity.

PURPOSE

The intertwined scatter and spectral effects within CBCT pose significant challenges to the quantitative performance of spectral imaging. In this work, we present the first attempt to develop a stationary spectral modulator with flying focal spot (SMFFS) technology as a promising, low-cost approach to accurately solving the x-ray scattering problem and physically enabling spectral imaging in a unified framework, and with no significant misalignment in data sampling of spectral projections.

METHODS

To deal with the intertwined scatter-spectral challenge, we propose a novel scatter-decoupled material decomposition (SDMD) method for SMFFS, which consists of four steps in total, including (1) spatial resolution-preserved and noise-suppressed multi-energy "residual" projection generation free from scatter, based on a hypothesis of scatter similarity; (2) first-pass material decomposition from the generated multi-energy residual projections in non-penumbra regions, with a structure similarity constraint to overcome the increased noise and penumbra effect; (3) scatter estimation for complete data; and (4) second-pass material decomposition for complete data by using a multi-material spectral correction method. Monte Carlo simulations of a pure-water cylinder phantom with different focal spot deflections are conducted to validate the scatter similarity hypothesis. Both numerical simulations using a clinical abdominal CT dataset, and physics experiments on a tabletop CBCT system using a Gammex multi-energy CT phantom and an anthropomorphic chest phantom, are carried out to demonstrate the feasibility of CBCT spectral imaging with SMFFS and our proposed SDMD method.

RESULTS

Monte Carlo simulations show that focal spot deflections within a range of 2 mm share quite similar scatter distributions overall. Numerical simulations demonstrate that SMFFS with SDMD method can achieve better material decomposition and CT number accuracy with fewer artifacts. In physics experiments, for the Gammex phantom, the average error of the mean values (E RMSE ROI $E^{\text{ROI}}_{\text{RMSE}}$ ) in selected regions of interest (ROIs) of virtual monochromatic image (VMI) at 70 keV is 8 HU in SMFFS cone-beam (CB) scan, and 19 and 210 HU in sequential 80/120 kVp (dual kVp, DKV) CB scan with and without scatter correction, respectively. For the chest phantom, theE RMSE ROI $E^{\text{ROI}}_{\text{RMSE}}$ in selected ROIs of VMIs is 12 HU for SMFFS CB scan, and 15 and 438 HU for sequential 80/140 kVp CB scan with and without scatter correction, respectively. Also, the non-uniformity among selected regions of the chest phantom is 14 HU for SMFFS CB scan, and 59 and 184 HU for the DKV CB scan with and without a traditional scatter correction method, respectively.

CONCLUSIONS

We propose a SDMD method for CBCT with SMFFS. Our preliminary results show that SMFFS can enable spectral imaging with simultaneous scatter correction for CBCT and effectively improve its quantitative imaging performance.

A model-based direct inversion network (MDIN) for dual spectral computed tomography

Objective. Dual spectral computed tomography (DSCT) is a very challenging problem in the field of imaging. Due to the nonlinearity of its mathematical model, the images reconstructed by the conventional CT usually suffer from the beam hardening artifacts. Additionally, several existing DSCT methods rely heavily on the information of the spectra, which is often not readily available in applications. To address this problem, in this study, we aim to develop a novel approach to improve the DSCT reconstruction performance.Approach. A model-based direct inversion network (MDIN) is proposed for DSCT, which can directly predict the basis material images from the collected polychromatic projections. The all operations are performed in the network, requiring neither the conventional algorithms nor the information of the spectra. It can be viewed as an approximation to the inverse procedure of DSCT imaging model. The MDIN is composed of projection pre-decomposition module (PD-module), domain transformation layer (DT-layer), and image post-decomposition module (ID-module). The PD-module first performs the pre-decomposition on the polychromatic projections that consists of a series of stacked one-dimensional convolution layers. The DT-layer is designed to obtain the preliminary decomposed results, which has the characteristics of sparsely connected and learnable parameters. And the ID-module uses a deep neural network to further decompose the reconstructed results of the DT-layer so as to achieve higher-quality basis material images.Main results. Numerical experiments demonstrate that the proposed MDIN has significant advantages in substance decomposition, artifact reduction and noise suppression compared to other methods in the DSCT reconstruction.Significance. The proposed method has a flexible applicability, which can be extended to other CT problems, such as multi-spectral CT and low dose CT.

A Locally Weighted Linear Regression Look-Up Table-Based Iterative Reconstruction Method for Dual Spectral CT

OBJECTIVE

Compared with traditional computed tomography (CT), dual spectral CT (DSCT) exhibits superior material distinguishability and thus has broad prospects in industrial and medical fields. In iterative DSCT algorithms, accurately modeling forward-projection functions is crucial, but it is very difficult to analytically provide accurate functions.

METHODS

In this article, we propose a locally weighted linear regression look-up table-based (LWLR-LUT) iterative reconstruction method for DSCT. First, the proposed method uses LWLR to establish LUTs for the forward-projection functions through calibration phantoms, achieving good local information calibration. Second, the reconstructed images can be iteratively obtained through the established LUTs. The proposed method not only does not require knowledge of the X-ray spectra and the attenuation coefficients, but also implicitly accounts for some scattered radiation while fitting locally the forward-projection functions in the calibration space.

RESULTS

Both numerical simulations and real data experiments demonstrate that the proposed method can achieve highly accurate polychromatic forward-projection functions and greatly improve the quality of the images reconstructed from scattering-free and scattering projections.

CONCLUSION

The proposed method is simple and practical, and achieves good material decomposition effects for objects with different complex structures through simple calibration phantoms.

Mineral quantitative characterization method based on basis material decomposition model by dual-energy computed tomography

BACKGROUND

Dual-energy computed tomography (DECT) can reconstruct electron density ρe and effective atomic number Zeff distribution for material discrimination. Image-domain basis material decomposition (IBMD) method is a widely used DECT method. However, IBMD method cannot be used for mineral identification directly due to limitations of complex basis material determination, beam hardening artifacts, and inherent errors caused by approximate empirical formulas.

OBJECTIVE

This study proposes an improved IBMD (IIBMD) method to overcome the above limitations.

METHODS

In IIBMD method, the composition of basis material is optimized to obtain accurate decomposition coefficients, which enables accurate ρe and Zeff distribution. Moreover, the thickness of basis material is optimized to reduce the effect of beam hardening. Furthermore, two formulas in place of empirical formulas are proposed to calculate ρe and Zeff. Finally, a threshold technique is applied to separate different mineral phases.

RESULTS

Numerical simulations and practical experiments using a photon-counting detector CT system are implemented to verify IIBMD method. Results show that the relative errors of ρe and Zeff for seven common minerals are down to 5%, lower than most of the existing DECT methods for rocks. Reasonable volume fraction results of mineral phases are thus obtained through threshold segmentation.

CONCLUSIONS

This study demonstrates that the proposed IIBMD method has high practical value in mineralogical identification.

An iterative reconstruction method based on monochromatic images for dual energy CT

PURPOSE

Dual-energy computed tomography (DECT) scans objects using two different X-ray spectra to acquire more information, which is also called dual spectral CT (DSCT) in some articles. Compared to traditional CT, DECT exhibits superior material distinguishability. Therefore, DECT can be widely used in the medical and industrial domains. However, owing to the nonlinearity and ill condition of DECT, studies are underway on DECT reconstruction to obtain high quality images and achieve fast convergence speed. Therefore, in this study, we propose an iterative reconstruction method based on monochromatic images (IRM-MI) to rapidly obtain high-quality images in DECT reconstruction.

METHODS

An IRM-MI is proposed for DECT. The proposed method converts DECT reconstruction problem from the basis material images decomposition to monochromatic images decomposition to significantly improve the convergence speed of DECT reconstruction by changing the coefficient matrix of the original equations to increase the angle of the high- and low-energy projection curves or reduce the condition number of the coefficient matrix. The monochromatic images were then decomposed into basis material images. Furthermore, we conducted numerical experiments to evaluate the performance of the proposed method.

RESULTS

The decomposition results of the simulated data and real data experiments confirmed the effectiveness of the proposed method. Compared to the extended algebraic reconstruction technique (E-ART) method, the proposed method exhibited a significant increase in the convergence speed by increasing the angle of polychromatic projection curves or decreasing the condition number of the coefficient matrix, when choosing the appropriate monochromatic images. Therefore, the proposed method is also advantageous in acquiring high quality and rapidly converged images.

CONCLUSIONS

We developed an iterative reconstruction method based on monochromatic images for the material decomposition for DECT. The numerical experiments using the proposed method validated its capability of decomposing the basis material images. Furthermore, the proposed method achieved faster convergence speed compared to the E-ART method.

Dual-Energy Computed Tomography for the Diagnosis of Mediastinal Lymph Node Metastasis in Lung Cancer Patients: A Preliminary Study

OBJECTIVE

This study explored the feasibility of dual-energy computed tomography (DECT) for the diagnosis of mediastinal lymph node (LN) metastasis in patients with lung cancer.

METHODS

Forty-two consecutive patients with lung cancer, who underwent DECT, were included in this retrospective study. The attenuation value (Hounsfield unit) in virtual monochromatic images and the iodine concentration in the iodine map were measured at mediastinal LNs. The slope of the spectral attenuation curve (K) and normalized iodine concentration (in thoracic aorta) were calculated. The measurement results were statistically compared using 2 independent samples t test. Receiver operating characteristic curve analysis, net reclassification improvement, and integrated discrimination improvement were used to evaluate the diagnostic performance of DECT for mediastinal LN metastasis.

RESULTS

A total of 74 mediastinal LNs were obtained, including 33 metastatic LNs and 41 nonmetastatic LNs. The attenuation value at the lower energy levels of virtual monochromatic images (40-90 keV), K, and normalized iodine concentration demonstrated a significant difference between metastatic LNs and nonmetastatic LNs. The attenuation value at 40 keV was the most favorable biomarker for the diagnosis of mediastinal LN metastasis (area under curve, 0.91; sensitivity, 0.94; specificity, 0.81), which showed a much better performance than the LN diameter-based evaluation method (area under curve, 0.72; sensitivity, 0.66; specificity, 0.82; net reclassification improvement, 0.359; integrated discrimination improvement, 0.330).

CONCLUSIONS

Dual-energy computed tomography is a promising diagnostic approach for the diagnosis of mediastinal LN metastasis in patients with lung cancer, which may help clinicians implement personalized treatment strategies.

A material decomposition method for dual-energy CT via dual interactive Wasserstein generative adversarial networks

PURPOSE

Dual-energy computed tomography (DECT) is highly promising for material characterization and identification, whereas reconstructed material-specific images are affected by magnified noise and beam-hardening artifacts. Although various DECT material decomposition methods have been proposed to solve this problem, the quality of the decomposed images is still unsatisfactory, particularly in the image edges. In this study, a data-driven approach using dual interactive Wasserstein generative adversarial networks (DIWGAN) is developed to improve DECT decomposition accuracy and perform edge-preserving images.

METHODS

In proposed DIWGAN, two interactive generators are used to synthesize decomposed images of two basis materials by modeling the spatial and spectral correlations from input DECT reconstructed images, and the corresponding discriminators are employed to distinguish the difference between the generated images and labels. The DECT images reconstructed from high- and low-energy bins are sent to two generators separately, and each generator synthesizes one material-specific image, thereby ensuring the specificity of the network modeling. In addition, the information from different energy bins is exploited through the feature sharing of two generators. During decomposition model training, a hybrid loss function including L₁ loss, edge loss, and adversarial loss is incorporated to preserve the texture and edges in the generated images. Additionally, a selector is employed to define the generator that should be trained in each iteration, which can ensure the modeling ability of two different generators and improve the material decomposition accuracy. The performance of the proposed method is evaluated using digital phantom, XCAT phantom, and real data from a mouse.

RESULTS

On the digital phantom, the regions of bone and soft tissue are strictly and accurately separated using the trained decomposition model. The material densities in different bone and soft-tissue regions are near the ground truth, and the error of material densities is lower than 3 mg/ml. The results from XCAT phantom show that the material-specific images generated by directed matrix inversion and iterative decomposition methods have severe noise and artifacts. Regarding to the learning-based methods, the decomposed images of fully convolutional network (FCN) and butterfly network (Butterfly-Net) still contain varying degrees of artifacts, while proposed DIWGAN can yield high quality images. Compared to Butterfly-Net, the root-mean-square error (RMSE) of soft-tissue images generated by the DIWGAN decreased by 0.01 g/ml, whereas the peak-signal-to-noise ratio (PSNR) and structural similarity (SSIM) of the soft-tissue images reached 31.43 dB and 0.9987, respectively. The mass densities of the decomposed materials are nearest to the ground truth when using the DIWGAN method. The noise standard deviation of the decomposition images reduced by 69%, 60%, 33%, and 21% compared with direct matrix inversion, iterative decomposition, FCN, and Butterfly-Net, respectively. Furthermore, the performance of the mouse data indicates the potential of the proposed material decomposition method in real scanned data.

CONCLUSIONS

A DECT material decomposition method based on deep learning is proposed, and the relationship between reconstructed and material-specific images is mapped by training the DIWGAN model. Results from both the simulation phantoms and real data demonstrate the advantages of this method in suppressing noise and beam-hardening artifacts.

An oblique projection modification technique (OPMT) for fast multispectral CT reconstruction

In x-ray multispectral (or photon-counting) computed tomography (MCT), the object of interest is scanned under multiple x-ray spectra, and it can acquire more information about the scanned object than conventional CT, in which only one x-ray spectrum is used. The obtained polychromatic projections are utilized to perform material-selective and energy-selective image reconstruction. Compared with the conventional single spectral CT, MCT has a superior material distinguishability. Therefore, it has wide potential applications in both medical and industrial areas. However, the nonlinearity and ill condition of the MCT problem make it difficult to get high-quality and fast convergence of images for existing MCT reconstruction algorithms. In this paper, we proposed an iterative reconstruction algorithm based on an oblique projection modification technique (OPMT) for fast basis material decomposition of MCT. In the case of geometric inconsistency, along the current x-ray path, the oblique projection modification direction not only relates to the polychromatic projection equation of the known spectrum, but it also comprehensively refers to the polychromatic projection equation information of the unknown spectra. Moreover, the ray-by-ray correction makes it applicable to geometrically consistent projection data. One feature of the proposed algorithm is its fast convergence speed. The OPMT considers the information from multiple polychromatic projection equations, which greatly speeds up the convergence of MCT reconstructed images. Another feature of the proposed algorithm is its high flexibility. The ray-by-ray correction will be suitable for any common MCT scanning mode. The proposed algorithm is validated with numerical experiments from both simulated and real data. Compared with the ASD-NC-POCS and E-ART algorithms, the proposed algorithm achieved high-quality reconstructed images while accelerating the convergence speed of them.

An effective sinogram inpainting for complementary limited-angle dual-energy computed tomography imaging using generative adversarial networks

Dual-energy computed tomography (DECT) provides more anatomical and functional information for image diagnosis. Presently, the popular DECT imaging systems need to scan at least full angle (i.e., 360°). In this study, we propose a DECT using complementary limited-angle scan (DECT-CL) technology to reduce the radiation dose and compress the spatial distribution of the imaging system. The dual-energy total scan is 180°, where the low- and high-energy scan range is the first 90° and last 90°, respectively. We describe this dual limited-angle problem as a complementary limited-angle problem, which is challenging to obtain high-quality images using traditional reconstruction algorithms. Furthermore, a complementary-sinogram-inpainting generative adversarial networks (CSI-GAN) with a sinogram loss is proposed to inpainting sinogram to suppress the singularity of truncated sinogram. The sinogram loss focuses on the data distribution of the generated sinogram while approaching the target sinogram. We use the simultaneous algebraic reconstruction technique namely, a total variable (SART-TV) algorithms for image reconstruction. Then, taking reconstructed CT images of pleural and cranial cavity slices as examples, we evaluate the performance of our method and numerically compare different methods based on root mean square error (RMSE), peak signal-to-noise ratio (PSNR) and structural similarity (SSIM). Compared with traditional algorithms, the proposed network shows advantages in numerical terms. Compared with Patch-GAN, the proposed network can also reduce the RMSE of the reconstruction results by an average of 40% and increase the PSNR by an average of 26%. In conclusion, both qualitative and quantitative comparison and analysis demonstrate that our proposed method achieves a good artifact suppression effect and can suitably solve the complementary limited-angle problem.

Multi-energy CT reconstruction using tensor nonlocal similarity and spatial sparsity regularization

Background

Multi-energy computed tomography (MECT) based on a photon-counting detector is an emerging imaging modality that collects projections at several energy bins with a single scan. However, the limited number of photons collected into the divided, narrow energy bins results in high quantum noise levels in reconstructed images. This study aims to improve MECT image quality by minimizing noise levels while retaining image details.

Methods

A novel MECT reconstruction method was proposed by exploiting the nonlocal tensor similarity among interchannel images and spatial sparsity in single-channel images. Similar patches were initially extracted from the interchannel images in spectral and spatial domains, then stacked into a new three-order tensor. Intrinsic tensor sparsity regularization that combined the Tuker and canonical polyadic (CP) low-rank decomposition techniques were applied to exploit the nonlocal similarity of the formulated tensor. Spatial sparsity in single-channel images was modeled by total variation (TV) regularization that utilizes the compressibility of gradient image. A new MECT reconstruction model was established by simultaneously incorporating the intrinsic tensor sparsity and TV regularizations. The iterative alternating minimization method was utilized to solve the reconstruction model based on a flexible framework.

Results

The proposed method was applied to the digital phantom and real mouse data to assess its feasibility and reliability. The reconstruction and decomposition results in the mouse data were encouraging and demonstrated the ability of the proposed method in noise suppression while preserving image details, not observed with other methods. Imaging data from the digital phantom illustrated this method as achieving the best intuitive reconstruction and decomposition results among all compared methods. They reduced the root mean square error (RMSE) by 89.75%, 50.75%, and 36.54% on the reconstructed images compared with analytic, TV-based, and tensor-based methods, respectively. This phenomenon was also observed with decomposition results, where the RMSE was also reduced by 97.96%, 67.74%, 72.05%, respectively.

Conclusions

In this study, we proposed a reconstruction method for photon counting detector-based MECT, using the intrinsic tensor sparsity and TV regularizations. Improvements in noise suppression and detail preservation in the digital phantom and real mouse data were validated by the qualitative and quantitative evaluations on the reconstruction and decomposition results, verifying the potential of the proposed method in MECT reconstruction.

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.