Skip to main content
  • Research article
  • Open access
  • Published:

Detection of stretch reflex onset based on empirical mode decomposition and modified sample entropy

Abstract

Background

Accurate spasticity assessment provides an objective evaluation index for the rehabilitation treatment of patients with spasticity, and the key is detecting stretch reflex onset. The surface electromyogram of patients with spasticity is prone to false peaks, and its data length is unstable. These conditions decrease signal differences before and after stretch reflex onset. Therefore, a method for detecting stretch reflex onset based on empirical mode decomposition denoising and modified sample entropy recognition is proposed in this study.

Results

The empirical mode decomposition algorithm is better than the wavelet threshold algorithm in denoising surface electromyogram signal. Without adding Gaussian white noise to the electromyogram signal, the stretch reflex onset recognition rate of the electromyogram signal before and after empirical mode decomposition denoising was increased by 56%. In particular, the recognition rate of stretch reflex onset under the optimal parameter of the modified sample entropy can reach up to 100% and the average recognition rate is 93%.

Conclusions

The empirical mode decomposition algorithm can eliminate the baseline activity of the surface electromyogram signal before stretch reflex onset and effectively remove noise from the signal. The identification of stretch reflex onset using combined empirical mode decomposition and modified sample entropy is better than that via modified sample entropy alone, and stretch reflex onset can be accurately determined.

Background

Muscle spasticity is an intermittent or persistent involuntary excessive movement of the skeletal muscle caused by an upper motor neuron injury [1]. Spasticity is clinically manifested by an increase in passive stretching resistance, i.e., an increase in muscle tension, and resistance increases with passive stretching speed [2]. Accurate spasticity assessment provides an objective evaluation index for the rehabilitation treatment of patients with spasticity [3]. The modified Ashworth scale is the most widely used method for assessing spasticity in the clinical setting. It is simple and easy to implement without the aid of instruments. However, the modified Ashworth scale exhibits strong subjectivity; thus, it cannot meet the requirements of accurate spasticity assessment [4].

The tonic stretch reflex threshold is currently recognized as the most effective and consistent value for assessing spasticity [5]. The stretch reflex threshold represents the joint angle at which antagonist muscles or motor neurons begin to contract when the affected limbs of the subject are stretched passively [6]. The starting point that corresponds to the beginning of contraction of antagonistic muscles or motor neurons is called stretch reflex onset [7]. Surface electromyogram (EMG) is a nonlinear and nonstationary signal obtained from electrodes attached on the skin surface of muscles [8]. The spasticity of a patient can be analyzed, and the nerve components during an increase in muscle tension can be distinguished by analyzing surface EMG signals. Stretch reflex onset detection using surface EMG signal is a prerequisite and basic step in biomedical research and clinical diagnosis, such as gait recognition, and automatic prosthetic control [9,10,11]. It is crucial for detecting the stretch reflex threshold of muscles when assessing spasticity [12,13,14].

At present, the commonly used stretch reflex onset detection method is based on experience or the standard deviation (SD) of surface EMG signal [1, 15]. Experience is too subjective [12], and the SD of surface EMG signal is susceptible to the baseline activity of such signal [16]. SD cannot achieve exact onset detection give the low signal-to-noise ratio (SNR) of surface EMG signal in patients with spasticity [17]. Detection methods based on time-domain analysis and statistical characteristics, such as the Teager–Kaiser energy operator and the maximum likelihood method, have been effectively applied in recent years. The Teager–Kaiser energy operator considers the amplitude and instantaneous frequency of signals to improve the recognition rate of onset detection; however, it requires an extremely high SNR [18]. The maximum likelihood method is commonly used for the onset detection of signals with low SNR (SNR = 1); however, it is only suitable for processing Gaussian noise [19]. Yang et al. combined Teager–Kaiser energy operator with an image enhancement technique and adopted morphological close operators to realize the accurate onset detection of pathological surface EMG signals with weak noise [20].

Sample Entropy is a method based on approximate entropy that is used to measure the complexity of time series. Sample entropy has been applied in evaluating the complexity of physiological time series and in diagnosing pathological state [21]. Zhou et al. used sample entropy in the onset detection of muscle activities and achieved good recognition effect in an environment with background burr noise [22]. However, Sample entropy is highly dependent on data length, and stability is hardly guaranteed when data length is short [23,24,25,26]. Compared with the standard sample entropy, the modified sample entropy exhibits stronger dependence on data length and has a smaller fluctuation in entropy value, making the former more suitable for processing short-term surface EMG signals [27]. Considering the poor quality of surface EMG signals in patients with spasticity, preprocessing should be performed before the stretch reflex onset detection of the modified sample entropy [28].

Empirical mode decomposition is the adaptive decomposition of a signal into several intrinsic mode functions in accordance with the time-scale characteristics of data without presetting any basis function [29]. The empirical mode decomposition method proposed by Huang et al. can be applied to any type of signal decomposition and exhibits evident advantages in dealing with nonlinear and nonstationary EMG signals [30, 31]. Zhou et al. combined empirical mode decomposition with soft thresholds to eliminate three types of common noise in surface EMG signal: power line interference, Gaussian white noise, and baseline drift. The empirical mode decomposition denoising has been proven to be better than those of other digital filters, such as traditional infinite impulse response filters [32].

