Identyfikacja osobnicza
Możliwość zastosowania techniki identyfikacji śladów w celu klasyfikowania indywidualnych gepardów zależy od dwóch czynników: użycia ujednoliconego protokołu gromadzenia śladów oraz nowego modelu statystycznego opartego na krzyżowo sprawdzonej analizie dyskryminacyjnej parami z analizą skupień Warda. Te czynniki są wspierane przez zintegrowany graficzny interfejs użytkownika do wizualizacji danych (Fig. 1). Wymagane jest minimalne wyposażenie, co czyni tę technikę opłacalną (Lista materiałów). Dane zebrane na podstawie śladów obejmowały liczbę gepardów, liczbę zebranych obrazów śladów, zakres śladów przypadających na jednego geparda, liczbę szlaków, zakres szlaków przypadających na jednego geparda oraz zakres wieku gepardów (Table 1).
Zebrano 781 śladów (kobiety:mężczyźni 395:386) należących do 110 tras, pochodzących od 38 osób, w celu utworzenia zestawu danych uczących. Tabela 1 zawiera podsumowanie zebranych danych. Przy użyciu okna ekstrakcji cech (Ryc. 2) z zestawu 25 punktów orientacyjnych wygenerowano 15 punktów pochodnych na każdym obrazie śladu. Z tych punktów orientacyjnych i pochodnych wygenerowano 136 zmiennych dla każdego śladu, obejmujących odległości, kąty i powierzchnie. Każdy wiersz w bazie danych reprezentował zatem 136 zmiennych wygenerowanych na podstawie pojedynczego śladu. Ślady przetwarzane były według tras. Różna liczba wierszy reprezentowała każdą trasę i była odpowiednio oznaczona.
Dane te zostały powielone w tabeli danych jako jednostka, a następnie określone jako wartość centroidu odniesienia (RCV), która służy do stabilizacji porównania parami śladów niezbędnego do klasyfikacji indywidualnej. Okno analizy parami (Ryc. 3) zostało zaprojektowane w celu ułatwienia weryfikacji danych i/lub testowania danych pochodzących z nieznanych populacji. Ryc. 4 przedstawia wynik porównania parami śladów pochodzących od tej samej osoby (A) oraz dwóch różnych osób (B) na podstawie dostosowanego modelu techniki identyfikacji odcisków stóp. Klasyfikator wbudowany w model opiera się na obecności lub braku nakładania się elips. Należy zauważyć, że analiza jest przeprowadzana dla każdego porównania parami w obecności trzeciej jednostki, tj. wartości centroidu odniesienia (RCV).
Za pomocą niezawodnej, wzajemnie sprawdzanej analizy dyskryminacyjnej połączonej z analizą skupień metodą Warda, opracowano algorytm umożliwiający skuteczną klasyfikację osobników. Algorytm techniki identyfikacji śladów opiera się na trzech regulowalnych elementach: liczbie wykorzystanych pomiarów, wielkości elipsy (zastosowanym przedziale ufności) oraz wartości progowej wyznaczającej granicę odcięcia dla skupień. Każdy z tych elementów jest dostosowywany w oprogramowaniu, aż do osiągnięcia najwyższej dokładności klasyfikacji dla zestawu treningowego zwierząt o znanej tożsamości. Tego samego algorytmu można następnie użyć do identyfikacji nieznanych gepardów. Na przykład rysunek 5a, b i c przedstawia dendrogram prób śladów siedmiu gepardów, pokazujący poprawne przewidywanie, gdy algorytm jest zoptymalizowany (a) oraz w przypadku nieoptymalnego algorytmu (b i c).
Próby z zatrzymaniem zostały przeprowadzone w celu zweryfikowania algorytmu opracowanego na podstawie zestawu treningowego „znanych” osobników. Przeprowadzono je kolejno, zmieniając udział gepardów w zestawach testowych i treningowych. Zamiast przydzielać gepardy do zestawów treningowych i testowych w sposób dowolny, analizy wykonywano sekwencyjnie, zwiększając rozmiar zestawu testowego. Dla każdego zestawu testowego przeprowadzono 10 iteracji, wybierając gepardy losowo w każdej iteracji. Dla każdego zestawu testowego umożliwiano obliczenie wartości średniej. Rysunek 6. przedstawia zmieniający się rozmiar zestawu testowego naniesiony na siebie (czerwony), a na osi y – przewidywaną wartość dla każdej iteracji rozmiaru testowego (zielony) oraz średnią przewidywaną wartość dla każdego rozmiaru testowego (niebieski). Wykres pokazuje, że nawet wtedy, gdy rozmiar zestawu testowego znacznie wzrasta (n=28) w porównaniu z rozmiarem zestawu treningowego (n=10), średnia przewidywana wartość jest zbliżona do wartości oczekiwanej.
Stosując kilka prób kontrolnych, dokładność indywidualnej identyfikacji była konsekwentna >90% zarówno dla przewidywanej liczby osobników, jak i – co równie ważne – klasyfikacji ścieżek, czyli, czy ślady pochodzące od tego samego osobnika (ślady własnego osobnika) oraz ślady pochodzące od różnych osobników (ślady nie-własne) zostały poprawnie sklasyfikowane. Przedstawiono dendrogram skupień obejmujący wszystkie 38 osobniki gepardaRyc. 7). Przeprowadzono 110 prób, uzyskując łącznie 5886 porównań parami. Wśród nich wystąpiło 46 błędnych klasyfikacji, co odpowiada dokładności 99% (Tabela 2).
| Liczba gepardów | Liczba obrazów śladów | Zakres śladów na jednego geparda | Liczba tras | Zakres tras na jednego geparda | Zakres wieku (lata) |
| Samice | 16 | 386 | 12 - 36 | 55 | 2 - 5 | 2,5 - 8,5 |
| Samce | 22 | 395 | 7 - 32 | 54 | 1 - 4 | 1 - 11 |
| Razem | 38 | 781 | 7 - 36 | 109 | 1 - 5 | 1 - 11 |
Tabela 1. Podsumowanie zebranych danych. Liczba gepardów, liczba zebranych obrazów śladów łap, zakres śladów na jednego geparda, liczba szlaków, zakres szlaków na jednego geparda oraz zakres wiekowy gepardów.
| Własne | Niewłasne | Łącznie | Błędne klasyfikacje |
| Własne (N) | 117 | 9 | 126 | 9 |
| Własne (%) | 93 | 7 | 100 | 7 |
| Niewłasne (N) | 37 | 5,723 | 5,760 | 37 |
| Niewłasne (%) | 1 | 99 | 100 | 1 |
| Łącznie (N) | - | - | 5,886 | 46 |
| Łącznie (%) | - | - | 100 | 1 |
Tabela 2. Wynik w oprogramowaniu metody identyfikacji śladów pokazuje klasyfikację ścieżek na podstawie porównań parami. „Samo” oznacza ścieżki pochodzące od tej samej osoby, a „nie-samo” – ścieżki od różnych osób. Każda ścieżka była porównywana z każdą inną za pomocą dostosowanego, odpornego modelu analizy dyskryminacyjnej z walidacją krzyżową. 110 ścieżek dało w efekcie 5886 porównań parami, a ogólna dokładność klasyfikacji wyniosła 99%.

