X img

官方微信

img

群号:冰川冻土交流群

QQ群:218834310

  • 中国百强报刊
  • 中国精品科技期刊
  • 中国国际影响力优秀学术期刊
  • 中文核心期刊要目总览
  • 中国科技核心期刊
  • 中国科学引文数据库源刊
  • CN 62-1072/P 
  • ISSN 1000-0240 
  • 创刊于1979年
  • 主管单位:中国科学院
  • 主办单位:中国科学院寒区旱区
  •                  环境与工程研究所
  •                  中国地理学会
高级检索
作者投稿 专家审稿 编辑办公 编委办公 主编办公

冰川冻土, 2020, 42(4): 1407-1416 doi: 10.7522/j.issn.1000-0240.2018.1144

冰冻圈技术

冻土CT图像中环形伪影和射线束硬化伪影的成因及校正方法

陈世杰,1,2, 马巍1,2, 李国玉,1,2

1.中国科学院 西北生态环境资源研究院 冻土工程国家重点实验室,甘肃 兰州 730000

2.中国科学院大学,北京 100049

Ring artifact and beam hardening artifact in CT image of frozen soil: causes and correction methods

CHEN Shijie,1,2, MA Wei1,2, LI Guoyu,1,2

1.State Key Laboratory of Frozen Soil Engineering,Northwest Institute of Eco-Environment and Resources,Chinese Academy of Sciences,Lanzhou 730000,China

2.University of Chinese Academy of Sciences,Beijing 100049,China

通讯作者: 李国玉, 研究员, 从事冻土工程与寒区环境研究. E-mail: guoyuli@lzb.ac.cn.

编委: 武俊杰

收稿日期: 2018-09-07   修回日期: 2018-11-28  

基金资助: 国家自然科学基金项目.  41801038.  41630636
国家重点研发计划项目.  2017YFC0405101
甘肃省科技计划项目.  18JR3RA376

Received: 2018-09-07   Revised: 2018-11-28  

作者简介 About authors

陈世杰(1984-),男,甘肃高台人,工程师,2010年在兰州交通大学获硕士学位,从事冻土CT扫描试验和研究工作.E-mail:chsj53719@lzb.ac.cn , E-mail:chsj53719@lzb.ac.cn

摘要

冻土CT图像中不可避免地会出现环形伪影和射线束硬化伪影。为解决这些问题, 分析了两类伪影的特征和形成机理, 提出了降低伪影的硬件优化方案。首先, 通过改变X射线管电压, 以此改变X射线光通量来降低冻土CT图像的环形伪影。其次, 利用不同厚度的滤过材料, 将X射线中的低能射线预先过滤, 使得穿过冻土试样的X射线能谱范围变窄, 以此来降低射线束硬化伪影对试样的影响。对比分析结果表明: 将扫描电压提高到120 kV既能有效降低同类型冻土试样CT图像中的环形伪影, 又能保证图像的质量; 选用材质为2A12、 厚度为18 mm的航空铝材, 在扫描电压为120 kV、 扫描电流为230 mA的条件下可以有效抑制射线束硬化所造成的伪影。

关键词: 冻土 ; 环形伪影 ; 射线束硬化伪影 ; CT图像

Abstract

Ring artifacts and beam hardening artifacts are unavoidable in CT image of frozen soil. In order to remove the two kinds of artifacts, in this paper the characteristics and formation mechanisms of the two kinds of artifacts are analyzed and a hardware optimization solution for reducing artifacts is proposed. Firstly, the ring artifacts are reduced by increasing the X-ray luminous flux to change the X-ray energy spectrum. Secondly, the low energy rays in X-rays are pre-filtered by the filter materials of different thickness, so that the X-ray energy spectrum range across the frozen soil specimen is narrowed to correct the beam hardening artifacts. The test results show that increasing the scanning voltage to 120 kV can not only effectively reduce the ring artifacts in CT images of frozen soil, but also ensure the quality of the images. At the same time, the aluminum material of 2A12 with 18 mm thickness can effectively correct the artifacts caused by the beam hardening under the scanning voltage of 120 kV and scanning current of 230 mA.

Keywords: frozen soil ; ring artifact ; beam hardening artifact ; CT image

