Rok / Year: 2014
Svazek / Volume: 16
Číslo / Number: 2
EEG biofeedback - návrh a realizace v prostředí LabVIEW EEG Biofeedback – design and implementation in LabVIEW Jiří Dvořák, Vratislav Čmiel
[email protected],
[email protected] Fakulta elektrotechniky a komunikačních technologií VUT v Brně
Abstrakt: Cílem této práce byl návrh a realizace pracoviště EEG biofeedback pro výukové i cvičné účely, a to za pomoci dostupných univerzitních prostředků. Vzhledem k požadavku práce v reálném čase byl návrh proveden pomocí akviziční jednotky MP 35 Biopac, univerzální měřící karty NI USB6221 a prostředí LabVIEW.
Abstract: The aim of this work was to design and implement an EEG Biofeedback training program for educational and training purposes, using university resources available. To meet the requirement of real-time processing, the system was made by means of the MP 35 Biopac acquisition unit, NI USB-6221 universal measuring card and LabVIEW.
VOL.16, NO.2, APRIL 2014
EEG biofeedback - návrh a realizace v prostředí LabVIEW Jiří Dvořák, Vratislav Čmiel Fakulta elektrotechniky a komunikačních technologií VUT v Brně Email:
[email protected],
[email protected]
Abstrakt – Cílem této práce byl návrh a realizace pracoviště EEG biofeedback pro výukové i cvičné účely, a to za pomoci dostupných univerzitních prostředků. Vzhledem k požadavku práce v reálném čase byl návrh proveden pomocí akviziční jednotky MP 35 Biopac, univerzální měřící karty NI USB-6221 a prostředí LabVIEW.
Z profesionálních aplikací je pak nejvíce rozšířený Deymed Diagnostic BRAINFEEDBACK nebo Mind WorkStation. Obecný princip jejich funkce je vždy stejný a lze ho odvodit z obrázku 1, které je zároveň blokovým schématem námi navrženého pracoviště EBFB. To se skládá z následujících hardwarových komponent: 1.
1 Úvod EEG biofeedback (dále jen EBFB) je terapeutická metoda léčby lehkých mozkových dysfunkcí (LMD) [1], využívající kontinuální vyhodnocování elektrických potenciálů šedé kůry mozkové. Jedná se o vysoce specifickou metodu pro posílení žádoucí aktivace nervové soustavy, především pro trénink pozornosti a soustředění, sebeovládání a sebekázně, zklidnění impulzivity a hyperaktivity, a zlepšení výkonů intelektu. Na rozdíl od farmakoterapie je účinek EEG biofeedbacku trvalý, protože se jedná o proces učení – co se mozek při tréninku naučí, už se nelze "odnaučit" [1, 2]. Metoda upravuje chyby v autonomní regulaci fyziologických funkcí a má princip ve snímání signálu EEG z těla pacienta a jeho následného předkládání zpět pacientovi, a to ve vizuální nebo zvukové podobě. Pacientovi je tak okamžitě poukázáno na přirozenou či patologickou činnost mozku. Ten se pak svou vůlí snaží odchylku korigovat [8]. Jedná se tedy o samostatnou regulaci mozkové aktivity za přispění cílené inhibice specifických frekvencí EEG a zesilování jiných. Obecné se provádí dvě formy tréninku: podle stimulačního neboli výcvikového modelu (trénink pozornost) a podle antistresového modelu (trénink relaxace) [8, 11]. Metoda ke své realizaci potřebuje hardwarové a současně i softwarové technické prostředky, což z ní dělá typickou úlohu biomedicínského inženýrství. Zpětná vazba bývá většinou -převážně v případě dětských pacientů- vytvořena pomocí interaktivní počítačové hry. Dítě dostává okamžitou informaci o stavu své mozkové aktivity pomocí měnícího se parametru zvoleného herního motivu (výška letícího letadla, hladina vody v řece apod.). Forma a parametry tréninku jsou nastaveny terapeutem vždy tak, aby bylo dítě motivováno zvyšujícím se herním skóre ke korekci odchylky správným směrem [8]. Bližší popis tréninku a vysvětlení principu terapie uvádí kapitola 5.
2.
3.
Akviziční jednotka MP35, sdružená s programovým balíkem Biopac Student Lab PRO [4], slouží ke sběru analogových dat nejrůznějšího biologického charakteru. Jejím úkolem je zesílení (50 000x) a filtrace pomocí horní propusti o mezní frekvenci 0,5 Hz (zájmové pásmo spektra EEG signálu je v našem případě 0,5 Hz až 35 Hz) [4]. Jednotka utlumí artefakty spojené s pomalým kolísáním nulové izolinie napětí, způsobené změnou impedance vlivem dýchacího cyklu či měnícím se půlčlánkovým napětím elektrod. Signál snímáme pomocí nalepovacích stíněných elektrod typu EL254S nebo EL502. Zapojení elektrod určuje MP35 akviziční jednotka jako bipolární. Pár aktivních diferentních elektrod je umístěno v oblasti hlavy podle aktuálního tréninkového modelu (většinou v oblasti Cz a Pz centrální části temene podle systému rozmístění elektrod 10/20) a další referenční (zemnicí) bývá umístěna do oblasti nejmenší elektrické aktivity (ušní lalůček) [5]. Měří se rozdíl napětí mezi dvěma aktivními elektrodami vztaženo k elektrodě zemnící. Měřící karta NI USB-6221 je využita jako A/D převodník signálu jdoucího s analogového výstupního portu akviziční jednotky MP35. Karta komunikuje s osobním počítačem přes USB sběrnici, využívá k tomu balík ovladačů skupiny NI-DAQmx. Osobní počítač se spuštěnou aplikací EBFB a dvěma monitory. Aplikace EBFB, vytvořená v prostředí LabVIEW, se stará o konečnou filtraci signálu, jeho rozdělení do EEG pásem a výpočet výkonu v těchto pásmech. Na základě tohoto výkonu je měněn herní parametr na Monitoru 2, který zobrazuje herní motiv pro pacienta. Monitor 1 naopak obsahuje ovládací panel, pomocí kterého terapeut sleduje vývoj snímaného EEG a nastavuje parametry tréninku.
3 Zpracování signálu
2 Schéma pracoviště
Konečná filtrace signálu EEG je provedena softwarově až v samotné aplikaci. První čtveřice pásmových propustí tvoří banku číslicových IIR filtrů a vymezuje jednotlivá EEG pásma
Dnes je dostupná celá řada zařízení pro domácí trénink EBFB (NEUROSET BrainGym, Laxball, BioGraph Infiniti).
73
VOL.16, NO.2, APRIL 2014
NI USB-6221 měřicí karta (digitalizace)
MP35
MONITOR 2
akviziční jednotka (zesílení EEG, potlačení offsetu, galvanické oddělení)
aplikace HRY.VI (vizuální zpětná vazba)
USB 2.0
MONITOR 1 aplikace EEG BIOFEEDBACK.VI (analýza EEG, řízení hry)
PC vývojové prostředí LabView 7.1
Obrázek 1: Blokové schéma hardwarové části pracoviště EBFB v pořadí Delta 0,5 - 4Hz, Theta 4 – 8 Hz, Alfa 8 – 13 Hz, Beta 13 – 22 Hz) [9]. Tato klasifikace frekvenční škály slouží čistě pro informativní zobrazení časových průběhů jednotlivých pásem na monitoru terapeuta a neovlivňují herní motiv. Proto byly použity eliptické filtry 5. řádu typu IIR. IIR filtry potřebují oproti FIR filtrům mnohem menší řád přenosové funkce k dosažení totožné odezvy. Klesají tak nároky na paměť při výpočtu koeficientů a stavových proměnných a tím i zpoždění při zpracování vstupních vzorků. Poslední pátý FIR filtr typu dolní propust s mezní frekvencí 40 Hz filtruje vstupní signál pro následné zpracování, vyhodnocení a výpočet úspěšností terapeutických tréninků. Z tohoto důvodu si nepřejeme velké zpoždění ve filtrační větvi a současně požadujeme minimální zkreslení signálu a hlavně stabilitu filtru. Typ filtru s konečnou impulsní charakteristikou je tedy zvolen zejména pro jeho stabilitu (neosciluje), možnost lineární fázové charakteristiky a tím i konstantního skupinového zpoždění (nemění tvar signálu). Sice FIR více zpožďuje vstupní signál oproti druhému typu, ale jelikož zde není -jak bude popsáno později- prioritou strmost přechodu z propustného pásma do pásma zádržného, vystačíme si s délkou impulsní charakteristiky rovnou padesáti vzorkům. Zpoždění signálu je proto rovno 25-ti vzorkům (polovina délky impulsní charakteristiky filtru), což je při vzorkovací frekvenci 250 Hz jedna desetina vteřiny. V našem případě je nejvíce směrodatným parametrem (hned vedle zpoždění výstupního signálu) plochost amplitudové frekvenční charakteristiky. Projeví-li se vlivem zvlnění plata modulu frekvenční charakteristiky nějaké diference mezi útlumem jednotlivých harmonických komponent, může to vést až k chybnému vyhodnocení úspěšnosti tréninku. Z tohoto důvodu je třeba vhodným váhováním impulsní charakteristiky posledního FIR filtru toto zvlnění eliminovat. Volíme takové váhové okno, které ve výsledku způsobí nejmenší zvlnění. Nejvhodnějším adeptem je Kaiserovo okno, které vykazuje největší útlum vedlejších laloků frekvenční charakteristiky. Výsledná amplitudová frekvenční charakteristika tvořena dílčími filtry použitými v aplikaci je znázorněna na obrázku 2. V aplikaci EEG Biofeedback realizuje výpočet individuálních výkonových spekter blok subVI s názvem Auto Power Spectrum.VI [3], což je v podstatě algoritmus výpočtu jednostranného výkonového spektra Sff pomocí součinu komplexního Fwi a komplexně sdruženého spektra *Fwi (1). Symbol wi
zde udává index příslušné spektrální realizace a N zase počet vzorků v jednom úseku vstupního signálu. Výsledná hodnota složek výkonového spektra je udávána v jednotce Vrms2. Umocnění N2 je zde ve významu jednostranného výkonového spektra. kΩ, nT a mT vyjadřují vzorky signálu ve frekvenční, resp. časové diskrétní oblasti. 2
𝑤 𝑖 𝑆𝑓𝑓
=
𝐹𝑤 𝑖 𝐹𝑤 𝑖 𝑘Ω ∙ 𝐹𝑤∗𝑖 𝑘Ω 𝑘Ω = = = 𝑁2 𝑁2 𝑁−1 −𝑗𝑘 Ω𝑛𝑇 𝑁−1 ∙ 𝑚 =0 𝑓𝑤 𝑖 (𝑚)𝑒 𝑗𝑘 Ω𝑚𝑇 𝑛=0 𝑓𝑤 𝑖 (𝑛)𝑒 𝑁2
(1)
Každý úsek signálu váhujeme obdélníkovým oknem o délce N=250 vzorků1. Ačkoliv lze zvyšováním délky okna dosáhnout většího frekvenčního rozlišení (časové rozlišení se vlivem překrývání oken stává na délce okna do značné míry invariantní), zhoršil by se dynamický rozsah a také rychlost odezev jednotlivých složek spektra na náhlé změny amplitud v odpovídajících kmitočtových pásmech (zejména pak těch vyšších). S přihlédnutím k tomuto omezení byla optimální délka okna zvolena podle minimálně požadovaného rozlišení frekvenční osy spekter, které je dle vztahu ∆f=fvz/N rovno 1 Hz [6]. Vstupní signál je vzorkován kmitočtem 250 Hz. Tato hodnota je při uvažování nejvyšší zájmové harmonické složky rovné 35-ti Hz více než dostačující a byla zvolena v zájmu zachování detailnějšího tvaru EEG signálu pro jeho záznam a pozdější vykreslení. Překrytí okna činí N/10 neboli 25 vzorků. To vše ve výsledku znamená, že desetkrát za vteřinu (∆t=100ms) je vypočítáno a současně vykresleno jednostranné výkonové spektrum z úseku signálu o délce 1s. Obnovovací frekvence 10 Hz by měla být dostatečná k vyvolání zpětnovazebních smyček v neuronové síti lidského mozku a tím i vytvoření efektivního procesu učení a autoregulace.
4 Realizace zpětné vazby Klíčovou podmínkou použitelnosti metody je správné navázání vizuálních prostředků zpětnovazební smyčky na sledovaný Délka signálu v jedné realizaci je shodná se vzorkovacím kmitočtem signálu. 1
74
1,2
60000
1
50000
0,8
40000
0,6
30000
0,4
20000
0,2
10000
0
AMP35 [-]
Adig.filtry [-]
VOL.16, NO.2, APRIL 2014
0 0
5
10
15
20
25
30
35
40
45
50
FIR [-] IIR Delta [-] IIR Theta [-] IIR Alpha [-] IIR Beta [-]
MP35 [-]
f [Hz] Obrázek 2: Amplitudové frekvenční charakteristiky filtrů parametr snímaného signálu. Tento parametr zde představuje okamžitá aktivita příslušných EEG pásem, kterou lze ze signálu získat několika rozdílnými způsoby. V této práci jsou diskutovány dvě níže uvedené metody.
ny pro získání reprezentace okamžitých výkonů pásem. Krátkodobé integrátory určují ekvivalentní časové okno odhadu, což umožňuje sledovat vývoj spektra s časem. Podělíme-li každý z výstupů schématu šířkou přiřazeného pásma, získáme normovaný údaj o výkonu v tomtéž pásmu [W.Hz-1]. Nespornou výhodou tohoto řešení je zlepšení časové odezvy, resp. zkrácení dopravního zpoždění při vyhodnocování aktivity vyšších frekvenčních pásem. Toho se dá docílit jednoduchým zkrácením časového okna příslušejících pásem. Úvaha je taková, že změní-li se např. v pásmu alfa skokově amplituda určité frekvenční složky, objeví se tato změna na výstupu systému nejrychleji právě při délce okna rovnajícího se jedné celé periodě uvažované složky. Při délce okna o hodnotě vyšší bude dána okamžitá úroveň spektrální složky průměrem všech těchto úrovní v daném okně obsažených a tím pádem se na výstupu změna projeví s větším zpožděním. V praxi by tedy pro zjednodušení postačilo vytvořit pro každé pásmo krátkodobý integrátor s takovou konstantou (délkou časového okna) rovnou periodě nejnižší frekvenční složky propuštěné příslušným filtrem. V našem případě, tj. při uvažování dolních mezních frekvencí EEG pásem fm podle kapitoly 3, překrytí oken rovnajícímu se 25-ti vzorkům a vzorkovací frekvenci 250-ti Hz, by vypadaly délky oken N tak, jak udává tabulka 1 (hodnoty jsou vždy zaokrouhleny na nejbližší celý vyšší násobek překrytí oken):
4.1 Průměrování spektrálních čar: Vhodným průměrováním všech výkonových spektrálních čar příslušejících svému EEG pásmu získáme vhodnou aproximaci aktivity v tomto pásmu, kterou pak lze použít k vyhodnocení úspěšnosti tréninku. Přesnost této metody stoupá s rozlišením frekvenční osy spektrogramu a klesá s prosakováním spektra (vliv volby použitého váhového okna). Principiálně tuto metodu vysvětluje obrázek 3. Index n v rovnici (2) zde uvádí pořadí spektrální složky v individuálním výkonovém spektru S.
𝑃𝐷𝐸𝐿𝑇𝐴
1 = 4
𝑃𝐴𝐿𝑃𝐻𝐴 =
1 5
4
𝑆 𝑛 ; 𝑃𝑇𝐻𝐸𝑇𝐴 𝑛=1
1 = 4
13
𝑆(𝑛) ; 𝑃𝐵𝐸𝑇𝐴 = 𝑛=9
1 22
8
𝑆 𝑛 ; 𝑛=5
(2)
35
𝑆(𝑛) 𝑛=14
Tabulka 1: Výsledné délky oken Pásmo DELTA THETA ALPHA SMR BETA Obrázek 3: Metoda průměrování spektrálních čar
Dolní mezní frekvence 0,5 Hz 4 Hz 8 Hz 12 Hz 13 Hz
Délka okna 500 vzorků 100 vzorků 50 vzorků 25 vzorků 25 vzorků
Příklad výpočtu pro pásmo THETA je následující:
4.2 Výpočet výkonu pásem pomocí banky filtrů:
𝑇𝑇𝐻𝐸𝑇𝐴 _𝑚𝑎𝑥 = 1 4 𝐻𝑧 = 0,250 𝑠 → 𝑁𝑇𝐻𝐸𝑇𝐴 = 0,250 𝑠 ∙ 𝑓𝑣𝑧 = = 0,250 𝑠 ∙ 250 = 62,5 𝑣𝑧 ≈ 100 𝑣𝑧
Principiální schéma metody je uvedené na obrázku 4. Vstupní signál je nejprve rozdělen na dílčí frekvenční pásma sadou paralelních filtrů. Takto získané složky jsou následně umocně-
75
VOL.16, NO.2, APRIL 2014
Toto řešení bohužel zavádí do systému problém oscilací všech harmonických složek v pásmu, jejichž perioda není rovna celému násobku délky okna. Kmity se projeví v případě delšího časového úseku signálu konstantní aktivity a rozsah vzniklého vlnění se nakonec ukázal být pro aplikaci EEG Biofeedback nežádoucí. Z toho důvodu byla pro aplikaci vybrána metoda první, čili metoda spektrálních čar.
načítání dat, zpracování dat a řízení zpětné vazby v reálném čase možnost volit kanály zahrnuté do tréninku, parametry tréninku a velikost rozhodovacích úrovní možnost záznamu zvolených kanálů snímaného EEG signálu do souboru a možnost zpětného prohlížení vytváření pacientské databáze zahrnující také informace o jednotlivých sezeních (délka, úspěšnost, trénovaný parametr,…) a vyhodnocení dlouhodobé úspěšnosti terapie jednoduché nastavení napěťového rozsahu a citlivosti jednotlivých vstupních kanálů přehledné a jednoduché uživatelské prostředí. výběr z většího počtu herních motivů možnost vkládání informačních popisků do snímaného průběhu
Obrázek 4: Výpočet výkonu pásem pomocí banky filtrů
5.1 Funkční popis
Výsledky srovnávacích testů obou metod udává obrázek 5. Vstupem testu byl napěťový skok signálů pěti harmonických složek o kmitočtu 3, 6, 10, 13 a 24 Hz. Tento skok by měl simulovat rázovou změnu aktivity (napětí) vybraných složek všech spektrálních pásem z hodnoty 1 uV na hodnotu 10 uV a zpět. Větší strmost náběžných a sestupných hran křivek druhé metody jasně potvrzuje lepší časové odezvy ve vyšších pásmech. Oproti nim je při první metodě dopravní zpoždění všech pásem konstantní, rovnající se dle vztahu ∆t=T*(N/M) jedné vteřině (M zde představuje velikost překrytí časového okna neboli 25 vzorků, N=250 délku okna a T=0,1s periodu vykreslování spektra). Z grafu lze dále vidět, že jediná kmitočtová složka 10 Hz má hladkou odezvu a to právě z důvodu celistvého násobku délky časového okna (50 vz).
Detailní blokové schéma navržené aplikace EEG biofeedback je zobrazeno v příloze 1. Aplikaci lze z funkčního hlediska rozdělit na okno pro nastavení parametrů tréninku, čili okno terapeuta (tvořené programem EEG biofeedback.vi; detail na obrázku 6 – 8), a okno s herním motivem (program Hry.vi; detail na obrázku 9). Účel ovládacích a signalizačních prvků obou oken, stejně jako účel jednotlivých částí blokového schématu, vyplývá z následujícího popisu: Volba režimu - Načítání signálu z měřící karty, načítání s pomocí vstupu zvukové karty (použito pří vývoji aplikace), simulované načítání EEG z externího datového souboru (pro účely výuky a prezentace), zobrazení záznamu EEG a ukončení programu. Záložky - Zobrazení údajů o probíhajícím tréninku, dat pacientské databáze nebo zaznamenaných úseků signálu. Průběh EEG - Vykreslení snímaného EEG signálu po průchodu pásmovou propustí 0,5 Hz až 35 Hz, tvořenou blokem FIR.VI.
1.
2.
5 Aplikace LabVIEW 3. Aplikace pro terapii metodou EEG biofeedback by měla zahrnovat následující funkce [1, 8, 9]: 2 Hz
6 Hz
Zpoždění [s]
0,02 0,01
6
5
4
3
2
1
Metoda průměrování spektrálních čar P [Vrms2.Hz-1]
Metoda odhadu výkonu pásem pomocí banky filtrů P[W.Hz-1]
6
5
0 0
6
5
4
0 3
0
0,03
4
0,01
0,04
3
5
0,05
2
0,02
3 2,5 2 1,5 1 0,5 0 1
0,03
10
2
0
24 Hz 0,04
1
0,02
Zpoždění [s]
13 Hz
0
0,04
Zpoždění [s]
15
Zpoždění [s]
0,06
0
6
5
0
6
5
0
4
0 3
0 2
5
1
0,01 0
5
4
0,02 10
3
10
10 Hz 0,05 15 0,04 0,03 10 0,02 5 0,01 0 0
2
0,03 15
1
15
Zpoždění [s]
Obrázek 5: Porovnání metod odhadu aktivity pásem
76
VOL.16, NO.2, APRIL 2014
1
9
13
3 2 10
12
4
7
5 8 11 6 Obrázek 6: Hlavní okno aplikace EBFB
4.
5.
6.
7.
Průběh pásem - Vykreslení výstupů z jednotlivých filtrů z banky IIR filtrů (blok IIR.VI). Jejich účel je pouze informativní Výkonové spektrum - Ukazatel individuálních výkonových spekter signálu, z nichž se vyhodnocuje aktivita jednotlivých trénovaných pásem a tím i úspěšnost tréninku. Dílčí spektrum se počítá vždy z 250 vzorků signálů uložených v softwarovém posuvném registru, tedy vždy z části signálu o délce jedné vteřiny. Spektrum je vykreslováno v rozsahu 1 – 35 Hz s rozlišením frekvenční osy 1 Hz. Průměrná pásmová aktivita - Ukazatele normované aktivity (výkon vztažený k jednomu hertzi) v jednotlivých pásmech. Jsou pro terapeuta stěžejním ukazatelem distribuce aktivity ve všech pěti základních EEG pásmech (plus dvou přídavných, daných vzájemným poměrem základních pásem [2, 9]). Matematické vyjádření jejich hodnot koresponduje s metodou průměrování spektrálních čar. Volba pásem - Určuje pásma, které budou zahrnuty do tréninkového kola. Jakou formu tréninku terapeut zvolí, to závisí na formě postižení mozkových funkcí pacienta. Obecně lze vybírat mezi relaxačním a výcvikovým tréninkem. Do relaxačních metod patří pásma Delta, Theta a Alfa, naopak do výcvikových metod pásma SMR a Beta. Herní motivy mohou být maximálně dvouparametrické (výsledek hry může být ovlivňován maximálně dvěma parametry). Častokrát je snahou tréninku současně pásmo z jedné skupiny stimulovat, zatímco pásmo ze skupiny opačné potlačit.
Nastavení prahů - Pro každé trénované pásmo nastavuje terapeut hodnotu aktivity, která odlišuje úspěch tréninku od neúspěchu. 9. Výběr hry - Volbou herního tématu se zároveň i spustí odpovídající hra na monitoru pacienta. 10. Tréninkové kolo - Současné spuštění tréninkového kola a záznamu úspěšnosti. Délku kola lze ukončit manuálně nebo s využitím časovače. 11. Úspěšnost tréninku - Slouží pro kontrolu okamžité a průměrné úspěšnosti tréninku. Okamžitá úspěšnost může nabývat stavů z oboru hodnot {0 %, 50 %, 100 %}. Průměrná úspěšnost je vyhodnocována od začátku kola až do jeho konce dle vztahu: 8.
𝑝𝑟ů𝑚ě𝑟𝑛á ú𝑠𝑝ěš𝑛𝑜𝑠𝑡 = 100% ∙
𝑝𝑛 𝑁 𝑛=1 𝑝 ∙𝑁 𝑐
(3)
N…celkový počet parametrů vybraných pro trénink pn…počet cyklů od začátku záznamu, kdy velikost n-tého parametru přesáhla rozhodovací práh pc…počet všech cyklů od začátku záznamu 12. Záznam signálu - Uložení snímaného signál do binárního souboru. Celý zaznamenaný soubor lze zpětně použít jak pro jeho zobrazení ze záznamu (obrázek 8), tak k režimu simulace EBFT (bod 1 kap. 5.1). 13. Vybraný pacient - Zobrazuje jméno a rodné číslo pacienta a umístění souboru s odpovídajícím záznamem EEG. 14. Databáze pacientů - Databáze pro správu pacientů. Uchovávány jsou základní osobní data, diagnóza a parametry jednotlivých terapeutických sezení.
77
VOL.16, NO.2, APRIL 2014
14
Obrázek 7: Databáze pacientů Ty obsahují i hodnoty dosažených výsledků. Z těchto dat je následně pořízen graf dosažené úspěšnosti v závislosti na pořadí tréninkového kola (obrázek 7).
ny ve hře by naopak rušily klidný a relaxovaný stav mysli pacienta. Při relaxační formě tréninku usilujeme o co nejdelší setrvání pásmové aktivity pod inhibičním prahem, respektive nad prahem vybuzení. Potlačují se vyšší pásma (Beta) a naopak se trénuje vzrůst aktivity pásem nižších (Alfa, Theta). Hry se „stimulačním“ soutěžním motivem volíme pro výcvikový model terapie (snaha o vybuzení vyšších pásem). Zde má požadovaný přechod pásmových aktivit přes rozhodovací úrovně charakter impulsů, což může klidně odpovídat skokovým změnám v odvíjejícím se herním motivu. Tyto změny totiž soustředění pacienta nijak nenarušují. Odměňování skokovým vzrůstem skóre naopak vytváří pozitivní zpětnou vazbu a přispívá tak k upevnění správných vzorců funkce pacientova mozku.
Obrázek 8: Prohlížeč záznamu signálu 5.2 Herní témata Celkem byly pro aplikaci vytvořeny 4 herní témata. Jejich popis je uveden v tabulce 2. Obecně lze volit pro každý tréninkový program libovolný herní motiv. Aplikace není v tomto ohledu nijak limitována. Z velmi zobecněné interpretace výsledků [5, 8, 9, 11] však vyplývá, že je dobré volit hry s „udržovacím“ motivem pro relaxační formy tréninku. Cílem je vyvolat a udržet stav uvolnění. Herní motiv na obrazovce musí proto plynout spojitě a pomalu. Rychlé a nečekané změ-
Obrázek 9: Herní motivy
78
VOL.16, NO.2, APRIL 2014 Tabulka 2: Popis herních motivů a variací tréninku Název hry KRAJINA (1 parametr) KVĚTINA (1 parametr) AUTA (1 parametr) PUZZLE (1 parametr)
Parametry 1. parametr Doporučení: 1. parametr Doporučení: 1. parametr Doporučení: 1. parametr Doporučení: 1. parametr
KRAJINA (2 parametry)
2. parametr Doporučení:
KVĚTINA (2 parametry)
1. parametr 2. parametr Doporučení: 1. parametr
AUTA (2 parametry)
2. parametr Doporučení:
PUZZLE (2 parametry)
1. parametr 2. parametr Doporučení:
Možné přidělené pásmo δ, θ, α, β, SMR, SMR/β, α/θ
Odměňování ve hře posun obrazu ve vodorovném i svislém směru ke středu obrazovky Vhodné pro relaxační model tréninku - aktivace nižších pásem či inhibice vyšších pásem; dlouhá tréninková kola δ, θ, α, β, SMR, SMR/β, α/θ rozkvět poupěte, růst stonku a zbarvení květu Vhodné pro relaxační model tréninku - aktivace nižších pásem či inhibice vyšších pásem; kratší tréninková kola δ, θ, α, β, SMR, SMR/β, α/θ předjetí červeného vozu Vhodné pro výcvikový model tréninku - aktivace vyšších pásem; středně dlouhá tréninková kola δ, θ, α, β, SMR, SMR/β, α/θ odkrytí dalšího dílku obrázku Vhodné pro výcvikový model tréninku - aktivace vyšších pásem; krátká tréninková kola α, β, SMR, SMR/β posun obrazu ve vodorovném směru ke středu obrazovky δ, θ, α/θ posun obrazu ve svislém směru ke středu obrazovky Vhodné pro relaxační model tréninku - stimulace nižších pásem za soustavné inhibice vyšších pásem; dlouhá tréninková kola α, β, SMR, SMR/β rozkvět poupěte a růst stonku δ, θ, α/θ zbarvení květu Vhodné pro relaxační model tréninku - stimulace nižších pásem za soustavné inhibice vyšších pásem; krátká tréninková kola α, β, SMR, SMR/β předjetí červeného vozu (podmíněno dopředným pohybem obou vozů) δ, θ, α/θ dopředný pohyb obou vozů Vhodné pro výcvikový model tréninku - stimulace vyšších pásem za soustavné inhibice nižších pásem; středně dlouhá tréninková kola α, β, SMR, SMR/β odkrytí dalšího dílku obrázku δ, θ, α/θ odkrytí dalšího dílku obrázku Vhodné pro výcvikový model tréninku - stimulace vyšších pásem za soustavné inhibice nižších pásem; krátká tréninková kola
Pokaždé, když se pacient dokáže natolik zkoncentrovat, že průměrná pásmová aktivita v hlídaném pásmu skokově překročí požadovanou mez, dojde na jeho monitoru k předjetí automobilu nebo odkrytí dalšího pole skládanky. Z výše uvedených důvodů byly hry KRAJINA a KVĚTINA (obrázek 9 nahoře) navrženy spíše pro pásma Delta, Theta, Alpha a hry AUTA a PUZZLE (obrázek 9 dole) pro pásma Beta a SMR. Na základě uvedených zdrojů jsou také navržena všechna doporučení uvedená v tabulce 2.
Program umožňuje práci jen s jedním snímacím kanálem. Vícekanálové zpracování by znamenalo znásobení nároků na výpočetní výkon. Pro základní trénink je však jednokanálové EEG zcela dostačující [8]. Lze jím provádět jak trénink individuálních EEG pásem, tak i jejich kombinací, čímž je zajištěna vysoká míra variability tréninku. Herní témata pro vizuální zpětnou vazbu úspěšnosti tréninku jsou oblastí aplikace, ve které se projevily negativní důsledky volby programovacího prostředí LabVIEW, které není navrženo pro profesionální práci s herní grafikou (alespoň v jeho základní verzi). Graficky jednodušší herní motivy udrží pacientovu pozornost po kratší dobu, a proto je třeba je i v rámci jednoho tréninkového kola častěji střídat [2]. Tato možnost byla v aplikaci EBFB ponechána. Zvýšení atraktivity her a tím i efektivity tréninku by mohlo být dosaženo doplněním jednotlivých herních motivů o zvukové podněty, které velkou mírou přispívají k upevnění autoregulačních vazeb mozkové aktivity [10, 5]. Aplikace je ve zjednodušené verzi využívána jako demonstrační program při výuce a laboratorních cvičení na FEKT VUT v Brně [7].
6 Závěr Zdrojový kód hotové aplikace realizovaný pomocí LabVIEW umožňuje názornou, leč plně funkční ukázku tréninku za pomoci terapeutické metody EEG biofeedback. Nabízí plnohodnotnou obsluhu stanoviště terapeuta a současně v sobě implementuje čtyři jednoduché zpětnovazební herní motivy. Součástí je i funkční pacientská databáze.
79
VOL.16, NO.2, APRIL 2014
Literatura
[7]
CMIEL, V., O. JANOUSEK a J. KOLAROVA. EEG biofeedback. Proceedings of the 4th International Symposium on Applied Sciences in Biomedical and Communication Technologies - ISABEL '11. New York, New York, USA: ACM Press, 2011, s. 1-5. DOI: 10.1145/2093698.2093752. Dostupné z: http://dl.acm.org/citation.cfm?doid=2093698.2093752
[1]
TYL Jiří; TYLOVÁ, Vendula. Lehké mozkové dysfunkce [Online]. Praha : Asociace pro aplikovanou psychofiziologii a biofeedback ČR - Biofeedback institut, 2002. 3, str. 22. [citováno 24. března 2009.] Dostupný z: http://zsvojanova.cz/doc/LMD.pdf
[2]
TYL, Jiří. Závěrečná zpráva o řešení grantu Interní grantové agentury Ministerstva zdravotnictví ČR.. [Online] Biofeedback institut, 14. červen 2003. [citováno 8. prosince 2008.] Dostupný z: http://www.eegbiofeedback.cz
[8]
KOPŘIVOVÁ, Jana, Martin BRUNOVSKÝ, Ján PRAŠKO a Jiří HORÁČEK. EEG biofeedback a jeho použití v klinické praxi: EEG Biofeedback and its Use in Clinical Practice. Psychiatrie. Praha: Grada, 2008, roč. 12, č. 1, s. 10-16.
[3]
Auto Power Spectrum.VI. National Instruments [online]. National Instruments Corporation, © 2014 [cit. 2014-02-16]. Dostupné z: http://zone.ni.com/reference/en-XX/help/371361H01/lvanls/auto_power_spectrum/
[9]
FRIEL, Patrick. EEG Biofeedback in the Treatment of Atttention Deficit / Hyperactivity Disorder. Alternative Medicine Review. 2007, roč. 12, č. 2, s. 146-151.
[10]
FREDERICK, Jon A., Joel F. LUBAR, Howard W. RASEY, Sheryl A. BRIM a Jared BLACKBURN. Effects of 18.5 Hz Auditory and Visual Stimulation on EEG Amplitude at the Vertex. Journal of Neurotherapy. 1999, vol. 3, 3-4, s. 23-28. DOI: http://dx.doi.org/10.1300/j184v03n03_03.
[11]
SOKHADZE, Estate “Tato” a Olga M. BAZANOVA. Peak Performance Training Using Prefrontal EEG Biofeedback. Biofeedback. 2012, vol. 40, issue 1, s. 71-93. DOI: http://dx.doi.org/10.1002/9781119959984.ch5.
[4]
BIOPAC Systems, Inc. BSL hardware guide. BIOPAC Systems, Inc. [Online] 23. červen 2004. [citováno 24. března 2009.] Dostupný z: http://www.biopac.com
[5]
LINDEN, Michael, Thomas HABIB a Vesna RADOJEVIC. A controlled study of the effects of EEG biofeedback on cognition and behavior of children with attention deficit disorder and learning disabilities. Biofeedback and Self-Regulation. 1996, vol. 21, issue 3, s. 297-297. DOI: http://dx.doi.org/10.1007/bf02214740.
[6]
Discrete Fourier Transform (DFT). The MathWorks, Inc. [online]. © 1994-2014 [cit. 2014-02-16]. Dostupné z: http://www.mathworks.com/help/matlab/math/discretefourier-transform-dft.html
80
1
81
Načtení záznamu
Načítání ze souboru (režim simulace)
1
1
Načítání ze zvukové karty
Načítání z měřicí karty
1
Přepínač režimů
1
Vykreslení záznamu
2
FIR.VI (Filtrace PP 0.5 Hz – 35 Hz)
Řídicí jednotka: výběr režimu výběr trénovaných pásem nastavení rozhodovacích prahů přepínač INHIBICE/STIMULACE volba délky tréninkového kola volba herního motivu volba pacienta editace pacientské databáze záznam REC/PAUSE/STOP
Vykreslení EEG pásem δ, θ, α, β
6
5
4
Výpočet a vykreslení výkonu jednotl. EEG pásem
Úprava a vykreslení ind.výkonového spektra
7, 8, 9, 10, 12, 13
Záznam načítaného signálu (ukládání do souboru)
Výpočet ind.výkonového spektra
IIR.VI (Banka filtrů pro EEG pásma δ, θ, α, β)
Vykreslení snímaného EEG signálu
3
EEG biofeedback.VI
11
14
PRUMERNA USPESNOST.VI (Výpočet okamžité a průměrné úspěšnosti)
DATA.TXT (Pacientská databáze)
Instrukce
Data
Vykreslení okamžité úspěšnosti ve formě hry na pacientově monitoru
Hry.VI
VOL.16, NO.2, APRIL 2014
Příloha 1: Blokové schéma softwarového řešení aplikace EBFB