SECTION 01一条时间记录对应许多可能的位置
设探头位于 rᵢ,候选像素位于 r;均匀、无明显折射的介质中,传播时间是两点距离除以声速。在三维空间,固定时间对应以探头为中心的球面;二维教学模型中则是一条圆。单个探头不能知道球面上的哪一点产生了信号,多探头约束交汇后才可能定位。这里传播的是从吸收体到探头的一程声波,不能套用超声脉冲回波中的往返距离。
原始声压常有正负双极结构。若直接把包络强度沿圆弧相加,得到的是一种直观定位图,并不自动成为初始压力的严格逆解。保留原始符号、时间轴和处理顺序,能帮助解释为何看似相近的实现给出不同的目标宽度与背景。
SECTION 02通用反投影包含的不只是延时叠加
Xu 与 Wang 的通用反投影研究针对平面、球面和圆柱探测几何,结合经过时间处理的信号与立体角权重恢复初始压力。立体角描述探测面元从候选位置看过去覆盖多少方向;空间位置变化后,同一个面元的贡献也会变化。因此,把所有通道以同样权重相加,通常不能替代与几何一致的公式。(来源 1)
学习时可以并排实现一个最简延时叠加与一个有明确推导来源的反投影版本,但必须给它们不同名称。还应注明采用二维还是三维传播模型、探测面是否满足推导条件、采样面是否截断。所谓通用并不意味着任意介质和任意不完整数据都能精确恢复。
SECTION 03一个可以手算的声速偏差例子
以下数字是教学假设,并非设备实测。若目标到探头的真实距离为 15 毫米,真实声速为 1500 米每秒,则到达时间为 10 微秒。重建误用 1540 米每秒,会把这个时刻解释为 15.4 毫米,单条传播路径上出现 0.4 毫米的距离偏差。多探头情况下,各条路径的误差方向不同,表现可能是失焦、环纹或形变,而不只是整幅图平移。
若采样率设为 40 MHz,时间间隔为 25 纳秒,对应一程距离间隔 37.5 微米。这只是采样距离,不是空间分辨率;实际分辨率还受带宽、孔径、角度覆盖与插值影响。计算采样索引后应使用明确定义的插值规则,并把超出记录时间的索引处理为无贡献,不能静默折返。
SECTION 04用四组变化拆开算法与采集误差
建议从一个居中的圆形初始压力目标开始,固定显示范围,依次只改变声速、探头坐标、通道数和时间零点。比较目标质心、横向与轴向剖面、背景能量,避免同时改动多个因素后只挑最好看的图。再把目标移到边缘,检查几何权重在整个视场内是否一致。
当需要考虑异质声速或更一般的探测几何时,时间反演是另一种重建路线。k-Wave 官方线阵示例展示了这种方法与单步 FFT 方法的区别;这提供了比较工具,而不是证明某一种方法在所有数据上更好。(来源 2)
- 保存未归一化的重建结果,再统一生成展示图。
- 记录滤波器频段、插值方法、时间零点和声速。
- 模拟与重建使用不同网格,并加入明确水平的噪声,避免结果过于理想。
SECTION 05什么时候不应继续靠调色解决问题
若改变声速后目标明显收缩,首先检查传播模型;若缺失边缘随阵列旋转而变化,首先检查视角覆盖;若所有目标幅值随通道位置变化,检查灵敏度和空间加权。颜色范围只能改变可见性,不能增加缺失的测量。对重建基线的要求,应是参数明确、误差可解释、重复运行可比较。
k-Wave 的圆弧传感器示例专门讨论了避免“逆犯罪”:正向模拟和反向重建使用完全相同的离散条件,可能掩盖真实误差。将这一提醒加入自己的反投影练习,比追求一张无噪声的完美图更有学习价值。(来源 3)
带走这三个要点
- 延时只是反投影的一部分,传播模型和空间权重同样重要。
- 采样距离间隔不能直接当作成像分辨率。
- 先用可控算例定位声速、坐标和零点误差,再比较复杂算法。
参考资料与继续阅读
以下资料用于支持或进一步核对本文涉及的概念。本文的教学算例与组织方式由本站编写;具体研究结论应以原文的实验条件和限制为准。
编辑说明:初始文章由生成式 AI 根据所列公开资料辅助撰写,尚未经过领域专家审稿。示意与假设算例不构成研究数据,涉及健康的内容不替代专业医疗判断。阅读完整说明 ↗