
在上一节中,我们详细探讨了工业CT图像处理领域中普遍存在的噪声类型、其成因以及相应的解决策略。继噪声问题之后,本节我们转向CT图像处理中的另一个常见且复杂的问题——伪影(Artifacts)。
就工业CT扫描而言,被测物品通常由金属材料组成,例如最基础的铝制品与钢制品,而金属由于其材料特性而产生的金属伪影,相较于医学CT扫描,对于图像处理方面则是一个更为典型与常见的问题。
总的来说,金属伪影(Metal artifacts)的产生与金属材料的高衰减特性紧密相关,其成因复杂,可根据其产生原理分别归因于:
射束硬化效应(Beam hardening)
部分容积效应(Partial volume effect)
光子饥饿效应(Photon starvation)
X射线的散射效应(Scatter artifacts)等
根据上述四种不同的成因,本文将就每一种成因分别展开,并逐一给出图片示例。(注:下列CT图片均为在RETOMO软件中打开的效果)
1. 射束硬化效应导致的杯状伪影
图1 金属中心部的典型杯状伪影
2. 部分容积效应导致的金属伪影
图2 金属处的条纹状与放射状伪影
3. 光子饥饿效应导致的条纹状伪影
当X射线穿过高衰减度材料如金属时,材料会大量吸收X射线光子,导致较少的光子到达探测器,从而降低了图像相应部分的信噪比。低信噪比会在图像中产生条纹状阴影,特别是沿X射线的入射方向更为质密的区域,所以该种条纹状阴影经常出现在材料最质密处或处在较为中心处的高衰减材料的边缘处,如下图3展示了一组较为典型、出现在铝制铸件上的条纹状阴影。
图3 由于光子饥饿效应产生的条纹状阴影
4. 高衰减度材料的强散射效应导致的散射伪影
该种伪影是由于X射线束穿透高衰减度材料时,发生的较强的康普顿散射或相干散射效应。因此散射光子被探测器误判为射线的原始(直线)路径,从而导致图像上出现 “云雾“状伪影,这种伪影通常出现在高衰减度材料的边界处,会严重影响该种材料的边界识别,如下图4展示的铝制活塞的钢销边界处的典型云雾状散射伪影。
图4 典型的云雾状散射伪影
回顾上篇中我们提到的CT图像成像过程,虽然金属伪影在工业CT中最为常见,但都可以归纳为由于被扫描物体其自身材料性质所产生的伪影,仅发生在射线透射被测物体并达到探测器这一过程。除此之外,在CT成像的其他阶段如在发射阶段使用的不同射线类型;探测阶段所使用的不同检测设备或检测方法;不同仪器厂家在投影重构阶段所采用的不同算法;抑或是整个测试过程中由于操作产生的误差等等,都会在不同阶段产生不同的伪影。这些伪影往往具有一定的相似性,因此在下文中,我们将对发生在其他阶段的较为典型的伪影,就其成因、表现形式等方面进行展开,并提供示例以说明每种伪影的种类:
采样不足导致的阶梯状伪影(Step Artifacts)
当投影数据量不足或投影角度间隔过宽,则会产生此种伪影(Step Artifacts)。这些伪影通常表现为重建图像难以精确地捕捉材料边缘或细节处的信息,通常会表现为视图混叠现象。
图5 典型的由于采样数量不足而导致的阶梯状伪影
莫尔纹伪影(Moiré Artifacts)
图6 周期性噪声所导致的典型莫尔纹伪影
环状伪影(Ring Artifacts)
图7 两处典型的环状伪影
解决方案
由于伪影产生的特殊性,在BETA CAE开发的RETOMO软件中可以通过以下操作,减少伪影在图像分割中的影响。
图8 拖动绝对阈值以适应材料边界
RETOMO中还为用户开发了一些较为便捷的半自动图像分割工具,比较典型的可用于处理伪影导致的图像分割错误的工具如:
图9 孤岛算法位置
除了上述两种方法,在RETOMO中还有其他几种工具,如形貌学过滤器(Morphological filter)或者基于局部种子生长逻辑的轨迹分割工具(Trace Segmentation)或法兰工具(Flange Tool)均可用于处理由于伪影所导致的分割错误问题。此外,RETOMO中还支持更为快捷的人工智能复制分割,能大幅提升用户的使用效率。
至此,工业CT的应用场景、原理、常见问题分析的介绍就完成了。关于BETA RETOMO软件在这一领域的更多应用,欢迎随时与我们联系。
部分参考文献:
1. De Man B. Iterative reconstruction for reduction of metal artifacts in computed tomography[J]. 2001.
2. 丁忠明,张立凯.工业CT扫描文物的伪影与校正[J].文物保护与考古科学, 2022, 34(6):35-44.
3. Li B, Mo Y, Zhang P, et al. Streak artifact correction in X-ray CT imaging[J]. Atomic Energy Science and Technology, 2017, 51(11): 2075-2078.
END
作者 | 高凌宇
BETA CAE 工程师