UltraScatter: Simulación Basada en Rayos de la Dispersión de Ultrasonidos

Información Básica

• ID del Artículo: 2510.10612
• Título: UltraScatter: Ray-Based Simulation of Ultrasound Scattering
• Autores: Felix Duelmer, Mohammad Farid Azampour, Nassir Navab
• Clasificación: physics.med-ph cs.CV
• Fecha de Publicación: 12 de octubre de 2025 (preimpresión en arXiv)
• Enlace del Artículo: https://arxiv.org/abs/2510.10612

Resumen

Los métodos tradicionales de simulación de ultrasonidos logran alta precisión resolviendo numéricamente la ecuación de ondas, pero con un costo computacional enorme. Los métodos alternativos rápidos basados en convolución de respuestas de impulsos precomputados siguen siendo relativamente lentos, requiriendo típicamente varios minutos para generar una imagen B-mode completa. Este artículo presenta UltraScatter, un marco probabilístico de trazado de rayos que modela eficientemente y de manera realista la dispersión de ultrasonidos. El tejido se representa como campos volumétricos de probabilidad de dispersión y amplitud de dispersión, con interacciones de rayos simuladas mediante delta tracking de vuelo libre. Los rayos dispersados se trazan hasta el sensor, integrando información de fase mediante un modelo lineal de tiempo de vuelo. Combinando formación de imagen de ondas planas y conformación de haces, la arquitectura de trazado de rayos paralelizada genera imágenes B-mode en segundos. La validación con datos de fantomas demuestra patrones de moteado y de inclusiones realistas, posicionando a UltraScatter como una alternativa escalable a los métodos de ecuación de ondas.

Antecedentes de Investigación y Motivación

Definición del Problema

La simulación de ultrasonidos es crucial para mejorar algoritmos de reconstrucción, optimizar diseños de sensores y entrenar modelos de aprendizaje automático en condiciones controladas. Sin embargo, los métodos existentes enfrentan un compromiso entre eficiencia computacional y realismo físico.

Limitaciones de Métodos Existentes

1. Solucionadores de ecuación de ondas (como FDTD, métodos pseudoespectrales k-space): Alta precisión pero costo computacional enorme, difícil de escalar
2. Métodos rápidos como Field II: Mediante convolución de respuestas espaciales de impulsos precomputados, aún requieren varios minutos para generar una única imagen B-mode
3. SIMUS: Basado en operaciones de suma retrasada de distancias geométricas de ida y vuelta, eficiencia limitada

Motivación de la Investigación

Se necesita un método que combine realismo físico con cálculo rápido, haciendo posible la simulación de ultrasonidos en tiempo real o interactiva. Las técnicas de renderizado físico (PBR) de gráficos por computadora proporcionan inspiración, particularmente algoritmos de trazado de rayos altamente optimizados acelerados por CUDA.

Contribuciones Principales

1. Marco modular de alto rendimiento: Modelado de atenuación, absorción y dispersión múltiple en medios participantes
2. Estrategia de muestreo del transmisor: Conecta cada interacción de escena con todos los elementos del sensor
3. Cadena completa de conformación de haces de transmisión-recepción: Convierte directamente ecos simulados en imágenes B-mode
4. Mejora de rendimiento significativa: Aproximadamente 70 veces más rápido que implementaciones de código de dominio de frecuencia tradicional

Detalles del Método

Definición de la Tarea

Comenzando desde un mapa etiquetado, se asignan características de dispersión para cada categoría de tejido, utilizando un esquema de trazado de rayos Monte Carlo para modelar emisión de ondas de presión, dispersión y atenuación, escribiendo ecos que regresan al sensor en búferes RF de elementos específicos, y luego generando la imagen B-mode final mediante un conformador de haces digital convencional.

Arquitectura del Modelo

Modelo Matemático Principal

La señal de presión que llega al elemento del sensor e en el tiempo t se define como:

$P(e, t) = \int_\Omega \int_A P_i(x, t, \omega_i) f_d(\omega_i) d\omega da$

donde $P_i(x, t, \omega_i)$ es la presión incidente proveniente de la ubicación de escena x, dirección $\omega_i$, y $f_d$ es la función de directividad.

La presión dispersada se modela como:

$P_{scattered}(x, t, \omega_o) = \int_\Omega a(x) p(\omega_i, \omega_o)P_i(x, t, \omega_i) d\omega_i$

donde $a(x)$ representa la amplitud de dispersión y $p(\omega_i, \omega_o)$ es la función de fase que controla la redistribución angular.

Algoritmo de Trazado de Rayos

