A tutorial review of electrical neuroimaging from group-average to single-trial event-related potentials.

This tutorial review details some of the recent advances in signal analyses applied to event-related potential (ERP) data. These "electrical neuroimaging" analyses provide reference-independent measurements of response strength and response topography that circumvent statistical and interpretational caveats of canonical ERP analysis methods while also taking advantage of the greater information provided by high-density electrode montages. Electrical neuroimaging can be applied across scales ranging from group-averaged ERPs to single-subject and single-trial datasets. We illustrate these methods with a tutorial dataset and place particular emphasis on their suitability for studies of clinical and/or developmental populations.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Assessment of quadriceps muscle inactivation with a new electrical stimulation paradigm.

INTRODUCTION: In this study we evaluated the validity of garment-based quadriceps stimulation (GQS) for assessment of muscle inactivation in comparison with femoral nerve stimulation (FNS). METHODS: Inactivation estimates (superimposed doublet torque), self-reported discomfort, and twitch and doublet contractile properties were compared between GQS and FNS in 15 healthy subjects. RESULTS: Superimposed doublet torque was significantly lower for GQS than for FNS at 20% and 40% maximum voluntary contraction (MVC) (P < 0.01), but not at 60%, 80%, and 100% MVC. Discomfort scores were systematically lower for GQS than for FNS (P < 0.05). Resting twitch and doublet peak torque were lower for GQS, and time to peak torque was shorter for GQS than for FNS (P < 0.01). CONCLUSIONS: GQS can be used with confidence for straightforward evaluation of quadriceps muscle inactivation, whereas its validity for assessment of contractile properties remains to be determined. Muscle Nerve 51: 117-124, 2015.

CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

Continuous assessment of electrical epileptic activity in acute stroke.

OBJECTIVE: To determine the incidence and risk factors of electrical seizures and other electrical epileptic activity using continuous EEG (cEEG) in patients with acute stroke. METHODS: One hundred consecutive patients with acute stroke admitted to our stroke unit underwent cEEG using 10 electrodes. In addition to electrical seizures, repetitive focal sharp waves (RSHWs), repetitive focal spikes (RSPs), and periodic lateralized epileptic discharges (PLEDs) were recorded. RESULTS: In the 100 patients, cEEG was recorded for a mean duration of 17 hours 34 minutes (range 1 hour 12 minutes to 37 hours 10 minutes). Epileptic activity occurred in 17 patients and consisted of RSHWs in seven, RSPs in seven, and PLEDs in three. Electrical seizures occurred in two patients. On univariate Cox regression analysis, predictors for electrical epileptic activity were stroke severity (high score on the National Institutes of Health Stroke Scale) (hazard ratio [HR] 1.12; p = 0.002), cortical involvement (HR 5.71; p = 0.021), and thrombolysis (HR 3.27; p = 0.040). Age, sex, stroke type, use of EEG-modifying medication, and cardiovascular risk factors were not predictors of electrical epileptic activity. On multivariate analysis, stroke severity was the only independent predictor (HR 1.09; p = 0.016). CONCLUSION: In patients with acute stroke, electrical epileptic activity occurs more frequently than previously suspected.

Use of anisotropic modelling in electrical impedance tomography: description of method and preliminary assessment of utility in imaging brain function in the adult human head.

Electrical Impedance Tomography (EIT) is an imaging method which enables a volume conductivity map of a subject to be produced from multiple impedance measurements. It has the potential to become a portable non-invasive imaging technique of particular use in imaging brain function. Accurate numerical forward models may be used to improve image reconstruction but, until now, have employed an assumption of isotropic tissue conductivity. This may be expected to introduce inaccuracy, as body tissues, especially those such as white matter and the skull in head imaging, are highly anisotropic. The purpose of this study was, for the first time, to develop a method for incorporating anisotropy in a forward numerical model for EIT of the head and assess the resulting improvement in image quality in the case of linear reconstruction of one example of the human head. A realistic Finite Element Model (FEM) of an adult human head with segments for the scalp, skull, CSF, and brain was produced from a structural MRI. Anisotropy of the brain was estimated from a diffusion tensor-MRI of the same subject and anisotropy of the skull was approximated from the structural information. A method for incorporation of anisotropy in the forward model and its use in image reconstruction was produced. The improvement in reconstructed image quality was assessed in computer simulation by producing forward data, and then linear reconstruction using a sensitivity matrix approach. The mean boundary data difference between anisotropic and isotropic forward models for a reference conductivity was 50%. Use of the correct anisotropic FEM in image reconstruction, as opposed to an isotropic one, corrected an error of 24 mm in imaging a 10% conductivity decrease located in the hippocampus, improved localisation for conductivity changes deep in the brain and due to epilepsy by 4-17 mm, and, overall, led to a substantial improvement on image quality. This suggests that incorporation of anisotropy in numerical models used for image reconstruction is likely to improve EIT image quality.

