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

Resumo Técnico: Gêmeos Digitais Articulados Específicos de Paciente a partir de uma Única Tomografia Computadorizada de Corpo Inteiro

Definição do Problema
Modelos anatômicos específicos de pacientes derivados de imagens médicas são fundamentais para o planejamento cirúrgico, intervenções guiadas por imagem e desenvolvimento de algoritmos. No entanto, os modelos existentes são predominantamente estáticos; eles capturam a anatomia do paciente apenas na configuração específica presente durante o exame. Essa limitação é particularmente crítica em imagens radiográficas, onde a aparência projetada é uma função tanto da geometria de imagem quanto da pose do paciente. Embora os modelos de nível populacional (ex: XCAT) permitam variação de pose, eles carecem de geometria específica do paciente. Por outro lado, as Radiografias Digitalmente Reconstruídas (DRRs) derivadas de TC preservam a geometria específica do paciente, mas estão presas à pose de aquisição, tornando-as "controláveis em termos de visualização", mas não "controláveis em termos de pose". Há uma necessidade de um gêmeo digital de paciente que seja simultaneamente individualizado e articulável em novas configurações corporais.

Metodologia
Os autores propõem um pipeline de três estágios para construir um gêmeo digital articulado específico de paciente a partir de uma única tomografia computadorizada (TC) de corpo inteiro:

1. Modelagem Anatômica do Paciente:

   • Segmentação: Um volume de TC de corpo inteiro é processado usando o TotalSegmentator para extrair máscaras semânticas do corpo, ossos e órgãos. Marching Cubes e suavização de Taubin são aplicados para gerar malhas de superfície.
   • Ajuste de SMPL: Um modelo de corpo humano paramétrico (SMPL) é ajustado à anatomia derivada da TC para criar um arcabouço cinemático. Para abordar a discrepância entre a pose supina da TC e a pose canônica em pé do SMPL (e para evitar o ajuste não confiável a deformações de tecidos moles), o SMPL é inicializado usando 25 pontos-chave esqueléticos anotados manualmente (estilo OpenPose).
   • Otimização: Uma otimização de múltiplas fases minimiza uma função de perda combinando distância de chamfer (consistência geométrica), perda de inclusão (garantindo que o esqueleto permaneça dentro do envelope corporal) e priors de forma e pose. O ajuste prioriza o envolvimento da anatomia esquelética sobre o ajuste estrito à superfície externa do corpo.

2. Vinculação Cinemática Consciente da Anatomia:

   • Uma árvore cinemática reduzida com raiz na pelve é definida.
   • Os ossos derivados da TC são vinculados a juntas ou grupos de juntas semanticamente correspondentes. Diferente da animação SMPL padrão, que usa linear blend skinning para deformação de superfície, este método aplica transformações rígidas a cada osso, preservando a rigidez esquelética.
   • Âncoras de Juntas: Em vez de usar as localizações de juntas do SMPL, âncoras de juntas específicas da TC são estimadas a partir das regiões de superfície mais próximas entre os ossos filho e pai para servirem como centros de rotação.
   • Órgãos e tecidos moles são atribuídos a segmentos anatômicos mais grosseiros (ex: tórax, pelve) e movem-se rigidamente com esses segmentos.

3. Retargeting de Pose:

   • Poses corporais alvo (especificadas como parâmetros SMPL) são redirecionadas (retargeted) para a anatomia do paciente.
   • O movimento relativo das juntas da pose de aquisição para a pose alvo é extraído da cadeia cinemática do SMPL e propagado através da árvore cinemática reduzida da TC via cinemática direta (forward kinematics).
   • Os ossos sofrem transformações rígidas baseadas em suas âncoras, enquanto os órgãos seguem o movimento de seus segmentos anatômicos atribuídos.

Principais Resultados
O método foi avaliado em três sujeitos de TC de corpo inteiro do conjunto de dados NMDID, com DRRs renderizadas usando DeepDRR.

• Precisão de Ajuste Corporal: O arcabouço SMPL ajustado alcançou uma distância de chamfer global de 15,8 ± 4,0 mm contra a superfície corporal derivada da TC. Crucialmente, ele envolveu 95,9 ± 1,8% do esqueleto derivado da TC, garantindo que a geometria esquelética permanecesse dentro do arcabouço cinemático.
• Fidelidade de Recomposição: Quando o gêmeo foi recomposto na pose de aquisição original, ele preservou as principais estruturas radiográficas. A comparação entre as DRRs da TC bruta e as DRRs do gêmeo recomposto rendeu um SSIM global de 0,872 ± 0,016 e PSNR de 18,5 ± 1,4 dB. Os autores observam que os valores mais baixos de PSNR devem-se principalmente ao preenchimento constante de Unidades Hounsfield (HU) usado para o envelope de tecido mole, em vez de erros estruturais.
• Articulação de Poses Não Vistas: Os gêmeos foram articulados em 28 poses não vistas derivadas de sequências de movimento AMASS. A autossustentabilidade estrutural foi mantida com uma taxa de envolvimento de 94,4 ± 0,4%, indicando que não houve penetração catastrófica do esqueleto-envelope ou quebra do rig. Os resultados qualitativos mostraram mudanças dependentes da pose nas radiografias sintéticas, mantendo a geometria esquelética específica do paciente.

Principais Contribuições e Significância
O artigo apresenta uma prova de conceito para estender a simulação de TC baseada em visualização, que é estática, para gêmeos anatômicos controláveis em termos de pose. A principal contribuição é um método que funde a anatomia esquelética derivada de TC com um envelope corporal paramétrico ajustado (SMPL) para permitir a articulação em novas configurações, preservando a geometria esquelética específica do paciente como corpos rígidos.

Os autores afirmam que esta abordagem preenche uma lacuna prática na simulação baseada em TC: ela permite a geração de radiografias sintéticas diversas em termos de pose a partir de um único exame estático. Essa capacidade é significativa para futuros estudos de posicionamento e imagem sintética, oferecendo uma maneira de testar procedimentos e algoritmos em variadas configurações de pacientes sem a necessidade de múltiplos exames.

Limitações
Os autores explicitamente reconhecem as seguintes limitações:

• A avaliação é limitada a três sujeitos, e as poses não vistas são avaliadas via autossustentabilidade estrutural, em vez de anatomia ou radiografias de referência (ground-truth).
• O método modela primariamente a articulação esquelética; a deformação de órgãos e tecidos moles é modelada via segmentos rígidos grosseiros, em vez de deformação não rígida.
• Como o SMPL é um modelo de superfície paramétrico, pode não capturar totalmente as formas corporais ou restrições anatômicas específicas do paciente, podendo levar a desajustes locais (ex: ao redor das mãos) ou autointerseções em poses extremas.