UltraScatter: Ray-Based Simulation of Ultrasound Scattering

基本信息

• 论文ID: 2510.10612
• 标题: UltraScatter: Ray-Based Simulation of Ultrasound Scattering
• 作者: Felix Duelmer, Mohammad Farid Azampour, Nassir Navab
• 分类: physics.med-ph cs.CV
• 发表时间: 2025年10月12日 (arXiv预印本)
• 论文链接: https://arxiv.org/abs/2510.10612

摘要

传统的超声仿真方法通过数值求解波动方程实现高精度，但计算成本巨大。基于预计算脉冲响应卷积的快速替代方法仍然相对缓慢，通常需要几分钟才能生成完整的B模式图像。本文提出UltraScatter，一个概率射线追踪框架，能够高效且逼真地建模超声散射。组织被表示为散射概率和散射幅度的体积场，射线相互作用通过自由飞行delta追踪进行仿真。散射射线被追踪到传感器，相位信息通过线性飞行时间模型整合。结合平面波成像和波束成形，并行化射线追踪架构在几秒钟内生成B模式图像。通过体模数据验证显示出逼真的斑点和包含物模式，将UltraScatter定位为波动方程方法的可扩展替代方案。

研究背景与动机

问题定义

超声仿真在增强重建算法、优化传感器设计和在受控条件下训练机器学习模型方面至关重要。然而，现有方法面临计算效率与物理真实性的权衡问题。

现有方法局限性

1. 波动方程求解器（如FDTD、k-space伪谱方法）：精度高但计算成本巨大，难以扩展
2. Field II等快速方法：通过卷积预计算空间脉冲响应，仍需数分钟生成单个B模式图像
3. SIMUS：基于几何往返距离的延迟求和操作，效率有限

研究动机

需要一种结合物理真实性与快速计算的方法，使实时或交互式超声仿真成为可能。计算机图形学中的物理渲染(PBR)技术为此提供了启发，特别是高度优化的CUDA加速射线追踪算法。

核心贡献

1. 模块化高性能框架：建模参与介质中的衰减、吸收和多重散射
2. 发射器采样策略：将每个场景相互作用与所有传感器元素连接
3. 完整的发射-接收波束成形链：直接将仿真回波转换为B模式图像
4. 显著性能提升：相比传统频域代码实现约70倍速度提升

方法详解

任务定义

从标签图开始，为每个组织类别分配散射特性，使用蒙特卡罗射线追踪方案建模压力波发射、散射和衰减，将返回到传感器的回波写入特定元素的RF缓冲区，然后通过常规数字波束成形器生成最终B模式图像。

模型架构

核心数学模型

到达传感器元素e在时间t的压力信号P定义为：

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

其中$P_i(x, t, \omega_i)$是从场景位置x、方向$\omega_i$传来的入射压力，$f_d$是方向性函数。

散射压力建模为：

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

其中$a(x)$表示散射幅度，$p(\omega_i, \omega_o)$是控制角度重分布的相位函数。

射线追踪算法

1. 自由飞行距离采样： $s = s_{min} - \frac{1}{\mu} \ln(1-\xi)$

其中$\mu$是整个介质中的最大消光值，$\xi \sim U(0,1)$是随机数。

2. 相互作用判定： $\xi' < \frac{\sigma_t(x)}{\mu} \Rightarrow \text{真实相互作用，否则为空相互作用}$

3. 发射器采样策略：每次散射事件对所有传感器元素进行采样，发射一组次级射线，每个指向一个传感器元素。

4. 远场近似处理：将每个传感器元素细分为$\nu$个相同的子元素： $\nu = \left\lceil \frac{2b}{\lambda_{min}} \right\rceil$

其中$b$是单个传感器元素的方位宽度，$\lambda_{min}$是传输脉冲中的最小波长。

技术创新点

1. Woodcock追踪：采用带空相互作用的自由飞行delta追踪，提高计算效率
2. 全元素发射器采样：相比选择性方法，产生更清晰的图像和更锐利的结构边界
3. GPU并行化：基于Mitsuba 3框架，利用NVIDIA OptiX实现高效并行计算
4. 相位信息整合：通过线性飞行时间模型保持相位一致性

实验设置

数据集

使用CIRS通用体模(Model 054GS)进行验证，根据制造商规格重建其内部结构，并调整缺失参数直到仿真参考与临床B模式匹配。

对比方法

1. 临床扫描仪：Siemens Acuson Juniper配备12L3线性阵列（192个元素，2.9-11.5 MHz）
2. SIMUS：频域时间谐波仿真器，基于远场和近轴假设

实现细节