PDF (1760KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陈世杰, 马巍, 李国玉. 冻土CT图像中环形伪影和射线束硬化伪影的成因及校正方法[J]. 冰川冻土, 2020, 42(4): 1407-1416 doi:10.7522/j.issn.1000-0240.2018.1144

CHEN Shijie, MA Wei, LI Guoyu. Ring artifact and beam hardening artifact in CT image of frozen soil: causes and correction methods[J]. Journal of Glaciology and Geocryology, 2020, 42(4): 1407-1416 doi:10.7522/j.issn.1000-0240.2018.1144

0 引言

20世纪80年代, 医用CT扫描技术首次被Vinegar等1学者引入到岩土领域, 之后迅速发展, 成为岩土体细观结构变化观测的有效手段之一。在国内, 冻土工程国家重点实验室首次利用CT扫描技术对冻土试样冻融前后和受力条件下内部细观结构演化规律进行研究2-3, 发现扫描得到的冻土CT图像中存在环形或圆弧形的带状区域, 这种带状区域被称为环形伪影4-5, 整个环形或圆弧形的伪影通常与冻土的剪切带非常相似, 会对CT图像的分析造成错误的解读。为了消除环形伪影, Jin等6通过比较预测的初始图像值与实际投影数据收集的CT值, 提出用OSEM算法对环形伪影区域图像进行校准, 使伪影图像的均匀性提高5%。Titarenko等7基于CT图像经过傅里叶变换后其正弦图像沿水平方向是光滑的假设, 找到一种解析公式对环形伪影进行抑制。此外, 一些学者试图通过处理水平方向伪影图像的正弦图, 利用低通滤波的方法滤出高频分量, 来对环形伪影进行消除8-11, 例如Yenumula等12将CT图像投影正弦图用hann滤波器进行滤波后重建, 以此来减少环形伪影状况。Sijbers等13、 张国强等14利用坐标转化将CT图像投影数据转化成极坐标之后, 通过低通滤波算法对伪影图像进行伪影消除。Sadi等15基于投影数据平均曲线中环伪影的特征, 提出一种迭代加权中值滤波算法, 计算了中值滤波的中心权重, 从均值曲线偏差计算的局部偏差分布确定了迭代算法的收敛判据, 进而修正环形伪影。另有一些学者基于特征识别的方式对环形伪影正弦图中的不一致因子进行识别, 之后通过分段多项式拟合的方式对图像进行校正, 以达到消除环形伪影的目的16-18。以上校正方法虽然对环形伪影有一定的抑制作用, 但是这些方法均是基于图像后处理方法对伪影区域进行模糊和均匀处理, 会对其内部结构做出匀质化处理, 掩盖复杂组分材料内部的真实细观结构。特别是对于冻土试样来说, 其内部各组分物质之间密度的差异较小, 使得其投影信号没有突出的非连续性, 这些误差产生的图像信息也就没有清晰的边界, 利用上述方法进行校准时很难识别伪影的位置, 使用滤波的方式处理图像伪影时其阈值、 子集等相关参数较难确定。

同时, 也发现冻土CT图像边缘位置比中心位置的图像要亮, 在图像的中心位置处形成一个低密度带, 通常这种现象是由射线束硬化伪影造成的19-20。这种射线束硬化伪影也必然会给后期CT图像的解读带来严重的误导, 如在冻土的水分迁移过程中, 这种硬化伪影若不消除, 会得出水分由中心位置向试样边缘迁移的结论, 这与水分实际迁移机理是相悖的。消除射线束硬化伪影的主要方法可归纳为多项式拟合法和后处理重建算法21, 其中多项式拟合法是通过楔形模型利用多项式拟合的方法, 将射线穿透距离和多能量投影值进行线性拟合, 来达到校正目的22-25。但是这种方法要获得扫描材料的矿物组分以及穿透距离, 对于冻土材料有很大的局限性。后处理重建算法是通过对CT图像的两次数据重建, 将原始的投影数据从多色非线性转换为X射线在不同距离内通过物体的衰减和距离之间的单色线性关系, 利用正弦图像的投影数据估计X射线上高密物质的分布, 用该分布来对投影数据进行校正26-31。但这种方法需要预先知道扫描材料的能谱分布和空间分布, 这对于冻土复杂的内部结构来说是不可能实现的。

为了能够抑制冻土CT图像中的环形伪影和射线束硬化伪影, 本文分析了冻土CT图像中环形伪影和射线束硬化伪影的特征及产生机理, 在此基础上, 提出了降低以上两类伪影的硬件优化方案, 通过对比试验结果, 获得了最优的试验方案, 从而达到提高冻土CT图像质量的目的。

1 测试仪器及试样制备

本次试验所用的CT设备为中国科学院冻土工程国家重点实验室引进的荷兰飞利浦Brilliance 16排多能量螺旋医用CT机, 如图1所示。其空间分辨率为24 lp·cm-1, 密度分辨率为0.3%, 扫描电压为90~140 kV, 最大扫描电流为500 mA。

图1

图1   CT试验设备

Fig.1   The CT equipment used in the test


本次试验用土采自连霍高速甘肃省永登县某取土场附近, 按照《土工试验方法标准》测得土样的颗粒组成及物性指标, 如表1所示。将土制备成初始含水量为13%(最优含水率)、 干密度为1.82 g⋅cm-3、 压实度为95%的ø61.8 mm×125 mm圆柱状试样, 放置在-30 ℃冷库中快速降温6 h, 之后将冻土样放置在温度为-6 ℃的恒温箱中恒温24 h后, 开始CT扫描试验, 取中间层进行分析。

表1   土样的基本物理指标

Table 1  Basic physical parameters of the soil specimen

颗粒组成(质量占比)/%液限/%塑限/%比重最优含水率/%最大干密度/(g⋅cm-3
>0.05 mm0.05~0.005 mm<0.005 mm
15661926182.7131.9

新窗口打开| 下载CSV


2 环形伪影

2.1 环形伪影的特征及产生机理

制备好的冻土试样在90 kV扫描电压下的CT图像如图2所示。由图2(a)中标示区域可以看到CT图像中沿中心位置分布的半圆状环形伪影, 为了便于进一步比较, 将伪影CT图像进行傅里叶变换得到傅里叶投影图[图2(b)], 可以看到环形伪影在变换后的图像中呈现出放射状, 并且亮度不均匀。

图2

图2   冻土CT图像中的环形伪影

Fig.2   Ring artifacts corrupted in a CT image of frozen soil: CT image at 90 kV scanning voltage (a) and its projection image after Fourier transform (b)


按照CT图像的成像原理, CT图像的重建一般采用反投影重建算法, 在连续情况下的图像中某一点的反投影重建公式为32

f(r,θ)=1π0πpφrcos(θ-φ)dφ

式中: (rθ)为空间任意位置的极坐标; φ为不同投射方向与直角坐标的夹角; frθ)为经过给定点(rθ)的所有射线投影在φ=0~π区间内累加所得到的(rθ)点的像素值; pφφ方向上的投影值。

式(1)可知, 扫描过程中射线束的硬化、 探测器的响应不一致等引起的误差均会对投影值pφ 带来影响, 同时由于反投影重建算法中存在积分项, 这会使极小误差值在反投影重建过程中被放大, 从而影响最后重建得到的CT图像中的各个像素值, 使像素值的误差增大, 从而影响最终的成像结果。

有研究表明, 环形伪影的形成通常是由探测器通道错误或是探测器通道误差过大造成的33-34。为了表征它们之间的关系, 将不同探测器通道投影误差值设为εi, 由CT成像原理可知, 反投影图像形成环形, 与发生错误的探测器通道的投影线相切, 其结果是环形伪影的强度沿着整个环形的圆周积分Qi, 它等于探测器通道投影误差值εi 与投影数n的积, 由此可以得到CT图像中总的环形伪影强度C

C=1iQi2πri=n1iεi2πri

式中: C为总的环形伪影强度; i为探测器所对应的环形伪影个数; ri 为环形伪影半径; Qi 为半径为ri 的环形伪影强度的圆周积分; n为投影数; εi 为探测器通道投影误差值。

式(2)可以进一步验证引起冻土CT环形伪影的原因有两个: ①探测器通道误差偏大, ②探测器通道的响应不一致。

2.2 环形伪影的校正

为了解决冻土CT图像中的环形伪影问题, 首先要判定环形伪影是由探测器通道误差偏大所造成的, 或是由探测器通道的响应不一致造成的, 为此利用CT机初始校准用的标准物质(水膜)来进行判断, 通过水膜扫描得到CT图像和中心剖面线位置曲线如图3所示。

图3

图3   水膜原始CT图像及中心剖面线CT值分布

Fig.3   Original CT image of water phantom (a) and its CT number profile along the central section line (b)


图3可以看到, 水膜的CT图像并未出现肉眼可观测到的环形伪影, 同时中心剖面线CT值分布也基本趋于一条直线, 并未出现较大误差波动。基于以上原因, 可以排除冻土CT图像中的环形伪影是由探测器通道误差偏大所造成的。

但是在随后对同等扫描条件下的冻土试样扫描时发现, 其环形伪影依然存在。εi 引起的标准差为

SDεi=i=1N(εi-εi¯)2

为了减少上述统计标准差SDεi, 大多数医用CT生产厂商使用与人体密度(约1.06 g⋅cm-3)较为接近的标准物质水膜(密度约1 g⋅cm-3)校正来补偿SDεi, 通过水膜平均得到增益校正值来降低SDεi, 同时对探测器响应不一致的通道通过暗场和增益进行校正。

但是不同于水膜这种单一物质, 冻土是由冰颗粒、 未冻水、 土颗粒、 空气组成的多相特殊土, 其密度范围在1.4~2.6 g⋅cm-3之间, 同时冻土试样的密度、 相对分子质量、 样品截面尺寸、 形态与纯水等几何校准模型不一致, 使得照射到探测器单元的X射线光通量不一致。另外, 不同于人体CT扫描, 在利用医用CT进行冻土扫描时为了获得更高的图像质量, 会适当修改扫描参数以提高放射剂量, 这均会导致在原有标准物质(水膜)下得到的校准值超出了医用CT常规剂量探测器单元的校准范围, 使得校准值并不能很好地补偿探测器增益, 进而导致一定统计水平下的光通量通过探测器通道而出现的响应不一致的状况, 这也是经过水膜校准后的冻土试样环形伪影依然存在的原因。要解决这一问题, 要改变穿过冻土试样的X射线探测器单元的光通量, 提高探测器的光子统计水平, 使其范围控制在医用CT的校准范围之内。

依据X射线物理学原理, 连续X射线有一个波长极限, 它是由光子一次碰撞就耗尽能量所产生的X射线, 只与扫描电压有关, 并且有如下关系式。

ε=eU=hcλ0
εp=hλ0

式中: ε为光子的能量; e为电子电荷; U为扫描电压; h为普朗克常数; c为X射线传播速率; λ0为X射线波长; εp为动量。

式(4)和式(5)可知, 动能越大, 其穿透能力越强, 同时X射线波长越短动能越大, 要改变穿过冻土试样的X射线探测器单元的光通量只需要提高扫描电压就可以了。

2.3 校正结果

依据上述机理, 本文通过适当提高扫描电压从而改变X射线波长极限, 对冻土CT图像中的环形伪影进行硬件校准, 使得穿过冻土试样的X射线探测器单元的光通量能够满足原有增益校正值的校正范围, 从而补偿SDεi, 进而解决探测器响应不一致所造成的环形伪影。为此, 针对同一个冻土样品开展了三组不同扫描电压参数的CT试验, 为了保证每次扫描时的球管曝光功率近似一致且不超过机器的额定功率, 在改变不同电压参数的同时也适当地调整扫描电流参数, 得到了90 kV、 120 kV、 140 kV下的冻土CT图像和傅里叶变换投影图(图4), 其扫描参数及感兴趣区域的CT均值、 标准差、 图像信噪比(signal noise ratio, SNR)评价等如表2所示。

图4

图4   不同扫描电压条件下的CT图像

Fig.4   CT images under scanning voltage of 90, 120 and 140 V (a, c, d) and their projection images after Fourier transform (b, d, f)


表2   不同扫描条件下的CT图像参数

Table 2  The CT image parameters under different scanning condition

扫描条件感兴趣区域CT均值标准差信噪比环形伪影情况图像质量评价
扫描电压/kV扫描电流/mA层厚/mm曝光时间/ms
9031030.751 443.5198.57.14严重较差
12023030.751 230.1111.811.11降低良好
14020030.751 147.6138.18.33降低适中

新窗口打开| 下载CSV


为了对不同扫描电压下的图像质量进行量化对比, 可以通过信噪比(SNR)来对图像质量进行评价。信噪比越大其图像质量越高, 信噪比越小其图像噪声越大, 图像质量越低。SNR定义为12

SNR=σROIμROI

式中: SNR为CT图像中感兴趣区域的信噪比; μROI为感兴趣区域的标准差; σROI为感兴趣区域的CT均值。

通过对不同电压下的CT图像对比发现, 在扫描电压为90 kV时冻土试样的CT图像中存在环形伪影, 该组电压值也是利用标准物质校准时用到的扫描电压。当校准扫描电压提高到120 kV和140 kV时, 其环形伪影得到了很好的抑制, 傅里叶变换投影图中的线状伪影也全部消除。但是对比上述两组傅里叶变换投影图发现, 120 kV傅里叶变换投影图中的明、 暗区域较140 kV的投影图更加均匀, 并且噪点较少。另外, 对比图4中所选择感兴趣区域的图像信噪比, 发现扫描电压为120 kV时其信噪比最大, 图像质量最高。因此, 在扫描电压为120 kV的这组参数下, 既能减低探测器响应不一致所造成的环形伪影, 又能保证CT图像的质量, 该组扫描参数是冻土CT图像抑制环形伪影保证图像质量的最优扫描条件。

3 射线束硬化伪影

3.1 射线束硬化伪影的特征及产生机理

冻土CT图像的射线束硬化伪影具体表现为CT图像中心部分的图像亮度比外部边缘处暗, 在边缘区域形成一个圆形亮带, 同时从中心剖面线分布可以看到中心部位的CT值要小于试样边缘处的CT值, 在边缘位置出现两个峰值, 整个剖面线分布呈凹型, 如图5所示。

图5

图5   冻土CT图像中的射线束硬化伪影及中心剖面线CT值分布

Fig.5   Beam hardening artifacts in a CT image of frozen soil (a) and its CT number profile along the central section line (b)


按照Lambert-Beer定理, 单能X射线穿透均匀介质后其入射线和出射线的强度呈指数衰减规律。

I=I0e- μΔx

式中: II0分别为X射线出射和入射强度; μ为物质的质量吸收系数; Δx为穿过物体的厚度。

式(7)可得

μΔx=- lnII0=p

式中: p为X射线方向的投影值。

式(8)可知, 单能X射线穿透均匀介质后沿着射线方向的投影值p与穿过物体的厚度Δx呈线性关系, 其关系如图6所示。

图6

图6   投影值与材料厚度关系

Fig.6   Relationship between projection value p and material thickness Δx


但所有CT设备产生的X射线均有很宽的能谱值及波长范围, 其中X射线的波长范围在10-8~10-12 m, 也就是说实际中X射线是多能的, 这是CT设备所固有的属性, 无法改变。对于非均匀物质来说, 有

μ=1n(μmiρi)Wi

其中

μmiKλiZ3

式中: μmi 为单位长度的X射线质量吸收系数; ρi 为物质密度; Wi 为元素质量百分比; K为常数; λi 为不同能量的X射线波长; Z为靶材的原子序数。

由式(9)~(10)可知, 对于同一种物质, 当穿过的X射线波长λi 不同, 其衰减系数μ也不同, 若仍然用单能衰减进行补偿, 则会出现同一投影值p1对应不同厚度Δx1和Δx2的现象。同时, 波长较长的低能“软射线”先被吸收, 波长较短的高能“硬射线”得以穿透, 使得中心位置表现出“低密度性”, X射线穿透均匀介质后沿着射线方向的投影值p与穿过物体的厚度Δx呈非线性关系, 如图6所示。

正是基于以上原因, 要根本避免产生射线束硬化伪影, 其理想状况是采用单能的X射线源, 这可以使非线性问题得以解决, 但是这种单能X射线放射源造价昂贵, 在实际应用中并未得到广泛普及, 尤其在医用CT机领域尚未出现单能X射线球管。这就使得未校准的冻土CT图像均存在射线束硬化伪影现象。

3.2 射线束硬化伪影的校正

如果对穿过物质的多能射线束进行滤波, 其波长覆盖范围必然变窄, 在窄波长范围内可以近似认为X射线是单能的33。基于以上思路, 可以在样品周围加装滤过装置, 如图7所示, 滤过装置会优先吸收和过滤波长较长的低能X射线, 使得射线束穿过样品之前的射线均是短波高能射线, 而低能射线在未穿过试样之前已被滤过装置所吸收, 这就使得最终穿过试样和探测器的X射线波长范围变窄。

图7

图7   滤波示意图

Fig.7   Schematic diagram of the filter process


考虑到医用CT机的结构及试样安装拆卸的便利性及经济性, 最终设计的X射线束滤过装置如图8所示。滤过装置设计为圆柱状, 冻土样品放置在滤过装置中心位置, 试样和滤过装置之间周围用吸收系数极低的保温塑料泡沫材料填充。考虑到试验所用的CT机为医用CT, 其激励电压及X线管功率相对较低, X射线穿透能力相对较弱, 若选用原子序数较大的滤过材料如铅、 铜等, 会使得X射线衰减过快影响成像。综合考虑上述因素, 本次试验滤过装置的材料选用吸收系数相对较小的铝材, 其标号为2A12, 化学成分为: Si 0.5%, Fe 0.5%, Cu 3.8%~4.9%, Mg 1.2%~1.8%, Mn 0.3%~0.9%, Zn 0.3%, Ni 0.1%, 其余均为Al。为了降低射线束硬化伪影, 同时还要保证CT图像质量, 滤过装置的壁厚选择要综合考虑图像的噪声、 CT值的误差范围等因素, 为了确定滤过装置的最优厚度和验证校准效果, 开展了5组同等扫描条件下不同厚度的CT试验。扫描条件为扫描电压120 kV, 扫描电流230 mA, 层厚3 mm, 曝光时间0.75 ms, 材料厚度分别为8 mm、 13 mm、 18 mm、 21 mm、 25 mm。试验选取试样的中间位置作为分析层, 其层位数为16, 为了保证每组CT断层图像的一致性, 试样所选用的样品均为同一冻土试样, 并且每次扫描时对同分析层也做了定位, 这样就保证了对比分析的图像为同一试样的同一层位, 避免由于样品差异和定位不同而造成的误差, 其定位图像如图9所示。

图8

图8   滤过装置示意图

Fig.8   Schematic diagram of the filter device


图9

图9   样品定位示意图

Fig.9   Schematic diagram of the specimen positioning


3.3 校正结果

图10为经过不同厚度滤过材料校正后的冻土试样的CT图像及中心剖面线CT值分布图。由图可知, 在未加滤过材料和滤过材料的厚度分别为8 mm、 13 mm时, 冻土试样的射线束硬化伪影虽然有所降低但并未消除, 中心剖面线所示边缘区域的CT值仍然大于中心区域(图中标示位置所示), 其滤波校正效果不明显。随着滤过材料厚度的增加, CT图像边缘区域和中心区域亮度逐渐均匀并趋于一致, 边缘位置高亮环形区域消失, 同时中心剖面线边缘到中心位置的CT值波动较小, 也未出现边缘CT值大于中心CT值的情况。但是对比18 mm、 21 mm、 25 mm厚度的CT图像发现, 随着滤过材料厚度的增加, 其对应的CT图像损失的细节越大, 图像噪声也越来越大。

图10

图10   射线束硬化伪影校正结果

Fig.10   Results after beam hardening correction: (a), (c), (e), (g), (i) and (k) are, respectively, the CT images under 0, 8, 13, 18, 21 and 25 mm filtering; (b), (d), (f), (h), (j) and (l) are, respectively, their CT number profiles along the central section line


为了进一步对射线束硬化校准结果的图像质量进行量化表征, 选用式(6)中的信噪比参数对CT图像质量进行评价, 材料厚度与图像质量评价参数SNR关系曲线如图11所示, 材料厚度、 感兴趣区域CT值及标准差等参数如表3所示。

图11

图11   材料厚度与图像质量评价参数SNR的关系

Fig.11   Relation between material thickness and image quality evaluation parameters (signal noise ratio, SNR)


表3   不同厚度滤过材料下的图像质量比较

Table 3  Comparison of image quality under different filter material thickness

滤过材料材质滤过材料厚度/mm感兴趣区域面积/cm²感兴趣区域CT均值标准差信噪比环形伪影情况图像质量评价
029.281 576.2124.612.65严重较差
2A12829.281 434.596.214.91严重较差
2A121329.281 369.184.416.22严重适中
2A121829.281 291.174.517.33降低良好
2A122129.281 230.786.214.28降低较差
2A122529.281 216.5109.811.01降低较差

新窗口打开| 下载CSV


对比表3图11中的SNR发现, 厚度为18 mm时图像质量的信噪比最大, 同时也极大地降低了射线束硬化带来的影响, 由此滤过材料的最优校准厚度确定为18 mm。此时既能有效抑制冻土试样的射线束硬化所造成的伪影, 又能很好保证CT图像质量。

4 结论

本文针对冻土CT图像环形伪影和射线束硬化伪影问题, 分析了两类伪影的特征及产生机理, 在此基础上提出了抑制伪影的硬件优化方案, 通过优化前后试验结果的对比分析, 得出以下结论:

(1) 冻土CT图像产生环形伪影的原因是冻土试样的密度、 分子量、 样品截面尺寸、 形态与标准校准物质(水膜)的几何校准模型不一致, 使得照射到探测器单元的X射线光通量不一致, 导致在原有标准物质下得到的校准值并不能很好地补偿探测器增益, 引起探测器通道的响应不一致, 进而产生环形伪影。通过对比试验测试, 适当提高扫描电压, 通过改变X射线波长可以降低探测器响应不一致所造成的环形伪影。对于冻土试样, 将扫描电压提高到120 kV既能有效降低CT图像中的环形伪影, 又能保证图像质量, 同时该组扫描参数也是抑制冻土试样CT图像环形伪影的最优扫描条件。

(2) 冻土CT图像产生射线束硬化伪影的原因是多能X射线在穿过试样时, 波长较长的射线(低能软射线)过早被吸收, 使得波长较短的射线(高能硬射线)比例增加, 从而产生射线束的硬化伪影。通过在样品周围加入一定厚度的滤过装置, 可以有效降低同等类型条件下冻土CT图像的射线束硬化伪影。通过对比测试, 在同等扫描条件下(扫描电压120 kV, 扫描电流230 mA, 层厚3 mm, 曝光时间0.75 ms), 选用材质为2A12, 厚度为18 mm的航空铝材可以有效抑制射线束硬化所造成的伪影。

参考文献

Vinegar H JWellington S L.

Tomographic imaging of three-phase flow experiments

[J]. Review of Scientific Instruments, 1987581): 96 - 107.

