📚 课程目录第1章超声成像物理与图像形成

1. 超声成像物理与图像形成

超声图像增强与伪影处理的所有工程手段，最终都要落到一个物理事实上：屏幕上看到的每一帧灰度图，都是声波在人体组织里往返传播、被界面反射或散射、再被阵列换能器接收并经过波束合成与后处理重建出来的结果。如果不清楚这个链条中哪一步引入了什么信息、哪一步又可能引入什么伪影，那么后续的降噪、增强、去伪影操作就很容易变成"凭感觉调参"。本章先把这条链路拆开，明确每个环节的物理量、工程实现和它与伪影处理的对应关系。

1.1 声波传播：决定图像对比度与衰减的物理基础

医用超声通常工作在 2～15 MHz 频段，属于机械纵波。声波在介质中传播时，两个参数最关键：声速 c 和声阻抗 Z。软组织平均声速约 1540 m/s，脂肪约 1450 m/s，骨骼可达 3000～4000 m/s。成像系统默认全程按 1540 m/s 计算回波到达时间并换算深度，因此当声波穿过脂肪层或软骨时，实际声速与假设值不一致，就会造成几何位置偏移——这是很多"看起来结构错位"类伪影的根源。

声阻抗定义为 Z = ρc（密度×声速）。当声波垂直入射到两种介质界面时，反射系数为：

R = (Z₂ − Z₁) / (Z₂ + Z₁)
I_r / I_i = R²

例如软组织（Z≈1.63 MRayl）与空气（Z≈0.0004 MRayl）界面，R≈−0.999，几乎全反射，这就是为什么耦合不良时图像会出现大片暗区或强反射伪影。而软组织与血液之间 Z 差异很小，R 只有约 0.01，所以血管腔在 B 模式下呈现为低回声暗区，需要靠增强或血流模式来观察。

另一个必须记住的量是衰减。组织对声波的衰减近似随频率线性增长，典型软组织衰减系数约 0.5 dB/(cm·MHz)。这意味着 5 MHz 探头在 8 cm 深度处单程衰减约 20 dB，往返约 40 dB。工程上通过 TGC（时间增益补偿）按深度逐段放大回波来补偿，但 TGC 曲线是分段折线，补偿不完美时就会在图像上留下水平方向的明暗条带——这属于典型的增益类伪影，处理时不能简单当作噪声滤除。

1.2 回波接收与波束合成：空间信息的来源

阵列换能器由几十到几百个阵元组成。发射时通过控制各阵元激励延时，使声波在焦点处同相叠加，形成聚焦声束；接收时则对每个阵元收到的回波施加相反的延时补偿，再求和，这一过程称为波束合成（Beamforming）。它决定了图像的横向分辨率、旁瓣电平和焦点位置。

波束合成的核心是延时计算。对位于 (x_i, 0) 的第 i 个阵元，若目标点在 (x_p, z_p)，则往返延时为：

τ_i = ( √((x_i − x_p)² + z_p²) + z_p ) / c

实际系统中，c 固定为 1540 m/s，而真实组织声速有偏差，于是延时误差会导致相位失配，表现为横向模糊、焦点外能量泄漏。此外，阵元间距 d 若大于半波长 λ/2，会出现栅瓣（grating lobe），在图像上形成与主瓣方向对称的虚假回波——这是相控阵和线阵中常见的旁瓣/栅瓣伪影，处理时需要在波束合成层面加窗（如 Hamming、Hanning）来抑制。

波束合成之后，回波包络经过对数压缩、检波、DSC（数字扫描转换）映射到显示坐标系，才形成我们看到的 B 模式图像。对数压缩的动态范围通常设为 40～60 dB，压缩过强会丢失弱回声细节，压缩过弱则强回声区域饱和，两者都会影响后续增强算法的输入质量。

1.3 图像模式与工程增强的关联

B 模式提供的是灰度解剖结构，M 模式沿单条扫描线记录时间变化，彩色多普勒叠加血流速度信息，能量多普勒只显示血流存在性而不显示方向。不同模式的数据来源和噪声特性不同，增强策略也不能通用。例如 B 模式的主要噪声是斑点噪声（speckle），它来自大量微小散射体的相干叠加，表现为颗粒状纹理，既包含组织信息又干扰观察；而彩色多普勒的主要问题是速度混叠（aliasing）和壁滤波器引入的低速信号丢失。

工程上处理斑点噪声时，常用 Lee 滤波、Frost 滤波或非局部均值（NLM）。以 NLM 为例，其权重由图像块相似度决定：

w(p,q) = exp( −‖v(p) − v(q)‖² / h² )
NL[v](p) = Σ w(p,q)·v(q) / Σ w(p,q)

其中 h 是平滑参数，控制衰减强度。h 过大则细节被抹平，h 过小则噪声残留。实际调试时，可以先用 5×5 或 7×7 搜索窗，在斑点区域观察方差下降比例，同时用边缘保持指标（如 EPI）确认结构边缘未被破坏。这类参数关系正是后续章节要展开的内容，但它们的物理依据都来自本章所述的声波传播与回波统计特性。

工程提示：当图像出现"不该有的亮线"时，先判断它是否与声束方向对称。若对称，优先怀疑栅瓣或旁瓣；若沿深度方向水平分布，优先怀疑 TGC 或衰减补偿；若与已知解剖结构位置不符，优先怀疑声速失配。这个判断顺序能避免把物理伪影误当成随机噪声去滤除。

1.4 从物理量到处理参数：一个可复现的观察示例

假设使用线阵探头，中心频率 7.5 MHz，阵元间距 0.3 mm，声速按 1540 m/s 计算，则波长 λ≈0.205 mm，λ/2≈0.103 mm。阵元间距 0.3 mm 明显大于 λ/2，理论上会出现栅瓣。若在仿体或均匀组织区域观察到与主瓣对称的虚假回波，可尝试在波束合成中加 Hamming 窗，观察栅瓣幅度下降情况。这里给出的数值属于教学示例参数，用于说明判断逻辑，不代表某一台具体设备的实测结果。

同样，若在脂肪层较厚的腹部图像中看到深部结构位置偏浅或偏深，可以尝试将系统声速假设从 1540 m/s 调整为 1470 m/s 左右重新计算延时，观察目标点是否回到预期深度。这个调整过程本质上是在补偿声速失配，而不是改变图像后处理参数。

理解这条从声波到图像的完整链路，是后续所有增强与伪影处理工作的前提。只有知道每个伪影在物理上"从哪里来"，才能在选择算法时判断它"能不能被去掉"以及"去掉之后会不会损伤真实结构"。