Vidéoradiographie biplanaire pour étudier le poignet et les articulations radio-ulnaires distales

Résumé

La vidéoradiographie biplanaire (BVR) est une technique d’imagerie avancée permettant de comprendre le mouvement tridimensionnel des os squelettiques et des implants. Combinant des volumes d’images basés sur la densité et des vidéoradiographies du membre supérieur distal, la BVR est utilisée pour étudier le mouvement in vivo du poignet et de l’articulation radio-ulnaire distale, ainsi que les arthroplasties articulaires.

Résumé

Une mesure précise de la cinématique squelettique in vivo est essentielle pour comprendre le fonctionnement normal des articulations, l’influence de la pathologie, la progression de la maladie et les effets des traitements. Les systèmes de mesure qui utilisent des marqueurs de surface de la peau pour déduire le mouvement squelettique ont fourni des informations importantes sur la cinématique normale et pathologique, cependant, une arthrocinématique précise ne peut pas être obtenue à l’aide de ces systèmes, en particulier lors d’activités dynamiques. Au cours des deux dernières décennies, les systèmes de vidéoradiographie biplanaire (BVR) ont permis à de nombreux chercheurs d’étudier directement la cinématique squelettique des articulations pendant les activités de la vie quotidienne. Pour mettre en œuvre des systèmes BVR pour le membre supérieur distal, des vidéoradiographies du radius distal et de la main sont acquises à partir de deux sources de rayons X calibrées pendant qu’un sujet effectue une tâche désignée. Les positions tridimensionnelles (3D) des corps rigides sont calculées à partir des vidéoradiographies via un enregistrement optimal des projections de modèles 3D sur chaque vue BVR. Les modèles 3D sont des volumes d’image basés sur la densité de l’os spécifique dérivés de données de tomodensitométrie acquises indépendamment. Utilisant des processeurs graphiques et des systèmes de calcul haute performance, cette approche de suivi basée sur un modèle s’avère rapide et précise dans l’évaluation de la biomécanique du poignet et de l’articulation radio-ulnaire distale. Dans cette étude, nous avons d’abord résumé les études précédentes qui ont établi la concordance submillimétrique et submillimétrique de la BVR avec un système de capture de mouvement optique in vitro dans l’évaluation de la cinématique du poignet et de l’articulation radio-ulnaire distale. De plus, nous avons utilisé la BVR pour calculer le comportement du centre de rotation de l’articulation du poignet, pour évaluer le modèle d’articulation des composants de l’implant les uns sur les autres et pour évaluer le changement dynamique de la variance ulnaire pendant la pronosupination de l’avant-bras. À l’avenir, les os du carpe pourraient être capturés plus en détail grâce à l’ajout de détecteurs de rayons X à écran plat, d’un plus grand nombre de sources de rayons X (c’est-à-dire la vidéoradiographie multiplanaire) ou d’algorithmes avancés de vision par ordinateur.

Introduction

Une mesure précise de la cinématique squelettique in vivo est essentielle pour comprendre la fonction articulaire saine et remplacée, l’influence de la pathologie, la progression de la maladie et les effets des traitements. La quantification de la cinématique squelettique de manière non invasive à la surface de l’articulation (arthrocinématique) est cruciale pour comprendre les pathologies et les maladies articulaires, telles que l’arthrose, mais elle est techniquement difficile. Auparavant, les techniques qui utilisent des marqueurs de surface cutanée pour déduire le mouvement squelettique ont fourni des informations importantes sur la cinématique saine et pathologique. Cependant, il n’est pas possible d’obtenir une arthrocinématique précise à l’aide de ces techniques, en particulier lors d’activités dynamiques telles que les activités de la vie quotidienne. La précision de ces systèmes optiques est intrinsèquement limitée en raison du mouvement de la peau par rapport aux os sous-jacents, principale source d’erreur dans l’analyse du mouvement humain ^1,2.

Les méthodes de pointe actuelles pour quantifier la cinématique squelettique tridimensionnelle (3D) sont le suivi basé sur l’image, à savoir la vidéoradiographie biplan (BVR)³ et la tomodensitométrie en série (TDM) volumes⁴ et l’imagerie par résonance magnétique (IRM)⁵. Bien que les technologies habituelles de tomodensitométrie et d’IRM 3D soient très précises et accessibles dans de nombreux hôpitaux à travers le monde, elles sont incapables de mesurer le mouvement dynamique des articulations. Des techniques d’imagerie telles que la tomodensitométrie 4D⁶ et l’IRM dynamique⁷ ont été développées ces dernières années pour remédier à cette lacune. Cependant, ces méthodes exposent les patients à une forte dose de rayonnement ou souffrent d’une faible résolution temporelle.

