Topographical map of innervation zones within single motor unitsmeasured with a grid surface electrode

A grid electrode with closely spaced contacts was developed to measure the spatial distribution of motor unit action potentials on a skin surface. This electrode was used to estimate the configuration of innervation zones. The derived action potentials showed waveforms originating from the middle length of the muscle belly and propagating bidirectionally along the muscle fibers. The position of the innervation zones was estimated from the source of the propagation. By isolating action potentials from single motor units, the innervation zones for individual motor units were defined. Studies were performed in the biceps brachii of three normal subjects. The innervation zones comprised one to four separate regions and spread over up to 20 mm along the length of muscle fibers. The number and the spreading area of innervation zones varied with the subject and with the motor unit     

Modeling of surface myoelectric signals. II. Model-based signalinterpretation

Experimental electromyogram (EMG) data from the human biceps brachii were simulated using the model described in [10] of this work. A multichannel linear electrode array, spanning the length of the biceps, was used to detect monopolar and bipolar signals, from which double differential signals were computed, during either voluntary or electrically elicited isometric contractions. For relatively low-level voluntary contractions (10%-30% of maximum force) individual firings of three to four-different motor units were identified and their waveforms were closely approximated by the model. Motor unit parameters such as depth, size, fiber orientation and length, location of innervation and tendonous zones, propagation velocity, and source width were estimated using the model. Two applications of the model are described. The first analyzes the effects of electrode rotation with respect to the muscle fiber direction and shows the possibility of conduction velocity (CV) over- and under-estimation. The second focuses on the myoelectric manifestations of fatigue during a sustained electrically elicited contraction and the interrelationship between muscle fiber CV, spectral and amplitude variables, and the length of the depolarization zone. It is concluded that a) surface EMG detection using an electrode array, when combined with a model of signal propagation, provides a useful method for understanding the physiological and anatomical determinants of EMG waveform characteristics and b) the model provides a way for the interpretation of fatigue plots.     

Investigation of optimum electrode locations by using an automatized surface electromyography analysis technique.

Identification of the innervation zone is widely used to optimize the accuracy and precision of noninvasive surface electromyography (EMG) signals because the EMG signal is strongly influenced by innervation zones. However, simply structured fusiform muscle, such as biceps brachii muscle, has been employed mainly due to the simplicity with which the propagation from raw EMG signals can be observed. In this study, the optimum electrode location (OEL), free from innervational influence, was investigated by the propagation pattern of action potentials for brachii muscles and more complicated deltoid muscle structures using an automatized signal analysis technique. The technique employed newly developed computer software with additional clinical uses and minimized subjective differences. EMG signals were recorded using surface array electrodes during voluntary isometric contractions obtained from 12 healthy male subjects. Peaks in EMG signals were detected and averaged for each muscle. The propagation patterns and OEL were examined from biceps brachii muscles for all subjects and from deltoid muscles for seven subjects. The estimated locations were partially confirmed by comparing the root mean squares of the EMG signals. These results show that propagation patterns and OEL could be estimated simply and automatically even from the surface EMG signals of deltoid muscles.     

On Modeling the Single Motor Unit Action Potential

A model has been proposed for the generation of single motor unit potentials routinely observed in the clinical EMG examination of the normal biceps brachii muscle. A dipole representation was chosen for the single fiber activity. The motor unit was constructed from a uniform random array of single fibers. Motor unit potentials generated by this array have been observed at various distances both inside and outside the array. The effects of single fiber dipole axial dispersions on the potentials observed at increasing distances from the array have also been investigated. Motor unit potentials generated by the model have been compared with existing data from multielectrode studies in the biceps brachii.     

A novel approach for estimating muscle fiber conduction velocity by spatial and temporal filtering of surface EMG signals

We describe a new method for the estimation of muscle fiber conduction velocity (CV) from surface electromyography (EMG) signals. The method is based on the detection of two surface EMG signals with different spatial filters and on the compensation of the spatial filtering operations by two temporal filters (with CV as unknown parameter) applied to the signals. The transfer functions of the two spatial filters may have different magnitudes and phases, thus the detected signals have not necessarily the same shape. The two signals are first spatially and then temporally filtered and are ideally equal when the CV value selected as a parameter in the temporal filters corresponds to the velocity of propagation of the detected action potentials. This approach is the generalization of the classic spectral matching technique. A theoretical derivation of the method is provided together with its fast implementation by an iterative method based on the Newton's method. Moreover, the lowest CV estimate among those obtained by a number of filter pairs is selected to reduce the CV bias due to nonpropagating signal components. Simulation results indicate that the method described is less sensitive than the classic spectral matching approach to the presence of nonpropagating signals and that the two methods have similar standard deviation of estimation in the presence of additive, white, Gaussian noise. Finally, experimental signals have been collected from the biceps brachii muscle of ten healthy male subjects with an adhesive linear array of eight electrodes. The CV estimates depended on the electrode location with positive bias for the estimates from electrodes close to the innervation or tendon regions, as expected. The proposed method led to significantly lower bias than the spectral matching method in the experimental conditions, confirming the simulation results.     

