$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
L’ensemble de données comprend un total de 1 190 tranches de scanner CT provenant de 110 cas, classées en trois catégories : normales (55 cas), bénignes (15 cas) et malignes (40 cas). Chaque cas comprend plusieurs tranches CT (allant d’environ 80 à 200 tranches par cas), offrant des vues axiales variées de la région thoracique. Les images CT ont été obtenues en DICOM à l’aide du scanner SOMATOM avec des paramètres d’imagerie standards : tension du tube = 120 kV, épaisseur de tranche = 1 mm, largeur de fenêtre = 350–1 200 unités de Hounsfield (HU), et valeurs du centre de la fenêtre = 50–600 HU, lors d’une suspension complète de l’inspiration.
Toutes les images ont été entièrement désidentifiées avant analyse afin de supprimer toute information personnellement identifiable (PII). L’ensemble de données a été approuvé sur le plan éthique par les comités d’examen institutionnels des centres médicaux participants, avec un consentement écrit renoncé par le comité de surveillance. Les cas concernaient des personnes issues de milieux variés, tels que des fonctionnaires gouvernementaux, des agriculteurs et des résidents de plusieurs provinces irakiennes (dont Bagdad, Wasit, Diyala, Salahuddin et Babylone) avec une diversité de genre, d’âge, de niveau d’éducation et de statut de vie.
Comme l’ensemble de données est accessible et désidentifié, aucune approbation éthique supplémentaire n’a été nécessaire pour son utilisation dans cette étude. Le jeu de données respecte les termes et conditions spécifiés par ses contributeurs originaux et la plateforme hôte.
La méthodologie de cette recherche utilise un processus en plusieurs étapes conçu pour créer un modèle de prédiction à l’aide du jeu de données IQ-OTH/NCCD sur le cancer dupoumon 21. Cette méthodologie établit un terrain solide pour les dispositifs de santé centrés sur le consommateur, où une latence réduite est requise, ce qui peut offrir une meilleure précision. La figure 2 illustre le mécanisme de fonctionnement du modèle proposé.

Figure 2 : Diagramme de flux de travail du modèle proposé illustrant les étapes de prétraitement, d’augmentation, de classification par Transformateur de Vision et d’évaluation. Veuillez cliquer ici pour voir une version agrandie de cette figure.
1. Description du jeu de données
L’ensemble de données IQ-OTH/NCCD sur le cancer du poumon a été collecté à l’Irak-Oncology Teaching Hospital/National Center for Cancer Diseases à l’automne 2019. La partie entraînement du jeu de données utilisée contenait 1 097 images dont la description est donnée dans le tableau 2, qui illustre le nombre d’exemples d’instances de différentes classes.
| Classe | Nombre d’images |
| Normal | 416 |
| Bénin | 120 |
| Malin | 561 |
Tableau 2 : Exemples d’images CT normales, bénignes et malignes issues de l’ensemble de données.
L’ensemble de données IQ-OTH/NCCD sur le cancer du poumon a été choisi pour sa représentation complète des scanners CT normaux, bénins et malins. Bien que d’autres ensembles de données, tels que LIDC-IDRI et NSCLC-Radiomics, soient disponibles, ils se concentrent principalement sur la détection des nodules ou manquent d’échantillons suffisants de cas bénins. L’ensemble de données IQ-OTH/NCCD est particulièrement adapté à notre étude car il permet au Transformateur de Vision d’apprendre les caractéristiques discriminantes dans les trois classes, permettant une classification robuste des motifs subtils dans les images de scanner pulmonaire.
2. Prétraitement des données
Le prétraitement est une étape importante pour améliorer la performance du modèle grâce à la cohérence et aux ajustements appropriés des données qui seront traitées via le réseau de neurones. Afin d’assurer la cohérence de la taille d’entrée du réseau de neurones, nous avons mis toutes les images à la même dimension de 256 x 256 pixels, et normalisé les données à une valeur de pixel comprise entre 0 et 1 dans l’espoir d’améliorer l’entraînement et la convergence des modèles. Le tableau 3 contient les valeurs des techniques d’augmentation.
| Technique | Valeur |
| Normalisation | - |
| Redimensionnement | image_size x image_size |
| RotationRandom | facteur = 0,02 |
| RandomZoom | height_factor=0,2, width_factor=0,2 |
Tableau 3 : Techniques d’augmentation appliquées avec plages de paramètres.
Pour renforcer la robustesse du modèle face au surapprentissage et améliorer sa capacité de généralisation, des techniques d’augmentation de données telles que des rotations aléatoires et des zooms sont appliquées, simulant divers angles et tailles de manifestations du cancer du poumon telles qu’elles pourraient apparaître chez différents patients. Des rotations aléatoires avec des angles échantillonnés uniformément à partir de 0,02 radian et des transformations de zoom aléatoires avec des facteurs de hauteur et de largeur variables ont été appliquées dans cette étude. Les images suivant l’augmentation sont affichées à la Figure 3.

