BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Samuli Siltanen

Samuli Siltanen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath.NAphysics.med-pheess.IVmath.FAmath.PRphysics.data-anmath-phmath.MPmath.STphysics.comp-phStatistics TheoryComputer VisionMachine Learningmath.OCmath.SPNumerical Analysis

Catalog footprint

What is connected

23works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Learning a microlocal prior for limited-angle tomography

Limited-angle tomography is a highly ill-posed linear inverse problem. It arises in many applications, such as digital breast tomosynthesis. Reconstructions from limited-angle data typically suffer from severe stretching of features along the central direction of projections, leading to poor separation between slices perpendicular to the central direction. A new method is introduced, based on machine learning and geometry, producing an estimate for interfaces between regions of different X-ray attenuation. The estimate can be presented on top of the reconstruction, indicating more reliably the true form and extent of features. The method uses directional edge detection, implemented using complex wavelets and enhanced with morphological operations. By using machine learning, the visible part of the wavefront set is first extracted and then extended to the full domain, filling in the parts of the wavefront set that would otherwise be hidden due to the lack of measurement directions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Random tree Besov priors -- Towards fractal imaging

We propose alternatives to Bayesian a priori distributions that are frequently used in the study of inverse problems. Our aim is to construct priors that have similar good edge-preserving properties as total variation or Mumford-Shah priors but correspond to well defined infinite-dimensional random variables, and can be approximated by finite-dimensional random variables. We introduce a new wavelet-based model, where the non zero coefficient are chosen in a systematic way so that prior draws have certain fractal behaviour. We show that realisations of this new prior take values in some Besov spaces and have singularities only on a small set $τ$ that has a certain Hausdorff dimension. We also introduce an efficient algorithm for calculating the MAP estimator, arising from the the new prior, in denoising problem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deep neural networks for inverse problems with pseudodifferential operators: an application to limited-angle tomography

We propose a novel convolutional neural network (CNN), called $Ψ$DONet, designed for learning pseudodifferential operators ($Ψ$DOs) in the context of linear inverse problems. Our starting point is the Iterative Soft Thresholding Algorithm (ISTA), a well-known algorithm to solve sparsity-promoting minimization problems. We show that, under rather general assumptions on the forward operator, the unfolded iterations of ISTA can be interpreted as the successive layers of a CNN, which in turn provides fairly general network architectures that, for a specific choice of the parameters involved, allow to reproduce ISTA, or a perturbation of ISTA for which we can bound the coefficients of the filters. Our case study is the limited-angle X-ray transform and its application to limited-angle computed tomography (LA-CT). In particular, we prove that, in the case of LA-CT, the operations of upscaling, downscaling and convolution, which characterize our $Ψ$DONet and most deep learning schemes, can be exactly determined by combining the convolutional nature of the limited angle X-ray transform and basic properties defining an orthogonal wavelet system. We test two different implementations of $Ψ$DONet on simulated data from limited-angle geometry, generated from the ellipse data set. Both implementations provide equally good and noteworthy preliminary results, showing the potential of the approach we propose and paving the way to applying the same idea to other convolutional operators which are $Ψ$DOs or Fourier integral operators.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Joint Reconstruction in Low Dose Multi-Energy CT

Multi-energy CT takes advantage of the non-linearly varying attenuation properties of elemental media with respect to energy, enabling more precise material identification than single-energy CT. The increased precision comes with the cost of a higher radiation dose. A straightforward way to lower the dose is to reduce the number of projections per energy, but this makes tomographic reconstruction more ill-posed. In this paper, we propose how this problem can be overcome with a combination of a regularization method that promotes structural similarity between images at different energies and a suitably selected low-dose data acquisition protocol using non-overlapping projections. The performance of various joint regularization models is assessed with both simulated and experimental data, using the novel low-dose data acquisition protocol. Three of the models are well-established, namely the joint total variation, the linear parallel level sets and the spectral smoothness promoting regularization models. Furthermore, one new joint regularization model is introduced for multi-energy CT: a regularization based on the structure function from the structural similarity index. The findings show that joint regularization outperforms individual channel-by-channel reconstruction. Furthermore, the proposed combination of joint reconstruction and non-overlapping projection geometry enables significant reduction of radiation dose.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sparse dynamic tomography. A shearlet-based approach for iodine perfusion in plant stems

