SECTION 011. 同一声压可以对应不同材料参数
在简化模型中 p₀(r,λ)=Γ(r)μₐ(r,λ)F(r,λ)。如果 Γ 已知而 F 未知,一个像素的声压仍包含两个相乘的未知量。把 μₐ 增大同时把 F 减小,可能得到同样的乘积。这个代数观察已经足够说明:仅对声学图像做尺度调整,不能自动分离当地照明与吸收。
更复杂的地方在于 F 不是任意的独立图像,它随整个样品的吸收、散射和边界照明变化。2010 年的多光源定量光声研究在已知背景光学性质等条件下讨论吸收扰动重建,明确写出了自己的先验前提(来源 1)。这些前提应随结果一起阅读。
SECTION 022. 一个四倍亮度差的假设例子
假设两个完全相同的吸收靶都具有 Γ=0.1、μₐ=100 m⁻¹,但局部 F 分别是 20 和 5 J/m²。初始声压分别为 200 和 50 Pa。四倍的幅度差来自照明,没有任何浓度差。这是为说明公式而设定的计算,不是某个组织深度的真实衰减曲线。
若人为将深处像素乘四,能让这个例子的幅度相等,却不能据此证明真实数据已经定量。因为真实靶的 F 未必恰好低四倍,接收灵敏度也可能随位置改变。经验补偿可以改善可见性,但应在图注中标明补偿函数,并保留未补偿数据供比较。
SECTION 033. 两层反问题需要分开检查
声学层从探头时间序列推断初始压力分布,需要声速、探测几何与接收响应;光学层再结合光传播模型,把初始压力中的吸收贡献分离出来。即使第一层图像边缘很清楚,第二层仍可能受照明误差影响,因此视觉锐利不能替代参数准确性。
2012 年的 MIPAT 原始研究直接拟合多照明条件下的超声通道数据,同时估计光学参数,并用数值例子研究其性质(来源 2)。这说明声学层与光学层也可以联合求解;它并不意味着未知照明、噪声与模型偏差从此自动消失。
SECTION 044. 多照明和多波长各增加了什么
移动光源位置,会改变同一目标经历的内部光照分布。把这些观测放在一起,有机会约束单次照明下难以区分的参数,但收益依赖照明是否提供了足够不同的信息、模型是否合适,以及先验是否成立。来源 1、2 的重建方式和已知条件不同,不能把结论合成无条件的唯一性保证。
改变波长则同时改变吸收谱与内部 F。若忽略 F 的波长依赖,输入光谱分解的曲线已经被照明“染色”。教学上可以先用已知混合比例的均匀仿体,再加入不同厚度的覆盖层,观察估计是否随覆盖层而漂移;这种对照比只展示一张彩色浓度图更能发现问题。
SECTION 055. 把定量目标写成可验证的实验
先明确输出是相对变化、吸收系数还是浓度,再为每种输出指定参考测量。记录入射脉冲能量有助于校正光源波动,但它测不到目标内部的 F。温度、接收增益、照明位置和探头位置也应固定或独立记录,否则同一图像值的变化可能有多个解释。
建议保留三组结果:原始或能量归一化后的声学量、采用的光通量模型、最终参数及其误差。对于无法独立核实的内部 F,应把它作为不确定度来源写出来。相关性、重建残差和参数准确度分别回答不同问题,三者都需要比“图像看起来合理”更明确的评价。
带走这三个要点
- 初始声压包含吸收与局部光通量的乘积,亮度差不能直接解释成浓度差。
- 多照明和多波长增加观测,但有效性取决于模型、先验与数据条件。
- 入射能量归一化、可见性增强和定量光通量校正应分别报告。
参考资料与继续阅读
以下资料用于支持或进一步核对本文涉及的概念。本文的教学算例与组织方式由本站编写;具体研究结论应以原文的实验条件和限制为准。
- Quantitative photoacoustic tomography with multiple optical sources ↗2010 · opg.optica.org
- Iterative algorithm for multiple illumination photoacoustic tomography (MIPAT) using ultrasound channel data ↗2012 · pmc.ncbi.nlm.nih.gov
编辑说明:初始文章由生成式 AI 根据所列公开资料辅助撰写,尚未经过领域专家审稿。示意与假设算例不构成研究数据,涉及健康的内容不替代专业医疗判断。阅读完整说明 ↗