In this work, the problem of detecting neuronal spikes or action potentials (AP) in noisy recordings from a Microelectrode Array (MEA) is investigated. In particular, the spike detection algorithms should be less complex and with low computational complexity so as to be amenable for real time applications. The use of the MEA is that it allows collection of extracellular signals from either a single unit or multiple (45) units within a small area. The noisy MEA recordings then undergo basic filtering, digitization and are presented to a computer for further processing. The challenge lies in using this data for detection of spikes from neuronal firings and extracting spatiotemporal patterns from the spike train which may allow control of a robotic limb or other neuroprosthetic device directly from the brain. The aim is to understand the spiking action of the neurons, and use this knowledge to devise efficient algorithms for Brain Machine Interfaces (BMIs). An effective BMI will require a realtime, computationally efficient implementation which can be carried out on a DSP board or FPGA system. The aim is to devise algorithms which can detect spikes and underlying spatio-temporal correlations having computational and time complexities to make a real time implementation feasible on a specialized DSP chip or an FPGA device. The time-frequency localization, multiresolution representation and analysis properties of wavelets make them suitable for analysing sharp transients and spikes in signals and distinguish them from noise resembling a transient or the spike. Three algorithms for the detection of spikes in low SNR MEA neuronal recordings are proposed: 1. A wavelet denoising method based on the Discrete Wavelet Transform (DWT) to suppress the noise power in the MEA signal or improve the SNR followed by standard thresholding techniques to detect the spikes from the denoised signal. 2. Directly thresholding the coefficients of the Stationary (Undecimated) Wavelet Transform (SWT) to detect the spikes. 3. Thresholding the output of a Teager Energy Operator (TEO) applied to the signal on the discrete wavelet decomposed signal resulting in a multiresolution TEO framework. The performance of the proposed three wavelet based algorithms in terms of the accuracy of spike detection, percentage of false positives and the computational complexity for different types of wavelet families in the presence of colored AR(5) (autoregressive model with order 5) and additive white Gaussian noise (AWGN) is evaluated. The performance is further evaluated for the wavelet family chosen under different levels of SNR in the presence of the colored AR(5) and AWGN noise. Chapter 1 gives an introduction to the concept behind Brain Machine Interfaces (BMIs), an overview of their history, the current state-of-the-art and the trends for the future. It also describes the working of the Microelectrode Arrays (MEAs). The generation of a spike in a neuron, the proposed mechanism behind it and its modeling as an electrical circuit based on the Hodgkin-Huxley model is described. An overview of some of the algorithms that have been suggested for spike detection purposes whether in MEA recordings or Electroencephalographic (EEG) signals is given. Chapter 2 describes in brief the underlying ideas that lead us to the Wavelet Transform paradigm. An introduction to the Fourier Transform, the Short Time Fourier Transform (STFT) and the Time-Frequency Uncertainty Principle is provided. This is followed by a brief description of the Continuous Wavelet Transform and the Multiresolution Analysis (MRA) property of wavelets. The Discrete Wavelet Transform (DWT) and its filter bank implementation are described next. It is proposed to apply the wavelet denoising algorithm pioneered by Donoho, to first denoise the MEA recordings followed by standard thresholding technique for spike detection. Chapter 3 deals with the use of the Stationary or Undecimated Wavelet Transform (SWT) for spike detection. It brings out the differences between the DWT and the SWT. A brief discussion of the analysis of non-stationary time series using the SWT is presented. An algorithm for spike detection based on directly thresholding the SWT coefficients without any need for reconstructing the denoised signal followed by thresholding technique as in the first method is presented. In chapter 4 a spike detection method based on multiresolution Teager Energy Operator is discussed. The Teager Energy Operator (TEO) picks up localized spikes in signal energy and thus is directly used for spike detection in many applications including R wave detection in ECG and various (alpha, beta) rhythms in EEG. Some basic properties of the TEO are discussed followed by the need for a multiresolution approach to TEO and the methods existing in literature. The wavelet decomposition and the subsampled signal involved at each level naturally lends it to a multiresolution TEO framework at the same time significantly reducing the computational complexity due the subsampled signal at each level. A wavelet-TEO algorithm for spike detection with similar accuracies as the previous two algorithms is proposed. The method proposed here differs significantly from that in literature since wavelets are used instead of time domain processing. Chapter 5 describes the method of evaluation of the three algorithms proposed in the previous chapters. The spike templates are obtained from MEA recordings, resampled and normalized for use in spike trains simulated as Poisson processes. The noise is modeled as colored autoregressive (AR) of order 5, i.e AR(5), as well as Additive White Gaussian Noise (AWGN). The noise in most human and animal MEA recordings conforms to the autoregressive model with orders of around 5. The AWGN Noise model is used in most spike detection methods in the literature. The performance of the proposed three wavelet based algorithms is measured in terms of the accuracy of spike detection, percentage of false positives and the computational complexity for different types of wavelet families. The optimal wavelet for this purpose is then chosen from the wavelet family which gives the best results. Also, optimal levels of decomposition and threshold factors are chosen while maintaining a balance between accuracy and false positives. The algorithms are then tested for performance under different levels of SNR with the noise modeled as AR(5) or AWGN. The proposed wavelet based algorithms exhibit a detection accuracy of approximately 90% at a low SNR of 2.35 dB with the false positives below 5%. This constitutes a significant improvement over the results in existing literature which claim an accuracy of 80% with false positives of nearly 10%. As the SNR increases, the detection accuracy increases to close to 100% and the false alarm rate falls to 0. Chapter 6 summarizes the work. A comparison is made between the three proposed algorithms in terms of detection accuracy and false positives. Directions in which future work may be carried out are suggested.

