O protocolo de imagem 3D PREFUL mais simples inclui duas sequências. Um localizador anatômico + uma medida com respiração livre utilizando aquisição em formato 3D stack-of-stars com ângulo dourado. A trajetória stack-of-stars é mostrada na Figura 1. Essa trajetória combina amostragem radial (no plano) com amostragem cartesiana linear (na direção z). Inicialmente, são adquiridas amostras na direção da partição kz, seguidas pela rotação do ângulo dourado para amostrar o plano kx-ky. O esquema de aquisição por ângulo dourado promove auto-gatilhamento com cobertura uniforme do espaço-k, facilitando a reconstrução de imagens com dados fortemente subamostrados. O exame completo não ultrapassa 10 minutos. Mais detalhes sobre a aquisição 3D PREFUL estão apresentados na Tabela de Materiais.

Figura 1: Trajetória em pilha rotacionada de estrelas com proporção áurea. Observe a amostragem radial aplicada ao longo da dimensão no plano (kx-ky) e a amostragem cartesiana aplicada ao longo da direção de fatia (kz). Este esquema resulta em uma cobertura cilíndrica do espaço-k adquirido. Neste exemplo, são mostradas três partições, incluindo 32 projeções radiais por partição. Clique aqui para visualizar uma versão maior desta figura.
Uma lista dos parâmetros típicos de RM utilizados para a aquisição de imagens 3D PREFUL em um scanner de 1,5 T é apresentada na Tabela 1B. Para o scanner de 3 T, são propostos os parâmetros de sequência listados na Tabela 2.
Após a aquisição e a transferência dos dados serem concluídas, inicia-se a etapa de reconstrução da imagem. Esse passo é realizado por um código MATLAB pertencente à propriedade intelectual (PI) de uma empresa (BioVisioneers GmbH) e não pode ser compartilhado publicamente. O script é totalmente automatizado assim que o caminho para os dados brutos é fornecido. Uma visão esquemática do procedimento de reconstrução é apresentada na Figura 2.

Figura 2: Visão esquemática do método de ressonância magnética funcional pulmonar tridimensional com resolução de fase (PREFUL). Inicialmente, os dados são adquiridos utilizando uma aquisição por ressonância magnética com respiração livre, de duração de 8 minutos, com trajetória tipo pilha-de-estrelas. Após a transferência dos dados do aparelho de RM para uma unidade de processamento, são reconstruídas imagens tridimensionais de baixa resolução com resolução temporal de aproximadamente 100 ms. O parênquima pulmonar de cada imagem é segmentado e o volume pulmonar é calculado. As informações de volume pulmonar são então utilizadas como sinal de sincronização. Com base na amplitude e na fase do sinal de sincronização, as projeções radiais são classificadas em compartimentos respiratórios que abrangem um ciclo respiratório completo. A etapa de classificação é seguida pela reconstrução dinâmica de imagens de resolução total, que posteriormente são registradas ao nível do final da inspiração. Após diversas etapas de pós-processamento, calcula-se o ciclo de ventilação regional (RVent) e realiza-se a análise das imagens, incluindo a extração de parâmetros. O ciclo RVent é avaliado mediante o cálculo do laço fluxo-volume (FVL) para cada vóxel. Os mapas de ventilação FVL-CM são extraídos por meio da avaliação do FVL de cada vóxel em relação a um FVL de referência saudável, sendo a similaridade avaliada utilizando a métrica de correlação cruzada. Para ambos RVent e FVL-CM, são quantificados os valores globais percentuais totais de defeito de ventilação (VDP). Além disso, a dinâmica do ciclo RVent é analisada por meio da análise de tempo-para-o-pico, resultando em um mapa de parâmetro de tempo-para-o-pico da ventilação (VTTP). Adicionalmente, o desvio do pico esperado de ventilação em 50% do ciclo RVent é quantificado no mapa VTTPDev. Clique aqui para visualizar uma versão ampliada desta figura.
Após o procedimento de reconstrução da imagem, todas as imagens são alinhadas espacialmente a um único estado fixo. O próprio registro é realizado na CPU utilizando o Advanced Normalization Tools (ANTs27) ou o pacote de registro Forsberg28. Embora o pacote ANTs permaneça o padrão ouro para tarefas de registro de imagens em ressonância magnética e tomografia computadorizada, o registro Forsberg permite um procedimento até 10 vezes mais rápido, com duração máxima de 9 minutos, obtendo resultados comparáveis19. O pacote de registro pode ser escolhido com base nas prioridades do usuário e nas unidades computacionais disponíveis. Após a conclusão do registro, as imagens morfológicas registradas são novamente reinvertidas, de modo que o pulmão apareça escuro na imagem em escala de cinza.
Na Figura 3, é apresentada a próxima etapa do fluxo de trabalho, que envolve a segmentação do parênquima pulmonar. Primeiro, o parênquima pulmonar é segmentado a partir da imagem no final da inspiração utilizando uma rede neural convolucional com arquitetura nnUnet29. A segmentação do parênquima pulmonar é seguida por um reconhecimento de vasos30, que é excluído da máscara final de segmentação.