In this paper we propose a motion-aware variational approach to reconstruct moving objects from sparse dynamic data. The motivation of this work stems from X-ray imaging of plants perfused with a liquid contrast agent, aimed at increasing the contrast of the images and studying the phloem transport in plants over time. The key idea of our approach is to deploy 3D shearlets as a space-temporal prior, treating time as the third dimension. The rationale behind this model is that a continuous evolution of a cartoon-like object is well suited for the use of 3D shearlets. We provide a basic mathematical analysis of the variational model for the image reconstruction. The numerical minimization is carried out with primal-dual scheme coupled with an automated choice of regularization parameter. We test our model on different measurement setups: a simulated phantom especially designed to resemble a plant stem, with spreading points to simulate a spreading contrast agent; a measured agarose gel phantom to demonstrate iodide diffusion and geometry prior to imaging living sample; a measured living tree grown \textit{in vitro} and perfused with a liquid sugar-iodine-mix. The results, compared against a 2D static model, show that our approach provides reconstructions that capture well the time dynamic of the contrast agent onset and are encouraging to develop microCT as a tool to study phloem transport using iodine tracer.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Analysis of regularized inversion of data corrupted by white Gaussian noise

Tikhonov regularization is studied in the case of linear pseudodifferential operator as the forward map and additive white Gaussian noise as the measurement error. The measurement model for an unknown function $u(x)$ is \begin{eqnarray*} m(x) = Au(x) + δ\hspace{.2mm}\varepsilon(x), \end{eqnarray*} where $δ>0$ is the noise magnitude. If $\varepsilon$ was an $L^2$-function, Tikhonov regularization gives an estimate \begin{eqnarray*} T_α(m) = \text{argmin}_{u\in H^r}\big\{\|A u-m\|_{L^2}^2+ α\|u\|_{H^r}^2 \big\}\end{eqnarray*} for $u$ where $α=α(δ)$ is the regularization parameter. Here penalization of the Sobolev norm $ \|u\|_{H^r}$ covers the cases of standard Tikhonov regularization ($r=0$) and first derivative penalty ($r=1$). Realizations of white Gaussian noise are almost never in $L^2$, but do belong to $H^s$ with probability one if $s<0$ is small enough. A modification of Tikhonov regularization theory is presented, covering the case of white Gaussian measurement noise. Furthermore, the convergence of regularized reconstructions to the correct solution as $δ\rightarrow 0$ is proven in appropriate function spaces using microlocal analysis. The convergence of the related finite-dimensional problems to the infinite-dimensional problem is also analysed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Posterior consistency and convergence rates for Bayesian inversion with hypoelliptic operators

Bayesian approach to inverse problems is studied in the case where the forward map is a linear hypoelliptic pseudodifferential operator and measurement error is additive white Gaussian noise. The measurement model for an unknown Gaussian random variable $U(x,ω)$ is \begin{eqnarray*} M(y,ω) = A(U(x,ω) )+ δ\hspace{.2mm}\mathcal{E}(y,ω), \end{eqnarray*} where $A$ is a finitely many times smoothing linear hypoelliptic operator and $δ>0$ is the noise magnitude. The covariance operator $C_U$ of $U$ is $2r$ times smoothing, self-adjoint, injective and elliptic pseudodifferential operator. If $\mathcal{E}$ was taking values in $L^2$ then in Gaussian case solving the conditional mean (and maximum a posteriori) estimate is linked to solving the minimisation problem \begin{eqnarray*} T_δ(M) = \text{argmin}_{u\in H^r} \big\{\|A u-m\|_{L^2}^2+ δ^2\|C_U^{-1/2}u\|_{L^2}^2 \big\}. \end{eqnarray*} However, Gaussian white noise does not take values in $L^2$ but in $H^{-s}$ where $s>0$ is big enough. A modification of the above approach to solve the inverse problem is presented, covering the case of white Gaussian measurement noise. Furthermore, the convergence of conditional mean estimate to the correct solution as $δ\rightarrow 0$ is proven in appropriate function spaces using microlocal analysis. Also the contraction of the confidence regions is studied.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The D-bar Method for Diffuse Optical Tomography: a computational study