Given the pathological symptoms of patients with spasticity, abnormal overexcitement leads to discharges of motor units, and consequently poor quality of surface EMG signals. The surface EMG of patients with spasticity is prone to unconscious and involuntary false peaks, which cause a slight difference in surface EMG signal before and after the stretch reflex onset. Moreover, a false EMG peak is mixed into the mutant signal after the stretch reflex onset [33]. Therefore, a stretch reflex onset detection method based on empirical mode decomposition denoising and modified sample entropy recognition is proposed in the current study. First, the denoising effects of the empirical mode decomposition and wavelet threshold algorithms are compared through various evaluation indices. Then, the surface EMG signal is decomposed via empirical mode decomposition, and the effective intrinsic mode function is extracted in accordance with the correlation coefficient between each order of the intrinsic mode function and the original signal. The soft threshold is set to denoise the signal on the basis of the surface EMG signal of the subjects in resting state. Lastly, stretch reflex onset is identified using the modified sample entropy. Gaussian white noise with SNRs of 0, 5, 10, 15 and 20 dB is added. The accuracy and robustness of this method and those of the modified sample entropy are compared.

Results

Analysis of empirical mode decomposition denoising effect

The semi-synthetic surface EMG signals of 10 patients with spasticity were selected and compared with those of empirical mode decomposition denoising via the wavelet threshold denoising method. Among these signals, the wavelet basis selected the sym4 basis function, which is relatively similar to the surface EMG signals. Three evaluation indices were cited to compare the denoising effects of the two methods, namely, SNR, root mean square error (R), and the correlation coefficient (P).

Where SNR, R, and P are expressed as follows:

$$ SNR=101g\left(\frac{\sum \limits_{i=1}^N{x}_i(t)}{\sum \limits_{i=1}^N{\left({x}_i(t)-{y}_i(t)\right)}^2}\right), $$
(1)
$$ R=\sqrt{\frac{1}{N}\sum \limits_{i=1}^N{\left({x}_i(t)-{y}_i(t)\right)}^2}, $$
(2)
$$ P=\frac{\operatorname{cov}\left({x}_i(t),{y}_i(t)\right)}{\sqrt{D\left({x}_i(t)\right).\kern0.5em D\left({y}_i(t)\right)}} $$
(3)

xi(t) is the original signal, yi(t) is the signal after de-noising, and N is the number of sampling points. Cov() is the covariance function, D() is the variance function.

When SNR is higher, the signal denoising effect is better; when the value of R is smaller, the coincidence degree between the denoised and original signals is higher; when the value of P is larger, the correlation between the denoised and original signals is stronger. The semi-synthetic surface EMG signals of the 10 subjects were selected, and the evaluation indices were obtained using the two denoising methods (Table 1).

Table 1 The evaluation indexes of empirical mode decomposition and wavelet threshold denoising

As shown in Table 1, the SNR and P of the empirical mode decomposition denoised signals are higher, R is smaller, the P value of wavelet threshold denoising fluctuates considerably, and denoising performance is unstable. The experimental results indicate that the empirical mode decomposition algorithm is better than the wavelet threshold algorithm in denoising surface EMG signal. The former exhibits stronger denoising ability and better effect, and it can be used as an effective surface EMG signal denoising method.

Stretch reflex onset detection after denoising of EMG signals

The semi-synthetic surface EMG signals of 25 patients with spasticity were selected, and the data of each patient were added with 0, 5, 10, 15, and 20 dB Gaussian white noise. The stretch reflex onset detection process based on empirical mode decomposition denoising and modified sample entropy recognition is described as follows.

  1. (1)

    Soft threshold denoising of empirical mode decomposition.

  2. (2)

    The denoised surface EMG signal is processed via frame processing with an active sliding window of fixed length. The frame is shifted to one point to calculate the modified sample entropy of each frame.

  3. (3)

    The entropy value is calculated. If the corresponding entropy value at a certain time is greater than the set soft threshold and 50 consecutive points are greater than the soft threshold, then the moment is judged as the stretch reflex onset.

The sliding window length is 32 ms under the optimal parameters of the modified sample entropy, and the threshold sensitivity value is 0.3–0.45. The sliding window length is 64 ms under the optimal parameters of the modified sample entropy, and the threshold sensitivity value is 0.5–0.6. The sliding window length is 96 ms under the optimal parameters of the modified sample entropy, and the threshold sensitivity value is 0.5–0.65. The stretch reflex onset recognition rate of the modified sample entropy in a noisy environment is provided in Table 2. The stretch reflex onset recognition rate of the modified sample entropy after empirical mode decomposition denoising of the surface EMG signals is presented in Table 3. After empirical mode decomposition denoising of the surface EMG signals, the stretch reflex onset recognition rate was obtained under the optimal parameter of the modified sample entropy is presented in Table 4.

