ISSN ONLINE(2320-9801) PRINT (2320-9798)

All submissions of the EM system will be redirected to Online Manuscript Submission System. Authors are requested to submit articles directly to Online Manuscript Submission System of respective journal.

Ocular Artifacts Reduction in EEG Using DWT And ANC For Portable Applications

R.Hemakumar, G.Banumathy
PG Student, M.E-Communication Systems, Valliammai Engineering College, Chennai, Tamilnadu, India
Related article at Pubmed, Scholar Google

Visit for more related articles at International Journal of Innovative Research in Computer and Communication Engineering

Abstract

A new model to remove ocular artifacts (OA) from electroencephalograms (EEGs) is presented. The model is based on discrete wavelet transformation (DWT) and adaptive noise cancellation (ANC). Using simulated and measured data, the accuracy of the model is compared with the accuracy of other existing methods based on stationary wavelet transforms and our previous work based on wavelet packet transform and independent component analysis. A particularly novel feature of the new model is the use of DWTs to construct an OA reference signal, using the three lowest frequency wavelet coefficients of the EEGs. The results show that the new model demonstrates an improved performance with respect to the recovery of true EEG signals and also has a better tracking performance. Because the new model requires only single channel sources, it is well suited for use in portable environments.

Keywords

Adaptive noise cancellation (ANC), electroencephalogram (EEG), ocular artifacts (OAs), signal processing.

INTRODUCTION

Every day we are exposed to several solicitations for purchasing products, voting or supporting particular politicians and even improving our life style. Such pressure has become usual, being mediated by the entire current media available, video, audio, and even internet. How and to what extent these messages could be detected and recognized by our brain is still not well understood. In fact, the study of brain responses to commercial and political announcements has been measured mainly by thermodynamic responses of the different brain areas, by using the functional Magnetic Resonance Imaging devices. However, both the stimuli and the relative brain responses have rapidly shifting characteristics that are not tracked by the evolution of the hemodynamic blood flow, which usually lasts 4–6 seconds. Different brain imaging tools, mainly EEG[1,8] and Magneto encephalography, exhibit a sufficient time resolution to follow the brain activity at an expense of a coarse level of spatial resolution with respect to the MRI. In fact, during those last ten years, the use of the high resolution EEG techniques has retrieved an increased amount of information related to the brain during activities related to complex cognitive tasks, such as memory, visual attention, short-term memory, and so forth.
Starting from the interesting characteristic of the high resolution EEG techniques for the tracking of brain activity, the present work would like to describe neuro electric based methodology for the assessment of the efficacy of commercial, politic, and Public Service Announcements (PSAs). The aim of the brain imaging techniques applied to the fruition of commercial advertising is to understand mechanisms underlying customer’s engagement with brand or company advertised. In particular, the issue is to explain how the exposure of subsequent film segments is able to trigger in the consumer mind persisting stimuli leading to interest, preference, purchase, and repurchase of a given product.

WORKING OF PROJECT

COMPONENT ANALYSIS

In recent years, many traditional approaches have been proposed to remove or attenuate such OAs from recorded EEG. Widely used methods for attenuating OAs are mostly based on time domain or frequency domain techniques. Principal component analysis (PCA) has also been used to remove artifacts from EEG[2] . However, PCA cannot remove the artifacts from EEG completely, because sometimes the waveforms of the OAs are smaller in amplitude with respect to the ongoing EEG. Later, independent component analysis (ICA) was proposed. ICA is developed with respect to blind source separation, where the aim is to obtain components that are approximately independent . However, to remove OAs, ICA needs a reference signal that requires tedious visual classification of the components .
Some research methods have used the modelling of OAs components based on improved support vector machine techniques to isolate them from the EEG. Discrete wavelet transform (DWT) is a method that neither relies upon the reference OAs nor visual inspection. As frequency precision is also improved, researchers have often used DWT to detect the OAs region and then selected the correct threshold to remove the interference.

ARTIFACTS

During the measurement EEG signals, other kinds of signals, such as ECG, EMG, 50Hzelectric interference, so called artifacts, are also captured by electrodes. These signals mixwith EEG signal and EEG signal are therefore contaminated. Artifacts are undesired signals that can introduce significant changes in neurological signalsand ultimately affect the neurological phenomenon . Some types of artifacts will increase or decrease alpha power which leads to mistakes in the alpha wave measurement. In the application of relaxation state detection, we will have wrong detections caused by artifacts.

Characteristics of Artifacts

There are many types of artifacts. Some of them will not influence the frequency band ofalpha waves. For instance, we see a sharp peak around 50 Hz.This is not a feature of EEG signals; however, it is caused by the artifacts of 50Hz electric interference. Luckily this peak is located far from the alpha frequency bands. Therefore,it does not influence the power change of the alpha waves. Artifacts that have relatively considerable influence on the alpha wave can be categorized as follows.

Ocular artifacts (OA)

Ocular artifacts are caused by eye movements, such as blinking eyes and rolling eyeballs. It shows a segment of EEG signals contaminated by eyes blinking. Here we showsa segment of EEG [3] signals contaminated by eyeballs rolling. In the frequency domain, ocular artifacts increase the power of EEG signals from 2Hz to 20 Hz. Unluckily,alpha waves locate between 8Hz and 12Hz. Hence, the presence of ocular artifacts in EEG signals will cause unreliable detection of alpha power.

