Directional-TV algorithm for image reconstruction from limited-angular-range data

One-step inverse generation network for sparse-view dual-energy CT reconstruction and material imaging

Objective.Sparse-view dual-energy spectral computed tomography (DECT) imaging is a challenging inverse problem. Due to the incompleteness of the collected data, the presence of streak artifacts can result in the degradation of reconstructed spectral images. The subsequent material decomposition task in DECT can further lead to the amplification of artifacts and noise.Approach.To address this problem, we propose a novel one-step inverse generation network (OIGN) for sparse-view dual-energy CT imaging, which can achieve simultaneous imaging of spectral images and materials. The entire OIGN consists of five sub-networks that form four modules, including the pre-reconstruction module, the pre-decomposition module, and the following residual filtering module and residual decomposition module. The residual feedback mechanism is introduced to synchronize the optimization of spectral CT images and materials.Main results.Numerical simulation experiments show that the OIGN has better performance on both reconstruction and material decomposition than other state-of-the-art spectral CT imaging algorithms. OIGN also demonstrates high imaging efficiency by completing two high-quality imaging tasks in just 50 seconds. Additionally, anti-noise testing is conducted to evaluate the robustness of OIGN.Significance.These findings have great potential in high-quality multi-task spectral CT imaging in clinical diagnosis.

Accurate Reconstruction of Multiple Basis Images Directly From Dual Energy CT Data

OBJECTIVE

We develop optimization-based algorithms to accurately reconstruct multiple ( 2) basis images directly from dual-energy (DE) data in CT.

METHODS

In medical and industrial CT imaging, some basis materials such as bone, metals, and contrast agents of interest are confined often spatially within regions in the image. Exploiting this observation, we develop an optimization-based algorithm to reconstruct, directly from DE data, basis-region images from which multiple ( 2) basis images and virtual monochromatic images (VMIs) can be obtained over the entire image array.

RESULTS

We conduct experimental studies using simulated and real DE data in CT, and evaluate basis images and VMIs obtained in terms of visual inspection and quantitative metrics. The study results reveal that the algorithm developed can accurately and robustly reconstruct multiple ( 2) basis images directly from DE data.

CONCLUSIONS

The developed algorithm can yield accurate multiple ( 2) basis images, VMIs, and physical quantities of interest from DE data in CT.

SIGNIFICANCE

The work may provide insights into the development of practical procedures for reconstructing multiple basis images, VMIs, and physical quantities from DE data in applications. The work can be extended to reconstruct multiple basis images in multi-spectral or photon-counting CT.

A directional total variation minimization algorithm for isotropic resolution in digital breast tomosynthesis

An optimization-based image reconstruction algorithm is developed for contrast enhanced digital breast tomosynthesis (DBT) using dual-energy scanning. The algorithm minimizes directional total variation (TV) with a data discrepancy and non-negativity constraints. Iodinated contrast agent (ICA) imaging is performed by reconstructing images from dual-energy DBT data followed by weighted subtraction. Physical DBT data is acquired with a Siemens Mammomat scanner of a structured breast phantom with ICA inserts. Results are shown for both directional TV minimization and filtered back-projection for reference. It is seen that directional TV is able to substantially reduce depth blur for the ICA objects.

Directional TV algorithm for image reconstruction from sparse-view projections in EPR imaging

Objective.Electron paramagnetic resonance (EPR) imaging is an advanced in vivo oxygen imaging modality. The main drawback of EPR imaging is the long scanning time. Sparse-view projections collection is an effective fast scanning pattern. However, the commonly-used filtered back projection (FBP) algorithm is not competent to accurately reconstruct images from sparse-view projections because of the severe streak artifacts. The aim of this work is to develop an advanced algorithm for sparse reconstruction of 3D EPR imaging.Methods.The optimization based algorithms including the total variation (TV) algorithm have proven to be effective in sparse reconstruction in EPR imaging. To further improve the reconstruction accuracy, we propose the directional TV (DTV) model and derive its Chambolle-Pock solving algorithm.Results.After the algorithm correctness validation on simulation data, we explore the sparse reconstruction capability of the DTV algorithm via a simulated six-sphere phantom and two real bottle phantoms filled with OX063 trityl solution and scanned by an EPR imager with a magnetic field strength of 250 G.Conclusion.Both the simulated and real data experiments show that the DTV algorithm is superior to the existing FBP and TV-type algorithms and a deep learning based method according to visual inspection and quantitative evaluations in sparse reconstruction of EPR imaging.Significance.These insights gained in this work may be used in the development of fast EPR imaging workflow of practical significance.

