Ivan Gligorijevic

Source Separation From Single-Channel Recordings by Combining Empirical-Mode Decomposition and Independent Component Analysis

In biomedical signal processing, it is often the case that many sources are mixed into the measured signal. The goal is usually to analyze one or several of them separately. In the case of multichannel measurements, several blind source separation techniques are available for decomposing the signal into its components [e.g., independent component analysis (ICA)]. However, only a few techniques have been reported for analyses of single-channel recordings. Examples are single-channel ICA (SCICA) and wavelet-ICA (WICA), which all have certain limitations. In this paper, we propose a new method for a single-channel signal decomposition. This method combines empirical-mode decomposition with ICA. We compare the separation performance of our algorithm with SCICA and WICA through simulations, and we show that our method outperforms the other two, especially for high noise-to-signal ratios. The performance of the new algorithm was also demonstrated in two real-life applications.

A Multichannel Integrated Circuit for Electrical Recording of Neural Activity, With Independent Channel Programmability

Since a few decades, micro-fabricated neural probes are being used, together with microelectronic interfaces, to get more insight in the activity of neuronal networks. The need for higher temporal and spatial recording resolutions imposes new challenges on the design of integrated neural interfaces with respect to power consumption, data handling and versatility. In this paper, we present an integrated acquisition system for in vitro and in vivo recording of neural activity. The ASIC consists of 16 low-noise, fully-differential input channels with independent programmability of its amplification (from 100 to 6000 V/V) and filtering (1-6000 Hz range) capabilities. Each channel is AC-coupled and implements a fourth-order band-pass filter in order to steeply attenuate out-of-band noise and DC input offsets. The system achieves an input-referred noise density of 37 nV/√Hz, a NEF of 5.1, a CMRR >; 60 dB, a THD <; 1% and a sampling rate of 30 kS/s per channel, while consuming a maximum of 70 μA per channel from a single 3.3 V. The ASIC was implemented in a 0.35 μm CMOS technology and has a total area of 5.6 × 4.5 mm2. The recording system was successfully validated in in vitro and in vivo experiments, achieving simultaneous multichannel recordings of cell activity with satisfactory signal-to-noise ratios.

Multi-structural Signal Recovery for Biomedical Compressive Sensing

Compressive sensing has shown significant promise in biomedical fields. It reconstructs a signal from sub-Nyquist random linear measurements. Classical methods only exploit the sparsity in one domain. A lot of biomedical signals have additional structures, such as multi-sparsity in different domains, piecewise smoothness, low rank, etc. We propose a framework to exploit all the available structure information. A new convex programming problem is generated with multiple convex structure-inducing constraints and the linear measurement fitting constraint. With additional a priori information for solving the underdetermined system, the signal recovery performance can be improved. In numerical experiments, we compare the proposed method with classical methods. Both simulated data and real-life biomedical data are used. Results show that the newly proposed method achieves better reconstruction accuracy performance in term of both L1 and L2 errors.

Time-frequency heart rate variability characteristics of young adults during physical, mental and combined stress in laboratory environment

The goal of this study was to evaluate the changes in heart rate variability (HRV) parameters due to a specific physical, mental or combined load. More specifically, the difference in effect between mental load and physical activity is studied. In addition, the effect of the combined physical and mental demand on the HRV parameters was examined and compared with the changes during the single task. In a laboratory environment, 28 subjects went through a protocol with different types of load (physical and/or mental), each followed by a period of rest. Continuous wavelet transformation was applied to create time series of instantaneous power and frequency in specified frequency bands (LF and HF). HF could distinguish the active conditions from the rest condition, meaning that HRV is sensitive to any change in mental or physical state. Differences in HRV parameters were observed between physical, mental and the combined load. In conclusion, we were able to distinguish between rest, physical and mental condition by combining different HRV characteristics. The addition of a mental load to a physical task had an extra effect on the HRV characteristics.

A new and fast approach towards sEMG decomposition

The decomposition of high-density surface EMG (HD-sEMG) interference patterns into the contribution of motor units is still a challenging task. We introduce a new, fast solution to this problem. The method uses a data-driven approach for selecting a set of electrodes to enable discrimination of present motor unit action potentials (MUAPs). Then, using shapes detected on these channels, the hierarchical clustering algorithm as reported by Quian Quiroga et al. (Neural Comput 16:1661–1687, 2004) is extended for multichannel data in order to obtain the motor unit action potential (MUAP) signatures. After this first step, more motor unit firings are obtained using the extracted signatures by a novel demixing technique. In this demixing stage, we propose a time-efficient solution for the general convolutive system that models the motor unit firings on the HD-sEMG grid. We constrain this system by using the extracted signatures as prior knowledge and reconstruct the firing patterns in a computationally efficient way. The algorithm performance is successfully verified on simulated data containing up to 20 different MUAP signatures. Moreover, we tested the method on real low contraction recordings from the lateral vastus leg muscle by comparing the algorithm’s output to the results obtained by manual analysis of the data from two independent trained operators. The proposed method showed to perform about equally successful as the operators.