Table 2 Stretch reflex onset recognition rate of modified sample entropy under noise environment
Table 3 Stretch reflex onset recognition rate of modified sample entropy after denoising
Table 4 Stretch reflex onset recognition rate under Optimal Parameters after denoising

As shown in Figs. 1 and 2, the semi-synthetic surface EMG signals of a patient with the modified Ashworth scale level of 1+ were selected, and Gaussian white noise that corresponded to an SNR of 20 dB was added. The modified sample entropy values before and after denoising were compared. The difference in the EMG signals before and after the stretch reflex onset increased after empirical mode decomposition denoising as shown by the signal time-domain and modified sample entropy diagrams. As indicated in Tables 2, 3, and 4, the stretch reflex onset recognition rate gradually increased. The average stretch reflex onset recognition rate of the modified sample entropy in a noisy environment is 35%, and the average stretch reflex onset recognition rate of the modified sample entropy after empirical mode decomposition denoising of the surface EMG signals is 66%. The average stretch reflex onset recognition rate under the optimal parameters of the modified sample entropy is 93% after empirical mode decomposition denoising of the surface EMG signals. Without adding Gaussian white noise to the EMG signals, the stretch reflex onset recognition rate of the EMG signals before and after empirical mode decomposition denoising was increased by 56%. The experimental results showed that the stretch reflex onset recognition rate was significantly increased after empirical mode decomposition denoising. In particular, the stretch reflex onset recognition rate can reach up to 100% and the average recognition rate is 93% under the optimal parameters of the modified sample entropy. However, one study [34] reported that the stretch reflex onset recognition rate of the Teager–Kaiser energy operator is 88%. For the signals with added Gaussian noise, the stretch reflex onset recognition rate was significantly increased after empirical mode decomposition denoising. This result indicates that the empirical mode decomposition algorithm can effectively remove noise interference from surface EMG signals and exhibits good anti-noise performance.

Fig. 1
figure 1

The modified sample entropy under noise (SNR = 20 dB) environment

Fig. 2
figure 2

The modified sample entropy after denoising

Discussion

Objective evaluation indices are provided for the rehabilitation treatment of patients with spasticity through the stretch reflex onset recognition method of combined empirical mode decomposition denoising and modified sample entropy recognition. The empirical mode decomposition algorithm can remove noise from surface EMG signals and retain useful information, and the soft threshold of the empirical mode decomposition algorithm expands the difference between the signals before and after the stretch reflex onset. Although the modified sample entropy can recognize stretch reflex onset, the poor quality of surface EMG signals in patients with spasticity leads to a slight difference in surface EMG signals before and after the stretch reflex onset. Moreover, false EMG peaks are mixed into the mutation signal after the stretch reflex onset. Thus, EMG signals must be preprocessed before the stretch reflex onset detection. For the surface EMG signals of patients with spasticity, the performance of the stretch reflex onset detection method using combined empirical mode decomposition and the modified sample entropy algorithm is better than that of the modified sample entropy algorithm alone. The former is more suitable for clinical applications.

Analysis of results

We evaluated the denoising effect of the empirical mode decomposition and wavelet threshold algorithms through three evaluation indicators. The SNR and P values after empirical mode decomposition denoising were higher than the results after wavelet threshold denoising. The root mean square error value was lower than the results after wavelet threshold denoising. The experimental results show that the empirical mode decomposition algorithm is better than the wavelet threshold algorithm in denoising surface EMG signals. The former exhibits stronger denoising ability and better effect, and it can be used as an effective surface EMG signal denoising method. By using the stretch reflex onset detection method of combined empirical mode decomposition denoising and modified sample entropy recognition, the stretch reflex onset recognition rate can reach up to 100% and the average recognition rate is 93% under the optimal parameters of the modified sample entropy. The stretch reflex onset recognition rate of combined empirical mode decomposition denoising and modified sample entropy recognition is considerably higher than that of modified sample entropy recognition alone. The combination of empirical mode decomposition and the modified sample entropy can be applied to detect stretch reflex onset in patients with spasticity.

Study limitations

This study only analyzed the stretch reflex onset of surface EMG signals. The quantitative assessment of upper limb flexor spasticity is determined on the basis of the tension stretch reflex threshold, which is determined on the basis of surface EMG and joint angle signals [13, 14]. The tensile stretch reflex threshold depends on speed. In this experiment, a certain velocity fluctuation occurs when the doctor uniformly stretches a patient’s limbs. We plan to design an isokinetic instrument that can measure spasticity and eliminate the influence of velocity fluctuation. Surface EMG signals are used in this study to analyze the stretch reflection threshold. We can also add a linear acceleration sensor in the experiment to determine the stretch reflection threshold through multiple parameters. The sample size used in this experiment is small. In the future, the sample size of patients with different spasticity grades will be increased, and experiments will be conducted at different stretch speeds to further comprehensively verify the reliability of the stretch reflex onset detection based on combined empirical mode decomposition and the modified sample entropy.

Clinical implications

