scholarly journals Noise-Robust Image Reconstruction Based on Minimizing Extended Class of Power-Divergence Measures

Entropy ◽  
2021 ◽  
Vol 23 (8) ◽  
pp. 1005
Author(s):  
Ryosuke Kasai ◽  
Yusaku Yamaguchi ◽  
Takeshi Kojima ◽  
Omar M. Abou Al-Ola ◽  
Tetsuya Yoshinaga

The problem of tomographic image reconstruction can be reduced to an optimization problem of finding unknown pixel values subject to minimizing the difference between the measured and forward projections. Iterative image reconstruction algorithms provide significant improvements over transform methods in computed tomography. In this paper, we present an extended class of power-divergence measures (PDMs), which includes a large set of distance and relative entropy measures, and propose an iterative reconstruction algorithm based on the extended PDM (EPDM) as an objective function for the optimization strategy. For this purpose, we introduce a system of nonlinear differential equations whose Lyapunov function is equivalent to the EPDM. Then, we derive an iterative formula by multiplicative discretization of the continuous-time system. Since the parameterized EPDM family includes the Kullback–Leibler divergence, the resulting iterative algorithm is a natural extension of the maximum-likelihood expectation-maximization (MLEM) method. We conducted image reconstruction experiments using noisy projection data and found that the proposed algorithm outperformed MLEM and could reconstruct high-quality images that were robust to measured noise by properly selecting parameters.

2022 ◽  
pp. 1-13
Author(s):  
Lei Shi ◽  
Gangrong Qu ◽  
Yunsong Zhao

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.


Author(s):  
Tomasz Rymarczyk ◽  
Konrad Kania ◽  
Michał Gołąbek ◽  
Jan Sikora ◽  
Michał Maj ◽  
...  

Purpose The purpose of this study is to develop a reconstruction and measurement system for data analysis using ultrasonic transmission tomography. The problem of reconstruction from the projection is encountered in practical implementation, which consists in reconstructing an image that is an estimation of an unknown object from a finite set of projection data. Reconstructive algorithms used in transmission tomography are based on linear mathematical models, which makes it necessary to process non-linear data into estimates for a finite number of projections. The application of transformation methods requires building a mathematical model in which the projection data forming the known and unknown quantities are functions with arguments from a continuous set of real numbers, determining the function describing the unknown quantities sought in the form of inverse relation and adapting it to operate on discrete and noisy data. This was done by designing a tomographic device and proprietary algorithms capable of reconstructing two-dimensional images regardless of the size, shape, location or number of inclusions hidden in the examined object. Design/methodology/approach The application consists of a device and measuring sensors, as well as proprietary algorithms for image reconstruction. Ultrasonic transmission tomography makes it possible to analyse processes occurring in an object without interfering with the examined object. The proposed solution uses algorithms based on ray integration, the Fermat principle and deterministic methods. Two applications were developed, one based on C and implemented on the embedded device, while the other application was made in Matlab. Findings Research shows that ultrasonic transmission tomography provides an effective analysis of tested objects in closed tanks. Research limitations/implications In the presented technique, the use of ultrasonic absorption wave has been limited. Nevertheless, the effectiveness of such a solution has been confirmed. Practical implications The presented solution can be used for research and monitoring of technological processes. Originality/value Author’s tomographic system consisting of a measuring system and image reconstruction algorithms.


Symmetry ◽  
2018 ◽  
Vol 10 (10) ◽  
pp. 449
Author(s):  
Chang Liu ◽  
Jun Qiu

In this paper, we propose the symmetric structure of the reconstructed points discretization model to partition and order the subsets of Ordered Subset Expectation Maximum (OSEM) algorithms for image reconstruction and then simplify the calculation of the projection coefficient matrix while satisfying the balancing properties of subsets. The reconstructed points discretization model was utilized to describe the forward and inverse relationships between the reconstructed points and the projection data according to the distance from the point to the ray rather than the intersection length between the square pixel and the ray. This discretization model provides new approaches for improving and constructing the reconstruction algorithms on the basis of the geometry of the model. The experimental results show that the OSEM algorithms based on the reconstructed points discretization model and its geometric symmetry structure can effectively improve the imaging speed and the imaging precision.


Author(s):  
Santosh Bhattacharyya