[本文引用: 1]

Chen ShijieZhao ShupingMa Weiet al.

Status and prospects of the frozen soil studies using CT technology

[J]. Journal of Glaciology and Geocryology, 2013351): 193 - 200.

[本文引用: 1]

陈世杰赵淑萍马巍.

利用CT扫描技术进行冻土研究的现状和展望

[J]. 冰川冻土, 2013351): 193 - 200.

[本文引用: 1]

Zhao YueHan QiaolingZhao Yandonget al.

The research status and analysis of frozen soil interior materials based on CT nondestructive scanning technology

[J]. Journal of Glaciology and Geocryology, 2017396): 1307 - 1315.

[本文引用: 1]

赵玥韩巧玲赵燕东.

基于CT无损扫描技术的冻土内部物质研究现状与分析

[J]. 冰川冻土, 2017396): 1307 - 1315.

[本文引用: 1]

Kyriakou YPrell DKalender W A.

Ring artifact correction for high-resolution micro CT

[J]. Physics in Medicine and Biology, 20095417): 385 - 391.

[本文引用: 1]

Anas E MLee S YHasan K.

Classification of ring artifacts for their effective removal using type adaptive correction schemes

[J]. Computers in Biology and Medicine, 2011416): 390 - 401.

[本文引用: 1]

