层面位置对西门子CT模拟机图像噪声的影响

蒋志玮1a,霍新颖1b,时飞跃1b,2,秦伟1b,林强1a,秦航1a

1. 南京医科大学附属南京医院(南京市第一医院)a. 医疗设备处;b. 肿瘤内科,江苏 南京,210006;2. 南京医科大学 医学物理研究中心,江苏 南京,210029

[摘 要] 目的 通过实验分析水模体的CT图像数据,探讨层面位置对西门子CT模拟机图像噪声的影响。方法 调用西门子Sensation Open CT模拟机的日常质量检测程序,使用水模体进行10次测量,得到两种情形(标称管电压120 kV和140 kV)相邻的六层(S1~S6)CT图像。通过手工勾画得到感兴趣区CT值的标准偏差值(图像噪声)数据。对不同层面位置的图像噪声数据进行统计分析和比较研究。结果 120 kV情形下,S1~S6层的图像噪声分别为(10.73±0.05)、(10.66±0.05)、(10.58±0.09)、(10.58±0.08)、(10.72±0.06)、(10.88±0.06)HU;140 kV情形下,S1~S6层的图像噪声分别为(10.09±0.09)、(10.02±0.04)、(9.90±0.07)、(9.93±0.05)、(10.1±0.09)、(10.17±0.12)HU。统计分析显示,不同层面的图像噪声数据总体上有统计学差异。结论 层面位置对西门子CT模拟机图像噪声有影响,中心层面位置的图像噪声小于两端层面的数值。

[关键词] CT模拟机;层面位置;图像噪声;管电压

引言

随着肿瘤精确放射治疗技术的开展和应用,西门子[1-2]、飞利浦[3-4]、GE[5]和东芝[5]等厂家生产等的大孔径CT模拟机在放疗部门的使用越来越广泛。精确放射治疗包含精确定位、精确计划设计和精确照射三个方面,因此对大孔径CT模拟机的图像质量及稳定性提出了更高的要求。美国医学物理学家学会(AAPM)的TG-66号报告[6-7],论述了放疗虚拟定位对CT模拟机质量保证工作中影像质量检验和定位精度检验的要求与方法。我们在前期工作中,对西门子CT模拟机的CT值和CT值线性的长期稳定性进行了研究[8-9]。CT图像成像质量评价体系,包含CT值、空间分辨率、密度分辨率、图像噪声、图像均匀性等,图像噪声评价是其中一个重要方面。高永华等[10]、张富利等[11]和李广喜[12]对CT设备的图像噪声进行了论述和研究。CT图像噪声,表征均匀物质的图像中某一区域内CT值偏离平均值的程度。图像噪声的大小与光子数量、物体大小、扫描层厚和滤波函数等有关,可用感兴趣区域(Region of Interest,ROI)内均匀物质CT值的标准偏差(Standard Deviation,SD)表示[13]。关于图像噪声上述影响因素的研究多有报道,然而层面位置对图像噪声的影响尚未见报道。本研究针对一台西门子Sensation Open CT模拟机,通过实验分析水模体的CT图像数据,探讨层面位置对西门子CT模拟机图像噪声的影响。

1 材料与方法

1.1 设备和模体

我院于2010年引进安装的一台西门子SOMATOM Sensation Open CT模拟机,机架孔径为82 cm[14]。该CT模拟机随机配备了包含水模体、断层厚度模体、电线模体和准线模体四个部分的质控模体组一套。其中,水模体的外直径为20 cm,玻璃壁厚约0.5 cm。主机工作站SIEMENS Syngo系统配置有日常质量检测程序Daily Quality Check(DQC)[10]

1.2 检测方法

调用DQC程序的质控检测使用旋转扫描方式,有效层厚4.8 mm,每次旋转产生相邻的6层图像(记为S1、S2、S3、S4、S5和S6),每次扫描时间0.5 s,FOV(重建视野)为250 mm,卷积核为S80f[15]。该程序对标称管电压120 kV和140 kV两种情形进行测量,相应的电流时间值分别为125 mAs和100 mAs。