1. Muestreo de Distancia de Vuelo Libre: $s = s_{min} - \frac{1}{\mu} \ln(1-\xi)$

donde $\mu$ es el valor máximo de extinción en todo el medio, $\xi \sim U(0,1)$ es un número aleatorio.

2. Determinación de Interacción: $\xi' < \frac{\sigma_t(x)}{\mu} \Rightarrow \text{interacción real, de lo contrario interacción nula}$

3. Estrategia de Muestreo del Transmisor: En cada evento de dispersión, se muestrea para todos los elementos del sensor, emitiendo un conjunto de rayos secundarios, cada uno dirigido a un elemento del sensor.

4. Manejo de Aproximación de Campo Lejano: Se subdivide cada elemento del sensor en $\nu$ subelementos idénticos: $\nu = \left\lceil \frac{2b}{\lambda_{min}} \right\rceil$

donde $b$ es el ancho azimutal de un elemento de sensor individual, $\lambda_{min}$ es la longitud de onda mínima en el pulso de transmisión.

Puntos de Innovación Técnica

1. Delta Tracking de Woodcock: Adopta delta tracking de vuelo libre con interacciones nulas, mejorando la eficiencia computacional
2. Muestreo de Transmisor de Todos los Elementos: Comparado con métodos selectivos, produce imágenes más claras y bordes de estructura más nítidos
3. Paralelización en GPU: Basada en el marco Mitsuba 3, utilizando NVIDIA OptiX para computación paralela eficiente
4. Integración de Información de Fase: Mantiene consistencia de fase mediante un modelo lineal de tiempo de vuelo

Configuración Experimental

Conjunto de Datos

Se utiliza el fantoma universal CIRS (Modelo 054GS) para validación, reconstruyendo su estructura interna según especificaciones del fabricante y ajustando parámetros faltantes hasta que la referencia de simulación coincida con el B-mode clínico.

Métodos de Comparación

1. Escáner Clínico: Siemens Acuson Juniper equipado con matriz lineal 12L3 (192 elementos, 2.9-11.5 MHz)
2. SIMUS: Simulador armónico temporal de dominio de frecuencia, basado en aproximaciones de campo lejano y paraxial

Detalles de Implementación

• Frecuencia central: 6.2 MHz
• Profundidad de imagen: 10 cm
• Rango dinámico: 60 dB
• Frecuencia de muestreo: 25 MHz
• Número F: 1.0
• 100,000 rayos emitidos por subelemento
• Hardware: CPU Intel Core i7-12700 + GPU NVIDIA RTX 4070 Ti

Resultados Experimentales

Resultados Principales

Comparación de Rendimiento Computacional:

• SIMUS: 634 ± 3 segundos
• UltraScatter: 9.3 ± 0.8 segundos
• Mejora de velocidad: aproximadamente 70 veces

Evaluación de Calidad de Imagen:

1. Estructura Anatómica General: Los tres métodos (clínico, UltraScatter, SIMUS) comparten estructura anatómica general similar
2. Estadísticas de Moteado de Fondo: UltraScatter produce estadísticas de moteado similares a los escaneos reales
3. Sombra Distal: Refleja correctamente los efectos de sombra del escaneo real

Análisis de Casos

Características de UltraScatter:

• Existe sobreiluminación en la región directamente debajo de la sonda, resultado de la sobreestimación residual de ganancia de campo cercano
• Produce bordes de inclusión más nítidos, atribuido a la estrategia de muestreo de transmisor de todos los elementos
• Mayor claridad en la resolución del cilindro profundo comparado con SIMUS

Características de SIMUS:

• Reproduce más fielmente la atenuación de intensidad axial del sistema clínico
• Definición de bordes de inclusión ligeramente menos clara

Trabajo Relacionado

Métodos Tradicionales de Simulación de Ultrasonidos

1. Solucionadores de Onda Completa: FDTD, métodos pseudoespectrales k-space - Exactos pero computacionalmente costosos
2. Métodos de Convolución: Field II - Basado en respuestas de impulsos precomputadas
3. Métodos Geométricos: SIMUS - Basado en operaciones de suma retrasada

Trazado de Rayos de Convolución (CRT)

Alternativa rápida que ha surgido recientemente, aproximando la propagación de ondas acústicas mediante el trazado de grandes conjuntos de rayos, con reflexión y refracción en límites macroscópicos, seguido de modelado de dispersión de escala fina mediante convolución de funciones de dispersión puntual separables.

Técnicas de Renderizado Físico

Técnicas PBR en gráficos por computadora, incluyendo marcha de rayos, mapeo de fotones y delta tracking, proporcionan implementaciones altamente optimizadas para manejar propagación de luz en medios heterogéneos.

