Andere Themen suchen…

Andere Themen suchen…

Beschleunigen Sie Ihre analytischen EEG-Zeitpläne mit schnell aufzubauenden, hochdichten kabellosen Arrays, die für den flexiblen Feldeinsatz optimiert sind.

Beschleunigen Sie Ihre analytischen EEG-Zeitpläne mit schnell aufzubauenden, hochdichten kabellosen Arrays, die für den flexiblen Feldeinsatz optimiert sind.

Die empirische Modenzerlegung (Empirical Mode Decomposition, EMD) wurde als Reaktion auf ein Missverhältnis zwischen den EEG-Verarbeitungswerkzeugen und den zu interpretierenden Daten entwickelt. Anstatt ein Signal in einen festen Satz von im Voraus ausgewählten Sinuswellen oder Wavelets zu zwingen, lässt die EMD die Daten ihre eigenen Bausteine definieren. Diese Bausteine werden als intrinsische Modenfunktionen (Intrinsic Mode Functions, IMFs) bezeichnet, und der Prozess, der sie erzeugt, erfordert keine Annahme darüber, ob das Signal stationär oder linear ist.

Dies macht die EMD konzeptionell attraktiv für das EEG, bei dem transiente, unregelmäßige Ereignisse oft genau die Merkmale sind, die ein Forscher erkennen möchte. Diese Attraktivität hat zu einer Reihe von angewandten Forschungsarbeiten in den Bereichen Anfallserkennung, Emotionsklassifikation, Artefaktentfernung und Gehirn-Computer-Schnittstellen geführt.

Beschleunigen Sie Ihre analytischen EEG-Zeitpläne mit schnell aufzubauenden, hochdichten kabellosen Arrays, die für den flexiblen Feldeinsatz optimiert sind.

Beschleunigen Sie Ihre analytischen EEG-Zeitpläne mit schnell aufzubauenden, hochdichten kabellosen Arrays, die für den flexiblen Feldeinsatz optimiert sind.

Warum EEG-Signale klassischen Signalverarbeitungswerkzeugen widerstehen

Das Kernproblem bei der Anwendung klassischer Methoden der Signalverarbeitung auf EEG ist struktureller Natur. Eine Fourier-basierte Analyse zerlegt ein Signal in eine feste Bibliothek von Sinus- und Kosinuswellen – ein Ansatz, der sauber funktioniert, wenn der zugrunde liegende Prozess sein statistisches Verhalten im Laufe der Zeit nicht ändert.

Das EEG erfüllt diese Bedingung nicht. Die Gehirnaktivität einer Person im ruhigen Ruhezustand sieht statistisch anders aus als die Aktivität derselben Person während eines Ausbruchs fokussierter Aufmerksamkeit oder, im klinischen Umfeld, während eines Krampfanfalls. Die Forschung führt diese nicht-stationäre Eigenschaft als primäre Rechtfertigung dafür an, sich der EMD anstelle von frequenzbasierten Werkzeugen zuzuwenden, die für stabile Signale entwickelt wurden.

Die praktische Konsequenz ist, dass eine Methode, die Stationarität voraussetzt, transiente neuronale Ereignisse verschmieren oder falsch darstellen kann, indem sie eine kurze Spitze so behandelt, als wäre sie gleichmäßig über ein längeres Fenster verteilt. Für Anwendungen wie die Epilepsiediagnostik, bei denen das Timing und die Form einer Spitze diagnostisches Gewicht haben, ist diese Art der Verschmierung eine echte Einschränkung.

Der Lösungsvorschlag von EMD besteht darin, die Annahme einer festen Basis vollständig zu vermeiden und die Zerlegung sich an die Form anpassen zu lassen, die das Signal tatsächlich annimmt.

Was die empirische Modenzerlegung (EMD) tut und wie sie funktioniert

Die EMD arbeitet mit einem wiederholten Prozess, der als Siebung (Sifting) bezeichnet wird. Der Algorithmus identifiziert die lokalen Maxima und Minima eines Signals, verbindet sie zu einer oberen und unteren Hüllkurve und subtrahiert den Durchschnitt dieser Hüllkurven vom Originalsignal.

Dieser Schritt wird so lange wiederholt, bis die resultierende Komponente zwei Bedingungen erfüllt:

  1. Ihre obere und untere Hüllkurve sind symmetrisch um Null

  2. Ihr Mittelwert bleibt über das gesamte Segment nahe bei Null.

