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

Sampsa Pursiainen

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

ResearcherUnclaimed source record
astro-ph.IMmath.NANumerical Analysisphysics.med-phastro-ph.EPeess.IVmath-phmath.APmath.MPmath.OC

Catalog footprint

What is connected

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

preprint2022arXiv

Conditionally Exponential Prior in Focal Near- and Far-Field EEG Source Localization via Randomized Multiresolution Scanning (RAMUS)

In this paper, we focus on the inverse problem of reconstructing distributional brain activity with cortical and weakly detectable deep components in non-invasive Electroencephalography. In particular, we aim to generalize the previously extensively used conditionally Gaussian prior (CGP) formalism to achieve distributional reconstructions with higher focality. For this purpose, we introduce as a hierarchical prior, a general exponential distribution, refered to as conditionally exponential prior (CEP). The first-degree CEP corresponds to focality enforcing Laplace prior that suffers from strong depth bias making the deep activity unrecoverable. We sample over multiple resolution levels via RAMUS to reduce this bias as it is known to depend on the resolution of the source space. Moreover, we introduce a procedure based on the physiological a priori knowledge of the brain activity to obtain the shape and scale parameters of the gamma hyperprior that steer the CEP. The posterior estimates are calculated using iterative statistical methods, expectation maximization and iterative alternating sequential algorithm, which we show to be algorithmically similar and to have a close resemblance to the iterative reweighting methods. The performance of CEP is compared with the recent sampling-based dipole localization method Sequential semi-analytic Monte Carlo estimation (SESAME) in numerical experiments of simulated somatosensory evoked potentials related to the human median nerve stimulation. Our results suggest that a hybrid of the first-degree CEP and RAMUS can achieve an accuracy comparable to the second-degree case (CGP) while being more focal. Further, the proposed hybrid is shown to be robust to noise effects and compare well to the dipole reconstructions obtained with SESAME.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Investigation of Wavelength-induced Uncertainties in Full-Wave Radar Tomography of High Contrast Domain: An Application to Small Solar System Bodies

This paper aims to reconstruct the internal structure of a two-dimensional test object via numerically simulated full-wave time domain radar tomography with the presence of wavelength-induced (WI) uncertainties, following from a complex domain structure, and domain diameters 21 or 64 times the wavelength of the signal propagating inside the target. In particular, we consider an application in planetary scientific studies of reconstructing the interior structure of an arbitrary high contrast small Solar System Body (SSSB), i.e., an asteroid, with a probing signal wavelength limited by the instrument and mission payload requirements. Our uncertainty reduction model finds the reconstruction via averaging multiple inverse solutions assuming that the WI deviations in the solutions correspond to random deviations, which we assume to be independent and identically distributed (IID). It incorporates error marginalisation via a randomised signal configuration, spatial-averaging of candidate solutions, frequency-based error marginalisation, and the truncated singular value decomposition (TSVD) filtering technique, based on our assumptions of the phase discrepancy of the signal, domain parameters, and the full-wave forward model. The numerical experiments are performed for 20 and 60 MHz centre frequencies proposed for CubeSat-based radars, the latter being the centre frequency of the Juventas Radar which will be aboard Hera mission to investigate the interior structure of asteroid Dimorphos. A benchmark reconstruction of the target was obtained with the spatial averaging, sparse point density and frequency randomised configuration for both 20 and 60 MHz frequency systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

L1-norm vs. L2-norm fitting in optimizing focal multi-channel tES stimulation: linear and semidefinite programming vs. weighted least squares

