Investigation of nonlinear pupil dynamics by recurrence quantification analysis.

Pupil is controlled by the autonomous nervous system (ANS). It shows complex movements and changes of size even in conditions of constant stimulation. The possibility of extracting information on ANS by processing data recorded during a short experiment using a low cost system for pupil investigation is studied. Moreover, the significance of nonlinear information contained in the pupillogram is investigated. We examined 13 healthy subjects in different stationary conditions, considering habitual dental occlusion (HDO) as a weak stimulation of the ANS with respect to the maintenance of the rest position (RP) of the jaw. Images of pupil captured by infrared cameras were processed to estimate position and size on each frame. From such time series, we extracted linear indexes (e.g., average size, average displacement, and spectral parameters) and nonlinear information using recurrence quantification analysis (RQA). Data were classified using multilayer perceptrons and support vector machines trained using different sets of input indexes: the best performance in classification was obtained including nonlinear indexes in the input features. These results indicate that RQA nonlinear indexes provide additional information on pupil dynamics with respect to linear descriptors, allowing the discrimination of even a slight stimulation of the ANS. Their use in the investigation of pathology is suggested.

Effects of transducer size on impedance spectroscopy measurements.

The response to an electric field of electrolytic solutions, gels, liquid crystals, and other soft materials is described by the drift-diffusion and Poisson equations. Existing models, used for the interpretation of experimental data, usually consider the system as one dimensional (1D), which is valid only for an infinite electrode size. Here we solve numerically the model equations in 2D, considering a circular electrode with a finite radius, and discuss the limit of validity of the 1D approximation.

Investigation of motor unit recruitment during stimulated contractions of tibialis anterior muscle.

This work investigated motor unit (MU) recruitment during transcutaneous electrical stimulation (TES) of the tibialis anterior (TA) muscle, using experimental and simulated data. Surface electromyogram (EMG) and torque were measured during electrically-elicited contractions at different current intensities, on eight healthy subjects. EMG detected during stimulation (M-wave) was simulated selecting the elicited MUs on the basis of: (a) the simulated current density distribution in the territory of each MU and (b) the excitation threshold characteristic of the MU. Exerted force was simulated by adding the contribution of each of the elicited MUs. The effects of different fat layer thickness (between 2 and 8mm), different distributions of excitation thresholds (random excitation threshold, higher threshold for larger MUs or smaller MUs), and different MU distributions within the muscle (random distribution, larger MU deeper in the muscle, smaller MU deeper) on EMG variables and torque were tested. Increase of the current intensity led to a first rapid increase of experimental M-wave amplitude, followed by a plateau. Further increases of the stimulation current determined an increase of the exerted force, without relevant changes of the M-wave. Similar results were obtained in simulations. Rate of change of conduction velocity (CV) and leading coefficient of the second order polynomial interpolating the force vs. stimulation level curve were estimated as a function of increasing current amplitudes. Experimental data showed an increase of estimated CV with increasing levels of the stimulation current (for all subjects) and a positive leading coefficient of force vs. stimulation current curve (for five of eight subjects). Simulations matched the experimental results only when larger MUs were preferably located deeper in the TA muscle (in line with a histochemical study). Marginal effect of MU excitation thresholds was observed, suggesting that MUs closer to the stimulation electrode are recruited first during TES regardless of their excitability.

Surface EMG: the issue of electrode location.

This paper contributes to clarifying the conditions under which electrode position for surface EMG detection is critical and leads to estimates of EMG variables that are different from those obtained in other nearby locations. Whereas a number of previous works outline the need to avoid the innervation zone (or the muscle belly), many authors place electrodes in the central part or bulge of the muscle of interest where the innervation zone is likely to be. Computer simulations are presented to explain the effect of the innervation zone on amplitude, frequency and conduction velocity estimates from the signal and the need to avoid placing electrodes near it. Experimental signals recorded from some superficial muscles of the limbs and trunk (abductor pollicis brevis, flexor pollicis brevis, biceps, upper trapezius, vastus medialis, vastus lateralis) were processed providing support for the findings obtained from simulations. The use of multichannel techniques is recommended to estimate the location of the innervation zone and to properly choose the optimal position of the detection point(s) allowing meaningful estimates of EMG variables during movement analysis.

Sensitivity of surface EMG-based conduction velocity estimates to local tissue in-homogeneities--influence of the number of channels and inter-channel distance.

