$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Il dataset include un totale di 1.190 sezioni di TAC provenienti da 110 casi, suddivise in tre categorie: normali (55 casi), benigni (15 casi) e maligni (40 casi). Ogni caso comprende più tacce CT (che variano da circa 80 a 200 fette per caso), offrendo diverse viste assiali della regione toracica. Le immagini TC sono state ottenute in DICOM utilizzando lo scanner SOMATOM con parametri di imaging standard: tensione del tubo = 120 kV, spessore della fetta = 1 mm, larghezza finestra = 350–1.200 unità di Hounsfield (HU) e valori centrici della finestra = 50–600 HU, durante una piena ricezione di ispirazione.
Tutte le immagini sono state completamente de-identificate prima dell'analisi per rimuovere qualsiasi informazione personale identificabile (PII). Il dataset è stato eticamente approvato dai comitati di revisione istituzionali dei centri medici partecipanti, con il consenso scritto rinunciato dal comitato di supervisione. I casi includevano persone di contesti vari come dipendenti pubblici, agricoltori e residenti di diverse province irachene (tra cui Baghdad, Wasit, Diyala, Salahuddin e Babilonia) con diversità di genere, età, livello educativo e status abitativo.
Poiché il dataset è pubblicamente disponibile e de-identificato, non è stata necessaria un'ulteriore approvazione etica per il suo utilizzo in questo studio. Il dataset rispetta i termini e le condizioni specificati dai suoi contributori originali e dalla piattaforma ospitante.
La metodologia di questa ricerca utilizza un processo multi-fasi progettato per creare un modello previsionale con l'aiuto del dataset IQ-OTH/NCCD sul cancropolmonare 21. Questa metodologia crea un terreno solido per dispositivi sanitari orientati al consumatore, dove è richiesta una latenza minore, il che può garantire una maggiore accuratezza. La Figura 2 mostra il meccanismo di funzionamento del modello proposto.

Figura 2: Diagramma del flusso di lavoro del modello proposto che mostra preprocessing, augmentation, classificazione Vision Transformer e passaggi di valutazione. Clicca qui per visualizzare una versione più grande di questa figura.
1. Descrizione del dataset
Il dataset IQ-OTH/NCCD sul cancro polmonare è stato raccolto presso l'Iraq-Oncology Teaching Hospital/National Center for Cancer Diseases nell'autunno del 2019. La parte di addestramento del dataset utilizzata conteneva 1.097 immagini la cui descrizione è fornita nella Tabella 2, che illustra il conteggio dei campioni di diverse classi.
| Classe | Numero di immagini |
| Normale | 416 |
| Benigno | 120 |
| Maligno | 561 |
Tabella 2: Esempio di immagini TC normali, benigne e maligne dal dataset.
Il dataset IQ-OTH/NCCD sul cancro polmonare è stato scelto per la sua rappresentazione completa delle TAC normali, benigne e maligne. Sebbene siano disponibili altri dataset, come LIDC-IDRI e NSCLC-Radiomics, essi si concentrano principalmente sulla rilevazione dei noduli o mancano di campioni sufficienti di casi benigni. Il dataset IQ-OTH/NCCD è particolarmente adatto al nostro studio poiché consente al Vision-Transformer di apprendere caratteristiche discriminative in tutte e tre le classi, consentendo una classificazione robusta di schemi sottili nelle immagini TC polmonari.
2. Preprocessing dei dati
La pre-elaborazione è un passaggio importante per migliorare le prestazioni del modello attraverso la coerenza e adeguate modifiche ai dati che verranno elaborati attraverso la rete neurale. Per garantire coerenza nella dimensione di input della rete neurale, abbiamo scalato tutte le immagini alla misura di 256 x 256 pixel e normalizzato i dati a un valore di pixel compreso tra 0 e 1 nella speranza di migliorare l'addestramento e la convergenza del modello. La Tabella 3 contiene i valori della tecnica di aumento.
| Tecnica | Valore |
| Normalizzazione | - |
| Ridimensionamento | image_size x image_size |
| RandomRotation | fattore=0,02 |
| RandomZoom | height_factor=0,2, width_factor=0,2 |
Tabella 3: Tecniche di aumento applicate con intervalli di parametri.
Per aumentare la robustezza del modello contro l'overfitting e migliorare la sua capacità di generalizzazione, vengono applicate tecniche di aumento dei dati come rotazioni casuali e zoom, che simulano vari angoli e dimensioni delle manifestazioni del tumore polmonare come potrebbero comparire in diversi pazienti. In questo studio sono state applicate rotazioni casuali con angoli campionati uniformemente da 0,02 radiani e trasformazioni di zoom casuali con fattori di altezza e larghezza variabili. Le immagini successive all'aumento sono mostrate nella Figura 3.

