4. 影像组学数据生命周期与标准化节点

影像组学研究的可重复性危机,很少源于算法本身,而更多来自数据在生命周期各环节中悄然发生的语义漂移。同一批 DICOM 图像,若在采集时层厚不同、重建核不同,在存储时被匿名化工具改写了私有标签,在标注时勾画边界相差两个像素,在特征提取时重采样插值方式不一致,最终得到的特征向量可能面目全非。因此,标准化不是某个单点动作,而是贯穿采集、存储、标注、特征提取、建模共享五个阶段的控制链。本节以一条真实可追踪的数据流为主线,说明每个阶段应固定什么、记录什么、校验什么。

一、采集阶段:把扫描协议当作实验条件来管理

影像组学特征对采集参数高度敏感。以 CT 为例,管电压(kVp)、管电流(mAs)、层厚、重建核(kernel)、像素间距、是否使用迭代重建,都会改变纹理特征的数值分布。标准化控制点包括:同一研究内尽量统一扫描协议;若无法统一,必须逐例记录关键参数,并在建模时将其作为协变量或分层因素。对于 MRI,场强、序列类型(T1WI/T2WI/DWI)、TR/TE、翻转角、是否压脂同样需要登记。

采集阶段还应完成设备与受试者标识的规范化。建议在扫描时即写入研究编号而非患者姓名,避免后续匿名化时因私有标签残留导致信息泄露。一个可操作的检查项是:导出 DICOM 头后,确认 PatientName、PatientID、StudyDate 等字段已按研究方案替换,且 BurnedInAnnotation 为 NO,防止像素中烧录了姓名。

二、存储阶段:DICOM 到研究数据集的受控转换

原始 DICOM 应保持只读归档,任何转换都生成新副本,确保可回溯。存储阶段的核心标准化动作是建立统一的目录结构与命名规则,例如按 研究编号/模态/序列/日期 分层,文件名使用 研究编号_模态_序列号.dcm。同时,将 DICOM 转换为 NIfTI 或 nrrd 时,必须固定方向约定。这里需要先明确参照系:本文及后续图示统一采用 RAS(+x 向右、+y 向前、+z 向上)作为世界坐标方向,若工具默认使用 LPS,应在转换脚本中显式翻转,并在元数据中记录。

存储阶段还要处理层厚不一致与缺失层问题。对于层厚大于 1 mm 的 CT,若后续要做三维纹理,需在提取前统一重采样到各向同性体素(如 1×1×1 mm³),并记录插值方法(线性或 B 样条)。这些参数一旦确定,应写入数据字典,而不是留给每位分析者自行选择。

影像组学数据生命周期标准化控制链 五个阶段 · 每阶段固定输入、处理参数与输出校验 ① 采集 kVp / mAs / 层厚 重建核 / 场强 / 序列 像素间距 / 压脂 控制点:协议统一 或逐例登记参数 ② 存储 DICOM 只读归档 匿名化 / 目录命名 DICOM→NIfTI 转换 控制点:RAS 方向 重采样 1mm³ 记录 ③ 标注 ROI 勾画 / 边界定义 双人独立 + 仲裁 Dice / 体积一致性 控制点:标注 SOP Dice ≥ 0.85 判合格 ④ 特征提取 重采样 / 灰度离散化 固定 bin 宽或 bin 数 IBSI 参数集 控制点:参数版本化 ICC 稳定性筛选 ⑤ 建模共享 特征矩阵 训练/验证划分 模型卡 / 数据字典 控制点:可复现 版本与许可 特征不稳定 → 回溯标注/采集参数复核 贯穿全流程的数据字典(示例字段) 字段名 示例值 标准化要求 slice_thickness 1.0 mm 统一或作为协变量 resample_spacing 1×1×1 mm³ 线性插值,记录方法 bin_width 25 HU 固定 bin 宽,全队列一致 roi_dice 0.88 ≥0.85 方可入组

三、标注阶段:把“勾画”变成可度量的过程

标注是影像组学中人为变异最大的环节。标准化不是要求所有人画出完全相同的边界,而是要求边界定义可复现、可评估。具体做法是:先制定标注 SOP,明确 ROI 是否包含瘤周、是否排除坏死区、是否沿肿瘤外缘还是内缘勾画;再由两名标注者独立勾画同一批 20~30 例,计算 Dice 相似系数与体积差异。若 Dice 低于 0.85,应回到 SOP 讨论边界定义,而不是直接取平均。

对于多中心研究,还需评估不同中心标注者之间的系统性偏移。例如,中心 A 倾向于把 ROI 画大 5%,中心 B 偏小 3%,这种偏移会直接进入特征值。此时可在建模前对特征做中心内标准化(z-score),或把中心作为批次效应变量纳入 ComBat 校正。标注文件本身也应标准化:统一使用 NIfTI 掩膜,标签值 1 表示 ROI,0 表示背景,避免同一研究内出现 1/2/3 多标签混用。

四、特征提取阶段:参数集必须版本化

特征提取的标准化核心是“固定参数集并记录版本”。以 PyRadiomics 为例,一个可复现的配置应至少固定:重采样体素尺寸、插值方式、灰度离散化方法(固定 bin 宽或固定 bin 数)、是否启用小波或 LoG 滤波、以及特征类别(一阶、形状、GLCM、GLRLM 等)。若使用固定 bin 宽,CT 常用 25 HU;若使用固定 bin 数,则需说明 bin 数(如 32 或 64)。两者不可混用,因为同一 ROI 在两种策略下得到的 GLCM 纹理值不可直接比较。

提取完成后,应做稳定性筛选。常用方法是计算 ICC(组内相关系数):对同一批病例重复提取两次(如不同重采样顺序或不同标注者掩膜),ICC ≥ 0.75 的特征保留,低于阈值的特征剔除或标记为不稳定。这一步能显著减少后续建模中的过拟合风险。此外,所有特征值应输出为带表头的 CSV,列名遵循 IBSI 命名规范,便于跨工具核对。

五、建模与共享阶段:让结果可被他人重放

建模阶段的标准化重点不在算法选择,而在数据划分与文档。训练集、验证集、测试集应按患者而非按图像划分,避免同一患者的多层图像同时进入训练和测试造成泄漏。若做交叉验证,应固定随机种子,并记录每一折的划分文件。特征矩阵、标签、划分索引、模型超参数、软件版本应一并归档。

共享阶段建议提供三样东西:数据字典(字段名、单位、缺失值编码)、模型卡(训练数据来源、特征数量、算法、性能指标及置信区间)、以及最小可复现脚本(从原始 DICOM 到最终预测的完整命令)。若数据涉及隐私,可共享去标识化后的特征矩阵与掩膜,而不共享原始图像。这样,其他研究者至少能验证特征与模型的一致性,而不必重新采集数据。

关键判断依据:当特征在不同中心间分布差异显著(如批次效应导致 AUC 下降超过 0.1),应优先回溯采集与重采样参数,而非直接更换分类器。标准化节点上的一个小偏差,往往比模型选择带来更大的性能波动。

综上,影像组学数据生命周期的标准化,本质上是把每个阶段的“默认选择”变成“显式记录”。采集时记录协议,存储时固定方向与重采样,标注时量化一致性,提取时版本化参数,建模共享时保留划分与文档。这五个控制点连成一条可追溯的链,才能让影像组学特征从“一次性的实验结果”变成“可积累的研究资产”。