En combinant de nouveaux algorithmes de vision par ordinateur et des systèmes à rayons X traditionnels, la BVR s’est avérée précise pour plusieurs articulations chez les animaux et les humains. résolu à l’aide d’algorithmes de suivi basés sur des marqueurs ou sur des modèles. Les approches basées sur des marqueurs suivent les billes de tantale insérées dans les os ou les tissus mous et sont optimales pour les tests sur les animaux et in vitro. Cependant, ils sont prohibitifs pour la recherche humaine in vivo . Heureusement, les améliorations apportées aux algorithmes de suivi basés sur des modèles offrent une alternative viable. Chez l’homme, les approches de suivi BVR basées sur un modèle consistent à préparer les ensembles d’images volumétriques acquises par TDM ou IRM dans une posture statique et à capturer les mouvements des intérêts dans le champ de vision de deux rayons X. La plupart des applications de suivi basées sur des modèles génèrent ensuite des radiographies reconstruites numériquement (RRC) de l’os ou de l’implant à partir des images statiques de TDM ou d’IRM et les associent à des vidéoradiographies améliorées à l’aide de mesures qui démontrent la similitude entre les RRM et les vidéoradiographies⁸. Ce processus s’appelle le « suivi » de l’os ou de l’implant.

Les principales variables de sortie du suivi des os ou des implants sont la cinématique des corps rigides, à partir de laquelle la cinématique des articulations, les allongements ligamentaires ^9,10, l’espacement des articulations comme substitut de l’épaisseur du cartilage¹¹, le contact articulaire^12,13 et d’autres biomarqueurs peuvent être calculés. Récemment, nous avons documenté la précision du suivi BVR basé sur un modèle dans le calcul de la biomécanique du poignet, de l’arthroplastie totale du poignet (TWA) et de l’articulation radio-ulnaire distale (DRUJ)^14,15. Dans la section suivante, un protocole détaillé de cette méthode validée pour l’étude du mouvement du poignet squelettique, de l’arthroplastie totale du poignet et de l’articulation radio-ulnaire distale lors de diverses tâches est présenté. Nous segmentons les volumes d’image basés sur la densité des os et des implants à partir des volumes d’images CT, suivons ces volumes d’image partiels dans les vidéoradiographies et déterminons des résultats tels que le centre de rotation, le modèle de contact et la variance ulnaire pour démontrer les forces et les limites de cette méthode.

Protocole

Cette étude a été approuvée par l’Institutional Review Board (IRB) de Lifespan - Rhode Island Hospital, un IRB accrédité par l’AAHRPP. Au total, 16 patients ont fourni un consentement éclairé signé conformément aux directives de l’établissement.

1. Acquisition de données

1. Tomodensitométrie (TDM)
   1. Préparez les échantillons ou les sujets pour le scanner.
      REMARQUE : Pour l’évaluation de la précision^14,15, 6 avant-bras intacts provenant de quatre spécimens de cadavres intacts (70,5 ± 12,3 ans, 4 droits, 2 femelles) et 1 avant-bras (49 ans, droit, femelle) avec une arthroplastie totale du poignet ont été utilisés. Pour l’étude in vivo ¹⁶, 10 volontaires avec des poignets sains (57,0 ± 5,2 ans, huit femmes, neuf dominants droitiers) et 6 patients (74,7 ± 5,6 ans, deux femmes, deux poignets droits, quatre poignets dominants) ayant subi une arthroplastie totale du poignet ont été recrutés.
   2. Pour imager chaque sujet⁴, acquérir une tomodensitométrie de l’avant-bras distal avec un diamètre de collecte de données de 25 cm et reconstruire avec un champ de vision de 20 cm à l’aide de l’algorithme de convolution Bone Plus à l’aide d’un scanner clinique à 16 coupes. Utilisez des réglages de tube de 80 kVp et 80 mA. Limitez l’exposition de l’avant-bras à environ 7 cm proximal de la surface du radius distal aux phalanges distales.
      REMARQUE : Dans cette étude, les dimensions du voxel de l’image étaient de 0,39 × 0,39 mm² dans le plan transversal de l’avant-bras et de 0,625 mm le long de l’axe long de l’avant-bras.
   3. Transférez les images CT acquises sur l’ordinateur du laboratoire à l’aide d’un protocole approuvé par HIPPA au format DICOM (Digital Imaging and Communications in Medicine).
