$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Der Datensatz umfasst insgesamt 1.190 CT-Scan-Schnitte aus 110 Fällen, die in drei Klassen eingeteilt sind: normal (55 Fälle), gutartig (15 Fälle) und bösartig (40 Fälle). Jeder Fall umfasst mehrere CT-Schnitte (mit etwa 80 bis 200 Schnitten pro Fall) und bieten vielfältige axiale Ansichten des Thoraxbereichs. Die CT-Bilder wurden in DICOM mit dem SOMATOM-Scanner mit Standardbildparametern aufgenommen: Röhrenspannung = 120 kV, Schnittdicke = 1 mm, Fensterbreite = 350–1.200 Hounsfield-Einheiten (HU) und Fenstermittelwerte = 50–600 HU, während eines vollständigen Atmungsanhaltens.
Alle Bilder wurden vor der Analyse vollständig deidentifiziert, um personenbezogene Daten (PII) zu entfernen. Der Datensatz wurde ethisch von den institutionellen Überprüfungsgremien der teilnehmenden medizinischen Zentren genehmigt, wobei die schriftliche Zustimmung vom Aufsichtsausschuss aufgehoben wurde. Die Fälle umfassten Personen aus unterschiedlichen Hintergründen, wie Regierungsangestellte, Landwirte und Einwohner mehrerer irakischer Provinzen (darunter Bagdad, Wasit, Diyala, Salahuddin und Babylon) mit unterschiedlichen Geschlechter, Alter, Bildungsniveau und Lebenshaltung.
Da der Datensatz öffentlich zugänglich und deidentifiziert ist, war für seine Verwendung in dieser Studie keine zusätzliche ethische Genehmigung erforderlich. Der Datensatz hält sich an die von den ursprünglichen Mitwirkenden und der Hosting-Plattform festgelegten Geschäftsbedingungen.
Die Methodik dieser Forschung verwendet einen mehrstufigen Prozess, der darauf ausgelegt ist, mit Hilfe des IQ-OTH/NCCD-Lungenkrebs-Datensatzes21 ein Vorhersagemodell zu erstellen. Diese Methodik schafft eine starke Grundlage für verbraucherzentrierte Gesundheitsgeräte, bei denen weniger Latenz erforderlich ist, was zu besserer Genauigkeit führen kann. Abbildung 2 zeigt den Arbeitsmechanismus des vorgeschlagenen Modells.

Abbildung 2: Workflow-Diagramm des vorgeschlagenen Modells, das Vorverarbeitung, Augmentation, Vision Transformer-Klassifikation und Bewertungsschritte darstellt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
1. Datensatzbeschreibung
Der IQ-OTH/NCCD-Lungenkrebs-Datensatz wurde im Herbst 2019 am Iraq-Oncology Teaching Hospital/National Center for Cancer Diseases gesammelt. Der Trainingsteil des verwendeten Datensatzes enthielt 1.097 Bilder, deren Beschreibung in Tabelle 2 angegeben ist, die die Anzahl der Beispielinstanzen verschiedener Klassen veranschaulicht.
| Klasse | Anzahl der Bilder |
| Normal | 416 |
| Harmlos | 120 |
| Bösartig | 561 |
Tabelle 2: Beispiele von normalen, gutartigen und bösartigen CT-Bildern aus dem Datensatz.
Der IQ-OTH/NCCD-Lungenkrebs-Datensatz wurde wegen seiner umfassenden Darstellung normaler, gutartiger und bösartiger CT-Scans ausgewählt. Obwohl andere Datensätze wie LIDC-IDRI und NSCLC-Radiomics verfügbar sind, konzentrieren sie sich entweder hauptsächlich auf die Entdeckung von Knoten oder es fehlen genügend Stichproben gutartiger Fälle. Der IQ-OTH/NCCD-Datensatz ist besonders geeignet für unsere Studie, da er es dem Vision Transformer ermöglicht, diskriminierende Merkmale in allen drei Klassen zu erlernen und so eine robuste Klassifikation subtiler Muster in Lungen-CT-Bildern zu ermöglichen.
2. Datenvorverarbeitung
Die Vorverarbeitung ist ein wichtiger Schritt zur Verbesserung der Modellleistung durch Konsistenz und angemessene Anpassungen der Daten, die im neuronalen Netzwerk verarbeitet werden. Um eine Konsistenz in der Eingabegröße für das neuronale Netzwerk sicherzustellen, haben wir alle Bilder auf die gleiche Größe von 256 x 256 Pixel skaliert und die Daten auf einen Pixelwert zwischen 0 und 1 normalisiert, in der Hoffnung, das Modelltraining und die Konvergenz zu verbessern. Tabelle 3 enthält die Werte der Augmentationstechnik.
| Technik | Wert |
| Normalisierung | - |
| Größenänderung | image_size x image_size |
| RandomRotation | faktor=0,02 |
| RandomZoom | height_factor=0,2, width_factor=0,2 |
Tabelle 3: Angewandte Augmentationstechniken mit Parameterbereichen.
Um die Robustheit des Modells gegen Overfitting zu erhöhen und seine Verallgemeinerungsfähigkeit zu verbessern, werden Datenerweiterungstechniken wie zufällige Rotationen und Zooms angewendet, die verschiedene Winkel und Größen von Lungenkrebsmanifestationen simulieren, wie sie bei verschiedenen Patienten auftreten können. In dieser Studie wurden zufällige Rotationen mit gleichmäßig abgeleiteten Winkeln von 0,02 Radianten und zufällige Zoom-Transformationen mit unterschiedlichen Höhen- und Breitenfaktoren angewendet. Die nachfolgenden Bilder sind in Abbildung 3 dargestellt.

