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

BACKGROUND AND OBJECTIVE: 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 a mathematical formulation for finding a current pattern that optimizes an L1-norm fit between a given focal target distribution and volumetric 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. METHODS: We present a linear programming approach that 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 from a pre-filtered set of candidates. RESULTS: The numerical simulation results obtained with both 8- and 20-channel electrode montages suggest that our hypothesis on the benefits of L1-norm data fitting is valid. Compared to an L1-norm regularized L2-norm fitting (L1L2) via semidefinite programming and weighted Tikhonov least-squares method (TLS), the L1L1 results were overall preferable for maximizing the focused current density at the target position, and the ratio between focused and nuisance current magnitudes. CONCLUSIONS: 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 while 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.

Reconstructing subcortical and cortical somatosensory activity via the RAMUS inverse source analysis technique using median nerve SEP data.

This study concerns reconstructing brain activity at various depths based on non-invasive EEG (electroencephalography) scalp measurements. We aimed at demonstrating the potential of the RAMUS (randomized multiresolution scanning) technique in localizing weakly distinguishable far-field sources in combination with coinciding cortical activity. As we have shown earlier theoretically and through simulations, RAMUS is a novel mathematical method that by employing the multigrid concept, allows marginalizing noise and depth bias effects and thus enables the recovery of both cortical and subcortical brain activity. To show this capability with experimental data, we examined the 14-30 ms post-stimulus somatosensory evoked potential (SEP) responses of human median nerve stimulation in three healthy adult subjects. We aim at reconstructing the different response components by evaluating a RAMUS-based estimate for the primary current density in the nervous tissue. We present source reconstructions obtained with RAMUS and compare them with the literature knowledge of the SEP components and the outcome of the unit-noise gain beamformer (UGNB) and standardized low-resolution brain electromagnetic tomography (sLORETA). We also analyzed the effect of the iterative alternating sequential technique, the optimization technique of RAMUS, compared to the classical minimum norm estimation (MNE) technique. Matching with our previous numerical studies, the current results suggest that RAMUS could have the potential to enhance the detection of simultaneous deep and cortical components and the distinction between the evoked sulcal and gyral activity.

Parametrizing the Conditionally Gaussian Prior Model for Source Localization with Reference to the P20/N20 Component of Median Nerve SEP/SEF.

In this article, we focused on developing the conditionally Gaussian hierarchical Bayesian model (CG-HBM), which forms a superclass of several inversion methods for source localization of brain activity using somatosensory evoked potential (SEP) and field (SEF) measurements. The goal of this proof-of-concept study was to improve the applicability of the CG-HBM as a superclass by proposing a robust approach for the parametrization of focal source scenarios. We aimed at a parametrization that is invariant with respect to altering the noise level and the source space size. The posterior difference between the gamma and inverse gamma hyperprior was minimized by optimizing the shape parameter, while a suitable range for the scale parameter can be obtained via the prior-over-measurement signal-to-noise ratio, which we introduce as a new concept in this study. In the source localization experiments, the primary generator of the P20/N20 component was detected in the Brodmann area 3b using the CG-HBM approach and a parameter range derived from the existing knowledge of the Tikhonov-regularized minimum norm estimate, i.e., the classical Gaussian prior model. Moreover, it seems that the detection of deep thalamic activity simultaneously with the P20/N20 component with the gamma hyperprior can be enhanced while using a close-to-optimal shape parameter value.

A mathematical model and iterative inversion for fluorescent optical projection tomography.

Solving the fluorophore distribution in a tomographic setting has been difficult because of the lack of physically meaningful and computationally applicable propagation models. This study concentrates on the direct modelling of fluorescence signals in optical projection tomography (OPT), and on the corresponding inverse problem. The reconstruction problem is solved using emission projections corresponding to a series of rotational imaging positions of the sample. Similarly to the bright field OPT bearing resemblance with the transmission x-ray computed tomography, the fluorescent mode OPT is analogous to x-ray fluorescence tomography (XFCT). As an improved direct model for the fluorescent OPT, we derive a weighted Radon transform based on the XFCT literature. Moreover, we propose a simple and fast iteration scheme for the slice-wise reconstruction of the sample. The developed methods are applied in both numerical experiments and inversion of fluorescent OPT data from a zebrafish embryo. The results demonstrate the importance of propagation modelling and our analysis provides a flexible modelling framework for fluorescent OPT that can easily be modified to adapt to different imaging setups.

A realistic, accurate and fast source modeling approach for the EEG forward problem.

The aim of this paper is to advance electroencephalography (EEG) source analysis using finite element method (FEM) head volume conductor models that go beyond the standard three compartment (skin, skull, brain) approach and take brain tissue inhomogeneity (gray and white matter and cerebrospinal fluid) into account. The new approach should enable accurate EEG forward modeling in the thin human cortical structures and, more specifically, in the especially thin cortices in children brain research or in pathological applications. The source model should thus be focal enough to be usable in the thin cortices, but should on the other side be more realistic than the current standard mathematical point dipole. Furthermore, it should be numerically accurate and computationally fast. We propose to achieve the best balance between these demands with a current preserving (divergence conforming) dipolar source model. We develop and investigate a varying number of current preserving source basis elements n (n=1,…,n=5). For validation, we conducted numerical experiments within a multi-layered spherical domain, where an analytical solution exists. We show that the accuracy increases along with the number of basis elements, while focality decreases. The results suggest that the best balance between accuracy and focality in thin cortices is achieved with n=4 (or in extreme cases even n=3) basis functions, while in thicker cortices n=5 is recommended to obtain the highest accuracy. We also compare the current preserving approach to two further FEM source modeling techniques, namely partial integration and St. Venant, and show that the best current preserving source model outperforms the competing methods with regard to overall balance. For all tested approaches, FEM transfer matrices enable high computational speed. We implemented the new EEG forward modeling approaches into the open source duneuro library for forward modeling in bioelectromagnetism to enable its broader use by the brain research community. This library is build upon the DUNE framework for parallel finite elements simulations and integrates with high-level toolboxes like FieldTrip. Additionally, an inversion test has been implemented using the realistic head model to demonstrate and compare the differences between the aforementioned source models.