$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Zbiór danych zawiera łącznie 1 190 przekrojów tomografii komputerowej z 110 przypadków, podzielonych na trzy kategorie: normalne (55 przypadków), łagodne (15 przypadków) i złośliwe (40 przypadków). Każda przypadek składa się z wielu przekrojów tomografii (od około 80 do 200 nacięć na przypadek), oferując różnorodne aksie, widoki na okolicę klatki piersiowej. Obrazy CT zostały uzyskane w DICOM za pomocą skanera SOMATOM ze standardowymi parametrami obrazowania: napięcie rurki = 120 kV, grubość przekroju = 1 mm, szerokość okna = 350–1 200 jednostek Hounsfielda (HU) oraz wartości środka okna = 50–600 HU, podczas pełnego wstrzymania oddechu.
Wszystkie obrazy zostały całkowicie zdeanonimowane przed analizą, aby usunąć wszelkie dane osobowe (PII). Zbiór danych został etycznie zatwierdzony przez instytucjonalne komisje rewizyjne uczestniczących ośrodków medycznych, a pisemna zgoda została zrzeczona przez komisję nadzorczą. Sprawy obejmowały osoby o różnym pochodzeniu, takie jak pracownicy rządowi, rolnicy oraz mieszkańcy kilku prowincji Iraku (w tym Bagdadu, Wasit, Diyala, Salahuddin i Babilonu), o zróżnicowanym płci, wieku, poziomie wykształcenia i statusie zamieszkania.
Ponieważ zbiór danych jest publicznie dostępny i bez identyfikacji, nie było potrzeby dodatkowej zgody etycznej do jego wykorzystania w tym badaniu. Zbiór danych spełnia warunki określone przez jego pierwotnych współtwórców oraz platformę hostingową.
Metodologia tych badań wykorzystuje wieloetapowy proces zaprojektowany do stworzenia modelu predykcyjnego z pomocą zestawu danych IQ-OTH/NCCD dotyczącego raka płuc21. Ta metodologia tworzy solidne podstawy dla urządzeń medycznych skoncentrowanych na konsumencie, gdzie wymagane są mniejsze opóźnienia, co może zapewnić lepszą dokładność. Rysunek 2 przedstawia mechanizm działania proponowanego modelu.

Rysunek 2: Diagram workflow proponowanego modelu przedstawiający etapy wstępnego przetwarzania, augmentacji, klasyfikacji Vision Transformer oraz ewaluacji. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
1. Opis zbioru danych
Zbiór danych IQ-OTH/NCCD dotyczący raka płuc został zebrany w Szpitalu Klinicznym Iracko-Onkologicznym/Narodowym Centrum Chorób Nowotworowych jesienią 2019 roku. Część treningowa użytego zbioru danych zawierała 1 097 obrazów, których opis przedstawiono w Tabeli 2, ilustrującej liczbę przykładów różnych klas.
| Klasa | Liczba obrazów |
| Normalne | 416 |
| Łagodne | 120 |
| Złośliwy | 561 |
Tabela 2: Przykładowe przypadki prawidłowych, łagodnych i złośliwych obrazów tomografii komputerowej z tego zbioru danych.
Zestaw danych IQ-OTH/NCCD dotyczący raka płuc został wybrany ze względu na kompleksową reprezentację prawidłowych, łagodnych i złośliwych tomografii komputerowej. Chociaż dostępne są inne zbiory danych, takie jak LIDC-IDRI i NSCLC-Radiomics, skupiają się one głównie na wykrywaniu guzków lub nie mają wystarczającej liczby próbek przypadków łagodnych. Zbiór danych IQ-OTH/NCCD jest szczególnie odpowiedni do naszych badań, ponieważ pozwala Transformatorowi Wzroku poznać cechy dyskryminacyjne we wszystkich trzech klasach, umożliwiając solidną klasyfikację subtelnych wzorców w obrazach tomografii płuc.
2. Wstępne przetwarzanie danych
Wstępne przetwarzanie jest ważnym krokiem w poprawie wydajności modelu poprzez spójność i odpowiednie dostosowania danych przetwarzanych przez sieć neuronową. Aby zapewnić spójność rozmiaru wejścia dla sieci neuronowej, przeskalowaliśmy wszystkie obrazy do tego samego wymiaru 256 x 256 pikseli i znormalizowaliśmy dane do wartości pikseli między 0 a 1, mając nadzieję na poprawę trenowania i zbieżności modeli. Tabela 3 zawiera wartości technik augmentacji.
| Technika | Wartość |
| Normalizacja | - |
| Zmiana rozmiaru | image_size x image_size |
| RandomRotation | czynnik=0,02 |
| RandomZoom | height_factor=0,2, width_factor=0,2 |
Tabela 3: Stosowane techniki augmentacji z zakresami parametrów.
Aby zwiększyć odporność modelu na nadmierne dopasowanie i poprawić jego zdolność do uogólniania, stosuje się techniki uzupełniania danych, takie jak losowe obroty i zoomy, symulujące różne kąty i rozmiary objawów raka płuc, tak jak mogą pojawić się u różnych pacjentów. W badaniu zastosowano losowe rotacje o kątach jednolicie pobieranych z 0,02 radiana oraz losowe transformacje powiększenia o zmiennych współczynnikach wysokości i szerokości. Obrazy po augmentacji przedstawiono na Rysunku 3.