2. Vidéoradiographie biplanaire (BVR)
   1. Dispositif expérimental (Figure 1). Pour permettre un diamètre de champ de vision capable de capturer le radius distal, le cubitus distal, les os du carpe et les métacarpiens ainsi qu’une perspective permettant de minimiser les occlusions osseuses dans diverses postures du poignet, positionnez les intensificateurs d’image (II) sous les avant-bras et réglez l’angle entre eux à 110° (sources de rayons X dirigées vers le bas vers les II correspondants à 55° par rapport au sol).
      1. Générez les rayons X en mode continu avec un réglage d’exposition compris entre 65 et 75 kV, en fonction de la visibilité des os/implants, et 80 mA. Ajustez les distances source-image des deux paires X-Ray-II à environ 130 cm et les distances source-main à 90 cm.
         REMARQUE : Dans cette étude, les sujets portaient un tablier de plomb de protection contre les rayons X, un collier de protection thyroïdienne et un dosimètre numérique en laboratoire. Tous les chercheurs se tenaient derrière un verre de protection en plomb pendant l’expérience.
      2. Placez un tréteau ou une structure similaire entre les II, et utilisez-le pour placer des objets (par exemple, une poignée de porte, un marteau, un pichet) pour les sujets.

Graphique 1. Dispositif expérimental. Veuillez cliquer ici pour voir une version agrandie de cette figure.

2. Grille d’absence de distorsion (Figure 2A). Étant donné que les champs électromagnétiques externes affectent le comportement des faisceaux de rayons X (c’est-à-dire les trajets d’électrons), les images radiographiques souffrent de distorsions non linéaires en coussinets et en forme de S¹⁷. Pour éliminer ces distorsions, utilisez l’image d’une pièce plate de tôle perforée de précision (c’est-à-dire une grille sans distorsion)^18,19. Placez des grilles sans distorsion sur les II et capturez leurs images avant et après les données de vidéoradiographie. Enregistrez ces images en tant qu’images TIFF 8 bits pour les deux sources de rayons X.
3. Étalonnage (Figure 2B). L’étalonnage des sources de rayons X est nécessaire pour calculer les coordonnées 3D des objets (os ou implants)²⁰. Utilisez des images d’un objet 3D à la géométrie connue, qui comporte des sphères radio-opaques régulièrement espacées, pour générer une relation entre la position et l’orientation de la source de rayons X dans l’espace 3D et l’espace de l’image (c’est-à-dire les radiographies). Un objet cubique (c’est-à-dire un cube d’étalonnage) fabriqué à partir de feuilles acryliques et d’entretoises avec des sphères d’acier est efficace et précis pour les réglages BVR¹⁹.
   1. Capturez des radiographies du cube^{d’étalonnage 18} dans différentes orientations 4 fois avant de réaliser l’expérience et 4 fois après l’expérience. Les étalonnages post-expérience servent à confirmer que la configuration expérimentale reste inchangée tout au long de l’expérience. Enregistrez toutes les images en tant qu’images TIFF 8 bits pour les deux sources de rayons X.

Graphique 2. A) Grille sans distorsion. B) Cube d’étalonnage et ses éléments de référence. Veuillez cliquer ici pour voir une version agrandie de cette figure.

4. Description des tâches. Définissez diverses tâches pour capturer le mouvement du poignet lors de la flexion-extension, de la déviation radiale-ulnaire et de la circumduction¹⁶. Décrire les tâches de mouvement DRUJ comme la rotation de l’avant-bras d’une posture de maintien à une pronation complète ou à une supination complète, définie comme la rotation d’une poignée en T dans le sens des aiguilles d’une montre et dans le sens inverse des aiguilles d’une montre²¹. Acquérir des vidéoradiographies du poignet dans sa pose neutre (c’est-à-dire une tâche statique neutre). Capturez la pose neutre avec le dos de la main coplanaire avec l’arrière de l’avant-bras.
   REMARQUE : Dans la posture de départ pour toutes les tâches, l’épaule était en adduction, le coude était soutenu au niveau de l’articulation et l’avant-bras était sans contrainte et dans sa rotation neutre. Les tâches de flexion-extension et de déviation radiale-ulnaire ont été décrites comme la rotation de la paume par rapport aux côtés palmaire-dorsal et radial-ulnaire de la main, respectivement. La circumduction a été décrite comme le mouvement du poignet qui visait à atteindre une amplitude de mouvement active maximale dans toutes les directions.
