Simultaneous algebraic reconstruction technique based on guided image filtering

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.

Image reconstruction method for incomplete CT projection based on self-guided image filtering

In some fields of medical diagnosis or industrial nondestructive testing, it is difficult to obtain complete computed tomography (CT) data due to the limitation of radiation dose or other factors. Therefore, image reconstruction of incomplete projection data is the focus of this paper. In this paper, a new image reconstruction model based on self-guided image filtering (SGIF) term is proposed for few-view and segmental limited-angle (SLA) CT reconstruction. Then the alternating direction method (ADM) is used to solve this model. For simplicity, we call it ADM-SGIF method. The key idea of ADM-SGIF method is to use the reconstructed image itself as a reference and utilize its structural features to guide CT reconstruction. This method can effectively preserve image structures and remove shading artifacts. To validate the effectiveness of the proposed reconstruction method, we conduct digital phantom and real CT data experiments. The results indicate that ADM-SGIF method outperforms competing methods, including total variation (TV), relative total variation (RTV), and L0-norm minimization solved by ADM (ADM-L0) methods, in both subjective and objective evaluations.

Sparse-view CT reconstruction based on group-based sparse representation using weighted guided image filtering

OBJECTIVES

In the past, guided image filtering (GIF)-based methods often utilized total variation (TV)-based methods to reconstruct guidance images. And they failed to reconstruct the intricate details of complex clinical images accurately. To address these problems, we propose a new sparse-view CT reconstruction method based on group-based sparse representation using weighted guided image filtering.

METHODS

In each iteration of the proposed algorithm, the result constrained by the group-based sparse representation (GSR) is used as the guidance image. Then, the weighted guided image filtering (WGIF) was used to transfer the important features from the guidance image to the reconstruction of the SART method.

RESULTS

Three representative slices were tested under 64 projection views, and the proposed method yielded the best visual effect. For the shoulder case, the PSNR can achieve 48.82, which is far superior to other methods.

CONCLUSIONS

The experimental results demonstrate that our method is more effective in preserving structures, suppressing noise, and reducing artifacts compared to other methods.

An image reconstruction method for transmission computed tomography with the constraint of the linear attenuation coefficients

For the reconstructed image of transmission computed tomography, the linear attenuation coefficients of the diagnosed object may improve the image quality by adding additional constraint besides the projection data. In the present work, an image reconstruction method with the constraint of the linear attenuation coefficients is developed and two models including a classical numerical Shepp-Logan model and a Monte Carlo model are used to show the corresponding benefits. The results indicate that the number of the projection angles is potentially decreased to 1/3 of itself while the quality of the reconstructed image is not deteriorated.

Sparse Angle CBCT Reconstruction Based on Guided Image Filtering

Cone-beam Computerized Tomography (CBCT) has the advantages of high ray utilization and detection efficiency, short scan time, high spatial and isotropic resolution. However, the X-rays emitted by CBCT examination are harmful to the human body, so reducing the radiation dose without damaging the reconstruction quality is the key to the reconstruction of CBCT. In this paper, we propose a sparse angle CBCT reconstruction algorithm based on Guided Image FilteringGIF, which combines the classic Simultaneous Algebra Reconstruction Technique(SART) and the Total p-Variation (TpV) minimization. Due to the good edge-preserving ability of SART and noise suppression ability of TpV minimization, the proposed method can suppress noise and artifacts while preserving edge and texture information in reconstructed images. Experimental results based on simulated and real-measured CBCT datasets show the advantages of the proposed method.

NeuRec: Incorporating Interpatient prior to Sparse-View Image Reconstruction for Neurorehabilitation

Medical imaging technologies such as computed tomography (CT) and magnetic resonance imaging (MRI) imaging are indispensable for contemporary neurorehabilitation diagnostics, intervention, and monitoring. It would be desirable to reconstruct images from sparse measurements to reduce the ionizing radiation and motion artifacts. Although recent coordinate-based representation methods have shown promise advances for sparse-view reconstruction, they overfit a single MLP on a single patient. In this work, we generalize it across many patients by incorporating an interpatient prior into the ill-posed inverse/reconstruction problem, which is the missing ingredient in the previous works. The experiment demonstrates that our method significantly improves image quality over the state-of-the-art both qualitatively and quantitatively. Thus, our method provides a powerful and principled means to deal with the measurement-scarce problem.

A Realistic Breast Phantom Proposal for 3D Image Reconstruction in Digital Breast Tomosynthesis