• 中心频率：6.2 MHz
• 成像深度：10 cm
• 动态范围：60 dB
• 采样率：25 MHz
• F数：1.0
• 每个子元素发射10万条射线
• 硬件：Intel Core i7-12700 CPU + NVIDIA RTX 4070 Ti GPU

实验结果

主要结果

计算性能对比：

• SIMUS：634 ± 3秒
• UltraScatter：9.3 ± 0.8秒
• 速度提升：约70倍

图像质量评估：

1. 整体解剖结构：三种方法（临床、UltraScatter、SIMUS）共享相似的整体解剖结构
2. 背景斑点统计：UltraScatter产生与真实扫描相似的斑点统计
3. 远端阴影：正确反映真实扫描的阴影效应

案例分析

UltraScatter特点：

• 在探头正下方区域存在过度照明，这是近场增益残余过估计的结果
• 产生更锐利的包含物边界，归因于全元素发射器采样策略
• 对深层圆柱体的分辨率比SIMUS更清晰

SIMUS特点：

• 更忠实地再现临床系统的轴向强度衰减
• 包含物边缘定义略不清晰

相关工作

传统超声仿真方法

1. 全波求解器：FDTD、k-space伪谱方法 - 精确但计算昂贵
2. 卷积方法：Field II - 基于预计算脉冲响应
3. 几何方法：SIMUS - 基于延迟求和操作

卷积射线追踪(CRT)

近期出现的快速替代方案，通过追踪大量射线集合近似声波传播，在宏观边界处反射和折射，然后通过卷积可分离点扩散函数建模细尺度散射。

物理渲染技术

计算机图形学中的PBR技术，包括射线行进、光子映射和delta追踪等，为处理异质介质中的光传播提供了高度优化的实现。

结论与讨论

主要结论

UltraScatter成功实现了射线追踪超声仿真，用蒙特卡罗路径采样替代频域求解器，在相同硬件上渲染典型CIRS体模视图约需9秒，而SIMUS需要超过10分钟，实现了显著的性能提升。

局限性

1. 斑点统计：由随机采样驱动，两次独立运行产生不同模式，探头横向移动时纹理去相关
2. 射线配置：当前原型从单一轴向线发射射线到立方网格，更现实的配置需要跨传感器全elevational孔径
3. 时间相干性：需要引入相关随机序列或共享种子图以在小探头运动中保持斑点稳定

未来方向

1. Elevational聚焦：通过动态接收延迟或投射额外射线实现
2. 增强物理建模：为每个体素分配独立的密度、声速和衰减
3. 非线性现象：建模宏观折射、相位像差甚至谐波生成等非线性现象
4. 动态场景支持：改进斑点一致性和elevational分辨率

深度评价

优点

1. 显著性能提升：70倍速度提升使实时仿真成为可能
2. 物理真实性：产生与临床扫描相似的图像质量和斑点模式
3. 模块化架构：基于Mitsuba 3的设计便于扩展和改进
4. GPU并行化：充分利用现代硬件的并行计算能力
5. 完整pipeline：从射线追踪到B模式图像生成的端到端解决方案

不足

1. 近场效应：存在近场增益过估计问题
2. 斑点一致性：缺乏时间相干性，不适合动态场景
3. 维度限制：当前仅支持2D成像，缺乏完整的3D elevational处理
4. 验证范围：仅在单一体模上验证，缺乏更广泛的验证数据

影响力

学术贡献：

• 首次将计算机图形学的射线追踪技术成功应用于超声仿真
• 为超声仿真领域提供了新的技术路径

实用价值：

• 使实时超声仿真成为可能，有助于医学教育和训练
• 为机器学习模型提供快速数据生成工具
• 支持传感器设计优化的快速迭代

可复现性：

• 基于开源Mitsuba 3框架
• 提供了详细的实现细节和参数设置
• 代码和数据的可获得性需要进一步确认

适用场景

1. 医学教育：实时超声仿真训练系统
2. 算法开发：快速生成训练数据用于深度学习模型
3. 设备设计：传感器参数优化的快速原型验证
4. 研究应用：需要大量仿真数据的科学研究项目

参考文献

论文引用了19篇重要参考文献，涵盖：

• 传统超声仿真方法（FDTD、k-Wave、Field II、SIMUS）
• 射线追踪超声仿真相关工作
• 计算机图形学渲染技术（Mitsuba 3、OptiX）
• Monte Carlo方法和delta追踪算法

总结：UltraScatter代表了超声仿真领域的重要技术突破，成功将计算机图形学的先进射线追踪技术引入医学超声仿真，实现了计算效率和物理真实性的良好平衡。尽管存在一些局限性，但其显著的性能提升和模块化设计为未来研究奠定了坚实基础。