An interactive method based on multi-objective optimization for limited-angle CT reconstruction

Objective. Limited-angle x-ray computed tomography (CT) is a typical ill-posed inverse problem, leading to artifacts in the reconstructed image due to the incomplete projection data. Most iteration CT reconstruction methods involve optimization for a single object. This paper explores a multi-objective optimization model and an interactive method based on multi-objective optimization to suppress the artifacts of limited-angle CT.Approach. The model includes two objective functions on the dual domain within the data consistency constraint. In the interactive method, the structural similarity index measure (SSIM) is regarded as the value function of the decision maker (DM) firstly. Secondly, the DM arranges the objective functions of the multi-objective optimization model to be optimized according to their absolute importance. Finally, the SSIM and the simulated annealing (SA) method help the DM choose the desirable reconstruction image by improving the SSIM value during the iteration process.Main results. Simulation and real data experiments demonstrate that the artifacts can be suppressed by the proposed method, and the results were superior to those reconstructed by the other three reconstruction methods in preserving the edge structure of the image.Significance. The proposed interactive method based on multi-objective optimization shows some potential advantages over classical single object optimization methods.

Directional TV algorithm for fast EPR imaging

Precise radiation guided by oxygen images has demonstrated superiority over the traditional radiation methods. Electron paramagnetic resonance (EPR) imaging has proven to be the most advanced oxygen imaging modality. However, the main drawback of EPR imaging is the long scan time. For each projection, we usually need to collect the projection many times and then average them to achieve high signal-to-noise ratio (SNR). One approach to fast scan is to reduce the repeating time for each projection. While the projections would be noisy and thus the traditional commonly-use filtered backprojection (FBP) algorithm would not be capable of accurately reconstructing images. Optimization-based iterative algorithms may accurately reconstruct images from noisy projections for they may incorporate prior information into optimization models. Based on the total variation (TV) algorithms for EPR imaging, in this work, we propose a directional TV (DTV) algorithm to further improve the reconstruction accuracy. We construct the DTV constrained, data divergence minimization (DTVcDM) model, derive its Chambolle-Pock (CP) solving algorithm, validate the correctness of the whole algorithm, and perform evaluations via simulated and real data. The experimental results show that the DTV algorithm outperforms the existing TV and FBP algorithms in fast EPR imaging. Compared to the standard FBP algorithm, the proposed algorithm may achieve 10 times of acceleration.

A nonlinear total variation based computed tomography (CT) image reconstruction method using gradient reinforcement

Compressed sensing-based reconstruction algorithms have been proven to be more successful than analytical or iterative methods for sparse computed tomography (CT) imaging by narrowing down the solution set thanks to its ability to seek a sparser solution. Total variation (TV), one of the most popular sparsifiers, exploits spatial continuity of features by restricting variation between two neighboring pixels in each direction as using partial derivatives. When the number of projections is much fewer than the one in conventional CT, which results in much less sampling rate than the minimum required one, TV may not provide satisfactory results. In this study, a new regularizer is proposed which seeks for a sparser solution by reinforcing the gradient of TV and empowering the spatial continuity of features. The experiments are done by using both analitical phantom and real human CT images and the results are compared with conventional, four-directional, and directional TV algorithms by using contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR) and Structural Similarity Index (SSIM) metrics. Both quantitative and visual evaluations show that the proposed method is promising for sparse CT image reconstruction by reducing the background noise while preserving the features and edges.

Prototyping optimization-based image reconstructions from limited-angular-range data in dual-energy CT

