Improved 3D Radial Phyllotaxis Trajectories for Uniform Density Distribution of Readout Directions and Sequential Binning

Resumen Técnico: Trayectorias de Filotaxia Radial 3D Mejoradas para una Distribución de Densidad Uniforme y Binning Secuencial

Planteamiento del Problema
Las trayectorias radiales tridimensionales, específicamente la filotaxia espiral, se utilizan ampliamente en la RM dinámica y con resolución de movimiento debido a su robustez ante el movimiento y sus capacidades de autonavigación. Sin embargo, la trayectoria de filotaxia espiral estándar sufre dos limitaciones primarias:

1. Densidad de Muestreo No Uniforme: La trayectoria estándar no logra una distribución esférica uniforme de la densidad de muestreo (SDSD, por sus siglas en inglés). Tiende a agrupar las direcciones de lectura cerca de la región ecuatorial, lo que provoca un muestreo de k-space anisotrópico, variaciones dependientes de la dirección en el campo de visión efectivo y la resolución espacial, y posibles artefactos de aliasing.
2. Artefactos de Binning Secuencial: Aunque la cobertura global del k-space es eficiente, muchas aplicaciones (por ejemplo, imágenes dinámicas, RM funcional, reconstrucción con resolución de movimiento) requieren la reconstrucción de imágenes a partir de subconjuntos temporalmente contiguos de lecturas (binning secuencial). En la filotaxia estándar, los entrelazados (interleaves) consecutivos exhiben un agrupamiento sistemático en los ángulos polares. Cuando estos subconjuntos consecutivos se agrupan (binning), crean huecos circulares en el k-space, lo que se manifiesta como artefactos de ringing (anillamiento) e inconsistencias temporales en las imágenes reconstruidas.

Metodología
Los autores proponen dos avances metodológicos para abordar estas limitaciones, evaluados in vivo en 10 voluntarios sanos utilizando un escáner de RM de 3T con dos secuencias de eco de gradiente (GRE ponderada en T1 y LIBRE).

1. Filotaxia Uniforme (UPhy):

   • Objetivo: Lograr una distribución de densidad de direcciones de lectura uniforme (SDSD uniforme).
   • Mecanismo: UPhy redefine la formulación del ángulo polar de la filotaxia estándar. Mientras que el modelo estándar utiliza $\theta_n = \frac{\pi}{2}\sqrt{\frac{n}{N}}$, UPhy utiliza $\theta_n = \arccos(1 - \frac{n}{N})$.
   • Efecto: Esta modificación asegura una distribución uniforme de puntos sobre la esfera, eliminando el agrupamiento ecuatorial observado en el modelo estándar, mientras mantiene el orden azimutal (ángulo áureo) original.

2. Filotaxia Flexible (FlexiPhy):

   • Objetivo: Mejorar la robustez ante un binning secuencial retrospectivo arbitrario.
   • Mecanismo: FlexiPhy se basa en UPhy pero introduce una permutación aleatoria ($\sigma$) de los entrelazados. Esto desacopla el orden del ángulo polar de la progresión del ángulo áureo azimutal. Específicamente, el ángulo polar para el $j$-ésimo punto en el $k$-ésimo entrelazado se determina mediante un índice permutado aleatoriamente $\sigma(k)$, en lugar del índice secuencial $k$.
   • Efecto: Esta aleatorización rompe el agrupamiento sistemático de los ángulos polares dentro de los bins temporales consecutivos. En consecuencia, los subconjuntos secuenciales de lecturas mantienen un patrón de muestreo bien distribuido, evitando la formación de huecos concéntricos en el k-space y los asociados artefactos de ringing.

Contribuciones Clave

• Uniformidad Analítica y Empírica: El artículo demuestra que UPhy proporciona una SDSD uniformemente probada analíticamente, corrigiendo el sesgo de densidad de la filotaxia estándar.
• Desacoplamiento para Binning: La introducción de FlexiPhy ofrece una estrategia novedosa para desacoplar el orden polar y el azimutal mediante la aleatorización, dirigida específicamente a mitigar la degradación de las propiedades de muestreo en subconjuntos temporales contiguos.
• Validación In Vivo: Las trayectorias fueron probadas utilizando dos secuencias distintas (GRE estándar y LIBRE con supresión de grasa) en escenarios de imagen cerebral estática e imagen ocular dinámica.

Resultados

• Imagen Estática: La inspección visual de las reconstrucciones de secuencias temporales consecutivas reveló que tanto la filotaxia estándar como la UPhy exhibieron visibles artefactos de ringing e inconsistencias temporales. FlexiPhy redujo sustancialmente estos artefactos, produciendo imágenes más estables temporalmente.
• Imagen Dinámica: En la imagen ocular dinámica, las trayectorias estándar y UPhy produjeron artefactos de ringing estructurados que se propagaron hacia la región orbital, oscureciendo detalles anatómicos. FlexiPhy redujo significativamente estos artefactos, mejorando la visualización del ojo y las estructuras circundantes.
• Métricas Cuantitativas:
  • SSIM y Error L2: FlexiPhy demostró valores de Índice de Similitud Estructural (SSIM) significativamente más altos y errores relativos L2 más bajos en comparación con la filotaxia estándar y UPhy para ambas secuencias GRE y LIBRE (Bonferroni-corrected $p < 0.05$).
  • UPhy vs. Estándar: Aunque UPhy mejoró la uniformidad angular teórica, no redujo significativamente los artefactos de ringing en el binning secuencial en comparación con la trayectoria estándar, excepto por una ligera mejora en el error L2 para la secuencia LIBRE. Esto confirma que la mejora de la SDSD global por sí sola es insuficiente para resolver los artefactos de binning secuencial.

Significancia y Reivindicaciones
El artículo concluye que FlexiPhy permite reconstrucciones de binning secuencial más fiables al reducir los artefactos de ringing inducidos por la trayectoria y las inconsistencias temporales. Su construcción aleatorizada no está ligada a una estrategia de binning específica, lo que la hace ampliamente compatible con los enfoques de binning retrospectivo utilizados en la RM dinámica y con resolución de movimiento.

Los autores enfatizan que, si bien UPhy corrige la distribución de densidad global, es el desacoplamiento del orden polar y azimutal en FlexiPhy lo que es crítico para mantener la calidad del muestreo dentro de los subconjuntos temporales. El trabajo sugiere que FlexiPhy preserva las ventajas prácticas de la filotaxia espiral (cobertura eficiente del k-space, transiciones de gradiente suaves) ofreciendo al mismo tiempo una robustez superior para aplicaciones que requieren resolución temporal o clasificación retrospectiva. El estudio reconoce limitaciones, señalando que trabajos futuros deberían investigar el impacto en tareas posteriores como la precisión de la estimación de movimiento y explorar la combinación de estas trayectorias con técnicas de adquisición de polo a polo.