$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Questo studio non ha coinvolto alcuna sperimentazione diretta su partecipanti umani o animali. Tutto il lavoro è stato condotto utilizzando il dataset pubblico di LC25000 anonimizzato di immagini istopatologiche, che non conteneva informazioni identificabili sui pazienti né manipolazione diretta di tessuti umani. Non era richiesta l'approvazione del Comitato di Revisione Istituzionale (IRB) o del Comitato per la Cura e l'Uso degli Animali Istituzionali (IACUC). Tutte le procedure rispettavano gli standard etici e rispettavano i termini di utilizzo del dataset per la ricerca accademica. La Figura 2 mostra i passaggi del diagramma del flusso di lavoro.

Figura 2: Flusso di lavoro del metodo proposto. Il flusso di lavoro include pre-elaborazione dei dati, aumentazione, addestramento del modello e valutazione. Clicca qui per visualizzare una versione più grande di questa figura.
Descrizione del dataset
Il dataset LC25000 è stato utilizzato per questo studio. Era composto da 25.000 immagini istopatologiche, tutte uniformmente formattate come file JPEG con una risoluzione di 768 × 768 pixel. Il dataset ha avuto origine da un set iniziale di 1.250 immagini originali, inclusi 250 immagini di tessuto polmonare benigno, 250 adenocarcinomi polmonari, 250 carcinomi squamosi polmonari, 250 immagini di tessuto colonico benigno e 250 adenocarcinomi del colon. Le immagini rimanenti sono state generate tramite aumento dei dati per ampliare il dataset, garantendo una collezione diversificata ed estesa per un addestramento e una validazione robusta del modello.
Le immagini ad alta risoluzione permisero un esame dettagliato delle strutture cellulari, essenziale per una classificazione accurata del cancro. L'aumento dei dati era necessario a causa del numero limitato di immagini originali; Tecniche come rotazioni, flip, zoom e spostamenti sono state impiegate per aumentare sinteticamente la dimensione e la diversità del dataset. Questo processo ha generato ulteriori 23.750 immagini sintetiche in modo bilanciato tra classi, risultando in un dataset finale di 25.000 immagini, con ogni classe contenente 5.000 immagini dopo l'augmentazione.
Preprocessing dei dati
Le immagini istopatologiche grezze all'interno del dataset avevano una risoluzione di 768 × 768 pixel. Per garantire la compatibilità con l'architettura EfficientNetB7 e massimizzare l'efficienza computazionale, ogni immagine veniva ridimensionata da 768 × 768 a 224 × 224 pixel. Questo processo di ridimensionamento rappresentava un compromesso tra la conservazione delle caratteristiche istologiche essenziali e la riduzione del sovraccarico computazionale. L'operazione di ridimensionamento fu rappresentata matematicamente come mostrato in
Equazione 1:
(1)
dove X indicava l'immagine originale, Xridimensionato era l'immagine ridimensionata, e h e w erano rispettivamente l'altezza e la larghezza desiderate. Tutti i valori dei pixel sono stati normalizzati in un intervallo da 0 a 1, il che ha stabilizzato e accelerato l'addestramento assicurando una distribuzione costante degli input tra le immagini. Questa normalizzazione è stata calcolata matematicamente come mostrato nell'Equazione 2. Questo passaggio era essenziale per migliorare la convergenza dei modelli durante l'addestramento.
(2)
L'aumento dei dati è stato impiegato come tecnica critica per migliorare la robustezza e l'applicabilità dei modelli di deep learning, specialmente nel contesto dell'imaging medico con dataset limitati e sbilanciati. In questo studio, sono state applicate diverse tecniche di aumento utilizzando il modulo ImageDataGenerator di TensorFlow per simulare la variazione in vivo nelle immagini istopatologiche. Queste tecniche includevano rotazioni casuali fino a 20 gradi per simulare la variabilità dell'orientamento dei tessuti, spostamenti orizzontali e verticali fino al 20% delle dimensioni dell'immagine per imitare la variabilità della posizione durante la preparazione delle diapositive, e zoom casuali fino al 20% per riflettere la variabilità nei livelli di ingrandimento.
Inoltre, è stato eseguito il flipping orizzontale per introdurre variazioni a immagine speculare, arricchendo ulteriormente il dataset di addestramento. Per gestire i pixel appena sintetizzati generati durante queste trasformazioni, è stata adottata una strategia di riempimento del vicino più vicino per garantire transizioni fluide e preservare l'integrità delle caratteristiche istologiche chiave. La matrice di rotazione 2D usata per l'aumento è rappresentata nell'Equazione 3, e l'Equazione 4 definisce una funzione con un cambiamento incrementale:
(3)
(4)
Queste tecniche di aumento compensavano sfide come lo squilibrio di classe, lievi variazioni istologiche e limitazioni del dataset permettendo al modello di apprendere le caratteristiche invarianti e riducendo il rischio di sovrafitting. Espandendo sinteticamente il dataset e incorporando variabilità controllata, l'aumento dei dati ha migliorato significativamente la capacità del modello di generalizzare a dati nuovi e mai visti, producendo così un sistema più forte e robusto per l'impiego clinico nella diagnosi del cancro ai polmoni e al colon.
Architettura del modello
Il modello è stato sviluppato sulla base dell'architettura EfficientNetB7, una rete di deep learning all'avanguardia rinomata per le sue prestazioni, scalabilità ed efficacia nei compiti di classificazione delle immagini. EfficientNetB7 è stato costruito utilizzando una serie di blocchi Mobile Inverted Bottleneck Convolution (MBConv), che incorporavano convoluzioni separabili in profondità. Questo design strutturale comprimeva efficacemente il modello, portando a riduzioni significative del numero di parametri con una perdita minima di prestazioni. Tale innovazione ha permesso al modello di raggiungere un'elevata accuratezza con costi computazionali inferiori, rendendolo particolarmente adatto ad applicazioni ad alta intensità computazionale come l'analisi istopatologica delle immagini. La Figura 3 mostra l'architettura del modello dell'EfficientNetB7 pre-addestrato.