Figure 3 : Images CT après augmentation démontrant la rotation, le zoom et la normalisation pour améliorer la généralisation du modèle. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Ces techniques de prétraitement sont nécessaires, et l’augmentation a été utilisée dans toutes les classes pour garantir la robustesse du modèle.
3. Architecture du modèle
Adaptant le nouveau cadre Vision Transformer (ViT) pour traiter efficacement des données d’images complexes, l’architecture du modèle vise à catégoriser les scanners CT en trois catégories : normales, bénignes et malignes. Avec cette méthode, des images d’entrée de 256 x 256 pixels sont traitées. Pour garantir la cohérence et favoriser un apprentissage plus efficace des modèles, chaque image est soumise à plusieurs opérations de prétraitement, telles que le redimensionnement et la normalisation. Pour améliorer la résilience du modèle aux changements d’échelle et d’orientation de l’image, la conception commence par une couche d’augmentation de données utilisant des méthodes telles que la normalisation, le redimensionnement, les rotations aléatoires de 0,02 radian, et des zooms aléatoires allant jusqu’à 20 %. Après augmentation, le modèle prend 16 patchs de 16 pixels sur chaque image, donc chaque image en a 256.
L’algorithme 1 définit le flux de travail du modèle proposé et sa construction, ainsi que l’ajustement fin effectué par rapport au modèle.
Algorithme 1 : Classification du cancer du poumon à l’aide de transformateurs de vision
Entrée : Ensemble d’images CT I issues de l’ensemble de données IQ-OTH/NCCD
Résultat : Résultats de classification C en catégories : Normal, bénin, Malin
Prétraitement:
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
Configuration du modèle :
Initialize Vision Transformer (ViT) Model
Set number of patches P = 16 x 16
Set number of heads H = 6 in multi-head attention
Formation:
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
Validation:
Split data into TrainingSet (80%) and ValidationSet (20%)
ComputeValidationMetrics(ValidationSet)
Évaluation:
C <- Classify(I)
ComputePerformanceMetrics(C)
Return C
Pour préserver les relations spatiales entre les patchs, ces patches sont aplatis (équation 1) en vecteurs de 768 dimensions (car chaque patch a 16 x 16 x 3 = 768) puis envoyés à travers une couche d’encodage de patch, qui projette chaque vecteur dans un espace de 64 dimensions en intégrant des embeddings positionnels. Le cœur de l’architecture comprend huit couches transformers, chacune dotée d’un mécanisme d’attention multi-têtes avec six têtes, permettant au modèle de se concentrer simultanément sur différentes parties de l’image et de capturer une large gamme de caractéristiques. Elle est représentée par les équations 2, 3 et 4.

Où μ est la moyenne et σ2 est la variance de l’entrée x, et ε est une petite constante pour éviter la division par zéro.

Où Q est la matrice de requête, K est la matrice clé, V est la matrice de valeurs, et dk est la dimension des vecteurs clés.


Où WiQ, WiK et WiV sont respectivement les matrices de poids apprises pour les requêtes, les clés et les valeurs, et WO est la matrice de poids de sortie.
Le schéma bloc du modèle proposé a été illustré à la Figure 4.

Figure 4 : Schéma bloc du modèle Vision Transformer avec tête MLP, détaillant les principales couches de traitement et l’architecture. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Chaque couche transformateur inclut une perception multi-couches (MLP) avec des unités cachées d’abord mises à l’échelle à 128 puis à 64, élargissant et comprimant efficacement le flux d’information pour capturer des dépendances complexes.
Le dropout est appliqué stratégiquement à un taux de 0,1 dans les mécanismes d’attention et les MLP pour éviter le sur-ajustement. La sortie de l’encodeur transformateur est traitée via une tête MLP composée de deux couches denses de 2 048 et 1 024 unités respectivement, utilisant des unités linéaires d’erreur gaussiennes (GELU) pour l’activation, ce qui a démontré de bonnes performances dans les applications d’apprentissage profond. Cela est réalisé à l’aide de l’équation 5.