Sobald eine Komponente diese Bedingungen erfüllt, wird sie als intrinsische Modenfunktion (IMF) bezeichnet und vom Originalsignal subtrahiert, damit der Prozess die nächste extrahieren kann.

Das Ergebnis ist ein Stapel von IMFs, geordnet von den schnellsten bis zu den langsamsten Schwingungen, gefolgt von einem endgültigen verbleibenden Trend, der den langsamen Drift erfasst. Da jede IMF direkt aus der lokalen Form der Daten extrahiert und nicht mit einer vorgegebenen Wellenform abgeglichen wird, kann die EMD im Prinzip unregelmäßige und sich verändernde Dynamiken darstellen, ohne dass der Analytiker im Voraus eine Basis auswählen muss.

Methode

Beste Bedingung

EMD

Starke Kontamination

ICA oder CCA

Leichtes Rauschen, einzelne Quelle

Wavelet-Transformation

Leichtes Rauschen, zwei Quellen

Wie EMD transiente Ereignisse wie Epilepsie-Spikes erfasst

Transiente Ereignisse im EEG, epileptiforme Spikes oder kurze Ausbrüche im Zusammenhang mit einer kognitiven Aufgabe, sind von herausragender klinischer und neurowissenschaftlicher Bedeutung, liegen jedoch oft unter der Hintergrundaktivität des Gehirns oder störendem Rauschen vergraben. Dies ist genau das Szenario, in dem die lokale, adaptive Zerlegung von EMD einen Vorteil gegenüber Methoden mit fester Basis bieten soll.

Eine Studie von Safieddine et al. testete dies direkt, indem sie epileptische Spike-ähnliche Signale simulierte, sie mit realen Muskelartefakten in verschiedenen Intensitäten kontaminierte und anschließend vier Korrekturmethoden verglich:

  1. EMD

  2. Unabhängige Komponentenanalyse (ICA)

  3. Kanonische Korrelationsanalyse (CCA)

  4. Wavelet-Transformation

Bei starker Kontamination übertraf EMD alle drei Alternativen. Dieser Vorteil war jedoch nicht absolut. Wenn das Signal weniger verrauscht war und die Spikes von einer einzigen kortikalen Quelle stammten, lieferten CCA und ICA bessere Ergebnisse. Wenn die Spikes aus zwei verschiedenen Quellen stammten, funktionierte bei starkem Rauschen entweder EMD oder ICA am besten, während die Wavelet-Transformation bei geringerem Rauschen besser abschnitt.

Mit anderen Worten: Das beste Werkzeug zur Artefaktbereinigung hing sowohl vom Rauschpegel als auch von der zugrunde liegenden Quellenkonfiguration ab, und der Vorteil der EMD war speziell im stark kontaminierten Fall am ausgeprägtesten.

Unterdessen gingen Forschungsstudien von Riaz et al. und Shafiul Alam et al. von der Artefaktkorrektur zur Klassifikation über und nutzten IMFs als Merkmalquelle zur Erkennung von Anfällen und Epilepsie.

Das Team von Riaz extrahierte aus jeder IMF zeitliche Momente, d. h. statistische Maße wie Varianz, Schiefe und Kurtosis, zusammen mit spektralen Merkmalen wie dem spektralen Schwerpunkt und der spektralen Schiefe. Diese Merkmale wurden in eine Support Vector Machine eingespeist, einen Standard-Klassifikator für überwachtes Lernen. Die resultierende Klassifikation von Epilepsiepatienten und Anfallssegmenten lieferte Ergebnisse, die von den Autoren als gut beschrieben wurden und sich im Vergleich zu anderen etablierten Methoden der Merkmalsextraktion behaupten konnten.

Das Team von Alam ging noch einen Schritt weiter und nutzte statistische Momente höherer Ordnung, die im EMD-Bereich berechnet wurden, und speiste sie in ein künstliches neuronales Netz ein. Bei mehreren Klassifikationsaufgaben, einschließlich der Unterscheidung von Anfallsaktivität und Nicht-Anfallsaktivität, meldete die Methode eine Genauigkeit, Sensitivität und Spezifität von bis zu 100 Prozent und lief dabei schneller als vergleichbare Zeit-Frequenz-Analyseverfahren.