Image reconstruction from data collected over full-angular range (FAR) in dual-energy CT (DECT) is well-studied. There exists interest in DECT with advanced scan configurations in which data are collected only over limited-angular ranges (LARs) for meeting unique workflow needs in certain practical imaging applications, and thus in the algorithm development for image reconstruction from such LAR data. The objective of the work is to investigate and prototype image reconstructions in DECT with LAR scans. We investigate and prototype optimization programs with various designs of constraints on the directional-total-variations (DTVs) of virtual monochromatic images and/or basis images, and derive the DTV algorithms to numerically solve the optimization programs for achieving accurate image reconstruction from data collected in a slew of different LAR scans. Using simulated and real data acquired with low- and high-kV spectra over LARs, we conduct quantitative studies to demonstrate and evaluate the optimization programs and their DTV algorithms developed. As the results of the numerical studies reveal, while the DTV algorithms yield images of visual quality and quantitative accuracy comparable to that of the existing algorithms from FAR data, the former reconstruct images with improved visualization, reduced artifacts, and also enhanced quantitative accuracy when applied to LAR data in DECT. Optimization-based, one-step algorithms, including the DTV algorithms demonstrated, can be developed for quantitative image reconstruction from spectral data collected over LARs of extents that are considerably smaller than the FAR in DECT. The theoretical and numerical results obtained can be exploited for prototyping designs of optimization-based reconstructions and LAR scans in DECT, and they may also yield insights into the development of reconstruction procedures in practical DECT applications. The approach and algorithms developed can naturally be applied to investigating image reconstruction from LAR data in multi-spectral and photon-counting CT.

Multiobjective optimization guided by image quality index for limited-angle CT image reconstruction

BACKGROUND

Due to the incomplete projection data collected by limited-angle computed tomography (CT), severe artifacts are present in the reconstructed image. Classical regularization methods such as total variation (TV) minimization, ℓ0 minimization, are unable to suppress artifacts at the edges perfectly. Most existing regularization methods are single-objective optimization approaches, stemming from scalarization methods for multiobjective optimization problems (MOP).

OBJECTIVE

To further suppress the artifacts and effectively preserve the edge structures of the reconstructed image.

METHOD

This study presents a multiobjective optimization model incorporates both data fidelity term and ℓ0-norm of the image gradient as objective functions. It employs an iterative approach different from traditional scalarization methods, using the maximization of structural similarity (SSIM) values to guide optimization rather than minimizing the objective function.The iterative method involves two steps, firstly, simultaneous algebraic reconstruction technique (SART) optimizes the data fidelity term using SSIM and the Simulated Annealing (SA) algorithm for guidance. The degradation solution is accepted in the form of probability, and guided image filtering (GIF) is introduced to further preserve the image edge when the degradation solution is rejected. Secondly, the result from the first step is integrated into the second objective function as a constraint, we use ℓ0 minimization to optimize ℓ0-norm of the image gradient, and the SSIM, SA algorithm and GIF are introduced to guide optimization process by improving SSIM value like the first step.

RESULTS

With visual inspection, the peak signal-to-noise ratio (PSNR), root mean square error (RMSE), and SSIM values indicate that our approach outperforms other traditional methods.

CONCLUSIONS

The experiments demonstrate the effectiveness of our method and its superiority over other classical methods in artifact suppression and edge detail restoration.

3D directional gradient L0 norm minimization guided limited-view reconstruction in a dual-panel positron emission mammography

BACKGROUND

Dual-panel PET is often used for local organ imaging, especially breast imaging, due to its simple structure, high sensitivity, good in-plane resolution, and straightforward fusion with other imaging modalities. Nevertheless, because of data loss caused by the dual-panel structure, using conventional image reconstruction methods results in limited-view artifacts and low image quality in dual-panel positron emission mammography (PEM), which may seriously affect the diagnosis. To mitigate the limited-view artifacts in the dual-panel PEM, we propose a 3D directional gradient L₀ norm minimization (3D-DL₀) guided reconstruction method.

METHODS

The detailed derivation and reasonable simplification of the 3D-DL₀ algorithm are given first. Using this algorithm, we then obtain a prior image with edge recovery but contrast loss. To limit the solution space, the 3D-DL₀ prior is introduced into the Maximum a Posteriori reconstruction. Meanwhile, a space-invariant point spread function is also implemented to restore image contrast and boundaries. Finally, the reconstructed images with limited-view artifact suppression are obtained. The proposed method was evaluated using the data acquired from physical phantoms and patients with breast tumors on a commercial dual-panel PET system.

RESULTS

The qualitative and quantitative studies for phantom data and the blind reader study for clinical data show that the proposed method is more effective in reaching a balance between artifact elimination and image contrast improvement compared with various limited-view reconstruction methods. In addition, the iteration process of the method is proved convergent numerically.

CONCLUSIONS

