2. 医学影像数据基础与多中心异质性

多中心医学影像模型验证的核心困难,往往不在模型结构本身,而在数据。同一套模型在 A 院表现良好、在 B 院明显下降,最常见的原因不是算法失效,而是两站数据在模态、协议、设备、人群四个维度上并不等价。理解这些差异如何进入像素、又如何转化为“站点效应”,是设计验证方案的前提。

一、影像模态:先明确“同一任务”是否可比

模态决定了图像所承载的物理信息。CT 反映 X 射线衰减系数(HU),MRI 反映质子密度与弛豫时间,超声反映声阻抗界面,PET 反映放射性示踪剂分布。多中心验证中,模态层面的异质性有三类:

  • 模态不同:如 A 站用增强 CT、B 站用 MRI 评估同一病灶,此时模型输入分布根本不同,通常不能直接合并验证,只能做分层评估。
  • 模态相同但序列/参数不同:如 MRI 的 T1WI、T2WI、DWI、FLAIR 属于不同对比度,混用会引入系统性偏差。
  • 模态相同但重建方式不同:如 CT 的滤波反投影(FBP)与迭代重建(IR),同一患者同一扫描会得到不同噪声纹理。

因此第一步不是“把数据放一起”,而是建立模态清单:每个站点、每个任务、每种序列的样本量、层厚、像素间距、窗宽窗位。只有清单对齐,后续比较才有意义。

二、采集协议:参数差异如何进入像素

采集协议是异质性最密集的来源。以 CT 为例,管电压(kVp)改变会改变 HU 与噪声;管电流(mAs)影响量子噪声;层厚(如 1 mm 与 5 mm)影响部分容积效应;重建核(如软组织核与骨核)影响边缘锐度。MRI 中,TR/TE 决定对比度,翻转角影响信号强度,并行采集因子影响信噪比与伪影。这些参数并非“设置项”,而是直接改变像素统计分布。

一个可操作的判断是:若两站同一模态的像素强度直方图与噪声功率谱差异显著,则模型在特征层面对应的响应也会偏移。验证前应至少记录:设备型号、扫描协议名称、kVp/mAs、层厚、重建核、TR/TE、场强、线圈类型。缺失这些字段,后续无法解释性能下降来自哪里。

三、设备差异:同一协议在不同机器上并不相同

即使协议名称一致,不同厂商、不同代际设备的探测器效率、散射校正、迭代重建强度、金属伪影抑制算法都不同。常见现象是:A 站 GE 设备与 B 站 Siemens 设备采用“相同”腹部协议,但 B 站图像更平滑、噪声更低,导致依赖纹理的模型在 B 站假阴性升高。

设备差异还体现在几何与标定上:像素间距、床高、重建视野(FOV)、层间距是否连续。若模型使用绝对尺寸特征(如病灶直径毫米数),FOV 与像素间距不一致会直接造成尺度误差。因此跨站验证必须统一到物理坐标,而不是像素坐标。

四、人群分布:病例构成改变先验概率

站点效应不仅来自机器,也来自人。不同医院的转诊结构、疾病谱、年龄分布、合并症比例不同。例如 A 站为肿瘤专科医院,B 站为综合医院,同一“肺结节”任务中,B 站良性结节比例更高,模型阈值不变时假阳性会上升。人群分布差异会改变患病率,进而改变阳性预测值,即使敏感度、特异度不变,临床可用性也会变化。

此外,扫描体位、呼吸配合、对比剂使用习惯、图像后处理流程也属于人群—流程耦合因素。验证报告中应给出各站点的年龄、性别、疾病亚型、病灶大小分布,用于判断性能差异是否可由病例构成解释。

五、站点效应:从“差异”到“可测量的偏移”

站点效应(site effect)指站点身份本身对模型输出或性能产生的系统性影响。它可由模态、协议、设备、人群中的一项或多项叠加形成。工程上常用两种方式观察:

  • 性能分层:按站点分别计算 AUC、敏感度、特异度,观察置信区间是否重叠。
  • 分布偏移检测:对模型中间特征或输入像素做站点间统计检验(如 KS 检验、MMD),定位偏移发生在输入层还是特征层。

若 A 站 AUC 为 0.91、B 站为 0.78,且 B 站图像噪声更高、层厚更厚,则优先怀疑协议与设备;若两站图像统计接近但 B 站阳性率显著更低,则优先怀疑人群分布。区分这两类原因,决定了后续是“做图像标准化”还是“做阈值/校准调整”。

多中心异质性来源与站点效应传导路径 影像模态 CT / MRI / PET / US 序列与对比度不同 采集协议 kVp / mAs / 层厚 TR / TE / 重建核 设备差异 厂商 / 代际 / 探测器 FOV / 像素间距 人群分布 疾病谱 / 患病率 年龄 / 合并症 像素统计分布改变 强度直方图 · 噪声功率谱 · 纹理 模型输入/特征偏移 特征分布漂移 · 决策边界错位 站点效应:AUC / 敏感度 / 特异度分层差异 例:A 站 AUC 0.91,B 站 AUC 0.78(示例) 偏移归因判断 图像统计差异大 → 协议 / 设备标准化 图像接近但阳性率低 → 阈值 / 校准调整 特征层漂移显著 → 域适应 / 重训练

六、验证前的数据对齐清单

把上述维度落到可执行层面,建议在建模前完成一份站点数据对齐表。以下为示例字段(数值为合理假设,用于说明格式):

站点模态/序列设备层厚(mm)kVp重建核样本量阳性率
A 院CT 平扫GE 64 排1.25120软组织4200.38
B 院CT 平扫Siemens 128 排3.0100B30f3600.22
C 院CT 增强Philips 64 排1.0120标准2800.45

该表一旦建立,就能快速判断:B 站层厚更厚、kVp 更低、阳性率更低,其性能下降可能同时来自协议与人群,需分别做图像重采样与阈值校准实验,而不是简单归因于“模型泛化差”。

实践提示:跨站验证时,先做“同站内交叉验证”作为上限参考,再做“跨站外部验证”。若同站内性能也低,问题在模型或标注;若同站内高、跨站低,问题在异质性,应回到本节的四个维度逐项排查。

总结来看,模态决定可比性,协议与设备决定像素分布,人群决定先验概率,三者叠加形成站点效应。验证工程的第一步不是训练,而是把数据差异说清楚、量出来、对齐好。只有输入层面对齐,后续的性能比较才具有解释力。