The D-bar method at negative energy is numerically implemented. Using the method we are able to numerically reconstruct potentials and investigate exceptional points at negative energy. Subsequently, applying the method to Diffusive Optical Tomography, a new way of reconstructing the diffusion coefficient from the associated Complex Geometrics Optics solution is suggested and numerically validated.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tomographic X-ray data of a lotus root filled with attenuating objects

This is the documentation of the tomographic X-ray data of a lotus root, filled with four different attenuating objects, of different sizes. Data are available at www.fips.fi/dataset.php, and can be freely used for scientific purposes with appropriate references to them, and to this document in http://arxiv.org/arXiv. The data set consists of (1) the X-ray sinogram of a single 2D slice of the lotus root with two different resolutions and (2) the corresponding measurement matrices modeling the linear operation of the X-ray transform. Each of these sinograms was obtained from a measured 360-projection fan-beam sinogram by down-sampling and taking logarithms. The original (measured) sinogram is also provided in its original form and resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Direct Reconstruction Method for Anisotropic Electrical Impedance Tomography

A novel computational, non-iterative and noise-robust reconstruction method is introduced for the planar anisotropic inverse conductivity problem. The method is based on bypassing the unstable step of the reconstruction of the values of the isothermal coordinates on the boundary of the domain. Non-uniqueness of the inverse problem is dealt with by recovering the unique isotropic conductivity that can be achieved as a deformation of the measured anisotropic conductivity by \emph{isothermal coordinates}. The method shows how isotropic D-bar reconstruction methods have produced reasonable and informative reconstructions even when used on EIT data known to come from anisotropic media, and when the boundary shape is not known precisely. Furthermore, the results pave the way for regularized anisotropic EIT. Key aspects of the approach involve D-bar methods and inverse scattering theory, complex geometrical optics solutions, and quasi-conformal mapping techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Hybrid Segmentation and D-bar Method for Electrical Impedance Tomography

The Regularized D-bar method for Electrical Impedance Tomography provides a rigorous mathematical approach for solving the full nonlinear inverse problem directly, i.e. without iterations. It is based on a low-pass filtering in the (nonlinear) frequency domain. However, the resulting D-bar reconstructions are inherently smoothed leading to a loss of edge distinction. In this paper, a novel approach that combines the rigor of the D-bar approach with the edge-preserving nature of Total Variation regularization is presented. The method also includes a data-driven contrast adjustment technique guided by the key functions (CGO solutions) of the D-bar method. The new TV-Enhanced D-bar Method produces reconstructions with sharper edges and improved contrast while still solving the full nonlinear problem. This is achieved by using the TV-induced edges to increase the truncation radius of the scattering data in the nonlinear frequency domain thereby increasing the radius of the low pass filter. The algorithm is tested on numerically simulated noisy EIT data and demonstrates significant improvements in edge preservation and contrast which can be highly valuable for absolute EIT imaging.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-resolution parameter choice method for total variation regularized tomography

A computational method is introduced for choosing the regularization parameter for total variation (TV) regularization. The approach is based on computing reconstructions at a few different resolutions and various values of regularization parameter. The chosen parameter is the smallest one resulting in approximately discretization-invariant TV norms of the reconstructions. The method is tested with X-ray tomography data measured from a walnut and compared to the S-curve method. The proposed method seems to automatically adapt to the desired resolution and noise level, and it yields useful results in the tests. The results are comparable to those of the S-curve method; however, the S-curve method needs a priori information about the sparsity of the unknown, while the proposed method does not need any a priori information (apart from the choice of a desired resolution). Mathematical analysis is presented for (partial) understanding of the properties of the proposed parameter choice method. It is rigorously proven that the TV norms of the reconstructions converge with any choice of regularization parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Positive-energy D-bar method for acoustic tomography: a computational study

