$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Este estudo foi aprovado pelo Conselho de Revisão Institucional do Hospital Geral da Universidade de Shenzhen (aprovação do IRB nº KYLL-2026-077-1). A exigência de consentimento informado por escrito foi dispensada devido ao desenho retrospectivo do estudo.
1. Coorte de estudo
Pacientes com AVC isquêmico agudo (AIS) que passaram por ressonância magnética multiparamétrica foram identificados retrospectivamente a partir do banco de dados clínico institucional. O fluxo de trabalho geral do estudo, incluindo particionamento de conjuntos de dados, desenvolvimento de modelos de imagem, construção de modelos clínicos e fusão multimodal, está ilustrado na Figura 1.

Figura 1: Fluxo de trabalho do framework multimodal de previsão de resultados. (A) Aquisição de dados. Pacientes com AVC isquêmico agudo (AIS) que atenderam aos critérios de inclusão pré-definidos foram identificados retrospectivamente. O conjunto de dados foi dividido em uma coorte de validação de treinamento (n = 250), uma coorte interna independente de teste (n = 50) e uma coorte externa de teste (n = 37). A validação cruzada estratificada de cinco vezes foi realizada dentro da coorte de validação de treinamento. (B) Desenvolvimento de modelos de imagem. Sequências de ressonância magnética multiparamétricas, incluindo imagens ponderadas por difusão (DWI), coeficiente de difusão aparente (ADC) e recuperação de inversão atenuada por fluido T2 (T2-FLAIR), foram processadas usando uma arquitetura híbrida de rede neural convolucional 3-Transformador de Visão (CNN-ViT) para previsão de resultados. (C) Desenvolvimento de modelos clínicos. Variáveis clínicas estruturadas foram usadas para treinar modelos de aprendizado de máquina para previsão do resultado funcional em 90 dias. (D) Estratégia de fusão multimodal. As previsões geradas pelos modelos de imagem e clínicos foram integradas usando um meta-aprendiz de regressão logística empilhada para produzir previsões de resultado final. Por favor, clique aqui para ver uma versão ampliada desta figura.
Exames de ressonância magnética foram realizados utilizando sistemas clínicos de ressonância magnética. Uma bobina de cabeça dedicada em matriz faseada era usada para aquisição de imagem. Todas as sequências foram adquiridas no plano axial usando posicionamento consistente dos cortes entre as modalidades. A imagem ponderada por difusão (DWI) foi adquirida usando uma sequência de imagem ecoplanar spin-echo de tiro único, com tempo de repetição (TR) de 3.000–5.000 ms e tempo de eco (TE) de 80–90 ms. A sensibilização por difusão foi aplicada usando valores b de 0 e 800–1.000 s/mm 2 em pelo menos 3 direções ortogonais. O campo de visão variava de 220 a 240 mm, com uma matriz de 128 × 128. A espessura da fatia era de 5–6 mm com uma folga entre a interfaia de 1–1,5 mm. Foram obtidas médias de dois a quatro sinais.
Mapas de coeficientes de difusão aparente (ADC) foram gerados automaticamente a partir de dados DWI na estação de trabalho scanner usando ajuste monoexponencial baseado nos valores b adquiridos. Os valores de ADC foram calculados voxel a par e exportados para análise quantitativa. Imagens de recuperação por inversão atenuada por fluido T2 (T2-FLAIR) foram adquiridas usando uma sequência de recuperação por inversão com TR de 8.000–10.000 ms, TE de 80–140 ms e tempo de inversão de 2.200–2.600 ms. O campo de visão variava de 220 a 240 mm, com tamanho de matriz de 192 × 192 a 256 × 256. A espessura da fatia era de 4–5 mm com uma folga entre a interfaia de 1–1,5 mm.
2. Pré-processamento de dados
Os dados clínicos foram importados de planilhas estruturadas contendo identificadores de casos, escores de desfechos do mRS, rótulos de divisão de conjunto de dados e variáveis demográficas e clínicas. Identificadores de casos foram padronizados para garantir consistência com nomes de arquivos de imagem. As variáveis clínicas incluíram características demográficas, gravidade do AVC medida usando a Escala de AVC dos Institutos Nacionais de Saúde (NIHSS), fatores de risco vasculares e condições comórbidas registradas na admissão.
O desfecho primário foi o status funcional 90 dias após o início do AVC, medido usando a Escala de Rankin (mRS) modificada. Um desfecho favorável foi definido como mRS ≤ 2 e um desfecho desfavorável como mRS > 2. O conjunto de dados foi dividido aleatoriamente em coortes de validação de treinamento e coortes internas independentes de teste, utilizando amostragem estratificada baseada na distribuição mRS para preservar o equilíbrio de resultados. Uma coorte externa de teste foi processada separadamente e não foi usada durante o desenvolvimento do modelo. Dentro da coorte de validação de treinamento, foi aplicada a validação cruzada estratificada k-fold para manter distribuições consistentes de resultados entre as dobras.
Todos os volumes DWI, ADC e T2-FLAIR foram pré-processados antes do treinamento do modelo para garantir consistência espacial e numérica entre as modalidades. O mesmo pipeline de pré-processamento foi aplicado às coortes de validação de treinamento, teste interno e teste externo sem modificações. As imagens foram reorientadas para a orientação canônica do RAS e convertidas em matrizes de ponto flutuante. A resolução espacial no plano foi reamostrada para 256 × 256 usando interpolação linear. A dimensão do plano passante foi padronizada para 20 fatias usando uma estratégia baseada no centro: volumes contendo mais de 20 fatias foram cortados no centro, enquanto volumes com menos de 20 fatias foram simetricamente preenchidos a zero. O volume final foi de 256 × 256 × 20.
A normalização da intensidade foi realizada independentemente para cada volume usando normalização do escore z:

onde x denota intensidade voxel, μ é a intensidade média do volume, e σ é o desvio padrão. Se σ = 0, a normalização não foi aplicada para evitar instabilidade numérica. Os volumes processados eram salvos no formato NIfTI com uma matriz afim padronizada para análise de aprendizado profundo a jusante.
3. Desenvolvimento de modelos clínicos
Variáveis clínicas estruturadas foram usadas para desenvolver modelos de aprendizado de máquina para previsão de resultados. Os preditores candidatos incluíam características demográficas, fatores de risco vasculares, etiologia do AVC e medidas clínicas de gravidade basal. Variáveis contínuas foram imputadas usando valores medianos e padronizadas. Variáveis categóricas foram imputadas usando a categoria mais frequente e codificadas usando codificação one-hot.
Múltiplos algoritmos de aprendizado de máquina foram avaliados, incluindo regressão logística, floresta aleatória, aumento de gradiente, máquina de vetores de suporte (SVM), aumento extremo de gradiente e modelos de máquinas de aumento de gradiente de luz. O desenvolvimento do modelo seguiu uma estrutura estratificada de validação cruzada cinco vezes dentro da coorte de validação de treinamento para estimar o desempenho de generalização. As previsões finais para as coortes de teste interna e externa foram geradas por meio da média das previsões de modelos treinados em cada dobra de validação cruzada. Nenhuma amostra externa de teste foi usada durante a validação cruzada ou seleção do modelo.
4. Arquitetura de modelos de deep learning
Uma arquitetura híbrida CNN-ViT tridimensional foi implementada para previsão de resultados a partir de volumes multimodais de ressonância magnética. A rede foi projetada para combinar extração local de características espaciais com modelagem contextual global dentro de um arcabouço unificado. Os volumes de entrada consistiam em imagens multicanais tridimensionais processadas de ponta a ponta.
A extração de características foi inicialmente realizada usando uma espinha dorsal convolucional hierárquica tridimensional composta por 4 estágios. Cada estágio incluía 2 camadas convolucionais com tamanho de grão de 3 × 3 × 3 e preenchimento de 1 voxel, seguida de normalização em lote e ativação linear retificada por unidades. A resolução espacial foi progressivamente reduzida usando camadas de max-pooling tridimensional aplicadas após os três primeiros estágios, enquanto a profundidade do canal de características aumentou em cada nível para capturar representações semânticas de nível superior. A regularização de dropout (taxa de dropout = 0,1) foi aplicada após o estágio convolucional final para reduzir o sobreajuste. A espinha dorsal convolucional transformava o volume de entrada de tamanho C × D × H × W em uma representação compacta de características de alto nível com dimensões espaciais reduzidas.
O mapa de características resultante foi remodelado em uma sequência de tokens ao achatar as dimensões espaciais de modo que o número de tokens correspondesse a:
N = D' x H' x W'
Relacionamentos contextuais globais entre tokens foram modelados usando camadas de codificadores Transformer compostas por redes de autoatenção e feedforward multi-cabeças. A autoatenção era calculada como:

onde Q, K e V denotam matrizes de consulta, chave e valor, respectivamente, e d representa a dimensão de incorporação. A normalização e o dropout de camadas foram aplicados dentro de cada camada de codificador para melhorar a estabilidade do treinamento. O módulo Transformer consistia em 3 camadas de codificador com 8 cabeças de atenção e uma dimensão de embedding de 256.
Após a codificação por transformador, a representação correspondente ao token de classificação era extraída e normalizada. Uma camada linear totalmente conectada produzia uma única saída logit para classificação binária.
A arquitetura híbrida proposta CNN-ViT foi intencionalmente projetada como um modelo leve e eficiente em termos de parâmetros para equilibrar a capacidade representacional e o risco de sobreajuste. O modelo compreendia 3,79 milhões de parâmetros treináveis (aproximadamente 14,4 MB na precisão fp32), incluindo 1,38 milhão na espinha dorsal convolucional e 2,37 milhões no codificador Transformer.
5. Treinamento de modelos de imagem
O treinamento do modelo de imagem foi realizado usando uma estrutura de validação cruzada estratificada de cinco vezes para preservar a distribuição de resultados entre as dobras, ao mesmo tempo em que melhora a robustez da estimativa de desempenho. O conjunto de dados foi dividido em uma coorte de validação de treinamento e uma coorte interna independente de teste, utilizando amostragem estratificada baseada na distribuição dos resultados. Dentro da coorte de validação de treinamento, foi aplicada a validação cruzada estratificada em cinco fases. Para cada dobra, o modelo de imagem foi treinado usando o subconjunto de treinamento e avaliado usando o subconjunto correspondente de validação, enquanto a coorte de teste mantida foi reservada exclusivamente para a avaliação final de desempenho.
A otimização do modelo foi realizada usando uma estrutura de deep learning em uma estação de trabalho equipada com GPU com o otimizador AdamW, taxa de aprendizado de 3 × 10⁻5 e decaimento de peso de 3 × 10⁻4. O treinamento era realizado em um lote de 8 para até 200 épocas. A entropia cruzada binária com logits era usada como função de perda.
Para corrigir o desequilíbrio de classes, um fator de ponderação de classe positiva foi calculado para cada dobra com base na razão de amostras negativas para positivas e incorporado à função de perda. O clipping de normas de gradiente com norma máxima de 0,5 foi aplicado para melhorar a estabilidade numérica durante a otimização. O treinamento automático de precisão mista foi possibilitado para melhorar a eficiência computacional.
A parada precoce foi implementada quando o desempenho de validação não melhorou em pelo menos 1 × 10⁻4 ao longo de 30 épocas consecutivas. O checkpoint de melhor desempenho do modelo de cada dobra foi mantido. Após a conclusão de todas as dobras, as previsões para as coortes de teste interna e externa foram geradas usando cada modelo específico de dobra, e as probabilidades finais foram obtidas por meio da média das previsões entre os 5 modelos para gerar resultados em conjunto.
6. Modelo de fusão multimodal
Uma estratégia de fusão empilhada foi implementada para integrar previsões derivadas de imagens com informações clínicas estruturadas. O modelo de imagem de aprendizado profundo e o modelo de predição clínica serviram como aprendizes base, e suas probabilidades previstas foram usadas como características de entrada para um meta-aprendiz de regressão logística. Recursos adicionais de interação, incluindo o produto e a diferença absoluta das probabilidades previstas, foram incorporados para capturar informações complementares entre imagens e previsões clínicas.
Para evitar vazamento de informações, o meta-aprendiz foi treinado usando probabilidades previstas fora de fold geradas pela coorte de validação de treinamento. Probabilidades fora de dobra validadas cruzadamente dos modelos de imagem e clínico foram combinadas por identificador de paciente para construir o conjunto de dados de treinamento do meta-aprendiz. Para as coortes de teste interna e externa, as probabilidades correspondentes do conjunto de teste dos modelos de imagem e clínico foram usadas como entradas para o meta-aprendiz treinado para gerar probabilidades fundidas. O meta-aprendiz treinado foi aplicado a ambas as coortes de teste sem reajuste.
7. Estudo de ablação
Experimentos de ablação foram realizados para avaliar as contribuições das sequências individuais de ressonância magnética e do componente Transformador Visual, utilizando configurações idênticas de treinamento e avaliação ao modelo principal. As ablações de sequência incluíam modelos de sequência única, modelos de deixar uma sequência fora e o modelo multiparamétrico completo. Para avaliar a contribuição do módulo Vision Transformer, a arquitetura híbrida proposta foi adicionalmente comparada com uma linha base apenas CNN, na qual o codificador Transformer foi removido, preservando a mesma espinha dorsal convolucional.
8. Análise estatística
O desempenho do modelo foi avaliado separadamente na coorte de teste independente interna e na coorte externa de teste, usando a área sob a curva característica operacional do receptor (AUC) como métrica principal de discriminação. Curvas de característica operacional do receptor (ROC) foram construídas usando probabilidades previstas geradas por cada modelo. Sensibilidade, especificidade e precisão geral foram calculadas adicionalmente para caracterizar o desempenho da classificação.
Resultados binários foram gerados usando limiares determinados de acordo com o índice de Youden. Sensibilidade, especificidade e precisão foram posteriormente calculadas no limiar ideal. Todas as análises estatísticas, desenvolvimento de modelos de aprendizado de máquina e treinamento de modelos de aprendizado profundo foram implementados usando pacotes padrão de computação científica e software de aprendizado de máquina.