4、光声信号特性：信号时域与频域特征、幅度与相位信息、噪声来源与信噪比概念
光声成像的核心在于将光能量转化为超声信号，再通过超声换能器进行采集与重建。理解光声信号的物理本质与数学表达，是设计成像系统、优化采集参数、提升图像质量的前提。本节围绕光声信号的时域与频域特征、幅度与相位信息、噪声来源与信噪比（SNR）概念展开，帮助学习者建立从信号产生到信号处理的全链路认知框架。
4.1 光声信号的时域特征
光声信号在时域上表现为典型的双极性脉冲，即“N形波”。当短脉冲激光照射吸收体时，吸收体瞬间升温并发生热弹性膨胀，产生初始声压脉冲。该脉冲在传播过程中，由于介质黏滞吸收和几何扩散，逐渐演化为一个正压峰后跟随一个负压尾波的波形。正压峰的幅度与吸收体的光吸收系数、局部光通量以及Grüneisen系数成正比，而负压尾波的幅度和持续时间则与吸收体的几何尺寸和声速相关。
时域特征的关键参数包括：
- 峰值幅度（pmax）：反映吸收体的光吸收强度，是定量成像的基础。
- 上升时间（tr）：从基线到正压峰的时间，与吸收体尺寸和声速有关，通常为纳秒量级。
- 脉冲宽度（τ）：正压峰的持续时间，决定了轴向分辨率，τ ≈ 0.5·d/c（d为吸收体直径，c为声速）。
- 负压尾波幅度：与吸收体的声阻抗失配程度相关，可用于区分不同组织类型。
在实际系统中，换能器的带宽和频率响应会改变接收到的波形形状。宽带宽换能器（如中心频率5 MHz、带宽80%）能较好地保留N形波特征，而窄带宽换能器则会使波形展宽、峰值降低。因此，时域波形的解读必须结合换能器的脉冲响应进行去卷积处理。
4.2 光声信号的频域特征
光声信号的频域特征由其时域波形决定。理想点状吸收体产生的光声信号为无限带宽的δ函数，但实际吸收体具有有限尺寸，其频谱呈现带通特性。对于直径为d的球形吸收体，其频谱的主瓣宽度约为 fc ≈ 0.66·c/d，且高频分量随频率增加而衰减。这意味着小尺寸吸收体产生更高频的信号，而大尺寸吸收体则主要辐射低频成分。
频域分析在光声成像中具有重要应用：
- 频域滤波：通过带通滤波器抑制低频漂移和高频噪声，保留有效信号频带。
- 频域重建：在傅里叶域进行重建算法（如F-K迁移），可提高重建速度和图像质量。
- 频谱定量分析：通过计算频谱斜率或中心频率偏移，可以估计吸收体的尺寸和光学参数。
| 吸收体直径 | 中心频率（估算） | 典型应用场景 |
|---|---|---|
| 10 μm（红细胞） | ~100 MHz | 微血管成像、光声显微 |
| 100 μm（小血管） | ~10 MHz | 小动物脑成像 |
| 1 mm（肿瘤） | ~1 MHz | 乳腺光声断层 |
4.3 幅度与相位信息
光声信号的幅度携带吸收体的光吸收定量信息，而相位则携带声波传播路径和介质声速的信息。在光声成像中，幅度信息主要用于重建吸收分布图，相位信息则用于校正声速不均匀性引起的畸变。
幅度信息的提取需要注意以下几点：
- 系统增益校准：不同通道的放大器增益可能存在差异，需进行一致性校准。
- 光通量归一化：激光能量在组织中的衰减导致不同深度处的光通量不同，需进行光通量补偿。
- 换能器响应校正：换能器的频率响应不均匀，需通过去卷积恢复真实信号幅度。
相位信息在光声成像中的应用相对较少，但在声速层析成像中具有重要意义。通过测量不同角度接收信号的相位差，可以反演介质的声速分布，进而校正重建图像中的几何畸变。此外，在光声相干成像中，相位信息用于实现合成孔径聚焦，提高横向分辨率。
4.4 噪声来源与信噪比概念
光声成像系统的噪声来源主要包括以下几类：
| 噪声类型 | 来源 | 特征 | 抑制方法 |
|---|---|---|---|
| 热噪声 | 换能器阻抗、前置放大器 | 白噪声，与频率无关 | 低噪声放大器、阻抗匹配 |
| 散粒噪声 | 光探测器的统计涨落 | 与信号幅度平方根成正比 | 提高激光能量、多次平均 |
| 电噪声 | 电源纹波、电磁干扰 | 50 Hz工频及谐波 | 屏蔽、滤波、差分采集 |
| 声学噪声 | 环境声、耦合介质气泡 | 低频宽带噪声 | 隔声处理、脱气水耦合 |
信噪比（SNR）定义为信号功率与噪声功率之比，通常以分贝（dB）表示：
SNR = 20 · log₁₀(A_signal / A_noise)其中A_signal为信号幅度，A_noise为噪声均方根幅度。在光声成像中，提高SNR的主要途径包括：
- 增加激光能量：但需注意安全阈值（ANSI标准限制最大允许曝光量）。
- 多次平均：N次平均可使SNR提升√N倍，但会降低成像帧率。
- 优化换能器设计：提高灵敏度、匹配声阻抗、降低热噪声。
- 信号滤波：匹配滤波或小波去噪可在不损失分辨率的前提下提升SNR。
4.5 案例拆解：小动物脑成像中的SNR优化
背景与约束：某实验室搭建了一套光声断层成像系统，用于小鼠脑部血管成像。系统采用10 MHz中心频率、70%带宽的环形换能器阵列（256阵元），激光波长750 nm，脉冲能量8 mJ/cm²（低于安全阈值）。初始图像中，深部血管（深度>5 mm）的信噪比仅为8 dB，无法清晰分辨小血管分支。
初始现象与目标：浅表血管（深度<3 mm）SNR约25 dB，成像清晰；但深部区域信号衰减严重，SNR降至8 dB以下。目标是将深部血管SNR提升至15 dB以上，同时保持成像帧率不低于2 fps。
分析判断过程：
- 首先检查系统噪声基底：关闭激光，采集纯噪声数据，测得噪声RMS为0.35 mV。
- 测量深部信号幅度：在深度5 mm处，信号峰值幅度约为1.2 mV，计算SNR = 20·log₁₀(1.2/0.35) ≈ 10.7 dB，与实测8 dB接近（考虑带宽内噪声积分）。
- 分析衰减原因：750 nm波长在脑组织中的有效衰减系数约0.8 cm⁻¹，5 mm深度处光通量衰减至表面的约67%。
- 评估换能器灵敏度：10 MHz换能器在5 mm深度处的接收灵敏度下降约30%。
实际操作步骤：
- 将激光能量从8 mJ/cm²提升至12 mJ/cm²（仍低于ANSI安全阈值），信号幅度提升至1.8 mV。
- 采用64次平均（原为16次平均），噪声RMS降至0.12 mV，SNR提升至20·log₁₀(1.8/0.12) ≈ 23.5 dB。
- 在信号处理链路中加入带通滤波器（5-15 MHz），滤除低频漂移和高频噪声，进一步将噪声RMS降至0.09 mV。
- 对深部区域进行光通量补偿算法校正，恢复幅度一致性。
结果对比：优化后深部血管SNR达到18.5 dB，小血管分支（直径约100 μm）清晰可辨。成像帧率从4 fps降至1.5 fps，但通过并行采集和GPU加速重建，最终维持在2 fps，满足实验需求。
风险与代价：激光能量提升50%增加了光热损伤风险，需实时监测皮肤温度；64次平均使采集时间延长4倍，小鼠呼吸运动可能导致轻微运动伪影，需采用门控采集或图像配准算法。
复盘经验：SNR优化需要系统性地平衡激光能量、平均次数、滤波策略和重建算法。单纯增加激光能量或平均次数并非最优解，应结合系统噪声基底分析和信号特征进行综合优化。此外，光通量补偿对于深部成像至关重要，否则即使SNR提升，重建图像的定量准确性仍受影响。