Zusammengenommen deuten diese Studien darauf hin, dass die EMD helfen kann, transiente epileptische Aktivität zu isolieren, aber ihre Leistung hängt empfindlich vom Rauschpegel, der Geometrie der Quelle und dem getesteten Datensatz ab, anstatt einen festen, garantierten Vorteil zu bieten.

Häufige Herausforderungen bei der empirischen Modenzerlegung

Die EMD ist flexibel, aber diese Flexibilität erfordert methodische Entscheidungen, die wissenschaftliche Schlussfolgerungen beeinflussen können. Kleine Störungen in einem Signal können die Extrema verändern und somit die Hüllkurven und die extrahierten Moden beeinflussen. Rauschen, Filterung, kurze Aufzeichnungen und scharfe Grenzen können die Zerlegung instabiler machen.

Gängige Überprüfungen helfen, diese Schwierigkeiten offenzulegen statt sie zu verbergen. Dazu gehören:

  • Messung des Rekonstruktionsfehlers nach Summierung der IMFs und des Residuums.

  • Prüfung der Empfindlichkeit gegenüber der Behandlung von Endpunkten und Siebungsschwellenwerten.

  • Vergleich von Zerlegungen über wiederholte Segmente oder Rauschensembles hinweg.

  • Prüfung, ob eine entfernte Komponente plausible neuronale Aktivität enthält.

Diese Prüfungen beseitigen Unklarheiten nicht, aber sie machen sie sichtbar. Das Mischen von Moden (Mode Mixing) kann eine Ensemble-Variante erfordern, während ein übermäßiger hochfrequenter Anteil eher auf eine Kontamination als auf einen nützlichen neuronalen Modus hindeuten kann. In jedem Fall sollte die Analyse die Entscheidungen dokumentieren, die zu der endgültigen Darstellung geführt haben.

Eine weitere Herausforderung ist die Kommunikation. Ein IMF-Index ist ein algorithmisches Label, keine Diagnose oder ein universell definiertes EEG-Band. Klare Abbildungen, Parameteraufzeichnungen und eine vorsichtige Sprache sind erforderlich, wenn Ergebnisse aus der explorativen Signalverarbeitung in die biomedizinische Interpretation übergehen.

Mode Mixing: Die Hauptschwachstelle der grundlegenden EMD

Die deutlichste strukturelle Schwachstelle der grundlegenden EMD ist ein Phänomen namens Mode Mixing (Modenmischung). Es äußert sich auf zwei verwandte Arten: Entweder werden Schwingungen, die auf sehr unterschiedlichen Zeitskalen ablaufen, in derselben IMF eingeschlossen, oder eine einzelne, kohärente rhythmische Komponente wird fragmentiert und über mehrere IMFs verstreut, anstatt in einer einzigen enthalten zu bleiben. Beide Varianten des Problems erschweren das Vertrauen darauf, dass eine IMF einem echten, einzigartigen physiologischen Prozess entspricht.

Das praktische Risiko ist interpretatorischer Natur. Wenn ein Forscher davon ausgeht, dass eine bestimmte IMF beispielsweise die Aktivität in einem bestimmten Frequenzband darstellt, diese IMF jedoch aufgrund von Mode Mixing tatsächlich eine Mischung aus Muskelartefakten und verbleibender Spike-Aktivität enthält, wird jede aus dieser IMF gezogene Schlussfolgerung unzuverlässig.

Complete Ensemble Empirical Mode Decomposition

Die Complete Ensemble Empirical Mode Decomposition, allgemein als CEEMD bezeichnet, ist eine rauschunterstützte Erweiterung der EMD. Sie fügt dem Signal sorgfältig kontrollierte Rauschrealisierungen hinzu, zerlegt das resultierende Ensemble und mittelt die entsprechenden Moden. Ziel ist es, die Zerlegung stabiler zu machen, wenn nahe beieinander liegende Schwingungsskalen schwer zu trennen sind.

Wie CEEMD die Standard-EMD verbessert

Die Standard-EMD kann unter Mode Mixing leiden, bei dem Schwingungen mit unterschiedlichen Skalen in einer IMF auftreten oder ähnliche Skalen über mehrere IMFs verteilt sind. CEEMD geht diese Instabilität statistisch an, anstatt das grundlegende Konzept der Hüllkurvensiebung zu ändern. Wiederholte Zerlegungen und Mittelwertbildungen können den Einfluss einer einzelnen Rauschrealisierung verringern, aber das Ergebnis hängt immer noch vom Ensemble-Design, der Rauschamplitude, der Anzahl der Versuche und Implementierungsdetails ab.