COMBINATION OF DWT AND ANC

We applied a model based on WPTand ICA to remove the OAs in a contaminated EEG signal. Although the model employed existing methods to compensatefor ICA’s uncertainties, the tracking performance of signals wasnot good and compensation for the uncertainties [10,11] of the ICAalgorithm had to be estimated, producing a distorted signal after denoising. Therefore, in this paper, we have applied a new modelbased on DWT and adaptive noise cancellation (ANC) to removethe OAs. The first and most important step of our new methodis the construction of a reference signal using DWT. With thisreferencesignal, a new model was established based on ANC; hence, a combination of DWT and ANC is employed. ANC isoften used to remove power line interference in EEG signals.
However, it does require a reference signal that has proveddifficult to obtain in biological systems. In this paper, we applyDWT to the contaminated EEG to derive a reference signalthat is used in a further stage of ANC processing. This allowsus to lever the biggest advantage of ANC: that it can followthe changes and automatically adjust its parameters to achieveoptimal performance of the filter when the statistical propertiesof the input signal are changing. This is a novel approach thatwe demonstrate to be effective even when the EEG signal hasonly one channel and is hence particularly suitable for portableapplications.

METHODOLOGY

The recorded EEG signals are contaminated by OAs, thiscontamination is considered to be an additive noise within theEEG signal. So, we can write the following expression
equation

Operations Involved

The first is the construction of the reference signal; the second isthe removal the OAs from the recorded EEG signal by applyingANC, based on a recursive least-squares (RLS) algorithm.OAs is mainly concentrated in the low frequency band, soDWT is used to construct the OAs in the frequency domain.DWT is a multi resolution representation of signals and images.It can be used to decompose signals into multiscale representations .
It is a commonly used tool for analyzing non stationary signals. The wavelets used in DWT are effectivein constructing both time- and frequency-domain information from time-varying and non stable EEG signals. There has been much research on the use of DWT to remove the artifacts in EEG signals. Researchers combine ICA and DWT to remove artifacts from EEG signals, employing the DWT to improve the contrast of the output after ICA. In DWT is used to detect the OA region and thenselect the correct threshold to remove interference.

INTERFERENCE CANCELLATION

As for the mathematical notation used throughout this section, all quantities are assumed to be real-valued. Scalar and vector quantities shall be indicated by lowercase (e.g., x) [4,5] and uppercase-bold (e.g., X) letters, respectively. We represent scalar and vector sequences or signals as x(n) and X(n), respectively, where n denotes the discrete time or discrete spatial index, depending on the application.
Herew represents the coefficients of the FIR filter tap weight vector, x(n) is the input vector samples, z -1 is a delay of one sample periods, y(n) is the adaptive filter output, d(n) is the desired echoed signal and e(n)is the estimation error.

Adaptive filter

The aim of an adaptive filter is to calculate the difference between the desired signal and the adaptive filter output, e(n).[6,9] This error signal is fed back into the adaptive filter and its coefficients are changed algorithmically in order to minimize a function of this difference, known as the cost function. In the case of acoustic echo cancellation, the optimal output of the adaptive filter is equal in value to the unwanted echoed signal. When the adaptive filter output is equal to desired signal the error signal goes to zero. In this situation the echoed signal would be completely cancelled and the far user would not hear any of their original speech returned to them.
This section examines adaptive filters and various algorithms utilized. The various methods used in this thesis can be dividedinto two groups based on their cost functions. The first class are known as Mean Square Error (MSE) adaptive filters, they aim to minimize a cost function equal to the expectation of the square of the difference between the desired signal d(n), and the actual output of the adaptive filter y(n) .
equation
Where k=1, 2, 3…n., k=1 corresponds to the time at which the RLS algorithm commences. Later we will see that in practice no tall previous values are considered; rather only the previous N (corresponding to the filter order) error signals are considered.
equation

RESULTS

In this paper, we have introduced a new model that combines DWT and ANC techniques to remove the OAs in contaminated EEG signals. We have demonstrated the effectiveness of our new model by using the model to process simulated and standard EEG data. Our new model is able to eliminate OAs in the low-frequency band even when their frequency is overlap-ping with that of the EEG signal.Fig 3 shows the original EEG signal. Fig 4 shows the EEG signal with artifacts. After wavelet decomposition and adaptive noise cancellation the coefficient table is updated based on the noise present in the EEG signal its shown in Fig 5. Finally the denoised EEG signal is extracted plotted in Fig 6. SNR value for the noisy EEG signal is 1.5935 and SNR for denoised EEG signal is 1.6265 is achieved.

DISCUSSION

In fig EEG signal was generated and the ocular artifacts were added to the EEG signal. Depends upon the noise or ocular artifacts the coefficient updation graph is obtained. Finally the noiseless EEG signal was extracted through MATLAB.Finally the SNR values for respected output is obtained. SNR value changes based on the noise level.

Figures at a glance

Figure 1 Figure 2 Figure 3
Figure 1 Figure 2 Figure 3
Figure 4 Figure 5 Figure 6
Figure 4 Figure 5 Figure 6
 

References