使用水模体,调用DQC程序,在83 min内进行10次测量。为了分析每一层(S1~S6)的图像噪声数据,本研究中参照程序自动勾画的ROI,在CT主机工作站Syngo的Viewing界面中,使用手动勾画ROI的方式,得到每一层图像上半径为4 cm(面积50.01 cm2,像素点数为20974)的ROI区域的SD值,即图像噪声值(图1)。使用手工勾画方式分别得到了120 kV和140 kV两种管电压情形S1~S6层的图像噪声值。

图1 手工勾画ROI区获取CT图像噪声值示例

1.3 统计学分析

使用WPS表格记录图像噪声数据,使用SPSS 18.0软件对S1~S6每一层图像噪声的最大值、最小值、平均值、标准差和极差值进行统计描述;对120 kV和140 kV不同层面的图像噪声数据,使用Kruskal-Wallis检验进行统计分析,α取0.05。

2 结果

管电压120 kV和140 kV两种情形下CT模拟机图像第1到第6层(S1~S6)的图像噪声值数据,分别见图2a和图2b。M1~M10分别代表第1次至第10次测量。

图2 两种情形CT图像噪声值的变化趋势

注:a. 120 kV管电压下图像噪声值;b. 140 kV管电压下图像噪声值。

管电压120 kV和140 kV情形下,S1~S6层图像噪声的统计数据(最大值、最小值、平均值、标准差和极差),分别见表1和表2。

表1 120 kV情形下S1~S6图像噪声值的统计数据(HU)

images/BZ_73_1284_1754_2238_1801.png最大值 10.8 10.7 10.7 10.7 10.8 11.0最小值 10.7 10.6 10.4 10.4 10.6 10.8平均值 10.73 10.66 10.58 10.58 10.72 10.88标准差 0.05 0.05 0.09 0.08 0.06 0.06极差 0.1 0.1 0.3 0.3 0.2 0.2

表2 140 kV情形下S1~S6图像噪声值的统计数据(HU)

参数 S1 S2 S3 S4 S5 S6最大值 10.2 10.1 10.0 10.0 10.2 10.3最小值 9.9 10.0 9.8 9.9 10.0 10.0平均值 10.09 10.02 9.90 9.93 10.10 10.17标准差 0.09 0.04 0.07 0.05 0.09 0.12极差 0.3 0.1 0.2 0.1 0.2 0.3

使用Kruskal-Wallis检验分析两种情形下S1~S6层图像噪声数据,结果显示:120 kV情形下,S1~S6层总体上有统计学差异(P<0.001);140 kV情形下,S1~S6层总体上有统计学差异(P<0.001)。S1和S6是两端的层面,S3和S4是中间层面。对S1和S3,S1和S4,S6和S3,S6和S4进行两两比较。120 kV情形下,P值分别为0.026、0.017、<0.001和<0.001;140 kV情形下,P值分别为0.003、0.014、<0.001和<0.001。上述P值均小于0.05,说明组间比较均有统计学差异。

以层面编号(代表层面位置)为横坐标,以图像噪声的平均值为纵坐标,得到两种情形图像噪声的平均值随层面位置的变化趋势,结果示于图3。

图3 图像噪声的平均值随层面位置的变化

为了对上述数据结果进行验证,我们使用本院体检中心的一台Emotion 16 CT模拟机进行了补充实验。该设备的DQC程序设置与Sensation Open CT模拟机略有不同。使用该设备的DQC程序,进行了3次数据测量。130 kV情形下,S1、S2和S3层面CT图像噪声的平均值分别为10.27、9.97和10.07 HU;110 kV情形下,S1、S2和S3层面CT图像噪声的平均值分别为14.73、14.37和14.43 HU。S1和S3代表两端层面,S2代表中间层面。

3 讨论

由图2、表1和表2可见,该CT模拟机的图像噪声值在10次测量中的变化较小。图像噪声值的变化范围,120 kV情形为10.4~11.0 HU,140 kV情形为9.8~10.3 HU。表1和表2中12个极差值均≤0.3 HU。由此可见,该CT模拟机图像噪声值具有较好的短期稳定性。