A new computational method for reconstructing a potential from the Dirichlet-to-Neumann map at positive energy is developed. The method is based on D-bar techniques and it works in absence of exceptional points -- in particular, if the potential is small enough compared to the energy. Numerical tests reveal exceptional points for perturbed, radial potentials. Reconstructions for several potentials are computed using simulated Dirichlet-to-Neumann maps with and without added noise. The new reconstruction method is shown to work well for energy values between $10^{-5}$ and $5$, smaller values giving better results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tomographic X-ray data of a walnut

This is the documentation of the tomographic X-ray data of a walnut made available at http://www.fips.fi/dataset.php . The data can be freely used for scientific purposes with appropriate references to the data and to this document in arXiv. The data set consists of (1) the X-ray sinogram of a single 2D slice of the walnut with three different resolutions and (2) the corresponding measurement matrices modeling the linear operation of the X-ray transform. Each of these sinograms was obtained from a measured 120-projection fan-beam sinogram by down-sampling and taking logarithms. The original (measured) sinogram is also provided in its original form and resolution. In addition, a larger set of 1200 projections of the same walnut was measured and a high-resolution filtered back-projection reconstruction was computed from this data; both the sinogram and the FBP reconstruction are included in the data set, the latter serving as a ground truth reconstruction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Estimation of conductivity changes in a region of interest with electrical impedance tomography

This paper proposes a novel approach to reconstruct changes in a target conductivity from electrical impedance tomography measurements. As in the conventional difference imaging, the reconstruction of the conductivity change is based on electrical potential measurements from the exterior boundary of the target before and after the change. In this paper, however, images of the conductivity before and after the change are reconstructed simultaneously based on the two data sets. The key feature of the approach is that the conductivity after the change is parameterized as a linear combination of the initial state and the change. This allows for modeling independently the spatial characteristics of the background conductivity and the change of the conductivity - by separate regularization functionals. The approach also allows in a straightforward way the restriction of the conductivity change to a localized region of interest inside the domain. While conventional difference imaging reconstruction is based on a global linearization of the observation model, the proposed approach amounts to solving a non-linear inverse problem. The feasibility of the proposed reconstruction method is tested experimentally and with a simulation which demonstrates a potential new medical application of electrical impedance tomography: imaging of vocal folds in voice loading studies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear Fourier analysis for discontinuous conductivities: computational results

Two reconstruction methods of Electrical Impedance Tomography (EIT) are numerically compared for nonsmooth conductivities in the plane based on the use of complex geometrical optics (CGO) solutions to D-bar equations involving the global uniqueness proofs for Calderón problem exposed in [Nachman; Annals of Mathematics 143, 1996] and [Astala and Päivärinta; Annals of Mathematics 163, 2006]: the Astala-Päivärinta theory-based "low-pass transport matrix method" implemented in [Astala et al.; Inverse Problems and Imaging 5, 2011] and the "shortcut method" which considers ingredients of both theories. The latter method is formally similar to the Nachman theory-based regularized EIT reconstruction algorithm studied in [Knudsen, Lassas, Mueller and Siltanen; Inverse Problems and Imaging 3, 2009] and several references from there. New numerical results are presented using parallel computation with size parameters larger than ever, leading mainly to two conclusions as follows. First, both methods can approximate piecewise constant conductivities better and better as the cutoff frequency increases, and there seems to be a Gibbs-like phenomenon producing ringing artifacts. Second, the transport matrix method loses accuracy away from a (freely chosen) pivot point located outside of the object to be studied, whereas the shortcut method produces reconstructions with more uniform quality.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear Inversion from Partial EIT Data: Computational Experiments

Electrical impedance tomography (EIT) is a non-invasive imaging method in which an unknown physical body is probed with electric currents applied on the boundary, and the internal conductivity distribution is recovered from the measured boundary voltage data. The reconstruction task is a nonlinear and ill-posed inverse problem, whose solution calls for special regularized algorithms, such as D-bar methods which are based on complex geometrical optics solutions (CGOs). In many applications of EIT, such as monitoring the heart and lungs of unconscious intensive care patients or locating the focus of an epileptic seizure, data acquisition on the entire boundary of the body is impractical, restricting the boundary area available for EIT measurements. An extension of the D-bar method to the case when data is collected only on a subset of the boundary is studied by computational simulation. The approach is based on solving a boundary integral equation for the traces of the CGOs using localized basis functions (Haar wavelets). The numerical evidence suggests that the D-bar method can be applied to partial-boundary data in dimension two and that the traces of the partial data CGOs approximate the full data CGO solutions on the available portion of the boundary, for the necessary small $k$ frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Novikov-Veselov Equation: Theory and Computation

