NeuroMotion: Open-source platform with neuromechanical and deep network modules to generate surface EMG signals during voluntary movement.

Neuromechanical studies investigate how the nervous system interacts with the musculoskeletal (MSK) system to generate volitional movements. Such studies have been supported by simulation models that provide insights into variables that cannot be measured experimentally and allow a large number of conditions to be tested before the experimental analysis. However, current simulation models of electromyography (EMG), a core physiological signal in neuromechanical analyses, remain either limited in accuracy and conditions or are computationally heavy to apply. Here, we provide a computational platform to enable future work to overcome these limitations by presenting NeuroMotion, an open-source simulator that can modularly test a variety of approaches to the full-spectrum synthesis of EMG signals during voluntary movements. We demonstrate NeuroMotion using three sample modules. The first module is an upper-limb MSK model with OpenSim API to estimate the muscle fibre lengths and muscle activations during movements. The second module is BioMime, a deep neural network-based EMG generator that receives nonstationary physiological parameter inputs, like the afore-estimated muscle fibre lengths, and efficiently outputs motor unit action potentials (MUAPs). The third module is a motor unit pool model that transforms the muscle activations into discharge timings of motor units. The discharge timings are convolved with the output of BioMime to simulate EMG signals during the movement. We first show how MUAP waveforms change during different levels of physiological parameter variations and different movements. We then show that the synthetic EMG signals during two-degree-of-freedom hand and wrist movements can be used to augment experimental data for regressing joint angles. Ridge regressors trained on the synthetic dataset were directly used to predict joint angles from experimental data. In this way, NeuroMotion was able to generate full-spectrum EMG for the first use-case of human forearm electrophysiology during voluntary hand, wrist, and forearm movements. All intermediate variables are available, which allows the user to study cause-effect relationships in the complex neuromechanical system, fast iterate algorithms before collecting experimental data, and validate algorithms that estimate non-measurable parameters in experiments. We expect this modular platform will enable validation of generative EMG models, complement experimental approaches and empower neuromechanical research.

A myoelectric digital twin for fast and realistic modelling in deep learning.

Muscle electrophysiology has emerged as a powerful tool to drive human machine interfaces, with many new recent applications outside the traditional clinical domains, such as robotics and virtual reality. However, more sophisticated, functional, and robust decoding algorithms are required to meet the fine control requirements of these applications. Deep learning has shown high potential in meeting these demands, but requires a large amount of high-quality annotated data, which is expensive and time-consuming to acquire. Data augmentation using simulations, a strategy applied in other deep learning applications, has never been attempted in electromyography due to the absence of computationally efficient models. We introduce a concept of Myoelectric Digital Twin - highly realistic and fast computational model tailored for the training of deep learning algorithms. It enables simulation of arbitrary large and perfectly annotated datasets of realistic electromyography signals, allowing new approaches to muscular signal decoding, accelerating the development of human-machine interfaces.

Fast Approximation of EEG Forward Problem and Application to Tissue Conductivity Estimation.

Bioelectric source analysis in the human brain from scalp electroencephalography (EEG) signals is sensitive to the conductivities of different head tissues. The conductivity of tissues is subject dependent, so non-invasive methods for conductivity estimation are necessary to fine tune EEG models. To do so, the EEG forward problem solution (so-called lead field matrix) must be computed for a large number of conductivity configurations. Computing a lead field requires a matrix inversion which is computationally intensive for realistic head models. Thus, the required time for computing a large number of lead fields can become impractical. In this work, we propose to approximate the lead field matrix for a set of conductivity configurations, using the exact solution only for a small set of support points in the conductivity space. Our approach accelerates the computation time, while controlling the approximation error. Our method is tested on simulated and measured EEG data for brain and skull conductivity estimation. This test demonstrates that the approximation does not introduce any bias and runs significantly faster than if exact lead field were to be computed.

Strategies for statistical thresholding of source localization maps in magnetoencephalography and estimating source extent.

BACKGROUND: Magnetoencephalography allows defining non-invasively the spatio-temporal activation of brain networks thanks to source localization algorithms. A major difficulty of MNE and beamforming methods, two classically used techniques, is the definition of proper thresholds that allow deciding the extent of activated cortex. NEW METHOD: We investigated two strategies for computing a threshold, taking into account the difficult multiple comparison issue. The strategies were based either on parametric statistics (Bonferroni, FDR correction) or on empirical estimates (local FDR and a custom measure based on the survival function). RESULTS: We found thanks to the simulations that parametric methods based on the sole estimation of H0 (Bonferroni, FDR) performed poorly, in particular in high SNR situations. This is due to the spatial leakage originating from the source localization methods, which give a 'blurred' reconstruction of the patch extension: the higher the SNR, the more this effect is visible. COMPARISON WITH EXISTING METHODS: Adaptive methods such as local FDR or our proposed 'concavity threshold' performed better than Bonferroni or classical FDR. We present an application to real data originating from auditory stimulation in MEG. CONCLUSION: In order to estimate source extent, adaptive strategies should be preferred to parametric statistics when dealing with 'leaking' source reconstruction algorithms.