No description available.

Work is aimed at the time-frequency decomposition of a signal application for monitoring the EKG trend progression. Goal is to create algorithm which would watch changes in the ST segment in EKG recording and its realization in the Matlab program. Analyzed is substance of the origin of EKG and its measuring. For trend calculations after reading the signal is necessary to preprocess the signal, it consists of filtration and detection of necessary points of EKG signal. For taking apart, also filtration and measuring the signal is used wavelet transformation. Source of the data is biomedicine database Physionet. As an outcome of the algorithm are drawn ST segment trends for three recordings from three different patients and its comparison with reference method of ST qualification. For qualification of the heart stability, as a system, where designed methods watching differences in position of the maximal value in two-zone spectrum and the Poincare mapping method. Realized method is attached to this thesis.

This thesis deals with a proposition of an algorithm for QRS complex and typical ECG waves boundaries detection. It incorporates a literature research focused on heart electrophysiology and commonly used methods for ECG fiducial points detection and delineation. Out of the presented methods an algorithm based on a continuous wavelet transform is implemented. Detection and delineation algorithm is tested on CSE standard signal database towards references determined both manually and automatically. Obtained results are compared to other congenerous methods. The diploma thesis is further concerned with an algorithm modification for experimental electrocardiograms of isolated rabbit hearts. Recording specifics of these data are introduced. Additionally, based on time and frequency analysis, particular modifications of the algorithm are proposed and realized. Due to a large extent of records functionality is verified on randomly selected database samples. Efficiency of the modified algorithm is evaluated through manually annotated references.

Signal compression is daily-used tool for memory capacities reduction and for fast data communication. Methods based on wavelet transform seem to be very effective nowadays. Signal decomposition with a suitable bank filters following with coefficients quantization represents one of the available technique. After packing quantized coefficients into one sequence, run length coding together with Huffman coding are implemented. This thesis focuses on compression effectiveness for the different wavelet transform and quantization settings.

This thesis deals with time-frequency analysis of electrograms measured on isolated guinea pig hearts perfused according to Langendorff. Time-frequency analysis is based on algorithms Matching Pursuit and Wigner-Ville Distribution. The theoretical part describes the basics of electrocardiography, measurement on isolated hearts, the theory of approximation method Matching Pursuit and its combination with the Wigner-Ville distribution spectrum showing the energy density of the signal. Also other common approaches of time-frequency analysis are presented including the theory of continuous wavelet transform. The presented algorithms were tested on a set of electrograms, on which were induced ischemia within measurement followed by reperfusion. The proposed method allows for the fast detection of ischemia without any a priori knowledge of the signal, and also serves as a tool for measurement of EG important points and intervals. In the conclusion efficacy of the method was presented and its possible uses has been discussed.