Conclusiones y Discusión

Conclusiones Principales

UltraScatter implementa exitosamente simulación de ultrasonidos mediante trazado de rayos, reemplazando solucionadores de dominio de frecuencia con muestreo de trayectorias Monte Carlo, renderizando vistas típicas de fantomas CIRS en aproximadamente 9 segundos en hardware idéntico, mientras que SIMUS requiere más de 10 minutos, logrando una mejora de rendimiento significativa.

Limitaciones

1. Estadísticas de Moteado: Impulsadas por muestreo aleatorio, dos ejecuciones independientes producen patrones diferentes, la textura se descorrelaciona cuando la sonda se mueve lateralmente
2. Configuración de Rayos: El prototipo actual emite rayos desde una única línea axial a una cuadrícula cúbica, configuraciones más realistas requieren apertura elevacional completa entre sensores
3. Coherencia Temporal: Requiere introducir secuencias aleatorias correlacionadas o mapas de semillas compartidas para mantener estabilidad de moteado con movimientos pequeños de sonda

Direcciones Futuras

1. Enfoque Elevacional: Mediante retrasos de recepción dinámicos o emisión de rayos adicionales
2. Modelado Físico Mejorado: Asignación de densidad, velocidad acústica y atenuación independientes para cada vóxel
3. Fenómenos No Lineales: Modelado de refracción macroscópica, aberración de fase e incluso generación armónica
4. Soporte de Escenas Dinámicas: Mejora de consistencia de moteado y resolución elevacional

Evaluación Profunda

Fortalezas

1. Mejora de Rendimiento Significativa: 70 veces más rápido hace posible la simulación en tiempo real
2. Realismo Físico: Produce calidad de imagen y patrones de moteado similares a escaneos clínicos
3. Arquitectura Modular: Diseño basado en Mitsuba 3 facilita extensión y mejora
4. Paralelización en GPU: Aprovecha plenamente la capacidad de computación paralela del hardware moderno
5. Pipeline Completo: Solución de extremo a extremo desde trazado de rayos hasta generación de imagen B-mode

Deficiencias

1. Efectos de Campo Cercano: Existe problema de sobreestimación de ganancia de campo cercano
2. Consistencia de Moteado: Falta coherencia temporal, inadecuado para escenas dinámicas
3. Limitación Dimensional: Actualmente solo soporta imagen 2D, carece de procesamiento elevacional 3D completo
4. Alcance de Validación: Validado solo en un único fantoma, carece de datos de validación más amplios

Impacto

Contribuciones Académicas:

• Primera aplicación exitosa de técnicas de trazado de rayos de gráficos por computadora a simulación de ultrasonidos
• Proporciona nueva ruta tecnológica para el campo de simulación de ultrasonidos

Valor Práctico:

• Hace posible simulación de ultrasonidos en tiempo real, beneficiando educación médica y entrenamiento
• Proporciona herramienta rápida de generación de datos para modelos de aprendizaje automático
• Soporta iteración rápida en optimización de diseño de sensores

Reproducibilidad:

• Basado en marco Mitsuba 3 de código abierto
• Proporciona detalles de implementación y configuración de parámetros detallados
• La disponibilidad de código y datos requiere confirmación adicional

Escenarios Aplicables

1. Educación Médica: Sistemas de entrenamiento de simulación de ultrasonidos en tiempo real
2. Desarrollo de Algoritmos: Generación rápida de datos de entrenamiento para modelos de aprendizaje profundo
3. Diseño de Dispositivos: Verificación de prototipo rápido para optimización de parámetros de sensores
4. Aplicaciones de Investigación: Proyectos de investigación científica que requieren grandes cantidades de datos de simulación

Referencias

El artículo cita 19 referencias importantes que cubren:

• Métodos tradicionales de simulación de ultrasonidos (FDTD, k-Wave, Field II, SIMUS)
• Trabajo relacionado con simulación de ultrasonidos mediante trazado de rayos
• Técnicas de renderizado de gráficos por computadora (Mitsuba 3, OptiX)
• Métodos Monte Carlo y algoritmos de delta tracking

Resumen: UltraScatter representa un avance tecnológico importante en el campo de la simulación de ultrasonidos, integrando exitosamente técnicas avanzadas de trazado de rayos de gráficos por computadora en simulación médica de ultrasonidos, logrando un buen equilibrio entre eficiencia computacional y realismo físico. A pesar de algunas limitaciones, su mejora de rendimiento significativa y diseño modular establecen una base sólida para investigación futura.