Para começar, adquira imagens de ressonância magnética de tempo de eco ultracurto do pulmão durante a respiração livre. Importe os dados e as trajetórias com espaçamento k para o MATLAB. Descarte as primeiras 1000 projeções para garantir que os dados atinjam uma magnetização de estado estacionário.
Em seguida, execute a reconstrução da imagem usando uma transformada de Fourier rápida não uniforme para um tamanho de matriz de 96 por 96 por 96. Use aproximadamente 200 projeções correspondentes a 0,6 a 0,8 segundos de dados. Em seguida, reconstrua e armazene imagens de todos os elementos da bobina, bem como a imagem combinada da bobina final.
Na imagem combinada da bobina, selecione uma fatia coronal que mostre claramente o diafragma. Uma vez que a fatia coronal é selecionada, visualize as imagens individuais da bobina para esta fatia e selecione um ou dois elementos da bobina que melhor mostram o diafragma. Agora reconstrua apenas os dados dos elementos da bobina usando uma janela deslizante para gerar imagens com resolução temporal de aproximadamente 0,5 segundo.
Use as primeiras 200 projeções para reconstruir uma imagem usando uma transformada rápida de Fourier não uniforme e armazene apenas a fatia do diafragma. Desloque 100 projeções e reconstrua uma imagem adicional armazenando a fatia do diafragma. Agora, selecione uma linha sobre o diafragma na primeira das imagens da janela deslizante.
Visualize o movimento respiratório visualizando este navegador respiratório para todas as projeções. Determine a localização do diafragma para todos os navegadores respiratórios e use esse local para rotular as projeções como pertencentes a um determinado compartimento respiratório. Em seguida, identifique o silo com o maior número de projeções correspondentes à expiração final e escolha-o para reconstrução.
Use um filtro exponencial para fornecer um peso de um para projeções dentro do compartimento primário e um peso drasticamente reduzido para projeções dentro de diferentes compartimentos respiratórios. Em seguida, use a caixa de ferramentas de reconstrução avançada de Berkeley para reconstruir uma imagem de alta resolução no compartimento respiratório desejado. Calcule os pesos de compensação de densidade usando a combinação de densidade iterativa.
Dimensione os pesos de compensação de densidade pelos pesos de passagem suave. Em seguida, dimensione os dados com base na compensação de densidade e nos pesos de passagem suave. Agora execute uma transformada de Fourier rápida não uniforme básica para facilitar a combinação da bobina.
Converta a imagem de transformada rápida de Fourier não uniforme em espaço k quadriculado para combinação de bobina. Em seguida, gere uma matriz de combinação de bobinas e use-a para combinar bobinas para os dados brutos e o espaço k em grade e estimar as sensibilidades da bobina. Em seguida, usando a compensação de densidade ponderada, os dados combinados da bobina e os mapas de sensibilidade da bobina, realizam a reconstrução do sentido comprimido de imagem paralela.
As imagens geradas na expiração final usando gating baseado em imagem e baseado em espaço k mostraram uma visualização clara do diafragma com gating baseado em imagem demonstrando compensação de movimento superior. O soft-gating aumentou a nitidez das imagens de inspiração, reduzindo os artefatos de subamostragem em comparação com o hard-gating. Tanto o gating baseado em imagem quanto o k-space detectaram com sucesso formas de onda respiratórias durante a respiração regular, com o gating baseado em imagem produzindo resultados mais claros em condições respiratórias irregulares.