Jin D S-CChou W-HChen J-C.

Application of sinogram-based moving window subtraction technique in OSEM algorithm to reduce ring artifacts for cone-beam micro-CT image quality assessment

[C]// Proceedings of the 2nd International Conference for Innovation in Biomedical Engineering and Life Sciences. SingaporeSpringer Singapore2018253 - 257.

[本文引用: 1]

Titarenko SWithers P JYagola A.

An analytical formula for ring artefact suppression in X-ray tomography

[J]. Applied Mathematics Letters, 20102312): 1489 - 1495.

[本文引用: 1]

Anas E M ALee S YHasan M K.

Removal of ring artifacts in CT imaging through detection and correction of stripes in the sinogram

[J]. Physics in Medicine and Biology, 20105522): 6911 - 6930.

[本文引用: 1]

Kowalski G.

Suppression of ring artefacts in CT fan-beam scanners

[J]. IEEE Transactions on Nuclear Science, 1978255): 1111 - 1116.

Raven C.

Numerical removal of ring artifacts in microtomography

[J]. Review of Scientific Instruments, 1998698): 2978 - 2980.

Münch BTrtik PMarone Fet al.

Stripe and ring artifact removal with combined wavelet: Fourier filtering

[J]. Optics Express, 20091710): 8567 - 8591.