Selective activation of distant nerve by surface electrode array

Neural prostheses for restoring lost functions can benefit from selective activation of nerves with limited number and density of electrodes. Here, we show by simulations and animal experiments that multipoint simultaneous stimulation with a surface electrode array can selectively activate nerves in a bundle at a desired location in between the array or at a desired depth, which are referred to as lateral or depth-wise gating stimulation, respectively. The stimulation broadly generates action potentials with cathodic source electrodes, and simultaneously blocks unnecessary propagation with downstream anodic gate electrodes. In general, stimulation with a small diameter electrode can affect a nearly hemispherical region, while a large electrode is effective at a more vertically compressed region, i.e., a surface of nerve bundle. The gating stimulation takes advantage of the size effects by utilizing an asymmetrical electrode array. The array of the lateral gating stimulation is designed to have four electrodes; a pair of large source electrodes and a pair of small gate electrodes. The depth-wise gating stimulation array consists of two electrodes; a large gate and small source electrodes. The simulation first demonstrated that appropriate combination of currents at the source and gate electrodes can change recruitment patterns of nerves with lateral or depth-wise selectivity as desired. We then applied the lateral gating stimulation on the rat spinal cords and obtained a preliminary support for the feasibility.     

The nerve-electrode interface of the cochlear implant: currentspread

One of the fundamental facets of the cochlear implant that must be understood to predict accurately the effect of an electrical stimulus on the auditory nerve is the nerve-electrode interface. One aspect of this interface is the degree to which current delivered by an electrode spreads to neurons distant from it. This paper reports a direct mapping of this current spread using recordings from single units from the cat auditory nerve. Large variations were seen in the degree to which the different units are selective in responding to electrodes at different positions within the scala tympani. Three types of units could be identified based on the selectiveness of their response to the different electrodes in a linear array. The first type of unit exhibited a gradual increase in threshold as the stimulating site was moved from more apical to more basal locations within the scala tympani. The second type of unit exhibited a sharp local minimum, with rapid increases in threshold in excess of 6 dB/mm in the vicinity of the minimum. At electrode sites distant from the local minima the rate of change of the threshold approached that of the first type of units. The final type of unit also demonstrated a gradual change in threshold with changing electrode position, however, two local minima, one apical and one basal, could be identified. These three types are hypothesized to correspond to units which originate apical to the electrode array, along the electrode array and basal to the electrode array     

Construction and validation of a plunge electrode array for three-dimensional determination of conductivity in the heart.

The heart's response to electrical shock, electrical propagation in sinus rhythm, and the spatiotemporal dynamics of ventricular fibrillation all depend critically on the electrical anisotropy of cardiac tissue. Analysis of the microstructure of the heart predicts that three unique intracellular electrical conductances can be defined at any point in the ventricular wall; however, to date, there has been no experimental confirmation of this concept. We report the design, fabrication, and validation of a novel plunge electrode array capable of addressing this issue. A new technique involving nylon coating of 24G hypodermic needles is performed to achieve nonconductive electrodes that can be combined to give moderate-density multisite intramural measurement of extracellular potential in the heart. Each needle houses 13 silver wires within a total diameter of 0.7 mm, and the combined electrode array gives 137 sites of recording. The ability of the electrode array to accurately assess conductances is validated by mapping the potential field induced by a point current source within baths of saline of varying concentration. A bidomain model of current injection in the heart is then used to test an approximate relationship between the monodomain conductivities measured by the array, and the full set of bidomain conductivities that describe cardiac tissue.     

Automatic Decomposition of the Clinical Electromyogram

