El protocolo de imagen 3D PREFUL más sencillo incluye dos secuencias: un localizador anatómico y una medición con respiración libre utilizando la adquisición 3D tipo pila de estrellas con ángulo dorado. La trayectoria de pila de estrellas se muestra en la Figura 1. Esta trayectoria combina muestreo radial (en el plano) con muestreo cartesiano lineal (en la dirección z). Inicialmente, se adquieren muestras en la dirección del particionado kz, seguidas de la rotación según el ángulo dorado para muestrear el plano kx-ky. El esquema de adquisición con ángulo dorado favorece la autogatillación con una cobertura uniforme del espacio k, facilitando la reconstrucción de imágenes a partir de datos altamente submuestreados. El examen completo no supera los 10 minutos. Más detalles sobre la adquisición 3D PREFUL se encuentran en la Tabla de Materiales.

Figura 1: Trayectoria de pila de estrellas rotada según la proporción áurea. Obsérvese el muestreo radial aplicado a lo largo de la dimensión en el plano (kx-ky) y el muestreo cartesiano aplicado a lo largo de la dirección de corte (kz). Este esquema da como resultado una cobertura cilíndrica del espacio-k adquirido. En este ejemplo, se muestran tres particiones, incluyendo 32 proyecciones radiales por partición. Haga clic aquí para ver una versión más grande de esta figura.
Se presenta una lista de los parámetros típicos de RM utilizados para la imagen 3D PREFUL en un escáner de 1,5 T en la Tabla 1B. Para el escáner de 3 T, se proponen los parámetros de secuencia que figuran en la Tabla 2.
Después de que se completan la adquisición de datos y la transferencia de datos, comienza la parte de reconstrucción de imágenes. Este paso se realiza mediante un código de MATLAB que pertenece a la propiedad intelectual (IP) de una empresa (BioVisioneers GmbH) y no puede compartirse abiertamente. El script es completamente automatizado una vez que se proporciona la ruta hacia los datos brutos. En Figura 2 se muestra un esquema general del procedimiento de reconstrucción.

Figura 2: Vista esquemática del método de resonancia magnética funcional pulmonar tridimensional con resolución de fase (PREFUL). En primer lugar, se adquieren los datos mediante una adquisición por resonancia magnética de 8 minutos de duración durante respiración libre, utilizando una trayectoria tipo «pila de estrellas». Una vez transferidos los datos desde el escáner de RM a una unidad de procesamiento, se reconstruyen imágenes tridimensionales de baja resolución con una resolución temporal de aproximadamente 100 ms. Se segmenta el parénquima pulmonar de cada imagen y se calcula el volumen pulmonar. Esta información de volumen pulmonar se utiliza posteriormente como señal de sincronización. En función de la amplitud y la fase de la señal de sincronización, las proyecciones radiales se clasifican en compartimentos respiratorios que abarcan un ciclo respiratorio completo. A esta clasificación le sigue la reconstrucción dinámica de imágenes de resolución completa, que posteriormente se registran al nivel de fin de la inspiración. Tras varios pasos de posprocesamiento, se calcula el ciclo de ventilación regional (RVent) y se realiza el análisis de la imagen, incluida la extracción de parámetros. El ciclo RVent se evalúa calculando el bucle flujo-volumen (FVL) para cada vóxel. Los mapas de ventilación FVL-CM se obtienen mediante la evaluación de cada FVL de vóxel frente a un FVL de referencia sano, utilizando como medida de similitud la métrica de correlación cruzada. Para RVent y FVL-CM, se cuantifican los valores globales del porcentaje total de defecto de ventilación (VDP). Además, la dinámica del ciclo RVent se analiza mediante el análisis del tiempo hasta el pico, obteniendo así un mapa de parámetros de tiempo hasta el pico de ventilación (VTTP). Asimismo, se cuantifica la desviación respecto al pico esperado de ventilación en el 50 % del ciclo RVent en el mapa VTTPDev. Haga clic aquí para ver una versión ampliada de esta figura.
Después del procedimiento de reconstrucción de la imagen, todas las imágenes se alinean espacialmente a un único estado fijo. El registro en sí se realiza en la CPU utilizando las Herramientas Avanzadas de Normalización (ANTs27) o el paquete de registro Forsberg28. Aunque el paquete ANTs sigue siendo el estándar de oro para tareas de registro de imágenes en resonancia magnética y tomografía computarizada, el registro Forsberg permite un procedimiento de registro hasta 10 veces más rápido, con tiempos de hasta 9 minutos, obteniendo resultados comparables19. El paquete de registro puede elegirse según las prioridades del usuario y las unidades de procesamiento disponibles. Una vez completado el registro, las imágenes morfológicas registradas se reinvierten nuevamente para que el pulmón aparezca oscuro en la imagen en escala de grises.
En la Figura 3 se muestra el siguiente paso del proceso, que consiste en la segmentación del parénquima pulmonar. Primero, el parénquima pulmonar se segmenta a partir de la imagen al final de la inspiración utilizando una red neuronal convolucional con arquitectura nnUnet29. A la segmentación del parénquima pulmonar sigue un reconocimiento de vasos30, que se excluye de la máscara final de segmentación.