The image quality improvement confirms the potential value of the proposed reconstruction algorithm to address the limited-view problem, and thus improve diagnostic accuracy in dual-panel PEM imaging.

4D-image reconstruction directly from limited-angular-range data in continuous-wave electron paramagnetic resonance imaging

OBJECTIVE

We investigate and develop optimization-based algorithms for accurate reconstruction of four-dimensional (4D)-spectral-spatial (SS) images directly from data collected over limited angular ranges (LARs) in continuous-wave (CW) electron paramagnetic resonance imaging (EPRI).

METHODS

Basing on a discrete-to-discrete data model devised in CW EPRI employing the Zeeman-modulation (ZM) scheme for data acquisition, we first formulate the image reconstruction problem as a convex, constrained optimization program that includes a data fidelity term and also constraints on the individual directional total variations (DTVs) of the 4D-SS image. Subsequently, we develop a primal-dual-based DTV algorithm, simply referred to as the DTV algorithm, to solve the constrained optimization program for achieving image reconstruction from data collected in LAR scans in CW-ZM EPRI.

RESULTS

We evaluate the DTV algorithm in simulated- and real-data studies for a variety of LAR scans of interest in CW-ZM EPRI, and visual and quantitative results of the studies reveal that 4D-SS images can be reconstructed directly from LAR data, which are visually and quantitatively comparable to those obtained from data acquired in the standard, full-angular-range (FAR) scan in CW-ZM EPRI.

CONCLUSION

An optimization-based DTV algorithm is developed for accurately reconstructing 4D-SS images directly from LAR data in CW-ZM EPRI. Future work includes the development and application of the optimization-based DTV algorithm for reconstructions of 4D-SS images from FAR and LAR data acquired in CW EPRI employing schemes other than the ZM scheme.

SIGNIFICANCE

The DTV algorithm developed may be exploited potentially for enabling and optimizing CW EPRI with minimized imaging time and artifacts by acquiring data in LAR scans.

A Simple but Universal Fully Linearized ADMM Algorithm for Optimization Based Image Reconstruction

Background and Objective

Optimization based image reconstruction algorithm is an advanced algorithm in medical imaging. However, the corresponding solving algorithm is challenging because the optimization model is usually large-scale and non-smooth. This work aims to devise a simple but universal solver for optimization models.

Methods

The alternating direction method of multipliers (ADMM) algorithm is a simple and effective solver of the optimization models. However, there always exists a sub-problem that has not closed-form solution. One may use gradient descent algorithm to solve this sub-problem, but the step-size selection via line search is time-consuming. Or, one may use fast Fourier transform (FFT) to get a closed-form solution if the system matrix and the sparse transform matrix are both of special structure. In this work, we propose a simple but universal fully linearized ADMM (FL-ADMM) algorithm that avoids line search to determine step-size and applies to system matrix and sparse transform of any structures.

Results

We derive the FL-ADMM algorithm instances for three total variation (TV) models in 2D computed tomography (CT). Further, we validate and evaluate one FL-ADMM algorithm and explore how the two important factors impact convergence rate. Also, we compare this algorithm with the Chambolle-Pock algorithm via real CT phantom reconstructions. These studies show that the FL-ADMM algorithm may accurately solve optimization models in image reconstruction.

Conclusion

The FL-ADMM algorithm is a simple, effective, convergent and universal solver of optimization models in image reconstruction. Compared to the existing ADMM algorithms, the new algorithm does not need time-consuming step-size line-search or special demand to system matrix and sparse transform. It is a rapid prototyping tool for optimization based image reconstruction.

Structure-guided computed tomography reconstruction from limited-angle projections

Limited-angle computed tomography (CT) imaging is one of the common imaging problems. The reconstructed images often encounter obvious artifacts and structure degradation. In recent years, the recoverability prior of image structure has been widely explored in limited-angle CT reconstruction, and the image quality has been greatly improved. However, the artifacts and structure degradation still exist. In this study, we establish a new reconstruction model based on weighted relative structure (wRS) determined by image gradients, which serves as weights to guide image reconstruction in order to reduce artifacts and preserve structures. Then, we develop an efficient algorithm using a surrogate function to solve this model. Moreover, this method is compared with some of other popular reconstruction methods, such as anisotropic total variation method and image gradient L0 norm minimization method and so on. Experiments on digital phantoms, real carved cheese and walnut projection are reported to demonstrate its superiority. Several quantitative indices including RMSE, PSNR, and SSIM of the reconstruction images from 90°-data of FORBILD head phantom are 0.0120, 43.52, and 0.9961. The experimental results indicate that the image obtained by our method is the closest to reference image. By comparing reconstruction images or their residual images, images reconstructed from real CT data, the experimental results of the residual images and the respective quantitative data analysis also demonstrate that the images reconstructed using our new method suffer from less artifacts and structure degradation.