由表1和表2,Kruskal-Wallis检验分析及图3可见,层面位置是影响CT设备图像噪声的因素之一。两种情形中,S3和S4的图像噪声值整体上明显低于两端(S1和S6)的值,中间和两端图像噪声平均值的差异约为0.2~0.3 HU。使用Emotion 16 CT机的补充实验数据,也初步验证了中心层面位置的噪声值低于两端层面噪声值的结果。

西门子CT模拟机在进出机架的方向上,排列着许多独有的超高速稀土陶瓷探测器。在一次旋转出束中,治疗床不动,产生了相邻的6层CT图像。S3和S4对应了中间位置的探测器,S1和S6对应了两端位置的探测器。扫描层面与扇形束X射线相对位置的不同,可能导致到达不同扫描层面的光子数量有差异,进而导致CT图像噪声有差异。中间扫描层面,理论上光子的数量最大,因而导致其图像噪声值小于两端层面的值。但是在实际工作中,影响因素很多,例如矫正曲线、扫描条件和阳极靶面角等,这些因素有的还会受环境温湿度、外部电压、球管温度等影响。因此,层面位置对该CT模拟机图像噪声影响的相关结果,其适用性还需要对其他具体的CT设备实际测量去验证。我们在前期工作中,发现CT图像均匀性也与层面位置相关[14]

影响CT图像噪声的因素有光子的数量(涉及管电压、mAs),物体的大小,扫描的层厚,滤波函数,矩阵的大小,散射线和电子噪声(探测器噪声)等。本研究的结果显示了层面位置也是影响CT设备图像噪声的因素之一。由表1和表2的数据比较可见,管电压120 kV和140 kV时的噪声水平整体上有区别,140 kV的图像噪声整体上低于120 kV的图像噪声。管电压对图像噪声的影响,许多文献都有提及[16-17]

郑庆增等[17]全面研究了CT图像噪声的各种影响因素,如管电压、管电流、层厚、不同算法重建以及扫描位置等。其中关于扫描位置方面,该研究采用管电压120 kV,管电流时间值300 mAs,层厚3 mm,重建算法为B20f,扫描位置位于孔径中心及偏上7 cm、下12 cm、左20 cm、右20 cm时,扫描头部电子密度模体。该研究的扫描位置时在同一层面的上下左右位置,本文研究的是进出机架方向的不同层面位置。

本研究也有不足之处。本研究仅对一台西门子CT模拟机的图像噪声进行了详细研究,而Emotion 16 CT机的数据也只是初步验证了层面位置对CT图像噪声影响的结论,层面位置对CT图像噪声影响的整体规律尚需进一步的实验进行深入研究。此外,本研究仅采用了DQC程序的扫描重建参数(电流时间、层厚、卷积核、FOV等),并未使用放疗工作常用的CT扫描重建参数。将来我们将针对上述不足,对更多CT设备,以及采用放疗工作常用的螺旋CT扫描重建参数,深入探究层面位置对放疗CT图像噪声的影响。

综上所述,本研究的结果表明,层面位置对本部门西门子CT模拟机图像噪声有影响,该西门子CT模拟机中心层面位置的图像噪声小于两端层面的值。

致谢

衷心感谢西门子公司郭强工程师关于层面位置对CT图像噪声影响原因分析的讨论。

[参考文献]

[1] 宋孟孟,王运来,黄祥,等.肺部肿瘤基于最大密度投影图像和四维CT勾画的内靶区差异性分析[J].中国医学物理学杂志,2018,35(2):135-140.

[2] 张国前,张书旭,余辉,等.放疗定位扫描条件对CT值的影响[J].中国医疗设备,2015,30(3):13-16.

[3] 廖雄飞,黎杰,王培.CT模拟定位机的扫描参数对放疗计划系统剂量计算的影响[J].肿瘤预防与治疗,2015,28(1):49-51.

[4] 欧阳斌,丘敏敏,钟嘉健,等.CT扫描及重建参数对放疗图像质量影响的研究[J].中国医药导报,2017,14(1):4-7.

[5] Tomic T,Papaconstadopoulosa P,Aldelaijana S,et al.Image quality for radiotherapy CT simulators with different scanner bore size[J].Phys Med,2018,45:65-71.

