2. 多模态医学影像数据获取与配准
多模态影像组学融合建模的第一步，是把来自不同设备、不同时间、不同物理机制的图像放到同一个空间参照系里。CT、MRI、PET 之所以能互补，正是因为它们对组织的“敏感点”完全不同：CT 反映电子密度与 X 射线衰减，骨皮质呈高密度、空气呈低密度，天然适合做解剖定位与衰减校正；MRI 依赖氢质子密度及 T1、T2、质子密度等弛豫特性，对软组织边界、脑灰白质、肿瘤浸润范围更敏感；PET 则探测正电子湮灭产生的 511 keV 光子对，反映葡萄糖代谢、受体表达等分子功能信息，但空间分辨率通常低于 CT 和 MRI。理解这些差异，才能理解为什么配准不是简单的“图像叠加”，而是要在解剖结构与功能信号之间建立可解释的空间对应。
成像原理差异带来的配准约束
CT 的像素值以 HU（Hounsfield Unit）表示，水约为 0 HU，骨可达 +1000 HU 以上，空气约 -1000 HU，具有较好的几何保真度，通常作为配准的参考空间。MRI 没有统一的绝对灰度标准，同一组织在不同序列、不同场强、不同线圈下信号强度可能显著变化，且存在几何畸变，尤其是 EPI 序列在脑边缘和鼻窦附近更明显。PET 图像本质是低分辨率、高噪声的功能图，病灶边界模糊，单独依赖灰度相似性容易失败。因此，跨模态配准通常采用“解剖引导功能”的策略：以 CT 或 T1 加权 MRI 作为固定参考，将 PET 或其它功能图变换到参考空间。
配准的数学本质是寻找一个空间变换 T，使浮动图像 I_F 经过 T 后与固定图像 I_M 在某种相似性度量下达到最优。变换可分为刚体、仿射、非线性形变三类。刚体变换包含 3 个平移和 3 个旋转参数，适合同一受试者头部或骨骼等刚性结构；仿射变换增加 3 个缩放和 3 个剪切参数，可校正不同扫描仪的尺度差异；非线性变换则用 B 样条或薄板样条等模型描述局部形变，适合腹部、肺部等受呼吸和器官运动影响的区域。相似性度量方面，同模态常用均方误差或归一化互相关，跨模态则更常用互信息（MI）或归一化互信息（NMI），因为它们不要求灰度值线性相关，只要求统计依赖性。
跨模态配准的实操路径
以脑部 CT-MRI-PET 三模态配准为例，一个可复现的流程如下。首先进行预处理：将 DICOM 序列转换为 NIfTI 格式，检查层厚、层间距、像素间距和方向矩阵，确保各模态的坐标系定义一致。若 PET 与 CT 来自同一台 PET/CT，通常已由设备完成刚性配准，可直接使用；若 PET 与 MRI 来自不同设备，则需要独立配准。接着进行粗配准：对头部数据，可用解剖标志点（如鼻根、外耳道、眼眶内缘）手动选取 4～6 对对应点，计算初始刚体变换；对腹部数据，可先用主轴对齐或基于分割的质心对齐。粗配准的作用是避免优化陷入局部极值。
精配准阶段，以 CT 为固定图像，MRI 为浮动图像，采用归一化互信息作为相似性度量，优化器可选 Powell 或 L-BFGS，多分辨率策略通常设置 3 层：先 4 倍下采样做全局搜索，再 2 倍下采样细化，最后全分辨率微调。变换模型先用刚体，若发现脑室、皮层等区域仍存在明显错位，再切换为 B 样条非线性变换，网格间距可设为 20～30 mm，正则化项权重约 0.01～0.1，防止形变过度。PET 配准到 CT 时，由于 PET 分辨率低，建议先对 PET 做高斯平滑（FWHM 约 4～6 mm），再以 CT 为参考做刚体或仿射配准，避免非线性变换把噪声拟合为形变。
重采样与 ROI 对齐
配准完成后，需要把浮动图像重采样到参考空间的统一网格。重采样包括三个关键参数：目标体素间距、插值方式和输出矩阵。对于 CT 参考空间，若原始层厚为 1.25 mm、像素间距 0.7 mm×0.7 mm，可重采样为 1 mm×1 mm×1 mm 各向同性体素，便于后续组学特征计算。插值方式上，CT 和 MRI 灰度图常用三线性插值或 B 样条插值，以保持灰度连续性；PET 可用三线性插值；而 ROI 掩膜、分割标签必须使用最近邻插值，否则会引入不存在的标签值。ROI 对齐时，需将配准变换矩阵同时应用于掩膜，再在参考空间重新生成二值掩膜，并检查掩膜体积变化是否在合理范围内，例如脑肿瘤 ROI 体积变化超过 15% 时应回查配准。
配准质量评估案例
假设一例脑胶质瘤患者的 T1 加权 MRI 与 CT 进行刚体配准，评估结果如下：视觉检查显示脑室边缘、颅骨内板、眼球轮廓在融合图上无明显重影；计算得到归一化互信息从初始的 0.18 提升到 0.42；以手动勾画的肿瘤 ROI 为参考，配准后 MRI ROI 与 CT ROI 的 Dice 系数为 0.86，Jaccard 为 0.75；选取 5 个解剖标志点计算 TRE，平均误差为 1.8 mm，最大误差 2.6 mm。综合判断该配准可用于后续影像组学特征提取。若 Dice 低于 0.7 或 TRE 大于 3 mm，则应检查标志点选取、预处理方向矩阵、相似性度量选择，必要时重新初始化并增加多分辨率层数。