Rysunek 3: Obrazy CT po augmentacji pokazujące rotację, powiększenie i normalizację w celu poprawy uogólnienia modelu. Proszę kliknąć tutaj, aby zobaczyć większą wersję tego rysunku.
Te techniki wstępnego przetwarzania są wymagane, a we wszystkich klasach zastosowano augmentację, aby zapewnić odporność modelu.
3. Architektura modelu
Dostosowując nowatorski framework Vision Transformer (ViT) do efektywnej obsługi skomplikowanych danych obrazowych, architektura modelu ma na celu podzielić tomografię komputerową na trzy kategorie: normalne, łagodne i złośliwe. Tą metodą przetwarzane są obrazy wejściowe o wymiarach 256 x 256 pikseli. Aby zagwarantować spójność i promować bardziej efektywne uczenie się modeli, każdy obraz przechodzi kilka operacji wstępnego przetwarzania, takich jak zmiana rozmiaru i normalizacja. Aby zwiększyć odporność modelu na zmiany skali i orientacji obrazu, projekt rozpoczyna się od warstwy augmentacji danych, która wykorzystuje metody takie jak normalizacja, zmiana rozmiaru, losowe obroty 0,02 radiana oraz losowe powiększenia do 20%. Po augmentacji model usuwa 16 x 16 pikseli z każdego zdjęcia, dzięki czemu każdy obraz ma 256 fragmentów.
Algorytm 1 definiuje przepływ pracy, proponowany model i sposób jego budowy, a także precyzyjne dostrojenie względem modelu.
Algorytm 1: Klasyfikacja raka płuc z użyciem transformatorów wzroku
Wejście wejściowe: Zestaw obrazów CT I z zestawu danych IQ-OTH/NCCD
Wyniki: Wyniki klasyfikacji C do kategorii: Normalne, Łagodne, Złośliwe
Wstępne przetwarzanie:
for each image i in I do
i <- Resize(i, 256 x 256) // Normalize and resize images
i <- Normalize(i)
i <- DataAugmentation(i) // Apply random rotations and zooms
end for
Konfiguracja modelu:
Initialize Vision Transformer (ViT) Model
Set number of patches P = 16 x 16
Set number of heads H = 6 in multi-head attention
Szkolenia:
for epoch = 1 to MaxEpochs do
for each batch b in I do
Patches <- ExtractPatches(b, P)
EncodedPatches <- PatchEncoder(Patches)
AttentionWeights <- MultiHeadAttention(EncodedPatches, H)
ClassLogits <- TransformerEncoder(AttentionWeights)
Loss <- ComputeLoss(ClassLogits, TrueLabels)
UpdateModelWeights(Loss)
end for
if EarlyStoppingCriteriaMet (Val_Loss No Improvement Patience = 5 Epochs) then
break
end if
end for
Weryfikacja:
Split data into TrainingSet (80%) and ValidationSet (20%)
ComputeValidationMetrics(ValidationSet)
Ocena:
C <- Classify(I)
ComputePerformanceMetrics(C)
Return C
Aby zachować relacje przestrzenne między łatkami, łaty te są spłaszczane (równanie 1) do 768-wymiarowych wektorów (ponieważ każdy łat ma 16 x 16 x 3 = 768), a następnie przesyłane przez warstwę kodującą łatki, która rzutuje każdy wektor do przestrzeni 64-wymiarowej poprzez integrację osadzeń pozycyjnych. Rdzeń architektury składa się z ośmiu warstw transformatorów, z których każda posiada mechanizm skupienia uwagi z wieloma głowami i sześcioma głowami, pozwalający modelowi jednocześnie skupić się na różnych częściach obrazu i uchwycić szeroki zakres cech. Jest reprezentowany za pomocą równań 2, 3 i 4.