Accurate Image Reconstruction in Dual-Energy CT with Limited-Angular-Range Data Using a Two-Step Method

Dual-energy CT (DECT) with scans over limited-angular ranges (LARs) may allow reductions in scan time and radiation dose and avoidance of possible collision between the moving parts of a scanner and the imaged object. The beam-hardening (BH) and LAR effects are two sources of image artifacts in DECT with LAR data. In this work, we investigate a two-step method to correct for both BH and LAR artifacts in order to yield accurate image reconstruction in DECT with LAR data. From low- and high-kVp LAR data in DECT, we first use a data-domain decomposition (DDD) algorithm to obtain LAR basis data with the non-linear BH effect corrected for. We then develop and tailor a directional-total-variation (DTV) algorithm to reconstruct from the LAR basis data obtained basis images with the LAR effect compensated for. Finally, using the basis images reconstructed, we create virtual monochromatic images (VMIs), and estimate physical quantities such as iodine concentrations and effective atomic numbers within the object imaged. We conduct numerical studies using two digital phantoms of different complexity levels and types of structures. LAR data of low- and high-kVp are generated from the phantoms over both single-arc (SA) and two-orthogonal-arc (TOA) LARs ranging from 14∘ to 180∘. Visual inspection and quantitative assessment of VMIs obtained reveal that the two-step method proposed can yield VMIs in which both BH and LAR artifacts are reduced, and estimation accuracy of physical quantities is improved. In addition, concerning SA and TOA scans with the same total LAR, the latter is shown to yield more accurate images and physical quantity estimations than the former. We investigate a two-step method that combines the DDD and DTV algorithms to correct for both BH and LAR artifacts in image reconstruction, yielding accurate VMIs and estimations of physical quantities, from low- and high-kVp LAR data in DECT. The results and knowledge acquired in the work on accurate image reconstruction in LAR DECT may give rise to further understanding and insights into the practical design of LAR scan configurations and reconstruction procedures for DECT applications.

Sam's Net: A Self-Augmented Multistage Deep-Learning Network for End-to-End Reconstruction of Limited Angle CT

