Extending the field of view in modulation-based X-ray phase microtomography

Título: Ampliación del campo de visión en la microtomografía de rayos X de fase basada en modulación

1. El Problema

Las técnicas de imagen de contraste de fase basadas en modulación (como la imagen basada en patrones estructurados o speckle) permiten obtener información cuantitativa sobre la densidad electrónica de los materiales, superando la limitación de la suposición de "material único" necesaria en otros métodos. Sin embargo, estas técnicas enfrentan desafíos críticos al intentar escanear objetos grandes (del orden de centímetros):

• Inestabilidad del haz: Requieren que el perfil del haz de rayos X permanezca constante durante todo el tiempo de medición. En las fuentes de sincrotrón, el perfil varía con el tiempo debido al movimiento de componentes o al control de la órbita del haz de electrones.
• Limitación del Campo de Visión (FOV): Las fuentes de sincrotrón de cuarta generación tienen perfiles de haz más pequeños, reduciendo el FOV. Para escanear objetos grandes, es necesario realizar múltiples escaneos y unirlos (stitching).
• Fallas en la corrección convencional: Los métodos actuales de corrección de campo plano (flat-field), como el enfoque de eigenflat (basado en componentes principales), requieren una región libre de muestra en la imagen de referencia para calibrar el desplazamiento del modulador. En los escaneos centrales de un objeto grande, el campo de visión está completamente ocupado por la muestra, haciendo imposible aplicar estas correcciones y dejando artefactos (patrones del modulador) que degradan la calidad y la precisión cuantitativa.
• Distorsiones no rígidas: La unión de múltiples escaneos suele asumir transformaciones rígidas (traslación/rotación), pero las variaciones temporales del haz y las distorsiones espaciales del detector requieren correcciones deformables más complejas.

2. Metodología

Los autores proponen una nueva técnica de procesamiento de imágenes que combina optimización de eigenflats con registro de imágenes deformables para permitir escaneos cuantitativos de alta resolución de objetos grandes. El flujo de trabajo se divide en las siguientes etapas:

• Adquisición de Escaneos Superpuestos: Se realizan múltiples tomografías parciales (ej. izquierda, central, derecha) con un pequeño solapamiento (aprox. 10% del FOV). Los escaneos exteriores tienen regiones libres de muestra, pero el central no.
• Registro Deformable (Deformable Image Registration):
  • Se utiliza un modelo de deformación libre (FFD) con interpolación B-Spline para alinear las imágenes de los escaneos exteriores con el escaneo central.
  • Esto corrige distorsiones no rígidas causadas por la variación del perfil del haz y la distorsión del detector, asegurando que las características de la muestra coincidan perfectamente en las zonas de solapamiento.
• Corrección de Campo Plano Dinámica Adaptada (Adapted Dynamic Flat-Field Correction):
  • Para el escaneo central (sin región libre de muestra), se plantea un problema de optimización inversa.
  • Se generan imágenes sintéticas de campo plano combinando componentes eigenflat con pesos variables.
  • Función de Costo: Se minimiza la Variación Total (Total Variation - TV) de la diferencia entre las imágenes de fase diferencial del escaneo central y el escaneo exterior (ya corregido y deformado) dentro de la zona de solapamiento.
  • La lógica es que, si la corrección es perfecta, la diferencia entre las imágenes en la zona de solapamiento debe ser homogénea (ruido aleatorio), eliminando cualquier patrón residual del modulador.
• Unión (Stitching) y Reconstrucción:
  • Una vez corregidos los escaneos, se unen mediante una suma ponderada por distancia en las zonas de solapamiento.
  • Las imágenes de fase diferencial unidas se integran (método de Fourier 2D) y se reconstruyen en un mapa 3D de densidad electrónica utilizando el algoritmo FDK (Feldkamp-Davis-Kress).
• Enfoque en Cascada: El método es escalable a n escaneos, propagando la información de referencia desde los escaneos exteriores hacia el interior.

3. Contribuciones Clave

• Eliminación de la restricción de región libre de muestra: Permiten aplicar correcciones de campo plano dinámicas a escaneos donde la muestra ocupa todo el campo de visión, algo imposible con métodos anteriores.
• Corrección de distorsiones complejas: Integran el registro deformable para manejar variaciones no rígidas del haz y del detector, mejorando la precisión de la unión de imágenes.
• Cuantificación robusta: Logran recuperar mapas de densidad electrónica cuantitativos para objetos de varios centímetros sin requerir estabilidad extrema del haz a lo largo del tiempo.
• Versatilidad del modulador: El método se valida tanto con iluminadores de matriz Talbot (rejillas) como con moduladores de papel de lija (sandpaper), demostrando su aplicabilidad general en técnicas de modulación.

4. Resultados

• Experimento: Se aplicó el método a un cerebro de rata de 15 mm de diámetro en la línea de luz P05 (PETRA III, DESY), que tiene un FOV horizontal limitado de 6 mm.
• Calidad de Imagen:
  • Se logró un mapa de densidad electrónica de todo el cerebro con resolución de múltiples micrómetros.
  • Reducción de Artefactos: Comparado con la corrección de campo plano promedio convencional, el método propuesto redujo la desviación estándar en un 8% y la variación total (TV) en un 22% en las zonas de prueba, eliminando casi por completo los artefactos de rejilla visibles.
  • Resolución: La resolución espacial se determinó en 10.9 µm (criterio de medio bit) y 12.4 µm (criterio de bit completo), homogénea en toda la imagen.
  • Uniformidad: No se observaron artefactos de anillos (típicos de errores de corrección) ni discontinuidades visibles en las uniones de los tres escaneos.
• Validación Adicional: Se demostró la viabilidad del enfoque en cascada para múltiples escaneos y su aplicación con papel de lija como modulador, reduciendo el ruido de fondo en un 33%.

5. Significado e Impacto

• Avance para Fuentes de Cuarta Generación: Esta técnica es crucial para aprovechar las nuevas fuentes de sincrotrón (como el ESRF-EBS o PETRA IV) que ofrecen haces más pequeños y brillantes, permitiendo escanear objetos grandes con alta resolución sin perder la capacidad cuantitativa.
• Aplicaciones Multidisciplinarias: Abre nuevas posibilidades en:
  • Biomedicina: Estudio de estructuras biológicas complejas (órganos completos) sin necesidad de cortes físicos, permitiendo correlaciones cuantitativas entre diferentes líneas de luz.
  • Ciencia de Materiales: Análisis de materiales heterogéneos y grandes.
  • Paleontología: Visualización de fósiles grandes con detalles internos.
• Cuantificación Real: Al eliminar la necesidad de la suposición de material único, permite obtener valores reales de densidad electrónica, facilitando la comparación directa de resultados entre diferentes experimentos e instalaciones.

En resumen, el artículo presenta una solución innovadora que supera las barreras de escala y estabilidad en la microtomografía de rayos X de fase, permitiendo la obtención de datos cuantitativos de alta resolución en objetos de tamaño centimétrico.