Gdzie μ to średnia, aσ 2 to wariancja wejścia x, a ε jest małą stałą, aby zapobiec dzieleniu przez zero.

Gdzie Q to macierz zapytania, K macierz klucza, V macierz wartości, a dk to wymiar wektorów klucza.


Gdzie WiQ, WiK i W iV to macierze wag nauczonych dla zapytań, kluczy i wartości, a WO to macierz wag wyjściowych.
Diagram blokowy proponowanego modelu został przedstawiony na Rysunku 4.

Rysunek 4: Schemat blokowy modelu Vision Transformer z głowicą MLP, przedstawiający główne warstwy przetwarzania i architekturę. Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.
Każda warstwa transformatora zawiera percepcję wielowarstwową (MLP) z ukrytymi jednostkami najpierw skalowanymi do 128, a następnie z powrotem do 64, skutecznie rozszerzając i kompresując przepływ informacji, aby uchwycić złożone zależności.
Rezygnacja jest strategicznie stosowana w tempie 0,1 w mechanizmach uwagi i MLP, aby zapobiec nadmiernemu dopasowaniu. Wyjście z enkodera transformatora jest przetwarzane przez głowicę MLP składającą się z dwóch gęstych warstw o odpowiednio 2 048 i 1 024 jednostkach, wykorzystując do aktywacji jednostki liniowe błędu Gaussa (GELU), które wykazują dobre wyniki w aplikacjach głębokiego uczenia. Osiąga się to za pomocą równania 5.

Gdy W1 i W2 to macierze wagowe, b1 i b2 to błędy, a GELU to używana funkcja aktywacji.
W warstwach transformatorów zastosowano wypadek 0,1, aby zapobiec nadmiernemu dopasowaniu bez utraty ważnych informacji o cechach. Funkcja aktywacji GELU została wykorzystana ze względu na jej gładkie, nieliniowe właściwości, które sprzyjają przepływowi gradientowemu i zbieżności w modelach opartych na transformatorach. Przerwanie jest regularizowane za pomocą równania 6.
Gdzie p to wskaźnik dropoutu (np. 0,1 lub 0,5), a warstwa dropout losowo ustawia ułamek p jednostek wejściowych na zero podczas treningu.
Wskaźnik wypadania 0,5 w tych warstwach dodatkowo pomaga w ograniczaniu nadmiernego dopasowania. Ostatnią warstwą modelu jest warstwa logitów (równanie 7), która generuje trzy klasy logitów odpowiadające kategoriom normalnym, łagodnym i złośliwym, wykorzystując rzadką kategoryczną funkcję strat entropii krzyżowej (równanie 8), odpowiednią dla tego wieloklasowego ustawienia klasyfikacji.