Figura 3: Resultados ejemplares de la segmentación basada en aprendizaje profundo para un sujeto masculino de 32 años. La fila superior muestra imágenes morfológicas de ocho cortes coronales representativos. En las filas segunda y tercera, se pueden observar la máscara correspondiente del parénquima pulmonar y la máscara final con la exclusión de los vasos, respectivamente. Haga clic aquí para ver una versión más grande de esta figura.
Un ejemplo de segmentación incorrecta del parénquima pulmonar se muestra en la Figura 4. Es fundamental inspeccionar visualmente las segmentaciones basadas en aprendizaje profundo y, si se consideran insatisfactorias, se deben considerar correcciones manuales para mejorar la precisión de la máscara final del parénquima pulmonar.

Figura 4: Resultados ejemplares de una segmentación incorrecta basada en aprendizaje profundo para un sujeto masculino de 57 años. La fila superior muestra imágenes morfológicas de ocho cortes coronales representativos. En las segunda y tercera filas, se pueden observar la máscara correspondiente del parénquima pulmonar y la máscara final con la exclusión de los vasos, respectivamente. Como es evidente, varias regiones fibroticas son reconocidas erróneamente como vasos o estructuras no pulmonares. Estas inexactitudes han sido corregidas manualmente, como se muestra en la cuarta fila. Haga clic aquí para ver una versión más grande de esta figura.
Después de varios pasos de filtrado descritos en el protocolo, se calculan los indicadores de ventilación. La resonancia magnética 3D PREFUL genera mapas cuantitativos de ventilación regional estática (RVent), métricas dinámicas de correlación flujo-volumen (FVL-CM) y dos parámetros basados en el análisis del tiempo hasta el pico de ventilación (VTTP), como se ilustra en la Figura 5. Esta figura muestra ocho cortes coronales de un hombre sano de 32 años. Obsérvese la distribución homogénea esperada de todos los parámetros de ventilación.

Figura 5: Resonancia magnética 3D PREFUL de un voluntario sano (hombre de 32 años de edad). Mapas representativos de morfología (fila superior) y parámetros de ventilación de resonancia magnética 3D PREFUL (segunda a quinta fila) para un voluntario sano (hombre de 32 años de edad). La ventilación regional estática se representa mediante la ventilación regional (RVent), mientras que la dinámica de la ventilación se evalúa utilizando la métrica de correlación del bucle flujo-volumen (FVL-CM), el tiempo hasta el pico de ventilación (VTTP) y la desviación del VTTP (VTTPDev), que son parámetros de ventilación que evalúan la dinámica respiratoria. Como era de esperar, se observan valores homogéneos de ventilación en todos los mapas de parámetros de ventilación. Haga clic aquí para ver una versión más grande de esta figura.
Para simplificar los mapas de ventilación, se derivan mapas de defecto de ventilación (VD) para RVent y FVL-CM, lo que permite una interpretación más rápida de los resultados. Los mapas de VD ejemplares se presentan en la Figura 6. Para VDRVent y VDFVL-CM, los valores de VDP fueron de 3,6 % y 3,0 %, respectivamente, dentro del rango normal saludable. Idealmente, y dependiendo de la edad de los voluntarios sanos, el valor de VDP no debería superar el 10 %. Los parámetros mencionados anteriormente (RVent, FVL-CM y sus mapas de VD) se han validado en varios estudios11,12,35,36 y son sensibles para la detección de enfermedades, así como para la detección de efectos inducidos por la terapia20,34,37,38,39.

Figura 6: Mapas representativos de RVent y FVL-CM, incluyendo sus mapas de defecto de ventilación. Se muestran los mapas de RVent (fila superior) y FVL-CM (segunda fila), incluyendo sus mapas de defecto de ventilación (tercera y cuarta fila) obtenidos mediante resonancia magnética 3D PREFUL en un voluntario sano (hombre de 32 años). En los mapas de defecto de ventilación (VD), las regiones sanas son verdes, y las áreas con déficit de ventilación se marcan en rojo. Haga clic aquí para ver una versión más grande de esta figura.
Se han publicado varios estudios de resonancia magnética pulmonar a intensidades de campo magnético de 1,5 T y 3 T. Aunque existe una ventaja teórica del 3 T debido al mayor cociente señal-ruido (SNR), esta ventaja puede verse superada por efectos más pronunciados de susceptibilidad magnética a 3 T. La influencia exacta de las intensidades del campo magnético sobre la calidad de imagen de los mapas de ventilación 3D PREFUL es actualmente desconocida. Aquí presentamos resultados de viabilidad (Figura 7) obtenidos en un voluntario sano (hombre de 35 años) utilizando un escáner de RM de 3 T. Obsérvese la apariencia más heterogénea de los parámetros de ventilación a 3 T en comparación con 1,5 T (Figura 5).