Objectives: Iterative (eg, simultaneous algebraic reconstruction technique [SART]) and analytical (eg, filtered back projection [FBP]) image reconstruction techniques have been suggested to provide adequate three-dimensional (3D) images of the breast for capturing microcalcifications in digital breast tomosynthesis (DBT). To decide on the reconstruction method in clinical DBT, it must first be tested in a simulation resembling the real clinical environment. The purpose of this study is to introduce a 3D realistic breast phantom for determining the reconstruction method in clinical applications. Methods: We designed a 3D realistic breast phantom with varying dimensions (64³-512³) mimicking some structures of a real breast such as milk ducts, lobules, and ribs using TomoPhantom software. We generated microcalcifications, which mimic cancerous cells, with a separate MATLAB code and embedded them into the phantom for testing and benchmark studies in DBT. To validate the characterization of the phantom, we tested the distinguishability of microcalcifications by performing 3D image reconstruction methods (SART and FBP) using Laboratory of Computer Vision (LAVI) open-source reconstruction toolbox. Results: The creation times of the proposed realistic breast phantom were seconds of 2.5916, 8.4626, 57.6858, and 472.1734 for 64³, 128³, 256³, and 512³, respectively. We presented reconstructed images and quantitative results of the phantom for SART (1-2-4-8 iterations) and FBP, with 11 to 23 projections. We determined qualitatively and quantitatively that SART (2-4 iter.) yields better results than FBP. For example, for 23 projections, the contrast-to-noise ratio (CNR) values of SART (2 iter.) and FBP were 2.871 and 0.497, respectively. Conclusions: We created a computationally efficient realistic breast phantom that is eligible for reconstruction and includes anatomical structures and microcalcifications, successfully. By proposing this breast phantom, we provided the opportunity to test which reconstruction methods can be used in clinical applications vary according to various parameters such as the No. of iterations and projections in DBT.

Computed Tomography as a Characterization Tool for Engineered Scaffolds with Biomedical Applications

The ever-growing field of materials with applications in the biomedical field holds great promise regarding the design and fabrication of devices with specific characteristics, especially scaffolds with personalized geometry and architecture. The continuous technological development pushes the limits of innovation in obtaining adequate scaffolds and establishing their characteristics and performance. To this end, computed tomography (CT) proved to be a reliable, nondestructive, high-performance machine, enabling visualization and structure analysis at submicronic resolutions. CT allows both qualitative and quantitative data of the 3D model, offering an overall image of its specific architectural features and reliable numerical data for rigorous analyses. The precise engineering of scaffolds consists in the fabrication of objects with well-defined morphometric parameters (e.g., shape, porosity, wall thickness) and in their performance validation through thorough control over their behavior (in situ visualization, degradation, new tissue formation, wear, etc.). This review is focused on the use of CT in biomaterial science with the aim of qualitatively and quantitatively assessing the scaffolds’ features and monitoring their behavior following in vivo or in vitro experiments. Furthermore, the paper presents the benefits and limitations regarding the employment of this technique when engineering materials with applications in the biomedical field.

Three-dimensional reconstruction of integrated implosion targets from simulated small-angle pinhole images

To indirectly evaluate the asymmetry of the radiation drive under limited measurement conditions in inertial confinement fusion research, we have proposed an integral method to approximate the three-dimensional self-radiation distribution of the compressed plasma core using only four pinhole images from a single laser entrance hole at a maximum projection angle of 10°. The simultaneous algebraic reconstruction technique (SART) that uses spatial constraints provided by the prior structural information and the central pinhole image is utilized in the simulation. The simulation results showed that the normalized mean square deviation between the original distribution and reconstruction results of the central radiation area of the simulated cavity was 0.4401, and the structural similarity of the cavity radiation distribution was 0.5566. Meanwhile, using more diagnostic holes could achieve better structural similarity and lower reconstruction error. In addition, the results indicated that our new proposed method could reconstruct the distribution of a compressed plasma core in a vacuum hohlraum with high accuracy.

Noise suppression-guided image filtering for low-SNR CT reconstruction

In practical computed tomography (CT) applications, projections with low signal-to-noise ratio (SNR) are often encountered due to the reduction of radiation dose or device limitations. In these situations, classical reconstruction algorithms, like simultaneous algebraic reconstruction technique (SART), cannot reconstruct high-quality CT images. Block-matching and 3D filtering (BM3D)-based iterative reconstruction algorithm (POCS-BM3D) has remarkable effect in dealing with CT reconstruction from noisy projections. However, BM3D may restrain noise with excessive loss of details in the case of low-SNR CT reconstruction. In order to achieve a preferable trade-off between noise suppression and edge preservation, we introduce guided image filtering (GIF) into low-SNR CT reconstruction, and propose noise suppression-guided image filtering reconstruction (NSGIFR) algorithm. In each iteration of NSGIFR, the output image of SART reserves more details and is used as input image of GIF, while the image denoised by BM3D serves as guidance image of GIF. Experimental results indicate that the proposed algorithm displays outstanding performance on preserving structures and suppressing noise for low-SNR CT reconstruction. NSGIFR can achieve more superior image quality than SART, POCS-TV and POCS-BM3D in terms of visual effect and quantitative analysis. Graphical abstract Block-matching and 3D filtering (BM3D)-based iterative reconstruction algorithm (POCS-BM3D) has remarkable effect in dealing with CT reconstruction from noisy projections. However, BM3D may restrain noise with excessive loss of details in the case of low-SNR CT reconstruction. In order to achieve a preferable trade-off between noise suppression and edge preservation, we introduce guided image filtering (GIF) into low-SNR CT reconstruction, and propose noise suppression-guided image filtering reconstruction (NSGIFR) algorithm.