Die Methode ist nicht automatisch für jeden Datensatz überlegen. Sie erhöht den Rechenaufwand und kann Restrauschen beibehalten oder eine Randempfindlichkeit einführen, wenn ihre Parameter schlecht gewählt sind. Für das EEG sollten Vergleiche daher die Rekonstruktionsqualität, die Stabilität über die Versuche hinweg, die Erhaltung transienter Ereignisse und die Frage untersuchen, ob die resultierenden Moden für die jeweilige Forschungsfrage interpretierbar bleiben.

Welche EEG-Anwendungen EMD bisher getestet haben

Die EMD wurde üblicherweise in drei großen Kategorien angewendet.

Merkmalsextraktion für die Klassifikation stellt den größten Anwendungsfall dar. Beispielsweise nutzte eine von Ning Zhuang geleitete Studie EMD zur Emotionserkennung, indem sie mehrdimensionale Merkmale aus IMFs extrahierte – darunter die erste Differenz der Zeitreihe, die erste Phasendifferenz und die normierte Energie – und diese dann an der öffentlich zugänglichen emotionalen DEAP-Datenbank testete. Die Studie ergab, dass die hochfrequente IMF1-Komponente einen signifikanten Einfluss auf die Unterscheidung emotionaler Zustände hatte, und der Gesamtansatz übertraf Vergleichsmethoden wie die fraktale Dimension, die Stichprobenentropie und die diskrete Wavelet-Transformation.

Artefaktbereinigung ist in Studien wie der von Safieddine et al. vertreten, in der die deterministische, hüllkurvenbasierte Zerlegung der EMD Muskelartefakte aus simuliertem epileptischem EEG effektiver entfernte als ICA, CCA oder Wavelet-basierte Methoden, insbesondere unter Bedingungen starker Kontamination.

Schließlich sind Mehrkanalanalysen und Gehirn-Computer-Schnittstellen in einer im Journal of Circuits, Systems and Computers veröffentlichten Arbeit zu finden. Hierbei wurde die EMD separat auf jeden Kanal einer mehrkanaligen EEG-Aufzeichnung angewendet, und die resultierenden IMFs wurden nach spektraler Ähnlichkeit gruppiert, um zu identifizieren, welche Komponenten mit einem externen Stimulus korrelierten.

Die Rekonstruktion des Signals nur aus diesen stimulusbezogenen Komponenten ermöglichte es den Forschern zu erkennen, wann eine Testperson aufmerksam auf einen bestimmten Stimulus reagierte – eine Fähigkeit, die direkt für BCI-Systeme relevant ist, die darauf angewiesen sind, beabsichtigte neuronale Signale vom Hintergrundrauschen zu unterscheiden.

Tutorial zur empirischen Modenzerlegung für Anfänger

Ein Anfänger kann sich die EMD wie das wiederholte Abschälen der schnellsten sichtbaren Schwingung von einer komplexen Wellenform vorstellen. Der Algorithmus beginnt nicht mit einer Liste von Frequenzbändern. Stattdessen beobachtet er, wo das Signal steigt und fällt, schätzt eine lokale Mittellinie und entfernt diese Mittellinie von der Kandidatenschwingung.

Eine sinnvolle Lernübung beginnt mit einem kurzen, sauberen synthetischen Signal, das zwei Schwingungen und einen langsamen Trend enthält. Das erwartete Ergebnis ist keine perfekte Trennung in jedem Fall, sondern ein Satz von Komponenten, deren Summe das Eingangssignal annähernd reproduziert. Das Auftragen des Originalsignals, jeder IMF und des Residuums auf einer gemeinsamen Zeitachse ist oft aufschlussreicher als die reine Betrachtung numerischer Werte.

Bei einer EEG-Aufzeichnung umfasst der Arbeitsablauf auch Entscheidungen, die vor der Zerlegung getroffen werden müssen:

  • Auswahl von Kanälen und Segmenten

  • Überprüfung von Abtastinformationen

  • Markierung fehlerhafter Intervalle

  • Dokumentation der Vorverarbeitung

