Análisis de la respuesta de detectores de Rayos X: Un Kit Completo para Radiofísicos

Introducción

Como radiofísicos, evaluar el funcionamiento de los detectores de equipos de rayos X es crucial para garantizar una calidad de imagen óptima y la optimización de la dosis del paciente. Hoy me complace compartir un kit de análisis integral que simplifica los cálculos complejos requeridos para la evaluación de la respuesta del detector según los estándares internacionales.

El Kit de Análisis de Imágenes de Rayos X es una aplicación web construida con Streamlit que proporciona a los radiofísicos una interfaz intuitiva para analizar los detectores a través de múltiples métricas de calidad.

Características y Capacidades Principales

1. Función Respuesta del Detector (MPV vs Kerma)

La característica principal del kit analiza la relación fundamental entre la respuesta del detector y el kerma incidente en aire. Este análisis es esencial para determinar el resto de parámetros esenciales a la hora de caracterizar un detector:

• Componentes del ruido
• Uniformidad de la imagen
• Espectro de potencia de ruido normalizado (NNPS)
• Resolución espacial (MTF)
• Umbral de sensibilidad

Fundamentos Matemáticos:

• Modelo lineal: MPV = a × Kerma + b
• Modelo logarítmico: MPV = a × ln(Kerma) + b
• Modelo polinomial: MPV = a × Kerma² + b × Kerma + c

El análisis determina automáticamente el modelo de mejor ajuste basado en el tipo de datos (archivos RAW vs STD) y proporciona validación estadística completa incluyendo valores R² y análisis de residuos.

2. Análisis de Uniformidad de Imagen

Siguiendo estándares internacionales, el kit realiza evaluaciones de uniformidad tanto globales como locales:

Métricas de Uniformidad Global:

• GU_PV (Uniformidad Global - Valor de Píxel): Desviación máxima de la media del ROI central en todos los ROIs móviles
• GU_SNR (Uniformidad Global - Relación Señal-Ruido): Evaluación de la variación de SNR

Métricas de Uniformidad Local:

• LU_PV (Uniformidad Local - Valor de Píxel): Desviación máxima entre cada ROI y sus 8 vecinos
• LU_SNR (Uniformidad Local - SNR): Evaluación de consistencia local de SNR

Parámetros de Análisis:

• ROI central: 80% del área total de la imagen para mediciones de línea base
• Tamaño de ROI móvil: Ventanas de 30mm × 30mm
• Tamaño de paso: Espaciado de cuadrícula de 15mm
• Conversión automática píxel-milímetro usando espaciado de píxeles DICOM

3. Análisis de Espectro de Potencia de Ruido (NPS)

El módulo NPS implementa la metodología estándar IEC 62220-1-2 para caracterización de ruido:

Características Clave:

• Promediado multi-imagen: Combina múltiples imágenes de campo plano para mejorar estadísticas
• Análisis basado en ROI: ROI central grande (125mm × 125mm) con ventanas de análisis más pequeñas
• Promediado radial: Cálculo de NPS 1D del espectro de potencia 2D
• Cálculo NNPS: Espectro de Potencia de Ruido Normalizado para comparación de detectores

Implementación Técnica:

• Eliminación automática de tendencias para remover variaciones sistemáticas
• Funciones de ventana para reducción de fuga espectral
• Validación estadística con mínimo 4 millones de píxeles totales de ROI

4. Función de Transferencia de Modulación (MTF)

Implementación del método de borde inclinado compatible con IEC 62220-1-1:2015:

Pipeline de Análisis:

1. Detección de Borde: Detección automática del ángulo del borde (óptimo 3-5° desde vertical/horizontal)
2. Extracción ESF: Función de Dispersión del Borde perpendicular a la dirección del borde
3. Cálculo LSF: Función de Dispersión de Línea vía diferenciación
4. Cálculo MTF: Transformada de Fourier con normalización adecuada

Correcciones Clave Aplicadas:

• Corrección de suavizado espectral para diferenciación de elemento finito
• Corrección de escalado del eje de frecuencia (1/cos α) para proyección oblicua
• Calibración automática de frecuencia espacial usando espaciado de píxeles

Métricas de Resolución Espacial:

• MTF50: Frecuencia espacial al 50% de modulación (métrica de resolución primaria)
• MTF10: Frecuencia espacial al 10% de modulación (resolución límite)
• MTF2: Frecuencia espacial al 2% de modulación (evaluación del piso de ruido)

5. Análisis de Contraste Umbral

Análisis avanzado de contraste que requiere datos tanto de MTF como de NPS:

Marco de Cálculo:

• Umbrales de contraste específicos del objeto
• Integración de parámetros geométricos del sistema
• Análisis de frecuencia Nyquist
• Evaluación combinada de rendimiento espacial y de ruido

Arquitectura Técnica

Soporte de Formatos de Archivo

• DICOM: Análisis completo de cabecera con extracción automática de espaciado de píxeles
• RAW/STD: Análisis directo de datos de píxeles con parámetros definidos por el usuario
• Procesamiento multi-archivo: Análisis en lote para mejora estadística

Pipeline de Procesamiento de Datos

• Validación de Calidad: Verificación automática de integridad de imagen
• Calibración: Conversión de píxel a unidad física
• Análisis Estadístico: Propagación completa de errores y estimación de incertidumbre
• Capacidades de Exportación: Exportación de resultados para análisis adicional y reporte

Aplicaciones Prácticas

Programas de Aseguramiento de Calidad

• Monitoreo rutinario del rendimiento del detector
• Pruebas de aceptación de nuevos sistemas de imagen
• Análisis de tendencias para mantenimiento predictivo

Aplicaciones de Investigación

• Estudios de comparación de detectores
• Investigación de optimización de calidad de imagen
• Desarrollo y validación de protocolos

Implementación Clínica

• Estudios de optimización de dosis usando respuesta cuantitativa del detector
• Evaluación de calidad de imagen para diferentes protocolos clínicos
• Reporte de cumplimiento regulatorio

Fundamentos Matemáticos

El kit implementa modelos matemáticos rigurosos basados en principios establecidos de física médica:

Cálculo de Uniformidad:

Término_Uniformidad = |Valor_ROI - Valor_Referencia| / Valor_Referencia × 100%

Normalización NPS:

NNPS(f) = NPS(f) × (Tamaño_Píxel)² / Valor_Medio_Píxel²

Análisis de Borde MTF: El método de borde inclinado sigue el pipeline completo de IEC con correcciones apropiadas para muestreo finito y efectos de ángulo del borde.

Conclusión

Este kit de análisis de imágenes de rayos X representa un avance significativo en software de física médica, proporcionando capacidades completas de análisis de detector en una interfaz web accesible. Al implementar estándares internacionales y mejores prácticas, permite a los físicos médicos realizar evaluaciones de calidad rigurosas de manera eficiente y consistente.

La combinación de análisis de respuesta del detector, evaluación de uniformidad, caracterización de ruido, y medición de resolución espacial proporciona una imagen completa del rendimiento del detector esencial para resultados óptimos de imagen clínica.

Ya sea que estés realizando QA de rutina, investigaciones de investigación, o pruebas de aceptación del sistema, este kit ofrece las capacidades analíticas necesarias para una evaluación exhaustiva del detector en imagen médica moderna.

---

Para documentación técnica detallada y código fuente, visita el repositorio de GitHub. El kit es de código abierto y mantenido activamente para la comunidad de física médica.