We describe a new, automatic signal-processing method (ADEMG) for extracting motor-unit action potentials (MUAP's) from the electromyographic interference pattern for clinical diagnostic purposes. The method employs digital filtering to select the spike components of the MUAP's from the background activity, identifies the spikes by template matching, averages the MUAP waveforms from the raw signal using the identified spikes as triggers, and measures their amplitudes, durations, rise rates, numbers of phases, and firing rates. Efficient new algorithms are used to align and compare spikes and to eliminate interference from the MUAP averages. In a typical 10-s signal recorded from the biceps brachii muscle using a needle electrode during a 20 percent-maximal isometric contraction, the method identifies 8-15 simultaneously active MUAP's and detects 30-70 percent of their occurrences. The analysis time is 90 s on a PDP-11/34A.     

Sensitivity- and effort-gain analysis: multilead ECG electrode array selection for activation time imaging

Methods for noninvasive imaging of electric function of the heart might become clinical standard procedure the next years. Thus, the overall procedure has to meet clinical requirements as an easy and fast application. In this paper, we propose a new electrode array which improves the resolution of methods for activation time imaging considering clinical constraints such as easy to apply and compatibility with routine leads. For identifying the body-surface regions where the body surface potential (BSP) is most sensitive to changes in transmembrane potential (TMP), a virtual array method was used to compute local linear dependency (LLD) maps. The virtual array method computes a measure for the LLD in every point on the body surface. The most suitable number and position of the electrodes within the sensitive body surface regions was selected by constructing effort gain (EG) plots. Such a plot depicts the relative attainable rank of the leadfield matrix in relation to the increase in number of electrodes required to build the electrode array. The attainable rank itself was computed by a detector criterion. Such a criterion estimates the maximum number of source space eigenvectors not covered by noise when being mapped to the electrode space by the leadfield matrix and recorded by a detector. From the sensitivity maps, we found that the BSP is most sensitive to changes in TMP on the upper left frontal and dorsal body surface. These sensitive regions are covered best by an electrode array consisting of two L-shaped parts of approximately 30 cm x 30 cm and approximately 20 cm x 20 cm. The EG analysis revealed that the array meeting clinical requirements best and improving the resolution of activation time imaging consists of 125 electrodes with a regular horizontal and vertical spacing of 2-3 cm.     

The Effect of Electromyographic Jitter on Single Motor Unit EMG Potentials

Variability between successive discharges of the single motor unit potential in the biceps brachii muscle, due to electromyographic (EMG) jitter, has been investigated. This jitter results from random arrival times of single fiber potentials at the electrode. A computer model has been used to generate single motor unit potentials incorporating the effects of EMG jitter. A computed variance peak was found in the fast rising edge of the motor unit potential for electrode sites outside of the motor unit territory. This peak was also observed in experimental data recorded from human subjects. The peak variance outside of the motor unit territoxy has also been mathematically related to the number of fibers in the motor unit, jitter, and the slope of the mean action potential at the center of the fast rising edge.     

Packaged AlGaAs waveguide modulator array at 830 nm wavelength

We report a fully packaged AlGaAs waveguide modulator array with four individually addressable elements operating at approximately 830 nm wavelength and a clock speed of 1 GHz. The modulators rely largely on the linear electro-optic effect for operation, and have been packaged with an E/D MESFET driver with complementary 3.5 V outputs, and a thick-film ceramic bias network. The device is compact, using multimode interference devices for on-chip splitters and combiners and has a 4 mm electrode length. Extinction ratios in excess of 10 dB have been demonstrated over a temperature range from room temperature to 143°C. The modulator array has been packaged with a remote high power (100 mW) diode laser using stable single mode input coupling, while the modulator output is packaged with a multimode fiber array of 52.5 μm core diameter     

基于任意麦克风阵列的声源二维DOA估计算法研究

对基于麦克风阵列的声源定位技术进行了研究,分析了基于麦克风阵列的远场信号模型,并结合子空间的方法推导出了声源二维(水平角和俯仰角)DOA估计——2D-MUSIC算法,该算法适用于任意拓扑结构的麦克风阵列。利用MATLAB仿真工具,对几种典型阵列结构进行了对比分析,提出了2种新型的三维麦克风阵列:均匀球面阵和三维均匀直线阵。仿真结果表明,提出的DOA估计算法在二维的均匀圆阵、三维的均匀球面阵和三维均匀直线阵中,均能得到较好的DOA估计效果。     

一种基于卡尔曼的红外相机动态范围自适应调整算法

卡尔曼算法在目标跟踪、风向预测、非均匀性矫正领域都有着广泛应用。提出了用卡尔曼算法进行线列相机扫描灰度值的预测。根据线列探测器的实时扫描,得出了灰度直方图的变化,并将这一过程归一化到卡尔曼公式中;用递推法快速测出了下一时间灰度的统计分布,间接预测了积分时间,最终达到了红外线列探测器动态范围自适应控制的目的。采用多幅8 bit图像为不同场景提供仿真环境,最终快速稳定地预测出了下一阶段的灰度分布统计图。结果为红外线列相机扫描提供了智能化解决方案。     

Polarity Correlator for Conduction Velocity Measurement

Techniques for the estimation of skeletal muscle fiber conduction velocity are of considerable interest. These techniques use, in general, some form of cross correlation or zero-crossing analysis. Cross correlation is a straightforward method of conduction velocity estimation, however, it is difficult to realize low-cost real-time processors. Polarity correlation offers a method which preserves the advantages of cross correlation and satisfies the speed-cost constraint. This paper describes the algorithm for implementation of a polarity cofrelator instrument using a microprocessor. The instrument is tested with deterministic and stochastic signals and used to estimate the conduction velocity of biceps brachii fibers of ten normal subjects. The performance of the instrument is compared to the zero-crossing technique described by Lynn.     

A silicon-based, three-dimensional neural interface: manufacturingprocesses for an intracortical electrode array

A method has been developed for the manufacture of a "three-dimensional" electrode array geometry for chronic intracortical stimulation. This silicon based array consists of a 4.2 x 4.2 x 0.12 mm thick monocrystalline substrate, from which project 100 conductive, silicon needles sharpened to facilitate cortical penetration. Each needle is electrically isolated from the other needles, and is about 0.09 mm thick at its base and 1.5 mm long. The sharpened end of each needle is coated with platinum to facilitate charge transfer into neural tissue. The following manufacturing processes were used to create this array. 1) Thermomigration of 100 aluminum pads through an n-type silicon block. This creates trails of highly conductive p+ silicon isolated from each other by opposing pn junctions. 2) A combination of mechanical and chemical micromachining which creates individual penetrating needles of the p+ silicon trails. 3) Metal deposition to create active electrode areas and electrical contact pads. 4) Array encapsulation with polyimide. The geometrical, mechanical, and electrical properties of these arrays should make them well suited as interfaces to cortical tissue.     