Ein allgemeiner Leitfaden zur Frequenzzerlegung kann helfen, die EMD zwischen Fourier-, Wavelet- und verwandten Zeit-Frequenz-Methoden einzuordnen. Diese Entscheidungen machen aus der EMD keinen klinischen Test; sie definieren die Bedingungen, unter denen ihre Ergebnisse bewertet werden können.

Empirische Modenzerlegung in Python: Ein Programmier-Leitfaden

Python-Implementierungen akzeptieren in der Regel ein eindimensionales Array und geben eine Sammlung von IMFs zurück, manchmal mit einem Residuum oder durch eine separate Restwertoperation. Ein grundlegender Arbeitsablauf lädt ein Signal, bestätigt dessen Abtastreihenfolge und den Status fehlender Werte, führt die Zerlegung aus und stellt die Komponenten grafisch dar. Eine reproduzierbare Analyse zeichnet auch die Bibliotheksversion und die Parameter auf, die das Siebungsverhalten und das Ensemble-Verhalten steuern.

Beliebte Python-Bibliotheken für EMD

Mehrere wissenschaftliche Python-Pakete bieten EMD oder verwandte Algorithmen an, aber ihre APIs und Standardwerte unterscheiden sich. Einige bieten Standard-EMD, Ensemble-Varianten, Abbruchschwellenwerte, Spline-Optionen und Visualisierungshilfen; andere erfordern mehr manuelle Einrichtung. Die zentrale Programmieraufgabe besteht nicht nur im Aufrufen einer Funktion, sondern im Überprüfen, was die zurückgegebenen Komponenten bedeuten und ob ihre Summe das Eingangssignal rekonstruiert.

Ein minimales Pseudocode-Muster sieht wie folgt aus:

signal = load_signal()
check_sampling_and_missing_values(signal)
imfs, residue = decompose(signal)
reconstruction = imfs.sum(axis=0) + residue
plot_components(signal, imfs, residue)
signal = load_signal()
check_sampling_and_missing_values(signal)
imfs, residue = decompose(signal)
reconstruction = imfs.sum(axis=0) + residue
plot_components(signal, imfs, residue)
signal = load_signal()
check_sampling_and_missing_values(signal)
imfs, residue = decompose(signal)
reconstruction = imfs.sum(axis=0) + residue
plot_components(signal, imfs, residue)
signal = load_signal()
check_sampling_and_missing_values(signal)
imfs, residue = decompose(signal)
reconstruction = imfs.sum(axis=0) + residue
plot_components(signal, imfs, residue)
signal = load_signal()
check_sampling_and_missing_values(signal)
imfs, residue = decompose(signal)
reconstruction = imfs.sum(axis=0) + residue
plot_components(signal, imfs, residue)

Der Code ist absichtlich schematisch gehalten, da sich Funktionsnamen und Rückgabewerte je nach Paket unterscheiden. Die Validierung sollte einen Rekonstruktionsfehler, Darstellungen von Extrema und Grenzen sowie Empfindlichkeitsprüfungen unter angemessenen Parameteränderungen umfassen.

Für das EEG spielen auch die Datendimensionen eine Rolle: Eine mehrkanalige Aufzeichnung erfordert möglicherweise einen kanalweisen Arbeitsablauf oder eine speziell multivariate Methode, anstatt alle Kanäle stillschweigend in ein einziges Array flachzudrücken.

Anwendung der empirischen Modenzerlegung auf EEG-Signale

Das EEG ist eine natürliche Umgebung für die empirische Modenzerlegung, da Spannungsänderungen überlagernde Rhythmen, transiente Ereignisse und nicht-stationäre Aktivitäten enthalten können. Die EMD kann auf einzelne Kanäle oder ausgewählte Segmente angewendet werden, um zu untersuchen, wie sich lokale Schwingungen entwickeln. Die relevante wissenschaftliche Frage kann die zeitliche Struktur, die Artefakttrennung, ereigniskorrelierte Änderungen oder Merkmale für ein nachfolgendes statistisches Modell betreffen.

Vor der Zerlegung erfordern die Signalqualität und der Aufzeichnungskontext Aufmerksamkeit. Die Platzierung der Elektroden, das Referenzschema, die Abtastrate, die Filterung, fehlende Abtastwerte und markierte Artefakte können alle die lokalen Extrema beeinflussen.

