Patient-Specific Articulated Digital Twins from a Single Full-Body CT Scan

Résumé technique : Jumeaux numériques articulés spécifiques au patient à partir d'un seul scanner CT corps entier

Énoncé du problème
Les modèles anatomiques spécifiques aux patients dérivés de l'imagerie médicale sont fondamentaux pour la planification chirurgicale, les interventions guidées par l'image et le développement d'algorithmes. Cependant, les modèles existants sont principalement statiques ; ils ne capturent l'anatomie du patient que dans la configuration spécifique présente lors du scanner. Cette limitation est particulièrement critique en imagerie radiographique, où l'apparence projetée est une fonction de la géométrie d'imagerie et de la pose du patient. Bien que les fantômes à l'échelle de la population (ex. : XCAT) permettent la variation de pose, ils manquent de géométrie spécifique au patient. Inversement, les radiographies reconstruites numériquement (DRR) dérivées du scanner CT préservent la géométrie spécifique au patient mais sont verrouillées à la pose d'acquisition, ce qui les rend « contrôlables par la vue » mais pas « contrôlables par la pose ». Il existe un besoin pour un jumeau numérique du patient qui soit à la fois individualisé et articulable dans de nouvelles configurations corporelles.

Méthodologie
Les auteurs proposent un pipeline en trois étapes pour construire un jumeau numérique articulé spécifique au patient à partir d'un seul scanner CT corps entier :

1. Modélisation anatomique du patient :

   • Segmentation : Un volume CT corps entier est traité à l'aide de TotalSegmentator pour extraire des masques sémantiques pour le corps, les os et les organes. Les algorithmes Marching Cubes et le lissage de Taubin sont appliqués pour générer des maillages de surface.
   • Ajustement SMPL : Un modèle corporel humain paramétrique (SMPL) est ajusté à l'anatomie dérivée du CT pour créer un échafaudage cinématique. Pour traiter l'écart entre la pose de décubitus dorsal du CT et la pose debout canonique du SMPL (et pour éviter un ajustement peu fiable aux déformations des tissus mous), le SMPL est initialisé à l'aide de 25 points clés squelettiques annotés manuellement (style OpenPose).
   • Optimisation : Une optimisation multiphasique minimise une fonction de perte combinant la distance de Chamfer (cohérence géométrique), la perte d'inclusion (garantissant que le squelette reste à l'intérieur de l'enveloppe corporelle) et des priorités sur la forme et la pose. L'ajustement privilégie l'encerclement de l'anatomie squelettique plutôt qu'une correspondance étroite avec la surface externe du corps.

2. Liaison cinématique sensible à l'anatomie :

   • Un arbre cinématique réduit, dont la racine est le pelvis, est défini.
   • Les os dérivés du CT sont liés à des articulations ou des groupes d'articulations sémantiquement correspondants. Contra-rentionnellement à l'animation SMPL standard qui utilise le linear blend skinning pour la déformation de surface, cette méthode applique des transformations rigides à chaque os, préservant la rigidité squelettique.
   • Ancres articulaires : Au lieu d'utiliser les localisations des articulations SMPL, des ancres articulaires spécifiques au CT sont estimées à partir des régions de surface les plus proches entre les os enfants et parents pour servir de centres de rotation.
   • Les organes et les tissus mous sont assignés à des segments anatomiques plus grossiers (ex. : thorax, pelvis) et se déplacent rigidement avec ces segments.

3. Retargeting de pose :

   • Les poses corporelles cibles (spécifiées par les paramètres SMPL $\theta$) sont re-projetées (retargeted) sur l'anatomie du patient.
   • Le mouvement relatif des articulations de la pose d'acquisition vers la pose cible est extrait de la chaîne cinématique SMPL et propagé à travers l'arbre cinématique réduit du CT via une cinématique directe.
   • Les os subissent des transformations rigides basées sur leurs ancres, tandis que les organes suivent le mouvement de leurs segments anatomiques assignés.

Résultats clés
La méthode a été évaluée sur trois sujets CT corps entier provenant du jeu de données NMDID, avec des DRR rendus via DeepDRR.

• Précision de l'ajustement du corps : L'échafaudage SMPL ajusté a atteint une distance de Chamfer globale de 15,8 ± 4,0 mm par rapport à la surface corporelle dérivée du CT. Crucialement, il a encerclé 95,9 ± 1,8 % du squelette dérivé du CT, garantissant que la géométrie squelettique restait à l'intérieur de l'échafaudage cinématique.
• Fidélité de la recomposition : Lorsque le jumeau était recomposé à la pose d'acquisition originale, il préservait les structures radiographiques majeures. La comparaison entre les DRR bruts du CT et les DRR du jumeau recomposé a produit un SSIM global de 0,872 ± 0,016 et un PSNR de 18,5 ± 1,4 dB. Les auteurs notent que les valeurs de PSNR plus faibles sont principalement dues au remplissage constant de l'unité Hounsfield (HU) utilisé pour l'enveloppe des tissus mous plutôt qu'à des erreurs structurelles.
• Articulation de poses non vues : Les jumeaux ont été articulés dans 28 poses non vues dérivées des séquences de mouvement AMASS. La cohérence structurelle de soi a été maintenue avec un taux d'encerclement de 94,4 ± 0,4 %, indiquant l'absence de pénétration catastrophique du squelette dans l'enveloppe ou de rupture de la structure (rig breakdown). Les résultats qualitatifs ont montré des changements de pose dépendants dans les radiographies synthétiques tout en conservant la géométrie squelettique spécifique au patient.

Contributions et significations clés
L'article présente une preuve de concept pour étendre la simulation CT statique et contrôlable par la vue vers des jumeaux anatomiques contrôlables par la pose. La principale contribution est une méthode qui fusionne l'anatomie squelettique dérivée du CT avec un enveloppe corporelle paramétrique ajustée (SMPL) pour permettre l'articulation dans de nouvelles configurations tout en préservant la géométrie squelettique spécifique du patient comme corps rigides.

Les auteurs affirment que cette approche comble une lacune pratique dans la simulation basée sur le CT : elle permet la génération de radiographies synthétiques diverses en termes de pose à partir d'un seul scan statique. Cette capacité est significative pour les futures études d'imagerie et de positionnement synthétiques, offrant un moyen de tester des procédures et des algorithmes dans des configurations de patients variées sans nécessiter de multiples scans.

Limites
Les auteurs reconnaissent explicitement certaines limites :

• L'évaluation est limitée à trois sujets, et les poses non vues sont évaluées via la cohérence structurelle de soi plutôt que par une anatomie ou des radiographies de référence (ground-truth).
• La méthode modélise principalement l'articulation squelettique ; la déformation des organes et des tissus mous est modélisée via des segments rigides grossiers plutôt que par une déformation non rigide.
• Comme le SMPL est un modèle de surface paramétrique, il peut ne pas capturer pleinement les formes corporelles ou les contraintes anatomiques spécifiques au patient, ce qui peut entraîner des écarts locaux (ex. : autour des mains) ou des auto-intersections lors de poses extrêmes.