Gdzie zi jest wektorem logitów dla klasy i, SoftMax(zi)) reprezentuje przewidywane prawdopodobieństwa, a yi to prawdziwa etykieta klasy.
Model jest kompilowany za pomocą optymalizatora AdamW (równania 9, 10, 11, 12 i 13), rozszerzenia optymalizatora Adama, które skuteczniej obsługuje spadek wagi, ustawiony z częstotliwością uczenia 0,001 i spadkiem wagi 0,0001, optymalizując zarówno szybkość uczenia, jak i stabilność modelu.





Gdzie mt i vt to pierwszy i drugi moment gradientów,
a
są momentami skorygowanymi przez polaryzację, η to szybkość uczenia się, a ε jest małą stałą zapobiegającą dzieleniu przez zero.
Model po wytrenowaniu wykorzystuje rozmiar partii 32, aby wykorzystać akcelerację GPU dla szybszych obliczeń i jest skonfigurowany do trenowania przez 100 epok.
Jednak trening ma mechanizm wczesnego zatrzymania straty walidacyjnej, aby ograniczyć trening, gdy model przestaje się poprawiać, oszczędzając zasoby obliczeniowe i unikając przetrenowania w 5 kolejnych epokach. Wydajność modelu śledzi się za pomocą wskaźników takich jak rzadka dokładność kategoryczna oraz dokładność top 2, które informują o zdolności modelu do klasyfikacji w najbardziej prawdopodobnych przewidywanych kategoriach. Callbacki do treningu modelu obejmują checkpointy, które zapisują najlepiej działający model według dokładności walidacji, aby zapewnić zachowanie najbardziej skutecznej wersji modelu po zakończeniu procesu treningowego. Ta solidna i dobrze zaplanowana architektura wykorzystuje możliwości transformatorów do przetwarzania danych obrazowania medycznego, oferując zaawansowane podejście do wykorzystania współczesnych metod AI w kluczowych zastosowaniach w diagnostyce opieki zdrowotnej.
4. Szkolenie i walidacja
Model był trenowany w środowisku Kaggle z użyciem karty NVIDIA P100, a podczas procesu treningowego zastosowano mini-batch gradient descent o rozmiarze partii 32. Optymalizator używany podczas treningu to AdamW o tempie uczenia 0,001 i spadku wagi 0,0001; To była dobra równowaga między szybką konwergencją a pewną regularizacją. Model był trenowany przez 100 epok, jednak stosowano wczesne zatrzymanie się na utracie walidacji z cierpliwością 5 epok. Mówiąc prosto, jeśli trening nie poprawiałby utraty walidacyjnej przez 5 ciągłych epok, trening zatrzymałby się z powodu nadmiernego dopasowania i efektywności obliczeniowej. W większości przypadków model ten przywracał wagi z epoki przy najmniejszych stratach walidacyjnych, co wskazywało, że działał optymalnie i nie musiał uruchamiać całych 100 epok. Łącznie około 32 minuty na wytrenowanie modelu.
Podczas szkolenia metryki obejmowały rzadką dokładność kategoryczną oraz dokładność z dwóch najwyższych ocen, monitorowane zarówno na zestawach treningowych, jak i walidacyjnych. Te metryki dają wgląd w to, jak często model przewiduje właściwą klasę oraz czy właściwa klasa znajduje się w dwóch najważniejszych prognozach, co jest istotne w zastosowaniach medycznych, ponieważ błędne klasyfikacje o wysokim stopniu pewności mogą mieć istotne konsekwencje. Krzywe uczenia się strat i dokładności przedstawiono na Rysunku 5.
.
Rysunek 5: Krzywe dokładności i strat treningu i walidacji na przestrzeni 50 epok pokazujące stabilne uczenie się i minimalne nadmierne dopasowanie. Proszę kliknąć tutaj, aby zobaczyć większą wersję tego rysunku.
W obecnym badaniu nie przeprowadzono żadnej walidacji krzyżowej. Takie rozwiązanie pozwalało na efektywne eksperymentowanie z architekturą modelu, w tym dostosowania warstw transformatorów i rozmiarów głowic MLP, jednocześnie zapewniając dobre uogólnienie modelu na podstawie niewidocznych danych walidacyjnych.
5. Analiza statystyczna
Wydajność modelu Vision Transformer została statystycznie przeanalizowana na podstawie zbioru danych IQ-OTH/NCCD dotyczących raka płuc. Precyzja, przypomnienie, wynik F1 i dokładność zostały obliczone za pomocą równań odpowiednio 14, 15, 16, 17. Te cztery są uznawane za konwencjonalne miary stosowane w zadaniach klasyfikacji, ponieważ mogą ujawniać informacje o wydajności modelu przy klasyfikowaniu i klasyfikowaniu na zasadzie każdej klasy.