Limited angle reconstruction is a typical ill-posed problem in computed tomography (CT). Given incomplete projection data, images reconstructed by conventional analytical algorithms and iterative methods suffer from severe structural distortions and artifacts. In this paper, we proposed a self-augmented multi-stage deep-learning network (Sam's Net) for end-to-end reconstruction of limited angle CT. With the merit of the alternating minimization technique, Sam's Net integrates multi-stage self-constraints into cross-domain optimization to provide additional constraints on the manifold of neural networks. In practice, a sinogram completion network (SCNet) and artifact suppression network (ASNet), together with domain transformation layers constitute the backbone for cross-domain optimization. An online self-augmentation module was designed following the manner defined by alternating minimization, which enables a self-augmented learning procedure and multi-stage inference manner. Besides, a substitution operation was applied as a hard constraint for the solution space based on the data fidelity and a learnable weighting layer was constructed for data consistency refinement. Sam's Net forms a new framework for ill-posed reconstruction problems. In the training phase, the self-augmented procedure guides the optimization into a tightened solution space with enriched diverse data distribution and enhanced data consistency. In the inference phase, multi-stage prediction can improve performance progressively. Extensive experiments with both simulated and practical projections under 90-degree and 120-degree fan-beam configurations validate that Sam's Net can significantly improve the reconstruction quality with high stability and robustness.

Image reconstruction from data over two orthogonal arcs of limited-angular ranges

PURPOSE

Computed tomography (CT) scanning over limited-angular ranges (LARs) is of practical interest in possible reduction of imaging dose and time and in design of non-standard scans. This work aims to investigate image reconstruction for two non-overlapping arcs of LARs, and to demonstrate that they may allow more accurate image reconstruction than may a single arc of LAR.

METHODS

We consider a configuration with two non-overlapping arcs of LARs α₁ and α₂ , whose centers are separated by 90°, and refer to it as a two-orthogonal-arc configuration. Data are generated from a chest phantom with two-orthogonal-arc configurations over total angular coverage α_τ = α₁ + α₂ ranging from 18° to 180°, and images are reconstructed subsequently by use of the directional-total-variation (DTV) algorithm. For comparison, we also consider image reconstruction for a single-arc configuration of angular range α_τ . Quantitative metrics such as the normalized root-mean-square-error (nRMSE) are used for evaluation of image reconstruction accuracy.

RESULTS

Visual inspection and quantitative analysis of images reconstructed reveal that a two-orthogonal-arc configuration generally yields more accurate image reconstruction than does its single-arc counterpart. As total angular range α_τ increases, the DTV algorithm yields image reconstruction with enhanced accuracy, as expected. Also, if α_τ remains constant, the two-orthogonal-arc configuration with α₁ = α₂ generally leads to image reconstruction more accurate than those of two-orthogonal-arc configurations with α₁ ≠ α₂ , as the nRMSE of the former can be lower than that of the latter for up to more than one order of magnitude.

CONCLUSIONS

Appropriately designed two-orthogonal-arc configurations may be exploited for improving image-reconstruction accuracy in CT imaging with reduced angular coverage. This study may yield insights into the design of innovative CT scans for lowering scan time and radiation dose, and/or for avoiding scan collision in, e.g., C-arm CT.

Exterior computed tomography image reconstruction based on anisotropic relative total variation in polar coordinates

In computed tomography (CT) image reconstruction problems, exterior CT is an important application in industrial non-destructive testing (NDT). Different from the limited-angle problem that misses part of the rotation angle, the rotation angle of the exterior problem is complete, but for each rotation angle, the projection data through the central region of the object cannot be collected, so that the exterior CT problem is ill-posed inverse problem. The results of traditional reconstruction methods like filtered back-projection (FBP) and simultaneous algebra reconstruction technique (SART) have artifacts along the radial direction edges for exterior CT reconstruction. In this study, we propose and test an anisotropic relative total variation in polar coordinates (P-ARTV) model for addressing the exterior CT problem. Since relative total variation (RTV) can effectively distinguish edges from noises, and P-ARTV with different weights in radial and tangential directions can effectively enhance radial edges, a two-step iteration algorithm was developed to solve the P-ARTV model in this study. The fidelity term and the regularization term are solved in Cartesian and polar coordinate systems, respectively. Numerical experiments show that our new model yields better performance than the existing state-of-the-art algorithms.

Ultra-limited-angle CT image reconstruction algorithm based on reweighting and edge-preserving

BACKGROUND

Ultra-limited-angle image reconstruction problem with a limited-angle scanning range less than or equal to π2 is severely ill-posed. Due to the considerably large condition number of a linear system for image reconstruction, it is extremely challenging to generate a valid reconstructed image by traditional iterative reconstruction algorithms.

OBJECTIVE

To develop and test a valid ultra-limited-angle CT image reconstruction algorithm.

METHODS

We propose a new optimized reconstruction model and Reweighted Alternating Edge-preserving Diffusion and Smoothing algorithm in which a reweighted method of improving the condition number is incorporated into the idea of AEDS image reconstruction algorithm. The AEDS algorithm utilizes the property of image sparsity to improve partially the results. In experiments, the different algorithms (the Pre-Landweber, AEDS algorithms and our algorithm) are used to reconstruct the Shepp-Logan phantom from the simulated projection data with noises and the flat object with a large ratio between length and width from the real projection data. PSNR and SSIM are used as the quantitative indices to evaluate quality of reconstructed images.

RESULTS

Experiment results showed that for simulated projection data, our algorithm improves PSNR and SSIM from 22.46db to 39.38db and from 0.71 to 0.96, respectively. For real projection data, our algorithm yields the highest PSNR and SSIM of 30.89db and 0.88, which obtains a valid reconstructed result.

CONCLUSIONS

Our algorithm successfully combines the merits of several image processing and reconstruction algorithms. Thus, our new algorithm outperforms significantly other two algorithms and is valid for ultra-limited-angle CT image reconstruction.

Dual-energy CT imaging with limited-angular-range data

In dual-energy computed tomography (DECT), low- and high-kVp data are collected often over a full-angular range (FAR) of 360^○. While there exists strong interest in DECT with low- and high-kVp data acquired over limited-angular ranges (LARs), there remains little investigation of image reconstruction in DECT with LAR data.Objective: We investigate image reconstruction with minimized LAR artifacts from low- and high-kVp data over LARs of ≤180^○by using a directional-total-variation (DTV) algorithm.Methods: Image reconstruction from LAR data is formulated as a convex optimization problem in which data-l₂is minimized with constraints on image's DTVs along orthogonal axes. We then achieve image reconstruction by applying the DTV algorithm to solve the optimization problem. We conduct numerical studies from data generated over arcs of LARs, ranging from 14^○to 180^○, and perform visual inspection and quantitative analysis of images reconstructed.Results: Monochromatic images of interest obtained with the DTV algorithm from LAR data show substantially reduced artifacts that are observed often in images obtained with existing algorithms. The improved image quality also leads to accurate estimation of physical quantities of interest, such as effective atomic number and iodine-contrast concentration.Conclusion: Our study reveals that from LAR data of low- and high-kVp, monochromatic images can be obtained that are visually, and physical quantities can be estimated that are quantitatively, comparable to those obtained in FAR DECT.Significance: As LAR DECT is of high practical application interest, the results acquired in the work may engender insights into the design of DECT with LAR scanning configurations of practical application significance.

Deep learning reconstruction of digital breast tomosynthesis images for accurate breast density and patient-specific radiation dose estimation

The two-dimensional nature of mammography makes estimation of the overall breast density challenging, and estimation of the true patient-specific radiation dose impossible. Digital breast tomosynthesis (DBT), a pseudo-3D technique, is now commonly used in breast cancer screening and diagnostics. Still, the severely limited 3rd dimension information in DBT has not been used, until now, to estimate the true breast density or the patient-specific dose. This study proposes a reconstruction algorithm for DBT based on deep learning specifically optimized for these tasks. The algorithm, which we name DBToR, is based on unrolling a proximal-dual optimization method. The proximal operators are replaced with convolutional neural networks and prior knowledge is included in the model. This extends previous work on a deep learning-based reconstruction model by providing both the primal and the dual blocks with breast thickness information, which is available in DBT. Training and testing of the model were performed using virtual patient phantoms from two different sources. Reconstruction performance, and accuracy in estimation of breast density and radiation dose, were estimated, showing high accuracy (density <±3%; dose <±20%) without bias, significantly improving on the current state-of-the-art. This work also lays the groundwork for developing a deep learning-based reconstruction algorithm for the task of image interpretation by radiologists.

Dual-energy CT imaging over non-overlapping, orthogonal arcs of limited-angular ranges

BACKGROUND

Interest exists in dual-energy computed tomography (DECT) imaging with scanning arcs of limited-angular ranges (LARs) for reducing scan time and radiation dose, and for enabling scan configurations of C-arm CT that can avoid possible collision between the rotating X-ray tube/detector and the imaged subject.

OBJECTIVE

In this work, we investigate image reconstruction for a type of configurations of practical DECT interest, referred to as the two-orthogonal-arc configuration, in which low- and high-kVp data are collected over two non-overlapping arcs of equal LAR α, ranging from 30° to 90°, separated by 90°. The configuration can readily be implemented, e.g., on CT with dual sources separated by 90° or with the slow-kVp-switching technique.

METHODS

The directional-total-variation (DTV) algorithm developed previously for image reconstruction in conventional, single-energy CT is tailored to enable image reconstruction in DECT with two-orthogonal-arc configurations.

RESULTS

Performing visual inspection and quantitative analysis of monochromatic images obtained and effective atomic numbers estimated, we observe that the monochromatic images of the DTV algorithm from LAR data are with substantially reduced LAR artifacts, which are observed otherwise in those of existing algorithms, and thus visually correlate reasonably well, in terms of metrics PCC and nMI, with their reference images obtained from full-angular-range data. In addition, effective atomic numbers estimated from LAR data of DECT with two-orthogonal-arc configurations are in reasonable agreement, with relative errors up to ∼  10%, with those estimated from full-angular-range data in DECT.

CONCLUSIONS

The results acquired in the work may yield insights into the design of LAR configurations of practical dual-energy application relevance in diagnostic CT or C-arm CT imaging.