Figura 3: Immagini CT dopo l'aumento che dimostrano rotazione, zoom e normalizzazione per migliorare la generalizzazione del modello. Clicca qui per visualizzare una versione più grande di questa figura.
Queste tecniche di preprocessing sono necessarie e l'aumento è stato utilizzato in tutte le classi per garantire la robustezza del modello.
3. Architettura del modello
Adattando il nuovo framework Vision Transformer (ViT) per gestire efficacemente dati di immagini complessi, l'architettura del modello è pensata per categorizzare le TAC in tre categorie: normale, benigna e maligna. Con questo metodo, vengono elaborate immagini di input da 256 x 256 pixel. Per garantire coerenza e promuovere un apprendimento più efficiente del modello, ogni immagine viene sottoposta a diverse operazioni di preelaborazione, come ridimensionamento e normalizzazione. Per migliorare la resilienza del modello ai cambiamenti di scala e orientamento dell'immagine, il progetto inizia con uno strato di aumento dati che utilizza metodi come normalizzazione, ridimensionamento, rotazioni casuali di 0,02 radiani e zoom casuali fino al 20%. Dopo l'augmentazione, il modello prende 16 patch da 16 pixel da ogni immagine, quindi ogni immagine ha 256 patch.
L'algoritmo 1 definisce il flusso di lavoro, il modello proposto e come viene costruito, nonché la fine tuning che viene eseguita rispetto al modello.
Algoritmo 1: Classificazione del cancro ai polmoni utilizzando Vision Transformers
Input: Set di immagini CT I dal dataset IQ-OTH/NCCD
Output: Risultati della classificazione C in categorie: Normale, Benigno, Maligno
Preelaborazione:
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
Configurazione dei modelli:
Initialize Vision Transformer (ViT) Model
Set number of patches P = 16 x 16
Set number of heads H = 6 in multi-head attention
Formazione:
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
Convalida:
Split data into TrainingSet (80%) and ValidationSet (20%)
ComputeValidationMetrics(ValidationSet)
Valutazione:
C <- Classify(I)
ComputePerformanceMetrics(C)
Return C
Per preservare le relazioni spaziali tra le patch, queste patch vengono appiattite (equazione 1) in vettori di 768 dimensioni (poiché ogni patch ha 16 x 16 x 3 = 768) e poi inviate attraverso uno strato di codifica patch, che proietta ogni vettore in uno spazio a 64 dimensioni integrando embedding posizionali. Il nucleo dell'architettura comprende otto strati trasformer, ciascuno dotato di un meccanismo di attenzione multi-testa con sei testine, permettendo al modello di mettere a fuoco su diverse parti dell'immagine contemporaneamente e catturare un'ampia gamma di caratteristiche. È rappresentato tramite le equazioni 2, 3 e 4.

Dove μ è la media e σ2 è la varianza dell'input x, e ε è una piccola costante per evitare la divisione per zero.

Dove Q è la matrice di interrogazione, K è la matrice chiave, V è la matrice dei valori e dk è la dimensione dei vettori chiave.


Dove WiQ, WiK e WiV sono rispettivamente le matrici di pesi apprese per query, chiavi e valori, e WO è la matrice di pesi in uscita.
Il diagramma a blocchi del modello proposto è stato mostrato nella Figura 4.

Figura 4: Diagramma a blocchi del modello Vision Transformer con testa MLP, che dettaglia i principali strati di elaborazione e l'architettura. Clicca qui per visualizzare una versione più grande di questa figura.
Ogni layer di trasformatore include una percezione multilivello (MLP) con unità nascoste inizialmente scalate fino a 128 e poi di nuovo a 64, espandendo e comprimendo efficacemente il flusso di informazioni per catturare dipendenze complesse.
Il dropout viene applicato strategicamente a un tasso di 0,1 nei meccanismi di attenzione e nelle MLP per prevenire il sovrafitting. L'uscita dall'encoder del trasformatore viene elaborata attraverso una testa MLP composta da due strati densi rispettivamente da 2.048 e 1.024 unità, impiegando Unità Lineari di Errore Gaussiano (GELU) per l'attivazione, che si è dimostrato efficace nelle applicazioni di deep learning. Questo si ottiene usando l'equazione 5.