The aim of this simulation study was to investigate the influence of local tissue in-homogeneities on the estimates of muscle fiber conduction velocity (CV) from surface EMG signals. A recently developed analytical surface EMG model was used to generate simulated surface EMG signals from a planar layered volume conductor, comprised of the muscle tissue and fat layer, with spheres (1 mm radius) in the fat layer of conductivity different from the surrounding tissue. CV was estimated with a maximum likelihood multi-channel approach, varying the number of channels and the inter-channel distance used for the estimate. The action potentials detected along the muscle fiber direction changed shape due to the presence of the in-homogeneities, thus affecting CV estimates. CV estimates were influenced by the location of the in-homogeneities with respect to the fiber and detection electrodes. The maximum percent variation of CV estimates due to the presence of in-homogeneities decreased with increasing number of channels and inter-channel distance: 19.6% (2 channels), 12.1% (3 channels), 6.4% (4 channels), for 5 mm inter-channel distance, and 12.0% (2 channels), 5.2% (3 channels), 2.4% (4 channels), for 10 mm inter-channel distance (for double differential detection). The results were in agreement and explained previous experimental findings. It was concluded that multi-channel methods for CV estimation significantly reduce the sensitivity of CV estimates to tissue in-homogeneities.

Advances in surface electromyographic signal simulation with analytical and numerical descriptions of the volume conductor.

Surface electromyographic (EMG) signal modelling is important for signal interpretation, testing of processing algorithms, detection system design and didactic purposes. Various surface EMG signal models have been proposed in the literature. This study focuses on the proposal of a method for modelling surface EMG signals, using either analytical or numerical descriptions of the volume conductor for space-invariant systems, and on the development of advanced models of the volume conductor by numerical approaches, accurately describing the volume conductor geometry and the conductivity, as mainly done in the past, but also the conductivity tensor of the muscle tissue. For volume conductors that are space-invariant in the direction of source propagation, the surface potentials generated by any source can be computed by one-dimensional convolutions, once the volume conductor transfer function has been derived (analytically or numerically). Conversely, more complex volume conductors require a complete numerical approach. In a numerical approach, the conductivity tensor of the muscle tissue should be matched with the fibre orientation. In some cases (e.g. multi-pinnate muscles), accurate description of the conductivity tensor can be very complex. A method for relating the conductivity tensor of the muscle tissue, to be used in a numerical approach, to the curve describing the muscle fibres is presented and applied to investigate representatively a bi-pinnate muscle with rectilinear and curvilinear fibres. The study thus proposes an approach for surface EMG signal simulation in space invariant systems, as well as new models of the volume conductor using numerical methods.

Multichannel surface EMG for the non-invasive assessment of the anal sphincter muscle.

BACKGROUND/AIMS: This work focuses on recording, processing and interpretation of multichannel surface EMG detected from the external anal sphincter muscle. The aim is to describe the information that can be extracted from signals recorded with such a technique. METHODS: The recording of many signals from different locations on a muscle allows the extraction of additional information on muscle physiology and anatomy with respect to that obtained by classic bipolar recordings. Multichannel EMG methods have been recently developed for the assessment of the external anal sphincter. An anal probe was used in this study to record signals at different depths within the anal canal during contractions at different effort levels. The plug is 150 mm in length and 14 mm in diameter, holding a circumferential array of 16 equally spaced silver bar electrodes, located at a distance of 20 mm from the probe tip and aligned with the probe axis. RESULTS: Information about localization of the innervation zone, fiber length, EMG amplitude, muscle fiber conduction velocity and single motor unit analysis can be obtained from the signals recorded with the circumferential array by means of innovative signal processing techniques. CONCLUSIONS: The type of information extracted from multichannel surface EMG signals cannot be obtained with other currently available techniques. The technological innovation described in this work is promising for a further insight into the investigation of pelvic floor pathologies and rehabilitation treatments.

Comparison of spatial filter selectivity in surface myoelectric signal detection: influence of the volume conductor model.

Spatial filters are used for increasing selectivity in surface EMG signal detection. The study investigated the importance of the description of the volume conductor to the inference of conclusions on comparing filter selectivity from simulation analyses. A cylindrical multi-layer description of the volume conductor was used for the simulation analysis. Different anatomies were analysed with this model, and results on filter selectivity were compared. The longitudinal single (LSD), double (LDD) and normal double differential (Laplacian, NDD) filters were investigated. Largely different conclusions could be drawn when comparing filter selectivity resulting from simulations with different volume conductor models. A filter that performed best with a particular anatomy could be the poorest with another anatomy. With a bone-muscle model and superficial fibres, the ratio between peak-to-peak values of the propagating and non-propagating signal components was approximately 220% for LDD and LSD and lower than for NDD (approximately 290%). With a bone-muscle-fat-skin model, LSD performed significantly worse (150%) than both LDD and NDD, which showed similar performances (approximately 300%). Similarly, if the lateral distance of the recording was increased by 10 degrees, the signal amplitude was reduced to 2% with LSD and LDD and to 4% with NDD. With another anatomy, LSD and LDD reduced signal amplitude to 20-25%, and NDD reduced it to 4%. Similar considerations could be drawn for other selectivity indexes. Thus, modelling should be used carefully to infer conclusions on spatial selectivity and to indicate particular choices of spatial filters.