Là où W1 et W2 sont des matrices de poids, b1 et b2 sont des biais, et GELU est la fonction d’activation utilisée.
Un coupure de 0,1 était utilisé dans les couches transformateurs pour éviter le surapprentissage sans perdre des informations importantes sur les caractéristiques. La fonction d’activation GELU a été utilisée pour ses caractéristiques non linéaires et lisses, qui favorisent l’écoulement du gradient et la convergence dans les modèles basés sur transformateurs. Le dropout est régularisé à l’aide de l’équation 6.
Où p est le taux d’abandon (par exemple, 0,1 ou 0,5), et la couche d’abandon règle aléatoirement une fraction p des unités d’entrée à zéro pendant l’entraînement.
Un taux de dégagement de 0,5 dans ces couches aide encore à atténuer le sur-ajustement. La couche finale du modèle est une couche logits (équation 7) qui produit trois logits de classes correspondant aux catégories normale, bénigne et maligne, en utilisant une fonction de perte croisée par entropie croisée catégorique clairsemée (équation 8) appropriée à ce contexte de classification multi-classes.


Où zi est le vecteur logits pour la classe i, SoftMax(zi)) représente les probabilités prédites, et yi est la véritable étiquette de classe.
Le modèle est compilé avec l’optimiseur AdamW (équations 9, 10, 11, 12 et 13), une extension de l’optimiseur Adam qui gère plus efficacement la décroissance de poids, réglée à un taux d’apprentissage de 0,001 et une décroissance de poids de 0,0001, optimisant à la fois la vitesse d’apprentissage et la stabilité du modèle.





Où mt et vt sont les premiers et seconds moments des gradients,
et
sont les moments corrigés par biais, η est le taux d’apprentissage, et ε est une petite constante pour éviter la division par zéro.
Le modèle, une fois entraîné, utilise une taille de lot de 32 pour tirer parti de l’accélération GPU afin de calculer plus rapidement et est configuré pour s’entraîner pendant 100 époques.
Mais l’entraînement dispose d’un mécanisme d’arrêt précoce sur la perte de validation pour limiter l’entraînement lorsque le modèle cesse de s’améliorer, afin d’économiser des ressources de calcul et d’éviter le surentraînement de 5 époques consécutives. La performance du modèle est suivie avec des indicateurs tels que la précision catégorielle clairsemée et la précision du top 2, qui déterminent la capacité de classification du modèle dans les catégories prédites les plus probables. Les rappels pour l’entraînement du modèle comprennent des points de contrôle qui sauvegardent le modèle le plus performant selon la précision de la validation afin de garantir que la version la plus réussie du modèle soit conservée une fois le processus d’entraînement terminé. Cette architecture robuste et bien conçue exploite les capacités des transformateurs pour traiter les données d’imagerie médicale avec une approche très avancée pour utiliser des méthodes d’IA contemporaines dans des applications majeures en diagnostic de santé.
4. Formation et validation
Le modèle a été entraîné dans un environnement Kaggle utilisant un GPU NVIDIA P100, et une mini-descente de gradient batch avec une taille de 32 a été utilisée lors du processus d’entraînement. L’optimiseur utilisé pendant l’entraînement était AdamW avec un taux d’apprentissage de 0,001 et une décroissance de poids de 0,0001 ; C’était un bon équilibre entre une convergence rapide et une certaine régularisation. Le modèle a été entraîné pendant 100 époques, cependant, un arrêt précoce lors de la perte de validation avec une patience de 5 époques a été utilisé. Pour faire simple, si l’entraînement n’améliorait pas la perte de validation pendant 5 époques consécutives, l’entraînement s’arrêtait pour surajustement et efficacité computationnelle. Dans la plupart des cas, ce modèle restaurait les poids de l’époque avec la plus faible perte de validation, indiquant qu’il fonctionnait de manière optimale et n’avait pas besoin d’exécuter les 100 époques complètes. Environ 32 minutes pour entraîner le modèle.
Pendant la formation, les métriques comprenaient une précision catégorielle clairsemée et une précision top 2, surveillées à la fois sur les ensembles d’entraînement et de validation. Ces indicateurs fournissent un aperçu de la fréquence à laquelle un modèle prédit la bonne classe et si la bonne classe se situe parmi les deux premières prédictions, ce qui est important dans les applications médicales car les erreurs de classification à haute confiance peuvent avoir des conséquences importantes. Les courbes d’apprentissage pour la perte et la précision sont visualisées à la Figure 5.
.
Figure 5 : Précision de l’entraînement et de la validation ainsi que des courbes de perte sur 50 époques montrant un apprentissage stable et un surajustement minimal. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Aucune validation croisée n’a été réalisée dans la présente étude. Cette configuration permettait des expérimentations efficaces avec l’architecture du modèle, y compris des ajustements des couches de transformateurs et des tailles des têtes MLP, tout en garantissant que le modèle se généralisait bien sur des données de validation invisibles.
5. Analyse statistique
La performance du modèle Vision Transformer a été analysée statistiquement à partir du jeu de données IQ-OTH/NCCD sur le cancer du poumon. La précision, le rappel, le score F1 et la précision ont été calculés à l’aide des équations 14, 15, 16, 17 respectivement. Ces quatre mesures sont considérées comme les méthodes conventionnelles pour les tâches de classification, car elles peuvent révéler des informations sur la performance du modèle lors de la classification et de la classification par classe pour chaque classe.




