3. PET与CT图像预处理
PET/CT融合影像处理的第一步并不是配准或融合，而是把两台设备采集到的原始断层数据整理成“可比较、可计算、可融合”的状态。PET与CT在物理机制、空间分辨率、噪声特性和灰度含义上差异极大：PET反映放射性核素分布，计数低、噪声服从泊松统计、分辨率通常为4～6 mm；CT反映X射线衰减系数，噪声近似高斯分布、分辨率可达0.5～1 mm，灰度以HU为单位且具有物理意义。预处理的目标，就是让二者在几何、灰度和伪影层面尽量对齐，为后续配准、融合和定量分析提供稳定输入。
去噪与平滑：先分清噪声来源
PET图像的主要噪声来自低计数统计涨落，表现为颗粒状斑点；CT噪声则来自量子噪声和电子学噪声，表现为细密纹理。去噪不能“一刀切”，否则会抹掉PET的病灶摄取细节或CT的骨小梁结构。工程上常采用分层策略：对PET先做高斯平滑或边缘保持滤波，再进入重建后处理；对CT则多用中值滤波或非局部均值抑制脉冲噪声。
以PET为例，若原始重建图像分辨率为4 mm，常用高斯核σ=1.5～2.0 mm进行轻度平滑，使噪声标准差下降约30%～50%，同时保持病灶对比度损失在可接受范围。若使用各向异性扩散滤波，则需控制迭代次数（如5～10次）和梯度阈值，避免过度平滑导致小病灶“消失”。平滑不是越强越好，判断依据是：在均匀本底区域噪声明显降低，而在已知病灶或高摄取区边缘仍能保持清晰边界。
重采样：统一空间网格
PET与CT的层厚、像素间距和矩阵尺寸往往不同。例如CT常为512×512、层厚1 mm，PET常为128×128或192×192、层厚3～5 mm。融合前必须把二者重采样到同一空间网格。重采样包括插值和坐标变换两步：先确定目标网格（通常以CT为参考，因为CT分辨率更高、几何更稳定），再把PET映射到该网格上。
插值方法选择直接影响结果：最近邻插值保留原始值但会产生块状效应，适合标签图；三线性插值平滑但可能轻微模糊；B样条插值精度高但计算量大。实际工程中，PET到CT的重采样常用三线性或B样条，目标像素间距取CT的原始间距，层厚取1～3 mm。重采样后要检查：PET的SUV值是否因插值发生系统性偏移，通常要求均匀区域均值变化小于2%。
灰度归一化：让两种模态可比
CT的HU值有明确物理意义（水为0 HU，空气约-1000 HU），一般不做全局归一化，但需要做窗宽窗位调整以突出不同组织。PET的SUV虽然也是定量指标，但受注射剂量、体重、采集时间影响，不同患者之间不可直接比较。因此，灰度归一化主要针对PET：将SUV映射到固定范围（如0～10或0～15），或按本底摄取做标准化。
若后续使用深度学习模型，常把CT截断到[-1000, 3000] HU后线性映射到[0,1]，PET则用SUVmax或SUVpeak做归一化。注意：归一化参数必须从训练集统计得到并固定，不能每例单独重算，否则会破坏模型的一致性。对于融合显示，则常用伪彩映射，把PET高摄取区映射为暖色，CT保持灰阶，便于医生判读。
金属伪影与截断伪影的初步处理
金属伪影主要出现在CT中，由高密度植入物（如髋关节、牙科材料）引起射线硬化、散射和光子 starvation，表现为亮暗条纹。初步处理包括：提高管电压、增加迭代重建权重、使用金属伪影校正算法（MAR）。在工程流水线中，若无法重新扫描，可先对金属区域做阈值分割，再用插值或正弦图修复抑制条纹。判断标准是：金属周围软组织HU值波动减小，条纹伪影不再干扰PET衰减校正。
截断伪影则常见于PET/CT视野边缘，当患者体型超出扫描FOV时，体表数据被截断，导致CT值错误、PET衰减校正偏差。初步处理包括：使用扩展FOV重建、外推体轮廓、或采用基于PET发射数据的衰减校正。若截断严重，应在报告中标注，避免定量结果被误读。
预处理流水线设计
把上述步骤串成一条可复现的流水线，是工程落地的关键。推荐顺序为：原始数据读取 → CT去噪/金属伪影校正 → PET去噪/平滑 → 以CT为参考重采样PET → PET灰度归一化 → 截断伪影检查与补偿 → 输出统一网格数据。每一步都应记录参数和中间结果，便于回溯和质控。
流水线设计中有三个容易忽视的细节。第一，重采样必须在去噪之后、归一化之前，否则插值会把噪声扩散到相邻体素，后续平滑难以区分真实信号与插值伪影。第二，金属伪影校正应尽量在CT域完成，因为PET衰减校正依赖CT的HU值，CT伪影会直接传递到PET。第三，所有参数（滤波核、插值方法、归一化范围）必须写入配置文件并随病例保存，否则不同操作者得到的结果不可比。
最后需要强调：预处理不是“越干净越好”，而是“在抑制干扰的同时保留定量信息”。PET的SUV、CT的HU都是后续诊断和放疗计划的重要依据，任何滤波、插值、归一化都可能引入偏差。因此，每一步都应有明确的参数依据和质控指标，并在报告中记录关键处理步骤，使结果可追溯、可复现。