Przypomnienie to miara zdolności modelu do identyfikacji wszystkich istotnych przypadków w klasie, natomiast precyzja to proporcja pozytywnych identyfikacji, które faktycznie były poprawne. Wynik F1 łączy przywołanie i precyzję, gdy oba są ważne.
Do tego zadania użyto miar wydajności Współczynnika Korelacji Matthewsa (MCC) 18 oraz Kappy 19 Cohena.


Gdzie P0 to zaobserwowane porozumienie, a Pe to oczekiwane porozumienie.
MCC jest wszechstronnym pomiarem wydajności klasyfikacji, uwzględniającym zarówno prawdziwie pozytywne, jak i negatywne wyniki oraz fałszywe pozytywy oraz fałszywe negatywy. Jej wartość może mieścić się w przedziale od -1 do 1, gdzie 1 to perfekcyjne przewidywanie, 0 to losowa, a -1 to całkowite niezgodność między przewidywaną a prawdziwą etykietą. MCC był cenny do porównań między różnymi wynikami modeli, gdy był stosowany do niezrównoważonych zbiorów danych oraz zapewniał miarę zrównoważoną niezależnie od różnic wielkości klas.
W ramach dodatkowej analizy statystycznej obliczono wynik F2, z priorytetem przypominania, zgodnie z równaniem 20.

Wydajność modelu była również kwantyfikowana za pomocą średniego błędu kwadratowego (MSE), rdzenia średniego błędu kwadratowego (RMSE) oraz średniego błędu bezwzględnego (MAE) – które zazwyczaj są pomiarami do zadań regresji, ale tutaj zmodyfikowanymi dla zadania klasyfikacyjnego, aby określić, jak duży był średnio błąd bez uwzględnienia kierunku.
Aby ustalić, czy praktyczne byłoby integrowanie ViT z urządzeniami medycznymi skierowanymi do konsumentów, przeprowadziliśmy szeroko zakrojone wyniki oceny wydajności, skupiając się wyłącznie na efektywności czasowej modelu. Stworzyliśmy funkcje raportujące czas potrzebny na przewidzenie pojedynczego lub po prostu wielu obrazów, ponieważ jest to szczególnie istotne w zastosowaniach w czasie rzeczywistym. Dla każdego obrazu, przed rozpoczęciem przewidywania, dane są najpierw wstępnie przetwarzane, aby miały kształt obrazu, jaki model chce. Czas rozpoczęcia prognozy i czas jej zakończenia zapisywaliśmy, aby uzyskać wyniki dotyczące czasu i opóźnień. Czas i średnie opóźnienie zostały obliczone odpowiednio za pomocą równania 21 i 22.


Gdzie:
- N to liczba obrazów (próbek).
- tend to czas zapisywany po przewidywaniu i-tego obrazu.
- Tstart to czas zapisany przed przewidywaniem i-tego obrazu.
Przeprowadzono dalsze eksperymenty, aby ocenić odporność modelu Vision Transformer, który zaproponowaliśmy, systematycznie wprowadzając dodatkowy szum i rozmycie do obrazów. Do każdego piksela dodawano szum gaussowski o współczynniku szumu 0,4, a wartości pikseli przycinano w zakresie [0,1]. Czynnik rozmycia został zmniejszony przez rozmycie Gaussa o rozmiarze jądra 45× 45. Eksperymenty te mają na celu kompleksową ocenę modelu przy symulowanych pogorszeniu jakości obrazu percepcyjnego.