EEG Data processing in GNU Octave

Resting EEG of control subjects was recorded in a quiet room with eyes open (5 min) and eyes closed (5 min)* using a 64-channel cap. The EEG data are in various formats, such as European Data Format (.edf), EEGLAB (.fdt, .set) or Biotrace Compressed Data (.bcd), and do not contain events.

The purpose of the programme is to perform the following tasks in GNU Octave or Matlab, which can run under GNU Octave:

1. transport the EEG data to GNU Octave
2. remove noise and other artefacts
3. provide the possibility to select certain positions in the extended 10-20 system from the EEG data (e.g. Cz, F3, F4)
4. apply Hanning windowing to the filtered signal at the selected position
5. use an elliptic bandpass filter [1.5-45] GHz (see attachment)
6. call the MATLAB elliptic filter function to extract theta, alpha, beta, beta-low, beta-high and gamma bands
7. plot the individual frequency bands
8. create an FFT and calculate and plot the amplitude, absolute power, relative power and power spectral density (PSD) of theta, alpha, beta-low, beta-high, gamma, alpha/beta, alpha/theta, theta/beta, alpha/(beta+gamma), (theta/beta+gamma).

*Not all data have this temporal structure, so it would be desirable for the time to be entered in some way, either automatically or manually.