[本文引用: 1]

Yenumula LAcharya R VKumar Uet al.

Ring artifact correction in gamma-ray process tomography imaging

[J]. Applied Radiation and Isotopes, 201712475 - 82.

[本文引用: 2]

Sijbers JPostnov A.

Reduction of ring artefacts in high resolution micro-CT reconstructions

[J]. Physics in Medicine and Biology, 20044914): 247 - 253.

[本文引用: 1]

Zhang GuoqiangZhou HuHe Youet al.

Ring artifacts correction of computerized tomography image based on polar-coordinate transform

[J]. Acta Optica Sinica, 2012325): 325 - 330.

[本文引用: 1]

张国强周虎和友.

基于极坐标变换去除计算机层析图像环形伪影

[J]. 光学学报, 2012325): 325 - 330.

[本文引用: 1]

Sadi FLee S YHasan M K.

Removal of ring artifacts in computed tomographic imaging using iterative center weighted median filter

[J]. Computers in Biology and Medicine, 2010401): 109 - 118.

[本文引用: 1]

Li BaoleiFeng YaLi Junjiang.

Research on beam hardening correction method for computed tomography imaging

[J]. Life Science Instruments, 200759): 34 - 38.

[本文引用: 1]

李保磊封亚李俊江.

CT成像射束硬化校正方法研究