Simultaneous correction of sensitivity and spatial resolution in projection-based magnetic particle imaging

PURPOSE

The purpose of this study was to develop a method to simultaneously correct the spatial resolution and inhomogeneous sensitivity of a receiving coil in projection-based magnetic particle imaging - and to investigate its efficacy through simulation and experimental studies.

METHODS

Magnetic particle imaging (MPI) images were reconstructed using the simultaneous algebraic reconstruction technique (SART), and simultaneous corrections to sensitivity and spatial resolution were performed by incorporating the sensitivity map of the receiving coil and the system function into the SART algorithm. After each SART update, the regularization method - with total variation (TV) minimization - was used to suppress noise amplification and artifact generation. For comparison, MPI images were also reconstructed using the filtered backprojection (FBP) method and the FBP-truncated singular value decomposition (TSVD) method, in which the system function was deconvolved from the projection data using TSVD. In simulation studies, the sensitivity map of a second-order, gradiometer-type receiving coil was generated using the Biot-Savart law, while the system function was obtained by calculating the MPI signals induced by magnetic nanoparticles at various distances from a field-free line (FFL), using a lock-in-amplifier model. The effects of a regularization parameter for TV minimization (α), number of iterations (N), and signal-to-noise ratio (SNR) of the MPI signals on the reconstructed MPI images of a numerical phantom were evaluated, using the image profiles and percent root mean square error (PRMSE). Experimental studies involved the calculation of the system function using a tube phantom. Projection data for an A-shaped phantom were acquired using our MPI scanner, and their MPI images were reconstructed from the projection data, as described above.

RESULTS

When both the sensitivity and spatial resolution were corrected (SART-SR), the quality of the reconstructed images was seen to have improved, compared to when the spatial resolution was not corrected - or when the FBP and FBP-TSVD methods were used. When SNR was low (20), a larger α value yielded a better image. The minimum PRMSE occurred at N ≈ 200-400 and increased with increasing N thereafter. When SNR was high (100), the image quality was generally not dependent on the α value within its studied range. The PRMSE decreased slowly with increasing N, and tended to converge to a constant value. The full width at half maximum (FWHM) of the profile was obtained from the A-shaped phantom, reconstructed using the SART-SR algorithm with α = 0.05 and N = 1000. The FWHM value of the tube (2 mm diameter) in the A-shaped phantom image was found to be 2.2 mm on average, whereas those calculated from the images obtained by the FBP and FBP-TSVD methods were 4.4 and 3.0 mm on average, respectively. Spatial resolution improved when using the FBP-TSVD method as compared to the FBP method but image distortion and artifacts were observed.

CONCLUSIONS

Although further studies are necessary to optimize the parameters used in the SART algorithm and in TV minimization, the present results suggest that the proposed method is useful for improving the image quality of projection-based MPI.

Three-dimensional visualization of microvasculature from few-projection data using a novel CT reconstruction algorithm for propagation-based X-ray phase-contrast imaging

Propagation-based X-ray phase-contrast imaging (PBI) is a powerful nondestructive imaging technique that can reveal the internal detailed structures in weakly absorbing samples. Extending PBI to CT (PBCT) enables high-resolution and high-contrast 3D visualization of microvasculature, which can be used for the understanding, diagnosis and therapy of diseases involving vasculopathy, such as cardiovascular disease, stroke and tumor. However, the long scan time for PBCT impedes its wider use in biomedical and preclinical microvascular studies. To address this issue, a novel CT reconstruction algorithm for PBCT is presented that aims at shortening the scan time for microvascular samples by reducing the number of projections while maintaining the high quality of reconstructed images. The proposed algorithm combines the filtered backprojection method into the iterative reconstruction framework, and a weighted guided image filtering approach (WGIF) is utilized to optimize the intermediate reconstructed images. Notably, the homogeneity assumption on the microvasculature sample is adopted as prior knowledge, and therefore, a prior image of microvasculature structures can be acquired by a k-means clustering approach. Then, the prior image is used as the guided image in the WGIF procedure to effectively suppress streaking artifacts and preserve microvasculature structures. To evaluate the effectiveness and capability of the proposed algorithm, simulation experiments on 3D microvasculature numerical phantom and real experiments with CT reconstruction on the microvasculature sample are performed. The results demonstrate that the proposed algorithm can, under noise-free and noisy conditions, significantly reduce the artifacts and effectively preserve the microvasculature structures on the reconstructed images and thus enables it to be used for clear and accurate 3D visualization of microvasculature from few-projection data. Therefore, for 3D visualization of microvasculature, the proposed algorithm can be considered an effective approach for reducing the scan time required by PBCT.