This study focuses on Multi-Channel Transcranial Electrical Stimulation, a non-invasive brain method for stimulating neuronal activity under the influence of low-intensity currents. We introduce mathematical formulation for finding a current pattern which optimizes a L1-norm fit between a given focal target distribution and volume current density inside the brain. L1-norm is well-known to favor well-localized or sparse distributions compared to L2-norm (least-squares) fitted estimates. We present a linear programming approach which performs L1-norm fitting and penalization of the current pattern (L1L1) to control the number of non-zero currents. The optimizer filters a large set of candidate solutions using a two-stage metaheuristic search in from a pre-filtered set of candidates. The numerical simulation results, obtained with both a 8- and 20-channel electrode montages, suggest that our hypothesis on the benefits of L1-norm data fitting is valid. As compared to L1-norm regularized L2-norm fitting (L1L2) via semidefinite programming and weighted Tikhonov least-squares method, the L1L1 results were overall preferable with respect to maximizing the focused current density at the target position and the ratio between focused and nuisance current magnitudes. We propose the metaheuristic L1L1 optimization approach as a potential technique to obtain a well-localized stimulus with a controllable magnitude at a given target position. L1L1 finds a current pattern with a steep contrast between the anodal and cathodal electrodes meanwhile suppressing the nuisance currents in the brain, hence, providing a potential alternative to modulate the effects of the stimulation, e.g., the sensation experienced by the subject.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Randomized Multiresolution Scanning in Focal and Fast E/MEG Sensing of Brain Activity with a Variable Depth

We focus on electromagnetoencephalography imaging of the neural activity and, in particular, finding a robust estimate for the primary current distribution via the hierarchical Bayesian model (HBM). Our aim is to develop a reasonably fast maximum a posteriori (MAP) estimation technique which would be applicable for both superficial and deep areas without specific a priori knowledge of the number or location of the activity. To enable source distinguishability for any depth, we introduce a randomized multiresolution scanning (RAMUS) approach in which the MAP estimate of the brain activity is varied during the reconstruction process. RAMUS aims to provide a robust and accurate imaging outcome for the whole brain, while maintaining the computational cost on an appropriate level. The inverse gamma (IG) distribution is applied as the primary hyperprior in order to achieve an optimal performance for the deep part of the brain. In this proof-of-the-concept study, we consider the detection of simultaneous thalamic and somatosensory activity via numerically simulated data modeling the 14-20 ms post-stimulus somatosensory evoked potential and field response to electrical wrist stimulation. Both a spherical and realistic model are utilized to analyze the source reconstruction discrepancies. In the numerically examined case, RAMUS was observed to enhance the visibility of deep components and also marginalizing the random effects of the discretization and optimization without a remarkable computation cost. A robust and accurate MAP estimate for the primary current density was obtained in both superficial and deep parts of the brain.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Tomographic inversion of gravity gradient field for a synthetic Itokawa model

This article investigates reconstructing the internal mass density of a numerical asteroid model using the gradient of a simulated gravity field as synthetic measurement data. Our goal is to advance the mathematical inversion methodology and find feasibility constraints for the resolution, noise and orbit selection for future space missions. We base our model on the shape of the asteroid Itokawa as well as on the recent observations and simulation studies which suggest that the internal density varies, increasing towards the center, and that the asteroid may have a detailed structure. We introduce randomized multiresolution scan algorithm which might provide a robust way to cancel out bias and artifact effects related to the measurement noise and numerical discretization. In this scheme, the inverse algorithm can reconstruct details of various sizes without fixing the exact resolution {\em a priori}, and the randomization minimizes the effect of discretization on the solution. We show that the adopted methodology provides an advantageous way to diminish the surface bias of the inverse solution. The results also suggest that a noise level below 80 Eotvos will be sufficient for the detection of internal voids and high density anomalies, if a sparse set of measurements can be obtained from a close-enough distance to the target.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Advanced boundary electrode modeling for tES and parallel tES/EEG