Abbildung 3: CT-Bilder nach der Augmentation, die Rotation, Zoom und Normalisierung zur Verbesserung der Modellverallgemeinerung demonstrieren. Bitte klicken Sie hier, um eine größere Version dieser Abbildung zu sehen.
Diese Vorverarbeitungstechniken sind erforderlich, und Augmentation wurde über alle Klassen hinweg eingesetzt, um die Robustheit des Modells sicherzustellen.
3. Modellarchitektur
Durch die Anpassung des neuartigen Vision Transformer (ViT)-Frameworks, um komplexe Bilddaten effektiv zu verarbeiten, soll die Modellarchitektur CT-Scans in drei Kategorien einteilen: normal, gutartig und bösartig. Mit dieser Methode werden 256 x 256 Pixel große Eingabebilder verarbeitet. Um Konsistenz zu gewährleisten und ein effizienteres Modelllernen zu fördern, wird jedes Bild mehreren Vorverarbeitungsoperationen unterzogen, wie z. B. Größenänderung und Normalisierung. Um die Widerstandsfähigkeit des Modells gegenüber Änderungen der Bildgröße und -ausrichtung zu verbessern, beginnt das Design mit einer Datenerweiterungsschicht, die Methoden wie Normalisierung, Größenänderung, zufällige Rotationen von 0,02 Radianten und zufällige Zooms von bis zu 20 % verwendet. Nach der Augmentation entnimmt das Modell 16 x 16 Pixel Patches aus jedem Bild, sodass jedes Bild 256 Patches hat.
Algorithmus 1 definiert den Workflow, das vorgeschlagene Modell und dessen Aufbau sowie die Feinabstimmung, die in Bezug auf das Modell durchgeführt wird.
Algorithmus 1: Lungenkrebsklassifikation mittels Vision Transformers
Eingabe: Satz von CT-Bildern I aus dem IQ-OTH/NCCD-Datensatz
Ausgabe: Klassifikationsergebnisse C in Kategorien: Normal, Gutartig, Bösartig
Vorverarbeitung:
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
Modellaufbau:
Initialize Vision Transformer (ViT) Model
Set number of patches P = 16 x 16
Set number of heads H = 6 in multi-head attention
Ausbildung:
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
Validierung:
Split data into TrainingSet (80%) and ValidationSet (20%)
ComputeValidationMetrics(ValidationSet)
Auswertung:
C <- Classify(I)
ComputePerformanceMetrics(C)
Return C
Um die räumlichen Beziehungen zwischen den Patches zu erhalten, werden diese Patches (Gleichung 1) in 768-dimensionale Vektoren abgeflacht (da jeder Patch 16 x 16 x 3 = 768 hat) und dann durch eine Patch-Codierungsschicht geschickt, die jeden Vektor durch Integration von Positionseinbettungen in einen 64-dimensionalen Raum projiziert. Der Kern der Architektur besteht aus acht Transformatorschichten, von denen jede einen Multi-Head-Attention-Mechanismus mit sechs Köpfen besitzt, sodass das Modell gleichzeitig auf verschiedene Bildteile fokussieren und eine breite Palette von Merkmalen erfassen kann. Sie wird durch die Gleichungen 2, 3 und 4 dargestellt.

