Inverse Problems & Imaging
April 2020 , Volume 14 , Issue 2
Select all articles
In this paper, we present a patch based weighted means filter for removing an impulse noise by adapting the fundamental idea of the non-local means filter to the random-valued impulse noise. Our approach is to give a weight to a pixel in order to evaluate the probability that the pixel is contaminated by the impulse noise, which we call Reliable Weight of the pixel. With the help of the Reliable Weights we introduce the similarity function to measure the similarity among patches of the image contaminated by a random impulse noise. It turns out that the similarity function has significant anti impulse noise interference ability. We then incorporate the Reliable Weights and the similarity function into a filter designed to remove the random impulse noise. Under suitable conditions, we establish two convergence theorems to demonstrate that our method is feasible. Simulation results confirm that our filter is competitive compared to recently proposed methods.
We consider the inverse problem of denoising an image where each point (pixel) is an element of a target set, which we refer to as a target-valued image. The target sets considered are either (ⅰ) a closed convex set of Euclidean space or (ⅱ) a closed subset of the sphere. The energy for the denoising problem consists of an
We establish the exact recovery guarantees for a class of Riemannian optimization methods based on the embedded manifold of low rank matrices for matrix completion. Assume
via the resampled Riemannian gradient descent initialization. The analysis of the new initialization procedure relies on an asymmetric restricted isometry property of the sampling operator and the curvature of the low rank matrix manifold. Numerical simulation shows that the algorithms are able to recover a low rank matrix from nearly the minimum number of measurements.
Low rank approximation has been extensively studied in the past. It is most suitable to reproduce rectangular like structures in the data. In this work we introduce a generalization using "shifted" rank-
These kind of shifts naturally appear in applications, where an object
The main difficulty of the above stated problem lies in finding a suitable shift vector
We validate our approach in several numerical experiments on different kinds of data. We compare the technique to shift-invariant dictionary learning algorithms. Furthermore, we provide examples from application including object segmentation in non-destructive testing and seismic exploration as well as object tracking in video processing.
Region-of-interest computed tomography (ROI CT) aims at reconstructing a region within the field of view by using only ROI-focused projections. The solution of this inverse problem is challenging and methods of tomographic reconstruction that are designed to work with full projection data may perform poorly or fail when applied to this setting. In this work, we study the ROI CT problem in the presence of measurement noise and formulate the reconstruction problem by relaxing data fidelity and consistency requirements. Under the assumption of a robust width prior that provides a form of stability for data satisfying appropriate sparsity-inducing norms, we derive reconstruction performance guarantees and controllable error bounds. Based on this theoretical setting, we introduce a novel iterative reconstruction algorithm from ROI-focused projection data that is guaranteed to converge with controllable error while satisfying predetermined fidelity and consistency tolerances. Numerical tests on experimental data show that our algorithm for ROI CT is competitive with state-of-the-art methods especially when the ROI radius is small.
In the context of international nuclear safeguards, the International Atomic Energy Agency (IAEA) has recently approved passive gamma emission tomography (PGET) as a method for inspecting spent nuclear fuel assemblies (SFAs). The PGET instrument is essentially a single photon emission computed tomography (SPECT) system that allows the reconstruction of axial cross-sections of the emission map of an SFA. The fuel material heavily self-attenuates its gamma-ray emissions, so that correctly accounting for the attenuation is a critical factor in producing accurate images. Due to the nature of the inspections, it is desirable to use as little a priori information as possible about the fuel, including the attenuation map, in the reconstruction process. Current reconstruction methods either do not correct for attenuation, assume a uniform attenuation throughout the fuel assembly, or assume an attenuation map based on an initial filtered back-projection reconstruction. We propose a method to simultaneously reconstruct the emission and attenuation maps by formulating the reconstruction as a constrained minimization problem with a least squares data fidelity term and regularization terms. Using simulated data, we show that our approach produces clear reconstructions which allow for a highly reliable classification of spent, missing, and fresh fuel rods.
Hyperspectral image (HSI) super-resolution is a technique to improve the spatial resolution of a HSI for better visual perception and down stream applications. This is a very ill-posed inverse problem and is often solved by fusing the low-resolution (LR) HSI with a high-resolution (HR) multispectral image (MSI). It is more challenging for blind HSI super-resolution, i.e., when the spatial degradation operators are completely unknown. In this paper, we propose a novel sparse tensor factorization model for the task of blind HSI super-resolution using the spatial non-local self-similarity and spectral global correlation of HSIs. Image clustering method is employed to collect some similar 3D cubes of HSIs which can be formed as some 4D image clusters with high correlation. We conduct cluster wise computation to not only save computation time but also to introduce a non-local regularity originated from the redundancy of cubes. By using the sparsity of tensor decomposition and the low-rank in non-local self-similarity direction underlying 4D similar clusters, we design a sparse tensor regularization term, which preserves the spatial-spectral structural correlation of HSIs. In addition, we present a proximal alternating direction method of multipliers (ADMM) based algorithm to efficiently solve the proposed model. Numerical experiments demonstrate that the proposed model outperforms many state-of-the-art HSI super-resolution methods.
This paper addresses the problem of identifying impenetrable obstacles in a Kirchhoff-Love infinite plate from multistatic near-field data. The Linear Sampling Method is introduced in this context. We firstly prove a uniqueness result for such an inverse problem. We secondly provide the classical theoretical foundation of the Linear Sampling Method. We lastly show the feasibility of the method with the help of numerical experiments.
We model electrical impedance tomography (EIT) based on the minimum energy principle. It results in a constrained minimization problem in terms of current density. The new formulation is proved to have a unique solution within appropriate function spaces. By characterizing its solution with the Lagrange multiplier method, we relate the new formulation to the so-called shunt model and the complete electrode model (CEM) of EIT. Based on the new formulation, we also propose a new numerical method to solve the forward problem of EIT. The new solver is formulated in terms of current. It was shown to give similar results to that of the traditional finite element method, with simulations on a 2D EIT model.
Add your name and e-mail address to receive news of forthcoming issues of this journal:
[Back to Top]