5. Acquisition de vidéoradiographies. Acquérez des vidéoradiographies à une fréquence de 200 Hz avec la vitesse d’obturation de la caméra réglée à 500 μs. Enregistrez deux secondes d’imagerie pour chaque tâche afin d’acquérir 2 répétitions. Stockez les vidéoradiographies acquises à partir de chaque source de rayons X sous forme de vidéo CINE dans un format 8 bits.
   REMARQUE : Dans cette étude, les radiographies avaient une taille d’image de 1760×1760 pixels et avaient une résolution d’environ 0,22 mm×0,22 mm par pixel.
   REMARQUE : Les sources de rayons X actuelles ont une option pour agrandir le champ de vision (FOV). Cependant, dans cette étude, le mode agrandi offert par le matériel n’a pas été utilisé. D’après la description de la tâche, cette fonctionnalité peut vous permettre d’acquérir des sorties de résolution plus élevée. Si cette fonction est utilisée, les étapes Grille de non-distorsion et Étalonnage doivent être capturées dans ce mode.

3. Capture optique de mouvement (OMC) - étude de précision
   1. Utilisez le protocole d’acquisition de données de capture de mouvement optique (OMC) décrit dans des études précédentes^14,15. Fixez directement deux grappes de marqueurs, chacune avec quatre sphères de marqueur rétroréfléchissantes de 9,5 mm de diamètre, au troisième métacarpien et au radius. Capturez le mouvement des marqueurs à l’aide de 8 caméras de capture de mouvement optiques. Synchronisez le démarrage de l’acquisition des données OMC avec le BVR par un déclencheur externe faible actif.

2. Traitement des données

1. Tomodensitométrie (TDM)
   1. Anonymisation DICOM. Anonymisez les DICOM et importez-les dans un logiciel de traitement d’images pour la modélisation 3D (par exemple, Mimics Materialise). Dans le logiciel Mimics Materialise, dans la barre d’outils Fichier , ouvrez la fenêtre Anonymiser . Lorsque la nouvelle fenêtre est ouverte, sélectionnez Anonymiser toutes les images dans l’onglet de gauche, puis appuyez sur le bouton Anonymiser en bas à droite.
   2. Segmentation (figure 3). Segmenter les os ou les implants d’intérêt en suivant une méthodologie décrite précédemment 16,22,23.
      1. En bref, dans la barre d’outils Segment , sélectionnez l’option Seuil , puis choisissez les ensembles de seuils prédéfinis pour l’os à l’ouverture de la fenêtre. Une fois le masque seuillé créé, utilisez l’option Masquer scindé et peignez manuellement les régions des os ou les composants de l’implant.
      2. Lorsque toutes les régions sont colorées, cliquez sur OK et attendez que le logiciel génère automatiquement les masques pour chaque os et le composant de l’implant. Vérifiez tous les masques pour vous assurer qu’ils couvrent complètement les régions prévues. Si ce n’est pas le cas, utilisez l’option Modifier les masques pour modifier les régions manuellement et une par une.
         REMARQUE : Pour étudier l’articulation du poignet, les 2e^et 3e^{métacarpiens} et les os du radius doivent être segmentés et séparés du reste des os. En option, le capitate et d’autres os peuvent être segmentés pour les étapes de post-traitement.
         REMARQUE : Pour étudier le DRUJ, le radius et le cubitus doivent être segmentés et séparés du reste des os.
         REMARQUE : Pour les implants, le composant radial et le composant carpien avec ses vis doivent être segmentés.

Graphique 3. Image par tomodensitométrie du poignet et modèles reconstruits du radius, du troisième métacarpien et du cubitus. Veuillez cliquer ici pour voir une version agrandie de cette figure.