Epilepsy affects more than 50 million people and is one of the most prevalent neurological disorders and has a high impact on the quality of life of those suffering from it. However, 70% of epilepsy patients can live seizure free with proper diagnosis and treatment. Patients are evaluated using scalp EEG recordings which is cheap and non-invasive. Diagnostic yield is however low and qualified personnel need to process large amounts of data in order to accurately assess patients. MindReader is an unsupervised classifier which detects spectral anomalies and generates a hypothesis of the underlying patient state over time. The aim is to highlight abnormal, potentially epileptiform states, which could expedite analysis of patients and let qualified personnel attest the results. It was used to evaluate 95 scalp EEG recordings from healthy adults and adult patients with epilepsy. Interictal Epileptiform discharges (IED) occurring in the samples had been retroactively annotated, along with the patient state and maneuvers performed by personnel, to enable characterization of the classifier’s detection performance. The performance was slightly worse than previous benchmarks on pediatric scalp EEG recordings, with a 7% and 33% drop in specificity and sensitivity, respectively. Electrode positioning and partial spatial extent of events saw notable impact on performance. However, no correlation between annotated disturbances and reduction in performance could be found. Additional explorative analysis was performed on serialized intermediate data to evaluate the analysis design. Hyperparameters and electrode montage options were exposed to optimize for the average Mathew’s correlation coefficient (MCC) per electrode per patient, on a subset of the patients with epilepsy. An increased window length and lowered amount of training along with an common average montage proved most successful. The Euclidean distance of cumulative spectra (ECS), a metric suitable for spectral analysis, and homologous L2 and L1 loss function were implemented, of which the ECS further improved the average performance for all samples. Four additional analyses, featuring new time-frequency transforms and multichannel convolutional autoencoders were evaluated and an analysis using the continuous wavelet transform (CWT) and a convolutional autoencoder (CNN) performed the best, with an average MCC score of 0.19 and 56.9% sensitivity with approximately 13.9 false positives per minute.