Dove W1 e W2 sono matrici di peso, b1 e b2 sono bias e GELU è la funzione di attivazione utilizzata.
Un dropout di 0,1 veniva utilizzato negli strati trasformatori per evitare il sovrafitting senza perdere informazioni importanti sulle caratteristiche. La funzione di attivazione GELU è stata utilizzata per le sue caratteristiche non lineari e lisce, che promuovono il flusso gradiente e la convergenza nei modelli basati su trasformatori. Il dropout viene regolarizzato usando l'equazione 6.
Dove p è il tasso di abbandono (ad esempio, 0,1 o 0,5), e lo strato di abbandono imposta casualmente una frazione p delle unità di input a zero durante l'addestramento.
Un tasso di dropout di 0,5 in questi strati aiuta ulteriormente a mitigare il sovrafitting. L'ultimo livello del modello è uno strato logits (equazione 7) che produce tre logit di classe corrispondenti alle categorie normale, benigna e maligna, utilizzando una funzione di perdita di entropia incrociata categorica sparsa (equazione 8) appropriata per questo contesto di classificazione multiclasse.


Dove zi è il vettore logits per la classe i, SoftMax(zi)) rappresenta le probabilità previste, e yi è la vera etichetta di classe.
Il modello è compilato con l'ottimizzatore AdamW (equazioni 9, 10, 11, 12 e 13), un'estensione dell'ottimizzatore Adam che gestisce in modo più efficace il decadimento del peso, impostato a un tasso di apprendimento di 0,001 e un decadimento del peso di 0,0001, ottimizzando sia la velocità di apprendimento che la stabilità del modello.





Dove mt e vt sono il primo e il secondo momento dei gradienti,
e
sono i momenti corretti per bias, η è il tasso di apprendimento, e ε è una piccola costante per evitare la divisione per zero.
Il modello, una volta addestrato, utilizza una dimensione di batch di 32 per sfruttare l'accelerazione GPU e un calcolo più rapido ed è configurato per addestrarsi per 100 epoche.
Ma l'addestramento ha un meccanismo di arresto precoce sulla perdita di validazione per limitare l'addestramento quando il modello smette di migliorare, per risparmiare risorse computazionali ed evitare il sovraaddestramento di 5 epoche consecutive. Le prestazioni del modello vengono monitorate con metriche come l'accuratezza categoriale scarsa e l'accuratezza delle prime 2, che informano sulla capacità di classificazione del modello nelle categorie previste più probabili. I callback per l'addestramento del modello comprendono checkpointing che salverà il modello con le prestazioni più alte in base all'accuratezza della validazione, per garantire che la versione più efficace del modello venga mantenuta una volta terminato il processo di addestramento. Questa architettura robusta e ben pianificata sfrutta le capacità dei transformer per elaborare dati di imaging medico con un approccio altamente avanzato all'utilizzo dei metodi IA contemporanei per applicazioni importanti nella diagnostica sanitaria.
4. Formazione e validazione
Il modello è stato addestrato in un ambiente Kaggle utilizzando una GPU NVIDIA P100, e durante il processo di addestramento è stato utilizzato un mini-batch gradient descent con una dimensione di 32 lotti. L'ottimizzatore utilizzato durante l'addestramento era AdamW con un tasso di apprendimento di 0,001 e un decadimento del peso di 0,0001; Questo era un buon equilibrio tra una rapida convergenza e una certa regolarizzazione. Il modello è stato addestrato per 100 epoche, tuttavia si è utilizzato l'arresto precoce sulla perdita di validazione con una pazienza di 5 epoche. Per dirla semplicemente, se l'addestramento non migliorava la perdita di validazione per 5 epoche continue, l'addestramento si fermava per sovrafitting e efficienza computazionale. Nella maggior parte dei casi, questo modello ripristinava i pesi dell'epoca con la minima perdita di validazione, indicando che stava funzionando in modo ottimale e non aveva bisogno di eseguire l'intera epoca di 100 epoche. Un totale di circa 32 minuti per addestrare il modello.
Durante l'addestramento, le metriche includevano un'accuratezza categoriale scarsa e una precisione top 2, monitorate sia su set di addestramento che su quelli di validazione. Queste metriche forniscono informazioni su quanto spesso un modello prevede la classe corretta e se la classe corretta si trova tra le prime due previsioni, cosa che nelle applicazioni mediche è importante perché le errate classificazioni ad alta fiducia possono avere conseguenze significative. Le curve di apprendimento per perdita e accuratezza sono visualizzate nella Figura 5.
.
Figura 5: Accuratezza dell'addestramento e della validazione e curve di perdita in 50 epoche che mostrano un apprendimento stabile e un overfitting minimo. Clicca qui per visualizzare una versione più grande di questa figura.
Non è stata effettuata alcuna validazione incrociata nel presente studio. Questa configurazione permise una sperimentazione efficiente dell'architettura del modello, inclusi aggiustamenti ai livelli dei transformer e alle dimensioni delle testine MLP, assicurando al contempo che il modello si generalizzasse bene su dati di validazione non visti.
5. Analisi statistica
Le prestazioni del modello Vision Transformer sono state analizzate statisticamente sulla base del dataset IQ-OTH/NCCD sul cancro ai polmoni. Precisione, richiamo, punteggio F1 e accuratezza sono stati calcolati utilizzando rispettivamente le equazioni 14, 15, 16 e 17. Queste quattro sono considerate le misure convenzionali da utilizzare per compiti di classificazione, in quanto possono rivelare informazioni sulle prestazioni del modello nella classificazione e classificazione per classe per ogni classe.