Dabei ist μ der Mittelwert und σ2 die Varianz der Eingabe x und ε eine kleine Konstante, um eine Division durch Null zu verhindern.

Dabei ist Q die Abfragematrix, K die Schlüsselmatrix, V die Wertmatrix und dk die Dimension der Schlüsselvektoren.


Dabei sind WiQ, WiK und WiV die gelernten Gewichtungsmatrizen für Abfragen, Schlüssel und Werte bzw. WO die Ausgabegewichtmatrix.
Das Blockdiagramm des vorgeschlagenen Modells ist in Abbildung 4 dargestellt.

Abbildung 4: Blockdiagramm des Vision Transformer-Modells mit MLP-Kopf, das die wichtigsten Verarbeitungsschichten und die Architektur darstellt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Jede Transformatorschicht enthält eine Mehrschicht-Wahrnehmung (MLP) mit versteckten Einheiten, die zunächst auf 128 und dann wieder auf 64 skaliert werden, wodurch der Informationsfluss effektiv erweitert und komprimiert wird, um komplexe Abhängigkeiten zu erfassen.
Dropout wird strategisch mit einer Rate von 0,1 in den Aufmerksamkeitsmechanismen und MLPs angewendet, um Überanpassungen zu verhindern. Der Ausgang des Transformatorencoders wird durch einen MLP-Kopf verarbeitet, der aus zwei dichten Schichten mit jeweils 2.048 bzw. 1.024 Einheiten besteht, wobei Gaussian Error Linear Units (GELU) zur Aktivierung verwendet werden, was sich als gut in Deep-Learning-Anwendungen erwiesen hat. Dies wird mit Gleichung 5 erreicht.

Dabei sind W1 und W2 Gewichtsmatrizen, b1 und b2 Verspannungen und GELU die verwendete Aktivierungsfunktion.
Ein Dropout von 0,1 wurde in den Transformatorschichten verwendet, um Überfitting zu verhindern, ohne wichtige Merkmalsinformationen zu verlieren. Die GELU-Aktivierungsfunktion wurde wegen ihrer glatten, nichtlinearen Eigenschaften genutzt, die in transformatorbasierten Modellen den Gradientenfluss und die Konvergenz fördern. Der Ausstieg wird mit Gleichung 6 regularisiert.
Dabei ist p die Dropout-Rate (z. B. 0,1 oder 0,5), und die Dropout-Schicht setzt während des Trainings zufällig einen Bruchteil p der Eingangseinheiten auf null.
Eine Aussetzrate von 0,5 in diesen Schichten trägt zusätzlich zur Abmilderung von Overfitting bei. Die letzte Schicht des Modells ist eine Logit-Schicht (Gleichung 7), die drei Klassen-Logits ausgibt, die den Kategorien normal, gutartig und bösartig entsprechen, wobei eine spärliche kategoriale Kreuzentropieverlustfunktion (Gleichung 8) verwendet wird, die für diese Mehrklassenklassifikation geeignet ist.


Dabei ist zi der Logitsvektor für die Klasse i, SoftMax(zi)) die vorhergesagten Wahrscheinlichkeiten und yi das wahre Klassenlabel.
Das Modell wird mit dem AdamW-Optimierer (Gleichungen 9, 10, 11, 12 und 13) kompiliert, einer Erweiterung des Adam-Optimierers, der den Gewichtsabfall effektiver handhabt, gesetzt mit einer Lernrate von 0,001 und einem Gewichtsabfall von 0,0001, wodurch sowohl die Lerngeschwindigkeit als auch die Modellstabilität optimiert werden.