The upper limb spasticity assessment system based on surface EMG signals is simple and convenient to use in spasticity assessment in the clinical setting. The empirical mode decomposition and modified sample entropy algorithm exhibit satisfactory clinical applicability for the pathological surface EMG signals of patients with spasticity. Accurate spasticity assessment provides an objective evaluation index for the rehabilitation treatment of patients with spasticity. On the one hand, it can objectively reflect the rehabilitation status of stroke survivors. On the other hand, it can assist in adjusting treatment programs and formulating personalized rehabilitation programs. The realization of clinically accurate stretch reflex onset and the measurement of stretch reflex intensity can help understand the spasticity and pathological mechanism of increased muscle tension, track the rehabilitation of patients with spasticity, and screen targeted reasonable treatment programs.

Conclusion

A stretch reflex onset detection method based on combined empirical mode decomposition and modified sample entropy was proposed to address the problem of EMG signals in patients with spasticity. On the basis of empirical mode decomposition denoising and modified sample entropy recognition experiments, the following conclusions were drawn from this study. 1) The empirical mode decomposition algorithm is better than the wavelet threshold algorithm in denoising surface EMG signals because it can retain useful information and effectively remove noise from surface EMG signals. 2) Gaussian white noise with different SNRs was added to the original signal, and the results showed that empirical mode decomposition exhibits good anti-noise performance and its denoising capability is helpful in the stretch reflex onset detection. The stretch reflex onset recognition rate of combined empirical mode decomposition denoising and modified sample entropy recognition is considerably higher than that of the modified sample entropy alone. The stretch reflex onset recognition rate under the optimal parameters of the modified sample entropy can reach up to 100% and the average recognition rate is 93%. In summary, the empirical mode decomposition algorithm can effectively remove noise from surface EMG signals, and the stretch reflex onset detection rate of combined empirical mode decomposition and the modified sample entropy is high. Thus, accurate and reliable stretch reflex onset detection can be realized, and a new method for the objective and accurate evaluation of spasticity is provided.

Methods

Experimental

The surface EMG signals of biceps brachii during the passive extension of the upper limbs were collected from patients with spasticity. The subject sat upright while the doctor passively stretched the upper limb of the affected side of the subject, allowing the upper limb to adapt to the traction speed. Sudden traction was prevented from causing tension in the subjects that might affect experimental data. At the end of the preparatory activity, the doctor collected the surface EMG signals of each subject at an appropriate uniform speed on the basis of experience, and the process was repeated four times at 2 min intervals. Reacquisition was performed after 1 day. The elbow joint was fixed with one hand, and the wrist joint was held with the other hand. The passive traction of the subject from the maximum bending angle to the maximum extension angle was completed. Spasticity is velocity-dependent; thus, the fluctuation of the stretching speed must be maintained within a certain range during the experiment.

The EMG signal sensor used in the experiment was based on the hardware circuit and software design of a single-chip microcomputer (STM32). The collected analog signals were A/D converted and then transmitted to the upper computer. The data sampling frequency was 1000 Hz. The positive and negative electrodes of the electrode sheet were attached to the biceps muscle at 2 cm intervals. The reference end was attached to the side of the biceps. The experimental procedure is illustrated in Fig. 3.

Fig. 3
figure 3

The evaluation process of Spasticity

A total of 25 eligible stroke survivors from the Anhui Provincial Hospital were selected for the experiment. The study was approved by the Anhui Provincial Hospital Ethics Committee, and all the participants provided informed consent. Among the subjects, 11 belonged to the modified Ashworth scale 1, 7 to the modified Ashworth scale 1+, and 7 to the modified Ashworth scale 2.

Inclusion criteria: patients with upper limb flexor spasticity caused by stroke; can sit independently; an elbow flexion and extension range of at least 90°; clear consciousness, without serious cognitive and audiovisual impairment, can cooperate with the examination to complete simple instructions.

Exclusion criteria: patients with central nervous system disorders, such as multiple sclerosis and other diseases that may cause spasticity of the limbs; patients with diseases that can affect the movement of the elbow joints of the upper limbs, such as fracture of the upper limbs; patients with the modified Ashworth scale grades of 0, 3, and 4.

Empirical mode decomposition

The empirical mode decomposition method proposed by N.E. Huang is based on the time-scale characteristics of data for signal decomposition without setting any basis function [31]. The original signal is decomposed into a number of narrow-band components, and each component is called intrinsic mode function. The decomposition result consists of several intrinsic mode functions and one residual signal.

$$ s(t)=\sum imf\kern0.5em {}_i(t)\kern0.5em +{r}_n(t), $$
(4)

Each intrinsic mode function should meet the following conditions [31]:

  1. (1)

    The number of extreme points and zero crossings in the entire dataset should not be more than one.

  2. (2)

    The mean values of the upper and lower envelopes should be zero at any point.