Zudem sollten den resultierenden IMFs physiologische Bezeichnungen nicht allein aufgrund ihrer Reihenfolge zugewiesen werden. Schnelle Komponenten können Muskel- oder Instrumentenrauschen enthalten, und langsamere Komponenten können sowohl Drift als auch sinnvolle Aktivität beinhalten.

Die Interpretation ist fundierter, wenn sie rekonstruierte Signale vergleicht, ein klares Protokoll der ausgeschlossenen Komponenten führt und testet, ob die Ergebnisse über Kanäle, Teilnehmer und Analyseentscheidungen hinweg stabil bleiben.

Empirische Modenzerlegung vs. Wavelet-Transformation

Sowohl die EMD als auch die Wavelet-Transformation unterstützen die Zeit-Frequenz-Analyse, organisieren die Signalstruktur jedoch unterschiedlich. Die EMD leitet Moden aus lokalen Extrema der beobachteten Daten ab, wohingegen die Wavelet-Analyse das Signal mit einer ausgewählten Familie von skalierten und verschobenen Basisfunktionen vergleicht. Die eine Methode passt sich den Daten an, die andere stellt ein vordefiniertes mehrskaliges Koordinatensystem bereit.

Der Unterschied wird deutlicher, wenn man die analytischen Prioritäten vergleicht:

Merkmal

Empirische Modenzerlegung

Wavelet-Transformation

Darstellung

Aus den Daten abgeleitete intrinsische Modenfunktionen

Koeffizienten auf ausgewählten Wavelet-Skalen

Basis

Adaptiv und signalabhängig

Vorgegebene Wavelet-Familie

Zeitliche Lokalisierung

Lokal, mit datenabhängigen Moden

Skalenabhängig und durch das Wavelet gesteuert

Hauptprobleme

Mode Mixing, Endpunkteffekte, Abbruchregeln

Basiswahl, Auflösungskompromisse, Randeffekte

Keine der beiden Methoden ist universell überlegen. Die EMD kann eine intuitive Zerlegung für sich verändernde Schwingungen bieten, während Wavelets einen konsistenten Rahmen für den Vergleich von Signalen liefern können, wenn die Basis und die Skalen wissenschaftlich angemessen sind.

Die Wahl sollte sich nach der Fragestellung, den erwarteten Transienten, dem Bedarf an Reproduzierbarkeit und der Validierungsstrategie richten und nicht nach der pauschalen Behauptung, eine Methode sei fortschrittlicher.

Was die Analyse von Gehirnsignalen durch datengestützte Zerlegung gewinnt

Die Gehirnaktivität verweilt selten lange genug, als dass klassische Frequenzwerkzeuge ihre aufschlussreichsten Momente erfassen könnten. Aus diesem Grund hat eine Methode, die es den Signalen erlaubt, ihre eigenen Bausteine zu definieren, große Aufmerksamkeit erregt.

Die empirische Modenzerlegung umgeht die Notwendigkeit fester mathematischer Vorlagen und liefert Ergebnisse, die sich bei der Isolierung von Epilepsie-Spikes unter starkem Rauschen und der Extraktion aussagekräftiger Merkmale für die Emotionserkennung als nützlich erwiesen haben. Dieser Übergang zu einer adaptiv gestalteten Analyse eröffnet einen Weg, transiente neuronale Ereignisse zu erfassen, die andernfalls verschmiert oder übersehen würden. Das praktische Ergebnis, wie es sich in den Studien zeigt, ist ein Werkzeug, das die Klassifikation und Artefaktbereinigung unterstützen kann, wenn die Bedingungen zu seinen Gunsten zusammenspielen.