Benchmarking project of the Department of Electrical Engineering of Lappeenranta University of Technology

Lappeenrannan teknillisen korkeakoulun sähkötekniikan osasto muutti 1.8.2005 sähkötekniikan tutkinnon kaksiportaiseksi ja vastaamaan näin Bologna-prosessia ja Suomen yliopistolainsäädäntöä. Tutkinnonuudistuksen myötä osasto haluaa varmistaa sähkötekniikan tutkintojen laadun ja vertailtavuuden sekä parantaa opiskelijoidensekä henkilökunnan liikkuvuutta. Tutkintojen laatu ja vertailtavuus osoitetaan sähkötekniikan osaston benchmark-projektilla, jossa kerätään tietoja maisteri- ja tohtorintutkintoa tarjoavista eurooppalaisista yliopistoista. Diplomityö käsittää BM-projektin kolmannen vaiheen suunnittelun ja toteutuksen sisältäen teoriaa benchmark-projekteille tyypillisistä toimintatavoista. Hyväksi havaitut menetelmiä, kuten kyselyt ja matriisit, on tässä työssä otettu soveltuvin osin sähkötekniikan osaston BM-projektin työkaluiksi. Diplomityössä analysoidaan työkalujen avulla BM-kumppaneilta kerättyjä tietoja sekä esitetään ratkaisuja, miten sähkötekniikan osastolla voidaan jatkaa parhaiden toimintatapojen löytämistä.

Electrical source imaging for presurgical focus localization in epilepsy patients with normal MRI.

PURPOSE: Patients with magnetic resonance (MR)-negative focal epilepsy (MRN-E) have less favorable surgical outcomes (between 40% and 70%) compared to those in whom an MRI lesion guides the site of surgical intervention (60-90%). Patients with extratemporal MRN-E have the worst outcome (around 50% chance of seizure freedom). We studied whether electroencephalography (EEG) source imaging (ESI) of interictal epileptic activity can contribute to the identification of the epileptic focus in patients with normal MRI. METHODS: We carried out ESI in 10 operated patients with nonlesional MRI and a postsurgical follow-up of at least 1 year. Five of the 10 patients had extratemporal lobe epilepsy. Evaluation comprised surface and intracranial EEG monitoring of ictal and interictal events, structural MRI, [(18)F]fluorodeoxyglucose positron emission tomography (FDG-PET), ictal and interictal perfusion single photon emission computed tomography (SPECT) scans. Eight of the 10 patients also underwent intracranial monitoring. RESULTS: ESI correctly localized the epileptic focus within the resection margins in 8 of 10 patients, 9 of whom experienced favorable postsurgical outcomes. DISCUSSION: The results highlight the diagnostic value of ESI and encourage broadening its application to patients with MRN-E. If the surface EEG contains fairly localized spikes, ESI contributes to the presurgical decision process.

Spatial distribution of motor units recruited during electrical stimulation of the quadriceps muscle versus the femoral nerve.

INTRODUCTION: In this study we investigated differences in the spatial recruitment of motor units (MUs) in the quadriceps when electrical stimulation is applied over the quadriceps belly versus the femoral nerve. METHODS: M-waves and mechanical twitches were evoked using over-the-quadriceps and femoral nerve stimulation of gradually increasing intensity from 22 young, healthy subjects. Spatial recruitment was investigated using recruitment curves of M-waves recorded from the vastus medialis (VM) and vastus lateralis (VL) and of twitches recorded from the quadriceps. RESULTS: At maximal stimulation intensity (Imax), no differences were found between nerve and over-the-quadriceps stimulation. At submaximal intensities, VL M-wave amplitude was higher for over-the-quadriceps stimulation at 40% Imax, and peak twitch force was greater for nerve stimulation at 60% and 80% Imax. CONCLUSIONS: For the VM, MU spatial recruitment during nerve and over-the-quadriceps stimulation of increasing intensity occurred in a similar manner, whereas significant differences were observed for the VL. Muscle Nerve, 2013.

Single acquisition electrical property mapping based on relative coil sensitivities: A proof-of-concept demonstration.