[J]. 生命科学仪器, 200759): 34 - 38.

[本文引用: 1]

Boin MHaibel A.

Compensation of ring artefacts in synchrotron tomographic images

[J]. Optics Express, 20061425): 12071 - 12075.

Wei YingWu WenjinYang Min.

A correction method for ring artifacts in CT image based on polynomial fitting

[J]. Journal of Beijing Institute of Machinery (Natural Science), 2010254): 45 - 48.

[本文引用: 1]

魏英吴文晋杨民.

基于多项式拟合的CT图像环状伪影校正

[J]. 北京信息科技大学学报(自然科学版), 2010254): 45 - 48.

[本文引用: 1]

Ketcham R ACarlson W D.

Acquisition, optimization and interpretation of X-ray computed tomographic imagery: applications to the geosciences

[J]. Computers and Geosciences, 2001274): 381 - 400.

[本文引用: 1]

Zöllner CRit SKurz Cet al.

Decomposing a prior-CT-based cone-beam CT projection correction algorithm into scatter and beam hardening components

[J]. Physics and Imaging in Radiation Oncology, 2017349 - 52.

[本文引用: 1]

Jovanović ZKhan FEnzmann Fet al.

Simultaneous segmentation and beam-hardening correction in computed microtomography of rock cores

[J]. Computers and Geosciences, 201356142 - 150.

