3. 超声图像退化建模与预处理
超声成像的物理链路决定了图像质量不可能“天然干净”。从探头发射脉冲、声波在组织中传播、界面反射、回波被接收,到波束合成、包络检波、对数压缩、扫描转换,每一步都会引入特定类型的退化。工程上把这些退化归纳为四类:斑点噪声、深度衰减、时间增益补偿失配、扫描转换伪影。理解它们的成因与数学形态,是设计抑制策略的前提。
一、四类退化的成因与图像表现
斑点噪声(speckle)并非电子噪声,而是相干成像的固有结果。当分辨率单元内存在大量随机分布的散射子时,各回波复振幅随机叠加,形成乘性干涉图样。其幅度近似服从瑞利分布,方差与均值平方成正比,因此表现为“亮暗颗粒”而非独立像素噪声。关键工程含义是:斑点与组织信号相乘,简单低通滤波会同时模糊边界。
深度衰减源于声波吸收与散射,回波幅度随往返深度近似指数下降。设组织平均衰减系数为 α(dB/cm/MHz),发射频率 f,深度 z,则单程衰减约 α·f·z dB。5 MHz 探头在 6 cm 深度、α=0.5 dB/cm/MHz 时,单程衰减约 15 dB,往返约 30 dB。若不补偿,深部组织将淹没在暗区中。
时间增益补偿(TGC)正是为抵消衰减而设的深度相关增益。它按深度线性或分段提升接收增益,使同质组织在不同深度呈现相近亮度。TGC 失配的典型表现是:增益曲线过陡导致深部过亮、近场过暗;过平则深部仍暗。临床上操作者调节的“近场/远场增益”旋钮,本质就是修改 TGC 曲线。
扫描转换伪影发生在极坐标→直角坐标重采样阶段。扇形扫描的极坐标数据被映射到矩形像素网格时,若插值核选择不当或网格过粗,会出现近场横向模糊、远场阶梯状锯齿、以及扇形边缘的“空洞”或重复条纹。此外,声影、混响、镜像等传播伪影也会在扫描转换后被放大。
二、斑点噪声的工程化抑制
斑点抑制的核心矛盾是“去噪”与“保边”。工程上常用三类方法。第一类是自适应空域滤波:Lee 滤波假设局部均值与方差可估计,按 I_out = I_mean + k·(I_in − I_mean) 加权,其中 k = 1 − (σ_n²/σ_local²),在均匀区强平滑、在边缘区弱平滑。Frost 滤波则用指数衰减核,核宽随局部变异系数自适应。第二类是各向异性扩散(SRAD),以偏微分方程迭代,扩散系数在梯度大处趋近于零,从而在平滑斑点的同时保留边界。第三类是多帧平均:对同一解剖位置采集 N 帧,斑点因去相关而平均,信噪比提升约 √N,但要求帧间运动小于半波长,否则会引入运动模糊。
实际工程中,建议先做对数压缩前的同态滤波:将乘性斑点转为加性,再用 Lee 或 NLM 处理,最后反变换。参数示例:Lee 滤波窗口 5×5,噪声方差由均匀区估计;SRAD 迭代 20 次,时间步长 0.1,扩散阈值取局部变异系数 0.8。这些参数需根据探头频率和显示动态范围微调。
三、衰减补偿与 TGC 标定
TGC 标定的目标是让同质体模在不同深度呈现一致亮度。操作流程如下:将组织等效体模(如 0.5 dB/cm/MHz)置于探头下,固定发射功率和总增益,从近场到远场逐段读取平均灰度;若远场灰度低于近场 10% 以上,则提升对应深度段的 TGC 增益,直至曲线平坦。记录各段增益值,形成该探头的 TGC 查找表。实际扫描时,系统按深度索引查表并实时施加增益。
需要注意,TGC 是深度的一维函数,无法补偿横向衰减差异。对于脂肪层厚、声影明显的病例,可结合深度增益补偿(DGC)和侧向增益补偿(LGC),后者沿阵元方向调整增益,抑制因声束偏转导致的边缘变暗。
四、扫描转换伪影的识别与抑制
扫描转换伪影的典型表现有三类。第一类是近场横向模糊:扇形扫描近场弧长小,极坐标采样密度低,重采样到直角网格时横向分辨率下降。抑制方法是提高近场采样密度或采用自适应插值核。第二类是远场阶梯锯齿:远场弧长大,若直角网格过粗,会出现明显的阶梯状边缘。工程上采用双三次插值或先上采样再插值,可将锯齿降低到 1 像素以内。第三类是扇形边缘空洞:极坐标数据在扇形边界外无定义,直接映射会产生黑色三角区。解决方法是边缘填充或使用扇形掩膜,避免插值核越界。
识别伪影时,可观察均匀体模图像:若近场出现横向条纹、远场出现周期性锯齿、边缘出现放射状暗带,即可判定为扫描转换伪影。此时应检查插值核类型、网格间距和扇形角度范围,必要时重新设计重采样网格。
五、预处理链的工程实现要点
一个可复用的预处理链通常包含:原始 RF 数据 → 正交解调 → 包络检波 → 对数压缩 → TGC 补偿 → 斑点滤波 → 扫描转换 → 显示。其中 TGC 补偿可在对数域或线性域进行,对数域实现更简单,但会放大低幅噪声;线性域更精确,但需注意数值范围。斑点滤波建议在对数域进行,因为乘性噪声在对数域变为加性,便于使用加性噪声模型。
参数配置示例:中心频率 5 MHz,采样率 40 MHz,TGC 分段 8 段,每段增益步进 2 dB;Lee 滤波窗口 5×5,噪声方差 0.02;扫描转换网格 0.2 mm×0.2 mm,双三次插值。这些参数需根据实际探头和组织类型调整,并通过体模实验验证。
最后,预处理不是“越干净越好”。过度滤波会抹除微小病灶的边界信息,过度增益补偿会放大深部噪声。工程上应建立客观评价指标,如斑点信噪比(SNR)、边缘保持指数(EPI)和对比度噪声比(CNR),在抑制退化和保留诊断信息之间取得平衡。