$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
De dataset bevat in totaal 1.190 CT-scansnedes uit 110 gevallen, onderverdeeld in drie klassen: normaal (55 gevallen), goedaardig (15 gevallen) en kwaadaardig (40 gevallen). Elke casus bestaat uit meerdere CT-sneden (variërend van ongeveer 80 tot 200 sneden per geval), die diverse axiale aanzichten van het thoracale gebied bieden. De CT-beelden werden verkregen in DICOM met behulp van de SOMATOM-scanner en standaard beeldparameters: buisspanning = 120 kV, slicedikte = 1 mm, vensterbreedte = 350–1.200 Hounsfield-units (HU), en venstermiddenwaarden = 50–600 HU, tijdens een volledige inademhaling.
Alle beelden werden volledig gedeïdenificeerd vóór de analyse om persoonlijke identificeerbare informatie (PII) te verwijderen. De dataset werd ethisch goedgekeurd door de institutionele beoordelingscommissies van de deelnemende medische centra, waarbij schriftelijke toestemming werd opgegeven door de toezichtscommissie. De gevallen betroffen personen uit verschillende achtergronden, zoals overheidsmedewerkers, boeren en inwoners uit verschillende Iraakse provincies (waaronder Bagdad, Wasit, Diyala, Salahuddin en Babylon) met diversiteit in geslacht, leeftijd, opleidingsniveau en levensstatus.
Omdat de dataset openbaar beschikbaar is en gedeidentificeerd, was er geen extra ethische goedkeuring nodig voor het gebruik ervan in deze studie. De dataset voldoet aan de voorwaarden die zijn gespecificeerd door de oorspronkelijke bijdragers en het hostingplatform.
De methodologie van dit onderzoek maakt gebruik van een meerfasenproces dat is ontworpen om een voorspellingsmodel te creëren met behulp van de IQ-OTH/NCCD Longkankerdataset21. Deze methodologie vormt een sterk fundament voor consumentgerichte zorgapparaten, waarbij minder latentie nodig is, wat voor betere nauwkeurigheid kan zorgen. Figuur 2 toont het werkingsmechanisme van het voorgestelde model.

Figuur 2: Workflowdiagram van het voorgestelde model met preprocessing, augmentatie, Vision Transformer-classificatie en evaluatiestappen. Klik hier om een grotere versie van deze figuur te bekijken.
1. Datasetbeschrijving
De IQ-OTH/NCCD-longkankerdataset werd verzameld in het Iraq-Oncology Teaching Hospital/National Center for Cancer Diseases in het najaar van 2019. Het trainingsgedeelte van de gebruikte dataset bevatte 1.097 afbeeldingen waarvan de beschrijving in Tabel 2 wordt gegeven, die het aantal voorbeelden van verschillende klassen illustreert.
| Klas | Aantal afbeeldingen |
| Normaal | 416 |
| Goedaardig | 120 |
| Kwaadaardig | 561 |
Tabel 2: Voorbeelden van normale, goedaardige en kwaadaardige CT-beelden uit de dataset.
De IQ-OTH/NCCD-longkankerdataset is gekozen vanwege de uitgebreide weergave van normale, goedaardige en kwaadaardige CT-scans. Hoewel andere datasets, zoals LIDC-IDRI en NSCLC-Radiomics, beschikbaar zijn, richten deze zich ofwel voornamelijk op noduledetectie of ontbreken ze voldoende monsters van goedaardige gevallen. De IQ-OTH/NCCD-dataset is bijzonder geschikt voor ons onderzoek omdat het de Vision Transformer in staat stelt om onderscheidende kenmerken te leren over alle drie de klassen, waardoor subtiele patronen in long-CT-beelden robuust kunnen worden geclassificeerd.
2. Gegevensvoorverwerking
Preprocessing is een belangrijke stap in het verbeteren van modelprestaties door consistentie en passende aanpassingen aan de data die via het neurale netwerk verwerkt zullen worden. Om consistentie in de invoergrootte voor het neurale netwerk te waarborgen, hebben we alle afbeeldingen geschaald tot dezelfde dimensie van 256 x 256 pixels, en de data genormaliseerd naar een pixelwaarde tussen 0 en 1 in de hoop modeltraining en convergentie te verbeteren. Tabel 3 bevat de waarden van de augmentatietechniek.
| Techniek | Waarde |
| Normalisatie | - |
| Resizing | image_size x image_size |
| RandomRotation | factor=0,02 |
| RandomZoom | height_factor=0,2, width_factor=0,2 |
Tabel 3: Toegepaste augmentatietechnieken met parameterbereiken.
Om de robuustheid van het model tegen overfitting te vergroten en het vermogen om te generaliseren te verbeteren, worden data-augmentatietechnieken zoals willekeurige rotaties en zooms toegepast, waarbij verschillende hoeken en groottes van longkankermanifestaties worden gesimuleerd zoals ze bij verschillende patiënten kunnen voorkomen. Willekeurige rotaties met hoeken die uniform werden bemonsterd uit 0,02 radialen en willekeurige zoomtransformaties met wisselende hoogte- en breedtefactoren werden in deze studie toegepast. De beelden na de augmentatie zijn weergegeven in Figuur 3.

