跳至正文
仇天赐.研究与工程笔记
搜索标题、摘要和标签 Esc 关闭
返回技术档案

CT 残影排查:先找到信号在哪一层出错

从投影序列、时间响应和重建切片出发,区分时间滞后、运动与其他伪影,建立可验证的 CT 残影排查路径。

回顾整理

本文归入 2024.11 主题档案,整理于 2026.11.05。月份用于组织阅读路径,不代表具体实验或项目的完成日期。

把「残影」拆成可以验证的问题

这篇回顾把 CT 残影处理里需要先想清楚的问题拆开。归档月份用于组织主题,不代表当月完成了对应实验。

「轮廓后面有拖尾」只描述外观,不能直接定位原因。时间滞后、运动、金属相关误差,甚至显示窗设置,都可能让人产生相近的观感。第一步应保留原始输入,写清异常出现的层级:连续投影、正弦图、重建体,还是导出的显示图。缺少哪一层数据,也要写进问题记录。

用证据区分几类候选原因

探测器时间滞后(Detector Lag)意味着先前输入的信号进入后续帧。针对非晶硅平板探测器的研究说明了这种现象及其在锥束 CT 中造成伪影的可能性;它不能证明任何一份未知数据都存在同样机制。Starman 等,2012

可先建立一个排查表,每项只记录「观察到什么」和「还缺什么证据」。

候选原因 优先查看的证据 需要避免的推断
时间滞后 输入停止后是否仍有随时间衰减的响应 看到双边缘就认定为滞后
运动 相邻投影中的结构位置与采集时序 把角度变化当成物体运动
金属相关误差 异常与高衰减路径的关系 所有条纹都归为同一原因
处理或显示问题 中间结果、数值范围与固定显示窗 只比较自动拉伸后的截图

这张表是排查顺序,不是仅凭外观完成分类的规则。多种因素也可能同时出现。

让时序检查具有可比性

若具备受控阶跃输入,可在固定感兴趣区域(ROI)中记录每帧均值 m_t。以前段稳定曝光均值 m_on 和独立背景估计 b 定义归一化尾迹:

r_t = (m_t - b) / (m_on - b)

只有在分母远离背景噪声、曝光结束时刻已知时,这条曲线才有解释价值。必须保留实际时间间隔,不能把不等间隔采集的帧编号直接当作时间。曝光条件、背景估计方式和 ROI 位置也应一并保存。

没有受控序列时,相邻帧差分仍能帮助定位异常,却不能单独鉴别探测器记忆效应,因为真实场景本来就在变化。

从投影域追到重建域

建议固定显示窗,把原始投影、候选校正投影和重建切片并排检查。投影差值适合发现校正影响了哪些位置;重建结果用于观察这些变化是否沿扫描几何扩散。

还要检查远离明显伪影的平稳区域。若目标边缘更干净,但平稳区域的均值明显漂移,方案可能只是重新分配误差。能使用同一份原始投影做对照时,应固定重建参数,避免把滤波核变化带来的平滑误认为校正收益。

每次排查至少留下异常区域、正常区域和边缘剖面。它们比一张最漂亮的结果图更容易暴露副作用。

失败时回到假设,而非继续加滤波

常见失败是先做空间平滑,再尝试估计时间响应:噪声变小了,细结构也被改写,原本要解释的现象反而不清楚。另一种失败是用单张重建切片反推完整采集机制,证据不足却给出确定结论。

更稳妥的输出是一份分层结论:哪些异常已经复现,哪些候选原因得到支持,下一次最小对照需要什么数据。如果只有重建图,应明确保留因果判断的不确定性。排查的目标是缩小问题范围,让下一步校正拥有可检验的前提。

参考与阅读边界

笔记至此,探索继续。交流与纠错