1. 多模态医学影像配准工程导论
多模态医学影像配准，是把来自不同成像设备、不同时间点或不同患者姿态的影像，变换到同一个空间坐标系下，使解剖结构在空间上对齐。它既不是单纯的图像处理算法题，也不是纯粹的临床操作，而是一项贯穿"成像物理—图像表示—空间变换—相似性度量—优化求解—临床验证"的系统工程。本课程的目标，是让学习者能够独立完成一次可复现、可评估、可交付的多模态配准任务，并理解每一步背后的物理与数学依据。
为什么临床与科研离不开配准
不同模态的影像回答的是不同问题。CT 以电子密度为基础，对骨性结构、钙化和植入物显示清晰，且具有天然的几何精度，常用于放疗剂量计算与手术导航；MRI 软组织对比度优异，T1、T2、DWI、FLAIR 等序列分别突出解剖、水肿与白质病变，是脑、前列腺、肝脏等软组织诊断的主力；PET 反映代谢与受体分布，能发现 CT/MRI 上形态尚不明显的功能异常，但其空间分辨率低、解剖边界模糊；超声实时、无辐射、可床旁操作，却视野受限、切面依赖操作者、且存在斑点噪声与几何畸变。
正因为单一模态无法同时提供"精确几何 + 优良软组织对比 + 功能代谢信息"，临床才需要把它们叠加起来。典型场景包括：PET/CT 融合让核医学医生把高代谢灶精确定位到某个肺段或淋巴结；MRI 与 CT 融合用于脑肿瘤放疗靶区勾画，把 MRI 上的水肿区映射到 CT 计划系统；术前 MRI 与术中超声配准用于肝肿瘤消融导航；随访时把新一次 MRI 与基线 MRI 对齐，量化肿瘤体积变化。科研中，配准是群体分析（VBM、纤维束追踪、影像组学）的前置步骤——没有空间对齐，跨被试的统计就失去意义。
四种模态的成像特性与配准难点
理解配准，必须先理解各模态"长什么样、错在哪"。下表给出工程上最需要关注的差异。
| 模态 | 物理对比来源 | 空间分辨率 | 典型配准难点 |
|---|---|---|---|
| CT | X 射线衰减（电子密度） | 0.5–1 mm 各向同性 | 金属伪影、呼吸相位差异、床板形变 |
| MRI | 质子密度、T1/T2 弛豫 | 1–3 mm，层厚常大于层内 | 几何畸变（EPI）、强度不均匀、序列间灰度非线性 |
| PET | 放射性示踪剂分布 | 4–6 mm，部分容积效应强 | 低分辨率、灰度与解剖无对应、SUV 归一化差异 |
| 超声 | 声阻抗界面反射 | 0.3–1 mm，但依赖切面 | 斑点噪声、阴影、探头压迫形变、无固定坐标系 |
这些差异直接决定了配准策略：CT–CT 或 MRI–MRI 属单模态，灰度关系近似一致，可用灰度差平方和（SSD）或归一化互相关（NCC）；CT–MRI、CT–PET 属多模态，灰度无线性对应，必须使用互信息（MI）或其归一化形式（NMI）这类基于统计依赖的度量；超声与 CT/MRI 配准则常需先做特征提取或采用基于物理模型的形变约束。
配准任务的分类维度
工程上描述一个配准任务，通常沿四个维度定位：
- 变换类型：刚体（6 自由度：3 平移 + 3 旋转）、仿射（12 自由度，含缩放与剪切）、形变（B 样条自由形变、Demons、光流等，自由度可达数千）。头骨内脑组织可用刚体近似，腹部因呼吸与脏器滑动必须用形变。
- 空间维度：2D–2D（如 DSA 与血管造影）、2D–3D（术中透视与术前 CT）、3D–3D（PET/CT、MRI/CT）。
- 配准依据：基于灰度（MI、NCC）、基于特征（点、线、面、SIFT/SURF）、基于物理模型（有限元、生物力学）。
- 交互方式：刚性手动、半自动（选点初始化 + 自动优化）、全自动。
下面这张图把"模态特性—配准依据—变换类型—典型应用"串成一条可对照的决策链，帮助你在拿到一对影像时快速判断该走哪条路。
从"能配"到"配得准"：工程闭环
一次完整的配准工程，通常包含六个环节：数据准备（DICOM 读取、层厚与像素间距校正、重采样到各向同性）、预处理（去噪、偏置场校正、颅骨剥离或掩膜生成）、初始化（质心对齐、主轴对齐或手动选点）、相似性度量与优化（选择 MI/NCC，设定 Powell、L-BFGS 或梯度下降，控制多分辨率层数与迭代次数）、变换应用与重采样（最近邻用于标签，线性/B 样条用于灰度）、评估与验证。
以 CT–MRI 脑部配准为例，一个可复现的示例参数是：将 MRI 重采样到 1 mm 各向同性，CT 保持原始 0.8 mm；使用 NMI 作为度量，多分辨率 3 层（4×、2×、1× 下采样），优化器为 L-BFGS，最大迭代 200；刚体 6 自由度；最终以互信息收敛曲线和棋盘格叠加图判断。若 TRE 超过 2 mm，则检查是否因颅骨剥离不一致或初始化偏差过大导致陷入局部极值。这些数值是教学示例，实际项目需根据数据重新标定。