Figuur 3: CT-beelden na augmentatie die rotatie, zoom en normalisatie tonen om de generalisatie van het model te verbeteren. Klik hier om een grotere versie van deze figuur te bekijken.
Deze preprocessingtechnieken zijn vereist en augmentatie is gebruikt over alle klassen om de robuustheid van het model te waarborgen.
3. Modelarchitectuur
Door het nieuwe Vision Transformer (ViT)-framework aan te passen om effectief met ingewikkelde beeldgegevens om te gaan, is de modelarchitectuur bedoeld om CT-scans in drie categorieën te categoriseren: normaal, goedaardig en kwaadaardig. Met deze methode worden invoerbeelden van 256 x 256 pixels verwerkt. Om consistentie te garanderen en efficiënter modelleren te bevorderen, wordt elke afbeelding door verschillende preprocessing-operaties gevoerd, zoals het vergroten van grootte en normalisatie. Om de veerkracht van het model tegen veranderingen in beeldschaal en -oriëntatie te verbeteren, begint het ontwerp met een data-augmentatielaag die methoden gebruikt zoals normalisatie, verkleining, willekeurige rotaties van 0,02 radialen en willekeurige zooms tot 20%. Na de augmentatie haalt het model 16 x 16 pixelpatches uit elke afbeelding, waardoor elke afbeelding 256 patches heeft.
Algoritme 1 definieert de workflow, het voorgestelde model, hoe het wordt opgebouwd, en de fijne afstemming die wordt uitgevoerd ten opzichte van het model.
Algoritme 1: Classificatie van longkanker met behulp van Vision Transformers
Invoer: Set CT-beelden I uit de IQ-OTH/NCCD-dataset
Output: Classificatieresultaten C in categorieën: Normaal, Goedaardig, Kwaadaardig
Preprocessing:
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
Modelopstelling:
Initialize Vision Transformer (ViT) Model
Set number of patches P = 16 x 16
Set number of heads H = 6 in multi-head attention
Opleiding:
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
Validering:
Split data into TrainingSet (80%) and ValidationSet (20%)
ComputeValidationMetrics(ValidationSet)
Evaluatie:
C <- Classify(I)
ComputePerformanceMetrics(C)
Return C
Om de ruimtelijke relaties tussen patches te behouden, worden deze patches afgevlakt (vergelijking 1) tot 768-dimensionale vectoren (omdat elke patch 16 x 16 x 3 = 768 heeft) en vervolgens door een patch-coderinglaag gestuurd, die elke vector projecteert in een 64-dimensionale ruimte door positionele embeddings te integreren. De kern van de architectuur bestaat uit acht transformatorlagen, elk met een multi-head aandachtmechanisme met zes koppen, waardoor het model gelijktijdig op verschillende delen van het beeld kan focussen en een breed scala aan kenmerken kan vastleggen. Het wordt weergegeven met behulp van vergelijkingen 2, 3 en 4.

Waarbij μ het gemiddelde is en σ2 de variantie van de invoer x, en ε een kleine constante is om deling door nul te voorkomen.

Waarbij Q de querymatrix is, K de sleutelmatrix, V de waardematrix en dk de dimensie van de sleutelvectoren.


Waarbij WiQ, WiK en WiV de geleerde gewichtmatrices zijn voor respectievelijk zoekopdrachten, sleutels en waarden, en WO de outputgewichtmatrix is.
Het blokdiagram van het voorgestelde model is weergegeven in Figuur 4.

Figuur 4: Blokdiagram van het Vision Transformer-model met MLP-kop, waarin de belangrijkste verwerkingslagen en architectuur worden beschreven. Klik hier om een grotere versie van deze figuur te bekijken.
Elke transformatorlaag bevat een multi-layer perception (MLP) met verborgen eenheden die eerst worden opgeschaald tot 128 en daarna weer terug naar 64, waardoor de informatiestroom effectief wordt uitgebreid en gecomprimeerd om complexe afhankelijkheden te vangen.
Drop-out wordt strategisch toegepast met een frequentie van 0,1 binnen de aandachtsmechanismen en MLP's om overfitting te voorkomen. De uitgang van de transformatorencoder wordt verwerkt via een MLP-kop bestaande uit twee dichte lagen met respectievelijk 2.048 en 1.024 eenheden, waarbij gebruik wordt gemaakt van Gaussian Error Linear Units (GELU) voor activatie, wat goed is aangetoond te presteren in deep learning-toepassingen. Dit wordt bereikt met behulp van vergelijking 5.