However, actual signals are complex, and thus they do not satisfy the aforementioned conditions. Therefore, Huang made the following assumptions [31]:

  1. (1)

    Any signal is composed of several intrinsic mode functions.

  2. (2)

    Each intrinsic mode function can be linear or nonlinear. The numbers of local zeros and extreme points of each intrinsic mode function are the same, and the upper and lower envelopes are locally symmetric with respect to the time axis.

  3. (3)

    Each signal can contain several intrinsic mode functions at any time. If the modal functions overlap with one another, then a composite signal is formed.

The empirical mode decomposition algorithm assumes that a signal is composed of several intrinsic mode functions. First, all the maximum and minimum points in the original signal s(t) are extracted. Then the upper envelope e+(t) and the lower envelope e−(t) are fitted by a cubic spline function. Lastly, the mean values of the upper and lower envelopes are calculated as follows:

$$ m(t)=\kern0.5em \frac{e_{+}(t)+e\_(t)}{2}, $$
(5)

The difference between the original signal and the mean envelope is calculated as follows:

$$ {h}_1^1(t)=s(t)={m}_1(t), $$
(6)

If h1(t) meets the preceding intrinsic mode function conditions, then h1(t) can be used as the first-order intrinsic mode function. If the preceding intrinsic mode function conditions are not met, then h1(t) is used as the input and the aforementioned steps are repeated. Assume that h1k(t) meets the preceding intrinsic mode function conditions after k times. Then, h1k(t) is used as the first-order intrinsic mode function of the original signal as follows:

$$ {c}_1(t)={h}_1^k, $$
(7)

The difference between the original signal s(t) and the first-order intrinsic mode function can be calculated as follows:

$$ {r}_1(t)=s(t)-{c}_1(t), $$
(8)

c1(t) is derived repeatedly to enable r1(t) to obtain a second-order component, i.e., c2(t). Empirical mode decomposition stops when the N-th order intrinsic mode function component cn(t) or the remaining amount rn(t) is less than a preset value, or when the residual component rn(t) is a monotonic function or a constant.

Lastly, empirical mode decomposition decomposes the original signal into x(t) as follows:

$$ x(t)=\sum \limits_{i=1}^n{c}_i(t)+{r}_i(t) $$
(9)

Surface EMG signal denoising based on empirical mode decomposition

Before the stretch reflex onset surface EMG signal of patients with spasticity exhibits involuntary muscle activities and is mixed with noise during the acquisition process. Consequently, the difference in EMG signal before and after the stretch reflex onset is reduced, and the probability of the stretch reflex onset misjudgment is increased. To eliminate the baseline activity and high-frequency noise components of surface EMG signal, the surface EMG signals of the subjects with non-autonomous activities in resting state were set as the soft threshold. Inspired by the soft threshold function in [12], the soft threshold setting rule is presented as follows:

$$ \eta \left({IMF}_{ij}\right)=\mathit{\operatorname{sign}}\left({IMF}_{ij}\right){\left(\left|{IMF}_{ij}\right|-{\lambda}_i\right)}_{+}, $$
(10)

where η (IMFij) denotes the j-th value in the i-th-order intrinsic mode function component after denoising, IMFi represents the i-th-order intrinsic mode function component that should be preprocessed, λi indicates that each intrinsic mode function is filtered with two SDs of the surface EMG signal at rest, and (∙)+ designates the positive portion.

The SD of surface EMG signals is twice larger than that of baseline signals in the stretch reflex onset [13]. The SD of IMFi after the stretch reflex onset is twice greater than that of the intrinsic mode function in resting state. Therefore, pretreatment can effectively eliminate involuntary muscle activities generated during spasticity evaluation, reduce the baseline activity and noise interference of the surface EMG signals before the stretch reflex onset, and expand signal difference before and after the stretch reflex onset.

First, the correlation between each intrinsic mode function component and the original signal s(t) is calculated, and the effective components are selected in accordance with the correlation results for filtering. The specific processes are as follows:

  1. (1)

    First, the original signal is normalized, and then the signal is decomposed via empirical mode decomposition to obtain IMFi.

  2. (2)

    The correlation between intrinsic mode function and the original signal of each order is calculated. In accordance with the magnitude of the correlation, the effective component is found with the soft threshold of one-tenth of the maximum correlation coefficient, and the pseudo-component is eliminated.

  3. (3)

    Filtering is performed on the effective component, and the signal is reconstructed after filtering to reduce the noise of the original signal.

Modified sample entropy

Sample entropy reflects the probability of generating new information in a nonlinear dynamic system. The definition of standard sample entropy, shows that any distance greater than the similar tolerance r is discarded, and only the number of distances less than or equal to the similarity tolerance r is retained. The similarity of vectors is based on the Heaviside function [21, 35,36,37].

The boundary between various types of signals in an actual physical environment is unclear; thus, determining whether the input sample completely belongs to a certain class is difficult. In addition, the distance fluctuation in an actual surface EMG signal sequence is extremely large, and the slight fluctuations of d and r may also cause violent fluctuations in the output entropy value. Therefore, the standard sample entropy is insufficient for processing surface EMG signals. From [24], the function is defined on the basis of the sigmoid function as follows:

$$ f\left(d,r\right)=\frac{1}{1+\exp \left[\left(d-r\right)/r\right]}, $$
(11)