3. Segmentation - composante radiale uniquement. En raison des artefacts de traînée basés sur la densité de la surface de la composante radiale dans les images CT, le modèle de composante radiale ne peut pas être segmenté à partir des images. Ainsi, utilisez un modèle CAO ou un scanner de surface 3D pour générer le modèle 3D du composant radial explanté. Enregistrez ce modèle dans la segmentation approximative de la composante radiale à l’aide de la méthode par défaut de Mimics Materialise.
   REMARQUE : Dans cette étude, nous avons utilisé un scanner de surface d’une résolution de 0,1 mm. L’implant a été placé sur une table et le scanner 3D a été déplacé manuellement autour de l’implant pour capturer toutes les surfaces externes.
4. Volumes d’images partiels. Le logiciel de suivi nécessite qu’un corps d’intérêt segmenté soit saisi en tant que volume d’image basé sur la densité. Pour créer ce volume partiel, générez un masque complet du volume d’image contenant toutes les valeurs grises en sélectionnant l’option Seuil dans la barre d’outils Segment . Ajustez les valeurs de seuil de la valeur minimale aux valeurs de gris maximales.
   1. Ensuite, utilisez l’option Opérateurs booléens et soustrayez le volume segmenté de l’os ou de l’implant d’intérêt du masque complet généré. Changez la couleur du masque de ce volume partiel en noir (RVB : 0.0/0.0/0.0) car le logiciel d’enregistrement n’accepte que les volumes d’images en niveaux de gris. Exportez le masque en tant que DICOM en sélectionnant Exporter DICOM... dans le logiciel. La sortie est un volume d’image de l’os isolé ou de l’implant.
5. Modèles tridimensionnels. Convertissez les masques segmentés en modèles de surface 3D maillés triangulés en sélectionnant l’option Calculer 3D , avec la résolution optimale. Exportez ces modèles au format VRML, contenant l’emplacement 3D des sommets (points) et leur ordre de connexions. Les modèles sont générés pour la construction du système de coordonnées et une visualisation plus poussée.

2. Prétraitement de la vidéoradiographie
   1. Importation XMALab¹⁸ . Pour chaque sujet, créez un nouvel essai dans le logiciel XMALab. Attribuez les chemins d’accès aux dossiers aux images de la grille de non-distorsion et à 8 ensembles d’images d’étalonnage dans la boîte de dialogue Nouvelle version d’essai .
      REMARQUE : Si les sources de rayons X sont inversées tout au long de la capture, l’option de retournement doit être cochée dans le logiciel XMALab.
   2. Dé-distorsion. Pour déformer les radiographies, une fois l’essai créé et le logiciel ayant ouvert les images de grille dans l’espace de travail Undistortion, cliquez sur le bouton Compute Undistortion . XMALab détermine et affiche automatiquement les images déformées des grilles. Une fois que les centres de tous les trous ont été détectés par le logiciel, remplacez l’espace de travail dans le menu en haut à gauche par l’espace de travail Étalonnage . Pour plus d’informations et pour obtenir de l’aide sur le débogage des problèmes potentiels, consultez la référence²⁴.
   3. Calibrage. À ce stade, le logiciel XMALab détermine les matrices d’étalonnage. Dans XMALab, les sources de rayons X sont modélisées comme des caméras à sténopé avec des paramètres intrinsèques et extrinsèques, tels que la distance focale, le centre de la source et les distances source-image, et l’étape d’étalonnage est un processus qui détermine ces métriques¹⁸. Pour calibrer, sélectionnez chacun des 4 éléments de référence placés dans le cube de calibration pour tous les ensembles d’images et pour chaque vue. Dans la configuration, ces références étaient Pyramid, Jack, Dome et Helix (Figure 2B). Si le logiciel ne trouve pas automatiquement les points, appuyez sur Compute Calibration. Une fois toutes les images d’étalonnage étalonnées, sélectionnez Optimiser la configuration de l’appareil photo.
   4. Exportation des matrices d’étalonnage. Après l’optimisation, triez les 8 jeux d’étalonnage (c’est-à-dire les 4 ensembles qui ont été capturés avant l’expérience et les 4 ensembles qui ont été capturés après l’expérience) et sélectionnez l’ensemble d’images avec la plus petite erreur non déformée. Trouvez cette valeur en appuyant sur l’option d’information dans le coin supérieur droit de chaque fenêtre d’étalonnage. Sur l’image avec l’erreur la plus faible, exportez les images en sélectionnant l’option MayaCams 2.0 dans le menu Fichier-Exportation . Pour plus d’informations et pour déboguer d’éventuels problèmes, consultez la référence²⁰.
   5. Vidéoradiographie sans distorsion. Changez l’espace de travail en Suivi des marqueurs et chargez les vidéoradiographies qui ont été capturées au format CNE en appuyant sur l’icône +. Une fois le processus de chargement terminé, exportez les images d’essai non déformées en sélectionnant Images d’essai non déformées dans le menu Exportation de fichiers .