Dabei sind mt und vt das erste und zweite Moment der Gradienten
und
die bias-korrigierten Momente, η die Lernrate und ε eine kleine Konstante, um eine Division durch Null zu verhindern.
Das Modell verwendet beim Training eine Batchgröße von 32, um die GPU-Beschleunigung für schnellere Berechnungen zu nutzen, und ist so konfiguriert, dass es für 100 Epochen trainiert.
Aber das Training hat einen frühen Stoppmechanismus bei Validierungsverlusten, um das Training zu begrenzen, wenn das Modell aufhört zu verbessern, um Rechenressourcen zu sparen und Übertraining von fünf aufeinanderfolgenden Epochen zu vermeiden. Die Modellleistung wird mit Kennzahlen wie der spärlichen kategorialen Genauigkeit und der Top-2-Genauigkeit erfasst, die über die Klassifikationsfähigkeit des Modells in den wahrscheinlichsten vorhergesagten Kategorien informieren. Rückrufe für das Modelltraining umfassen Checkpointing, das das leistungsstärkste Modell entsprechend der Validierungsgenauigkeit speichert, um sicherzustellen, dass die erfolgreichste Version des Modells nach Abschluss des Trainingsprozesses erhalten bleibt. Diese robuste und gut geplante Architektur nutzt die Fähigkeiten von Transformatoren, um medizinische Bilddaten mit einem hochentwickelten Ansatz zur Nutzung moderner KI-Methoden für bedeutende Anwendungen in der Gesundheitsdiagnostik zu verarbeiten.
4. Ausbildung und Validierung
Das Modell wurde in einer Kaggle-Umgebung mit einer NVIDIA P100 GPU trainiert, und während des Trainingsprozesses wurde ein Mini-Batch-Gradientenabstieg mit einer Batchgröße von 32 verwendet. Der während des Trainings verwendete Optimierer war AdamW mit einer Lernrate von 0,001 und einem Gewichtsabfall von 0,0001; Das war ein gutes Gleichgewicht zwischen einer schnellen Konvergenz und etwas Regularisierung. Das Modell wurde für 100 Epochen trainiert, jedoch wurde frühes Stoppen bei Validierungsverlust mit einer Geduld von 5 Epochen verwendet. Einfach ausgedrückt: Wenn das Training den Validierungsverlust für 5 aufeinanderfolgende Epochen nicht verbessert, würde das Training wegen Überanpassung und Recheneffizienz stoppen. In den meisten Fällen stellte dieses Modell Gewichte aus der Epoche mit dem geringsten Validierungsverlust wieder her, was darauf hindeutet, dass es optimal arbeitete und nicht die gesamten 100 Epochen durchlaufen musste. Insgesamt etwa 32 Minuten, um das Modell zu trainieren.
Während des Trainings umfassten die Kennzahlen eine knappe kategoriale Genauigkeit und Top-2-Genauigkeit, die sowohl auf Trainings- als auch auf Validierungssets überwacht wurden. Diese Kennzahlen geben Einblicke darin, wie oft ein Modell die richtige Klasse vorhersagt und ob die richtige Klasse zu den beiden obersten Vorhersagen gehört, was in medizinischen Anwendungen von großer Bedeutung ist, da Fehlklassifikationen mit hoher Zuverlässigkeit erhebliche Folgen haben können. Die Lernkurven für Verlust und Genauigkeit sind in Abbildung 5 dargestellt.
.
Abbildung 5: Trainings- und Validierungsgenauigkeit sowie Verlustkurven über 50 Epochen, die stabiles Lernen und minimale Überanpassung zeigen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
In der vorliegenden Studie wurde keine Kreuzvalidierung durchgeführt. Dieses Setup ermöglichte effiziente Experimente mit der Modellarchitektur, einschließlich Anpassungen an Transformatorschichten und den Größen der MLP-Köpfe, während das Modell auf unsichtbaren Validierungsdaten gut verallgemeinert wurde.
5. Statistische Analyse
Die Leistung des Vision Transformer-Modells wurde statistisch auf Basis des IQ-OTH/NCCD-Lungenkrebs-Datensatzes analysiert. Präzision, Rückruf, F1-Wert und Genauigkeit wurden mit den Gleichungen 14, 15, 16 und 17 berechnet. Diese vier gelten als konventionelle Maße für Klassifikationsaufgaben, da sie Informationen über die Leistung des Modells bei der Klassifizierung und Klassifizierung pro Klasse für jede Klasse offenbaren können.