Proposal of a new method for narrowing and moving the stimulated region of cochlear implants: animal experiment and numerical analysis

We have proposed the tripolar electrode stimulation method (TESM) for narrowing the stimulation region and continuously moving the stimulation site for cochlear implants. The TESM stimulates the auditory nerve array using three adjacent electrodes which are selected among the electrodes of an electrode array within the lymphatic fluid. Current is emitted from each of the two lateral electrodes and received by the central electrode. The current received by the central electrode is made equal to the sum of the currents emitted from the lateral electrodes. In this paper, we evaluate whether or not TESM works according to a theory which is based on numerical analysis using an electrical equivalent circuit model of the auditory nerve fibers. In this simulation, the sums of the excited model fibers are compared to the compound action potentials (CAP's) which we obtained through animal experiments. To identify the main parameter while maintaining the amplitude of the CAP (the sum of the fired fibers), we assumed the presence of some parameters from the radial current density profile. In the case of the width value among the parameters being kept constant, the amplitude of the CAP was almost constant; thus, the number of the fired fibers was also almost constant. The width value equals the distance between the points at which the profile of the radial current density of the electrode array and the line of the radial threshold current density of the electrode array intersect. It is possible to determine the measure of the stimulation region or site by controlling the width value and the ratios of the currents emitted from the lateral electrodes. As a result, we succeeded in narrowing the stimulation region by controlling the sum of the currents emitted from the two lateral electrodes. Also we succeeded in continuously moving the stimulation site by modifying the currents emitted from the two lateral electrodes.     

Measurements of the Uptake Area of Small-Size Electromyographic Electrodes

Extracellular action potentials from 76 different muscle fibers in the human brachial biceps were recorded, with a 14 lead multielectrode, each leading-off surface being 25 ?m in diameter. Volume conduction of these action potentials was calculated by representing the low-pass filter characteristics of the muscle tissue by a transfer function with one time constant and an attenuation factor. The radial decline of the action potentials was calculated in steps of 5 ?m, and the pickup radius of the electrode, defined as the distance at which the peak-to-peak amplitude of the action potentials declines to 200 ?V, was computed. The pickup radius for the 25 ?m diameter electrode, assuming an average action potential peak-to-peak amplitude of 6.2 mV was 292 Mm. With this uptake area, a fiber density in the brachial biceps of 1.37 fibers/uptake area (the average number of fibers belonging to one motor unit and included within the electrode uptake area) and a fiber radius of 28 ?m, the electrode "sees" 34 different motor units.     

Motion Artifact from Spot and Band Electrodes During Impedance Cardiography

We have modified impedance cardiography for monitoring cardiac output during stress tests. Employing an off-line microcomputer, our instrument ensemble averaged impedance signals to minimize the effect of motion artifacts. We proposed a new four-spot electrode array and replaced the usual encircling band electrode array with it. We tested ten normal subjects and compared the signal-to-noise ratio (SNR) from our spot electrode array to that from a typical band electrode array at rest and during four levels of exercise on a treadmill. The average of the sighal-to-noise ratios for ten subjects from our spot electrode array was 13.6-45.5 percent larger than that from the band electrode array at rest and during four levels of exercise. Thus, it is desirable to replace band electrodes with spot electrodes in impedance cardiography for exercising subjects.