A joint deep learning model to recover information and reduce artifacts in missing-wedge sinograms for electron tomography and beyond

We present a joint model based on deep learning that is designed to inpaint the missing-wedge sinogram of electron tomography and reduce the residual artifacts in the reconstructed tomograms. Traditional methods, such as weighted back projection (WBP) and simultaneous algebraic reconstruction technique (SART), lack the ability to recover the unacquired project information as a result of the limited tilt range; consequently, the tomograms reconstructed using these methods are distorted and contaminated with the elongation, streaking, and ghost tail artifacts. To tackle this problem, we first design a sinogram filling model based on the use of Residual-in-Residual Dense Blocks in a Generative Adversarial Network (GAN). Then, we use a U-net structured Generative Adversarial Network to reduce the residual artifacts. We build a two-step model to perform information recovery and artifacts removal in their respective suitable domain. Compared with the traditional methods, our method offers superior Peak Signal to Noise Ratio (PSNR) and the Structural Similarity Index (SSIM) to WBP and SART; even with a missing wedge of 45°, our method offers reconstructed images that closely resemble the ground truth with nearly no artifacts. In addition, our model has the advantage of not needing inputs from human operators or setting hyperparameters such as iteration steps and relaxation coefficient used in TV-based methods, which highly relies on human experience and parameter fine turning.

Super-resolution CT Image Reconstruction Based on Dictionary Learning and Sparse Representation

In this paper, a single-computed tomography (CT) image super-resolution (SR) reconstruction scheme is proposed. This SR reconstruction scheme is based on sparse representation theory and dictionary learning of low- and high-resolution image patch pairs to improve the poor quality of low-resolution CT images obtained in clinical practice using low-dose CT technology. The proposed strategy is based on the idea that image patches can be well represented by sparse coding of elements from an overcomplete dictionary. To obtain similarity of the sparse representations, two dictionaries of low- and high-resolution image patches are jointly trained. Then, sparse representation coefficients extracted from the low-resolution input patches are used to reconstruct the high-resolution output. Sparse representation is used such that the trained dictionary pair can reduce computational costs. Combined with several appropriate iteration operations, the reconstructed high-resolution image can attain better image quality. The effectiveness of the proposed method is demonstrated using both clinical CT data and simulation image data. Image quality evaluation indexes (root mean squared error (RMSE) and peak signal-to-noise ratio (PSNR)) indicate that the proposed method can effectively improve the resolution of a single CT image.

Combine TV-L1 model with guided image filtering for wide and faint ring artifacts correction of in-line x-ray phase contrast computed tomography

In practice, mis-calibrated detector pixels give rise to wide and faint ring artifacts in the reconstruction image of the In-line phase-contrast computed tomography (IL-PC-CT). Ring artifacts correction is essential in IL-PC-CT. In this study, a novel method of wide and faint ring artifacts correction was presented based on combining TV-L1 model with guided image filtering (GIF) in the reconstruction image domain. The new correction method includes two main steps namely, the GIF step and the TV-L1 step. To validate the performance of this method, simulation data and real experimental synchrotron data are provided. The results demonstrate that TV-L1 model with GIF step can effectively correct the wide and faint ring artifacts for IL-PC-CT.

Edge-preserving reconstruction from sparse projections of limited-angle computed tomography using ℓ0-regularized gradient prior

Accurate images reconstructed from limited computed tomography (CT) data are desired when reducing the X-ray radiation exposure imposed on patients. The total variation (TV), known as the l₁-norm of the image gradient magnitudes, is popular in CT reconstruction from incomplete projection data. However, as the projection data collected are from a sparse-view of the limited scanning angular range, the results reconstructed by a TV-based method suffer from blocky artifact and gradual changed artifacts near the edges, which in turn make the reconstruction images degraded. Different from the TV, the ℓ₀-norm of an image gradient counts the number of its non-zero coefficients of the image gradient. Since the regularization based on the ℓ₀-norm of the image gradient will not penalize the large gradient magnitudes, the edge can be effectively retained. In this work, an edge-preserving image reconstruction method based on l₀-regularized gradient prior was investigated for limited-angle computed tomography from sparse projections. To solve the optimization model effectively, the variable splitting and the alternating direction method (ADM) were utilized. Experiments demonstrated that the ADM-like method used for the non-convex optimization problem has better performance than other classical iterative reconstruction algorithms in terms of edge preservation and artifact reduction.