Recent progress in the theory and computation for the Novikov-Veselov (NV) equation is reviewed with initial potentials decaying at infinity, focusing mainly on the zero-energy case. The inverse scattering method for the zero-energy NV equation is presented in the context of Manakov triples, treating initial data of conductivity type rigorously. Special closed-form solutions are presented, including multisolitons, ring solitons, and breathers. The computational inverse scattering method is used to study zero-energy exceptional points and the relationship between supercritical, critical, and subcritical potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exceptional circles of radial potentials

A nonlinear scattering transform is studied for the two-dimensional Schrodinger equation at zero energy with a radial potential. First explicit examples are presented, both theoretically and computationally, of potentials with nontrivial singularities in the scattering transform. The singularities arise from non-uniqueness of the complex geometric optics solutions that define the scattering transform. The values of the complex spectral parameter at which the singularities appear are called exceptional points. The singularity formation is closely related to the fact that potentials of conductivity type are critical in the sense of Murata.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reconstructing conductivities with boundary corrected D-bar method

The aim of electrical impedance tomography is to form an image of the conductivity distribution inside an unknown body using electric boundary measurements. The computation of the image from measurement data is a non-linear ill-posed inverse problem and calls for a special regularized algorithm. One such algorithm, the so-called D-bar method, is improved in this work by introducing new computational steps that remove the so far necessary requirement that the conductivity should be constant near the boundary. The numerical experiments presented suggest two conclusions. First, for most conductivities arising in medical imaging, it seems the previous approach of using a best possible constant near the boundary is sufficient. Second, for conductivities that have high contrast features at the boundary, the new approach produces reconstructions with smaller quantitative error and with better visual quality.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Novikov-Veselov Equation and the Inverse Scattering Method, Part I: Analysis

The Novikov-Veselov (NV) equation is a (2+1)-dimensional nonlinear evolution equation that generalizes the (1+1)-dimensional Korteweg-deVries (KdV) equation. Solution of the NV equation using the inverse scattering method has been discussed in the literature, but only formally (or with smallness assumptions in case of nonzero energy) because of the possibility of exceptional points, or singularities in the scattering data. In this work, absence of exceptional points is proved at zero energy for evolutions with compactly supported, smooth and rotationally symmetric initial data of the conductivity type: $q_0=γ^{-1/2}Δγ^{1/2}$ with a strictly positive function $γ$. The inverse scattering evolution is shown to be well-defined, real-valued, and preserving conductivity-type. There is no smallness assumption on the initial data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Determining electrical and heat transfer parameters using coupled boundary measurements

Let $Ω\subset\R^n$, $n\ge 3$, be a smooth bounded domain and consider a coupled system in $Ω$ consisting of a conductivity equation $\nabla \cdot γ(x) \nabla u(t,x)=0$ and an anisotropic heat equation $κ^{-1}(x)\partial_tψ(t,x)=\nabla\cdot (A(x)\nabla ψ(t,x))+(γ\nabla u(t,x))\cdot \nabla u(t,x), \quad t\ge 0$. It is shown that the coefficients $γ$, $κ$ and $A=(a_{jk})$ are uniquely determined from the knowledge of the boundary map $u|_{\partialΩ}\mapsto ν\cdot A\nabla ψ|_{\partialΩ}$, where $ν$ is the unit outer normal to $\partialΩ$. The coupled system models the following physical phenomenon. Given a fixed voltage distribution, maintained on the boundary $\partialΩ$, an electric current distribution appears inside $Ω$. The current in turn acts as a source of heat inside $Ω$, and the heat flows out of the body through the boundary. The boundary measurements above then correspond to the map taking a voltage distribution on the boundary to the resulting heat flow through the boundary. The presented mathematical results suggest a new hybrid diffuse imaging modality combining electrical prospecting and heat transfer-based probing.

0 citations0 reviewsNo code linkedOpen access