Figura 7: Resonancia magnética PREFUL 3D de un voluntario sano (hombre de 35 años de edad). Mapas representativos de los parámetros morfológicos (fila superior) y de ventilación de la resonancia magnética PREFUL 3D (filas segunda a quinta) para un voluntario sano (hombre de 35 años de edad) a 3T. La ventilación regional estática se representa mediante la ventilación regional (RVent), mientras que la dinámica de la ventilación se evalúa utilizando la métrica de correlación del bucle flujo-volumen (FVL-CM), el tiempo hasta el pico de ventilación (VTTP) y la desviación del VTTP (VTTPDev), parámetros que evalúan la dinámica de la ventilación. Haga clic aquí para ver una versión más grande de esta figura.
Actualmente no se ha publicado ninguna comparación o validación de la ventilación tridimensional obtenida mediante PREFUL con mediciones más directas. Existen varios estudios que informan una correlación positiva entre los valores de VDP y la superposición espacial de los defectos de ventilación entre la técnica bidimensional de PREFUL y la RM con 129Xe11,12,14 y la RM con 19F35. Una comparación reciente entre la ventilación derivada del PREFUL tridimensional y mediciones directas de ventilación mediante RM con 19F en pacientes con enfermedad pulmonar obstructiva crónica (EPOC), asma y voluntarios sanos reveló una correlación moderada a fuerte a nivel global40. A continuación se presenta un ejemplo de comparación (Figura 8) en un paciente diagnosticado con enfermedad pulmonar obstructiva crónica (mujer de 54 años, VEF1 = 42% del valor predicho), evaluado mediante PREFUL tridimensional y RM de lavado con 19F.

Figura 8: Resonancia magnética 3D PREFUL de una paciente femenina de 54 años con EPOC. Se presentan imágenes morfológicas (fila superior), parámetros de ventilación y mapas respectivos de defectos de ventilación de una paciente femenina de 54 años con EPOC (VEF1 = 42 % pred., CVF = 102 % pred.) obtenidos mediante resonancia magnética 3D PREFUL (segunda y tercera fila) y resonancia magnética con 19F (cuarta y quinta fila), junto con una comparación de los mapas de defectos de ventilación de ambos métodos (última fila). En general, en todas las secciones, el coeficiente de Sørensen-Dice fue del 25,5 % en áreas defectuosas y del 80,6 % en regiones sanas, lo que resulta en una superposición espacial total del 69,2 %. Destacable es la correlación visual entre las áreas de ventilación coincidentes en regiones sanas y defectuosas (representadas en verde oscuro). Haga clic aquí para ver una versión más grande de esta figura.
La resonancia magnética PREFUL 3D puede utilizarse para medir las respuestas regionales al tratamiento20. Como ejemplo, la Figura 9 muestra una comparación de tres cortes del pulmón de un sujeto con fibrosis quística (mujer de 43 años) en la línea base (a la izquierda, FEV1 = 94 % del valor predicho) y después del tratamiento con moduladores del CFTR (a la derecha, FEV1 = 112 % del valor predicho). Dado que ambas mediciones están espacialmente alineadas, la técnica PREFUL 3D permite el análisis regional de la respuesta al tratamiento, como se ilustra en la parte inferior de la Figura 9. Obsérvense los valores aumentados de la métrica de correlación del bucle flujo-volumen, así como la reducción de los defectos de ventilación tras el tratamiento.

Figura 9: Mediciones 3D de PREFUL en una paciente de 43 años con fibrosis quística (FQ). Mapas representativos de marcadores de ventilación en las mediciones 3D de PREFUL antes del tratamiento (izquierda) y después del tratamiento (derecha) en una paciente de 43 años con FQ. El VDPFVL-CM disminuyó del 18,0 % (línea base) al 3,6 % (post-tratamiento). FEV1 % pred. línea base: 94 %, FEV1 % pred. post-tratamiento: 112 %. LCI línea base: 10,48, LCI post-tratamiento: 9,39. Los mapas correspondientes de respuesta al tratamiento se muestran en la parte inferior. Obsérvese las áreas verdes que indican la ventilación resuelta tras la terapia. Haga clic aquí para ver una versión ampliada de esta figura.
Tabla 1: Lista de parámetros típicos de RM utilizados para la adquisición del localizador y de PREFUL 3D en un escáner de 1,5 T. (A) Localizador y (B) adquisición de PREFUL 3D. Haga clic aquí para descargar esta tabla.
Tabla 2: Parámetros propuestos de RM para la adquisición 3D PREFUL en un escáner de 3T. Haga clic aquí para descargar esta tabla.
Tabla 3: Informe estadístico ejemplar de los parámetros 3D PREFUL en un voluntario sano (hombre de 35 años) Haga clic aquí para descargar esta tabla.