where d is less than r; and the closer to zero, the closer the output value is to 1. d is greater than r; and the greater the difference, the closer the output value is to zero. The modified sample entropy is described as follows [21, 24, 38].

For a given N-point time series x(n) = {x (1), x (2), …, x(N)}, the modified sample entropy algorithm is calculated as follows.

The continuous m values in each frame signal sequence constitute an m-dimensional vector, where Xm(i) = [x(i); x(i + 1); …; x(i + m − 1)], 1 ≤ i ≤ N–m + 1.

The distance \( {d}_{ij}^m \) between the two vectors Xm(i) and Xm(j) (1 ≤ j ≤ N − m, j ≠ i) is defined as follows:

$$ {d}_{ij}^m=d\left[{X}^m(i),{X}^m(j)\right]=\underset{k\in \left[0,m\kern0.5em 1\right]}{\max}\kern0.5em \left\{\left|x\left(i+k\right)-\left(j+k\right)\right|\right\}, $$
(12)

The similarity \( {D}_{ij}^m \) between Xm(i) and Xm(j) can be calculated by providing a similar tolerance r as follows:

$$ {D}_{ij}^m=f\left({d}_{ij}^m,r\right)=\frac{1}{1+\exp \left[\left({d}_{ij}^m\right)-r/r\right]} $$
(13)

The function \( {B}_r^m \) can be calculated as follows:

$$ {B}_r^m(i)=\frac{1}{N-m-1}\sum \limits_{j=1,j\ne i}^{N-m}{D}_{ij}^m, $$
(14)
$$ {B}_r^m=\frac{1}{N-m}\sum \limits_{i=1}^{N-m}{B}_r^m(i), $$
(15)

Similarly, the function \( {A}_r^m \) can be determined by changing the vector dimension m to m + 1 as follows:

$$ {A}_r^m(i)=\frac{1}{N-m-1}\sum \limits_{j=1,j\ne i}^{N-m}{D}_{ij}^{m+1}, $$
(16)
$$ {A}_r^m=\frac{1}{N-m}\sum \limits_{i=1}^{N-m}{A}_r^m(i), $$
(17)

Lastly, the modified sample entropy (mSampEn) can be defined as

$$ mSampEn\left(m,r,N\right)=-1n\left[\frac{A_r^m}{B_r^m}\right] $$
(18)

Stretch reflex onset recognition based on the modified sample entropy

The modified sample entropy depends on N, m, and r. In accordance with the empirical value [22, 24], m is 2 and r is generally (0.1–0.25) SD(X). When X is the data of each frame, r is the local similarity tolerance. When X is the total data, r is the global similarity tolerance. Studies have confirmed that global similarity tolerance is more suitable for surface EMG signal mutation point detection than local similarity tolerance [22]. A previous study [39] verified that the result was optimal when r was 0.25 × SD(X). In conclusion, the current study sets r as the global similar tolerance, and takes 0.25 × SD(X). The sliding window length N ranges from 32 ms to 96 ms, with an interval of 32 ms.

The stretch reflex onset of the surface EMG signal of biceps brachii is determined using the modified sample entropy. The surface EMG signal is divided into frames by using a sliding window, and each frame is shifted to one point. The modified sample entropy of each frame signal is calculated as mSampEn. The adaptive threshold Th is set in accordance with Formula (19). The value of mSampEn lower than Th is set as zero and the value of mSampEn larger than Th is retained. The moment is judged as the stretch reflex onset when the value of mSampEn is greater than zero at a certain time and 50 consecutive mSampEn values are all greater than zero. The threshold sensitivity value α ranges from 0.3 to 0.65, with an interval of 0.05.

$$ Th=\min (mSampEn)+a\left[\max (mSampEn)-\min (mSampEn)\right], $$
(19)

The semi-synthetic surface EMG signal is selected. That is, the stretch reflex onset of the surface EMG signal is known, and thus the stretch reflex onset detection capability of combining empirical mode decomposition denoising with the modified sample entropy is verified. The semi-synthetic surface EMG signals consist of two types. The first type is the surface EMG signal of subjects in resting state, and the second type is the surface EMG signal of subjects with spasticity. The signal duration of each group is 3000 ms. The interval [− 50, + 50] is selected on the basis of the stretch reflex onset that is set in advance. If the detection result is within the interval, then it is considered true detection (TD). If it is beyond the interval, then it is considered false detection (FD). The recognition rate is defined as:

$$ recognition\_ rate=\frac{TD}{TD+ FD} $$
(20)

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request, without breaching participant confidentiality.

Abbreviations

EMG:

Electromyography

SD:

Standard deviation

SNR:

Signal-to-noise ratios

