📚 课程目录第4章配准数据采集与预处理工程

4. 配准数据采集与预处理工程

术前术中配准的精度上限，往往在数据进入配准算法之前就已经被决定了。影像读取方式、去噪强度、偏置场校正是否充分、重采样插值选择、感兴趣区域（ROI）裁剪范围，以及多源数据的时间戳对齐，都会直接影响后续相似性测度和优化过程的稳定性。本节围绕一条可复现的预处理流水线展开：从 DICOM/NIfTI 读取开始，到形成一份几何一致、灰度可比、时间可对齐的配准输入数据。

一、DICOM 与 NIfTI 的读取与几何信息解析

DICOM 是术中影像（C 臂、CBCT、超声）与术前 CT/MR 的主要载体，其几何信息分散在多个标签中：ImagePositionPatient (0020,0032) 给出切片左上角体素中心在患者坐标系中的位置，ImageOrientationPatient (0020,0037) 给出图像行、列方向余弦，PixelSpacing (0028,0030) 给出面内像素间距，SliceThickness (0018,0050) 与相邻切片位置差共同决定层间距。读取时必须按 InstanceNumber 或投影后的位置排序，否则会出现层序错乱。

NIfTI 则把几何信息集中在头信息中：sform/qform 矩阵描述体素索引到世界坐标的仿射变换，pixdim 给出各轴体素尺寸。工程上建议统一以 NIfTI 的仿射矩阵作为内部几何表示，把 DICOM 序列重建为带仿射的体数据，这样后续重采样、ROI 裁剪和配准变换都在同一坐标系下运算。一个常见的坑是 DICOM 的 RescaleSlope/Intercept 未应用，导致 CT 值不是 HU；正确做法是 HU = pixel * slope + intercept，再进入后续处理。

# 示例：DICOM 序列读取与 HU 转换（伪代码）
slices = read_dicom_series(folder)
slices.sort(key=lambda s: s.ImagePositionPatient[2])
volume = stack([s.pixel_array * s.RescaleSlope + s.RescaleIntercept
                for s in slices])
spacing = (slices[0].PixelSpacing[1], slices[0].PixelSpacing[0],
           slices[1].ImagePositionPatient[2] - slices[0].ImagePositionPatient[2])
affine = build_affine(slices[0], spacing)   # 由方向余弦与首层位置构造

二、去噪与偏置场校正

术中 C 臂和 CBCT 图像受量子噪声影响明显，MR 则存在射频场不均匀导致的低频偏置场。去噪的目标不是把图像“磨平”，而是在抑制噪声的同时保留血管、骨缘、电极等配准依赖的高频结构。工程上常用三种策略：高斯平滑（快但模糊边缘）、中值滤波（对椒盐噪声有效）、非局部均值或双边滤波（保边性好，计算量较大）。示例参数：对 0.5 mm 层厚的 CBCT，可采用 3×3×3 中值滤波去除孤立噪点，再用 σ=0.8 体素的高斯核轻度平滑；若用于血管配准，则改用双边滤波，空间 σ=1.5、灰度 σ=30 HU，以保留对比剂边缘。

偏置场校正主要针对 MR。常用 N4ITK 算法：先对图像取对数，用 B 样条拟合低频场，迭代估计并除以该场。示例配置为 4 级分辨率、每级 50 次迭代、B 样条网格 4×4×4，收敛容差 1e-4。校正后应检查白质区域灰度是否趋于均匀——若仍存在明显亮暗梯度，说明网格过粗或迭代不足。需要注意：偏置场校正应在去噪之后、重采样之前进行，避免插值引入的伪影被误当作偏置场拟合。

三、重采样与 ROI 提取

重采样解决两个问题：一是把各向异性体素（如 CT 常见的 0.7×0.7×3.0 mm）转为各向同性，二是把不同模态数据统一到同一网格。推荐目标间距 1.0 mm 各向同性，插值方式按用途选择：CT 骨结构用线性插值即可，MR 与超声建议用 B 样条插值以保留灰度连续性。重采样后必须同步更新仿射矩阵，否则后续配准的物理距离会整体缩放。

ROI 提取的目的是减少无关背景对相似性测度的干扰。以脊柱导航为例：先在 CT 上用 200 HU 阈值提取骨结构，做 3×3×3 形态学闭运算填补骨小梁间隙，再取最大连通域并向外膨胀 10 mm，得到包含椎体与部分软组织的配准区域。这样既保留了骨性标志，又避免把手术床、器械等非解剖结构纳入优化。ROI 一旦确定，应同时作用于固定图像和浮动图像，保证两者空间范围一致。

四、数据同步与时间戳对齐

术中超声、电磁定位、光学跟踪等设备各自带有时钟，若直接按采集顺序配对，会出现几十毫秒到数百毫秒的错位。工程做法是：以系统主时钟为基准，记录每帧数据的采集时间戳（如 t_us = 1234567890.123 s），再插值到统一的 30 Hz 时间栅格上。对于跟踪数据，可用线性插值；对于超声图像，则按最近邻帧选取，避免图像插值产生虚假结构。

对齐后应做一致性检查：计算相邻帧时间间隔的均值和标准差，若标准差超过 50 ms，说明存在丢帧或时钟漂移，需要重新采集或对该段数据降权。最终输出的配准输入应满足三个条件：几何一致（同一仿射、同一间距）、灰度可比（去噪与偏置场校正完成）、时间对齐（公共时间栅格、超差帧已剔除）。

工程提示：预处理顺序建议固定为“读取与几何解析 → 去噪 → 偏置场校正 → 重采样 → ROI 提取 → 时间对齐”。其中去噪与偏置场校正的顺序不可颠倒，重采样应放在 ROI 提取之前，以便 ROI 边界按物理距离而非体素索引定义。

常见错误：① 未应用 RescaleSlope/Intercept，导致 CT 值范围错误；② 重采样后忘记更新仿射矩阵；③ 对超声图像做时间插值，引入不存在的解剖结构；④ ROI 只作用于浮动图像，造成相似性测度偏置。

完成上述步骤后，可得到一份可直接送入配准算法的数据集。建议在流水线末端保存一份预处理日志，记录每步的参数（滤波核、迭代次数、目标间距、ROI 边界、时间栅格）和输出数据的校验值（尺寸、间距、仿射矩阵、时间戳范围），便于术中复现和问题回溯。