Tworzenie bazy danych
Stworzyliśmy multimodalny zbiór danych do detekcji upadków i rozpoznawania aktywności ludzkiej, o nazwie UP-Fall Detection21. Dane zostały zebrane w ciągu czterech tygodni w School of Engineering na Universidad Panamericana (Meksyk, Meksyk). Scenariusz testowy wybrano, biorąc pod uwagę następujące wymagania: (a) przestrzeń, w której badani mogli komfortowo i bezpiecznie symulować upadki oraz wykonywać czynności, oraz (b) środowisko wewnętrzne z naturalnym i sztucznym oświetleniem, odpowiednie dla ustawień sensorów multimodalnych.
Dostępne są próbki danych od 17 badanych, którzy wykonali 5 typów upadków oraz 6 różnych prostych czynności w trakcie 3 prób. Wszystkie informacje zebrano za pomocą wewnętrznego systemu akwizycji danych, składającego się z 5 czujników ubieralnych (trójosiowy akcelerometr, żyroskop i czujnik natężenia światła), 1 kasku elektroencefalograficznego, 6 czujników podczerwieni jako czujników otoczenia oraz 2 kamer umieszczonych w punktach widzenia bocznego i przedniego. Rysunek 1 przedstawia schemat rozmieszczenia czujników w środowisku i na ciele. Częstotliwość próbkowania całego zbioru danych wynosi 18 Hz. Baza danych zawiera dwa zestawy danych: skonsolidowany zestaw danych surowych (812 GB) oraz zestaw danych cech (171 GB). Wszystkie bazy danych zostały zapisane w chmurze w celu zapewnienia publicznego dostępu: https://sites.google.com/up.edu.mx/har-up/. Więcej szczegółów dotyczących akwizycji, wstępnego przetwarzania, konsolidacji i przechowywania tej bazy danych, a także informacje o synchronizacji i spójności danych można znaleźć w pracy Martínez-Villaseñor et al.21.
W przypadku tej bazy danych wszystkimi badanymi byli zdrowi młodzi ochotnicy (9 mężczyzn i 8 kobiet) bez żadnych dysfunkcji, w wieku od 18 do 24 lat, ze średnim wzrostem 1,6 m i średnią masą ciała 6,8 kg. Podczas gromadzenia danych badacz odpowiedzialny za kwestie techniczne nadzorował, aby wszystkie czynności były wykonywane przez badanych prawidłowo. Osoby badane wykonały pięć rodzajów upadków, z każdego po 10 s, obejmujących: upadek do przodu z użyciem dłoni (1), upadek do przodu z użyciem kolan (2), upadek w tył (3), siadanie na pustym krześle (4) oraz upadek na bok (5). Wykonano również sześć czynności dnia codziennego trwających po 60 s każda, z wyjątkiem skakania (30 s): chodzenie (6), stanie (7), podnoszenie przedmiotu (8), siedzenie (9), skakanie (10) i leżenie (1). Chociaż symulowane upadki nie mogą odwzorować wszystkich rodzajów upadków w życiu rzeczywistym, istotne jest włączenie przynajmniej reprezentatywnych typów upadków, co umożliwia stworzenie lepszych modeli detekcji upadków. Istotne jest również zastosowanie ADL (czynności dnia codziennego), a w szczególności aktywności, które mogą być zazwyczaj pomylone z upadkami, takich jak podnoszenie przedmiotu. Rodzaje upadków i ADL wybrano po przeanalizowaniu pokrewnych systemów detekcji upadków21. Przykładowo, Rycina 2 przedstawia sekwencję obrazów z jednej próby, w której badany upada na bok.
Wyekstrahowano 12 cech czasowych (średnia, odchylenie standardowe, amplituda maksymalna, amplituda minimalna, wartość skuteczna RMS, mediana, liczba przejść przez zero, skośność, kurtoza, pierwszy kwartyl, trzeci kwartyl oraz autokorelacja) i 6 cech częstotliwościowych (średnia, mediana, entropia, energia, częstotliwość główna oraz centroid widmowy)21 z każdego kanału czujników ubieralnych i środowiskowych, co łącznie dało 756 cech. Obliczono również 40 cech wizualnych21 dla każdej kamery, dotyczących względnego ruchu pikseli pomiędzy dwoma sąsiednimi obrazami w materiałach wideo.
Analiza danych pomiędzy podejściami unimodalnymi a multimodalnymi
Wykorzystując bazę danych UP-Fall Detection, przeanalizowano dane w celu porównania podejść unimodalnych i multimodalnych. W tym celu porównano siedem różnych kombinacji źródeł informacji: wyłącznie czujniki podczerwieni (IR); wyłącznie czujniki ubieralne (IMU); czujniki ubieralne i kask (IMU+EEG); czujniki podczerwieni, ubieralne oraz kask (IR+IMU+EEG); wyłącznie kamery (CAM); czujniki podczerwieni i kamery (IR+CAM); oraz czujniki ubieralne, kask i kamery (IMU+EEG+CAM). Dodatkowo porównano trzy różne rozmiary okien czasowych z 50% nakładaniem: jedną sekundę, dwie sekundy i trzy sekundy. Dla każdego segmentu wybrano najużyteczniejsze cechy, stosując selekcję i ranking cech. Przy zastosowaniu tej strategii wykorzystano jedynie 10 cech na każdą modalność, z wyjątkiem modalności IR, w której użyto 40 cech. Ponadto porównanie przeprowadzono z wykorzystaniem czterech znanych klasyfikatorów uczenia maszynowego: RF, SVM, MLP i KNN. Do trenowania modeli uczenia maszynowego zastosowano 10-krotną walidację krzyżową, z podziałem zbiorów danych na 70% do trenowania i 30% do testowania. Tabela 1 przedstawia wyniki tego benchmarku, raportując najlepszą wydajność uzyskaną dla każdej modalności w zależności od modelu uczenia maszynowego i najlepszej konfiguracji długości okna. Metryki oceny obejmują dokładność (accuracy), precyzję (precision), czułość (sensitivity), swoistość (specificity) oraz wskaźnik F1 (F1-score). Rysunek 3 przedstawia te wyniki w formie graficznej w odniesieniu do wskaźnika F1.
Z Tabeli 1 wynika, że podejścia multimodalne (czujniki podczerwieni i urządzenia ubieralne oraz kask, IR+IMU+EEG; oraz urządzenia ubieralne, kask i kamery, IMU+EEG+CAM) uzyskały najwyższe wartości wskaźnika F1 w porównaniu z podejściami unimodalnymi (wyłącznie podczerwień, IR; oraz wyłącznie kamery, CAM). Zauważyliśmy również, że zastosowanie samych urządzeń ubieralnych (IMU) pozwoliło uzyskać wydajność zbliżoną do podejścia multimodalnego. W tym przypadku zdecydowaliśmy się na podejście multimodalne, ponieważ różne źródła informacji mogą niwelować ograniczenia wynikające z innych źródeł. Na przykład uciążliwość kamer może zostać zrekompensowana przez czujniki ubieralne, a brak wszystkich czujników ubieralnych może zostać uzupełniony przez kamery lub czujniki środowiskowe.
W odniesieniu do punktu odniesienia dla modeli opartych na danych, eksperymenty przedstawione w Tabeli 1 wykazały, że RF osiąga najlepsze wyniki w niemal wszystkich badaniach; tymczasem MLP i SVM nie wykazywały spójnej wydajności (np. odchylenie standardowe w tych technikach wykazuje większą zmienność niż w RF). Jeśli chodzi o rozmiary okna, nie odnotowano między nimi żadnej znaczącej poprawy. Należy zauważyć, że eksperymenty te zostały przeprowadzone w celu klasyfikacji upadków oraz aktywności człowieka.
Rozmieszczenie czujników i najlepsza kombinacja multimodalna
Z drugiej strony, dążyliśmy do wyznaczenia najlepszej kombinacji urządzeń multimodalnych do detekcji upadków. W tej analizie ograniczyliśmy źródła informacji do pięciu czujników noszonych oraz dwóch kamer. Urządzenia te zapewniają największy komfort w ramach zastosowanego podejścia. Ponadto rozważono dwie klasy: upadek (każdy rodzaj upadku) lub brak upadku (jakakolwiek inna aktywność). Wszystkie modele uczenia maszynowego oraz rozmiary okien pozostały takie same jak w poprzedniej analizie.
Dla każdego czujnika ubieralnego zbudowano niezależny model klasyfikatora dla każdej długości okna. Model trenowano z wykorzystaniem 10-krotnej walidacji krzyżowej, z podziałem danych na zbiór treningowy w 70% i testowy w 30%. W Tabeli 2 podsumowano wyniki rankingu czujników ubieralnych według wydajności klasyfikatora na podstawie wskaźnika F1-score. Wyniki te zostały uporządkowane w kolejności malejącej. Jak widać w Tabeli 2, najlepszą wydajność uzyskano przy zastosowaniu pojedynczego czujnika umieszczonego w pasie, na szyi lub w prawej kieszeni spodni (obszar zacieniony). Ponadto czujniki ubieralne umieszczone na kostce i lewym nadgarstku wykazały najgorsze wyniki. Tabela 3 przedstawia preferowaną długość okna dla każdego czujnika ubieralnego w celu uzyskania najlepszej wydajności w każdym klasyfikatorze. Z wyników wynika, że czujniki w pasie, na szyi i w prawej kieszeni spodni, przy zastosowaniu klasyfikatora RF oraz rozmiaru okna 3 s z 50% nakładaniem, są najbardziej odpowiednimi czujnikami ubieralnymi do detekcji upadków.
Przeprowadziliśmy podobną analizę dla każdej kamery w systemie. Dla każdego rozmiaru okna zbudowaliśmy niezależny model klasyfikatora. Do szkolenia zastosowano 10-krotną walidację krzyżową z podziałem danych na 70% zestawu treningowego i 30% zestawu testowego. Tabela 4 przedstawia ranking najlepszych punktów widzenia kamer dla każdego klasyfikatora na podstawie miary F1-score. Zaobserwowano, że widok boczny (kamera 1) zapewnił najskuteczniejszą detekcję upadków. Ponadto RF uzyskał lepsze wyniki w porównaniu z pozostałymi klasyfikatorami. Tabela 5 przedstawia preferowaną długość okna dla każdego punktu widzenia kamery. Z wyników wynika, że najlepszą lokalizacją kamery jest widok boczny przy zastosowaniu RF z rozmiarem okna 3 s i 50% nakładaniem.
Na koniec wybrano dwa możliwe miejsca rozmieszczenia czujników noszonych (tj. pas oraz ciasna prawa kieszeń), które miały zostać połączone z kamerą o punkcie widzenia bocznym. Po przeprowadzeniu tej samej procedury treningowej uzyskano wyniki przedstawione w Tabeli 6. Jak pokazano, klasyfikator modelu RF osiągnął najlepszą dokładność oraz wartość wskaźnika F1 w obu multimodalnościach. Ponadto kombinacja pasa i kamery 1 zajęła pierwsze miejsce, osiągając 98,72% dokładności i 95,7% wskaźnika F1.