PURPOSE: All methods presented to date to map both conductivity and permittivity rely on multiple acquisitions to compute quantitatively the magnitude of radiofrequency transmit fields, B1+. In this work, we propose a method to compute both conductivity and permittivity based solely on relative receive coil sensitivities ( B1-) that can be obtained in one single measurement without the need to neither explicitly perform transmit/receive phase separation nor make assumptions regarding those phases. THEORY AND METHODS: To demonstrate the validity and the noise sensitivity of our method we used electromagnetic finite differences simulations of a 16-channel transceiver array. To experimentally validate our methodology at 7 Tesla, multi compartment phantom data was acquired using a standard 32-channel receive coil system and two-dimensional (2D) and 3D gradient echo acquisition. The reconstructed electric properties were correlated to those measured using dielectric probes. RESULTS: The method was demonstrated both in simulations and in phantom data with correlations to both the modeled and bench measurements being close to identity. The noise properties were modeled and understood. CONCLUSION: The proposed methodology allows to quantitatively determine the electrical properties of a sample using any MR contrast, with the only constraint being the need to have 4 or more receive coils and high SNR. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Comparison of electrical nerve stimulation, electrical muscle stimulation and magnetic nerve stimulation to assess the neuromuscular function of the plantar flexor muscles.

INTRODUCTION: As it might lead to less discomfort, magnetic nerve stimulation (MNS) is increasingly used as an alternative to electrical stimulation methods. Yet, MNS and electrical nerve stimulation (ENS) and electrical muscle stimulation (EMS) have not been formally compared for the evaluation of plantar flexor neuromuscular function. METHODS: We quantified plantar flexor neuromuscular function with ENS, EMS and MNS in 10 volunteers in fresh and fatigued muscles. Central alterations were assessed through changes in voluntary activation level (VAL) and peripheral function through changes in M-wave, twitch and doublet (PS100) amplitudes. Discomfort associated with 100-Hz paired stimuli delivered with each method was evaluated on a 10-cm visual analog scale. RESULTS: VAL, agonist and antagonist M-wave amplitudes and PS100 were similar between the different methods in both fresh and fatigued states. Potentiated peak twitch was lower in EMS compared to ENS, whereas no difference was found between ENS and MNS for any parameter. Discomfort associated with MNS (1.5 ± 1.4 cm) was significantly less compared to ENS (5.5 ± 1.9 cm) and EMS (4.2 ± 2.6 cm) (p < 0.05). CONCLUSION: When PS100 is used to evaluate neuromuscular properties, MNS, EMS and ENS can be used interchangeably for plantar flexor neuromuscular function assessment as they provide similar evaluation of central and peripheral factors in unfatigued and fatigued states. Importantly, electrical current spread to antagonist muscles was similar between the three methods while discomfort from MNS was much less compared to ENS and EMS. MNS may be potentially employed to assess neuromuscular function of plantar flexor muscles in fragile populations.

Extra Forces induced by wide-pulse, high-frequency electrical stimulation: Occurrence, magnitude, variability and underlying mechanisms.

OBJECTIVE: In contrast to conventional (CONV) neuromuscular electrical stimulation (NMES), the use of "wide-pulse, high-frequencies" (WPHF) can generate higher forces than expected by the direct activation of motor axons alone. We aimed at investigating the occurrence, magnitude, variability and underlying neuromuscular mechanisms of these "Extra Forces" (EF). METHODS: Electrically-evoked isometric plantar flexion force was recorded in 42 healthy subjects. Additionally, twitch potentiation, H-reflex and M-wave responses were assessed in 13 participants. CONV (25Hz, 0.05ms) and WPHF (100Hz, 1ms) NMES consisted of five stimulation trains (20s on-90s off). RESULTS: K-means clustering analysis disclosed a responder rate of almost 60%. Within this group of responders, force significantly increased from 4% to 16% of the maximal voluntary contraction force and H-reflexes were depressed after WPHF NMES. In contrast, non-responders showed neither EF nor H-reflex depression. Twitch potentiation and resting EMG data were similar between groups. Interestingly, a large inter- and intrasubject variability of EF was observed. CONCLUSION: The responder percentage was overestimated in previous studies. SIGNIFICANCE: This study proposes a novel methodological framework for unraveling the neurophysiological mechanisms involved in EF and provides further evidence for a central contribution to EF in responders.

Two-dimensional probabilistic inversion of plane-wave electromagnetic data: Methodology, model constraints and joint inversion with electrical resistivity data

