ˇ e vysok´e uˇcen´ı technick´e v Praze Cesk´ Fakulta elektrotechnick´a Katedra poˇc´ıtaˇcov´e grafiky a interakce
Bakal´aˇrsk´a pr´ace
Vyhodnocov´ an´ı biometrick´ ych mˇ eˇ ren´ı laboratoˇ re ch˚ uze Petr Stolaˇr
Vedouc´ı pr´ace: Ing. Jan Buri´anek
Studijn´ı program: Softwarov´e technologie a management, Bakal´aˇrsk´ y Obor: Web a multimedia 27. kvˇetna 2010
iv
v
Podˇ ekov´ an´ı Chtˇel bych pˇredevˇs´ım podˇekovat sv´emu vedouc´ımu pr´ace Ing. Janu Buri´ankovi za nab´ıdnut´ı moˇznosti zpracov´ an´ı tohoto t´ematu, za jeho konstruktivn´ı kritiku, rady a pˇripom´ınky. D´ ale bych chtˇel podˇekovat rodiˇc˚ um a pˇr´atel˚ um za jejich podporu.
vi
vii
Prohl´ aˇ sen´ı Prohlaˇsuji, ˇze jsem pr´ aci vypracoval samostatnˇe a pouˇzil jsem pouze podklady uveden´e v pˇriloˇzen´em seznamu. Nem´am z´ avaˇzn´ y d˚ uvod proti uˇzit´ı tohoto ˇskoln´ıho d´ıla ve smyslu §60 Z´akona ˇc. 121/2000 Sb., o pr´ avu autorsk´em, o pr´ avech souvisej´ıc´ıch s pr´avem autorsk´ ym a o zmˇenˇe nˇekter´ ych z´akon˚ u (autorsk´ y z´ akon).
V Praze, dne 27. 5. 2010
.............................................................
viii
Abstract This document deals with evaluation of biometric data from gait laboratory. It introduces to readers a laboratory environment and hardware and software equipment. It describes how the subject position is determined in capture volume of laboratory and how is the position stored. It analyzes used output formats of captured data and it mentions current possibilites of visualization. The result of this work is functional aplication prototyp, which allows visualization of captured data from gait laboratory.
Abstrakt ˇ Tato pr´ ace se zab´ yv´ a vyhodnocov´an´ım biometrick´ ych dat laboratoˇre ch˚ uze. Cten´ aˇre seznamuje s prostˇred´ım laboratoˇre a s jeho hardwarov´ ym a softwarov´ ym vybaven´ım. Popisuje jak je urˇcena pozice subjektu v sn´ıman´em prostoru laboratoˇre a jak se tato pozice zaznamen´av´ a. Analyzuje pouˇz´ıvan´e v´ ystupn´ı form´aty mˇeˇren´ ych dat a zmiˇ nuje st´avaj´ıc´ı moˇznosti jejich vizualizace. V´ ysledkem t´eto pr´ ace je funkˇcn´ı prototyp aplikace umoˇzn ˇuj´ıc´ı vizualizovat data poˇr´ızen´ a syst´emy v laboratoˇri ch˚ uze.
ix
x
Obsah ´ 1 Uvod 1.1 Motivace zad´ an´ı . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.2 Struktura pr´ ace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1 1 1
2 Anal´ yza pohybu ch˚ uze 2.1 Souˇradn´ y syst´em doln´ı poloviny lidsk´eho tˇela . . . . . . . . . . . . . . 2.1.1 P´ anev . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.1.1.1 Definice p´anevn´ıho technick´eho souˇradn´eho syst´emu . 2.1.1.2 V´ ypoˇcet stˇred˚ u kyˇceln´ıch kloub˚ u (HJC) . . . . . . . . 2.1.1.3 Definice p´anevn´ıho anatomick´eho souˇradn´eho syst´emu 2.1.2 Stˇred kolenn´ıho kloubu (KJC) . . . . . . . . . . . . . . . . . . 2.1.3 Stehno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.1.4 Doln´ı ˇc´ ast doln´ı konˇctiny . . . . . . . . . . . . . . . . . . . . . 2.1.4.1 Holeˇ n . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.1.4.2 Zadn´ı odd´ıl chodidla . . . . . . . . . . . . . . . . . . . 2.1.4.3 Pˇredn´ı odd´ıl chodidla . . . . . . . . . . . . . . . . . . 2.2 V´ ystupn´ı parametry . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.2.1 P´ anev . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.2.2 Kyˇcel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.2.3 Koleno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2.2.4 Chodidlo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . .
3 4 4 4 5 5 6 6 7 9 10 10 13 15 15 16 16
3 Hardwarov´ e vybaven´ı laboratoˇ re ch˚ uze 3.1 Motion Capture . . . . . . . . . . . . . . 3.1.1 Vyuˇzit´ı . . . . . . . . . . . . . . 3.1.2 Optick´e syst´emy . . . . . . . . . 3.2 Tlakov´ a ploˇsina AMTI OR6-7 . . . . . . 3.2.1 Vyuˇzit´ı . . . . . . . . . . . . . . 3.2.2 Popis zaˇr´ızen´ı . . . . . . . . . . . 3.3 Sn´ım´ an´ı EMG . . . . . . . . . . . . . . . 3.3.1 Co je EMG? . . . . . . . . . . . 3.3.2 Vyuˇzit´ı v laboratoˇri ch˚ uze . . . . 3.3.3 Vys´ılaˇc TeleMyo 2400T G2 . . . 3.3.4 Pˇrij´ımaˇc TeleMyo 2400R G2 . . 3.4 Digit´ aln´ı video kamery . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
19 20 20 20 21 21 21 22 22 23 23 23 24
xi
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
. . . . . . . . . . . .
xii
4 Softwarov´ e vybaven´ı laboratoˇ re ch˚ uze 4.1 Vicon Nexus . . . . . . . . . . . . . . 4.2 Vicon Polygon . . . . . . . . . . . . . 4.2.1 Muskuloskelet´ aln´ı syst´em . . . 4.2.2 Tˇeˇziˇstˇe . . . . . . . . . . . . .
OBSAH
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
25 25 27 27 28
5 Postup sn´ıman´ı pohybu v laboratoˇ ri ch˚ uze 5.1 Nasn´ım´ an´ı dat . . . . . . . . . . . . . . . . . . . . . 5.1.1 Pˇr´ıprava subjektu . . . . . . . . . . . . . . . 5.1.2 Sn´ım´ an´ı a zpracov´ an´ı dynamick´eho z´aznamu 5.1.3 Pˇrezkoum´ an´ı a v´ yplˇ n mezer . . . . . . . . . . 5.1.4 Export dat . . . . . . . . . . . . . . . . . . . 5.2 Zpracov´ an´ı vyexportovan´ ych dat . . . . . . . . . . . 5.2.1 Popis uˇzivatelsk´eho rozhran´ı . . . . . . . . . 5.2.1.1 Zobrazen´ı a skryt´ı mnoˇziny kost´ı . . 5.2.1.2 Panel zobrazen´ı (View pane) . . . . 5.2.1.3 Zobrazen´ı grafu (Graph View) . . . 5.3 Zhodnocen´ı aplikac´ı . . . . . . . . . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
. . . . . . . . . . .
29 30 30 31 31 31 32 32 33 34 34 35
6 Zobrazovaˇ c namˇ eˇ ren´ ych dat laboratoˇ re ch˚ uze 6.1 Rozbor v´ ystupn´ıho form´ atu biometrick´ ych dat 6.1.1 Z´ akladn´ı struktura C3D souboru . . . . 6.1.2 Popis C3D souboru . . . . . . . . . . . . 6.2 Knihovna C3DServer SDK . . . . . . . . . . . . 6.3 Implementace . . . . . . . . . . . . . . . . . . . 6.3.1 Z´ aklad aplikace . . . . . . . . . . . . . . 6.3.2 Naˇcten´ı C3D souboru . . . . . . . . . . 6.3.3 Zobrazovac´ı smyˇcka . . . . . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
. . . . . . . .
37 37 38 39 40 41 42 42 43
. . . . . . . .
. . . . . . . .
. . . . . . . .
7 Testov´ an´ı aplikace na re´ aln´ ych datech
45
8 Z´ avˇ er
47
Literatura
49
A Obrazov´ a pˇ r´ıloha
51
B Seznam pouˇ zit´ ych zkratek
55
C Obsah pˇ riloˇ zen´ eho CD
57
Seznam obr´ azk˚ u 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 2.10 2.11 2.12 2.13
Technick´ y souˇradn´ y syst´em p´anve. Zdroj: [1] . . . . . . . . . . . . . Anatomick´ y souˇradn´ y syst´em p´anve. Zdroj: [1] . . . . . . . . . . . . Vlevo: Stˇred kolenn´ıho kloubu, vpravo: anatomick´ y souˇradn´ y syst´em Zdroj: [1] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Struktura chodidla. Zdroj [1]. . . . . . . . . . . . . . . . . . . . . . . Tˇri rigidn´ı segmenty. Zdroj: [1]. . . . . . . . . . . . . . . . . . . . . . Um´ıstˇen´ı marker˚ u na holeni. Zdroj: [1]. . . . . . . . . . . . . . . . . . Anatomick´ y souˇradn´ y syst´em holenˇe. Zdroj: [1]. . . . . . . . . . . . . Zadn´ı ˇc´ ast chodidla. Zdroj: [1]. . . . . . . . . . . . . . . . . . . . . . Anatomick´ y souˇradn´ y syst´em zadn´ıho odd´ılu chodidla. Zdroj: [1]. . . Rozm´ıstˇen´ı marker˚ u na pˇredn´ım odd´ılu chodidla. Zdroj: [1]. . . . . . Anatomick´ y souˇradn´ y syst´em pˇredn´ıho odd´ılu chodidla. Zdroj: [1]. . Blokov´ y diagram parametr˚ u modelu. . . . . . . . . . . . . . . . . . . Definice kinematick´ ych veliˇcin. . . . . . . . . . . . . . . . . . . . . .
3.1 3.2 3.3 3.4
Laboratoˇr ch˚ uze. Zdroj [3]. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 Z´ akladn´ı Vicon MX syst´em s osmi kamerami. Zdroj: [11]. . . . . . . . . . . . 21 Souˇradn´ y syst´em ploˇsiny AMTI. Zdroj: [12]. . . . . . . . . . . . . . . . . . . . 22 Vlevo: vys´ılaˇc TeleMyo 2400T G2, vpravo: pˇrij´ımaˇc TeleMyo 2400R G2. Zdroj:[13] 23
4.1 4.2
Okno aplikace Vicon Nexus. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 Okno aplikace Vicon Polygon s popsan´ ym uˇzivatelsk´ ym rozhran´ım. Zdroj: [10] 28
5.1 5.2 5.3 5.4 5.5
30 31 33 33
5.6
Zv´ yraznˇen´e povinn´e hodnoty antropometrick´ ych parametr˚ u. . . . . . . . . . . Zrekonstruovan´ y z´ aznam. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Okno aplikace Vicon Polygon s v´ıce otevˇren´ ymi dialogy. . . . . . . . . . . . . Data Object View. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Kontextov´ a nab´ıdka v 3D pracovn´ı ploˇse. Zde zobrazena moˇznost pro skryt´ı mnoˇziny kost´ı. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Dialog pro smaz´ an´ı komponenty. . . . . . . . . . . . . . . . . . . . . . . . . .
6.1 6.2 6.3 6.4
Obsah hlaviˇcky vypsan´ y aplikac´ı VBC3DEditor. . . . Sekce parametru vypsan´a aplikac´ı VBC3DEditor. . . . Sekce 3D dat vypsan´ a aplikac´ı VBC3DEditor. . . . . . Sn´ımek obrazovky aplikace GDV pro vizualizace dat. .
39 40 40 42
xiii
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . . . . . . . . stehna. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . .
. . . .
. . . .
. . . .
. . . .
4 6 7 8 8 9 9 11 11 12 12 13 14
34 34
´ U ˚ SEZNAM OBRAZK
xiv
A.1 A.2 A.3 A.4 A.5
Sn´ımek z v´ ysledn´e aplikace. . . . . . . . . . . . . . . . . . . . Sn´ımek z v´ ysledn´e aplikace. . . . . . . . . . . . . . . . . . . . Sn´ımek z v´ ysledn´e aplikace. . . . . . . . . . . . . . . . . . . . Prostˇred´ı brnˇensk´e laboratoˇre ch˚ uze. Zdroj: [2] . . . . . . . . Zdravotn´ı sestra Aneta Fedrov´ a z brnˇensk´e laboratoˇre ch˚ uze nohy pacientky retroreflexn´ı markery. Zdroj: [4] . . . . . . . . A.6 Prov´adˇen´ı dynamick´eho z´ aznamu. Zdroj: [4] . . . . . . . . . . A.7 Prov´adˇen´ı dynamick´eho z´ aznamu. Zdroj: [4] . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . pˇripevˇ nuje na . . . . . . . . . . . . . . . . . . . . . . . . . . .
51 52 52 53
C.1 Seznam pˇriloˇzen´eho CD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
57
53 54 54
Seznam tabulek 2.1 2.2 2.3 2.4 2.5 2.6
Technick´ y souˇradn´ y syst´em p´anve . . . . . . . . . . . . . Anatomick´ y souˇradn´ y syst´em p´anve . . . . . . . . . . . Anatomick´ y souˇradn´ y syst´em stehna . . . . . . . . . . . Anatomick´ y souˇradn´ y syst´em holenˇe . . . . . . . . . . . Anatomick´ y souˇradn´ y syst´em zadn´ıho odd´ılu chodidla . Anatomick´ y souˇradn´ y syst´em pˇredn´ıho odd´ılu chodidla .
. . . . . .
5 6 7 10 11 12
6.1
Z´ akladn´ı struktura C3D souboru. . . . . . . . . . . . . . . . . . . . . . . . . .
39
xv
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
xvi
SEZNAM TABULEK
Kapitola 1
´ Uvod 1.1
Motivace zad´ an´ı
Pro potˇreby t´eto pr´ ace bylo vyuˇzito prostˇred´ı laboratoˇre ch˚ uze ve Fakultn´ı dˇetsk´e nemocnici J. G. Mendela v Brnˇe (viz. poˇrad PORT [17]). Tato laboratoˇr je souˇc´ast´ı Kliniky dˇetsk´e chirurgie, ortopedie a traumatologie. Vyuˇz´ıv´a se k upˇresnˇen´ı indikac´ı chirurgick´e l´eˇcby pohybov´eho apar´ atu dˇet´ı postiˇzen´ ych dˇetskou mozkovou obrnou. Laboratoˇr vybaven´ a nejnovˇejˇs´ımi technologiemi, poskytuje l´ekaˇr˚ um flexibiln´ı, rychl´e a dostateˇcnˇe pˇresn´e prostˇredky k urˇcen´ı diagn´oz a postupu l´eˇcby. Kombinuj´ı se zde optick´e sn´ımaˇce pohybu, tlakov´e ploˇsiny a EMG syst´emy. Takov´ato zaˇr´ızen´ı b´ yvaj´ı ˇcasto technicky velmi sloˇzit´ a a kladou na obsluhuj´ıc´ı person´al znaˇcn´e n´aroky. Jelikoˇz se jedn´a pˇrev´aˇznˇe o lidi vzdˇelan´e v medic´ınsk´em oboru, je d˚ uleˇzit´e, aby software ˇr´ıd´ıc´ı tyto syst´emy poskytoval snadn´e ovl´ ad´ an´ı. Extr´emnˇe d˚ uleˇzit´a je v´ ysledn´a vizualizace namˇeˇren´ ych biometrick´ ych dat. Mus´ı striktnˇe odpov´ıdat potˇreb´ am diagnostikuj´ıc´ıch l´ekaˇr˚ u, kteˇr´ı na z´akladˇe tˇechto dat urˇc´ı n´asleduj´ıc´ı postup l´eˇcby. C´ılem t´eto pr´ ace je analyzovat st´avaj´ıc´ı pouˇz´ıvan´e syst´emy umoˇzn ˇuj´ıc´ı z´aznam pacientovi ch˚ uze a v´ yslednou vizualizaci. Pro spr´avn´e pochopen´ı vˇsech souvislost´ı je zapotˇreb´ı nastudovat celkovou problematiku lidsk´e ch˚ uze a tak´e form´aty dat, kter´e slouˇz´ı k uchov´ an´ı nasn´ıman´e informace. K tomu poslouˇzili dokumenty od R. B. Davise a spol. [5] a M. P. Kadaby a spol. [6]. Na z´ akladˇe t´eto anal´ yzy byl implementov´an vizualizaˇcn´ı syst´em, jehoˇz vstupem jsou nasn´ıman´ a data z optick´ ych sn´ımaˇc˚ u.
1.2
Struktura pr´ ace
Struktura bakal´ aˇrsk´e pr´ ace je takov´ato. V n´asleduj´ıc´ı druh´e kapitole je pops´ano, jak je definov´ano rozm´ıstˇen´ı retroreflexn´ıch marker˚ u na tˇele pacienta. D´ale jak se urˇc´ı v´ ysledn´a poloha jednotliv´ ych ˇc´ ast´ı tˇela v 3D prostoru laboratoˇre. Kapitola 3 se zab´ yv´a hardwarov´ ym vybaven´ı laboratoˇre. N´ asleduje kapitola 4. V n´ı je zm´ınka o softwarov´em vybaven´ı laboratoˇre. Zde jsou pops´ any programy Vicon Nexus a Vicon Polygon jako kl´ıˇcov´e prvky pro z´aznam a anal´ yzu pacientovi ch˚ uze. V kapitole 5 je pops´an cel´ y postup pr´ace vedouc´ı k nasn´ım´an´ı
1
2
´ KAPITOLA 1. UVOD
dynamick´eho z´ aznamu pomoc´ı programu Vicon Nexus. Je zde tak´e nast´ınˇena pr´ace s programem Vicon Polygon, kter´ y umoˇzn ˇuje v´ yslednou vizualizaci vˇsech dat. V kapitole 6 pˇrich´az´ı na ˇradu rozbor zobrazovaˇce namˇeˇren´ ych dat spoleˇcnˇe s popisem v´ ystupn´ıch form´at˚ u a v´ yslednou implementac´ı. V kapitole 7 je zm´ınˇeno testov´an´ı aplikace na re´aln´ ych datech. Pot´e n´asleduje z´ avˇer, ve kter´em je shrnut v´ ysledek pr´ace.
Kapitola 2
Anal´ yza pohybu ch˚ uze Na zaˇc´ atku tohoto dokumentu je zapotˇreb´ı nast´ınit nutnou teorii, na kter´e je cel´a problematika postavena. Anal´ yza pohybu, jak ji dnes zn´ame, se zaˇcala vyv´ıjet na poˇc´atku osmdes´at´ ych let dvac´ at´eho stolet´ı. V˚ ubec prvn´ı laboratoˇr´ı, kter´a se zab´ yvala myˇslenkou vyuˇzit´ı sn´ım´an´ı pohybu v medic´ınˇe, byla laboratoˇr v Dˇetsk´e nemocnici Newington (Hartford, Connecticut, USA). Nˇekolika ˇclenn´e t´ ymy pracovali na spoleˇcn´em c´ıli. V´ ysledkem toho vˇseho bylo vyvinut´ı biomechanick´eho modelu, kter´ y popsali M. P. Kadaba a spol. v textu [6] a R. B. Davis a spol. v textu [5].
Na z´ akladˇe tohoto biomechanick´eho modelu, naz´ yvan´eho tak´e jako konvenˇcn´ı model ch˚ uze, vznikl matematick´ y model, kter´ y umoˇzn ˇuje sledovat chov´an´ı syst´emu ovlivnˇen´eho urˇcit´ ymi vstupy a parametry. Model funguje tak, ˇze z trajektorie fyzick´ ych marker˚ u dopoˇc´ıt´a trajektorii virtu´ aln´ıch marker˚ u, kter´e slouˇz´ı k urˇcen´ı kinematick´e a kinetick´e veliˇciny (´ uhly, momenty, atd.) a reprezentaci segment˚ u tˇela. Pˇred samotn´ ym mˇeˇren´ım je tˇreba pomˇeˇrit antropometrick´e parametry pacienta (hmotnost a v´ yˇska osoby, d´elka nohy, ˇs´ıˇrka kolena a kotn´ıku). D´ale se definuj´ı rigidn´ı segmenty. Kaˇzd´ y segment je pops´ an sv´ ym ortogon´aln´ım souˇradn´ ym syst´em vyuˇz´ıvaj´ıc´ım pravidla prav´e ruky. Model ch˚ uze je zaloˇzen na n´asleduj´ıc´ıch pˇredpokladech: • hypot´eza rigidn´ıho tˇela – ˇc´ asti lidsk´eho tˇela jsou povaˇzov´any za rigidn´ı a nedeformovateln´e. • sn´ım´ an´ı rigidn´ıch segment˚ u – fyzick´e markery jsou pˇripevnˇeny na povrchu tˇela osoby; virtu´aln´ı markery jsou generov´ any z fyzick´ ych marker˚ u a z fyzick´ ych parametr˚ u osoby. Z v´ yˇse uveden´eho je patrn´e, ˇze kvalita, pˇresnost a konzistence pozice fyzick´ ych a virtu´ aln´ıch marker˚ u m´ a velk´ y vliv na v´ ysledn´ y v´ ystup.
3
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
4
2.1
Souˇ radn´ y syst´ em doln´ı poloviny lidsk´ eho tˇ ela
Jak bylo v´ yˇse zm´ınˇeno, kaˇzd´ a ˇc´ ast doln´ı poloviny lidsk´eho tˇela m´a sv˚ uj vlastn´ı ortogon´aln´ı souˇradn´ y syst´em. K urˇcen´ı tˇechto syst´em˚ u n´am slouˇz´ı z´ıskan´e pozice retroreflexn´ıch marker˚ u. Tyto fyzick´e markery n´ am umoˇzn ˇuj´ı vytvoˇrit tzv. technick´ y souˇradn´ y syst´em, kter´ y s pomoc´ı virtu´aln´ıch marker˚ u umoˇzn´ı vytvoˇren´ı anatomick´eho souˇradn´eho syst´emu segmentu. D˚ uleˇzit´ y je smˇer ch˚ uze subjektu v glob´ aln´ım souˇradn´em syst´emu. Vypoˇc´ıt´a se z pˇredchoz´ı a n´asledn´e pozice LASI markeru. Pˇritom se porovn´ av´a posunut´ı po ose X a Y. Podle toho kter´e posunut´ı m´a vˇetˇs´ı hodnotu, subjekt se pohybuje po t´e ˇci on´e ose.
2.1.1
P´ anev
Vytvoˇren´ı p´anevn´ıho anatomick´eho souˇradn´eho syst´emu ovlivˇ nuj´ı n´asleduj´ıc´ı kroky: 2.1.1.1
Definice p´ anevn´ıho technick´ eho souˇ radn´ eho syst´ emu
P´anevn´ı technick´ y souˇradn´ y syst´em 2.1 je urˇcen markery pˇripevnˇen´ ymi na p´anev subjektu. Markery maj´ı oznaˇcen´ı LASI, RASI, LPSI, RPSI nebo alternativnˇe LASI, RASI, SACR. Z´aleˇz´ı na um´ıstˇen´ı marker˚ u na p´ anvi. Kdyˇz se pouˇzij´ı PSI markery, vytvoˇr´ı se mezi nimi virtu´aln´ı marker pouˇz´ıvan´ y pro kalkulaci. Jinak kdyˇz je pouˇzit SACR marker, bude pouˇzit pro kalkulaci on.
Obr´ azek 2.1: Technick´ y souˇradn´ y syst´em p´anve. Zdroj: [1]
ˇ ´ SYSTEM ´ DOLN´I POLOVINY LIDSKEHO ´ ˇ 2.1. SOURADN Y TELA
Co Poˇc´atek 1. definovan´ a osa 2. definovan´ a osa
Y Z
3. definovan´ a osa
X
Jak (LASI + RASI)/2 RASI → LASI Kolm´ a k rovinˇe definovan´e markery LASI, RASI, LPSI, RPSI nebo LASI, RASI, SACR Vektorov´ y souˇcin mezi jednotkov´ ym vektorem Y a Z
5
Popis Horizont´aln´ı, zprava doleva Vertik´aln´ı, ukazuj´ıc´ı vzh˚ uru
Ukazuj´ıc´ı dopˇredu
Tabulka 2.1: Technick´ y souˇradn´ y syst´em p´anve 2.1.1.2
V´ ypoˇ cet stˇ red˚ u kyˇ celn´ıch kloub˚ u (HJC)
Problematice nalezen´ı stˇredu kyˇceln´ıch kloub˚ u v p´anevn´ım technick´em souˇradn´em syst´emu, se vˇenovali Davis a spol. v dokumentu [5]. Vyuˇz´ıv´a se regresivn´ıho pˇr´ıstupu k nalezen´ı vektoru spojuj´ıc´ıho p˚ uvodn´ı p´ anevn´ı technick´ y souˇradn´ y syst´em s HJC. Regrese je pˇritom zaloˇzena na vzd´ alenosti mezi anterior superior iliaca spina (ASIS - trn kyˇceln´ı pˇredn´ı, horn´ı) a d´elkou nohy. Lok´ aln´ı souˇradnice stˇredu lev´eho kyˇceln´ıho kloubu v p´anevn´ım technick´e souˇradn´e syst´emu jsou (pro prav´ y HJC je potˇreba invertovat Y souˇradnici): LHJCx = C * cos(θ) * sin(β) - (ASIS-TrochanterDistance + mm) * cos(β) LHJCy = -(C * sin(θ) - aa) LHJCz = -C * cos(θ) * cos(β) - (ASIS-TrochanterDistance + mm) * sin(β) Kde: θ = 0.5 rad β = 0.314 rad ASIS-TrochanterDistance = 0.1288 * delkaNohy - 48.56 C = 0.115 * delkaNohy - 15.3 aa = (vzdalenost(LASI,RASI))/2 mm = polomˇer markeru 2.1.1.3
Definice p´ anevn´ıho anatomick´ eho souˇ radn´ eho syst´ emu
Jakmile je lokalizov´ an HJC a urˇcen p´anevn´ı technick´ y souˇradn´ y syst´em, m˚ uˇze b´ yt definov´ an anatomick´ y p´ anevn´ı souˇradn´ y syst´em. Smˇer p´anevn´ıho anatomick´eho souˇradn´eho syst´emu je stejn´ y jako technick´ y souˇradn´ y syst´em, ale poˇc´atek je posunut do stˇredu LHJC a RHJC definovan´ ych v 2.1.1.2.
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
6
Obr´ azek 2.2: Anatomick´ y souˇradn´ y syst´em p´anve. Zdroj: [1] Co Poˇc´atek 1. definovan´a osa 2. definovan´a osa
Y Z
3. definovan´a osa
X
Jak (LHJC + RHJC)/2 RASI → LASI Kolm´ a k rovinˇe definovan´e markery LASI, RASI, LPSI, RPSI nebo LASI, RASI, SACR Vektorov´ y souˇcin mezi jednotkov´ ym vektorem Y a Z
Popis Horizont´aln´ı, zprava doleva Vertik´aln´ı, ukazuj´ıc´ı vzh˚ uru
Ukazuj´ıc´ı dopˇredu
Tabulka 2.2: Anatomick´ y souˇradn´ y syst´em p´anve
2.1.2
Stˇ red kolenn´ıho kloubu (KJC)
Vypoˇc´ıtat stˇred kolenn´ıho kloubu lze pomoc´ı dvou metod. Prvn´ı pro urˇcen´ı KJC vyuˇz´ıv´a speci´aln´ı kolenn´ı marker - Knee Aligment Device (KAD). V brnˇensk´e laboratoˇri ch˚ uze se nepouˇz´ıv´a, proto nebude d´ ale zmiˇ nov´ an. Druh´a metoda urˇcuje pozici KJC d´ıky znalosti glob´aln´ı pozice HJC, pozice stehenn´ıho markeru THI, kolenn´ıho markeru KNE a offsetu kolena KO (KO = (ˇs´ıˇrka kolena + pr˚ umˇer markeru)/2). KJC je ve vzd´ alenosti KO od KNE, kolmo na rovinu tvoˇrenou markery KNE,THI a HJC. V´ıce obr. 2.3 vlevo.
2.1.3
Stehno
Anatomick´ y souˇradn´ y syst´em stehna vyuˇz´ıv´a pouze jeden fyzick´ y marker THI. K popisu rigidn´ı ˇc´asti tˇela v prostoru potˇrebujeme alespoˇ n tˇri body. Vyuˇzijeme proto poznatku o um´ıstˇen´ı KJC definovan´eho v 2.1.2 a HJC definovan´eho v 2.1.1.2. V´ıce obr. 2.3 vpravo.
ˇ ´ SYSTEM ´ DOLN´I POLOVINY LIDSKEHO ´ ˇ 2.1. SOURADN Y TELA
7
Obr´azek 2.3: Vlevo: Stˇred kolenn´ıho kloubu, vpravo: anatomick´ y souˇradn´ y syst´em stehna. Zdroj: [1]
Co Poˇc´atek 1. definovan´ a osa 2. definovan´ a osa
Z X
3. definovan´ a osa
Y
Jak RKJC RKJC → RHJC Kolm´ a k rovinˇe definovan´e RHJC, RKJC, THI Vektorov´ y souˇcin mezi jednotkov´ ym vektorem Z a X
Popis Vertik´aln´ı, ukazuj´ıc´ı vzh˚ uru Ukazuj´ıc´ı dopˇredu Horizont´aln´ı, zprava doleva
Tabulka 2.3: Anatomick´ y souˇradn´ y syst´em stehna
2.1.4
Doln´ı ˇ c´ ast doln´ı konˇ ctiny
Jako kl´ıˇcov´ a se pˇri anal´ yze ch˚ uze jev´ı potˇreba velmi detailn´ıho nasn´ım´an´ı pohyb˚ u pacientovi doln´ı konˇcetiny a hlavnˇe chodidla. To vyˇzaduje um´ıstˇen´ı vˇetˇs´ıho mnoˇzstv´ı marker˚ u na povrch doln´ı konˇcetiny. Samotn´ı pacienti - dˇeti, pˇredstavuj´ı urˇcit´e v´ yzvy. A to pˇredevˇs´ım jejich mal´ y povrch chodidla a vetˇs´ı promˇenlivost pohybu chodidla. Pro detailn´ı anal´ yzu byl vytvoˇren Oxford Foot Model (OFM) vyvinut´ y Nuffieldsk´ ym ortopedick´ ym centrem v Oxfordu, pˇredn´ım svˇetov´ ym centrem pro klinickou a v´ yzkumnou anal´ yzu ch˚ uze, ve spolupr´ aci s Oxfordskou universitou. Samotn´e chodidlo je velice sloˇzit´a struktura 2.4. Chodidlo utv´aˇr´ı: • 26 kost´ı, • 33 kloub˚ u,
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
8
• 107 vaz˚ ua • 19 sval˚ u.
Obr´ azek 2.4: Struktura chodidla. Zdroj [1].
Z v´ yˇse uveden´eho je patrn´e, ˇze k dostateˇcn´emu popisu chodidla nestaˇc´ı pouze jeden segment, ale je zapotˇreb´ı chodidlo rozdˇelit do v´ıce segment˚ u. I v tomto pˇr´ıpadˇe se vyuˇz´ıv´a pˇredpokladu rigidn´ıho tˇela. Doln´ı ˇc´ ast doln´ı konˇcetiny je rozdˇelena do tˇr´ı segment˚ u - holeˇ n, zadn´ı odd´ıl chodidla a pˇredn´ı odd´ıl chodidla, obr. 2.5. Jako doplˇ nkov´a ˇc´ast m˚ uˇze b´ yt do v´ ypoˇctu zahrnut hallux (palec).
Obr´ azek 2.5: Tˇri rigidn´ı segmenty. Zdroj: [1].
ˇ ´ SYSTEM ´ DOLN´I POLOVINY LIDSKEHO ´ ˇ 2.1. SOURADN Y TELA
2.1.4.1
9
Holeˇ n
OFM definuje um´ıstˇen´ı fyzick´ ych marker˚ u na holeni, viz. obr. 2.6. Anatomick´ y souˇradn´ y syst´em holenˇe je definovan´ y stˇredem kotn´ıku (AJC = (ANK + MMA)/2) a virtu´aln´ım markerem KJC, viz. obr. 2.7.
Obr´ azek 2.6: Um´ıstˇen´ı marker˚ u na holeni. Zdroj: [1].
Obr´ azek 2.7: Anatomick´ y souˇradn´ y syst´em holenˇe. Zdroj: [1].
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
10
Co Poˇc´atek 1. definovan´a osa 2. definovan´a osa
Y X
3. definovan´a osa
Z
Jak AJC = (ANK + MMA)/2 AJC → KJC Kolm´ a k rovinˇe definovan´e KJC, ANK, MMA Vektorov´ y souˇcin mezi jednotkov´ ym vektorem Y a X
Popis Vertik´aln´ı, ukazuj´ıc´ı vzh˚ uru Ukazuj´ıc´ı dopˇredu Horizont´aln´ı, zleva doprava
Tabulka 2.4: Anatomick´ y souˇradn´ y syst´em holenˇe
2.1.4.2
Zadn´ı odd´ıl chodidla
Anatomick´ y souˇradn´ y syst´em zadn´ıho odd´ılu chodidla popisuje sagit´aln´ı rovina kosti patn´ı, tvoˇren´a markery HEE, PCA a (LCA + STL)/2, viz obr. 2.8 a obr. 2.9. 2.1.4.3
Pˇ redn´ı odd´ıl chodidla
Pˇredn´ı odd´ıl chodidla je opatˇren markery P1M, P5M, D5M, TOE a D1M, viz. obr. 2.10. Transvers´aln´ı rovina je tvoˇren´ a metatrs´ aln´ımi markery D1M, D5M a P5M. V´ ysledn´ y anatomick´ y souˇradn´ y syst´em zn´ azorˇ nuje obr. 2.11.
ˇ ´ SYSTEM ´ DOLN´I POLOVINY LIDSKEHO ´ ˇ 2.1. SOURADN Y TELA
11
Obr´ azek 2.8: Zadn´ı ˇc´ast chodidla. Zdroj: [1].
Obr´ azek 2.9: Anatomick´ y souˇradn´ y syst´em zadn´ıho odd´ılu chodidla. Zdroj: [1].
Co Poˇc´atek 1. definovan´ a osa 2. definovan´ a osa
Y X
3. definovan´ a osa
Z
Jak HEE HEE → KJC Rovnobˇeˇzn´a s podlahou, HEE → (STL+LCA)/2 Vektorov´ y souˇcin mezi jednotkov´ ym vektorem Z a X
Popis Vertik´aln´ı, ukazuj´ıc´ı vzh˚ uru Ukazuj´ıc´ı dopˇredu Horizont´aln´ı, zleva doprava
Tabulka 2.5: Anatomick´ y souˇradn´ y syst´em zadn´ıho odd´ılu chodidla
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
12
Obr´azek 2.10: Rozm´ıstˇen´ı marker˚ u na pˇredn´ım odd´ılu chodidla. Zdroj: [1].
Obr´azek 2.11: Anatomick´ y souˇradn´ y syst´em pˇredn´ıho odd´ılu chodidla. Zdroj: [1]. Co Poˇc´atek 1. definovan´a osa
X
2. definovan´a osa
Y
3. definovan´a osa
Z
Jak (P5M + P1M)/2 Pˇr´ımka z poˇc´atku do TOE prom´ıtnut´ a do transvers´aln´ı roviny metatars´ aln´ıch ˇcl´ank˚ u Kolm´ a k transvers´aln´ıch rovinˇe metatars´ aln´ıch ˇcl´ank˚ u Vektorov´ y souˇcin mezi jednotkov´ ym vektorem Y a X
Popis Ukazuj´ıc´ı dopˇredu
Vertik´aln´ı, ukazuj´ıc´ı vzh˚ uru Horizont´aln´ı, zleva doprava
Tabulka 2.6: Anatomick´ y souˇradn´ y syst´em pˇredn´ıho odd´ılu chodidla
´ ´I PARAMETRY 2.2. VYSTUPN
2.2
13
V´ ystupn´ı parametry
Je - li zn´ ama pozice jednotliv´ ych marker˚ u a anatomick´e souˇradn´e syst´emy jednotliv´ ych segment˚ u, m˚ uˇzeme je spoleˇcnˇe s antropometrick´ ymi parametry aplikovat na matematick´ y model ch˚ uze, kter´ y generuje na v´ ystupu kloubn´ı kinematiku a kinetiku. Viz. obr. 2.12.
Obr´ azek 2.12: Blokov´ y diagram parametr˚ u modelu.
• Kinematika ´ – Uhly kloub˚ u ∗ Relativn´ı u ´hly mezi dvˇema rigidn´ımi segmenty. Vˇzdy tvoˇren´e flex´ı, abdukc´ı a rotac´ı. ∗ Absolutn´ı u ´hly mezi rigidn´ım segmentem a fixn´ım prostˇred´ım laboratoˇre. Vˇzdy tvoˇren´e rotacemi kolem glob´aln´ıch os X, Y a Z. – Kosti ∗ Mnoˇzina ˇctyˇr virtu´aln´ıch bod˚ u asociovan´ ych s kaˇzd´ ym modelovan´ ym segmentem tˇela. · jmenoSegmentuO - poˇc´atek segmentu jmenoSegmentu · jmenoSegmentuL - virtu´aln´ı bod ve smˇeru pˇr´ıˇcn´e (later´aln´ı) osy segmentu jmenoSegmentu · jmenoSegmentuA - virtu´aln´ı bod ve smˇeru pod´eln´e osy segmentu jmenoSegmentu · jmenoSegmentuP - virtu´aln´ı bod ve smˇeru vertik´aln´ı osy segmentu jmenoSegmentu • Kinetika – S´ıly ∗ P˚ usob´ıc´ı s´ıly vyj´ adˇren´e v lok´aln´ım referenˇcn´ım syst´emu kaˇzd´e rigidn´ı ˇc´ asti tˇela. ∗ Jednotky [N/Kg]. – Momenty ∗ Je - li pˇripojena tlakov´a ploˇsina, kalkuluje se s produkovan´ ymi momenty prodlouˇzen´ı.
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
14
∗ Vˇzdy seˇrazen´e jako flexe, abdukce, rotace. ∗ Jednotky [Nmm/Kg]. – V´ ykony ∗ Vypoˇc´ıtaj´ı se jako skal´ arn´ı souˇcin mezi momentem a u ´hlovou rychlost´ı. ∗ Mohou b´ yt vyj´ adˇreny jako skal´ar nebo jako tˇri separovan´e smˇerov´e komponenty. ∗ Jednotky [W/Kg].
Kinematick´ e promˇ enn´ e Jedn´a se o veliˇciny, kter´e se z´ısk´ avaj´ı bˇehem cyklu mˇeˇren´ı. Viz. obr. 2.13
Obr´ azek 2.13: Definice kinematick´ ych veliˇcin.
V n´asleduj´ıc´ım textu se zmiˇ nuje nˇekolik term´ın˚ u, kter´e je nejprve tˇreba objasnit. Transverz´aln´ı roviny jsou takov´e roviny, kter´e proch´ azej´ı tˇelem v horizont´aln´ım smˇeru, front´aln´ı roviny vych´azej´ı ze stˇredu tˇela smˇerem k hlavˇe, sagit´aln´ı roviny jsou takov´e roviny, kter´e jsou kolm´e na front´aln´ı i transvers´ aln´ı roviny a rozdˇeluj´ı tˇelo na pravou a levou ˇc´ast.
´ ´I PARAMETRY 2.2. VYSTUPN
2.2.1
15
P´ anev
Sklon p´ anve (Pelvic Tilt) Sklon p´ anve se bˇeˇznˇe poˇc´ıt´ a kolem transvers´aln´ı osy laboratoˇre. Bl´ıˇz´ı-li se vˇsak dopˇredn´ y pohyb pacienta sagit´ aln´ı ose laboratoˇre, sklon p´anve se mˇeˇr´ı podle n´ı. Sagit´aln´ı osa p´anve, kter´a leˇz´ı na transvers´ aln´ı rovinˇe p´anve, se prom´ıtne do sagit´aln´ı roviny laboratoˇre. Skon p´anve je urˇcen jako u ´hel sevˇren´ y mezi prom´ıtnutou sagit´aln´ı rovinou a sagit´aln´ı rovinou laboratoˇre. Kladn´ yu ´hel odpov´ıd´ a norm´aln´ı situaci, kdy PSIS je v´ yˇse neˇz ASIS. Rotace p´ anve (Pelvic Rotation) Rotace p´ anve je poˇc´ıt´ ana kolem front´aln´ı osy p´anevn´ıho souˇradn´eho syst´emu. Jedn´a se o u ´hel sevˇren´ y mezi sagit´ aln´ı osou p´ anve a sagit´aln´ı osou laboratoˇre prom´ıtnut´ y do transvers´ aln´ı roviny p´ anve. Obliquita p´ anve (Pelvic Obliquity) Obliquita p´ anve je u ´hel, kter´ y se mˇeˇr´ı mezi transvers´aln´ı osou laboratoˇre a front´aln´ı osou p´anve. Jin´ a neˇz nulov´ a hodnota obliquity p´anve pˇredstavuje situaci, ve kter´e je jedna strana p´anve n´ıˇze neˇz druh´ a. V´ ysledn´ a pozice p´ anve V´ ysledn´ a rotace p´ anve je urˇcena sklonem p´anve kolem transvers´aln´ı osy, hodnotu p´anevn´ı obliquity kolem jej´ı sagit´ aln´ı osy a rotac´ı kolem sv´e front´aln´ı osy. P´anev je pak v poloze urˇcen´e tˇemito tˇremi u ´hly - sklon, rotace, obliquita p´anve.
2.2.2
Kyˇ cel
Kyˇ celn´ı flexe/extenze (Hip Flexion/Extension) Kyˇceln´ı flexe je poˇc´ıt´ ana jako osa kolm´a k transvers´aln´ı ose p´anve, kter´a proch´az´ı stˇredem kyˇceln´ıho kloubu. Sagit´ aln´ı osa stehna je prom´ıtnuta na rovinu kolmou k ose kyˇceln´ı flexe. Kyˇceln´ı flexe je pak u ´hel mezi prom´ıtnutou sagit´aln´ı osu stehna a sagit´aln´ı osou p´anve. Kladn´ yu ´hel (flexe) odpov´ıd´ a situaci, kdy koleno je pˇred tˇelem. Kyˇ celn´ı rotace (Hip Rotation) Kyˇceln´ı rotace se mˇeˇr´ı kolem pod´eln´e osy stehna a je poˇc´ıt´ana mezi sagit´aln´ı osou stehna a sagit´aln´ı osou p´ anve, prom´ıtnut´e do roviny kolm´e na pod´elnou osu stehna. Kladn´a kyˇceln´ı rotace odpov´ıd´ a vnitˇrn´ı rotaci stehna. Kyˇ celn´ı abdukce/addukce (Hip Ab/Adduction) ´ Kyˇceln´ı addukce se zjist´ı v rovinˇe osy kyˇceln´ı flexe a stˇredu kolenn´ıho kloubu. Uhel je poˇc´ıt´ an mezi pod´elnou osou stehna a front´aln´ı osou p´anve, prom´ıtnut´e do t´eto roviny. Kladn´a hodnota odpov´ıd´ a otoˇcen´ı nohy smˇerem dovnitˇr (addukce).
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
16
V´ ysledn´ a pozice kyˇ cle Stehno jehoˇz pod´eln´ a osa byla rovnobˇeˇzn´a s front´aln´ı osou p´anve, a ve kter´e osa flexe kolena byla rovnobˇeˇzn´ a s transvers´ aln´ı osou p´ anve, by mˇela b´ yt v neutr´aln´ı poloze (vˇsechny 3 u ´hly jsou nulov´e). Zmˇena z t´eto neutr´ aln´ı pozice do aktu´aln´ı pozice stehna je urˇcen´a tˇremi u ´hly. Jedn´a se o u ´hel kyˇceln´ı flexe, u ´hle kyˇceln´ı addukce a u ´hel rotace stehna kolem pod´eln´e osy stehna.
2.2.3
Koleno
Kolenn´ı flexe/extenze (Knee Flexion/Extension) Sagit´aln´ı osa holenˇe se prom´ıtne do roviny kolm´e k ose flexe kolena. Flexe kolena je u ´hel mezi t´ımto prom´ıtnut´ım a sagit´ aln´ı osou stehna. Kladn´a hodnota u ´hlu odpov´ıd´a flexi kolena. Rotace kolena (Knee Rotation) Rotace kolena se mˇeˇr´ı kolem pod´eln´e osy holenˇe. Jedn´a se o u ´hel mezi sagit´aln´ı osou holenˇe a sagit´aln´ı osou stehna, prom´ıtnutou do roviny kolm´e k pod´eln´e ose holenˇe. Kladn´ y u ´hel odpov´ıd´a vnitˇrn´ı rotaci. Valgus/varus kolena (Valgus/Varus) Jedn´a se o u ´hel mezi pod´elnou osou holenˇe a pod´elnou osou stehna prom´ıtnutou do roviny osy kolenn´ı flexe a stˇredu kotn´ıku. Kladn´e hodnotˇe odpov´ıd´a varus (zevn´ı otoˇcen´ı kolena). V´ ysledn´ a pozice kolena Neutr´aln´ı pozice holenˇe je takov´ a, ˇze holeˇ n je ve stejn´e linii jako pod´eln´a osa stehna a osa flexe kotn´ıku je rovnobˇeˇzn´ a s osou flexe kolena. Z t´eto pozice se pomoc´ı u ´hl˚ u kolenn´ı flexe, otoˇcen´ım o hodnotu valgus / varus a u ´hlu rotace kolena, dostane do aktu´aln´ı pozice popsan´e tˇemito tˇremi u ´hly.
2.2.4
Chodidlo
Progrese chodidla (Foot Progression) Progrese chodidla je u ´hel mezi vektorem chodidla (prom´ıtnut´ y do transvers´aln´ı roviny laboratoˇre) a sagit´ aln´ı osou laboratoˇre. Kladn´a hodnota odpov´ıd´a vnitˇrn´ı rotaci chodidla. Dorsi / plant´ arn´ı flexe kotn´ıku (Ankle Dorsi/Plantar Flexion) ´ Vektor chodidla je prom´ıtnut´ y do sagit´ aln´ı roviny chodidla. Uhel mezi vektorem chodidla a sagit´aln´ı osou holenˇe se naz´ yv´ a dorsi / plant´arn´ı flexe chodidla. Kladn´a hodnota u ´hlu odpov´ıd´a dorsiflexi.
´ ´I PARAMETRY 2.2. VYSTUPN
17
Rotace chodidla (Foot Rotation) Mˇeˇr´ı se kolem osy kolm´e k vektoru chodidla a ose flexe kotn´ıku. Jedn´a se o u ´hel mezi vektorem chodidla a sagit´ aln´ı osou holenˇe, prom´ıtnutou do transvers´aln´ı roviny chodidla. Kladn´ a hodnota odpov´ıd´ a vnitˇrn´ı rotaci.
18
´ ˚ KAPITOLA 2. ANALYZA POHYBU CHUZE
Kapitola 3
Hardwarov´ e vybaven´ı laboratoˇ re ch˚ uze V t´eto kapitole jsou uvedeny hardwarov´e technologie, kter´e se obecnˇe v laboratoˇri ch˚ uze vyskytuj´ı. Jak zhruba laboratoˇr vypad´a? Velikost m´ıstnosti je zvolena tak, aby pacient mohl prov´est nˇekolik krok˚ u v pˇr´ım´em smˇeru. Za nutnost se povaˇzuje izolace m´ıstnosti od otˇres˚ u z okoln´ıho prostˇred´ı. I nepatrn´e chvˇen´ı zemˇe m´a vliv na v´ ystupn´ı hodnoty tlakov´e ploˇsiny. Kromˇe tlakov´e ploˇsiny jsou zde jeˇstˇe rozm´ıstˇen´e motion capturov´e kamery, pˇrij´ımaˇc EMG sign´ al˚ u a PC stanice slouˇz´ıc´ı k zaznamen´an´ı a n´asledn´e anal´ yze biometrick´ ych dat. Na n´asleduj´ıc´ım obr´ azku 3.1 je vidˇet sch´ema typick´eho rozvrˇzen´ı vybaven´ı laboratoˇre ch˚ uze. N´asleduje charakteristika jednotliv´ ych technologi´ı.
Obr´ azek 3.1: Laboratoˇr ch˚ uze. Zdroj [3].
19
´ VYBAVEN´I LABORATORE ˇ CHUZE ˚ KAPITOLA 3. HARDWAROVE
20
3.1 3.1.1
Motion Capture Vyuˇ zit´ı
Motion capture, nˇekdy tak´e naz´ yv´ an jako motion tracking nebo zkr´acenˇe MoCap, je technologie slouˇz´ıc´ı k sn´ım´ an´ı pohyb˚ u lid´ı, zv´ıˇrat, popˇr´ıpadˇe jin´ ych objekt˚ u. Sn´ıman´e pohyby se digitalizuj´ı a pˇred´ avaj´ı k dalˇs´ımu zpracov´an´ı. Tato technologie se v praxi vyuˇz´ıv´a pˇredevˇs´ım ve filmov´em a hern´ım pr˚ umyslu. D´ ale sv´e uplatnˇen´ı najde ve sportu, ale tak´e v arm´adˇe. V posledn´ı dobˇe se rozr˚ ust´ a jej´ı uˇzit´ı i v medic´ınˇe, kde l´ekaˇri mohou pomˇernˇe jednoduˇse, rychle a pˇredevˇs´ım detailnˇe zaznamenat pohyby lid´ı a potom je analyzovat. Existuj´ı r˚ uzn´e syst´emy sn´ım´ an´ı pohyb˚ u. Podle toho na jak´em principu jsou zaloˇzeny, je m˚ uˇzeme rozdˇelit do dvou z´ akladn´ıch skupin. Jedn´a se o optick´e syst´emy a neoptick´e syst´emy. Mezi neoptick´e patˇr´ı mechanick´e, magnetick´e a setrvaˇcn´e syst´emy. Jelikoˇz se v laboratoˇri ch˚ uze vyuˇz´ıv´a optick´ ych syst´emu sn´ım´ an´ı pohybu, budu v n´asleduj´ıc´ım textu popisovat pouze je.
3.1.2
Optick´ e syst´ emy
Optick´e syst´emy vynikaj´ı svoj´ı velkou pˇresnost´ı, v´ ykonem a praktiˇcnost´ı. Z´akladem je soustava kamer rozm´ıstˇen´ ych v prostoru a retroreflexn´ı kuliˇcky, tzv. markery, pˇripevnˇen´e ke sledovan´emu objektu - v naˇsem pˇr´ıpadˇe k pacientovi. Kamery jsou opatˇreny skupinou LED, kter´e vyzaˇruj´ı svˇetlo ve smˇeru pohledu kamery. K tomuto u ´ˇcelu se vyuˇz´ıvaj´ı LED emituj´ıc´ı ˇcerven´e svˇetlo, svˇetlo bl´ızk´e infraˇcerven´emu p´asmu a svˇetlo infraˇcerven´e. Od retroreflexn´ıch marker˚ u se generovan´e z´ aˇren´ı odraz´ı zpˇet do kamery. Zde proch´az´ı optick´ ym filtrem, kter´ y propouˇst´ı jenom svˇetlo se stejnou charakteristikou, jakou m´a svˇetlo vyzaˇrovan´e. Moˇznost pˇrizp˚ usobit citlivost kamery, n´ am umoˇzn ˇuje sn´ımat pouze retroreflexn´ı markery a ignorovat odlesky pokoˇzky a ostatn´ı neˇz´ adouc´ı odrazy. Syst´em pouˇz´ıvan´ y v laboratoˇri ch˚ uze dodala firma Vicon [14]. Jedn´a se konkr´etnˇe o variantu Vicon MX. Z´ akladn´ı Vicon MX syst´em sn´ım´an´ı pohyb˚ u (obr. 3.2) se sest´av´a z osmi kamer rozm´ıstˇen´ ych v prostoru. Synchronizace a komunikace mezi jednotliv´ ymi kamerami je zajiˇstˇena prvkem zvan´ ym MX Ultranet. Ten slouˇz´ı tak´e jako zdroj elektrick´e energie pro kamery. D´ıky velk´e sn´ımkov´e rychlosti kamer vznik´a velk´e mnoˇzstv´ı objemn´ ych dat. PC stanici komunikuj´ıc´ı s kamerami je proto tˇreba propojit pˇres MX Ultranet vysokorychlostn´ım gigabitov´ ym ethernetem. V praxi existuj´ı i syst´emy s vˇetˇs´ım mnoˇzstv´ım kamer. U takto rozs´ahl´ ych syst´em˚ u je tˇreba ˇreˇsit jejich spr´avnou synchronizaci a komunikaci pˇrid´an´ım nov´ ych prvk˚ u do s´ıt’ov´e topologie. Pˇred samotn´ ym mˇeˇren´ı je zapotˇreb´ı kamery zkalibrovat. Proces kalibrace zahrnuje identifikaci vnitˇrn´ıch (ohniskov´ a vzd´ alenost a obrazov´e zkreslen´ı) a vnˇejˇs´ıch (pozice a orientace) parametr˚ u kamery. Ke kalibraci se vyuˇz´ıv´a speci´aln´ı h˚ ulka ve tvaru p´ısmene T s pˇeti markery. Ta se um´ıst´ı do prostoru ve kter´em se pot´e bude sn´ımat pohyb a pomoc´ı softwaru v poˇc´ıtaˇci se provede kalibrace. Po kalibraci maj´ı kamery zmapov´anu oblast ve kter´e se nach´azej´ı a znaj´ı svoj´ı vz´ajemnou polohu. Pˇri samotn´em sn´ım´an´ı je pak d˚ uleˇzit´e, aby kaˇzd´ y marker vidˇely alespoˇ n dvˇe kamery. D´ıky tomu lze urˇcit polohu marker˚ u v prostoru. Zn´ame - li polohu vˇsech marker˚ u ve sc´enˇe, uˇz n´ am nic nebr´an´ı v zah´ajen´ı mˇeˇren´ı a dalˇs´ımu zpracov´an´ı dat ve specializovan´ ych softwarech.
´ PLOSINA ˇ 3.2. TLAKOVA AMTI OR6-7
21
Obr´ azek 3.2: Z´ akladn´ı Vicon MX syst´em s osmi kamerami. Zdroj: [11].
3.2 3.2.1
Tlakov´ a ploˇ sina AMTI OR6-7 Vyuˇ zit´ı
Tlakov´ a ploˇsina AMTI OR6-7 je n´astroj umoˇzn ˇuj´ıc´ı detekovat s´ıly p˚ usob´ıc´ı na jej´ı povrch. Vzniku t´eto ploˇsiny pˇredch´ azelo mnoho let v´ yzkumu a v´ yvoje. V´ ysledkem je zaˇr´ızen´ı s vysokou citlivost´ı a pˇresnost´ı. Proto se vyuˇz´ıv´a tak´e v brnˇensk´e laboratoˇri ch˚ uze, kde se zaznamen´ av´ a p˚ usoben´ı pacientova doˇslapu na ploˇsinu. V´ ysledn´e namˇeˇren´e hodnoty maj´ı nemal´ y vliv na v´ yslednou anal´ yzu pacientovi ch˚ uze.
3.2.2
Popis zaˇ r´ızen´ı
Ploˇsina vyˇzaduje instalaci v m´ıstech odst´ınˇen´ ych od okoln´ıch vibrac´ı, kter´e by mohly ovlivnit v´ ysledky mˇeˇren´ı. Dostateˇcn´e upevnˇen´ı je samozˇrejmost´ı. Jak rad´ı v´ yrobce, ploˇsina byla pˇriˇsroubov´ ana k p´ aru kolejnic, kter´e byly dod´any spoleˇcnˇe s ploˇsinou. Tyto kolejnice jsou zapuˇstˇeny do podlahy laboratoˇre a napevno pˇrilepeny dvousloˇzkov´ ym epoxy lepidlem. D´ıky tomu deska ploˇsiny nevystupuje nad u ´roveˇ n zemˇe laboratoˇre. Ploˇsina mˇeˇr´ı s´ıly a momenty p˚ usob´ıc´ı na jej´ı povrch. Zaˇr´ızen´ı umoˇzn ˇuje tyto s´ıly a momenty rozloˇzit do tˇr´ıdimenzion´ aln´ıho souˇradn´eho syst´emu podle os X, Y a Z. Na v´ ystupu je tak ˇsest sloˇzek odpov´ıdaj´ıc´ı tˇrem moment˚ um Mx , My a Mz a sloˇzk´am s´ıly Fx , Fy a Fz . Ploˇsina AMTI OR6-7 vyuˇz´ıv´ a souˇradn´eho syst´emu s pravidlem prav´e ruky. Viz. obr´azek 3.3. Podle pravidla prav´e ruky smˇeˇruje kladn´a osa X ve smˇeru palce, ukazov´ak ukazuje kladn´ y smˇer osy Y a ohnut´ y prostˇredn´ık kladn´ y smˇer osy Z. V tomto pˇr´ıpadˇe je kladn´a osa Z ploˇsiny smˇeˇrov´ ana dol˚ u, kladn´ a osa Y je orientov´ana v opaˇcn´em smˇeru neˇz je um´ıstˇen konektor pro pˇripojen´ı a kladn´ a osa X smˇeˇruje pˇri pohledu kladn´e osy Y doleva. Uvnitˇr ploˇsiny jsou um´ıstˇeny tenzometry. Tenzometr je pˇr´ıstroj umoˇzn ˇuj´ıc´ı mˇeˇren´ı tenze, coˇz je napˇet´ı nebo tlak vyvolan´ y vnˇejˇs´ım p˚ usoben´ım. K mˇeˇren´ı se vyuˇz´ıv´a ˇcidlo, kter´e zaznamen´ av´ a mechanick´e prodlouˇzen´ı tenk´ ych kovov´ ych dr´atk˚ u, ze kter´ ych se tenzometr skl´ ad´ a.
22
´ VYBAVEN´I LABORATORE ˇ CHUZE ˚ KAPITOLA 3. HARDWAROVE
Obr´ azek 3.3: Souˇradn´ y syst´em ploˇsiny AMTI. Zdroj: [12].
Zmˇeny mechanick´eho prodlouˇzen´ı jsou tenzometrem pˇrev´adˇeny na zmˇeny elektrick´eho odporu. Dojde - li k takov´e zmˇenˇe, je vybuzen velmi mal´ y elektrick´ y sign´al, kter´ y je vysl´an do s´ıt’ov´eho mostu. Vˇsechny tenzometry jsou t´ımto mostem propojeny. V´ ystupn´ı analogov´ y sign´al z tenzometrov´eho mostu mus´ı b´ yt zes´ılen, aby byl sign´al dostateˇcnˇe siln´ y pro dalˇs´ı zpracov´an´ı. Takov´ yto sign´ al je digitalizov´ an pomoc´ı A/D pˇrevodn´ıku. Je obecnˇe v´ yhodn´e m´ıt takto z´ıskan´a data ve vhodn´ ych fyzik´ aln´ıch jednotk´ach. Je proto potˇreba aplikovat na z´ıskan´a data spr´avn´e mˇeˇr´ıtko. Po aplikov´ an´ı vhodn´eho mˇeˇr´ıtka je moˇzno data zaˇc´ıt uchov´avat a zobrazovat fin´ aln´ı v´ ysledky.
3.3 3.3.1
Sn´ım´ an´ı EMG Co je EMG?
Pˇri ˇcinnosti sval˚ u se vytv´ aˇr´ı mal´e mikrovoltov´e elektrick´e sign´aly. Tyto elektrick´e sign´aly je moˇzno mˇeˇrit procesem odbornˇe naz´ yvan´ ym jako elektromyografie neboli EMG. Jedn´a se o vyˇsetˇrovac´ı techniku, kter´ a napom´ ah´ a l´ekaˇr˚ um hodnotit funkˇcn´ı stav pohybov´eho syst´emu. Principem EMG je mˇeˇren´ı rychlosti veden´ı vzruchu ve stimulovan´em nervu a velikost elektrick´e odpovˇedi na stimulace ve svalu. Vyˇsetˇren´ı se prov´ad´ı bud’ jehlovou nebo konduktivn´ı technikou. Techniky se liˇs´ı ve zp˚ usobu z´ısk´an´ı informace o svalu. Jehlov´a elektroda se zav´ad´ı vpichem pˇr´ımo do svalu a z´ısk´ av´ a informaci o jednotliv´ ych svalov´ ych vl´aknech. Oproti tomu kunduktivn´ı proces z´ısk´ av´ a informaci o celkov´em stavu svalu. Vyuˇz´ıv´a pˇritom elektrod um´ıstˇen´ ych na povrchu lidsk´e k˚ uˇze. Pro pacienta je v´ yhodou neinvazivnost a bezbolestnost t´eto metody. V´ ysledkem elektromyografie je EMG kˇrivka, kter´a l´ekaˇr˚ um slouˇz´ı jako n´astroj pro detailn´ı anal´ yzu.
´ ´I EMG 3.3. SN´IMAN
3.3.2
23
Vyuˇ zit´ı v laboratoˇ ri ch˚ uze
Pˇri zkoum´ an´ı lidsk´e ch˚ uze bylo zjiˇstˇeno, ˇze EMG poskytuje informace, kter´e je tˇreba pˇri fin´aln´ım rozboru zohlednit. Tento vyˇsetˇrovac´ı postup se proto vyuˇz´ıv´a i v brnˇensk´e laboratoˇri ch˚ uze. Syst´emy, poskytuj´ıc´ı moˇznost tohoto mˇeˇren´ı, dodala firma Noraxon [13], kter´a se zab´ yv´ a EMG a senzorov´ ymi syst´emy. Pro svobodn´ y pohyb pacienta byla zvolena bezdr´atov´ a technologie, kdy pacient m´a na sobˇe pˇripevnˇen vys´ılaˇc TeleMyo 2400T G2 (obr. 3.4) sn´ımaj´ıc´ı EMG konduktivn´ı metodou. Vys´ılan´a data jsou zachycena pˇrij´ımaˇcem TeleMyo 2400R G2 (obr. 3.4) k dalˇs´ımu zpracov´an´ı.
3.3.3
Vys´ılaˇ c TeleMyo 2400T G2
TeleMyo 2400T G2 (obr. 3.4) je bezdr´atov´ y pˇrenositeln´ y vys´ılaˇc nejnovˇejˇs´ı generace, umoˇzn ˇuj´ıc´ı sn´ım´an´ı povrchov´eho EMG a jin´ ych analogov´ ych sign´al˚ u v re´aln´em ˇcase. Vys´ılan´e sign´ aly maj´ı dosah aˇz sto metr˚ u, ˇc´ımˇz zaˇr´ızen´ı zaruˇcuje pˇr´ıjem na dlouh´e vzd´alenosti a pacienta neomezuje v pohybu pˇr´ıpadnou kabel´aˇz´ı. Syst´em zahrnuje ˇctyˇr, osmi, dvan´acti a ˇsestn´ acti kan´alovou vys´ılac´ı jednotku, aktivn´ı pˇredzes´ılen´e elektrody, nab´ıjeˇcku bateri´ı, pouzdro a opasek pro snadnˇejˇs´ı pˇripevnˇen´ı k pacientovi. Standardn´ı syst´em je vybaven 4, 8, 12 nebo 16 aktivn´ımi EMG senzory. Syst´em je navrˇzen tak, ˇze pˇrid´an´ım dalˇs´ıho vys´ılaˇce rozˇs´ıˇr´ı st´avaj´ıc´ı syst´em aˇz na 32 kan´ al˚ u kombinuj´ıc´ıch EMG a jin´e analogov´e sign´aly. Prakticky lze pˇripojit jak´ ykoliv analogov´ y senzor, kter´ y je nap´ajen +/- 5 V.
Obr´azek 3.4: Vlevo: vys´ılaˇc TeleMyo 2400T G2, vpravo: pˇrij´ımaˇc TeleMyo 2400R G2. Zdroj:[13]
3.3.4
Pˇ rij´ımaˇ c TeleMyo 2400R G2
Pˇrenositeln´ y vys´ılaˇc TeleMyo 2400T G2 je tˇreba doplnit odpov´ıdaj´ıc´ım bezdr´atov´ ym pˇrij´ımaˇcem sign´al˚ u. Pro potˇreby laboratoˇre ch˚ uze byl z nˇekolika moˇznost´ı vybr´an pˇrij´ımaˇc TeleMyo 2400R G2 (obr. 3.4). Stanice poskytuje konverzi digit´aln´ıch telemetrick´ ych dat na analogov´ y v´ ystupn´ı sign´ al. Kromˇe analogov´eho v´ ystupu mohou b´ yt data z´ısk´ana v digit´aln´ı formˇe pˇripojen´ım k PC pˇres USB rozhran´ı. Zaˇr´ızen´ı nab´ız´ı 32 kan´al˚ u pro r˚ uzn´e mˇeˇr´ıc´ı pˇr´ıstroje
24
´ VYBAVEN´I LABORATORE ˇ CHUZE ˚ KAPITOLA 3. HARDWAROVE
vys´ılaj´ıc´ı telemetrick´ y sign´ al. Tyto kan´ aly slouˇz´ı pro pˇr´ıjem EMG a jin´ ych analogov´ ych dat, z´ıskan´ ych napˇr. z tlakov´e ploˇsiny nebo pˇr´ıstroj˚ u pro anal´ yzu pohybu. Uvnitˇr stanice je synchronizaˇcn´ı trigger poskytuj´ıc´ı ˇcasovou synchronizaci zaˇc´atku mˇeˇren´ı dat vys´ılaˇce TeleMyo 2400T G2 s ostatn´ımi zaˇr´ızen´ımi v laboratoˇri ch˚ uze. V´ yhodou tohoto pˇrij´ımaˇce je moˇznost pr´ace s vys´ılaˇcem TeleMyo 2400T G2 i bez nutnosti pˇripojen´ı k poˇc´ıtaˇci.
3.4
Digit´ aln´ı video kamery
V laboratoˇri ch˚ uze se tak´e vyskytuj´ı dvˇe digit´aln´ı video kamery. Jsou zde um´ıstˇeny kv˚ uli re´aln´emu z´aznamu pacientovi ch˚ uze. Tento referenˇcn´ı z´aznam je pˇri z´avˇereˇcn´e anal´ yze takt´eˇz zohlednˇen. Kaˇzd´ a kamera sn´ım´ a pr˚ ubˇeh jit´ı z jin´e strany. Jedna zaznamen´av´a ch˚ uzi z boˇcn´ıho pohledu, druh´ a z ˇceln´ıho.
Kapitola 4
Softwarov´ e vybaven´ı laboratoˇ re ch˚ uze Cel´ y proces zpracov´ an´ı a vyhodnocov´an´ı biometrick´ ych dat obstar´avaj´ı dvˇe hlavn´ı aplikace. Jedn´a se o Vicon Nexus a Vicon Polygon. Tyto aplikace bˇeˇz´ı na operaˇcn´ım syst´emu Microsoft Windows XP. Vicon Nexus by se dal pˇrirovnat k rozhran´ı, umoˇzn ˇuj´ıc´ı obsluhuj´ıc´ımu person´ alu ovl´ adat hardwarov´e vybaven´ı laboratoˇre. Pˇritom poskytuje prostˇredky k z´aznamu informace z´ıskan´e z motion capturov´ ych kamer, tlakov´e ploˇsiny, EMG a pomocn´ ych DV kamer. Oproti tomu Vicon Polygon slouˇz´ı jako interaktivn´ı vizualizaˇcn´ı a reportovac´ı n´astroj, vhodn´ y pro anal´ yzu a vizualizaci zaznamenan´ ych 3D dat. Polygon um´ı ˇc´ıst data vytvoˇren´ a Vicon Nexusem, kter´ a obsahuj´ı modelov´a data vytvoˇren´a biomechanick´ ym modelaˇcn´ım softwarem (Plug-in Gait 4.1). V´ ysledn´a vizualizace a z´avˇereˇcn´a zpr´ava spoleˇcnˇe dovoluj´ı prezentovat v´ ysledky nasn´ıman´ ych dat, ze kter´ ych se pot´e stanov´ı pacientova diagn´oza a n´asleduj´ıc´ı postup l´eˇcby. Polygon umoˇzn ˇuje sjednotit vˇsechny z´ıskan´e informace (grafy, ilustraˇcn´ı videa, 3D renderovanou grafiku, pozn´ amky, atd.) do jednoho reportu.
4.1
Vicon Nexus
Vicon Nexus je specifick´ y software pouˇz´ıvan´ y pro v´ yzkumnou ˇcinnost v oblasti sn´ım´ an´ı pohybu. Vicon Nexus byl od zaˇc´ atku vyv´ıjen pro vˇedeck´e u ´ˇcely. D´ıky tomu umoˇzn ˇuje snadnˇe a flexibilnˇe prov´ adˇet kaˇzdodenn´ı rutinn´ı u ´koly. Nexus je jednoduch´ y na nauˇcen´ı a pouˇz´ıv´an´ı. Byl vyv´ıjen podle modern´ıch princip˚ u uˇzivatelsk´eho rozhran´ı. Pouˇz´ıvan´e n´astroje a nastaven´ı jsou jasnˇe a logicky rozm´ıstˇeny a d´ıky tomu mohou b´ yt nalezeny na m´ıstech, kde se jejich v´ yskyt oˇcek´av´a. Nic nen´ı skryto v nejasn´ ych menu nebo dialoz´ıch. Jedn´ a se o real time aplikaci, kter´a dok´aˇze zobrazovat v re´aln´em ˇcase informaci o markerech, analogov´ a data, referenˇcn´ı videoz´aznam, stejnˇe jako souˇcasn´ y stav zaˇr´ızen´ı. Velice zaj´ımavou schopnost´ı Vicon Nexusu je moˇznost pˇrekryt´ı videoz´aznamu z DV kamery 3D informac´ı. Do videoz´aznamu tak m˚ uˇze b´ yt prom´ıtnut silov´ y vektor doˇslapu sn´ıman´e osoby nebo pozice sn´ıman´ ych marker˚ u.
25
26
´ VYBAVEN´I LABORATORE ˇ CHUZE ˚ KAPITOLA 4. SOFTWAROVE
Obr´ azek 4.1: Okno aplikace Vicon Nexus.
Pˇri prvn´ım spuˇstˇen´ı Vicon Nexusu je zapotˇreb´ı prov´est konfiguraci vˇsech syst´em˚ u. Jelikoˇz tento dokument nen´ı prim´ arnˇe zamˇeˇren na optick´e syst´emy a jejich nastaven´ı, budou nast´ınˇeny pouze podstatn´e parametry: • Kvalitn´ı nasn´ım´ an´ı pozic retroreflexn´ıch marker˚ u ovlivˇ nuje nastaven´ı Grayscale Mode, kter´e ˇr´ık´ a jak´ a data budou posl´ ana z kamery do poˇc´ıtaˇce. • Maximum Blob Height ud´ av´ a jak´ y je maxim´aln´ı poˇcet vertik´aln´ıch pixel˚ u markeru, kter´ y senzor kamery jeˇstˇe pˇrijme. • Camera Strobe Intensity mˇen´ı mnoˇzstv´ı emitovan´eho svˇetla kamery. • Camera Threshold urˇcuje jak jasn´ y odraz od markeru je pˇrijateln´ y kamerou. D´ale je zapotˇreb´ı v Nexusu pˇridat a nakonfigurovat tlakovou ploˇsinu a EMG zaˇr´ızen´ı. U tlakov´e ploˇsiny je tˇreba nastavit vlastn´ı identifikaci - unik´atn´ı jm´eno nebo jin´ y identifik´ator, zdroj - specifikuj´ıc´ı fyzick´e pˇripojen´ı zaˇr´ızen´ı, pozici geometrick´eho stˇredu ve vztahu k poˇc´atku sn´ıman´e oblasti a pˇriloˇzit kalibraˇcn´ı soubor, dod´avan´ y s ploˇsinou. Po pˇrid´an´ı nov´eho EMG zaˇr´ızen´ı do Vicon Nexusu se nastav´ı jeho identifikace a zdroj stejnˇe jako u tlakov´e ploˇsiny, nav´ıc se jeˇstˇe urˇc´ı poˇcet jeho kan´al˚ u. Jeˇstˇe zb´ yv´ a pˇrid´ an´ı a nastaven´ı DV kamer, slouˇz´ıc´ıch k re´aln´emu n´ahledu na ch˚ uzi pacienta. Pokud to kamera umoˇzn ˇuje, lze nastavit sn´ımkovou frekvenci, ohniskovou vzd´alenost optiky DV kamer a jin´ a nastaven´ı. Po proveden´ı vˇsech potˇrebn´ ych nastaven´ı, je nutnost´ı kalibrace syst´emu. V´ yrobce doporuˇcuje kalibraci na zaˇc´ atku kaˇzd´eho dne, kdy se prov´ad´ı mˇeˇren´ı a to jak motion capturov´ ych kamer, tak i DV kamer. Kalibrace motion capturov´ ych kamer byla nast´ınˇena v sekci
4.2. VICON POLYGON
27
3.1.2. Kalibrace DV kamer se prov´ad´ı podobn´ ym zp˚ usobem. Takt´eˇz se pouˇz´ıv´a h˚ ulka ve tvaru p´ısmene T s pˇeti markery. Je-li syst´em nakonfigurov´an, nic nebr´an´ı v z´aznamu dat.
Plug-in Gait Souˇc´ast´ı Vicon Nexusu je Plug-in Gait. Plug-in Gait je biomechanick´ y model implementuj´ıc´ı konvenˇcn´ı model ch˚ uze (viz. text z kapitoly 2). Ten umoˇzn ˇuje v´ ypoˇcet kloubn´ı kinematiky a kinetiky na z´ akladˇe pozice marker˚ u a antropometrick´ ych mˇeˇren´ı pacienta tak, jak je pops´ ano v kapitole 2.
4.2
Vicon Polygon
Vicon Polygon je vizualizaˇcn´ı a reportovac´ı n´astroj slouˇz´ıc´ı k anal´ yze nasn´ıman´ ych 3D dat. Poskytuje jednotn´e prostˇred´ı, ve kter´em uˇzivatel´e mohou zobrazovat pohyb, pouˇz´ıvat grafy, ps´at strukturovan´e texty, pˇrehr´ avat videa a vyuˇz´ıvat dalˇs´ı poskytovan´e funkce k u ´pln´e anal´ yze hybnosti pacienta. Pˇritom Polygon dok´aˇze v´ ysledky reprodukovat v elektronick´e nebo tiˇstˇen´e podobˇe.
4.2.1
Muskuloskelet´ aln´ı syst´ em
Kl´ıˇcov´ ym prvkem Vicon Polygonu je Musculoskeletal Modeling Module. Tento modul poskytuje funkce pro anal´ yzu muskuloskelet´aln´ıho syst´emu, nˇekdy naz´ yvan´eho jako lokomoˇcn´ı syst´em. V z´ asadˇe tento syst´em m˚ uˇzeme rozdˇelit na syst´em sval˚ u, kter´ y je vlastn´ım vykonavatelem pohybu a syst´em kost´ı, kam kromˇe kost´ı m˚ uˇzeme zaˇradit klouby, vazy a ˇslachy. Tento komplexn´ı syst´em je velice d˚ uleˇzit´ y pro spr´avn´e drˇzen´ı tˇela a samotn´ y pohyb ˇclovˇeka. Polygon je proto velmi cenn´ ym n´astrojem pro diagnostikov´an´ı poruch muskuloskelet´aln´ıho syst´emu pacient˚ u. Skelet´ aln´ı model Velmi d˚ uleˇzitou souˇc´ ast´ı kaˇzd´eho muskuloskelet´aln´ıho modulu je vlastn´ı geometrie a hierarchie kostry. Muskuloskelet´ aln´ı modely Polygonu (a soubory kter´e je definuj´ı) obsahuj´ı popis kostry jako sekvenci kost´ı spojenou klouby. Svaly Svaly jsou nejd˚ uleˇzitˇejˇs´ı ˇc´ ast´ı muskuloskelet´aln´ıho modelu. Jejich pozice, d´elky a ostatn´ı charakteristiky jsou urˇcen´e 3D daty nasn´ıman´ ymi Vicon Nexusem, jeˇz d´avaj´ı nahl´ednout do jejich biomechanick´eho chov´ an´ı. Modelovan´e svaly jsou pops´ any jejich napojen´ım ke kost´ım kostry v normalizovan´em souˇradn´em syst´emu segmentu. Jedn´a se o standardn´ı postup kdy pozice napojen´ı jsou dˇeleny d´elkou kosti tak, aby tyto pozice mohly b´ yt mapov´any na kosti odliˇsn´ ych d´elek a tud´ıˇz i na r˚ uzn´e pacienty.
´ VYBAVEN´I LABORATORE ˇ CHUZE ˚ KAPITOLA 4. SOFTWAROVE
28
Obr´azek 4.2: Okno aplikace Vicon Polygon s popsan´ ym uˇzivatelsk´ ym rozhran´ım. Zdroj: [10]
Sval je pops´ an sv´ ym poˇc´ atkem (muscle origin) a koncem (muscle insertion). Muscle origin je napojen´ı na kost, kter´e se pˇri kontrakci svalu neh´ ybe, zat´ımco muscle insertion odpov´ıd´a konci svalu pˇripojen´eho k pohybuj´ıc´ı se kosti, pˇri jehoˇz kontrakci se mˇen´ı jeho poloha. Svaly mohou b´ yt pops´ any svalov´ ym modelem, kter´ y stanovuje jak se mohou mˇenit s´ıly kter´e sval produkuje, v z´ avislosti na jeho d´elce a rychlosti zkr´acen´ı. K tˇemto kalkulac´ım se nejˇcastˇeji vyuˇz´ıv´ a Hill˚ uv model svalu. Kromˇe vlastn´ıho modelu, m´a kaˇzd´ y sval tak´e unik´atn´ı parametry, kter´e pˇrizp˚ usobuj´ı mˇeˇr´ıtko modelu podle velikosti, d´elky a typu svalu. V´ıce informac´ı o t´eto problematice se m˚ uˇzete doˇc´ıst v textu od F. E. Zajaka [16].
4.2.2
Tˇ eˇ ziˇ stˇ e
Vicon Polygon disponuje moˇznost´ı zobrazen´ı tˇeˇziˇstˇe (center of mass) a vykreslen´ı souˇradnic jednotliv´ ych kost´ı segmentu nebo cel´eho tˇela. To dovoluje analyzovat rovnov´ahu, posturiku a energii lidsk´eho pohybu.
Kapitola 5
Postup sn´ıman´ı pohybu v laboratoˇ ri ch˚ uze Nejprve je tˇreba objasnit a ˇcten´ aˇre sezn´amit s typy soubor˚ u, kter´e se poˇz´ıvaj´ı k vytv´ aˇren´ı model˚ u a zaznamen´ an´ı informac´ı o pacientovi. Jedn´a se pˇredevˇs´ım o soubory s koncovkami .vst, .vsk, .mkr a .mp. N´ asleduje jejich popis: Vicon Skeleton Template (.vst) - ˇsablona kter´a specifikuje obecn´e vztahy mezi segmenty tˇela a klouby. Zahrnuje: • jm´ena jednotliv´ ych marker˚ u (napˇr. ANK, KNE), • kostern´ı strukturu sn´ıman´eho subjektu (napˇr. spojen´ı p´anve a stehna, stehna a holenˇe), • vztah mezi markery a strukturou kostry (napˇr. LASI marker patˇr´ı do segmentu p´ anve), • typ kloub˚ u spojuj´ıc´ı jednotliv´e segmenty (napˇr. kulovit´e, v´alcovit´e,...), • dalˇs´ı vlastnosti segment˚ u, marker˚ u a kloub˚ u (napˇr. barva kterou jsou v 3D pracovn´ı ploˇse segmenty zv´ yraznˇeny, atd.) Vicon Skeleton (.vsk) - vznik´ a aplikov´an´ım mˇeˇr´ıtka konkr´etn´ıho pacienta na Vicon Skeleton Template. Tento soubor se vytv´aˇr´ı individu´alnˇe pro kaˇzd´eho pacienta. Vicon Marker File (.mkr) - soubor, kter´ y specifikuje oznaˇcen´ı pouˇz´ıvan´e pro markery. Nexus tento soubor vytv´ aˇr´ı automaticky, kdyˇz se vytv´aˇr´ı Vicon Skeleton Template. Vicon Model Parameters File (.mp) - soubor, kter´ y uchov´av´a statick´e informace o pacientovi (v´ yˇska, hmotnost,...).
29
30
5.1
ˇ CHUZE ˚ KAPITOLA 5. POSTUP SN´IMAN´I POHYBU V LABORATORI
Nasn´ım´ an´ı dat
Jak bylo zm´ınˇeno v pˇredchoz´ı kapitole, k sn´ım´an´ı pacientovi ch˚ uze se vyuˇz´ıv´a aplikace zvan´a Vicon Nexus. Pokud byla provedena u ´spˇeˇsn´a kalibrace syst´emu, m˚ uˇze b´ yt provedeno mˇeˇren´ı. Pracovn´ı postup sn´ım´ an´ı ch˚ uze pomoc´ı Vicon Nexus s Plug-in Gait se skl´ad´a z n´asleduj´ıc´ıch ˇc´ast´ı:
5.1.1
Pˇ r´ıprava subjektu
V prv´e ˇradˇe je zapotˇreb´ı vybrat spr´avn´ y Vicon Skeleton Template. Vicon Nexus v sobˇe implementuje nˇekolik ˇsablon pro sn´ım´ an´ı lidsk´e ch˚ uze. Pro potˇreby laboratoˇre ch˚ uze z Fakultn´ı dˇetsk´e nemocnice J.G. Mendela v Brnˇe byla za vyhovuj´ıc´ı shled´ ana ˇsablona OxfordFootModel BILATERAL.vst. Ta obsahuje popis doln´ı ˇc´asti lidsk´eho tˇela zaloˇzen´eho na OFM popsan´eho v dˇr´ıve zm´ınˇen´em textu 2.1.4. V syst´emu je kaˇzd´ y pacient reprezentov´an stromem, jehoˇz koˇren obsahuje jm´eno pacienta. Strom je d´ale tvoˇren podstromy obsahuj´ıc´ımi informaci o markerech, segmentech tˇela, kloubech a v´ ystupn´ıch hodnot´ach, kter´e generuje Plug-in Gait 4.1. Po vytvoˇren´ı nov´eho stromu, potˇrebuje Nexus zn´at pacientovi antropometrick´e parametry. Ty se jednoduˇse vypln´ı do vstupn´ıch pol´ı v z´aloˇzce pacienta. V t´eto f´azi m˚ uˇze obsluhuj´ıc´ı person´ al pˇripevnit retroreflexn´ı markery na pacienta zp˚ usobem popsan´ ym v kapitole 2. Jsou-li pˇripevnˇeny markery k pacientovi, vytvoˇr´ı se statick´ y z´aznam a zrekonstruuj´ı se markery v syst´emu. Tento krok je nutn´ y k vytvoˇren´ı 3D reprezentace Obr´azek 5.1: Zv´ yraznˇen´e pacienta ve Vicon Nexusu. Pˇri vytv´aˇren´ı statick´eho z´aznamu povinn´e hodnoty antropo- by mˇel pacient st´at v sn´ıman´e oblasti v neutr´aln´ı statick´e metrick´ ych parametr˚ u. p´ oze, jako je p´oza ve tvaru p´ısmene T. Po vytvoˇren´ı statick´eho z´ aznamu se mus´ı oznaˇcit rekonstruovan´e 3D markery, aby se asociovaly s tˇemi definovan´ ymi v ˇsablonˇe OxfordFootModel BILATERAL.vst. Na v´ ybˇer je bud’ automatick´e nebo manu´ aln´ı oznaˇcen´ı. Pot´e se jeˇstˇe spust´ı kalibraˇcn´ı proces, kter´ y pˇrizp˚ usob´ı Vicon Skeleton Template pro pr´avˇe sn´ıman´eho pacienta. Kdyˇz uˇz je pacient zkalibrov´an, Nexus vytvoˇr´ı pacient˚ uv specifick´ y Vicon Skeleton (.vsk) popisuj´ıc´ı sadu jeho marker˚ u. Po uloˇzen´ı pacientova Vicon Skeletonu m˚ uˇzeme uskuteˇcnit tolik dynamick´ ych z´aznam˚ u kolik potˇrebujeme.
´ ´I DAT 5.1. NASN´IMAN
5.1.2
31
Sn´ım´ an´ı a zpracov´ an´ı dynamick´ eho z´ aznamu
Pˇri sn´ım´ an´ı a zpracov´ an´ı dynamick´eho z´aznamu pacienta se prov´ad´ı n´ asleduj´ıc´ı: Sn´ım´ an´ı dynamick´ eho z´ aznamu - jedn´a se o sn´ım´an´ı dat bˇehem kter´eho se pacient, na nˇemˇz jsou rozm´ıstˇeny markery podle specifikace ve Vicon Skeleton (.vsk), proch´ az´ı ve sn´ıman´em prostoru. Rekonstrukce a automatick´ e oznaˇ cen´ı marker˚ u - Vicon Nexus pouˇz´ıv´ a .vsk soubor k automatick´emu oznaˇcen´ı marker˚ u dynamick´eho z´ aznamu. Pokud automatick´e oznaˇcen´ı marker˚ u selˇze, mus´ı se pouˇz´ıt oznaˇcen´ı manu´aln´ı. D´ıky rekonstrukci marker˚ u mohou b´ yt spoˇc´ıt´any lok´aln´ı souˇradn´e syst´emy jednotliv´ ych segment˚ u (viz. obr. 5.2). Oznaˇcen´e markery potom slouˇz´ı ke kalkulaci kloubn´ı kinematiky a kinetiky, kterou m´a na starosti Plug-in Gait. D´ ale je zapotˇreb´ı na ˇcasov´e ose oznaˇcit sn´ımek, kdy nastal pacient˚ uv doˇslap a opuˇstˇen´ı ˇspiˇcky chodidla od zemˇe. Vhodn´e je tak´e oˇr´ıznut´ı ˇcasov´e osy z´aznamu zleva, popˇr´ıpadˇe zprava tak, aby na z´aznamu nebyly okamˇziky, Obr´azek 5.2: Zrekonstruovan´ y kter´e nejsou pro anal´ yzu potˇrebn´e. z´aznam.
5.1.3
Pˇ rezkoum´ an´ı a v´ yplˇ n mezer
Bˇehem vytv´ aˇren´ı dynamick´eho z´aznamu se m˚ uˇze st´at, ˇze nˇekter´e markery jsou na kr´ atk´ y okamˇzik zast´ınˇeny a ˇz´ adn´ a z kamer jejich pozici nezachyt´ı. Tento okamˇzik je ve v´ ysledn´e trajektorii oznaˇcov´ an jako trhlina nebo mezera, kter´a m´a trv´an´ı nˇekolika sn´ımk˚ u. Vicon Nexus poskytuje prostˇredky jak tyto mezery dopoˇc´ıtat. Opˇet m´ame na v´ ybˇer mezi automatick´ ym nebo manu´ aln´ım postupem. Pˇri automatick´em dopoˇc´ıt´an´ı se stanov´ı maxim´ aln´ı velikost mezer, kter´e se maj´ı dopoˇc´ıtat. Tato operace je zaloˇzena na kubick´e interpolaci kˇrivky. Je proto vhodn´e omezit maxim´aln´ı velikost mezer na pomˇernˇe n´ızk´ y poˇcet sn´ımk˚ u. Typicky 10% sn´ımkov´e frekvence. Pˇri vyˇsˇs´ım poˇctu sn´ımk˚ u doch´az´ı k nepˇresnosti v´ ypoˇctu. Neodstran´ı-li automatick´ y postup vˇsechny mezery, mus´ı se jednotliv´e mezery vyplnit bud’ pomoc´ı nˇejak´eho vzoru, tzn. stejnou ˇc´ast´ı podobn´e trajektorie bez mezer (ˇcasto se vyuˇz´ıv´ a u marker˚ u p´ anve, kter´e maj´ı shodnou trajektorii) nebo pro velk´e mezery pouˇz´ıt spline kˇrivku.
5.1.4
Export dat
Po u ´pravˇe nasn´ıman´ ych dat m˚ uˇzeme prov´est jejich export. Vyexportovan´a biometrick´a data jsou n´ aslednˇe k dispozici pro Vicon Polygon, kde se d´ale analyzuj´ı. Vicon Nexus poskytuje nˇekolik operac´ı pro import a export dat z/do souboru. Pro naˇse u ´ˇcely jsou podstatn´e n´asleduj´ıc´ı operace: Export 3D Overlay Video - tato operace exportuje referenˇcn´ı video z DV kamer ve form´ atu AVI, pˇriˇcemˇz jeˇstˇe m˚ uˇze video obraz doplnit o n´ahled 3D pracovn´ı plochy. M˚ uˇze tak b´ yt do video nahr´ avky prom´ıtnut napˇr. vektor s´ıly doˇslapu pacienta.
32
ˇ CHUZE ˚ KAPITOLA 5. POSTUP SN´IMAN´I POHYBU V LABORATORI
Export C3D - jedn´ a se o export dat nasn´ıman´ ych motion capturov´ ymi kamerami (tzn. markery a jejich trajektorie). Zde jsou tak´e v´ ystupn´ı parametry (kinetika a kinematika) spoˇcten´e Plug-in Gaitem. Vˇse je uloˇzeno do C3D form´atu. Export MP - exportuje parametry vyˇsetˇrovan´eho pacienta (v´ yˇska, v´aha, d´elka nohy,...) do .MP souboru. Export XCP - umoˇzn ˇuje export kalibraˇcn´ıch informac´ı proveden´eho z´aznamu do kalibraˇcn´ıho souboru .XCP.
5.2
Zpracov´ an´ı vyexportovan´ ych dat
Soubory vyexportovan´e z Vicon Nexusu, se vyuˇz´ıvaj´ı v aplikaci Vicon Polygon pro biomechanick´e vizualizace. Pˇredevˇs´ım C3D a VST/VSK soubory, d´ale MKR soubor s popisem vyuˇz´ıvan´ ych marker˚ u. Jak bylo zm´ınˇeno v pˇredchoz´ı kapitole, Vicon Polygon poskytuje funkce pro vypracov´an´ı z´avˇereˇcn´e zpr´ avy. Ta m˚ uˇze obsahovat texty, videoz´aznamy z DV kamer, grafy, 3D data, jeˇz lze zobrazit v 3D pracovn´ı ploˇse a dalˇs´ı poloˇzky. Vˇetˇsina obsahu je tvoˇrena pr´avˇe daty z Vicon Nexusu. ˇ Velice v´ yhodn´ a se jev´ı moˇznost vytvoˇren´ı vlastn´ı ˇsablony reportu. Sablona definuje hieˇ rarchickou strukturu a obsah reportu. Sablona obsahuje typy soubor˚ u a odkazy, kter´e maj´ı b´ yt zahrnuty do vˇsech novˇe vznikaj´ıc´ıch report˚ u.
5.2.1
Popis uˇ zivatelsk´ eho rozhran´ı
Pohodln´e a intuitivn´ı vytv´ aˇren´ı z´ avˇereˇcn´e zpr´avy z´avis´ı na kvalitˇe uˇzivatelsk´eho rozhran´ı. Na obr´azku 4.2 si m˚ uˇzete prohl´ednout uˇzivatelsk´e rozhran´ı aplikace Vicon Polygon. Jedn´a se o klasickou okenn´ı Windows aplikaci. V horn´ı ˇc´asti okna je menu liˇsta s n´ astrojovou liˇstou. V doln´ı ˇc´ asti aplikace je ˇcasov´ a osa a stavov´ a liˇsta. V lev´e ˇc´asti okna je navigace, kter´a poskytuje snadn´e listov´ an´ı strukturou cel´eho reportu. D´a se ˇr´ıci, ˇze toto rozm´ıstˇen´ı je standardn´ı pro takov´ yto typ aplikac´ı. Aplikace Vicon Polygon se stala podnˇetem pro tuto pr´aci. Z´akladn´ı myˇslenka tohoto programu je zcela spr´ avn´ a. Poskytnout l´ekaˇr˚ um n´astroj, kter´ y jim umoˇzn´ı prezentovat a analyzovat biometrick´ a data a vytv´ aˇret a archivovat posudkov´e zpr´avy. Probl´em nen´ı ani tak ve funkˇcnosti, jako v samotn´em ovl´ ad´an´ı a moˇznostech vizualizace. Na trhu existuje velk´e mnoˇzstv´ı komerˇcn´ıch i nekomerˇcn´ıch aplikac´ı. Nˇekter´e z nich disponuj´ı vˇetˇs´ı, nˇekter´e menˇs´ı m´ırou pouˇzitelnosti (usability). Vicon Polygon svoj´ı pouˇzitelnost´ı patˇr´ı do podpr˚ umˇeru dneˇsn´ıch program˚ u. Ovl´ad´an´ı je pomˇernˇe krkolomn´e a obˇcas i frustruj´ıc´ı. Aby si ˇcten´ aˇr udˇelal pˇredstavu, budou uvedeny na n´asleduj´ıc´ıch obr´azc´ıch spoleˇcnˇe s popisem nˇekter´e bˇeˇznˇe pouˇz´ıvan´e operace. Nejprve obr´azek 5.3 toho, jak nepˇrehlednˇe m˚ uˇze tak´e Vicon Polygon vypadat:
´ ´I VYEXPORTOVANYCH ´ 5.2. ZPRACOVAN DAT
33
Obr´ azek 5.3: Okno aplikace Vicon Polygon s v´ıce otevˇren´ ymi dialogy.
5.2.1.1
Zobrazen´ı a skryt´ı mnoˇ ziny kost´ı
Pro zobrazen´ı kost´ı v 3D pracovn´ı ploˇse je tˇreba v Data Object View (obr. 5.4) dvakr´ at kliknout lev´ ym tlaˇc´ıtkem myˇsi na poloˇzku Bones, pˇr´ıpadnˇe prav´ ym tlaˇc´ıtkem vyvolat kontextovou nab´ıdku a vybrat Show Objects. Co m´a ale uˇzivatel dˇelat, potˇrebuje - li mnoˇzinu v 3D ploˇse skr´ yt? Vˇetˇsina uˇzivatel˚ u by oˇcek´ avala, ˇze opˇetovn´ ym dvojklikem na poloˇzku Bones nebo moˇznost´ı v kontextov´em menu, bude dosaˇzeno k´ yˇzen´eho efektu. Ovˇsem marnˇe. V kontextov´em menu je pouze moˇznost jejich zobrazen´ı. Skryt´ı t´eto mnoˇziny je moˇzn´e pouze po kliknut´ı lev´ ym a pak prav´ ym tlaˇc´ıtkem myˇsi na jednu z kost´ı v 3D pracovn´ı ploˇse a zvolen´ı moˇznosti Display Set > Remove z kontextov´e nab´ıdky (obr. 5.5). Zde Obr´azek 5.4: Data Object View. je tˇreba zd˚ uraznit, ˇze nestaˇc´ı nam´ıˇrit kurzor myˇsi na jednu z kost´ı a pak kliknout prav´ ym tlaˇc´ıtkem myˇsi. Mus´ı se nejprve kost vybrat lev´ ym tlaˇc´ıtkem, jinak moˇznost Display Set nen´ı zobrazena. Na stejn´em principu funguje i zobrazov´ an´ı/skr´ yv´ an´ı ostatn´ıch prvk˚ u, jako jsou markery, bounding boxy a dalˇs´ı.
34
ˇ CHUZE ˚ KAPITOLA 5. POSTUP SN´IMAN´I POHYBU V LABORATORI
Obr´azek 5.5: Kontextov´ a nab´ıdka v 3D pracovn´ı ploˇse. Zde zobrazena moˇznost pro skryt´ı mnoˇziny kost´ı.
5.2.1.2
Panel zobrazen´ı (View pane)
Na obr´ azku 5.3 si m˚ uˇzete vˇsimnout, jak m˚ uˇze vypadat zobrazen´ı v´ıce komponent v panelu zobrazen´ı. Zde m˚ uˇze b´ yt v jednom okamˇziku zobrazena 3D pracovn´ı plocha, ˇ adn´a z komponent n´ ahled video nahr´avky, grafy atd. Z´ v sobˇe neobsahuje ikonku pro pˇr´ıpadnou minimalizaci nebo u ´pln´e uzavˇren´ı komponenty. Uzav´ır´an´ı tˇechto komponent se ad´ı kliknut´ım na pˇr´ısluˇsnou ikonu v panelu n´astroj˚ u. Obr´azek 5.6: Dialog pro smaz´ an´ı prov´ T´ ım se zobraz´ ı dialogov´ e okno s v´ y zvou vybr´ a n´ ı komponenty, komponenty. kter´ a m´ a b´ yt odstranˇena (obr. 5.6). Aˇz pot´e se zvol´ı kompo’ nenta k uzavˇren´ı. Nejedn´ a se o u ´plnˇe nejˇst astnˇejˇs´ı ˇreˇsen´ı. Pro uzavˇren´ı takov´eto komponenty by mˇelo staˇcit pouze jedno kliknut´ı.
5.2.1.3
Zobrazen´ı grafu (Graph View)
Pr´ace s grafem je takt´eˇz velice neintuitivn´ı. Chce - li uˇzivatel napˇr´ıklad zoomovat zobrazen´ y graf, mus´ı kliknout na jednu z os grafu prav´ ym tlaˇc´ıtkem myˇsi a zvolit jednu z moˇznost´ı z kontextov´e nab´ıdky Zoom Out / Zoom In. Zde je tˇreba podotknout, ˇze klikne - li se prav´ ym tlaˇc´ıtkem myˇsi do plochy grafu, zobraz´ı se u ´plnˇe jin´a kontextov´a nab´ıdka neˇz ta na os´ach grafu. Co ale uˇzivateli m˚ uˇze nejv´ıce vadit, je pˇrech´azen´ı mezi v procentech normalizovan´em grafu a grafem odpov´ıdaj´ıc´ım celkov´emu trv´ an´ı z´aznamu ch˚ uze. Pˇri prvn´ım zobrazen´ı je vykreslen normalizovan´ y graf. Klikne-li se dvakr´ at na kˇrivku grafu, vykresl´ı se graf urˇcen´ y celkov´ ym ˇcasem z´aznamu. Jak´e je pˇrekvapen´ı, kdyˇz opˇetovn´ y dvojklik nevr´at´ı zpˇet normalizovan´ y graf. Opˇet je to ˇreˇseno na jin´em m´ıstˇe a to ve zmiˇ novan´em kontextov´em menu osy grafu.
5.3. ZHODNOCEN´I APLIKAC´I
5.3
35
Zhodnocen´ı aplikac´ı
Vicon Nexus je ˇspiˇckov´ y komplexn´ı program pro z´ısk´an´ı dat z r˚ uzn´ ych typ˚ u zaˇr´ızen´ı. Jeho podpora pro export dat ve vhodn´ ych form´atech je excelentn´ı. Sv´ ymi vlastnostmi se ide´ alnˇe hod´ı pro potˇreby laboratoˇre ch˚ uze. Vicon Polygon m´ a sice skvˇel´e vyhodnocovac´ı prostˇredky, ale kvalitu shazuje jeho nepˇr´ıvˇetiv´e uˇzivatelsk´e rozhran´ı. Pˇritom tento n´astroj je pro vyhodnocen´ı pacientovi ch˚ uze kl´ıˇcov´ y. Probl´em je ten, ˇze uˇz se d´ ale nevyv´ıj´ı. V d˚ usledku toho vznik´a i tato pr´ace, kter´a by mˇela poloˇzit z´ aklad pro nov´ y a lepˇs´ı syst´em vizualizace biometrick´ ych dat.
36
ˇ CHUZE ˚ KAPITOLA 5. POSTUP SN´IMAN´I POHYBU V LABORATORI
Kapitola 6
Zobrazovaˇ c namˇ eˇ ren´ ych dat laboratoˇ re ch˚ uze Jak jiˇz bylo pops´ ano v pˇredchoz´ı kapitole 5, aplikace Vicon Polygon m´a znaˇcn´e nedostatky v oblasti uˇzivatelsk´eho rozhran´ı. Z´aroveˇ n zastaven´ı v´ yvoje tohoto softwaru nepˇrin´aˇs´ı nadˇeji, ˇze by mohlo doj´ıt ke zmˇen´ am a pˇr´ıpadn´ ym vylepˇsen´ım nebo rozˇs´ıˇren´ı aplikace. Po nabyt´ı znalost´ı uveden´ ych v pˇredchoz´ıch kapitol´ach a z´ısk´an´ı pˇrehledu o funkc´ıch dostupn´ ych n´ astroj˚ u, pˇriˇslo na ˇradu vytvoˇren´ı funkˇcn´ıho prototypu aplikace, kter´a se stane z´akladem pro n´ asledn´ y v´ yvoj a koneˇcn´e nahrazen´ı aplikace Vicon Polygon. Aplikace dostala n´azev Gait Data Viewer, d´ ale v textu zmiˇ nov´ana jako GDV.
6.1
Rozbor v´ ystupn´ıho form´ atu biometrick´ ych dat
Existuje nˇekolik zp˚ usob˚ u jak uchov´avat nasn´ıman´a data motion capturov´ ymi syst´emy. Tom´ aˇs Sol´ar ve sv´e diplomov´e pr´ aci [9] zmiˇ nuje moˇzn´e varianty. Mezi nejrozˇs´ıˇrenˇejˇs´ı patˇr´ı form´ aty BVA, BVH, ACM nebo C3D. Vicon Nexus vˇsak poskytuje uchov´an´ı dat pouze v C3D form´atu. C3D (Coordinate 3D) form´at je podmnoˇzinou v´ıce obecn´eho ADTech souborov´eho form´atu. Byl navrˇzen Andrewem Dainisem k uchov´an´ı 3D koordin´at a analogov´ ych dat spoleˇcnˇe s informacemi popisuj´ıc´ımi takto uloˇzen´a data. To vˇse v jednom souboru. Schopnost ukl´ adat velk´e mnoˇzstv´ı informac´ı o datech, je v´ yhoda, kter´a odliˇsuje C3D form´at od ostatn´ıch biomechanick´ ych form´ at˚ u. V d˚ usledku toho je v C3D souboru mal´e mnoˇzstv´ı obecnˇe spoleˇcn´ ych parametr˚ u, kter´e popisuj´ı 3D data a pak libovoln´e mnoˇzstv´ı uˇzivatelem nebo laboratoˇr´ı definovan´ ych, vytvoˇren´ ych a uloˇzen´ ych datov´ ych poloˇzek. Vzhledem k tomu, ˇze C3D form´ at nen´ı v´ az´ an na konkr´etn´ıho v´ yrobce, m˚ uˇzeme libovolnˇe pˇrid´avat nov´e parametry a pˇrizp˚ usobit ukl´ ad´ an´ı informace podle naˇsich potˇreb. V dneˇsn´ı dobˇe nach´ az´ı C3D form´at uplatnˇen´ı pr´avˇe v takov´ ych laboratoˇr´ıch, jako je zmiˇ novan´ a laboratoˇr ch˚ uze v Brnˇe. Form´at m´a tak´e ˇsirok´e vyuˇzit´ı i v jin´ ych oblastech, ˇ pˇredevˇs´ım z´ abavn´ım a filmov´em pr˚ umyslu. Rada 3D model´aˇr˚ u podporuje jeho naˇcten´ı. Jmenujme napˇr´ıklad komerˇcn´ı Autodesk Maya, Autodesk 3ds Max nebo open source model´ aˇr Blender.
37
38
6.1.1
ˇ NAME ˇ REN ˇ ´ ˇ CHUZE ˚ KAPITOLA 6. ZOBRAZOVAC YCH DAT LABORATORE
Z´ akladn´ı struktura C3D souboru
V uˇzivatelsk´e pˇr´ıruˇcce [7] se uv´ ad´ı, ˇze C3D form´at m´a tˇri z´akladn´ı ˇc´asti: Data - na t´eto u ´rovni je C3D soubor jednoduch´ y bin´arn´ı soubor, kter´ y obsahuje surov´e 3D a analogov´e informace. Standard Parameters - jedn´ a se o implicitn´ı informace o surov´ ych 3D a analogov´ ych datech, kter´e jsou potˇrebn´e k pˇr´ıstupu k u ´daj˚ um. Custom Parameters - zde jsou informace specifikuj´ıc´ı konkr´etn´ı software, kter´ y C3D soubor vytvoˇril, pˇr´ıpadnˇe subjekt, kter´ y byl sn´ım´an. Jedn´ım z c´ıl˚ u pˇri n´ avrhu C3D form´ atu bylo udˇelat ho jednoduch´ y pro listov´an´ı, prohled´av´an´ı popˇr´ıpadˇe modifikov´ an´ı parametr˚ u v C3D souboru. D´ıky tˇemto moˇznostem vynik´a C3D form´at schopnost´ı uchov´ avat velk´e mnoˇzstv´ı informace o datech. C3D form´at tak kombinuje funkce pro ukl´ ad´ an´ı dat a z´ aroveˇ n nese charakteristiku datab´aze. Dalˇs´ım c´ılem pro C3D form´ at byla minimalizace poˇzadavk˚ u na u ´loˇzn´ı prostor, minimalizace poˇctu soubor˚ u potˇrebn´ ych k uchov´an´ı dat a poskytnut´ı dostateˇcn´e rychlosti a jednoduch´eho pˇr´ıstupu k obsahu souboru. C3D form´at umoˇzn ˇuje opatˇrit kaˇzd´ y parametr unik´atn´ım oznaˇcen´ım tak, aby jednotliv´e funkce kaˇzd´eho parametru mohli b´ yt dokumentov´any v samotn´em C3D souboru. To poskytuje flexibilitu pro uˇzivatele, kter´ y tak m˚ uˇze uchov´avat velk´e mnoˇzstv´ı rozd´ıln´ ych dat a souvisej´ıc´ıch parametr˚ u v r´amci jednoho C3D souboru, pˇri zachov´ an´ı m´ıry intern´ı dokumentace, kter´a chyb´ı ve vˇetˇsinˇe ostatn´ıch form´at˚ u. C3D form´at nem´ a ˇz´ adn´ a limituj´ıc´ı krit´eria na ukl´adan´a data, typy nebo poˇcet parametr˚ u. Form´at je lehce rozˇsiˇriteln´ y pro pˇrid´ an´ı dalˇs´ıch parametr˚ u a dat. Soubor C3D je v bin´ arn´ım form´ atu. Bin´arn´ı soubory, oproti ASCII soubor˚ um (textov´ ym), jsou efektivnˇejˇs´ı z hlediska ukl´ ad´ an´ı dat a pˇr´ıstupu k obsahu. Nicm´enˇe bin´arn´ı pˇr´ıstup k dat˚ um obecnˇe vyˇzaduje specifickou aplikaci, kter´a nese detailn´ı znalosti o struktuˇre souboru a o tom jak jsou hodnoty dat v souboru uloˇzeny. Navzdory t´eto komplikaci, je u ´ˇcinnost a rychlost pˇr´ıstupu k bin´ arn´ımu souboru jednoznaˇcnˇe pˇrevaˇzuj´ıc´ı v´ yhodou. C3D form´at se d´ıky tomu stal jasnou volbou pro uchov´ an´ı biomechanick´ ych dat a parametr˚ u. Form´at C3D zach´ az´ı s informac´ı, jako by patˇrila do jedn´e ze dvou tˇr´ıd: Fyzik´ aln´ı mˇeˇren´ı a Parametry, kter´e se t´ ykaj´ı fyzik´ aln´ıho mˇeˇren´ı. Fyzik´ aln´ı mˇ eˇ ren´ı C3D specifikace poˇc´ıt´ a za fyzik´ aln´ı mˇeˇren´ı bud’ poziˇcn´ı informaci (3D souˇradnice) nebo numerick´a data (digitalizovan´ a analogov´ a informace). Kaˇzd´ y 3D koordin´ at je uloˇzen jako surov´ y (raw) X, Y a Z z´aznam s informac´ı o vzorku (tj. pˇresnost a pˇr´ıspˇevek jednotliv´ ych kamer). Kaˇzd´ y vzorek numerick´ ych dat m˚ uˇze obsahovat digitalizovanou analogovou informaci ze zdroje jako je EMG a tlakov´ a ploˇsina. Tyto vzorky jsou synchronizov´any se vzorky 3D souˇradnic, takˇze je snadn´e urˇcit ˇc´ıseln´e hodnoty dat, kter´e se vztahuj´ı ke kaˇzd´e souˇradnici v r´amci souboru.
´ ´IHO FORMATU ´ ´ 6.1. ROZBOR VYSTUPN BIOMETRICKYCH DAT
39
V d˚ usledku toho, mnoho C3D soubor˚ u obsahuje analogov´a a 3D data synchronizov´ ana sn´ımek po sn´ımku. To je velk´ a v´ yhoda oproti ostatn´ım form´at˚ um, kter´e obvykle uchov´ avaj´ı informace oddˇelenˇe ve v´ıce souborech. Uchov´av´an´ı vˇsech informac´ı v jednom souboru d´ av´ a vˇetˇs´ı d˚ uvˇeru v data. Je snazˇs´ı z´ıskat relevantn´ı data a takov´eto uchov´an´ı dat d´av´a jistotu, ˇze data z v´ıce zdroj˚ u, jako jsou motion capturov´e kamery, tlakov´e ploˇsiny a EMG zaˇr´ızen´ı jsou synchronizov´ ana v ˇcase a 3D prostoru. Parametry Kromˇe 3D a numerick´e informace obsahuje C3D soubor tak´e informace o datech, jako mˇern´e jednotky, popisy marker˚ u atd. Nicm´enˇe, na rozd´ıl od jin´ ych form´at˚ u, umoˇzn ˇuje C3D form´ at uloˇzen´ı takov´e informace jako jm´eno pacienta, diagn´ozy a dalˇs´ı poloˇzky, kter´e mohou b´ yt specifick´e pro n´ asledn´e vyhodnocen´ı.
6.1.2
Popis C3D souboru
Data v C3D souboru jsou uloˇzena v 16-bitov´em signed integer form´atu nebo volitelnˇe ve form´atu s plovouc´ı desetinou ˇc´ arkou. Form´at dat m˚ uˇze b´ yt urˇcen v hlaviˇcce souboru na bin´arn´ı u ´rovni bez vytv´ aˇren´ı jak´ ychkoliv dohad˚ u o form´atu dat. Pro kompatibilitu s r˚ uzn´ ymi operaˇcn´ımi syst´emy by mˇeli b´ yt vˇsechny C3D soubory povaˇzov´any za konzistentn´ı bloky, kter´e maj´ı d´elku 512 byt˚ u (nebo 256 byt˚ u pro 16-bit slova). Kaˇzd´ y C3D soubor mus´ı obsahovat minim´alnˇe tˇri sekce informac´ı: Sekce hlaviˇcky (jeden 512 bytov´ y blok). Sekce parametr˚ u (jeden nebo v´ıce 512 bytov´ ych blok˚ u). Sekce 3D dat (jeden nebo v´ıce 512 bytov´ ych blok˚ u). Tabulka 6.1: Z´akladn´ı struktura C3D souboru.
Hlaviˇ cka souboru (Header section) Prvn´ı sekc´ı C3D souboru je hlaviˇcka, kter´a vˇzdy zaˇc´ın´a blokem 1 (prvn´ı blok souboru). Prvn´ı slovo C3D souboru urˇcuje pozici sekce parametr˚ u, kter´a zase obsahuje ukazatel na zaˇc´ atek sekce 3D dat, stejnˇe jako detailn´ı informace nezbytn´e k ˇcten´ı sekce s 3D daty a interpretaci obsahu. Nav´ıc hlaviˇcka d´ ale obsahuje duplicitn´ı z´aznamy z ˇrad parametr˚ u uloˇzen´ ych v sekci parametr˚ u a omezen´e mnoˇzstv´ı o ud´alostech. Pozice sekce 3D dat, stejnˇe jako ostatn´ı d˚ uleˇzit´e atributy, by mˇely b´ yt vˇzdy naˇc´ıt´any ze sekce parametr˚ u v C3D souboru. D˚ uvod proˇc pozice sekce s 3D daty a ostatn´ıch parametr˚ u je duplicitnˇe v hlaviˇcce souboru je ten, Obr´azek 6.1: Obsah hlaviˇcky ˇze umoˇzn ˇuje jednoduch´ ym utilit´ am pˇristupovat k 3D dat˚ um vypsan´ y aplikac´ı VBC3DEditor. bez nutnosti naˇc´ıt´ an´ı a dek´ odov´ an´ı cel´e sekce parametr˚ u.
40
ˇ NAME ˇ REN ˇ ´ ˇ CHUZE ˚ KAPITOLA 6. ZOBRAZOVAC YCH DAT LABORATORE
Sekce parametr˚ u (Parametr section) Sekce parametr˚ u obvykle zaˇc´ın´a v C3D souboru blokem ˇc´ıslo 2, avˇsak nen´ı to nutnost´ı a mus´ıme poˇc´ıtat s t´ım, ˇze ne kaˇzd´ y C3D soubor je takto stavˇen. C3D specifikace vyˇzaduje aby sekce parametr˚ u zaˇc´ınala na zaˇc´atku 512ti bytov´em bloku, urˇcen´em ukazatelem v hlaviˇcce souboru. Sekce parametr˚ u m˚ uˇze b´ yt r˚ uznˇe dlouh´a, ale obvykle m´a d´elku osm aˇz deset blok˚ u. Na obr´azku 6.2 je vidˇet zv´ yraznˇen´a sekce parametr˚ u C3D souboru vyexportovan´eho Vicon Nexusem. D´ale si m˚ uˇzete vˇsimnout, ˇze je vybr´ana skupina POINT, kter´a popisuje 3D data. V prav´e ˇc´asti obr´azku stoj´ı za zm´ınku seznam popsan´ ych marker˚ u (at’ uˇz fyzick´ ych nebo virtu´aln´ıch) a jim pˇriˇrazen´ a unik´atn´ı ˇc´ısla. Obr´azek 6.2: Sekce parametru vypsan´a aplikac´ı VBC3DEditor.
Sekce 3D dat (3D data section) C3D soubor obsahuje 3D souˇradnice marker˚ u a analogov´a data v sekci, kter´a je obvykle um´ıstˇena na nˇejak´em m´ıstˇe za sekc´ı parametr˚ u. Specifikace C3D ˇr´ık´ a, ˇze sekce 3D dat zaˇc´ın´a 512 bytov´ ym blokem, na kter´ y ukazuje ukazatel POINT:DATA START v C3D parametrick´e sekci. 3D i analogov´a data maj´ı promˇenlivou d´elku podle mnoˇzstv´ı uloˇzen´ ych dat. C3D soubor m˚ uˇze obsahovat kombinaci 3D dat a analogov´e informace. Stejnˇe tak vˇsak m˚ uˇze obsahovat pouze jednu z nich. Na obr´azku 6.3 je vidˇet sekce 3D dat, obsahuj´ıc´ı koordin´aty jednotliv´ ych marker˚ u. V prav´e ˇc´asti obr´azku jsou vidˇet X, Y a Z souˇradnice ym sn´ımk˚ um. Obr´azek 6.3: Sekce 3D dat vypsan´ a ap- odpov´ıdaj´ıc´ı jednotliv´ likac´ı VBC3DEditor.
6.2
Knihovna C3DServer SDK
Aby program´atoˇri mohli jednoduˇse pracovat s daty uloˇzen´ ymi v bin´arn´ım C3D souboru, bylo zapotˇreb´ı vytvoˇrit n´ astroj umoˇzn ˇuj´ıc´ı k nim pˇristupovat. Takov´ ym n´astrojem je C3DServer Software Development Kit (SDK) vyvinut´ y firmou Motion Lab Systems. C3DServer umoˇzn ˇuje pˇristupovat k dat˚ um v C3D souboru a pˇr´ıpadnˇe je editovat. C3DServer je navrˇzen´ y k pr´aci na 32-bitov´em operaˇcn´ım syst´emu Microsoft Windows prostˇrednictv´ım technologie COM,
6.3. IMPLEMENTACE
41
pˇres jednoduch´e, dobˇre definovan´e rozhran´ı. C3DServer eliminuje potˇrebu program´ator˚ u, ps´at vlastn´ı k´ od podporuj´ıc´ı C3D form´at na bin´arn´ı u ´rovni. V referenˇcn´ım manu´ alu [8] se p´ıˇse, ˇze jeho pouˇzit´ı je moˇzn´e v programech psan´ ych v jazyce Java, C/C++ ˇci Visual Basic. Pˇri implementaci GDV vˇsak nebyl ˇz´adn´ y probl´em vyuˇz´ıt C3DServer SDK v aplikaci psan´e v jazyce C# a to d´ıky technologii COM od spoleˇcnosti Microsoft. D´ıky funkc´ım kter´e C3DServer poskytuje, se program´ator m˚ uˇze soustˇredit pouze na zpracov´ an´ı dat z C3D souboru a nemus´ı ˇreˇsit jak jsou data v souboru uchov´ana. Jeˇstˇe je tˇreba zm´ınit, ˇze C3DServer SDK existuje ve dvou verz´ıch a to pro komerˇcn´ı a nekomerˇcn´ı vyuˇzit´ı. Implementace GDV vyuˇz´ıv´a nekomerˇcn´ı verze, kter´a se od komerˇcn´ı verze liˇs´ı pouze rychlost´ı vykon´ an´ı pˇr´ıkaz˚ u. Tom´aˇs Sol´ar ve sv´e pr´aci [9] p´ıˇse, ˇze se jedn´ ao 5ms zpoˇzdˇen´ı pˇri vol´ an´ı kaˇzd´e funkce.
6.3
Implementace
Pro implementaci GDV byla zvolena technologie Microsoft XNA, spoleˇcnˇe s programovac´ım jazykem C#. Microsoft XNA je sada n´astroj˚ u, kter´a je urˇcen´a k v´ yvoji poˇc´ıtaˇcov´ ych her pro menˇs´ı v´ yvojov´e t´ ymy. Programuj´ı-li se hry napˇr´ıklad pomoc´ı technologie OpenGL nebo Microsoft DirectX, mus´ı se napsat kv˚ uli ban´aln´ım vˇecem pomˇernˇe velk´e mnoˇzstv´ı k´ odu. Framework XNA je postaven nad technologi´ı DirectX a obsahuje spoustu funkc´ı, kter´e by program´ atoˇri jinak museli sami programovat. V´ yvoj´aˇri her se tak mohou soustˇredit na samotnou hern´ı logiku, tedy na to, co m´a hra dˇelat ve skuteˇcnosti. Pˇritom se nemus´ı zaj´ımat na jak´em hardwaru hra pobˇeˇz´ı. Hry napsan´e v XNA mohou bˇeˇzet na operaˇcn´ım syt´emu Windows XP a novˇejˇs´ım, pˇr´ıpadnˇe na platformˇe Xbox 360. N´astroje kter´e XNA poskytuje v´ yvoj´aˇr˚ um her, se osvˇedˇcili i pˇri v´ yvoji GDV. Co je zapotˇreb´ı k tomu, aby v˚ ubec tento framework mohl b´ yt pouˇz´ıv´an? Framework XNA [15] je obsaˇzen ve v´ yvojov´em bal´ıku XNA Game Studio. Nejnovˇejˇs´ı verze nese oznaˇcen´ı 3.1 a je kompatibiln´ı i s nejnovˇejˇs´ım operaˇcn´ım syst´emem od Microsoftu, Windows 7. K jeho u ´spˇeˇsn´e instalaci je zapotˇreb´ı m´ıt nainstalovan´e Microsoft Visual C# 2008 Express Edition, jehoˇz souˇc´ast´ı je i Microsoft .NET Framework 3.5.
Pozn´ amka k implementaci V r´amci implementace GDV bylo zvoleno jen zobrazov´an´ı fyzick´ ych a virtu´aln´ıch marker˚ u, kter´e pˇredstavuj´ı klouby sn´ıman´eho pacienta. Jen pro pˇripomenut´ı, v´ ypoˇcet pozice tˇechto marker˚ u m´ a na starosti Plug-in Gait, popsan´ y v sekci 4.1. Prvotn´ı snahou bylo vykreslovat model kostry, jej´ıˇz pohyb by byl d´an znalost´ı trajektorie kloub˚ u z C3D souboru. Takov´ yto typ u ´lohy se ˇreˇs´ı pomoc´ı inversn´ı kinematiky. Ta se vˇsak nepodaˇrila zav´est do v´ ysledn´eho k´odu aplikace v d˚ usledku nedostatku hlubˇs´ı znalosti a tak´e nezkuˇsenosti v oblasti implementace. V r´amci ˇcasov´eho harmonogramu nezbyl ˇcas pro doplnˇen´ı t´eto znalosti a tak byla zvolena alternativn´ı cesta. Kaˇzd´a kost je tedy realizov´ ana pomoc´ı pˇr´ımky mezi virtu´ aln´ımi markery pˇredstavuj´ıc´ımi klouby. V´ ysledkem t´eto pr´ ace je tedy vizualizace fyzick´ ych a virtu´aln´ıch marker˚ u zn´azorˇ nuj´ıc´ıch klouby a naznaˇcen´ı kost´ı. Pro pˇredstavu n´asleduje sn´ımek (6.4) z fin´aln´ı aplikace Gait Data Viewer.
42
ˇ NAME ˇ REN ˇ ´ ˇ CHUZE ˚ KAPITOLA 6. ZOBRAZOVAC YCH DAT LABORATORE
Obr´ azek 6.4: Sn´ımek obrazovky aplikace GDV pro vizualizace dat.
6.3.1
Z´ aklad aplikace
V t´eto ˇc´asti jsou pops´ any z´ akladn´ı metody, kter´e XNA framework vyˇzaduje v kaˇzd´e sv´e aplikaci. Pˇritom bude rovnou pops´ ano co v aplikaci GDV obsahuj´ı. V konstruktoru tˇr´ıdy GaitDataViewer se vytv´ aˇr´ı struktury GraphicsDeviceManager (handle pro konfiguraci a spr´avu grafick´ ych zaˇr´ızen´ı) a ContentManager (zajiˇst’uje spr´avu naˇc´ıt´an´ı extern´ıch bin´arn´ıch soubor˚ u). Po samotn´em spuˇstˇen´ı aplikace jsou po konstruktoru jeˇstˇe vol´any dvˇe funkce. Prvn´ı funkce Initialize() inicializuje glob´ aln´ı promˇenn´e jako transformaˇcn´ı matice kamery atd. Tou druhou metodou je LoadContent(), kde se naˇc´ıtaj´ı modely z extern´ıch soubor˚ u, jako je model retroreflexn´ıho markeru. Cel´a myˇslenka frameworku XNA je zaloˇzena na metod´ach Update() a Draw(). Tyto metody jsou pˇri spuˇstˇen´e aplikaci vol´ any ve smyˇcce. D´ıky tomu mohou b´ yt v metodˇe Update() neust´ale aktualizov´ any vˇsechny v´ ypoˇcty a zpracov´any vstupy od uˇzivatele. V metodˇe Draw() jsou pak aktualizovan´e promˇenn´e prom´ıtnuty do fin´aln´ıho vykreslen´ı sn´ımku na obrazovku.
6.3.2
Naˇ cten´ı C3D souboru
Pro naˇc´ıt´an´ı bin´ arn´ıch C3D soubor˚ u byl vyuˇzit C3DServer SDK, popsan´ y v pˇredch´azej´ıc´ı ˇc´asti 6.2. Naˇc´ıt´ an´ı i pr´ ace s C3D souborem prob´ıh´a ve tˇr´ıdˇe C3dData. Zde se naˇc´ıt´a tak´e .vsk soubor (5). Jedn´ a se o XML soubor, ve kter´em jsou mino jin´e uvedeny n´azvy fyzick´ ych
6.3. IMPLEMENTACE
43
marker˚ u a jejich barva, kterou maj´ı b´ yt zv´ yraznˇeny. D´ıky znalosti tˇechto marker˚ u mohou b´ yt vytˇeˇzeny z C3D souboru jejich pozice v jednotliv´ ych sn´ımc´ıch animace a pot´e b´ yt vykresleny. D´ale bylo zapotˇreb´ı z´ıskat n´ azvy virtu´aln´ıch marker˚ u pˇredstavuj´ıc´ı klouby. Tyto markery jsou spoleˇcnˇe s velk´ ym poˇctem dalˇs´ıch uvedeny v .mkr souboru (5). Jedn´a se o textov´ y soubor, kter´ y vˇsak nem´ a XML strukturu. Je zde vˇzdy ˇr´adek po ˇr´adku uveden n´azev markeru s jeho popisem, pˇr´ıpadnˇe do jak´eho patˇr´ı segmentu. Musel by b´ yt naimplementov´an vlastn´ı zp˚ usob naˇc´ıt´an´ı takov´ehoto souboru, k ˇcemuˇz nebyl prostor. Z toho d˚ uvodu byla vytvoˇrena funkce kde se pozice tˇechto marker˚ u zjiˇst’uj´ı podle jejich ruˇcnˇe vypsan´ ych n´azv˚ u. Ano, nejedn´a se o u ´plnˇe nejlepˇs´ı ˇreˇsen´ı, ale pro v´ ysledn´ y prototyp je to v tuto chv´ıli postaˇcuj´ıc´ı. Jeˇstˇe mus´ı b´ yt zm´ınˇena jedna ne zcela standardn´ı z´aleˇzitost, kter´a se v tˇr´ıdˇe C3dData objevuje. Jelikoˇz aplikace GDV pouˇz´ıv´a Free verzi C3DServeru, je start aplikace dosti zdˇ adovˇe 2-3 minuty. Z´ louhav´ y. R´ aleˇz´ı ovˇsem na velikosti C3D souboru. Je to zp˚ usobeno t´ım, ˇze pˇri startu aplikace si do jednorozmˇern´eho pole aplikace ukl´ad´a pozice vˇsech marker˚ u ve vˇsech sn´ımc´ıch z´ aznamu. Pˇritom vyuˇz´ıv´a funkc´ı z SDK C3DServer a jak bylo zm´ınˇeno dˇr´ıve, daˇ n za bezplatnou verzi je pomal´ a rychlost ˇcten´ı dat z C3D souboru. Obch´az´ı se t´ım vˇsak vol´ an´ı tˇechto metod pˇri samotn´e vizualizaci, kter´a by d´ıky t´eto vlastnosti neodpov´ıdala poˇr´ızen´emu z´aznamu.
6.3.3
Zobrazovac´ı smyˇ cka
Jak bylo dˇr´ıve zm´ınˇeno, zobrazovac´ı smyˇcka je realizov´ana metodou Draw(). V t´eto metodˇe se vol´a metoda pro vykreslen´ı pˇr´ımek zn´azorˇ nuj´ıc´ıch kostru. Kaˇzd´a kost je zakonˇcena virtu´ aln´ım kloubem. Mezi klouby pak staˇc´ı vykreslit pˇr´ımku pomoc´ı n´asleduj´ıc´ıho pˇr´ıkazu: game.GraphicsDevice.DrawUserPrimitives
( PrimitiveType.LineList, // definuje jaky typ primitiv se bude vykreslovat vertices.ToArray(), // pole pozic kloubu 0, // zacneme prvnim kloubem v poli vertices.Count / 2); // skoncime poslednim Pot´e je na ˇradˇe vykreslen´ı vˇsech marker˚ u, at’ uˇz fyzick´ ych nebo virtu´aln´ıch. Pro vykreslen´ı markeru se pouˇz´ıv´ a model koule vyexportovan´ y z model´aˇre Autodesk Maya do form´atu FBX. Pro naˇcten´e modely z extern´ıho souboru byla vytvoˇrena tˇr´ıda ModelObject, kter´a v sobˇe uchov´av´ a informace specifick´e pro kaˇzd´ y objekt. Je zde metoda pro vykreslen´ı modelu, kter´ a je vol´ana z hlavn´ı smyˇcky Draw(). Zobrazen´ı modelu realizuje n´asleduj´ıc´ı k´od: public void DrawModelObject() { foreach (ModelMesh mesh in model.Meshes) { foreach (BasicEffect effect in mesh.Effects) { // povoleni defaultniho osvetleni effect.EnableDefaultLighting();
ˇ NAME ˇ REN ˇ ´ ˇ CHUZE ˚ KAPITOLA 6. ZOBRAZOVAC YCH DAT LABORATORE
44
// pokud to graficke zarizeni umoznuje, // melo by byt pouzito per-pixel osvetleni effect.PreferPerPixelLighting = true; effect.World = Matrix.CreateScale(scale) * Matrix.CreateTranslation(modelPosition);// effect.View = cameraViewMatrix; // effect.Projection = cameraProjectionMatrix; // if (this.isDiffuseChange) effect.DiffuseColor = diffuseColor; // } mesh.Draw();
svetova matice pohledova matice projekcni matice nastaveni difusni barvy
// vykreslit mesh
} } Protoˇze k struktur´ am marker˚ u je tˇreba pˇristupovat z v´ıce m´ıst v programu, byla vytvoˇrena tˇr´ıda podle n´avrhov´eho vzoru singleton. Motivac´ı t´eto tˇr´ıdy je uˇzit´ı pouze jedn´e instance v cel´em programu. Nach´ az´ı se zde datov´ a struktura Dictionary. Prvky vkl´adan´e do tohoto slovn´ıku jsou typu ModelObject a jsou identifikov´any kl´ıˇcem, jehoˇz hodnota odpov´ıd´a indexu markeru v C3D souboru. Pˇri vykreslov´an´ı jednotliv´ ych marker˚ u z hlavn´ı Draw() metody, se prolistuje cel´ y slovn´ık a nad kaˇzd´ ym markerem se zavol´a dˇr´ıve uveden´a metoda DrawModelObject(). D´ ale je uvedena ˇc´ast k´odu hlavn´ı metody Draw(): // vykresleni kosti this.boneLines.DrawBones(cameraViewMatrix, cameraProjectionMatrix); // vykreselni vsech modelu ktere chceme vykreslit foreach (int i in DrawableModelsSingleton.Instance.ModelsDictionary.Keys) { DrawableModelsSingleton.Instance.ModelsDictionary[i]. CameraProjectionMatrix = cameraProjectionMatrix; DrawableModelsSingleton.Instance.ModelsDictionary[i]. CameraViewMatrix = cameraViewMatrix; DrawableModelsSingleton.Instance.ModelsDictionary[i].DrawModelObject(); }
Kapitola 7
Testov´ an´ı aplikace na re´ aln´ ych datech Naˇ cten´ı dat Aplikace Gait Data Viewer byla otestov´ana na mnoˇzstv´ı vybran´ ych z´aznam˚ u, kter´e byly poˇr´ızeny bˇehem instalace syst´em˚ u v brnˇensk´e laboratoˇri. GDV aplikaci byl vˇzdy pˇredloˇzen C3D soubor s odpov´ıdaj´ıc´ım .vsk souborem. Jedn´a se o soubory, kter´e generuje Vicon Nexus (viz. sekce 5.1.4). Nejprve jsou naˇcteny jm´ena fyzick´ ych marker˚ u z .vsk souboru. D´ıky znalosti jmen tˇechto marker˚ u pak mohou b´ yt naˇcteny z C3D souboru jejich pozice v jednotliv´ ych sn´ımc´ıch animace. Pot´e uˇz staˇc´ı naˇcten´a data pouze zobrazit.
Zobrazen´ı dat Aplikace GDV zobrazuje, trajektorie fyzick´ ych marker˚ u pˇripevnˇen´ ych k pacientovi bˇehem prov´adˇen´ı dynamick´eho z´ aznamu. D´ale zobrazuje virtu´aln´ı markery pˇredstavuj´ıc´ı klouby, urˇcen´e na z´ akladˇe OFM 2.1.4. Mezi virtu´aln´ımi markery jsou vykresleny pˇr´ımky zn´azorˇ nuj´ıc´ı kosti segment˚ u tˇela. D´elka trv´ an´ı i rychlost pˇrehr´av´an´ı animace byla porovn´ana s v´ ysledky, kter´e prezentuje Vicon Polygon a byly shled´any za totoˇzn´e.
Funkˇ cnost aplikace Aplikace GDV m˚ uˇze b´ yt povaˇzov´ ana za stabiln´ı. Jej´ı funkˇcnost nen´ı nijak limitov´ana velikost´ı vstupn´ıho C3D souboru. Jedin´ a pot´ıˇz je v pomal´em startu programu, kdy se naˇc´ıt´a informace z C3D souboru. To je ovˇsem zp˚ usobeno nekomerˇcn´ı verz´ı knihovny C3DServer SDK, kter´ a je limitov´ ana rychlost´ı prov´ adˇen´ ych funkc´ı. Start aplikace tak zabere ˇr´adovˇe 2 aˇz 5 min v z´avislosti na velikosti vstupn´ıho C3D souboru.
45
46
´ ´I APLIKACE NA REALN ´ YCH ´ KAPITOLA 7. TESTOVAN DATECH
Kapitola 8
Z´ avˇ er V reˇserˇsn´ı ˇc´ asti bakal´ aˇrsk´e pr´ ace bylo pops´ano jak se urˇcuje v´ ysledn´a pozice doln´ı poloviny ˇ aˇr byl sezn´amen s hardwarov´ lidsk´eho tˇela v sn´ıman´e oblasti. Cten´ ymi a softwarov´ ymi syst´emy k zaznamen´ an´ı pacientovi ch˚ uze, pouˇz´ıvan´e v laboratoˇri ch˚ uze ve Fakultn´ı dˇetsk´e nemocnici J. G. Mendela v Brnˇe. Byl zm´ınˇen postup sn´ım´an´ı dat v laboratoˇri a rozebr´any souborov´e form´aty, kter´e slouˇz´ı k uchov´ an´ı nasn´ıman´ ych dat. Poznatky o st´avaj´ıc´ıch syst´emech a souˇcasn´ ych moˇznostech zobrazov´an´ı dat, byly vyuˇzity v z´avˇereˇcn´e implementaci.
Zhodnocen´ı pr´ ace Problematika lidsk´e ch˚ uze spoleˇcnˇe s jej´ım zaznamen´an´ım, vizualizac´ı a vyhodnocov´an´ım je velmi sloˇzit´ a a obs´ ahl´ a. D˚ ukazem toho je mnohalet´ y v´ yzkum v t´eto oblasti, kter´ y zapoˇcal v 80-t´ ych letech 20. stolet´ı. D´ıky t´eto pr´aci jsem nahl´edl do problematiky anal´ yzy lidsk´e ch˚ uze. Sezn´ amil jsem se s prostˇred´ım laboratoˇre ch˚ uze a z´ıskal pˇrehled o moˇznostech sn´ım´ an´ı ´ eˇsnˇe jsem implementoval prototyp aplikace Gait Data Viewer, kter´a umoˇzn pohybu. Uspˇ ˇuje vizualizaci nasn´ıman´ ych dat. Funkˇcnost aplikace jsem ovˇeˇril naˇcten´ım a zobrazen´ı re´aln´ ych dat, namˇeˇren´ ych v laboratoˇri ch˚ uze. Bohuˇzel se do posledn´ı chv´ıle nepodaˇrilo pˇresvˇedˇcit veden´ı brnˇensk´e laboratoˇre ch˚ uze o prospˇeˇsnosti t´eto pr´ ace. V d˚ usledku toho mi nebyl povolen pˇr´ıstup do laboratoˇre a nemohl jsem se proto sezn´ amit se syst´emy v provozu. Nemohl jsem ani konzultovat pˇr´ıpadn´e nejasnosti s obsluhuj´ıc´ım person´ alem, pˇr´ıpadnˇe samotn´ ymi l´ekaˇri. Byl jsem odk´az´an pouze na dostupnou literaturu. Pˇresto si mysl´ım, ˇze tento handicap kvalitu pr´ace neovlivnil, pouze pozdrˇzel proces v´ yvoje, kdy jsem si musel z v´ıce zdroj˚ u doplnit znalosti, pˇredevˇs´ım z oblasti anatomie. Z v´ yˇse uveden´ ych d˚ uvod˚ u tak nezbylo mnoho ˇcasu na v´ yslednou implementaci aplikace GDV. Aplikaci, kterou spoleˇcnˇe s t´ımto textem pˇredkl´ad´am, jsem se proto snaˇzil koncipovat tak, aby se mohlo pokraˇcovat v jej´ım v´ yvoji. Mrz´ı mˇe, ˇze se nepodaˇrila naimplementovat inversn´ı kinematika modelu kostry, kter´a by jeˇstˇe l´epe ilustrovala skuteˇcn´ y pohyb pacienta. I pˇresto se vˇsak moje ˇreˇsen´ı bl´ıˇz´ı poˇzadovan´emu v´ ysledku.
47
48
´ ER ˇ KAPITOLA 8. ZAV
Budouc´ı vylepˇ sen´ı V r´amci t´eto pr´ ace byl poloˇzen z´ aklad aplikace pro vizualizaci biometrick´ ych dat. V dalˇs´ı f´azi v´ yvoje je v pl´ anu implementace kostry pomoc´ı inversn´ı kinematiky. D´ale bych chtˇel plnˇe vyuˇz´ıt vˇsech informac´ı obsaˇzen´ ych v souborech produkovan´ ych Vicon Nexusem. Jako zaj´ımav´a myˇslenka se mi jev´ı prezentace z´avˇereˇcn´e zpr´avy na webu. L´ekaˇri by tak nemuseli analyzovat data v prostor´ ach nemocnice, ale napˇr´ıklad v pohodl´ı domova nebo na jak´emkoliv jin´em m´ıstˇe. Dalˇs´ı moˇznost prezentace je prostˇrednictv´ım mobiln´ıho telefonu. Zn´ı to moˇzn´a troˇsku jako sci-fi, ale mobiln´ı telefony s operaˇcn´ı syst´emem Windows Phone 7 Series podporuj´ı technologii XNA. Takˇze cesty pro v´ yvoj jsou otevˇren´e.
Literatura [1] Vicon User Group Meeting - THE OXFORD FOOT MODEL. Motion presentation, 2008. [2] Fakultn´ı nemocnice brno, Naposledy navˇst´ıveno 27.5.2010. http://www.fnbrno.cz. [3] TeleMyo 2400R G2 Analog Output Receiver, Naposledy navˇst´ıveno 27.5.2010. http://www.noraxon.com/downloads/documents/instruments/ TeleMyo-2400R-Receiver-G2.pdf. [4] blesk.cz. Veronika m´ a mozkovou obrnu: Pom˚ uˇze j´ı nov´a laboratoˇr?, Naposledy navˇst´ıveno 27.5.2010. http://www.blesk.cz/clanek/zpravy-pribehy/133058/ veronika-ma-mozkovou-obrnu-pomuze-ji-nova-laborator.html. [5] R. B. Davis, S. Ounpuu, D. Tyburski, and J. Gage. A gait analysis data collection and reduction technique. Human Movement Sciences, pages 575 – 587, 1991. [6] M. P. Kadaba, H. K. Ramakrishnan, M. E. Wootten, J. Gainey, G. Gorton, and G. V. B. Cochran. Repeatability of Kinematic, Kinetic, and Electromyographic Data in Normal Adult Gait. Jurnal of Orthopaedic Research, pages 849 – 860, 1989. [7] Motion Lab Systems. C3D User Manual, 2008. [8] Motion Lab Systems. User Reference Manual C3DServer, 2010. ˇ e vysok´e uˇcen´ı technick´e, [9] T. Sol´ ar. Diplomov´ a pr´ ace - 3D zobrazovaˇc znakov´e ˇreˇci. Cesk´ Fakulta elektrotechnick´ a, 2007. [10] Vicon Motion Systems. Polygon System Reference, revision 1.2 edition, 2005. [11] Vicon Motion Systems. Vicon MX Hardware System Reference, revision 1.6 (draft 4) edition, 2006. [12] AMTI Multiaxis Force Platform, Naposledy navˇst´ıveno 27.5.2010. http://amti.biz/ Platformhome/PlatformHome.htm. [13] Noraxon EMG and sensor systems, Naposledy navˇst´ıveno 27.5.2010. http://www. noraxon.com. [14] Vicon Peak, Naposledy navˇst´ıveno 27.5.2010. http://www.vicon.com. [15] Xna framework, Naposledy navˇst´ıveno 27.5.2010. http://xna.com.
49
50
LITERATURA
[16] F. E. Zajak. Muscle and tendon: properties, models, scaling and application to biomechanics and motor control, pages 359–411. CRC Crit. Rev. Biomed., 1989. ˇ a Televize. [17] Cesk´ Poˇrad PORT - laboratoˇr ch˚ uze, Naposledy navˇst´ıveno 27.5.2010. http://www.ceskatelevize.cz/program/port/17.02.2010-17:30-2/ tema/528-laborator-chuze/.
Pˇ r´ıloha A
Obrazov´ a pˇ r´ıloha
Obr´ azek A.1: Sn´ımek z v´ ysledn´e aplikace.
51
52
ˇ ´ILOHA A. OBRAZOVA ´ PR ˇ ´ILOHA PR
Obr´ azek A.2: Sn´ımek z v´ ysledn´e aplikace.
Obr´ azek A.3: Sn´ımek z v´ ysledn´e aplikace.
53
Obr´ azek A.4: Prostˇred´ı brnˇensk´e laboratoˇre ch˚ uze. Zdroj: [2]
Obr´azek A.5: Zdravotn´ı sestra Aneta Fedrov´a z brnˇensk´e laboratoˇre ch˚ uze pˇripevˇ nuje na nohy pacientky retroreflexn´ı markery. Zdroj: [4]
54
ˇ ´ILOHA A. OBRAZOVA ´ PR ˇ ´ILOHA PR
Obr´ azek A.6: Prov´ adˇen´ı dynamick´eho z´aznamu. Zdroj: [4]
Obr´ azek A.7: Prov´ adˇen´ı dynamick´eho z´aznamu. Zdroj: [4]
Pˇ r´ıloha B
Seznam pouˇ zit´ ych zkratek A/D Analogovˇe digit´ aln´ı pˇrevodn´ık EMG Elektromyografie LED Light Emitting Diodes USB Universal Serial Bus HJC Hip Joint Centers ASIS Anterior Superior Iliaca Spina KJC Knee Joint Center THI Thigh marker KNE Knee marker AJC Ankle Joint Centre TIB Tibia marker OFM Oxford Foot Model VST Vicon Skeleton Template VSK Vicon Skeleton MKR Vicon Marker File MP Vicon Model Parameters File C3D Coordinate 3D XML Extensible Markup Language
55
56
ˇ ´ILOHA B. SEZNAM POUZIT ˇ YCH ´ PR ZKRATEK
Pˇ r´ıloha C
Obsah pˇ riloˇ zen´ eho CD
Obr´ azek C.1: Seznam pˇriloˇzen´eho CD
57