Three dimensional microscopic structures play an important role in the understanding of various biological and physiological phenomena. Structural details of neurons, such as the density, caliber and volumes of dendrites, are important in understanding physiological and pathological functioning of nervous systems. Even so, many of the widely used stains in biology and neurophysiology are absorbing stains, such as horseradish peroxidase (HRP), and yet most of the iterative, constrained 3D optical image reconstruction research has concentrated on fluorescence microscopy. It is clear that iterative, constrained 3D image reconstruction methodologies are needed for transmitted light brightfield (TLB) imaging as well. One of the difficulties in doing so, in the past, has been in determining the point spread function of the system.We have been developing several variations of iterative, constrained image reconstruction algorithms for TLB imaging. Some of our early testing with one of them was reported previously. These algorithms are based on a linearized model of TLB imaging.


2020 ◽  
Vol 28 (6) ◽  
pp. 829-847
Author(s):  
Hua Huang ◽  
Chengwu Lu ◽  
Lingli Zhang ◽  
Weiwei Wang

AbstractThe projection data obtained using the computed tomography (CT) technique are often incomplete and inconsistent owing to the radiation exposure and practical environment of the CT process, which may lead to a few-view reconstruction problem. Reconstructing an object from few projection views is often an ill-posed inverse problem. To solve such problems, regularization is an effective technique, in which the ill-posed problem is approximated considering a family of neighboring well-posed problems. In this study, we considered the {\ell_{1/2}} regularization to solve such ill-posed problems. Subsequently, the half thresholding algorithm was employed to solve the {\ell_{1/2}} regularization-based problem. The convergence analysis of the proposed method was performed, and the error bound between the reference image and reconstructed image was clarified. Finally, the stability of the proposed method was analyzed. The result of numerical experiments demonstrated that the proposed method can outperform the classical reconstruction algorithms in terms of noise suppression and preserving the details of the reconstructed image.


2017 ◽  
Vol 2017 ◽  
pp. 1-10
Author(s):  
Hsuan-Ming Huang ◽  
Ing-Tsung Hsiao

Background and Objective. Over the past decade, image quality in low-dose computed tomography has been greatly improved by various compressive sensing- (CS-) based reconstruction methods. However, these methods have some disadvantages including high computational cost and slow convergence rate. Many different speed-up techniques for CS-based reconstruction algorithms have been developed. The purpose of this paper is to propose a fast reconstruction framework that combines a CS-based reconstruction algorithm with several speed-up techniques.Methods. First, total difference minimization (TDM) was implemented using the soft-threshold filtering (STF). Second, we combined TDM-STF with the ordered subsets transmission (OSTR) algorithm for accelerating the convergence. To further speed up the convergence of the proposed method, we applied the power factor and the fast iterative shrinkage thresholding algorithm to OSTR and TDM-STF, respectively.Results. Results obtained from simulation and phantom studies showed that many speed-up techniques could be combined to greatly improve the convergence speed of a CS-based reconstruction algorithm. More importantly, the increased computation time (≤10%) was minor as compared to the acceleration provided by the proposed method.Conclusions. In this paper, we have presented a CS-based reconstruction framework that combines several acceleration techniques. Both simulation and phantom studies provide evidence that the proposed method has the potential to satisfy the requirement of fast image reconstruction in practical CT.


2013 ◽  
Vol 2013 ◽  
pp. 1-14
Author(s):  
Joshua Kim ◽  
Huaiqun Guan ◽  
David Gersten ◽  
Tiezhi Zhang

Tetrahedron beam computed tomography (TBCT) performs volumetric imaging using a stack of fan beams generated by a multiple pixel X-ray source. While the TBCT system was designed to overcome the scatter and detector issues faced by cone beam computed tomography (CBCT), it still suffers the same large cone angle artifacts as CBCT due to the use of approximate reconstruction algorithms. It has been shown that iterative reconstruction algorithms are better able to model irregular system geometries and that algebraic iterative algorithms in particular have been able to reduce cone artifacts appearing at large cone angles. In this paper, the SART algorithm is modified for the use with the different TBCT geometries and is tested using both simulated projection data and data acquired using the TBCT benchtop system. The modified SART reconstruction algorithms were able to mitigate the effects of using data generated at large cone angles and were also able to reconstruct CT images without the introduction of artifacts due to either the longitudinal or transverse truncation in the data sets. Algebraic iterative reconstruction can be especially useful for dual-source dual-detector TBCT, wherein the cone angle is the largest in the center of the field of view.


Sign in / Sign up

Export Citation Format

Share Document