Referenzen

  1. Safieddine, D., Kachenoura, A., Albera, L., Birot, G., Karfoul, A., Pasnicu, A., ... & Merlet, I. (2012). Removal of muscle artifact from EEG data: comparison between stochastic (ICA and CCA) and deterministic (EMD and wavelet-based) approaches. EURASIP Journal on Advances in Signal Processing, 2012(1), 127. https://doi.org/10.1186/1687-6180-2012-127

  2. Riaz, F., Hassan, A., Rehman, S., Niazi, I. K., & Dremstrup, K. (2015). EMD-based temporal and spectral features for the classification of EEG signals using supervised learning. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 24(1), 28-35. https://doi.org/10.1109/TNSRE.2015.2441835

  3. Alam, S. S., & Bhuiyan, M. I. H. (2013). Detection of seizure and epilepsy using higher order statistics in the EMD domain. IEEE journal of biomedical and health informatics, 17(2), 312-318. https://doi.org/10.1109/JBHI.2012.2237409

  4. Zhuang, N., Zeng, Y., Tong, L., Zhang, C., Zhang, H., & Yan, B. (2017). Emotion recognition from EEG signals using multidimensional information in EMD domain. BioMed research international, 2017(1), 8317357. https://doi.org/10.1155/2017/8317357

  5. Rutkowski, T. M., Mandic, D. P., Cichocki, A., & Przybyszewski, A. W. (2010). EMD approach to multichannel EEG data—the amplitude and phase components clustering analysis. Journal of Circuits, Systems, and Computers, 19(01), 215-229. https://doi.org/10.1142/S0218126610006037

Häufig gestellte Fragen (FAQ)

Was ist die empirische Modenzerlegung (EMD) und warum ist sie für das EEG nützlich?

Die EMD ist eine adaptive Methode der Signalverarbeitung, die eine Wellenform in oszillatorische Komponenten, sogenannte intrinsische Modenfunktionen, zerlegt, ohne anzunehmen, dass das Signal stationär ist. Sie ist für das EEG nützlich, da sich die Gehirnaktivität kontinuierlich verändert und transiente Ereignisse enthält, deren genaue Darstellung klassischen Werkzeugen mit fester Basis schwerfällt.

Wie funktioniert der Siebprozess (Sifting) bei der EMD?

Der Siebprozess identifiziert wiederholt die lokalen Maxima und Minima eines Signals, konstruiert eine obere und untere Hüllkurve und subtrahiert deren Durchschnitt vom Signal. Dieser Schritt wird so lange wiederholt, bis die verbleibende Komponente symmetrische Hüllkurven und einen Mittelwert nahe Null aufweist. An diesem Punkt wird sie zu einer intrinsischen Modenfunktion.

Was sind intrinsische Modenfunktionen (IMFs)?

IMFs sind die datengestützten Bausteine, die durch die EMD erzeugt werden, geordnet von den schnellsten bis zu den langsamsten Schwingungen. Da sie direkt aus der lokalen Form des Signals extrahiert und nicht an eine vorgegebene Wellenform angepasst werden, können sie unregelmäßige und sich verändernde neuronale Dynamiken erfassen.

Warum ist Mode Mixing ein Problem bei der EMD?

Mode Mixing tritt auf, wenn Schwingungen unterschiedlicher Zeitskalen in derselben IMF eingeschlossen werden oder wenn ein einzelner Rhythmus über mehrere IMFs aufgeteilt wird. Diese Unklarheit macht es schwer darauf zu vertrauen, dass eine IMF einen einzelnen physiologischen Prozess darstellt, wodurch Rauschen oder Artefakte mit echter Gehirnaktivität vermischt werden können.

Wie geht die EMD mit transienten Ereignissen wie Epilepsie-Spikes um?

Die EMD isoliert kurze Spikes, indem sie diese in hochfrequente IMFs extrahiert, wodurch die Hintergrundaktivität abgetrennt wird. Bei starker Muskelkontamination übertraf die EMD andere Methoden zur Artefaktbereinigung, wenngleich ihr Vorteil bei saubereren Signalen oder anderen Quellenkonfigurationen abnahm.

Kann die EMD traditionelle Methoden wie die Fourier-Transformation bei allen EEG-Analysen ersetzen?

Die EMD ist kein universeller Ersatz; sie ist am wertvollsten bei nicht-stationären und transienten Signalen, bei denen Fourier-Methoden mit fester Basis an ihre Grenzen stoßen. Ihre Zuverlässigkeit variiert jedoch mit dem Rauschpegel und ist noch nicht für alle EEG-Aufzeichnungsbedingungen nachgewiesen.

Erfordert die EMD die Annahme, dass das Signal stationär ist?

Nein, die EMD setzt keine Stationarität voraus, was einen entscheidenden Vorteil gegenüber der klassischen Fourier-Analyse für das EEG darstellt. Die Zerlegung passt sich direkt dem sich ändernden statistischen Charakter der Daten an, sodass sie sich verschiebende Gehirndynamiken und kurze Ereignisse ohne Vorbehandlung nachverfolgen kann.

Kann die EMD mehrkanaliges EEG direkt analysieren?

