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

Ville Kolehmainen

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

ResearcherUnclaimed source record
math.NAphysics.med-pheess.IVNumerical Analysisphysics.comp-phmath.AP

Catalog footprint

What is connected

10works
6topics
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

10 published item(s)

preprint2026arXiv

Spotlight inversion by orthogonal projections

Many computational problems involve solving a linear system of equations, although only a subset of the entries of the solution are needed. In inverse problems, where the goal is to estimate unknown parameters from indirect noisy observations, it is not uncommon that the forward model linking the observed variables to the unknowns depends on variables that are not of primary interest, often referred to as nuisance parameters. In this article, we consider linear problems, and propose a novel projection technique to eliminate, or at least mitigate, the contribution of the nuisance parameters in the model. We refer to this approach as spotlight inversion, as it allows to focus on only the portion of primary interest of the unknown parameter vector, leaving the uninteresting part in the shadow. The viability of the approach is illustrated with two computed examples, one where it works as model reduction for a finite element approximation of an elliptic PDE, the other amounting to local fanbeam X-ray tomography, spotlighting the region of interest that is part of the full target.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Embedded quantitative MRI T1rho mapping using non-linear primal-dual proximal splitting

Quantitative MRI (qMRI) methods allow reducing the subjectivity of clinical MRI by providing numerical values on which diagnostic assessment or predictions of tissue properties can be based. However, qMRI measurements typically take more time than anatomical imaging due to requiring multiple measurements with varying contrasts for, e.g., relaxation time mapping. To reduce the scanning time, undersampled data may be combined with compressed sensing reconstruction techniques. Typical CS reconstructions first reconstruct a complex-valued set of images corresponding to the varying contrasts, followed by a non-linear signal model fit to obtain the parameter maps. We propose a direct, embedded reconstruction method for T1rho mapping. The proposed method capitalizes on a known signal model to directly reconstruct the desired parameter map using a non-linear optimization model. The proposed reconstruction method also allows directly regularizing the parameter map of interest, and greatly reduces the number of unknowns in the reconstruction. We test the proposed model using a simulated radially sampled data from a 2D phantom and 2D cartesian ex vivo measurements of a mouse kidney specimen. We compare the embedded reconstruction model to two CS reconstruction models, and in the cartesian test case also iFFT. The proposed, embedded model outperformed the reference methods on both test cases, especially with higher acceleration factors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

3D EIT Reconstructions from Electrode Data using Direct Inversion D-bar and Calderon Methods

The first numerical implementation of a D-bar method in 3D using electrode data is presented. Results are compared to Calderón's method as well as more common TV and smoothness regularization-based methods. D-bar methods are based on tailor-made non-linear Fourier transforms involving the measured current and voltage data. Low-pass filtering in the non-linear Fourier domain is used to stabilize the reconstruction process. D-bar methods have shown great promise in 2D for providing robust real-time absolute and time-difference conductivity reconstructions but have yet to be used on practical electrode data in 3D, until now. Results are presented for simulated data for conductivity and permittivity with disjoint non-radially symmetric targets on spherical domains and noisy voltage data. The 3D D-bar and Calderón methods are demonstrated to provide comparable quality to their 2D CGO counterparts, and hold promise for real-time reconstructions.

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

Parameter selection in dynamic contrast-enhanced magnetic resonance tomography

In this work we consider the image reconstruction problem of sparsely sampled dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI). DCE-MRI is a technique for acquiring a series of MR images before, during and after intravenous contrast agent administration, and it is used to study microvascular structure and perfusion. To overcome the ill-posedness of the related spatio-temporal inverse problem, we use regularization. In regularization one of the main problems is how to determine the regularization parameter which controls the balance between data fitting term and regularization term. Most methods for selecting this parameter require the computation of a large number of estimates even in stationary problems. In dynamic imaging, the parameter selection is even more time consuming since separate regularization parameters are needed for the spatial and temporal regularization functionals. In this work, we study the possibility of using the S-curve with DCE-MR data. We select the spatial regularization parameter using the S-curve, leaving the temporal regularization parameter as the only free parameter in the reconstruction problem. In this work, the temporal regularization parameter is selected manually by computing reconstructions with several values of the temporal regularization parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Temporal Huber regularization for DCE-MRI

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is used to study microvascular structure and tissue perfusion. In DCE-MRI a bolus of gadolinium based contrast agent is injected into the blood stream and spatiotemporal changes induced by the contrast agent flow are estimated from a time series of MRI data. Sufficient time resolution can often only be obtained by using an imaging protocol which produces undersampled data for each image in the time series. This has led to the popularity of compressed sensing based image reconstruction approaches, where all the images in the time series are reconstructed simultaneously, and temporal coupling between the images is introduced into the problem by a sparsity promoting regularization functional. We propose the use of Huber penalty for temporal regularization in DCE-MRI, and compare it to total variation, total generalized variation and smoothness based temporal regularization models. We also study the effect of spatial regularization to the reconstruction and compare the reconstruction accuracy with different temporal resolutions due to varying undersampling. The approaches are tested using simulated and experimental radial golden angle DCE-MRI data from a rat brain specimen. The results indicate that Huber regularization produces similar reconstruction accuracy with the total variation based models, but the computation times are significantly faster.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Approximate Marginalization of Absorption and Scattering in Fluorescence Diffuse Optical Tomography

In fluorescence diffuse optical tomography (fDOT), the reconstruction of the fluorophore concentration inside the target body is usually carried out using a normalized Born approximation model where the measured fluorescent emission data is scaled by measured excitation data. One of the benefits of the model is that it can tolerate inaccuracy in the absorption and scattering distributions that are used in the construction of the forward model to some extent. In this paper, we employ the recently proposed Bayesian approximation error approach to fDOT for compensating for the modeling errors caused by the inaccurately known optical properties of the target in combination with the normalized Born approximation model. The approach is evaluated using a simulated test case with different amount of error in the optical properties. The results show that the Bayesian approximation error approach improves the tolerance of fDOT imaging against modeling errors caused by inaccurately known absorption and scattering of the target.

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
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
preprint2006arXiv

The inverse conductivity problem with an imperfectly known boundary in three dimensions

We consider the inverse conductivity problem in a strictly convex domain whose boundary is not known. Usually the numerical reconstruction from the measured current and voltage data is done assuming the domain has a known fixed geometry. However, in practical applications the geometry of the domain is usually not known. This introduces an error, and effectively changes the problem into an anisotropic one. The main result of this paper is a uniqueness result characterizing the isotropic conductivities on convex domains in terms of measurements done on a different domain, which we call the model domain, up to an affine isometry. As data for the inverse problem, we assume the Robin-to-Neumann map and the contact impedance function on the boundary of the model domain to be given. Also, we present a minimization algorithm based on the use of Cotton--York tensor, that finds the pushforward of the isotropic conductivity to our model domain, and also finds the boundary of the original domain up to an affine isometry. This algorithm works also in dimensions higher than three, but then the Cotton--York tensor has to replaced with the Weyl--tensor.

0 citations0 reviewsNo code linkedOpen access