[6] Mutic S,Palta JR,Butker EK,et al.Quality assurance for computed-tomography simulators and the computedtomography-simulation process: report of the AAPM Radiation Therapy Committee Task Group No. 66[J].Med Phys,2003,30(10):2762-2792.

[7] 邓小武,黄邵敏,祁振宇.CT模拟机的质量控制和质量保证检验[J].中国肿瘤,2004,13(9):546-550.

[8] 时飞跃,任军,吴正参,等.西门子Sensation Open CT模拟机CT值稳定性分析[J].中华放射肿瘤学杂志,2017,26(12):1407-1410.

[9] 任军,时飞跃,柏正璐.西门子CT模拟机CT值线性的稳定性分析[J].中国医疗设备,2018,33(6):33-35.

[10] 高华永,晁勇,刘帅,等.CT设备的质量控制分析与探讨[J].中国医疗设备,2014,29(6):8-10.

[11] 张富利,陈建平,郑明民.PHILIPS大孔径CT模拟机噪声测试研究[J].中国医疗设备,2009,24(1):123-124.

[12] 李广喜.CT设备质量检测存在的问题与质控措施分析[J].中国卫生标准管理,2017,8(14):124-126.

[13] JJG 961-2017,医用诊断螺旋计算机断层摄影装置(CT)X射线辐射源[S].

[14] 时飞跃,柏正璐,秦航,等.西门子CT模拟机图像均匀性的长期稳定性分析[J].中国医学物理学杂志,2018,35(9):1035-1038.

[15] 时飞跃,柏正璐,秦伟.西门子CT模拟机水的CT值、图像噪声和管电压的长期稳定性分析[J].中国医学物理学杂志,2017,34(12):1232-1235.

[16] 亓恒涛,秦伟昌,王巍,等.64层螺旋CT噪声测试及其影响因素分析[J].中华放射医学与防护杂志,2007,27(2):195-198.

[17] 郑庆增,鞠忠建,邵莹,等.模拟定位CT的机型及扫描参数对CT值和噪声的影响[J].中国医学物理学杂志,2017,34(5):439-444.

Influence of Slice Position on Image Noise for a SIEMENS CT Simulator

JIANG Zhiwei1a, HUO Xinying1b, SHI Feiyue1b,2, QIN Wei1b, LIN Qiang1a, QIN Hang1a
1. a. Department of Medical Equipment, b. Department of Oncology, Nanjing First Hospital, Nanjing Medical University, Nanjing Jiangsu 210006, China; 2. Medical Physics Center, Nanjing Medical University, Nanjing Jiangsu 210029, China

Abstract: Objective To study influence of the slice position on image noise for SIEMENS CT simulator by analyzing the experimental CT image data using the water phantom. Methods The daily quality check program for the SIEMENS Sensation Open CT simulator was utilized for 10 measurements with a water phantom, and the CT images of 6 neighboring slices (S1-S6) were obtained under two conditions of the nominal tube voltages of 120 kV and 140 kV, respectively. Then the values of image noise were obtained by manual delineation with the region of interest. Statistical analysis and comparative study were conducted on the image noise data for different slice positions. Results Under the case with the nominal tube voltages of 120 kV, the values of image noise for the first to sixth slices were (10.73±0.05), (10.66±0.05), (10.58±0.09), (10.58±0.08), (10.72±0.06), and (10.88±0.06) HU,respectively. Under the case with the nominal tube voltages of 140 kV, the corresponding values were (10.09±0.09), (10.02±0.04),(9.90±0.07), (9.93±0.05), (10.1±0.09), (10.17±0.12) HU, respectively. Statistical analysis showed that overall there was statistical difference in image noise at different slice positions. Conclusion The slice position has influence on the image noise of the SIEMENS CT simulator. The values of image noise in the middle slice position are smaller than those in the two ends.

Key words: CT simulator; slice position; image noise; tube voltage

收稿日期:2019-08-20

通信作者:时飞跃,物理学博士,主要研究方向为肿瘤放射物理。

通信作者邮箱:shifeiyue2013@126.com

[中图分类号] R814.42

[文献标识码] A

doi:10.3969/j.issn.1674-1633.2020.03.016

[文章编号] 1674-1633(2020)03-0064-03