采样、噪声与伪影
当理想 Radon 模型不再成立
角度欠采样、有限角、噪声、截断、束硬化和散射伪影。
1 伪影分类
重构伪影可以按来源分为三类:
| 类别 | 典型原因 | 常见表现 |
|---|---|---|
| 采样问题 | 角度过少、探测器间距过大 | 放射状条纹、混叠、分辨率损失 |
| 数据缺失 | 有限角、截断、坏像素 | 方向性模糊、边缘延伸、环状伪影 |
| 模型失配 | 多能谱、散射、运动 | 杯状伪影、暗带、重影、定量偏差 |
2 角度欠采样
当投影角数量不足时,反投影路径无法在图像域均匀覆盖,强边缘会产生星状条纹。经验上,直径为 \(N\) 个像素的圆形视野通常需要数量级为
\[ N_\theta \gtrsim \frac{\pi}{2}N \tag{1}\]
的独立角度,才能使角向采样与探测器采样大致匹配。式 1 是工程估计而非对所有物体都严格成立的定理。
3 有限角问题
如果只采集 \(\theta\in[\theta_0,\theta_1]\) 且角度范围小于 \(180^\circ\),二维频域中将出现未被径向切片覆盖的楔形区域。其结果通常是沿缺失方向的边缘无法稳定恢复。
4 噪声传播
光子计数近似服从 Poisson 分布:
\[ R\sim\operatorname{Poisson}(\bar R). \tag{2}\]
负对数变换使噪声变为信号相关;Ramp 滤波又会增强高频噪声。因此低剂量条件下常使用窗化 FBP、统计迭代重构或学习型先验。
5 探测器截断
当物体超出探测器视野时,投影两端被截断,傅里叶切片条件所需的完整线积分不再可用。直接 FBP 会在视野边缘产生亮边或暗带,并造成全局偏差。
6 束硬化
多能谱成像模型 不满足严格线性线积分。低能光子优先被吸收后,穿过厚区域的平均能量升高,使投影与路径长度之间呈非线性关系。均匀物体常表现为中心偏暗的杯状伪影,两个高衰减物体之间可能出现暗带。
7 环状伪影
探测器通道增益校正误差在每个角度中位于近似固定的 \(s\) 坐标。旋转重构后,这种固定列误差会形成以旋转中心为圆心的环。
重要诊断原则
不要把所有伪影都归因于重构算法。先检查原始投影和正弦图:固定列、角度突变、饱和区、截断边界和非线性趋势往往能直接暴露问题来源。