Il ricordo è una misura della capacità di un modello di identificare tutte le istanze rilevanti in una classe, mentre la precisione è la proporzione di identificazioni positive che erano effettivamente corrette. Il punteggio F1 bilancia richiamo e precisione quando entrambe sono importanti.
Le misure di performance utilizzate per questo compito sono state l'equazione 18 del Coefficiente di Correlazione di Matthews (MCC) e l'equazione di Kappa 19 di Cohen.


Dove P0 è l'accordo osservato mentre Pe è l'accordo atteso.
Il MCC è una misura onnicomprensiva delle prestazioni nella classificazione, tenendo conto sia dei veri positivi che dei negativi, dei falsi positivi e dei falsi negativi. Il suo valore può essere compreso tra -1 e 1, con 1 che rappresenta la previsione perfetta, 0 la previsione casuale e -1 il completo disaccordo tra le etichette prevista e vera. Il MCC era utile per confronti tra diverse prestazioni di modelli quando applicato a dataset sbilanciati e fornendo una misura bilanciata indipendentemente dalle diverse dimensioni delle classi.
Come parte di un'ulteriore analisi statistica, il punteggio F2, che dà priorità al richiamo, è stato calcolato come mostrato dall'equazione 20.

Le prestazioni del modello sono state anche quantificate utilizzando l'Errore Quadratico Medio (MSE), l'Errore Quadratico Medio Radice (RMSE) e l'Errore Assoluto Medio (MAE) - che sono tipicamente misurazioni per compiti di regressione ma modificate qui per il compito di classificazione per quantificare quanto fosse grande in media l'errore senza considerare la direzione.
Per determinare se fosse pratico integrare i ViT nei dispositivi sanitari orientati al consumatore, abbiamo effettuato ampi risultati di valutazione delle prestazioni concentrandoci esclusivamente sull'efficienza temporale del modello. Abbiamo creato funzioni per riportare il tempo necessario per prevedere una singola o semplicemente più immagini, poiché questo diventa particolarmente rilevante per applicazioni in tempo reale. Per ogni immagine, prima di iniziare a prevedere, i dati vengono prima pre-elaborati per assumere la forma dell'immagine di input desiderata dal modello. L'ora in cui la previsione è iniziata e quella completata la previsione li abbiamo registrati per ottenere i risultati di tempo e latenza. Il tempo e la latenza media sono stati calcolati utilizzando rispettivamente l'equazione 21 e l'equazione 22.


Dove:
- N è il numero di immagini (campioni).
- Tend è il tempo registrato dopo la previsione dell'i-esima immagine.
- Tstart è il tempo registrato prima di prevedere la i-esima immagine.
Ulteriori esperimenti sono stati condotti per valutare la robustezza del modello Vision Transformer che abbiamo proposto, introducendo sistematicamente ulteriore rumore e sfocatura nelle immagini. Il rumore gaussiano veniva aggiunto a ogni pixel con un fattore di rumore di 0,4 mentre i valori dei pixel erano ritagliati nell'intervallo di [0,1]. Il fattore di sfocatura è stato ridotto tramite la sfocatura gaussiana con una dimensione del chicco di 45× 45. Questi esperimenti sono progettati per valutare il modello in modo completo sotto degradazioni simulate della qualità percettiva dell'immagine.