Rysunek 1: Układ czujników ubieralnych (po lewej) i środowiskowych (po prawej) w bazie danych UP-Fall Detection. Czujniki ubieralne są rozmieszczone na czole, lewym nadgarstku, szyi, w pasie, w prawej kieszeni spodni oraz na lewej kostce. Czujniki środowiskowe obejmują sześć par czujników podczerwieni do wykrywania obecności badanych oraz dwie kamery. Kamery są ustawione w widoku bocznym oraz widoku przednim, oba w odniesieniu do upadku człowieka. Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.

Rysunek 2: Przykład nagrania wideo wyekstrahowanego z bazy danych UP-Fall Detection. W górnej części znajduje się sekwencja obrazów przedstawiająca obiekt upadający na bok. W dolnej części znajduje się sekwencja obrazów reprezentująca wyekstrahowane cechy wizualne. Cechami tymi jest wzglendny ruch pikseli pomiędzy dwoma sąsiadującymi obrazami. Białe piksele reprezentują szybszy ruch, natomiast czarne piksele reprezentują ruch wolniejszy (lub bliski zeru). Sekwencja ta jest ułożona chronologicznie od lewej do prawej. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 3: Wyniki porównawcze przedstawiające najlepszą wartość miary F1 dla każdej modalności w zależności od modelu uczenia maszynowego oraz optymalnej długości okna. Słupki przedstawiają średnie wartości miary F1. Tekst w punktach danych przedstawia średnią oraz odchylenie standardowe w nawiasie. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Modalność | Model | Dokładność (%) | Precyzja (%) | Czułość (%) | Swoistość (%) | Wynik F1 (%) |
| IR | RF (3 sek.) | 67.38 ± 0.65 | 36.45 ± 2.46 | 31.26 ± 0.89 | 96.63 ± 0.07 | 32.16 ± 0.99 |
| SVM (3 sek.) | 65.16 ± 0.90 | 26.77 ± 0.58 | 25.16 ± 0.29 | 96.31 ± 0.09 | 23.89 ± 0.41 |
| MLP (3 sek.) | 65.69 ± 0.89 | 28.19 ± 3.56 | 26.40 ± 0.71 | 96.41 ± 0.08 | 25.13 ± 1.09 |
| kNN (3 s) | 61.79 ± 1.47 | 30.04 ± 1.44 | 27.55 ± 0.97 | 96.05 ± 0.16 | 27.89 ± 1.13 |
| IMU (jednostka pomiaru inercyjnego) | RF (1 s) | 95.76 ± 0.18 | 70.78 ± 1.53 | 66.91 ± 1.28 | 99.59 ± 0.02 | 68.35 ± 1.25 |
| SVM (1 s) | 93.32 ± 0.23 | 66.16 ± 3.33 | 58.82 ± 1.53 | 99.32 ± 0.02 | 60.00 ± 1.34 |
| MLP (1 s) | 95.48 ± 0.25 | 73.04 ± 1.89 | 69.39 ± 1.47 | 99.56 ± 0.02 | 70.31 ± 1.48 |
| kNN (1 s) | 94.90 ± 0.18 | 69.05 ± 1.63 | 64.28 ± 1.57 | 99.50 ± 0.02 | 66.03 ± 1.52 |
| IMU+EEG | RF (1 s) | 95.92 ± 0.29 | 74.14 ± 1.29 | 66.29 ± 1.66 | 99.59 ± 0.03 | 69.03 ± 1.48 |
| SVM (1 s) | 90.77 ± 0.36 | 62.51 ± 3.34 | 52.46 ± 1.19 | 99.03 ± 0.03 | 53.91 ± 1.16 |
| MLP (1 s) | 93.33 ± 0.55 | 74.10 ± 1.61 | 65.32 ± 1.15 | 99.32 ± 0.05 | 68.13 ± 1.16 |
| kNN (1 s) | 92.12 ± 0.31 | 66.86 ± 1.32 | 58.30 ± 1.20 | 98.89 ± 0.05 | 60.56 ± 1.02 |
| IR+IMU+EEG | RF (2 s) | 95.12 ± 0.36 | 74.63 ± 1.65 | 66.71 ± 1.98 | 99.51 ± 0.03 | 69.38 ± 1.72 |
| SVM (1 s) | 90.59 ± 0.27 | 64.75 ± 3.89 | 52.63 ± 1.42 | 99.01 ± 0.02 | 53.94 ± 1.47 |
| MLP (1 s) | 93.26 ± 0.69 | 73.51 ± 1.59 | 66.05 ± 1.11 | 99.31 ± 0.07 | 68.19 ± 1.02 |
| kNN (1 s) | 92.24 ± 0.25 | 67.33 ± 1.94 | 58.11 ± 1.61 | 99.21 ± 0.02 | 60.36 ± 1.71 |
| CAM | RF (3 s) | 32.33 ± 0.90 | 14.45 ± 1.07 | 14.48 ± 0.82 | 92.91 ± 0.09 | 14.38 ± 0.89 |
| SVM (2 s) | 34.40 ± 0.67 | 13.81 ± 0.22 | 14.30 ± 0.31 | 92.97 ± 0.06 | 13.83 ± 0.27 |
| MLP (3 s) | 27.08 ± 2.03 | 8.59 ± 1.69 | 10.59 ± 0.38 | 92.21 ± 0.09 | 7.31 ± 0.82 |
| kNN (3 sek) | 34.03 ± 1.11 | 15.32 ± 0.73 | 15.54 ± 0.57 | 93.09 ± 0.11 | 15.19 ± 0.52 |
| IR+CAM | RF (3 s) | 65.00 ± 0.65 | 33.93 ± 2.81 | 29.02 ± 0.89 | 96.34 ± 0.07 | 29.81 ± 1.16 |
| SVM (3 sek) | 64.07 ± 0.79 | 24.10 ± 0.98 | 24.18 ± 0.17 | 96.17 ± 0.07 | 22.38 ± 0.23 |
| MLP (3 sek) | 65.05 ± 0.66 | 28.25 ± 3.20 | 25.40 ± 0.51 | 96.29 ± 0.06 | 24.39 ± 0.88 |
| kNN (3 sek.) | 60.75 ± 1.29 | 29.91 ± 3.95 | 26.25 ± 0.90 | 95.95 ± 0.11 | 26.54 ± 1.42 |
| IMU+EEG+CAM | RF (1 s) | 95.09 ± 0.23 | 75.52 ± 2.31 | 66.23 ± 1.11 | 99.50 ± 0.02 | 69.36 ± 1.35 |
| SVM (1 s) | 91.16 ± 0.25 | 66.79 ± 2.79 | 53.82 ± 0.70 | 99.07 ± 0.02 | 55.82 ± 0.77 |
| MLP (1 s) | 94.32 ± 0.31 | 76.78 ± 1.59 | 67.29 ± 1.41 | 99.42 ± 0.03 | 70.44 ± 1.25 |
| kNN (1 s) | 92.06 ± 0.24 | 68.82 ± 1.61 | 58.49 ± 1.14 | 99.19 ± 0.02 | 60.51 ± 0.85 |
Tabela 1: Wyniki porównawcze przedstawiające najlepszą wydajność każdej modalności w odniesieniu do modelu uczenia maszynowego oraz najlepszej długości okna (w nawiasie). Wszystkie wartości wydajności reprezentują średnią i odchylenie standardowe.
| # | Typ IMU |
| RF | SVM | MLP | KNN |
| 1 | (98,36) Talia | (83,30) Prawa kieszeń | (57,67) Prawa kieszeń | (73,19) Prawa kieszeń |
| 2 | (95,7) Szyja | (83,2) Talia | (44,93) Szyja | (68,73) Talia |
| 3 | (95,35) Prawa kieszeń | (83,1) Szyja | (39,54) Talia | (65,06) Szyja |
| 4 | (95,06) Kostka | (82,96) Kostka | (39,06) Lewy nadgarstek | (58,26) Kostka |
| 5 | (94,6) Lewy nadgarstek | (82,82) Lewy nadgarstek | (37,56) Kostka | (51,63) Lewy nadgarstek |
Tabela 2: Ranking najlepszych czujników ubieralnych dla każdego klasyfikatora, posortowany według wyniku F1-score (w nawiasach). Zacienione obszary reprezentują trzech najlepszych klasyfikatorów w wykrywaniu upadków.
| typ IMU | Długość okna |
| RF | SVM (maszyna wektorów wspierających) | MLP (wielowarstwowy perceptron) | K-najbliższych sąsiadów (KNN) |
| Lewa kostka | 2 s | 3 sek. | 1 s | 3 sek. |
| Talia | 3 sek. | 1 s | 1 s | 2 s |
| Szyja | 3 sek. | 3 sek. | 2 sek. | 2 s |
| Prawa kieszeń | 3 sek. | 3 sek. | 2 sek. | 2 s |
| Lewy nadgarstek | 2 s | 2 sek. | 2 s | 2 sek. |
Tabela 3: Preferowana długość okna czasowego w czujnikach ubieralnych dla każdego klasyfikatora.
| # | Widok z kamery |
| RF | SVM (Maszyna Wektorów Wspierających) | MLP (wielowarstwowy perceptron) | K-najbliższych sąsiadów (KNN) |
| 1 | (62,27) Widok boczny | (24.25) Widok z boku | (13.78) Widok z przodu | (41.52) Widok boczny |
| 2 | (5,71) Widok z przodu | (0,20) Widok z przodu | (5.51) Widok boczny | (28.13) Widok z przodu |
Tabela 4: Ranking najlepszych punktów widzenia kamery dla każdego klasyfikatora, posortowany według wartości F1-score (w nawiasach). Obszary zaznaczone cieniem reprezentują najlepszy klasyfikator do wykrywania upadków.
| Kamera | Długość okna |
| RF | SVM (maszyna wektorów wspierających) | MLP (wielowarstwowy perceptron) | k-najbliższych sąsiadów |
| Widok z boku | 3 sek. | 3 sek. | 2 sek. | 3 sek. |
| Widok z przodu | 2 s | 2 s | 3 sek. | 2 s |
Tabela 5: Preferowana długość okna czasowego w widokach kamery dla każdego klasyfikatora.
| Tryb multimodalny | Klasyfikator | Dokładność (%) | Precyzja (%) | Czułość (%) | Wskaźnik F1 (%) |
Pas
+
Widok boczny | RF | 98.72 ± 0.35 | 94.01 ± 1.51 | 97.63 ± 1.56 | 95.7 ± 1.15 |
| SVM | 95.59 ± 0.40 | 10 | 70.26 ± 2.71 | 82.51 ± 1.85 |
| MLP | 7.67 ± 1.04 | 3.73 ± 1.69 | 37.1 ± 26.74 | 29.81 ± 12.81 |
| KNN | 91.71 ± 0.61 | 77.90 ± 3.3 | 61.64 ± 3.68 | 68.73 ± 2.58 |
Prawa kieszeń
+
Widok boczny | RF | 98.41 ± 0.49 | 93.64 ± 1.46 | 95.79 ± 2.65 | 94.69 ± 1.67 |
| SVM | 95.79 ± 0.58 | 10 | 71.58 ± 3.91 | 83.38 ± 2.64 |
| MLP | 84.92 ± 2.98 | 5.70 ± 1.36 | 48.29 ± 25.1 | 45.21 ± 14.19 |
| KNN | 91.71 ± 0.58 | 73.63 ± 3.19 | 68.95 ± 2.73 | 71.13 ± 1.69 |
Tabela 6: Wyniki porównawcze połączonego sensora ubieralnego i punktu widzenia kamery przy zastosowaniu okna o długości 3 s. Wszystkie wartości przedstawiają średnią i odchylenie standardowe.