3. Suivi basé sur un modèle (Figure 4)
   1. Amélioration des caractéristiques des radiographies/RRM (Figure 4A-B). L’amélioration des caractéristiques de l’image à l’aide de filtres est une étape nécessaire pour obtenir une correspondance précise. Définissez les paramètres des quatre filtres, y compris le contraste (détection de l’intensité), le Sobel (détection des bords), le gaussien (flou/plus lisse) et la netteté (mise en gras des bords) dans le logiciel. Ces filtres peuvent être sélectionnés en cliquant avec le bouton droit de la souris sur les objets Rad Renderer ou DRR Renderer affichés dans le widget en haut à gauche.
      REMARQUE : Dans cette étude, les bords des images des os sur les radiographies ont été améliorés à l’aide d’un filtre de Sobel avec un facteur d’échelle de 3 et une valeur de mélange de 0,4. De plus, un filtre de contraste avec un alpha (pour le contraste de l’image) de 2,5 et un bêta (pour la luminosité de l’image) de 0,9 a été utilisé pour améliorer encore la qualité des radiographies. Pour faire correspondre les DRR avec la radiographie, une valeur d’intensité de rayon de 0,35 a été utilisée avec un filtre Sobel avec une valeur de mélange de 0,1 et un facteur d’échelle de 1,7.
   2. Initialisation manuelle. Déplacez la RRC de chaque os et composant à l’aide de son point de pivot pour obtenir la meilleure correspondance possible entre la RRC et les radiographies. Le point de pivot aide les utilisateurs à faire pivoter et à traduire la DRR sélectionnée à partir du widget en bas à gauche. Cette étape est nécessaire car la position de l’os ou de l’implant doit être initialisée au moins une fois, et au moins pour la première image, avant d’effectuer l’optimisation. Appuyez sur W pour déplacer la DRR en mode traduction et appuyez sur E pour faire pivoter la DRR.
      REMARQUE : Il n’est généralement pas nécessaire de modifier l’emplacement du point de pivot, mais cela est possible à l’aide de l’option Déplacer le pivot de la barre d’outils.
   3. Optimisation. Appuyez sur C ou cliquez sur le bouton Suivre le courant dans la barre d’outils pour trouver la meilleure correspondance entre la DRR et les radiographies. La procédure par défaut du logiciel utilise la fonction de coût de corrélation croisée normalisée (NCC) et la méthode d’optimisation de l’essaim de particules et les procédures simplex descendantes^26,27. Ce processus est également appelé l’étape de correspondance ou de suivi. La procédure d’optimisation peut également être exécutée via la boîte de dialogue de suivi. Pour ce faire, appuyez sur l’option Boîte de dialogue de suivi dans la barre d’outils. Lorsque la nouvelle fenêtre est ouverte, remplacez les options À partir de l’image et À l’image sur 0, puis appuyez sur OK.
      REMARQUE : La NCC est une fonction qui mesure la similitude entre les radiographies et les DRR, et une valeur NCC plus proche de 0 a démontré une meilleure correspondance. Cependant, l’obtention de cette valeur dépend entièrement des filtres d’image et des qualités de capture. La valeur de NCC peut être observée dans la « Fenêtre du terminal ».
   4. Processus de jumelage - Premier passage. Une fois que la première image est automatiquement optimisée, ouvrez la boîte de dialogue de suivi, remplacez l’image De l’image par 0 et remplacez l’image de destination par le dernier numéro de l’image (situé en bas à droite). Le processus de correspondance par défaut du logiciel utilise la trame actuelle comme position initiale de la DRR pour le processus de correspondance. Le processus d’optimisation utilise la méthode d’optimisation par essaim de particules (PSO) dans une région de 3 mm et 3° de la position initiale.
   5. Processus d’appariement - Passage final (figure 4C). Étant donné que le processus de correspondance automatique peut échouer en raison d’une initialisation ou de filtres médiocres, vérifiez manuellement toutes les trames pour vous assurer que l’optimisation a été exécutée correctement. Si la RRC et les radiographies ne sont pas correctement appariées, répétez les étapes décrites aux sections 2.3.2 et 2.3.3.
      REMARQUE : Le suivi du poignet a d’abord été effectué avec un modèle qui incluait les deuxième et troisième métacarpiens (MC2-MC3) pour réduire les effets de chevauchement osseux et améliorer la précision de l’appariement. Après avoir localisé la position et la rotation combinées de la RRC dans les radiographies, les emplacements MC2-MC3 appariés et optimisés ont été transformés en volume d’image partiel MC3, et les positions initiales du MC3 ont été amorcées. En utilisant ces positions de tête, le MC3 a été suivi. À l’instar de MC3, le rayon, le cubitus, la composante radiale et la composante carpienne dans les vidéoradiographies pendant toutes les tâches ont également été suivis.

