Kim Knudsen

Simon Hubmer, Kim Knudsen, Changyou Li et al.

This paper considers the reconstruction problem in Acousto-Electrical Tomography, i.e., the problem of estimating a spatially varying conductivity in a bounded domain from measurements of the internal power densities resulting from different prescribed boundary conditions. Particular emphasis is placed on the limited angle scenario, in which the boundary conditions are supported only on a part of the boundary. The reconstruction problem is formulated as an optimization problem in a Hilbert space setting and solved using Landweber iteration. The resulting algorithm is implemented numerically in two spatial dimensions and tested on simulated data. The results quantify the intuition that features close to the measurement boundary are stably reconstructed and features further away are less well reconstructed. Finally, the ill-posedness of the limited angle problem is quantified numerically using the singular value decomposition of the corresponding linearized problem.

Mirza Karamehmedović, Adrian Kirkeby, Kim Knudsen

We consider the multi-frequency inverse source problem for the scalar Helmholtz equation in the plane. The goal is to reconstruct the source term in the equation from measurements of the solution on a surface outside the support of the source. We study the problem in a certain finite dimensional setting: From measurements made at a finite set of frequencies we uniquely determine and reconstruct sources in a subspace spanned by finitely many Fourier-Bessel functions. Further, we obtain a constructive criterion for identifying a minimal set of measurement frequencies sufficient for reconstruction, and under an additional, mild assumption, the reconstruction method is shown to be stable. Our analysis is based on a singular value decomposition of the source-to-measurement forward operators and the distribution of positive zeros of the Bessel functions of the first kind. The reconstruction method is implemented numerically and our theoretical findings are supported by numerical experiments.

Rasmus Dalgas Kongskov, Yiqiu Dong, Kim Knudsen

In inverse problems, prior information and a priori-based regularization techniques play important roles. In this paper, we focus on image restoration problems, especially on restoring images whose texture mainly follow one direction. In order to incorporate the directional information, we propose a new directional total generalized variation (DTGV) functional, which is based on total generalized variation (TGV) by Bredies \textit{et al}. [SIAM J. Imaging Sci., 3 (2010)]. After studying the mathematical properties of DTGV, we utilize it as regularizer and propose the L$^2$-DTGV variational model for solving image restoration problems. Due to the requirement of the directional information in DTGV, we give a direction estimation algorithm, and then apply a primal-dual algorithm to solve the minimization problem. Experimental results show the effectiveness of the proposed method for restoring the directional images. In comparison with isotropic regularizers like total variation and TGV, the improvement of texture preservation and noise removal is significant.

Henrik Garde, Kim Knudsen

This paper focuses on prior information for improved sparsity reconstruction in electrical impedance tomography with partial data, i.e. data measured only on subsets of the boundary. Sparsity is enforced using an $\ell_1$ norm of the basis coefficients as the penalty term in a Tikhonov functional, and prior information is incorporated by applying a spatially distributed regularization parameter. The resulting optimization problem allows great flexibility with respect to the choice of measurement boundaries and incorporation of prior knowledge. The problem is solved using a generalized conditional gradient method applying soft thresholding. Numerical examples show that the addition of prior information in the proposed algorithm gives vastly improved reconstructions even for the partial data problem. The method is in addition compared to a total variation approach.

Henrik Garde, Kim Knudsen

In electrical impedance tomography the electrical conductivity inside a physical body is computed from electro-static boundary measurements. The focus of this paper is to extend recent result for the 2D problem to 3D. Prior information about the sparsity and spatial distribution of the conductivity is used to improve reconstructions for the partial data problem with Cauchy data measured only on a subset of the boundary. A sparsity prior is enforced using the $\ell_1$ norm in the penalty term of a Tikhonov functional, and spatial prior information is incorporated by applying a spatially distributed regularization parameter. The optimization problem is solved numerically using a generalized conditional gradient method with soft thresholding. Numerical examples show the effectiveness of the suggested method even for the partial data problem with measurements affected by noise.