This paper explores advanced electrode modeling in the context of separate and parallel transcranial electrical stimulation (tES) and electroencephalography (EEG) measurements. We focus on boundary condition based approaches that do not necessitate adding auxiliary elements, e.g. sponges, to the computational domain. In particular, we investigate the complete electrode model (CEM) which incorporates a detailed description of the skin-electrode interface including its contact surface, impedance and normal current distribution. The CEM can be applied for both tES and EEG electrodes which is advantageous when a parallel system is used. In comparison to the CEM, we test two important reduced approaches: the gap model (GAP) and the point electrode model (PEM). We aim to find out the differences of these approaches for a realistic numerical setting based on the stimulation of the auditory cortex. The results obtained suggest, among other things, that GAP and GAP/PEM are sufficiently accurate for the practical application of tES and parallel tES/EEG, respectively. Differences between CEM and GAP were observed mainly in the skin compartment, where only CEM explains the heating effects characteristic to tES.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electroencephalography (EEG) Forward Modeling via H(div) Finite Element Sources with Focal Interpolation

The goal of this study is to develop focal, accurate and robust finite element method (FEM) based approaches which can predict the electric potential on the surface of the computational domain given its structure and internal primary source current distribution. While conducting an EEG evaluation, the placement of source currents to the geometrically complex grey matter compartment is a challenging but necessary task to avoid forward errors attributable to tissue conductivity jumps. Here, this task is approached via a mathematically rigorous formulation, in which the current field is modeled via divergence conforming H(div) basis functions. Both linear and quadratic functions are used while the potential field is discretized via the standard linear Lagrangian (nodal) basis. The resulting model includes dipolar sources which are interpolated into a random set of positions and orientations utilizing two alternative approaches: the position based optimization (PBO) and the mean position/orientation (MPO) method. These results demonstrate that the present dipolar approach can reach or even surpass, at least in some respects, the accuracy of two classical reference methods, the partial integration (PI) and St. Venant (SV) approach which utilize monopolar loads instead of dipolar currents.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electromagnetic 3D subsurface imaging with source sparsity for a synthetic object

This paper concerns electromagnetic 3D subsurface imaging in connection with sparsity of signal sources. We explored an imaging approach that can be implemented in situations that allow obtaining a large amount of data over a surface or a set of orbits but at the same time require sparsity of the signal sources. Characteristic to such a tomography scenario is that it necessitates the inversion technique to be genuinely three-dimensional: For example, slicing is not possible due to the low number of sources. Here, we primarily focused on astrophysical subsurface exploration purposes. As an example target of our numerical experiments we used a synthetic small planetary object containing three inclusions, e.g. voids, of the size of the wavelength. A tetrahedral arrangement of source positions was used, it being the simplest symmetric point configuration in 3D. Our results suggest that somewhat reliable inversion results can be produced within the present a priori assumptions, if the data can be recorded at a specific resolution. This is valuable early-stage knowledge especially for design of future planetary missions in which the payload needs to be minimized, and potentially also for the development of other lightweight subsurface inspection systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Orbiter-to-orbiter tomography: a potential approach for small planetary objects

The goal of this paper is to advance mathematical and computational methodology for orbiter-to-orbiter radio tomography of small planetary objects. In this study, an advanced full waveform forward model is coupled with a total variation based inversion technique. We use a satellite formation model in which a single unit receives a signal that is transmitted by one or more transponder satellites. Numerical results for a 2D domain are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sparse source travel-time tomography of a laboratory target: accuracy and robustness of anomaly detection

This study concerned conebeam travel-time tomography. The focus was on a sparse distribution of signal sources that can be necessary in a challenging in situ environment such as in asteroid tomography. The goal was to approximate the minimum number of source positions needed for robust detection of refractive anomalies, e.g., voids within an asteroid or a casting defects in concrete. Experimental ultrasonic data were recorded utilizing as a target a 150 mm plastic cast cube containing three stones with diameter between 22 and 41 mm. A signal frequency of 55 kHz (35 mm wavelength) was used. Source counts from one to six were tested for different placements. Based on our statistical inversion approach and analysis of the results, three or four sources were found to lead to reliable inversion. The source configurations investigated were also ranked according to their performance. Our results can be used, for example, in the planning of planetary missions as well as in material testing.

0 citations0 reviewsNo code linkedOpen access