[本文引用: 1]

Ketcham R AHanna R D.

Beam hardening correction for X-ray computed tomography of heterogeneous natural materials

[J]. Computers and Geosciences, 20146749 - 61.

[本文引用: 1]

Zhang QuanhongLu HongnianYang Min.

Beam hardening correction using polynomial fitting based on reprojection

[J]. Optical Technique, 2005314): 633 - 635.

张全红路宏年杨民.

基于重投影的多项式拟合校正射束硬化

[J]. 光学技术, 2005314): 633 - 635.

Zeng GangYu ZhongqiangYan Yonglian.

Beam hardening correction based on Monte Carlo simulation

[J]. High Energy Physics and Nuclear Physics, 2006302): 178 - 182.

曾钢郁忠强阎永廉.

基于蒙特卡罗模拟的射束硬化校正方法

[J]. 高能物理与核物理, 2006302): 178 - 182.

Zhang YihaiZhang CuiPan Xiaodonget al.

The methods for beam hardening correction of cone-beam industrial computed tomography

[J]. Nondestructive Testing, 2017396): 8 - 12.

[本文引用: 1]

张益海张催潘小东.

锥束工业CT射束硬化校正方法

[J]. 无损检测, 2017396): 8 - 12.

[本文引用: 1]

Joseph P MSpital R D.