Die grundlegende EMD wird üblicherweise auf eindimensionale Signale angewendet. Eine Mehrkanalanalyse erfordert eine bewusste kanalweise oder multivariate Strategie, da das einfache Kombinieren von Kanälen räumliche und zeitliche Bedeutungen verschleiern kann.

Beschleunigen Sie Ihre analytischen EEG-Zeitpläne mit schnell aufzubauenden, hochdichten kabellosen Arrays, die für den flexiblen Feldeinsatz optimiert sind.

Beschleunigen Sie Ihre analytischen EEG-Zeitpläne mit schnell aufzubauenden, hochdichten kabellosen Arrays, die für den flexiblen Feldeinsatz optimiert sind.

Emotiv ist ein führender Anbieter von Neurotechnologie, der die neurowissenschaftliche Forschung mit zugänglichen EEG- und Gehirndaten-Tools vorantreibt.

Christian Burgos

Das Neueste von uns

Leistungsspektraldichte im EEG

Die spektrale Leistungsdichte (Power Spectral Density, kurz PSD) ist das Werkzeug, das EEG-Signale entmischt und Ihnen zeigt, wie viel Energie jede dieser Geschwindigkeiten oder Frequenzen zur Gesamtaufzeichnung beiträgt. Sobald Sie ein PSD-Diagramm lesen können, lesen Sie eine Art Rhythmus-Partitur für das Gehirn – ein Diagramm, das zeigt, welche Tempi dominieren und welche im Hintergrund verblassen.

Artikel lesen

Die diskrete Kosinustransformation

Ein Elektroenzephalogramm (EEG) erzeugt über lange Zeiträume hinweg riesige Mengen an kontinuierlichen Daten über Dutzende von Kanälen. Dieses Volumen belastet den begrenzten Speicher tragbarer Headsets, schränkt Telemedizin-Netzwerke ein und verlangsamt Echtzeit-Gehirn-Computer-Schnittstellen (BCIs). Folglich müssen EEG-Rohdaten reduziert werden, um effizient verarbeitet werden zu können.

Die diskrete Kosinustransformation (DCT), die als mathematische Grundlage für die JPEG-Komprimierung dient, löst diese Herausforderung. Ebenso wie sie Bilder komprimiert und gleichzeitig die Erkennbarkeit beibehält, reduziert die DCT die EEG-Signalgröße, während sie dessen Gesamtform bewahrt. Das Verständnis ihrer Funktionsweise und Grenzen hilft bei der Entscheidung, wann die DCT geeignet ist oder wann alternative Transformationen vorzuziehen sind.

Artikel lesen

Ein Leitfaden zur Kurzzeit-Fourier-Transformation für EEG

Eine grundlegende Fourier-Transformation, das klassische mathematische Werkzeug, um ein Signal in seine einzelnen Frequenzkomponenten zu zerlegen, kann Ihnen sagen, welche Gehirnrhythmen irgendwo während einer gesamten Aufzeichnung vorhanden waren. Was sie Ihnen nicht sagen kann, ist, wann sie aufgetreten sind. Ein kurzer Ausbruch von Alpha-Aktivität, der in dem Moment auftritt, in dem jemand die Augen schließt, ist ein bedeutsames, zeitlich synchronisiertes Ereignis.

Gemittelt zu einer einzigen Frequenzzusammenfassung über eine zehnminütige Aufzeichnung wird dieser Ausbruch zu einer Statistik, die ununterscheidbar vom Hintergrundrauschen ist. Die zeitliche Rekonstruktion dieser Ereignisse ist der Ausgangspunkt für jede ernsthafte EEG-Forschung, und es ist genau das Problem, für dessen Lösung die Kurzzeit-Fourier-Transformation (STFT) entwickelt wurde.

Artikel lesen

Die Fourier-Transformation

Die Fourier-Transformation verändert die Art und Weise, wie bestehende Informationen dargestellt werden, indem sie ein Signal, das sich im Laufe der Zeit verändert, in eine Beschreibung der rhythmischen Komponenten umwandelt, die bereits darin verborgen sind. Diese Transformation als Linse und nicht als Gerät zu verstehen, ist der notwendige erste Schritt, bevor ein Gespräch über Gehirnrhythmen, Frequenzbänder oder spektrale Muster überhaupt Sinn ergeben kann.

Artikel lesen