Les érosions osseuses sont une caractéristique pathologique importante de la polyarthrite rhumatoïde. L’objectif de ce travail est de mettre en place un outil de formation pour fournir aux utilisateurs des conseils sur l’identification des cassures corticales pathologiques sur des images de tomodensitométrie quantitative périphérique à haute résolution pour l’analyse de l’érosion.
Les érosions osseuses sont une caractéristique pathologique de plusieurs formes d’arthrite inflammatoire, y compris la polyarthrite rhumatoïde (PR). L’augmentation de la présence et de l’ampleur des érosions est associée à de mauvais résultats, à une mauvaise fonction articulaire et à une mauvaise progression de la maladie. La tomodensitométrie quantitative périphérique à haute résolution (HR-pQCT) offre une visualisation in vivo inégalée des érosions osseuses. Cependant, à cette résolution, des discontinuités dans l’enveloppe corticale (cassures corticales) associées à des processus physiologiques normaux et à une pathologie sont également visibles. L’étude grouP for xtrEme Computed Tomography in Rheumatoid Arthritis a précédemment utilisé un processus consensuel pour développer une définition de l’érosion pathologique dans la HR-pQCT : une cassure corticale détectée dans au moins deux coupes consécutives, dans au moins deux plans perpendiculaires, de forme non linéaire, avec une perte osseuse trabéculaire sous-jacente. Cependant, malgré l’existence d’une définition consensuelle, l’identification de l’érosion est une tâche exigeante qui pose des problèmes de variabilité inter-évaluateurs. L’objectif de ce travail est de mettre en place un outil de formation pour fournir aux utilisateurs des conseils sur l’identification des cassures corticales pathologiques sur les images HR-pQCT pour l’analyse de l’érosion. Le protocole présenté ici utilise un module construit sur mesure (Bone Analysis Module (BAM) - Training), implémenté en tant qu’extension d’un logiciel de traitement d’images open-source (3D Slicer). À l’aide de ce module, les utilisateurs peuvent s’entraîner à identifier les érosions et comparer leurs résultats aux érosions annotées par des rhumatologues experts.
Les érosions osseuses se produisent lorsque l’inflammation provoque une perte osseuse localisée à la surface de l’os cortical. Ces érosions s’étendent dans la région osseuse trabéculaire sous-jacente. Ils sont une caractéristique pathologique de plusieurs formes d’arthrite inflammatoire, y compris la polyarthrite rhumatoïde (PR)1. La présence et l’ampleur de l’érosion sont associées à de mauvais résultats, à la fonction du patient et à la progression de la maladie 2,3,4,5. Alors que la radiographie simple reste la norme clinique pour l’évaluation de l’érosion, la tomodensitométrie quantitative périphérique à haute résolution (HR-pQCT) fournit des images 3D et une sensibilité et une spécificité supérieures pour la détection de l’érosion 6,7. Pour l’arthrite inflammatoire, telle que la PR, la HR-pQCT est généralement effectuée sur les 2e et 3e articulations métacarpo-phalangiennes - les articulations les plus touchées de la main8. Étant donné que les images HR-pQCT ont une résolution spatiale élevée, des interruptions physiologiques de la surface corticale sont également observées chez les personnes en bonne santé sans RA9. Ces interruptions corticales sont souvent associées à des canaux vasculaires ou à un foramen nutritif traversant l’os10. Ainsi, le défi consiste à distinguer les interruptions corticales associées à un processus pathologique (c’est-à-dire les érosions pathologiques) des caractéristiques non pathologiques.
La définition consensuelle d’une érosion osseuse pathologique a été publiée par l’étude GROUP for xtrEme Computed Tomography in Rheumatoid Arthritis (SPECTRA) comme la présence d’une interruption définie dans la couche corticale de l’os qui s’étend sur au moins deux coupes consécutives et est détectable dans deux plans perpendiculaires ou plus11. De plus, l’interruption doit être de forme non linéaire et s’accompagner d’une perte dans la région trabéculaire. Des exemples visuels d’interruptions corticales qui répondent ou non aux critères d’érosion sont présentés dans Klose-Jensen et al.12.
Cependant, toutes les interruptions corticales qui répondent aux critères ci-dessus ne sont pas classées comme érosions. Les interruptions sont parfois causées par des processus physiologiques tels que les canaux vasculaires (Figure 1). Ceux-ci peuvent être identifiés et différenciés des érosions en raison de leurs emplacements anatomiques prévisibles, de leurs marges parallèles et droites et de leur taille submillimétrique13. Les kystes sont une autre forme d’interruption corticale qui n’est pas considérée comme une érosion. Ils ont souvent une structure trabéculaire arrondie avec une paroi kystique claire 13. En contraste avec les arêtes vives et la structure trabéculaire ouverte montrée par les érosions. Cependant, il est possible que des érosions se forment dans les sites kystiques, ce qui rend ambiguë la délimitation du volume de perte osseuse causée par les érosions et non par les kystes. Bien que le but de cette étude ne soit pas de résoudre cette ambiguïté avec d’autres critères, il est nécessaire de fournir des exemples complets d’érosion pathologique et d’interruptions corticales physiologiques.
Figure 1 : Exemple d’interruptions corticales qui n’ont pas été causées uniquement par des érosions. (A) Un dessin illustrant l’emplacement commun des canaux vasculaires à la base de la tête métacarpienne. Exemples de canaux vasculaires dans les plans (B) coronal, (C) sagittal, (D) et (E) axial. (F) Exemple d’une interruption corticale causée par un kyste. (G) Exemple d’un volume vide dans la région trabéculaire de l’os impliquant à la fois des kystes et des érosions. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Malgré les défis liés à l’identification de l’érosion, il n’existe actuellement aucun outil de formation en place pour fournir aux utilisateurs moins expérimentés des conseils sur l’interprétation des images HR-pQCT pour l’analyse de l’érosion. Récemment, un module open source pour l’analyse de l’érosion appelé module d’analyse osseuse (BAM) - Erosion Volume a été développé, implémenté en tant qu’extension d’un logiciel de traitement d’images open source pour permettre la visualisation de l’érosion et les analyses volumétriques14. Le protocole présenté ici décrit l’utilisation d’un module de formation ajouté à BAM (BAM - Training), qui compare les tentatives d’identification de l’érosion d’un utilisateur en comparant l’identification de l’érosion avec les érosions annotées par des rhumatologues experts. Cet outil de formation fournit aux utilisateurs une rétroaction sur l’identification de l’érosion afin de guider les améliorations dans l’analyse de l’érosion. Les instructions d’installation du logiciel sont fournies à l’étape 1. Pour l’acquisition de nouvelles données, reportez-vous aux étapes 3 à 5.3. Pour l’utilisation du module de formation uniquement, reportez-vous à l’étape 2.
Toutes les méthodes de ce protocole suivent les lignes directrices établies par le Conjoint Health Research Ethics Board de l’Université de Calgary (REB19-0387).
1. Installez 3D Slicer 15 et les modules d’analyse osseuse
Figure 2 : Exemple de fenêtre de paramètres après l’ajout de modules d’analyse de structure à une installation de 3D Slicer. L’image montre une capture d’écran de la fenêtre des paramètres avec les modules mis en surbrillance dans la zone rouge. Veuillez cliquer ici pour voir une version agrandie de cette figure.
2. Module de formation
Figure 3 : menu déroulant de 3D Slicer. Le menu déroulant permet de trouver les modules d’analyse osseuse et de sélectionner le module d’entraînement. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Figure 4 : Identification du volume à l’intérieur de la surface périostée de l’os. (A) Exemple d’un masque. Le masque est visualisé sous la forme d’une image binaire. (B) Exemple de segmentation. La segmentation fait référence à la visualisation de l’image binaire superposée à l’image en niveaux de gris. Ces distinctions sont faites par 3D Slicer. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Figure 5 : Exemple de capture d’écran d’un module d’entraînement dans 3D Slicer. (A) Cliquez pour ajouter de nouveaux points d’amorçage. (B) Cliquez pour calculer les volumes d’érosion. (C) Cliquez pour importer des images. (D) Cliquez pour afficher les points de départ placés par les experts. Veuillez cliquer ici pour voir une version agrandie de cette figure.
3. Acquisition et exportation d’images à utiliser dans l’outil d’analyse de l’érosion
4. Conversion de fichiers et génération de masques osseux
REMARQUE : Selon le format d’image, suivez l’étape 4.1 pour les images AIM (format d’image propriétaire HR-pQCT), MHA (format ITK MetaImage), nii (NIfTI - Neuroimaging Informatics Technology Initiative), NRRD (Nearly Raw Raster Data) ou l’étape 4.2 pour les images DICOM.
5. Identification des érosions
6. Statistiques sur l’érosion
À l’aide de l’outil de formation, les utilisateurs peuvent s’entraîner à identifier les sites d’érosion tout en recevant une rétroaction sur leurs résultats. Cette boucle de rétroaction peut améliorer la capacité de l’utilisateur à identifier les érosions et éventuellement utiliser les modules BAM pour identifier les érosions sur ses propres images. Le retour d’information après le placement du point d’amorçage est basé sur les critères suivants. 1) Si le nombre de points d’amorçage placés ne correspond pas au nombre d’érosions de référence, l’utilisateur est alors invité à supprimer ou à ajouter le nombre approprié de points d’amorçage. 2) S’il n’est pas possible d’apparier l’emplacement du point d’amorçage à une érosion de référence, une rétroaction indiquant qu’il n’y a pas d’érosion à l’emplacement de ce point d’amorçage s’affiche. 3) Si un point d’ancrage est associé à une interruption corticale pathologique/physiologique de référence telle qu’un kyste ou un canal vasculaire, l’utilisateur est informé du type d’interruption corticale qu’il a tenté d’identifier comme une érosion et est invité à retirer le point d’amorçage. 4) Si l’emplacement du point d’amorçage chevauche une érosion de référence, l’algorithme peut toujours ne pas détecter l’érosion. Cela peut se produire lorsque le point d’amorçage n’a pas été centré dans l’érosion. Dans ce cas, l’utilisateur est invité à ajuster la position du point d’amorçage. 5) Si un point d’amorçage est placé trop loin de toute érosion, l’utilisateur est informé de son mauvais placement et encouragé à réessayer. 6) Lorsqu’un emplacement de point d’amorçage correspond à l’érosion de référence, une invite s’affiche pour informer l’utilisateur de sa tentative réussie d’identifier l’érosion à ce point d’amorçage spécifique.
La section suivante illustre des exemples de fonctionnement du module en fonction de différentes entrées. Les entrées correctes et incorrectes seront démontrées dans les exemples suivants. La figure 6A montre l’emplacement du point d’amorçage qui se trouve à l’intérieur de l’érosion. Il n’existe qu’une seule érosion dans cette image, donc le calcul des érosions avec le point de départ donnera les résultats attendus. La figure 6B montre l’invite qui s’affiche aux utilisateurs lorsque leur tentative d’identification des érosions correspond à l’image annotée par des experts. Le module affiche également les résultats sous forme de segmentations sur l’image en niveaux de gris (Figure 6C). Si l’utilisateur plaçait un point d’amorçage à un emplacement sans érosion, tel que la figure 7A, le module affichait une invite d’erreur (Figure 7B) indiquant qu’il n’y avait pas d’érosion à cet emplacement et suggérant à l’utilisateur de déplacer/supprimer les points d’amorçage.
Figure 6 : Exemple d’identification correcte de l’érosion. (A) Exemple d’un utilisateur qui place correctement un point d’amorçage dans le site d’érosion. (B) Exemple d’invite de rétroaction lorsque toutes les érosions ont été identifiées correctement. (C) Exemple de segmentation de l’érosion affichée lorsqu’une érosion est calculée correctement. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Figure 7 : Exemple d’identification incorrecte de l’érosion. (A) Exemple d’un point d’amorçage placé là où il n’y a pas d’érosion. (B) Exemple d’une invite d’erreur lorsqu’un point d’amorçage est placé à un endroit où il n’y a pas d’érosion. Veuillez cliquer ici pour voir une version agrandie de cette figure.
L’emplacement de tous les kystes et canaux vasculaires sur les images d’entraînement fournies a été identifié par un expert. Par conséquent, il est possible de détecter lorsqu’un utilisateur tente d’identifier de manière incorrecte un kyste ou un canal vasculaire. La figure 8A illustre une tentative d’identification d’un kyste en y plaçant un point de pérage. La figure 8B est l’invite d’erreur qui s’affiche par la suite.
Figure 8 : Exemple d’identification d’un kyste. (A) Exemple d’un point d’amorçage placé sur un kyste. (B) Exemple d’une invite d’erreur lorsqu’un point d’amorçage est placé sur un kyste. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Le module informera également l’utilisateur s’il dispose du bon nombre de points d’amorçage. Si l’utilisateur a placé un nombre incorrect de points d’amorçage, le module l’informerait du nombre exact de points d’amorçage manquants ou supplémentaires pour identifier toutes les érosions sur l’image. Le module donne également un retour d’information pour chaque point d’amorçage placé. Par conséquent, l’utilisateur sait quelles actions effectuer pour chaque point d’amorçage individuel. La figure 9 illustre un exemple où un utilisateur n’a placé qu’un seul point d’amorçage alors que deux étaient attendus.
Figure 9 : Exemple d’érosions calculées en l’absence d’un point d’amorçage. L’exemple montre un exemple où l’utilisateur n’a placé qu’un seul point d’amorçage alors que deux étaient attendus. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Si un utilisateur a de la difficulté à trouver une partie ou la totalité des érosions, il a la possibilité de révéler les emplacements annotés par des experts en appuyant sur un bouton intitulé Révéler les points de départ corrects (Figure 5D). Lorsqu’il est enfoncé, ce bouton charge les points d’amorçage corrects dans la fenêtre actuelle du segment 3D.
En résumé, cela montre que le module logiciel peut évaluer la justesse de la tentative de l’utilisateur d’identifier les érosions dans les images sélectionnées en comparant l’érosion calculée avec les érosions annotées par des experts. De plus, le module fournit un retour d’information basé sur chaque point d’amorçage placé par l’utilisateur pour le guider vers l’emplacement du point d’amorçage attendu et les paramètres d’entrée.
ID de numérisation | Interruption corticale | Os | Étiquette | Emplacement du centroïde | Volume (mm3) | Superficie (mm2) | Rondeur | Nombre de voxels (voxels) |
3_Training.nii | Érosion | Métacarpien | SEEDS_28-1 | 210, 108, 242 | 3.321668853 | 14.46818378 | 0.74411491 | 14853 |
3_Training.nii | Érosion | Métacarpien | SEEDS_28-3 | 179, 100, 241 | 1.100739562 | 7.121231239 | 0.7239659 | 4922 |
Tableau 1 : Exemple d’un fichier de sortie généré (format csv) décrivant les érosions calculées et leurs statistiques.
Cet outil de formation permet d’apprendre à identifier les érosions à l’aide du module d’analyse osseuse. L’utilisation ultérieure de cet outil d’analyse de l’érosion au-delà de l’entraînement nécessite l’accès à des images de bonne qualité, avec peu ou pas d’artefact de mouvement. La définition de l’érosion HR-pQCT basée sur la littérature décrit les caractéristiques anatomiques associées aux érosions pathologiques qui peuvent être rapportées avec une reproductibilité raisonnable11,20. Cependant, cette définition ne tient pas compte des localisations anatomiques communes des canaux vasculaires, ce qui peut entraîner leur classification erronée en tant qu’érosions osseuses10.
Les étapes critiques de ce protocole sont la génération du masque osseux, la mise en place des points de départ et la génération du volume d’érosion. Bien que des méthodes automatisées soient mises en œuvre pour générer les masques et le volume d’érosion, les masques nécessitent souvent une correction manuelle pour garantir des résultats satisfaisants. Une description complète des outils disponibles pour effectuer les corrections manuelles est fournie. Le placement des points d’amorçage est guidé par les exemples d’entraînement fournis par le module BAM-Training.
Sur la base des données utilisées à ce jour, ce protocole fournit des suggestions de dépannage lorsque le module d’analyse de l’érosion ne produit pas les résultats escomptés. Dans les travaux futurs, l’accès à des données d’entraînement supplémentaires sera fourni. Une étude antérieure a montré que les volumes d’érosion évalués avec cette méthode sont comparables aux méthodes existantes 14,21,22. La fourniture de données d’apprentissage permettra de les comparer aux nouveaux outils d’analyse de l’érosion au fur et à mesure de leur élaboration23.
L’outil de formation présenté ici aide principalement à l’identification de l’érosion ; Cependant, la méthode est actuellement limitée par l’absence de consensus sur la définition de l’étendue d’une érosion de l’os trabéculaire. Néanmoins, les modules BAM sont open source, de sorte que, à mesure que les définitions futures de l’étendue de l’érosion changent, d’autres chercheurs ont la possibilité de modifier les modules pour répondre à leurs besoins.
Au fur et à mesure que l’utilisation de la HR-pQCT dans la recherche rhumatologique se développe, l’outil de formation fournit aux utilisateurs inexpérimentés des conseils sur l’identification des interruptions corticales pathologiques sur les images HR-pQCT pour l’analyse de l’érosion. Cet outil s’appliquera aux chercheurs, quelle que soit la méthode choisie pour l’analyse de l’érosion. Bien qu’une identification entièrement automatisée de l’érosion soit souhaitable pour améliorer la reproductibilité et la rapidité de l’analyse, de grands ensembles de données de référence/référence avec des annotations précises sont nécessaires pour entraîner les modèles d’apprentissage automatique. En tant qu’outil open source, ce module offre la possibilité de développer collectivement de grands ensembles de données annotées pour une utilisation future dans l’apprentissage automatique. L’utilisation de cet outil de formation permettra à un plus grand nombre de chercheurs d’inclure l’analyse de l’érosion dans leurs recherches sur les RH-pQCT.
Tous les auteurs n’ont aucun conflit d’intérêts à signaler.
Les auteurs tiennent à remercier les organismes de financement suivants qui ont appuyé ce travail. La SLM est financée par la Société de l’arthrite (STAR-18-0189) et la subvention de planification et de diffusion des Instituts de recherche en santé du Canada. JJT est titulaire d’une bourse de recherche des IRSC.
Name | Company | Catalog Number | Comments |
3DSlicer | Open Source | N/A | Download at https://www.slicer.org/ |
BAM Erosion Analysis Modules | Open Source | N/A | Version used in manuscript: download at https://doi.org/10.5281/zenodo.7943007 |
XtremeCTII | Scanco Medical | N/A |
Demande d’autorisation pour utiliser le texte ou les figures de cet article JoVE
Demande d’autorisationExplorer plus d’articles
This article has been published
Video Coming Soon