Figura 3: Resultados exemplares da segmentação baseada em aprendizado profundo para um indivíduo do sexo masculino de 32 anos. A primeira linha exibe imagens morfológicas de oito cortes coronais representativos. Nas segunda e terceira linhas, podem-se observar, respectivamente, a máscara correspondente do parênquima pulmonar e a máscara final com a exclusão dos vasos. Clique aqui para visualizar uma versão ampliada desta figura.
Um exemplo de segmentação incorreta do parênquima pulmonar é mostrado na Figura 4. É essencial inspecionar visualmente as segmentações baseadas em aprendizado profundo, e, caso sejam consideradas insatisfatórias, devem-se considerar correções manuais para aumentar a precisão da máscara final do parênquima pulmonar.

Figura 4: Resultados exemplares de segmentação incorreta baseada em aprendizado profundo para um sujeito masculino de 57 anos. A primeira linha exibe imagens morfológicas de oito cortes coronais representativos. Nas segunda e terceira linhas, podem-se observar, respectivamente, a máscara do parênquima pulmonar correspondente e a máscara final com a exclusão dos vasos. Como evidente, diversas regiões fibroticas são erroneamente reconhecidas como vasos ou estruturas não pulmonares. Essas imprecisões foram corrigidas manualmente, conforme demonstrado na quarta linha. Clique aqui para visualizar uma versão ampliada desta figura.
Após várias etapas de filtragem descritas no protocolo, são calculados os parâmetros indiretos da ventilação. A ressonância magnética 3D PREFUL gera mapas quantitativos para ventilação regional estática (RVent), métricas dinâmicas de correlação fluxo-volume (FVL-CM) e dois parâmetros baseados na análise do tempo para pico da ventilação (VTTP), conforme ilustrado na Figura 5. Esta figura apresenta oito cortes coronais de um homem saudável de 32 anos. Observe a distribuição homogênea esperada de todos os parâmetros de ventilação.

Figura 5: Ressonância magnética 3D PREFUL de um voluntário saudável (homem de 32 anos). Mapas representativos dos parâmetros morfológicos (linha superior) e de ventilação por ressonância magnética 3D PREFUL (segunda à quinta linha) de um voluntário saudável (homem de 32 anos). A ventilação regional estática é representada pela ventilação regional (RVent), enquanto a dinâmica da ventilação é avaliada utilizando a métrica de correlação do laço fluxo-volume (FVL-CM), o tempo para pico da ventilação (VTTP) e o desvio do VTTP (VTTPDev), parâmetros que avaliam a dinâmica da ventilação. Conforme esperado, observam-se valores homogêneos de ventilação em todos os mapas de parâmetros ventilatórios. Clique aqui para visualizar uma versão maior desta figura.
Para simplificar os mapas de ventilação, são derivados mapas de defeito de ventilação (VD) para RVent e FVL-CM, o que permite uma interpretação mais rápida dos resultados. Os mapas exemplares de VD são apresentados na Figura 6. Para ambos VDRVent e VDFVL-CM, os valores de VDP foram de 3,6% e 3,0%, respectivamente, situando-se dentro da faixa normal saudável. Idealmente, e dependendo da idade dos voluntários saudáveis, o valor de VDP não deve exceder 10%. Os parâmetros mencionados acima (RVent, FVL-CM e seus mapas de VD) foram validados em diversos estudos11,12,35,36 e são sensíveis na detecção de doenças, bem como na detecção de efeitos induzidos por terapia20,34,37,38,39.

Figura 6: Mapas representativos de RVent e FVL-CM, incluindo seus mapas de defeito de ventilação. São apresentados os mapas de RVent (linha superior) e FVL-CM (segunda linha), incluindo seus mapas de defeito de ventilação (terceira e quarta linhas), derivados da ressonância magnética 3D PREFUL de um voluntário saudável (homem de 32 anos). Nos mapas de defeito de ventilação (VD), as regiões saudáveis são verdes e as áreas com déficit de ventilação são marcadas em vermelho. Clique aqui para visualizar uma versão maior desta figura.
Vários estudos de ressonância magnética pulmonar foram relatados em campos magnéticos de 1,5 T e 3 T. Embora exista uma vantagem teórica do campo de 3 T devido à maior relação sinal-ruído (SNR), essa vantagem pode ser superada pelos efeitos mais acentuados de suscetibilidade magnética em 3 T. A influência exata das intensidades do campo magnético na qualidade de imagem dos mapas de ventilação 3D PREFUL é atualmente desconhecida. Aqui, apresentamos resultados de viabilidade (Figura 7) obtidos com um voluntário saudável (homem de 35 anos) utilizando um scanner de ressonância magnética de 3 T. Observe a aparência mais heterogênea dos parâmetros de ventilação em 3 T quando comparada à de 1,5 T (Figura 5).