A method for correcting bone induced artifacts in computed tomography scanners

[J]. Journal of Computer Assisted Tomography, 197821): 100 - 108.

[本文引用: 1]

Meagher JMote CSkinner H.

CT image correction for beam hardening using simulated projection data

[J]. IEEE Transactions on Nuclear Science, 1990374): 1520 - 1524.

Remeysen KSwennen R.

Beam hardening artifact reduction in microfocus computed tomography for improved quantitative coal characterization

[J]. International Journal of Coal Geology, 2006671/2): 101 - 111.

Krumm MKasperl SFranz M.

Reducing non-linear artifacts of multi-material objects in industrial 3D computed tomography

[J]. Independent Nondestructive Testing and Evaluation, 2008414): 242 - 251.

Chen HuijuanPan Jinxiao.

A beam-hardening correction method based on poly-energetic statistics in X-ray CT

[J]. Computerized Tomography Theory and Applications, 2010191): 21 - 27.

陈慧娟潘晋孝.

一种基于多能统计的射束硬化校正方法

[J]. CT理论与应用研究, 2010191): 21 - 27.

Rasoulpour NKamali-Asl AHemmati H.

A new approach for beam hardening correction based on the local spectrum distributions

[J]. Nuclear Instruments and Methods in Physics Research: Section A Accelerators, Spectrometers, Detectors and Associated Equipment, 2015794177 - 184.

[本文引用: 1]

Zhuang Tiange. The principle and algorithm of CT[M]. ShanghaiShanghai Jiao Tong University Press1992.

[本文引用: 1]

庄天戈. CT原理与算法[M]. 上海上海交通大学出版社1992.

[本文引用: 1]

Yousuf M AAsaduzzaman M.

An efficient ring artifact reduction method based on projection data for micro-CT images

[J]. Journal of Scientific Research, 200921): 37 - 45.

[本文引用: 2]

Li JunjiangHu ShaoxingLi Baoleiet al.

Ring artifact correction for industrial CT images

[J]. Journal of Beijing University of Aeronautics and Astronautics, 20073311): 1378 - 1382.

[本文引用: 1]

李俊江胡少兴李保磊.

CT图像环状伪影校正方法

[J]. 北京航空航天大学学报, 20073311): 1378 - 1382.

[本文引用: 1]

/