Waar W1 en W2 gewichtsmatrices zijn, b1 en b2 biases, en GELU de gebruikte activatiefunctie is.
Een dropout van 0,1 werd gebruikt in de transformatorlagen om overfitting te voorkomen zonder belangrijke feature-informatie te verliezen. De GELU-activatiefunctie werd gebruikt vanwege zijn gladde niet-lineaire eigenschappen, die gradiëntstroom en convergentie bevorderen in transformator-gebaseerde modellen. De uitval wordt geregulariseerd met behulp van vergelijking 6.
Waarbij p het uitvalpercentage is (bijvoorbeeld 0,1 of 0,5), en de uitvallaag een fractie p van de invoereenheden willekeurig op nul zet tijdens de training.
Een uitvalpercentage van 0,5 in deze lagen helpt verder bij het beperken van overfitting. De laatste laag van het model is een logitslaag (vergelijking 7) die drie klasse-logits produceert die overeenkomen met de categorieën normaal, goedaardig en kwaadaardig, met behulp van een spaar categorische kruisentropieverliesfunctie (vergelijking 8) die geschikt is voor deze multi-klasse classificatie.


Waarbij zi de logitsvector is voor klasse i, SoftMax(zi)) de voorspelde kansen voorstelt, en yi het ware klasselabel.
Het model wordt gecompileerd met de AdamW-optimizer (vergelijkingen 9, 10, 11, 12 en 13), een uitbreiding van de Adam-optimizer die gewichtsafname effectiever aankan, ingesteld met een leersnelheid van 0,001 en een gewichtsafname van 0,0001, waarmee zowel de leersnelheid als de stabiliteit van het model worden geoptimaliseerd.





Waarbij mt en vt het eerste en tweede moment van de gradiënten zijn,
en
de bias-gecorrigeerde momenten, η de leersnelheid is, en ε een kleine constante is om deling door nul te voorkomen.
Het model gebruikt bij training een batchgrootte van 32 om gebruik te maken van GPU-versnelling voor snellere berekeningen en is geconfigureerd om te trainen voor 100 epochs.
Maar de training heeft een vroegtijdig stopmechanisme bij validatieverlies om de training te beperken wanneer het model stopt met verbeteren, om rekenkracht te besparen en overtraining van 5 opeenvolgende epochs te voorkomen. De prestaties van het model worden gevolgd met meetwaarden zoals spaarzame categorische nauwkeurigheid en top 2 nauwkeurigheid, die informeren over de classificatiecapaciteit van het model over de meest waarschijnlijke voorspelde categorieën. Callbacks voor modeltraining omvatten checkpointing waarbij het best presterende model wordt opgeslagen op basis van validatienauwkeurigheid, zodat de meest succesvolle versie van het model behouden blijft zodra het trainingsproces is afgerond. Deze robuuste en goed geplande architectuur maakt gebruik van de mogelijkheden van transformatoren om medische beeldvormingsdata te verwerken met een zeer geavanceerde benadering om hedendaagse AI-methoden te gebruiken voor belangrijke toepassingen in de gezondheidszorgdiagnostiek.
4. Training en validatie
Het model werd getraind in een Kaggle-omgeving met een NVIDIA P100 GPU, en tijdens het trainingsproces werd een mini-batch gradient descent met een batchgrootte van 32 gebruikt tijdens het trainingsproces. De optimizer die tijdens de training werd gebruikt was AdamW met een leersnelheid van 0,001 en een gewichtsafname van 0,0001; Dit was een goede balans tussen een snelle convergentie en enige regularisatie. Het model werd getraind voor 100 epochs, maar vroegtijdig stopte bij validatieverlies met een geduld van 5 epochs. Simpel gezegd: als de training validatieverlies niet verbetert voor 5 opeenvolgende epochs, zou de training stoppen vanwege overfitting en rekenefficiëntie. In de meeste gevallen herstelde dit model gewichten uit het tijdperk met het minste validatieverlies, wat aangeeft dat het optimaal presteerde en niet de volledige 100 epochen hoefde te lopen. In totaal ongeveer 32 minuten om het model te trainen.
Tijdens de training omvatten de meetwaarden een schaarse categorische nauwkeurigheid en top-2 nauwkeurigheid, gemonitord op zowel trainings- als validatiesets. Deze metrics geven inzicht in hoe vaak een model de juiste klasse voorspelt en of de juiste klasse tot de top twee voorspellingen behoort, wat in medische toepassingen belangrijk is omdat misclassificaties met hoge betrouwbaarheid aanzienlijke gevolgen kunnen hebben. De leercurves voor verlies en nauwkeurigheid zijn zichtbaar in Figuur 5.
.
Figuur 5: Trainings- en validatienauwkeurigheid en verliescurves over 50 tijdperken, wat stabiel leren en minimale overfitting laat zien. Klik hier om een grotere versie van deze figuur te bekijken.
In de huidige studie werd geen kruisvalidatie uitgevoerd. Deze opstelling maakte efficiënt experimenteren met de modelarchitectuur mogelijk, inclusief aanpassingen aan transformatorlagen en de afmetingen van MLP-koppen, terwijl het model goed gegeneraliseerd werd op onzichtbare validatiegegevens.
5. Statistische analyse
De prestaties van het Vision Transformer-model werden statistisch geanalyseerd op basis van de IQ-OTH/NCCD-longkankerdataset. Precisie, herinnering, F1-score en nauwkeurigheid werden respectievelijk berekend met vergelijkingen 14, 15, 16, 17. Deze vier worden beschouwd als de conventionele maatstaven voor classificatietaken, omdat ze informatie kunnen onthullen over de prestaties van het model bij classificatie en classificatie per klasse voor elke klasse.