Figura 3: Architettura del modello di EfficientNetB7. Questa figura mostra la struttura dettagliata del modello EfficientNetB7, illustrando i suoi principali livelli e i blocchi funzionali utilizzati per la classificazione delle immagini. Clicca qui per visualizzare una versione più grande di questa figura.
Il modello accettava immagini istopatologiche pre-processate che erano state ridimensionate dalle loro dimensioni originali di 768 × 768 pixel a 224 × 224 pixel per soddisfare i requisiti di input di EfficientNetB7. Ogni immagine ha subito diversi processi all'interno della rete, tra cui la normalizzazione batch, le funzioni di attivazione Swish e le operazioni di convoluzione. La normalizzazione batch è stata calcolata come mostrato nell'Equazione 5. La funzione di attivazione Swish ha introdotto la non linearità e ha permesso un addestramento più rapido ed efficiente rispetto alle funzioni tradizionali come ReLU. Questi livelli sono stati progettati per estrarre e migliorare le caratteristiche gerarchiche, permettendo al modello di apprendere progressivamente da pattern di basso livello a quelli di alto livello man mano che le immagini si propagavano nella rete.
(5)
La spina dorsale convoluzionale di EfficientNetB7, pre-addestrata sul dataset ImageNet, forniva un insieme altamente generalizzabile di rivelatori di caratteristiche efficaci su un'ampia gamma di domini di immagine. La base pre-addestrata ha svolto un ruolo chiave nel catturare schemi sottili ma informativi nelle immagini istopatologiche, spesso caratterizzate da un'elevata variabilità intra-classe e da caratteristiche a grana fine. L'output dei livelli convoluzionali veniva passato a un livello globale di pooling medio, che riduceva le dimensioni spaziali a un singolo vettore di caratteristiche per ogni immagine. Questo vettore, che racchiude le caratteristiche più rilevanti apprese dagli strati convoluzionali, fu poi appiattito in un array unidimensionale per la compatibilità con i livelli completamente connessi. L'algoritmo per una rilevazione migliorata del cancro ai polmoni e al colon è fornito nel File Supplementare 1.
L'architettura di rete successiva includeva due strati densi (completamente connessi) con rispettivamente 128 e 64 unità. Entrambi i livelli utilizzavano la funzione di attivazione Rettificata Linear Unit (ReLU) per introdurre ulteriore non linearità e per perfezionare le caratteristiche estratte, permettendo al modello di riconoscere schemi e relazioni più precise all'interno dei dati. L'ultimo livello del modello era un livello di output SoftMax con cinque neuroni, ciascuno rappresentante una delle cinque classi del dataset. La funzione SoftMax produceva una distribuzione di probabilità tra le classi, che rifletteva la previsione del modello per ogni immagine di input, come descritto nell'Equazione 6:
(6)
La combinazione di tecniche convoluzionali all'avanguardia, una regolarizzazione robusta e strategie di ottimizzazione efficaci ha portato a un modello di deep learning ben adatto alla classificazione delle immagini mediche ad alto rischio. Il modello ha raggiunto precisione, scalabilità ed efficienza operativa, dimostrando un forte potenziale per l'applicazione clinica nella diagnosi del cancro ai polmoni e al colon. All'interno di questo contesto, la regolarizzazione è stata realizzata principalmente tramite potenziamento avanzato dei dati e l'uso di stop precoce; non è stato applicato alcun caduto esplicito o decadimento del peso sugli strati densi.
Formazione e validazione
La pipeline di validazione e preaddestramento della rete neurale convoluzionale (CNN) per LC25000 classificazione istopatologica delle immagini è stata progettata per ottimizzare le prestazioni del modello e garantire una buona generalizzazione. La pipeline iniziò con una pre-elaborazione completa, che includeva la riduzione delle immagini originali da 768 × 768 pixel a 224 × 224 pixel per soddisfare i requisiti di input dell'architettura EfficientNetB7. Le intensità dei pixel sono state normalizzate all'intervallo [0, 1] per stabilizzare e accelerare l'allenamento, e sono state applicate tecniche di aumento dei dati, incluse rotazioni casuali, spostamenti, zoom e flip orizzontali -- per introdurre variabilità e mitigare l'overfitting. Queste operazioni di pre-elaborazione garantivano che il modello apprese le caratteristiche invarianti e si generalizzasse bene a nuovi dati.
Il modello EfficientNetB7, pre-addestrato sul dataset ImageNet, è stato utilizzato come base, fornendo un robusto set di rilevatori di caratteristiche in grado di identificare sottili schemi istologici. Il modello fu successivamente affinato con due strati densi (128 e 64 unità) attivati dalla funzione Rectified Linear Unit (ReLU), seguiti da uno strato di output SoftMax per la classificazione in cinque classi. Questa architettura era ottimizzata sia per la precisione che per l'efficienza computazionale, rendendola particolarmente adatta all'analisi delle immagini mediche. La regola di aggiornamento dell'ottimizzatore di Adam per θ è mostrata nell'Equazione 7.
(7)
L'addestramento del modello è stato eseguito in lotti di 128 immagini utilizzando l'ottimizzatore Adam, che ha scalato in modo adattivo il tasso di apprendimento per promuovere una convergenza efficace. Un criterio di arresto precoce è stato implementato per monitorare la perdita di validazione e interrompere l'addestramento in caso di non miglioramento, prevenendo così l'overfitting e ottimizzando le risorse computazionali. Durante l'addestramento, le metriche chiave di prestazione — tra cui accuratezza, perdita, precisione e richiamo — sono state monitorate utilizzando un set di validazione di riserva che comprende il 20% dei dati. Il gruppo di riserva ha fornito approfondimenti fondamentali sulle capacità di generalizzazione del modello e ha mitigato il rischio di sovraadattamento a dati non visti.
L'integrazione di strategie avanzate di preprocessing, aumento dei dati e ottimizzazione ha portato a un'eccellente accuratezza e robustezza del modello, dimostrando un forte potenziale per l'applicazione clinica nella diagnosi del cancro ai polmoni e al colon. Tutti gli esperimenti sono stati condotti su una GPU NVIDIA Tesla P100 con 16 GB di memoria e 64 GB di RAM di sistema. L'intera pipeline di addestramento — inclusi preprocessing, augmentation e ottimizzazione del modello — richiedeva circa tre ore per completare dieci epoche con un batch di 128 minuti. Questa infrastruttura computazionale è stata scelta per gestire in modo efficiente il grande dataset aumentato, mantenendo al contempo tempi di addestramento ragionevoli e adeguati agli ambienti di ricerca clinica.
Per la riproducibilità, tutti i parametri critici sono stati esplicitamente specificati. Il modello è stato addestrato utilizzando l'ottimizzatore Adam con un tasso di apprendimento costante di 0,001. L'addestramento veniva svolto per un massimo di 10 epoche, con arresti precoci basati sulla perdita di validazione e una pazienza di tre epoche per prevenire il sovrafitting. Un lotto di 128 è stato utilizzato sia per la fase di addestramento che per quella di validazione. Per garantire risultati coerenti, il seed casuale è stato impostato a 42 in tutte le fasi del caricamento dei dati e dell'addestramento del modello. Durante l'addestramento non è stato applicato alcun programma di decadimento del tasso di apprendimento.
Il protocollo combinava preprocessing sistematico, iperparametri ben calibrati, schemi di validazione rigorosi e un ambiente computazionale aperto per promuovere sia accuratezza che riproducibilità. Il flusso di lavoro, integrando un'architettura robusta di deep learning, un'ottimizzazione efficace e un disegno sperimentale riproducibile, ha stabilito una solida base per la classificazione istopatologica delle immagini. Oltre a dimostrare forti prestazioni nella rilevazione del cancro ai polmoni e al colon, il protocollo ha definito un quadro scalabile che potrebbe essere esteso a sfide simili di imaging biomedico in future ricerche.
Durante lo sviluppo del modello, sono state implementate diverse modifiche pratiche per ottenere risultati affidabili e riproducibili. Sono stati utilizzati arresti precoci e aumenti dati avanzati per affrontare lo squilibrio di classe e l'overfitting. L'elevato consumo di memoria durante l'addestramento è stato ridimensionato le immagini a 224 × 224 pixel e utilizzando una dimensione batch di 128. Quando si incontrava una convergenza lenta, la velocità di apprendimento veniva empiricamente impostata a 0,001. La verifica manuale della struttura delle directory delle immagini veniva effettuata per evitare errori di caricamento dei dati, e venivano corretti semi casuali per garantirne la riproducibilità. Queste strategie possono servire come guida pratica per i ricercatori che affrontano sfide analoghe nell'addestramento di modelli di deep learning con grandi dataset di immagini istopatologiche.
Analisi statistica
La modellazione statistica delle prestazioni del modello è stata condotta utilizzando un insieme completo di metriche per valutare accuratezza, stabilità e generalizzabilità. Le principali misure di valutazione includevano precisione, richiamo, punteggio F1 e accuratezza per quantificare la performance nella classificazione. La matrice di confusione è stata utilizzata per riportare metriche di errore specifiche per classe, calcolate secondo le Equazioni 8-10:
(8)
(9)
(10)
Inoltre, sono state valutate misure di errore come l'Errore Quadratico Medio (MSE), l'Errore Quadratico Medio Radice (RMSE) e l'Errore Assoluto Medio (MAE) per valutare l'accuratezza della previsione, calcolata utilizzando le Equazioni 11-13:
(11)
(12)
(13)
La curva Receiver Operating Characteristic (ROC) e l'Area Under the Curve (AUC) sono state impiegate per valutare la capacità del modello di discriminare tra classi a varie soglianze. L'AUC è stato calcolato come mostrato nell'Equazione 14:
(14)
Queste analisi statistiche hanno fornito una valutazione approfondita delle prestazioni di classificazione, della robustezza e della generalizzazione dei dati mai visti del modello.
Studio dell'ablazione
Per valutare ulteriormente il contributo di ciascun componente all'interno del quadro di classificazione, è stato condotto uno studio di ablazione omettendo o modificando selettivamente specifici elementi architettonici del modello finale. In particolare, sono state valutate due configurazioni sperimentali alternative: una includeva sia gli strati Global Average Pooling che quelli Flatten, seguita da due strati densi (rispettivamente 128 e 64 unità); l'altra configurazione ometteva sia lo strato Flatten che lo strato denso di 128 unità, lasciando solo lo strato di Global Average Pooling e un unico strato denso di 64 unità.
L'architettura completa, come riportato nei risultati principali, ha raggiunto un'accuratezza di validazione di circa il 96%. Quando il livello Flatten è stato escluso e sono stati utilizzati solo gli strati densi di 128 e 64 unità dopo il Global Average Pooling, l'accuratezza della validazione è scesa all'81,6%. Al contrario, quando è stato escluso lo strato denso di 128 unità e lo strato Flatten è stato mantenuto, l'accuratezza della validazione è diminuita all'88%. Questi risultati indicarono la necessità di includere sia lo strato extra denso sia i componenti architettonici appropriati per ottenere prestazioni ottimali di classificazione. I risultati degli esperimenti di ablazione hanno dimostrato costantemente che la combinazione di Global Average Pooling, Appiattimento e due strati completamente connessi permetteva una rappresentazione delle caratteristiche più espressiva, risultando in una classificazione istopatologica delle immagini più robusta e accurata.