skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A regularized weighted least gradient problem for conductivity imaging
We propose and study a method for imaging an approximate electrical conductivity from the magnitude of one interior current density field without any knowledge of the boundary voltage potential. Solely from this interior data, the exact conductivity is impossible to recover as non-unique solutions exist. We propose a method to recover a minimum residual type solution. The method is based on a weighted least gradient problem in the subspace of functions of bounded variations with square integrable traces. We prove existence and uniqueness for a nearby problem, and study the continuous dependence data for a regularized problem. The computational effectiveness and numerica  more » « less
Award ID(s):
1818882
PAR ID:
10093091
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Inverse problems
Volume:
35
Issue:
4
ISSN:
1361-6420
Page Range / eLocation ID:
045006(20pp)
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We consider an inverse shape problem coming from electrical impedance tomography with a Robin transmission condition. In general, a boundary condition of Robin type models corrosion. In this paper, we study two methods for recovering an interior corroded region from electrostatic data. We consider the case where we have small volume and extended regions. For the case where the region has small volume, we will derive an asymptotic expansion of the current gap operator and prove that a MUSIC-type algorithm can be used to recover the region. In the case where one has an extended region, we will show that the regularized factorization method can be used to recover said region. Numerical examples will be presented for both cases in two dimensions in the unit circle. 
    more » « less
  2. Recovering material properties of the subsurface using ground-penetrating radar (GPR) data of finite bandwidth with missing low frequencies and in the presence of strong attenuation is a challenging problem. We have adopted three nonlinear inverse methods for recovering electrical conductivity and permittivity of the subsurface by joining GPR multioffset and electrical resistivity (ER) data acquired at the surface. All of the methods use ER data to constrain the low spatial frequency of the conductivity solution. The first method uses the envelope of the GPR data to exploit low-frequency content in full-waveform inversion and does not assume structural similarities of the material properties. The second method uses cross gradients to manage weak amplitudes in the GPR data by assuming structural similarities between permittivity and conductivity. The third method uses the envelope of the GPR data and the cross gradient of the model parameters. By joining ER and GPR data, exploiting low-frequency content in the GPR data, and assuming structural similarities between the electrical permittivity and conductivity, we are able to recover subsurface parameters in regions where the GPR data have a signal-to-noise ratio close to one. 
    more » « less
  3. Electrical Impedance Tomography (EIT) is a well-known imaging technique for detecting the electrical properties of an object in order to detect anomalies, such as conductive or resistive targets. More specifically, EIT has many applications in medical imaging for the detection and location of bodily tumors since it is an affordable and non-invasive method, which aims to recover the internal conductivity of a body using voltage measurements resulting from applying low frequency current at electrodes placed at its surface. Mathematically, the reconstruction of the internal conductivity is a severely ill-posed inverse problem and yields a poor quality image reconstruction. To remedy this difficulty, at least in part, we regularize and solve the nonlinear minimization problem by the aid of a Krylov subspace-type method for the linear sub problem during each iteration. In EIT, a tumor or general anomaly can be modeled as a piecewise constant perturbation of a smooth background, hence, we solve the regularized problem on a subspace of relatively small dimension by the Flexible Golub-Kahan process that provides solutions that have sparse representation. For comparison, we use a well-known modified Gauss–Newton algorithm as a benchmark. Using simulations, we demonstrate the effectiveness of the proposed method. The obtained reconstructions indicate that the Krylov subspace method is better adapted to solve the ill-posed EIT problem and results in higher resolution images and faster convergence compared to reconstructions using the modified Gauss–Newton algorithm. 
    more » « less
  4. null (Ed.)
    In this paper, we develop a conceptually unified approach for characterizing and determining scattering poles and interior eigenvalues for a given scattering problem. Our approach explores a duality stemming from interchanging the roles of incident and scattered fields in our analysis. Both sets are related to the kernel of the relative scattering operator mapping incident fields to scattered fields, corresponding to the exterior scattering problem for the interior eigenvalues and the interior scattering problem for scattering poles. Our discussion includes the scattering problem for a Dirichlet obstacle where duality is between scattering poles and Dirichlet eigenvalues, and the inhomogeneous scattering problem where the duality is between scattering poles and transmission eigenvalues. Our new characterization of the scattering poles suggests a numerical method for their computation in terms of scattering data for the corresponding interior scattering problem. 
    more » « less
  5. Abstract In this paper, we consider the inverse scattering problem associated with an anisotropic medium with a conductive boundary. We will assume that the corresponding far–field pattern is known/measured and we consider two inverse problems. First, we show that the far–field data uniquely determines the boundary coefficient. Next, since it is known that anisotropic coefficients are not uniquely determined by this data we will develop a qualitative method to recover the scatterer. To this end, we study the so–called monotonicity method applied to this inverse shape problem. This method has recently been applied to some inverse scattering problems but this is the first time it has been applied to an anisotropic scatterer. This method allows one to recover the scatterer by considering the eigenvalues of an operator associated with the far–field operator. We present some simple numerical reconstructions to illustrate our theory in two dimensions. For our reconstructions, we need to compute the adjoint of the Herglotz wave function as an operator mapping intoH1of a small ball. 
    more » « less