Recall is een maat voor het vermogen van een model om alle relevante gevallen in een klasse te identificeren, terwijl precisie het aandeel positieve identificaties is dat daadwerkelijk correct was. De F1-score balanceert herinnering en precisie wanneer beide belangrijk zijn.
De prestatiematen die voor deze opdracht werden gebruikt, waren de Matthews Correlatiecoëfficiënt (MCC)-vergelijking 18 en Cohen's Kappa-vergelijking 19.


Waarbij P0 de waargenomen overeenkomst is en Pe de verwachte overeenkomst.
De MCC is een allesomvattende maatstaf voor classificatieprestaties, waarbij zowel ware positieven als negatieven, vals-positieven en vals-negatief worden meegenomen. De waarde kan tussen -1 en 1 liggen, waarbij 1 perfecte voorspelling is, 0 willekeurige voorspelling, en -1 volledige onenigheid tussen de voorspelde en ware labels. De MCC was waardevol voor vergelijkingen tussen verschillende modelprestaties wanneer toegepast op onevenwichtige datasets en biedt een maat die in balans is, ongeacht verschillende klassengroottes.
Als onderdeel van aanvullende statistische analyse werd de F2-score, met prioriteit voor terugroepen, berekend zoals weergegeven in vergelijking 20.

De prestaties van het model werden ook gekwantificeerd met behulp van Mean Squared Error (MSE), Root Mean Squared Error (RMSE) en Mean Absolute Error (MAE) – dit zijn typisch metingen voor regressietaken, maar hier aangepast voor de classificatietaak om te kwantificeren hoe groot de gemiddelde fout was, zonder rekening te houden met richting.
Om te bepalen of het praktisch zou zijn om ViTs te integreren in consumentgerichte zorgapparatuur, hebben we uitgebreide prestatie-evaluatieresultaten uitgevoerd die zich uitsluitend richtten op de temporele efficiëntie van het model. We hebben functies gemaakt om de tijd te rapporteren die nodig is om één of meerdere afbeeldingen te voorspellen, omdat dit vooral relevant wordt voor realtime toepassingen. Voor elke afbeelding wordt de data vóór het beginnen met voorspellen eerst voorbewerkt zodat deze de vorm heeft van het door het model gewenste invoerbeeld. Het tijdstip waarop de voorspelling begon en het moment waarop de voorspelling voltooid was, registreerden we om de resultaten van tijd en latentie te verkrijgen. Tijd en gemiddelde latentie werden berekend met behulp van respectievelijk vergelijking 21 en vergelijking 22.


Waar:
- N is het aantal afbeeldingen (samples).
- Tend is de tijd die wordt vastgelegd na het voorspellen van het i-de beeld.
- tstart is de tijd die wordt vastgelegd voordat het i-de beeld wordt voorspeld.
Verdere experimenten werden uitgevoerd om de robuustheid van het Vision Transformer-model dat we voorstelden te evalueren door systematisch extra ruis en vervaging aan de beelden toe te voegen. Gaussische ruis werd toegevoegd aan elke pixel met een ruisfactor van 0,4 terwijl pixelwaarden in het bereik van [0,1] werden afgeknipt. De vervagingsfactor werd verminderd door Gaussische blur met een kernelgrootte van 45× 45. Deze experimenten zijn ontworpen om het model volledig te evalueren onder gesimuleerde waarneembare beeldkwaliteitsverlagingen.