Linear and Nonlinear Inverse Problems with Practical Applications

Jennifer L. Mueller and Samuli Siltanen

This page contains the computational Matlab files related to the book, which you can purchase in the SIAM Bookstore.

Using the codes requires some familiarity with the Matlab software. Also, some routines make use of Matlab toolboxes (image processing toolbox in the case of X-ray tomography codes, optimization toolbox for total variation regularization for the 1D deconvolution, and PDE toolbox in the D-bar method files), which you should have available in those cases.

These Matlab files will help you get started with computational inversion.

Part I: Linear Inverse Problems

We consider several examples of linear inverse problems: X_ray tomography, one-dimensional deconvolution, image deblurring, and backward heat equation. For each example, we show how to construct measurement matrices and how to reconstruct unknown functions from indirect data corrupted by simulated noise.

One-dimensional deconvolution

The 1D convolution model has a lot of applications. We show how to simulate data and invert it using regularized methods on the following page:

One-dimensional deconvolution

Backward heat equation

Going back in time in heat propagation is a notoriously ill-posed problem. The following page contains algorithms for robust computational inversion:

Backward heat equation

X-ray tomography

We start with quite explicit computations based on actually constructing the measurement matrix A for a tomographic application:

X-ray tomography with matrices

Once you are familiar with the above matrix-based approach, it's time to get more practical. In large-scale imaging, constructing the measurement matrix A is out of the question because it is too large. Instead, we implement reconstruction methods that only need to apply the matrix A and its transpose to given vectors. This in turn can be done using Matlab's Radon.m routine.

Matrix-free X-ray tomography

Now the stuff gets really interesting: we consider tomographic problems with sparsely sampled data! You can modify the files further to examine limited-angle tomography as well.

Matrix-free X-ray tomography with sparse data

Image deblurring

You can apply Tikhonov regularization and approximate total variation regularization to two-dimensional deconvolution problems using the routines on the following page:

Image deblurring page

Part II: Nonlinear Inverse Problems

The guiding example of Part II is electrical impedance tomography. We show how to simulate the Dirichlet-to-Neumann map using the Finite Element Method and how to recover the conductivity using the D-bar method.

We discuss the following two examples:

EIT with the D-bar method: smooth and radial case

EIT with the D-bar method: discontinuous heart-and-lungs phantom

Donate · Contact Us · Site Map · Join SIAM · My Account
Facebook Twitter Youtube linkedin google+