Graphique 4. A) Radiographie capturée d’une source de rayons X avec radiographies reconstruites numériquement (RRC) des os. B) Radiographie améliorée (filtrée) et RRC. C) TRRC appariées après le processus d’optimisation. Veuillez cliquer ici pour voir une version agrandie de cette figure.

6. Réduction des données. En fonction de la fréquence d’images et de la qualité des images radiographiques, des approches de filtrage telles que la méthode de la moyenne mobile (avec une portée de 5 images) ont pu être utilisées. Pour ce faire, dans le menu de la barre d’outils, sélectionnez Suivi-Boîte de dialogue avancée, puis cliquez sur le bouton Lisser . Cette approche a été utilisée pour toutes les captures de cette étude.
7. Exportation des positions osseuses/implantaires. Une fois que chaque os est adapté aux radiographies et que toutes les images sont optimisées, exportez les matrices de transformation 4x4 des DRR dans le monde de la radiologie en cliquant sur le bouton Enregistrer le suivi et en sélectionnant la sortie du format de matrice. Cette matrice de transformation définit la rotation et la translation du volume d’image partiel de l’os/de l’implant dans le système de coordonnées du monde à rayons X calibré et peut être utilisée pour un traitement supplémentaire et le calcul du mouvement. Le système de coordonnées du volume partiel de l’image doit être transféré au système de coordonnées anatomiques pour une analyse biomécanique significative.

3. Analyse des données

1. Systèmes de coordonnées (Figure 5). Pour rendre compte de la cinématique, construisez des systèmes de coordonnées (CS) pour les os et les composants de l’implant en fonction de leurs repères anatomiques et géométriques. Une fois ces CS construits, enregistrez les modèles qui se trouvent dans le nouveau CS dans les modèles de surface qui ont été générés dans la Section 2.1.4 à l’aide de l’algorithme itératif du point le plus proche. Ce recalage un-à-un calcule une nouvelle matrice de transformation 4x4 qui sera utilisée pour les calculs cinématiques.
   REMARQUE : Dans cette étude, le CS du radius, du cubitus, du troisième métacarpien, du composant carpien, du composant radial, du chapeau de polyéthylène et du capitate ont été créés à l’aide de méthodes qui ont été décrites précédemment 16,22,23.

Graphique 5. Systèmes de coordonnées des os et des composants de l’implant. Veuillez cliquer ici pour voir une version agrandie de cette figure.

2. Cinématique. Dans cette étude, on définit le mouvement du poignet comme le mouvement du MC3 par rapport au rayon, et les données cinématiques ont été rapportées par rapport à la position neutre du poignet.
   1. Pour l’analyse DRUJ, calculez les paramètres d’axe de mouvement hélicoïdal (HAM)²⁸ pour décrire la cinématique comme le mouvement relatif du rayon CS dans le CS ulnaire, par rapport à la pose neutre.
   2. Utilisez l’analyse HAM pour déterminer le centre de rotation (COR) du poignet en flexion-extension et en déviation radiale-ulnaire.
   3. Pour décrire le décalage proximal-distal unidimensionnel du COR, définissez le COR projeté pour les poignets sains comme un point sur l’axe central de l’os capité qui avait la distance la plus courte de l’axe de rotation de la vis¹⁶.
   4. Décomposez la rotation HAM dans le système de coordonnées ulnaires à l’aide de l’axe de vis HAM, et calculez la pronation (positive) et la supination (négative) de l’avant-bras. La pose neutre a été définie lorsque la pronation-supination (PS) était nulle.
   5. Enfin, définissez la variance ulnaire (UV) comme la distance entre la fovéa ulnaire et l’encoche sigmoïde sur le rayon.