Probabilistic inversion methods based on Markov chain Monte Carlo (MCMC) simulation are well suited to quantify parameter and model uncertainty of nonlinear inverse problems. Yet, application of such methods to CPU-intensive forward models can be a daunting task, particularly if the parameter space is high dimensional. Here, we present a 2-D pixel-based MCMC inversion of plane-wave electromagnetic (EM) data. Using synthetic data, we investigate how model parameter uncertainty depends on model structure constraints using different norms of the likelihood function and the model constraints, and study the added benefits of joint inversion of EM and electrical resistivity tomography (ERT) data. Our results demonstrate that model structure constraints are necessary to stabilize the MCMC inversion results of a highly discretized model. These constraints decrease model parameter uncertainty and facilitate model interpretation. A drawback is that these constraints may lead to posterior distributions that do not fully include the true underlying model, because some of its features exhibit a low sensitivity to the EM data, and hence are difficult to resolve. This problem can be partly mitigated if the plane-wave EM data is augmented with ERT observations. The hierarchical Bayesian inverse formulation introduced and used herein is able to successfully recover the probabilistic properties of the measurement data errors and a model regularization weight. Application of the proposed inversion methodology to field data from an aquifer demonstrates that the posterior mean model realization is very similar to that derived from a deterministic inversion with similar model constraints.

Methods of Design and Testing of Real Time Systems

Teollusuussovelluksissa vaaditaan nykyisin yhä useammin reaaliaikaista tiedon käsittelyä. Luotettavuus on yksi tärkeimmistä reaaliaikaiseen tiedonkäsittelyyn kykenevän järjestelmän ominaisuuksista. Sen saavuttamiseksi on sekä laitteisto, että ohjelmisto testattava. Tämän työn päätavoitteena on laitteiston testaaminen ja laitteiston testattavuus, koska luotettava laitteistoalusta on perusta tulevaisuuden reaaliaikajärjestelmille. Diplomityössä esitetään digitaaliseen signaalinkäsittelyyn soveltuvan prosessorikortin suunnittelu. Prosessorikortti on tarkoitettu sähkökoneiden ennakoivaa kunnonvalvontaa varten. Uusimmat DFT (Desing for Testability) menetelmät esitellään ja niitä sovelletaan prosessorikortin sunnittelussa yhdessä vanhempien menetelmien kanssa. Kokemukset ja huomiot menetelmien soveltuvuudesta raportoidaan työn lopussa. Työn tavoitteena on kehittää osakomponentti web -pohjaiseen valvontajärjestelmään, jota on kehitetty Sähkötekniikan osastolla Lappeenrannan teknillisellä korkeakoululla.

Architectural characterization of a delta-front reservoir analogue combining ground penetrating radar and electrical resistivity tomography: Roda Sandstone (Lower Eocene, Graus-Tremp basin, Spain)

Three-dimensional reconstruction of reservoir analogues can be improved combining data from different geophysical methods. Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT) data are valuable tools, since they provide subsurface information from internal architecture and facies distribution of sedimentary rock bodies, enabling the upgrading of depositional models and heterogeneity reconstruction. The Lower Eocene Roda Sandstone is a well-known deltaic complex widely studied as a reservoir analogue that displays a series of sandstone wedges with a general NE to SW progradational trend. To provide a better understanding of internal heterogeneity of a 10m-thick progradational delta-front sandstone unit, 3D GPR data were acquired. In addition, common midpoints (CMP) to measure the sandstone subsoil velocity, test profiles with different frequency antennas (25, 50 and 100MHz) and topographic data for subsequent correction in the geophysical data were also obtained. Three ERT profiles were also acquired to further constrain GPR analysis. These geophysical results illustrate the geometry of reservoir analogue heterogeneities both depositional and diagenetic in nature, improving and complementing previous outcrop-derived data. GPR interpretation using radar stratigraphy principles and attributes analysis provided: 1)tridimensional geometry of major stratigraphic surfaces that define four units in the GPR Prism, 2) image the internal architecture of the units and their statistical study of azimuth and dips, useful for a quick determination of paleocurrent directions. These results were used to define the depositional architecture of the progradational sandbody that shows an arrangement in very-high-frequency sequences characterized by clockwise paleocurrent variations and decrease of the sedimentary flow, similar to those observed at a greater scale in the same system. This high-frequency sequential arrangement has been attributed to the autocyclic dynamics of a supply-dominated delta- front where fluvial and tidal currents are in competition. The resistivity models enhanced the viewing of reservoir quality associated with cement distribution caused by depositional and early diagenetic processes related to the development of transgressive and regressive systems tracts in igh-frequency sequences.