Combining EMD with ICA for extracting independent sources from single channel and two-channel data

Blind Source Separation (BSS) techniques are frequently needed in the processing of biomedical signals. This need comes from the fact that these signals are often composed of many different sources, which are mixed in the measured signal. However, we are usually only interested in examining one or a limited set of sources of interest separately. A variety of algorithms exist for separating multichannel mixtures into its independent sources (e.g. different Independent Component Analysis (ICA) techniques). These techniques only work if the number of channels is larger than, or equal to the number of sources present in the signal. On the other hand, only a few algorithms have been reported for the analysis of single channel sources, or other mixtures where the number of sources is higher than the number of channels. In this work we show a new technique which combines Empirical Mode Decomposition (EMD) and Independent Component Analysis (ICA). We will show that this technique is capable in separating independent sources when the number of these sources is higher than the number of channels available. We show the performance in single channel and two-channel biosignal processing.

Multi-sparse signal recovery for compressive sensing

Signal recovery is one of the key techniques of compressive sensing (CS). It reconstructs the original signal from the linear sub-Nyquist measurements. Classical methods exploit the sparsity in one domain to formulate the L0 norm optimization. Recent investigation shows that some signals are sparse in multiple domains. To further improve the signal reconstruction performance, we can exploit this multi-sparsity to generate a new convex programming model. The latter is formulated with multiple sparsity constraints in multiple domains and the linear measurement fitting constraint. It improves signal recovery performance by additional a priori information. Since some EMG signals exhibit sparsity both in time and frequency domains, we take them as example in numerical experiments. Results show that the newly proposed method achieves better performance for multi-sparse signals.

Neuronal activity in the bed nucleus of the stria terminalis in a rat model for obsessive–compulsive disorder

In search of a new potential target for deep brain stimulation in patients with obsessive-compulsive disorder (OCD), we evaluated the single-cell activity of neurons in the bed nucleus of the stria terminalis (BST) in urethane-anesthetized rats in an animal model for OCD, the schedule-induced polydipsia (SIP) model, and compared this to the BST activity in control rats and to a third group of rats which were introduced in the model but did not develop the SIP, and thus were considered resistant. We compared the firing rate and firing pattern of BST neurons between these groups, between hemispheres and made a correlation of the firing rate and firing pattern to the position in the BST. The variability of BST neurons in SIP rats was lower and the randomness higher than BST neurons in control rats or resistant rats. The firing rate of BST neurons in SIP rats was significantly higher and the burst index lower than BST neurons in resistant rats but not in control rats. Also, neurons from the right hemisphere in the SIP group had a higher burst index than neurons from the left hemisphere. However, this is opposite in the resistant and control group. Third, we found a higher bursting index with increasing (more ventral) depth of recording. These findings suggest that schedule-induced polydipsia, which models compulsive behavior in humans, induces a change in firing behavior of BST neurons.

Automated way to obtain motor units' signatures and estimate their firing patterns during voluntary contractions using HD-sEMG

A new, automated way to obtain signatures of active motor units (MUs) from high density surface EMG recordings during voluntary contractions is presented. It relies on clustering of repetitive shapes corresponding to different MU action potentials (MUAPs) present. The number of clusters and the mean shapes of the MUAPs as observed on the electrode grid, are estimated in a fast way without user interaction. The algorithm is tested on simulated signals mimicking a small muscle. Our results show that at least 8 MUAPs can be reliably reconstructed and their MU mean firing frequencies can be estimated.

Towards continuous and real-time attention monitoring at work: reaction time versus brain response

Abstract Continuous and objective measurement of the user attention state still represents a major challenge in the ergonomics research. Recently available wearable electroencephalography (EEG) opens new opportunities for objective and continuous evaluation of operators’ attention, which may provide a new paradigm in ergonomics. In this study, wearable EEG was recorded during simulated assembly operation, with the aim to analyse P300 event-related potential component, which provides reliable information on attention processing. In parallel, reaction times (RTs) were recorded and the correlation between these two attention-related modalities was investigated. Negative correlation between P300 amplitudes and RTs has been observed on the group level (p < .001). However, on the individual level, the obtained correlations were not consistent. As a result, we propose the P300 amplitude for accurate attention monitoring in ergonomics research. On the other hand, no significant correlation between RTs and P300 latency was found on group, neither on individual level. Practitioner Summary: Ergonomic studies of assembly operations mainly investigated physical aspects, while mental states of the assemblers were not sufficiently addressed. Presented study aims at attention tracking, using realistic workplace replica. It is shown that drops in attention could be successfully traced only by direct brainwave observation, using wireless electroencephalographic measurements.