3. Modèle d’articulation
   1. Calculez les modèles de contact entre le capuchon en polyéthylène et le composant radial à l’aide de la cinématique du poignet et des champs de distance spécifiques au composant²⁹.
   2. Calculez les champs de distance pour chaque modèle de surface 3D sous la forme d’un tableau volumétrique de distances signées par rapport à la surface.
   3. Calculez les valeurs de proximité des champs de distance sur la surface du capuchon en polyéthylène et de la composante radiale pour chaque posture.
   4. Pour obtenir la surface de contact entre les composants, ajustez ces valeurs de proximité à la résolution du système d’acquisition (le seuil d’exclusion de distance a été sélectionné à 0,70 mm dans cette étude).
4. Analyse statistique
   1. Utilisez l’analyse de Bland-Altman pour évaluer l’exactitude en démontrant le biais (différences moyennes) et la limite de concordance à 95 % (biais ± 1,96 écart-type) entre le BVR et l’OMC^{14, 30}.
   2. Utilisez des modèles linéaires généralisés (GLM) pour décrire le comportement du centre de rotation du poignet sain et pour trouver une relation entre la variation des UV et de la PS (UV = p₁×PS² + p₂×PS), où p₁ et p₂ étaient des coefficients optimisés à l’aide de la régression linéaire. Utilisez la métrique d’erreur quadratique moyenne (RMSE) pour évaluer la précision du modèle.

Résultats

Le choix d’un logiciel de recalage d’images 2D-3D pour le suivi basé sur le modèle dépend en partie de l’accès à l’unité de processeur graphique (GPU) et aux systèmes de calcul haute performance (HPC). Ces programmes ont des pipelines différents, et à l’heure actuelle, il n’y a pas de méthodologie commune entre les programmes. Dans cette étude, nous utilisons Autoscoper, un programme open source de recalage d’images 2D-à-3D développé à l’Université Brown

Discussion

La vidéoradiographie biplanaire (BVR) est une méthode basée sur l’image qui peut être utilisée pour mesurer le mouvement de l’os et de l’implant dans le poignet et l’articulation radio-ulnaire distale avec une précision submillimétrique et submillimétrique. Dans les études que nous avons décrites ici, la BVR a été utilisée pour identifier un modèle précis de COR projeté pour un poignet sain ainsi que des modèles de contact TWA. De tels résultats peuvent éclaire...

matériels

Name	Company	Catalog Number	Comments
3D Surface Scanner	Artec 3D	Artec Space SpiderTM	Luxembourg
Autoscoper	Brown University	https://simtk.org/projects/autoscoper	https://doi.org/10.1016/j.jbiomech.2019.05.040
CT Scanner	General Electric (GE)	Lightspeed 16	Milwaukee, WI, USA
Geomagic Wrap 3D	3DSystems	Version 2017	Rock Hill, SC, USA
Graphics Processing Unit (GPU)	Nvidia	GeForce GTX 1080	CUDA-enabled GPU
High-speed Video Cameras	Phantom	Version 10	Vision Research, Wayne, NJ, USA
Image Intensifier	Dunlee	40 cm diameter	Aurora, IL, USA
ImageJ	Open-source (Brown University)	https://imagej.net/Fiji	https://doi.org/10.1038/nmeth.2019
Matlab	The MathWorks, Inc.	R2017a to R2020a	Natick, MA, USA
Mimics	Materialise	Version 19.0 to 22.0	Leuven, Belgium
Motion Capture Cameras	Qualisys	Oqus 5+	Gothenburg, Sweden
Pulsed X-ray Generators	EMD Technologies	EPS 45–80	Saint-Eustache, Quebec, QC, Canada
Undistortion Grid	McMaster-Carr	9255T641	Steel Perforated Sheet Staggered Holes, 0.048" Thk, 0.125" Hole Dia, 36" X 40"
Wrist Implant (In-vitro Study)	Integra LifeSciences	Universal 2	Plainsboro, NJ, USA
Wrist Implant (In-vivo Study)	Integra LifeSciences	Freedom	Plainsboro, NJ, USA
WristViz	Open-source (Brown University)	https://github.com/DavidLaidlaw/WristVisualizer/tree/master	Open-source software
X-ray Tubes	Varian Medical Systems	Model G-1086	Palo Alto, CA, USA
XMALab	Open-source (Brown University)	https://www.xromm.org/xmalab/	https://doi.org/10.1242/jeb.145383