Feature Extraction of Eeg Signals Using Power Spectral Entropy
Feature Extraction of Eeg Signals Using Power Spectral Entropy
Authorized licensed use limited to: UNIVERSITY PUTRA MALAYSIA. Downloaded on April 21,2010 at 00:13:55 UTC from IEEE Xplore. Restrictions apply.
session. The subject sat in a relaxing chair with
armrest. The task was to control a feedback bar by
means of imaginary left or right hand movement. The
experimental is described in detail in [10, 11]. The aim
of signal processing was to distinguish the mental tasks 1 2 3
5 cm
experiment consists of 7 runs with 40 trials each. All C3 1 Cz 2 C4 3
Power spectral entropy is information entropy that is 1) The Discrete Fourier Transform X (ωi ) of signal
able to quantify the spectral complexity of an uncertain
system. It is defined in [12]. For an uncertain system, can be obtained by FFT; where ωi is the frequency
let’s assume a random variable X as states of the point of the number i .
system, the value of X is 2) Calculate its power spectral density (PSD)
1
X = {x1 , x2 ,", xn } Pˆ (ωi ) =
2
(n ≥ 1) (1) X (ωi )
N
The corresponding probability is ˆ
3) Normalize P (ωi ) , and obtain power spectral
P = {p1, p2 ," , pn } 0 ≤ pi ≤ 1 , i = 1,2," , n (2) density distribution function
Under constraints Pˆ (ωi )
pi =
n
∑ Pˆ (ω )
∑ pi = 1
i
(3) i
i =1 4) Use equation (4), the PSE can be worked out.
Therefore, the information entropy of the system can PSE can be interpreted as a measurement of the time
be expressed as uncertainty in frequency domain. In power spectrum of
n EEG signals, when the spectrum peak is narrow, its
H = −∑ pi ln pi (4) entropy value is small. It indicates that the signal has
i =1 an obvious concussive rhythm, that is to say, if wave is
The time-series signals become power spectrum by orderliness, its complexity is small; whereas, the
FFT transform, and the information entropy of power spectrum peak is more smoothness, and its entropy
spectrum is called power spectral entropy. The value is more greatness. Therefore, the PSE can reflect
algorithm can be summarized as: the spectra structure of EEG signals.
1360
436
Authorized licensed use limited to: UNIVERSITY PUTRA MALAYSIA. Downloaded on April 21,2010 at 00:13:55 UTC from IEEE Xplore. Restrictions apply.
2.3. Classification
1.25
w = Sw
−1
(µ1 − µ 2 ) (7) 1.16
y = w x − 0.5( w, µ1 + w, µ 2
T
) (8) 1.12
1.1
To improve classification accuracy of time t0 , a
1.08
time-variable linear classifier is proposed in this paper 1.06
0.75
0.65
EEG data set which was used in the BCI competition 0.55
1361
437
Authorized licensed use limited to: UNIVERSITY PUTRA MALAYSIA. Downloaded on April 21,2010 at 00:13:55 UTC from IEEE Xplore. Restrictions apply.
generally a circumscribed “event-related desynchron- Science Foundation of China (No. 30670529) and the
ization” (ERD) is characteristic for the contralateral, Chunhui Plan of Ministry of Education of China
and an “event-related synchronization” (ERS) for the (Z2005-2-62007)
ipsilateral hemisphere [13]. This fact is a foundation for
distinguish left and right hand movements. References
Fig.2 displays the PSE on channels C3 and C4
corresponding to imagine left hand movement. It [1] R. Herbert, M. G. Johannes, and G. Pfurtscheller,
shows difference in PSE distribution between the two “Optimal Spatial Filtering of Signal Trial EEG during
tasks. For channel C3 [Fig.2 (a)], the values of PSE by Imagined Hand Movement”, IEEE Trans .Rehab. Eng. , Dec.
imaginary left hand is lower than that of right hand. 2000, vol. 8(4), pp. 441-446.
Similar conclusion can be drawn from Fig.2 (b) for
channel C4. The time change regulation of PSE is [2] B. O. Peters, G. Pfurscheller, and H. Flyvjberg,
“Automatic differentiation of multi-channel EEG signals” ,
accordance with ERD/ERS. And the PSE has a good IEEE Trans. Biomed. Eng. , 2001, vol. 48(1), pp.111-116.
separability for the two tasks. So the PSE on channels
C3 and C4 can be used as feature of EEG signals [3] D. Coyle, G. Prasad, and T. M. McGinnity, “A time-
during imagined right and left hand movements. series prediction approach for feature extraction in a Brain-
According to the analysis above, the classification Computer Interface”, IEEE Trans. Neur. Syst. Rehabil. Eng. ,
accuracy can be obtained with the feature of PSE and Dec. 2005, vol.13, pp. 461-467.
time-variable linear classifier. For training data set, the
feature vectors were extracted at every sampling point, [4] V. Bostanov, “BCI Competition 2003-Data sets Ib and
with the sliding window size being 128 samples and IIb: Feature extraction from Event-Related brain potentials
with the continuous wavelet transform and t-Value
the interval is one sample. The time-variable linear scalogram”, IEEE Trans. Biomed. Eng. , June 2004, vol.51,
classifier was trained with these feature vectors. Then pp. 1057-1061.
the feature vectors of testing data set was calculated by
similar method and classified by trained classifier. [5] S. M. Pincus, and A. L. Goldberger, “Physiological time-
From Fig.2, it shows obviously difference in PSE of series analysis: what does regularity quantify?”, American
imaginary left and right hand movements during 4.5s Journal of Physiology,1994, 1643-1656.
to 7s. So the accumulation of information was started
from 4.5s in the algorithm of classification. The time [6] H. Z. Qi, and L. L. Cheng, “Study on Information
course of classification accuracy is shown in Fig.3. It Entropy of EEG Signals during Motor Imaginary”, Chinese
Journal Biomedical Engineering, Feb. 2007, vol.26 (1), pp.
shows the maximal accuracy achieves 90% around 74-77.
t=6.7s. Subject was quiet at first 2s, and at t=3s
beginning to imagine the mental tasks. Because of the [7] H. J. Wu, H. Zhang, and C. X. Cheng, “Application of
accumulation of information, the accuracy achieves power spectral entropy to the noninvasive detection of focal
max around t=6.7s. ischemic cerebral injury”, Journal of Biomedical Engineer-
ing, 2003, vol. 20 (2), pp.229-232.
4. Conclusion
[8] H. Liu, W. X. He, and X. P. Chen, “EEG time-series
analysis using nonlinear dynamics for sleep monitoring”,
This paper proposes new feature extraction and Journal of Jiangsu University: Natural Science, 2005, vol.26
classification methods for EEG signals during imaged (2), pp.174-177.
left and right hand movements. The Experimental
results have shown that the PSE can be used as one of [9] D. Z. Zhang, “Relationship between age and power
the feature of mental tasks successfully. Furthermore, spectral entropy of blood pressure waveform”, Chinese
we applied time-variable classifier for classify and Journal of Medical Physics, 2005, vol.22 (4), pp.606-607.
gained a good result. The result demonstrates that the
proposed algorithm is reasonable and effective for [10] A. Schlögl, K. Lugger, and G.Pfurtscheller, “Using
classification of imaginary left and right hand Adaptive Autoregressive Parameters for a Brain-Computer
movements. High classification accuracy and quick Interface Experiment”, Proceedings of the 19th Annual
calculation speed of the algorithm make it becoming a International Conference IEEE/EMBS, 1997, pp.1533-1535.
feasible choice for a practical BCI system.
[11] G. Pfurtscheller, C. Neuper, A. Schlögl, et al, “Separabi-
lity of EEG signals recorded during right and left motor
Acknowledgment imagery using adaptive autoregressive parameters”, IEEE
Trans. Rehab. Eng ., 1998, vol. 6(3), pp.316-25.
This work was supported by the National Natural
1362
438
Authorized licensed use limited to: UNIVERSITY PUTRA MALAYSIA. Downloaded on April 21,2010 at 00:13:55 UTC from IEEE Xplore. Restrictions apply.
[12] J. N. Kapur, and H.K. Kesavan, Entropy optimization
principles with applications, Academic Press, San Diego,
1992.
1363
439
Authorized licensed use limited to: UNIVERSITY PUTRA MALAYSIA. Downloaded on April 21,2010 at 00:13:55 UTC from IEEE Xplore. Restrictions apply.