Figura 7: Ressonância magnética 3D PREFUL de um voluntário saudável (homem de 35 anos). Mapas representativos de parâmetros morfológicos (linha superior) e de ventilação por ressonância magnética 3D PREFUL (segunda à quinta linha) para um voluntário saudável (homem de 35 anos) em 3T. A ventilação regional estática é representada pela ventilação regional (RVent), enquanto a dinâmica da ventilação é avaliada utilizando a métrica de correlação do laço fluxo-volume (FVL-CM), o tempo para pico da ventilação (VTTP) e o desvio do VTTP (VTTPDev), parâmetros que avaliam a dinâmica da ventilação. Clique aqui para visualizar uma versão maior desta figura.
A comparação/validação da ventilação tridimensional derivada do PREFUL com medições mais diretas ainda não foi publicada. Existem diversos estudos que relatam uma correlação positiva entre os valores de VDP e a sobreposição espacial de defeitos de ventilação entre a técnica bidimensional de PREFUL e a ressonância magnética com 129Xe11,12,14 e com 19F35. Uma comparação recente entre a ventilação derivada do PREFUL tridimensional e medições diretas de ventilação usando ressonância magnética com 19F em pacientes com doença pulmonar obstrutiva crônica (DPOC), asma e voluntários saudáveis revelou uma correlação moderada a forte em nível global40. Apresenta-se um exemplo de comparação (Figura 8) de um paciente diagnosticado com doença pulmonar obstrutiva crônica (mulher de 54 anos, VEF1 = 42% do valor previsto), que foi examinado usando PREFUL tridimensional e ressonância magnética com lavagem de 19F.

Figura 8: RNM PREFUL 3D de uma paciente com DPOC de 54 anos. São apresentadas imagens morfológicas (primeira linha), parâmetros de ventilação e respectivos mapas de defeito de ventilação de uma paciente com DPOC de 54 anos (VEF1 = 42 % pred., CVF = 102% pred.) obtidos por meio de RNM PREFUL 3D (segunda e terceira linha) e RNM de 19F (quarta e quinta linha), juntamente com uma comparação dos mapas de defeito de ventilação provenientes de ambos os métodos (última linha). No geral, em todos os cortes, o coeficiente de Sørensen-Dice foi de 25,5% nas áreas defeituosas e de 80,6% nas regiões saudáveis, resultando em uma sobreposição espacial total de 69,2%. Notavelmente, há uma correlação visual entre as áreas de ventilação correspondentes nas regiões saudáveis e defeituosas (representadas em verde-escuro). Clique aqui para visualizar uma versão maior desta figura.
A ressonância magnética 3D PREFUL pode ser usada para medir respostas regionais à terapia20. Como exemplo, Figura 9 mostra uma comparação de três cortes do pulmão de um indivíduo com fibrose cística (mulher de 43 anos) no momento basal (à esquerda, VEF1 = 94% do valor previsto) e após terapia com modulador de CFTR (à direita, VEF1 = 112% do valor previsto). Como ambas as medições são espacialmente correspondentes, a técnica 3D PREFUL permite a análise da resposta regional ao tratamento, conforme ilustrado na parte inferior da Figura 9. Observe os valores aumentados da métrica de correlação do laço fluxo-volume, bem como a redução nos defeitos de ventilação após a terapia.

Figura 9: Medições 3D PREFUL de uma paciente com fibrose cística (FC) de 43 anos de idade. Mapas exemplares de marcadores de ventilação nas medições 3D PREFUL pré-tratamento (esquerda) e pós-tratamento (direita) de uma paciente com FC de 43 anos de idade. O VDPFVL-CM diminuiu de 18,0% (pré-tratamento) para 3,6% (pós-tratamento). FEV1 % pred. pré-tratamento: 94%, FEV1 % pred. pós-tratamento: 112%. LCI pré-tratamento: 10,48, LCI pós-tratamento: 9,39. Os mapas correspondentes de resposta ao tratamento são mostrados na parte inferior. Observe as áreas verdes indicando a resolução da ventilação após a terapia. Clique aqui para visualizar uma versão maior desta figura.
Tabela 1: Lista dos parâmetros típicos de RM utilizados para a aquisição do localizador e do PREFUL 3D em um scanner de 1,5 T. (A) Localizador e (B) aquisição de PREFUL 3D. Clique aqui para baixar esta tabela.
Tabela 2: Parâmetros propostos de RNM para aquisição 3D PREFUL em um scanner de 3T. Clique aqui para baixar esta tabela.
Tabela 3: Relatório estatístico exemplar dos parâmetros 3D PREFUL para um voluntário saudável (homem de 35 anos) Clique aqui para baixar esta tabela.