References

  1. Pandyan AD, Gregoric M, Barnes MP, Wood D, Wijck FV, Burridge J, Hermen H, Johnson GR. Spasticity: clinical perceptions, neurological realities and meaningful measurement. Disabil Rehabil. 2005;27(1–2):2–6.

    Article  Google Scholar 

  2. Huang HW, Ju MS, Lin CCK. Flexor and extensor muscle tone evaluated using the quantitative pendulum test in stroke and parkinsonian patients. J Clin Neurosci. 2016;27:48–52.

    Article  Google Scholar 

  3. Valerie S, Diane P. Neurological rehabilitation and the management of spasticity. Med. 2012;40(9):513–7.

    Article  Google Scholar 

  4. Fleuren JFM, Voerman GE, Erren-Wolters CV, Snoek GJ, Rietman JS, Hermens HJ, Nene AV. Stop using the Ashworth scale for the assessment of spasticity. J Neurol Neurosurg Psychiatry. 2010;81(1):46–52.

    Article  Google Scholar 

  5. Germanotta M, Taborri J, Rossi S, Frascarelli F, Palermo E, Cappa P, Castelli E, Petrarca M. Spasticity measurement based on tonic stretch reflex threshold in children with cerebral palsy using the PediAnklebot. Front Hum Neurosci. 2017;11:1–11.

    Article  Google Scholar 

  6. Park HS, Kim J, Damiano DL. Development of a haptic elbow spasticity simulator (HESS) for improving accuracy and reliability of clinical assessment of spasticity. IEEE Trans Neur Sys Rehabil Eng. 2012;20(3):361–70.

    Article  Google Scholar 

  7. Hu BH, Zhang XF, Mu JS, Wu M, Zhu ZJ, Liu ZS, Wang Y. Spasticity measurement based on the HHT marginal Spectrum entropy of sEMG using a portable system: a preliminary study. IEEE Trans Neur Sys Rehabil Eng. 2018;26(7):1424–34.

    Article  Google Scholar 

  8. Lei M, Wang Z, Feng Z. Detecting nonlinearity of action surface EMG signal. Phys Lett A. 2001;290(5–6):297–303.

    Article  Google Scholar 

  9. Wentink EC, Beijen SI, Hermens HJ, Rietman JS, Veltink PH. Intention detection of gait initiation using EMG and kinematic data. Gait Posture. 2013;37(2):223–8.

    Article  Google Scholar 

  10. Vannozzi G, Conforto S, D’Alessio T. Automatic detection of surface EMG activation timing using a wavelet transform based method. J Electromyogr Kinesiol. 2010;20(4):767–72.

    Article  Google Scholar 

  11. Dalley SA, Varol HA, Goldfarb M. An improved approach for measuring the tonic stretch reflex response of spastic muscles. IEEE Trans Neur Sys Rehabil Eng. 2012;20(1):58–67.

    Article  Google Scholar 

  12. Silva MB, Silva AN, Naves ELM, Palomari ET, Soares AB. An improved approach for measuring the tonic stretch reflex response of spastic muscles. Comput Biol Med. 2017;80(1):166–74.

    Article  Google Scholar 

  13. Calota A, Feldman AG, Levin MF. Spasticity measurement based on tonic stretch reflex threshold in stroke using a portable device. Clin Neurophysiol. 2008;119(10):2329–37.

    Article  Google Scholar 

  14. Kim KS, Seo JH, Song CG. Portable measurement system for the objective evaluation of the spasticity of hemiplegic patients based on the tonic stretch reflex threshold. Med Eng Phys. 2011;33(1):62–9.

    Article  Google Scholar 

  15. Young RR. Spasticity: a review. Neurology. 1994;44(11):S12–20.

    Google Scholar 

  16. Lee AS, Cholewicki J, Reeves NP. The effect of background muscle activity on computerized detection of sEMG onset and offset. J Biomech. 2007;40(15):3521–6.

    Article  Google Scholar 

  17. Hodges PW, Bui BH. A comparison of computer-based methods for the determination of onset of muscle contraction using electromyography. Electroencephalogr Clin Neurophysiol. 1996;101(6):511–9.

    Google Scholar 

  18. Solnik S, Rider P, Steinweg K, DeVita P, Hortobagyi T. Teager-Kaiser energy operator signal conditioning improves EMG onset detection. Eur J Appl Physiol. 2010;110(3):489–98.

    Article  Google Scholar 

  19. Xu Q, Quan YZ, Yang L, He J. An adaptive algorithm for the determination of the onset and offset of muscle contraction by EMG signal processing. IEEE Trans Neur Sys Rehabil Eng. 2013;21(1):65–73.

    Article  Google Scholar 

  20. Yang DP, Zhang HJ, Gu YK, Liu H. Accurate EMG onset detection in pathological, weak and noisy myoelectric signals. Biom Sig Proce Cont. 2017;33:306–15.

    Article  Google Scholar 

  21. Richman JS, Moorman JR. Physiological time-series analysis using approximate entropy and sample entropy. Am J Physiol Heart Circ Physiol. 2000;278(6):H2039–49.

    Article  Google Scholar 

  22. Zhang X, Zhou P. Sample entropy analysis of surface EMG for improved muscle activity onset detection against spurious background spikes. J Electromyogr Kinesiol. 2012;22(6):901–7.

    Article  Google Scholar 

  23. Xie HB, Guo JY, Zheng YP. Fuzzy approximate entropy analysis of chaotic and natural complex systems: detecting muscle fatigue using electromyography signals. Ann Biomed Eng. 2010;38(4):1483–96.

    Article  Google Scholar 

  24. Xie HB, He WX, Liu H. Measuring time series regularity using nonlinear similarity-based sample entropy. Phys Lett A. 2008;372(48):7140–6.

    Article  Google Scholar 

  25. Chen WT, Wang ZZ, Xie HB, Yu W. Characterization of surface EMG signal based on fuzzy entropy. IEEE Trans Neur Sys Rehabil Eng. 2007;15(2):266–72.

    Article  Google Scholar 

  26. Xie HB, Chen WT, He WX, Liu H. Complexity analysis of the biomedical signal using fuzzy entropy measurement. Appl Soft Comput. 2011;11(2):2871–9.

    Article  Google Scholar 

  27. Chen WT, Zhuang J, Yu WX, Wang ZZ. Measuring complexity using FuzzyEn, ApEn, and SampEn. Med Eng Phys. 2009;31(1):61–8.

    Article  Google Scholar 

  28. Molina-picó A, Cuesta-frau D, Aboy M, Crespo C, Pau MM, Sandra OC. Comparative study of approximate entropy and sample entropy robustness to spikes. Artif Intell Med. 2011;53(2):97–106.

    Article  Google Scholar 

  29. Islam MK, Rastegarnia R, Yang Z. Methods for artifact detection and removal from scalp EEG: a review. Neurophysiol Cliniq Clin Neurophysiol. 2016;46(4–5):287–305.

    Article  Google Scholar 

  30. Andrade AO, Nasuto S, Kyberd P, Sweeney-Reed CM, Van Kanijn FR. EMG signal filtering based on empirical mode decomposition. Biom Sigl Proce Cont. 2006;1(1):44–55.

    Article  Google Scholar 

  31. Huang NE, Shen Z, Long SR, Wu MC, Shih HH, Zheng Q, Yen NC, Tung CC, Liu HH. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc R Soc Lond. 1998;454:903–95.

    Article  MathSciNet  Google Scholar 

  32. Zhang X, Zhou P. Filtering of surface EMG using ensemble empirical mode decomposition. Med Eng Phys. 2013;35(4):537–42.

    Article  Google Scholar 

  33. Hu BH, Zhang XF, Mu JS, Wu M, Wang Y. Spasticity assessment based on the Hilbert–Huang transform marginal spectrum entropy and the root mean square of surface electromyography signals: a preliminary study. Biomed Eng Online. 2018;17:27.

    Article  Google Scholar 

  34. Malone A, Meldrum D, Gleeson J, Bolger C. Reliability of surface electromyography timing parameters in gait in cervical spondylotic myelopathy. J Electromyogr Kinesiol. 2011;21(6):1004–10.

    Article  Google Scholar 

  35. Silva LEV, Murta LOJ. Evaluation of physiologic complexity in time series using generalized sample entropy and surrogate data analysis. Chao Interdiscipl Jour Nonl Scie. 2012;22(4):043105.

    Article  MathSciNet  Google Scholar 

  36. Lake DE, Richman JS, Griffin MP, Moorman JR. Sample entropy analysis of neonatal heart rate variability. Am J Phys Regul Integr Comp Phys. 2002;283(3):R789–97.

    Google Scholar 

  37. Ramdani S, Bouchara F, Lagarde J. Influence of noise on the sample entropy algorithm. Chao Interdiscipl Jour Nonl Scie. 2009;19(1):013123.

    Article  MathSciNet  Google Scholar 

  38. Marwaha P, Sunkaria RK. Complexity quantification of cardiac variability time series using improved sample entropy (I-SampEn). Austral Physi Eng Scien Med. 2016;39(3):755–63.

    Article  Google Scholar 

  39. Alcaraz R, Rieta JJ. A non-invasive method to predict electrical cardioversion outcome of persistent atrial fibrillation. Med Biol Eng Comput. 2008;46(7):625–35.

    Article  Google Scholar 

Download references

Acknowledgements

The authors are very grateful to the participants and staff who participated in this project.

Funding

This work was supported by the National Natural Science Foundations of China (No. U1713210).

Author information

Authors and Affiliations

Authors

Contributions

MD and BH implemented the algorithm, analyzed the data and interpreted the results. MD wrote the first draft of the manuscript. MD, FX and YW designed the study and performed all stages of the study including data collection, analysis, interpretation, and substantial revision of the manuscript. MW and ZZ gave valuable comments about rehabilitation medicine. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yong Wang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Anhui Provincial Hospital (Approval Number: 2015-Ethical Approval-08) and all the participants gave their informed consent.

Consent for publication

All participants in the study verbally agreed to publish personal data.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Du, M., Hu, B., Xiao, F. et al. Detection of stretch reflex onset based on empirical mode decomposition and modified sample entropy. BMC biomed eng 1, 23 (2019). https://doi.org/10.1186/s42490-019-0023-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42490-019-0023-y

Keywords