Rycina 1. Okno głównego menu startowego w technice identyfikacji odcisków śladów. Jest to dodatek do oprogramowania wizualizacji danych, zaprojektowany do klasyfikowania odcisków śladów według jednostki, płci i klasy wiekowej na podstawie pomiarów morfometrycznych. Interfejs graficzny użytkownika umożliwia płynną nawigację między różnymi opcjami. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 2. Okno ekstrakcji cech. Możliwości obejmują przeciąganie i upuszczanie obrazów, automatyczne dopasowanie rozmiaru do okna, obracanie obrazów w celu standaryzacji, uwzględnianie głębokości podłoża, itp. Wstępnie przypisane punkty orientacyjne są ręcznie umieszczane i generują serię pochodnych punktów wyznaczonych przez skrypt, umożliwiając ekstrakcję metryk w postaci odległości, kątów i powierzchni. Wynik ma postać wiersza danych zawierającego współrzędne x,y oraz metryki. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Ryc. 3. Okno analizy danych parami w technice identyfikacji śladów. Po utworzeniu bazy danych pomiarów, okno analizy parami służy do weryfikacji danych i/lub testowania danych pochodzących z nieznanych populacji. Analiza opiera się na niestandardowym modelu obejmującym stałą wartość centroidu odniesienia (RCV), który porównuje pary szlaków sekwencyjnie16,17. Końcowym wynikiem jest dendrogram skupień, który pozwala przewidzieć liczbę osobników oraz relacje między szlakami. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 4. Porównania parami. Na rycinie przedstawiono wyniki porównania parami torów pochodzących od tego samego osobnika (A) oraz od dwóch różnych osobników (B), na podstawie spersonalizowanego modelu w oprogramowaniu do wizualizacji danych. Klasyfikator wbudowany w model opiera się na obecności lub braku nakładania się elips. Należy zauważyć, że analizę przeprowadza się dla każdego porównania parami w obecności trzeciego podmiotu, tj. wartości centroidu odniesienia (RCV). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5. Dendrogram próbki tras siedmiu gepardów, pokazujący poprawne przewidywanie, gdy algorytm jest zoptymalizowany (a) oraz w przypadku algorytmu nieoptymalnego (b i c). d przedstawia wynik dla 18 zmiennych, przy przesunięciu suwaka w jedną stronę, co pokazuje, że szansa dla dziesięciu gepardów jest mniejsza niż 50%. Algorytm opiera się na trzech regulowalnych parametrach: liczbie wykorzystanych pomiarów, wielkości elipsy (używanym przedziale ufności) oraz wartości progowej, która określa wartość odcięcia dla skupień. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 6. Badanie z odłożeniem przeprowadzone sekwencyjnie poprzez zmianę proporcji gepardów w zestawach testowych i uczących. Zamiast arbitralnie przydzielać gepardy do zestawów uczących i testowych, analizę przeprowadzono sekwencyjnie zwiększając rozmiar zestawu testowego. Dla każdego zestawu testowego wykonano dziesięć iteracji, przy czym gepardy były wybierane losowo w każdej iteracji. Dla każdego zestawu testowego umożliwiono obliczenie wartości średniej. Wykres przedstawia zmieniający się rozmiar zestawu testowego naniesiony na siebie (czerwony), a na osi y wartość przewidywaną dla każdej iteracji rozmiaru zestawu testowego (zielony) oraz średnią wartość przewidywaną dla każdego rozmiaru zestawu testowego (niebieski). Wykres pokazuje, że nawet wtedy, gdy rozmiar zestawu testowego jest znacznie zwiększony (n=28) w porównaniu do rozmiaru zestawu uczącego (n=10), średnia wartość przewidywana jest zbliżona do wartości oczekiwanej. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 7. Dendrogram przedstawiający przewidywane wyniki, gdy wszystkie 110 tras z 38 gepardów są uwzględnione w analizie. Zwróć uwagę na spójność tras tworzących klastry. Ciekawe, że wiele błędnych klasyfikacji wystąpiło między rodzeństwem, np. gepardami Letotse/Duma oraz Vincent/Bonsai. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.