Recall ist ein Maß für die Fähigkeit eines Modells, alle relevanten Instanzen einer Klasse zu identifizieren, während Präzision der Anteil der positiven Identifikationen ist, die tatsächlich korrekt waren. Der F1-Score balanciert Erinnerung und Präzision, wenn beides wichtig ist.
Die für diese Aufgabe verwendeten Leistungsmaße waren die Matthews-Korrelationskoeffizient (MCC)-Gleichung 18 und Cohens Kappa-Gleichung 19.


Wobei P0 die beobachtete Übereinstimmung ist, während Pe die erwartete Übereinstimmung ist.
Der MCC ist eine umfassende Messung der Klassifikationsleistung, die sowohl wahre positive als auch negative sowie falsch-positive und falsch-negative Ergebnisse berücksichtigt. Sein Wert kann zwischen -1 und 1 liegen, wobei 1 die perfekte Vorhersage, 0 die zufällige Vorhersage und -1 die vollständige Diskrepanz zwischen den vorhergesagten und wahren Labels darstellt. Das MCC war für Vergleiche verschiedener Modellleistungen wertvoll, wenn es auf unausgewogene Datensätze angewandt wurde, und bot ein Maß, das unabhängig von unterschiedlichen Klassengrößen ausgewogen ist.
Im Rahmen der zusätzlichen statistischen Analyse wurde der F2-Wert mit Priorität auf Abruf berechnet, wie in Gleichung 20 gezeigt.

Die Leistung des Modells wurde außerdem mittelmäßig quadratischer Fehler (), Root Mean Squared Error (RMSE) und Mean Absolute Error (MAE) quantifiziert – typischerweise Messungen für Regressionsaufgaben, die hier aber für die Klassifikationsaufgabe modifiziert werden, um die durchschnittliche Größe des Fehlers ohne Rücksicht auf die Richtung zu quantifizieren.
Um festzustellen, ob es praktikabel wäre, ViTs in verbraucherorientierte Gesundheitsgeräte zu integrieren, führten wir umfangreiche Leistungsbewertungsergebnisse durch, die sich ausschließlich auf die zeitliche Effizienz des Modells konzentrierten. Wir haben Funktionen erstellt, um die Zeit zu berichten, die benötigt wird, um ein einzelnes oder mehrere Bilder vorherzusagen, da dies besonders für Echtzeitanwendungen relevant wird. Für jedes Bild werden vor Beginn der Vorhersage die Daten zunächst so vorverarbeitet, dass sie die Form des vom Modell gewünschten Eingabebildes haben. Die Zeit, zu der die Vorhersage begann, und die Zeit, zu der die Vorhersage abgeschlossen war, haben wir aufgezeichnet, um die Ergebnisse von Zeit und Latenz zu erhalten. Zeit und durchschnittliche Latenz wurden mit Gleichung 21 bzw. 22 berechnet.


Wo:
- N ist die Anzahl der Bilder (Proben).
- Tend ist die Zeit, die nach der Vorhersage des i-ten Bildes aufgezeichnet wird.
- tstart ist die aufgezeichnete Zeit, bevor das i-te Bild vorhergesagt wird.
Weitere Experimente wurden durchgeführt, um die Robustheit des von uns vorgeschlagenen Vision Transformer-Modells zu bewerten, indem systematisch zusätzliches Rauschen und Unschärfe in die Bilder eingeführt wurden. Gaußsches Rauschen wurde jedem Pixel mit einem Rauschfaktor von 0,4 hinzugefügt, während Pixelwerte im Bereich von [0,1] abgeschnitten wurden. Der Unschärfefaktor wurde durch Gaußsche Unschärfe mit einer Kerngröße von 45× 45 verringert. Diese Experimente sind darauf ausgelegt, das Modell unter simulierten wahrnehmbaren Bildqualitätsverschlechterungen umfassend zu bewerten.