La mémoire est une mesure de la capacité d’un modèle à identifier toutes les instances pertinentes dans une classe, tandis que la précision est la proportion d’identifications positives qui étaient effectivement correctes. Le score F1 équilibre rappel et précision lorsque les deux sont importants.
Les mesures de performance utilisées pour ce devoir étaient l’équation 18 du coefficient de corrélation de Matthews (MCC) et l’équation de Kappa 19 de Cohen.


Où P0 est l’accord observé tandis que Pe est l’accord attendu.
Le MCC est une mesure globale de la performance dans la classification, prenant en compte à la fois les vrais positifs et négatifs, les faux positifs et les faux négatifs. Sa valeur peut être comprise entre -1 et 1, 1 étant la prédiction parfaite, 0 la prédiction aléatoire, et -1 le désaccord total entre les étiquettes prédite et vraie. Le MCC était utile pour des comparaisons entre différentes performances de modèles lorsqu’il était appliqué à des ensembles de données déséquilibrés et fournissait une mesure équilibrée indépendamment des tailles de classes différentes.
Dans le cadre d’une analyse statistique supplémentaire, le score F2, qui priorise le rappel, a été calculé comme le montre l’équation 20.

La performance du modèle a également été quantifiée à l’aide de l’erreur quadratique moyenne (MSE), de l’erreur quadratique moyenne (RMSE) et de l’erreur absolue moyenne (MAE) – qui sont généralement des mesures pour les tâches de régression mais modifiées ici pour la tâche de classification afin de quantifier la taille moyenne de l’erreur sans tenir compte de la direction.
Afin de déterminer s’il serait pratique d’intégrer les ViT dans des dispositifs de santé destinés aux consommateurs, nous avons réalisé des résultats d’évaluation de performance approfondis en nous concentrant uniquement sur l’efficacité temporelle du modèle. Nous avons créé des fonctions pour rapporter le temps nécessaire pour prédire une seule ou simplement plusieurs images, car cela devient particulièrement pertinent pour les applications en temps réel. Pour chaque image, avant de commencer à prédire, les données sont d’abord prétraitées pour avoir la forme de l’image d’entrée désirée par le modèle. Nous avons enregistré le moment où la prédiction a commencé et le moment où la prédiction a été complétée, afin d’obtenir les résultats du temps et de la latence. Le temps et la latence moyenne ont été calculés à l’aide respectivement de l’équation 21 et de l’équation 22.


Où:
- N est le nombre d’images (échantillons).
- tend est le temps enregistré après avoir prédit la i-ème image.
- Tstart est le temps enregistré avant de prédire la i-ème image.
D’autres expériences ont été réalisées pour évaluer la robustesse du modèle Vision Transformer que nous avons proposé en introduisant systématiquement du bruit et du flou supplémentaires dans les images. Le bruit gaussien était ajouté à chaque pixel avec un facteur de bruit de 0,4 tout en découpant des valeurs de pixels dans la plage de [0,1]. Le facteur de flou a été diminué via le flou gaussienne avec une taille de grain de 45× 45. Ces expériences visent à évaluer le modèle de manière exhaustive sous des dégradations simulées de la qualité de l’image perceptive.