The current research work aims at the development of Biomedical Neuroengineering tools (Biotechnologies) for the in-depth functional study, rapid diagnosis, continuous monitoring and well-timed management of acute and chronic brain disorders, of individuals that are subjected to or suffer from any kind of systemic hypoxaemia or more localized brain hypoxia; as well as the functional assessment and continuous control of adaptability during the training of “altinauts” and generally of individuals that practice activities and function within environments of increased visual-cognitive-motor response demands (a type of brain “stress test”). For this purpose, we subject the entire visuocognitive system, from the elementary sensory to the most complex cognitive level, to an experimental test of categorical discrimination of complex visuocognitive stimuli, following ultra-rapid visual stimulation that leads to a motor response upon categorization of targets (images of animals elicit productive responses) and to its suppression upon categorization of nontargets (images of nonanimals elicit inhibitory responses). The oscillatory electro-physiological responses that are concurrently recorded at the occipital-temporal-parietal brain areas are analyzed in the time-domain (<20Hz) and in the joint time-frequency domain broadband (1-60Hz) with the Continuous Wavelet Transform that optimizes the multiresolution analysis of the high frequency (≥20Hz) γ-band oscillatory activity. This visuocognitive categorization test takes place in normoxaemic as well as hypoxaemic conditions (monitored reduction in the blood oxygen saturation from ≥97% to around 80% under conditions of hypobaric hypoxia within a hypobaric chamber), in order to assess electrophysiological markers that can detect and capture in the most sensitive and dynamic way even so transient, short-living and rather mild changes in brain function. The statistical parametric analysis of the time-frequency maps and the generalized, statistically safer, method of analysis of variance have established as the most sensitive and reliable the following markers: the major deflections of the evoked potentials, the phase-coherence factor of the oscillations across single-trials and the elicited energy of the evoked/phase-locked and the induced/total oscillatory activity. These electrophysiological markers in conjunction with psychometric tests allow for the investigation of the stages/levels of the decline as well as of the compensatory reserves in the visual-perceptive and cognitive-mental brain functions in order to determine the functional sensitivity thresholds of different brain functions to hypoxia. They open up the way for the functional characterization, the diagnosis and monitoring of brain insults or other acute and chronic pathological brain conditions. / Η παρούσα ερευνητική εργασία στοχεύει στην ανάπτυξη εργαλείων Βιοϊατρικής Νευρομηχανικής (Βιοτεχνολογίες) για την σε βάθος λειτουργική μελέτη, ταχεία διάγνωση, συνεχή παρακολούθηση και έγκαιρη αντιμετώπιση οξέων και χρόνιων εγκεφαλικών διαταραχών, ατόμων που υπόκεινται σε ή πάσχουν από οιαδήποτε μορφή συστηματικής υποξαιμίας ή πιο εντοπισμένης εγκεφαλικής υποξίας, καθώς και για την λειτουργική αξιολόγηση και το συνεχή έλεγχο της προσαρμοστικότητας κατά την εξάσκηση των «υψιβατών», και γενικότερα ατόμων που ασκούν δραστηριότητες και λειτουργούν μέσα σε περιβάλλοντα αυξημένων οπτικο-γνωστικο-κινητικών απαιτήσεων (ένα είδος «στρες τεστ» για τον εγκέφαλο). Για το σκοπό αυτό υποβάλλουμε ολόκληρο το οπτικογνωστικό σύστημα, από το στοιχειώδες αισθητηριακό έως το πιο πολύπλοκο νοητικό επίπεδο, σε μια πειραματική δοκιμασία κατηγορικής διάκρισης σύνθετων οπτικογνωστικών ερεθισμάτων, μετά από υπερταχεία οπτική διέγερση που οδηγεί στην έκλυση κινητικής απάντησης κατά την κατηγοριοποίηση στόχων (εικόνες «ζώων» εκλύουν παραγωγικές αποκρίσεις) και στην καταστολή της κατά την κατηγοριοποίηση μη-στόχων (εικόνες «μη-ζώων» εκλύουν ανασταλτικές αποκρίσεις). Οι ταλαντωτικές ηλεκτροφυσιολογικές αποκρίσεις που συγχρόνως καταγράφονται στις ινιακές-κροταφικές-βρεγματικές περιοχές του εγκεφάλου αναλύονται στο πεδίο του χρόνου (<20Hz) και στο συζευγμένο χρονοφασματικό πεδίο ευρυζωνικά (1-60Hz) με το συνεχή μετασχηματισμό του κυματίου που βελτιστοποιεί την πολυφασματική ανάλυση της υψίσυχνης (≥20Hz) γ-ταλαντωτικής δραστηριότητας. Αυτή η δοκιμασία οπτικογνωστικής κατηγοριοποίησης λαμβάνει χώρα τόσο σε νορμοξαιμικές όσο και υποξαιμικές συνθήκες (ελεγχόμενη μείωση στον κορεσμό του αίματος σε οξυγόνο από ≥97% γύρω στο 80% για 15 λεπτά κάτω από συνθήκες υποβαρικής υποξίας μέσα σε υποβαρικό θάλαμο), προκειμένου να ελέγξουμε ηλεκτροφυσιολογικούς δείκτες που μπορούν να ανιχνεύσουν και να συλλάβουν με τον πιο ευαίσθητο και δυναμικό τρόπο ακόμη και τόσο βραχύβιες και σχετικά ήπιες μεταβολές της εγκεφαλικής λειτουργίας. Η στατιστική παραμετρική ανάλυση των χρονοφασματικών χαρτών και η γενικευμένη, στατιστικά πιο ασφαλής, μέθοδος ανάλυσης των διακυμάνσεων ανέδειξαν ως πλέον ευαίσθητους και αξιόπιστους τους ακόλουθους δείκτες: τις κύριες αιχμές των προκλητών δυναμικών, τον παράγοντα φασικής συνάφειας των ταλαντώσεων μεταξύ των μοναδιαίων καταγραφών και την εκλυόμενη ενέργεια των προκλητών/φασικά-κλειδωμένων και επαγόμενων/ολικών ταλαντώσεων. Οι ηλεκτροφυσιολογικοί αυτοί δείκτες σε συνδυασμό με ψυχομετρικές δοκιμασίες επιτρέπουν τη διερεύνηση των σταδίων/επιπέδων κάμψης καθώς και των αποθεμάτων αντιρρόπησης των οπτικο-αντιληπτικών και γνωστικών-νοητικών λειτουργιών του εγκεφάλου για τον καθορισμό των λειτουργικών ουδών ευαισθησίας διάφορων εγκεφαλικών λειτουργιών στην υποξία. Ανοίγουν μάλιστα το δρόμο. για το λειτουργικό χαρακτηρισμό, τη διάγνωση και την παρακολούθηση εγκεφαλικών προσβολών ή άλλων οξέων και χρόνιων παθολογικών καταστάσεων του εγκεφάλου.