高级搜索

一种基于毛细管X光透镜的微型锥束CT扫描仪

,  ,  ,  ,  ,  , 

downloadPDF

上一篇

下一篇

周腊珍, 夏文静, 许倩倩, 陈赞, 李坊佐, 刘志国, 孙天希. 一种基于毛细管X光透镜的微型锥束CT扫描仪[J]. 物理学报, 2022, 71(9): 090701-1. doi: 10.7498/aps.71.20212195
引用本文: 周腊珍, 夏文静, 许倩倩, 陈赞, 李坊佐, 刘志国, 孙天希. 一种基于毛细管X光透镜的微型锥束CT扫描仪[J]. 物理学报, 2022, 71(9): 090701-1. doi: 10.7498/aps.71.20212195
shu
La-Zhen Zhou, Wen-Jing Xia, Qian-Qian Xu, Zan Chen, Fang-Zuo Li, Zhi-Guo Liu, Tian-Xi Sun. Micro cone-beam CT scanner based on X-ray polycapillary optics[J]. Acta Physica Sinica, 2022, 71(9): 090701-1. doi: 10.7498/aps.71.20212195
Citation: La-Zhen Zhou, Wen-Jing Xia, Qian-Qian Xu, Zan Chen, Fang-Zuo Li, Zhi-Guo Liu, Tian-Xi Sun. Micro cone-beam CT scanner based on X-ray polycapillary optics[J]. Acta Physica Sinica, 2022, 71(9): 090701-1. doi: 10.7498/aps.71.20212195
shu

一种基于毛细管X光透镜的微型锥束CT扫描仪

  • 1.

    赣南医学院医学信息工程学院, 赣州　341000

  • 2.

    赣南医学院, 组织工程材料与生物制造江西省重点实验室, 赣州　341000

  • 3.

    北京师范大学核科学与技术学院, 射线束技术教育部重点实验室, 北京　100875

    • 通讯作者: E-mail: lfz880920@163.com; 

  • 中图分类号: 07.85.Fv, 41.50.+h, 87.57.C-, 87.59.-e

Micro cone-beam CT scanner based on X-ray polycapillary optics

  • 1.

    School of Medical Information Engineering, Gannan Medical University, Ganzhou 341000, China

  • 2.

    Key Laboratory of Biomaterials and Biofabrication in Tissue Engineering of Jiangxi Province, Gannan Medical University, Ganzhou 341000, China

  • 3.

    Key Laboratory of Beam Technology Ministry of Education, College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China

    • Corresponding author: E-mail: lfz880920@163.com; 

  • MSC: 07.85.Fv, 41.50.+h, 87.57.C-, 87.59.-e

  • 摘要: 活体小动物成像系统是疾病研究、新药开发的重要组成部分. 其中, X射线微型锥束计算机断层成像(X-ray micro cone-beam computed tomography, Micro-CBCT)能实现数十至数百微米空间分辨率的解剖结构成像研究. Micro-CBCT成像仪的一个关键挑战是其空间和对比度分辨率主要取决于X射线源焦斑大小、探测器分辨率和系统几何结构等因素. 为提高Micro-CBCT的空间分辨率、对比度分辨率和成像均一性, 本文基于毛细管X光透镜研制了一款能够调控照射X射线束斑孔径的Micro-CBCT扫描仪, 用于小动物成像研究. 此系统由微焦斑X射线源、非晶硅平板探测器、旋转工作台和控制电脑组成, 并采用Feldkamp-Daivs-Kress算法重建投影图像. 对该系统的性能进行了测试, 结果表明, 系统在10%调制传递函数下的空间分辨率为9.1 lp/mm, 提高了1.35倍. 同时, 由于毛细管X光透镜对低能X射线的吸收和散射抑制作用, 实现了2倍以上的对比度增强, 减轻了多色X射线束硬化效应引起的图像均一性恶化问题. 应用该Micro-CBCT系统对麻醉小鼠进行了活体成像, 验证了该系统在小动物成像研究中的实用性.
    Abstract: In-vivo small animal imaging system is an important part of disease research and new drug development. It is essential for living small animal imaging system to be able to provide the anatomical structure, molecular and functional information. The X-ray micro cone-beam computed tomography (micro-CBCT) can perform longitudinal study with a resolution of tens-to-hundreds of microns in a short imaging time at a relatively low cost. Furthermore, it is easy to combine with other modalities to provide abundant information about small animals. A key challenge to the micro-CBCT scanner is that its spatial and contrast resolution determined primarily by the X-ray focal spot size, the detector element size, and the system geometry. Aiming to improve the spatial resolution, contrast resolution, and imaging uniformity of the micro-CBCT system, we use the X-ray polycapillary optics for adjusting the X-ray source. A micro-CBCT based on X-ray polycapillary optics with a large field of view is constructed for the small animal imaging study. The micro-CBCT system is composed of microfocus X-ray tube with an attached polycapillary focusing X-ray lens, amorphous silicon-based flat panel detector, rotation stage, and controlling PC. The Feldkamp-Daivs-Kress (FDK) algorithm is adopted to reconstruct the image. The system performances are evaluated. The magnification of this micro-CBCT system is 1.97. The results show that the spatial resolution of the system at 10% modulation transfer function (MTF) is 9.1 lp/mm, which is 1.35 times higher than that in the case of no optics. The image uniformity deterioration caused by hardening effect is effectively alleviated by filtrating the low energy X-rays with the X-ray polycapillary optics and the contrast enhancement is more than twice. The anesthetic rats are imaged with this micro-CBCT system in vivo and the practicability of the system in small animal imaging research is verified.
  • 加载中

  • 图 1  X射线源焦斑大小、SOD和SDD共同决定了Micro-CBCT的有效探测器孔径大小　(a) 焦点大小与有效探测器孔径α成正比; (b) SOD/SDD比值与有效探测器孔径α成正比

    Figure 1.  X-ray tube focal spot size, SOD and SDD jointly determine the effective detector aperture size of the micro-CT system: (a) Focal spot size is proportional to the effective detector aperture (α); (b) ratio of SOD/SDD is proportional to the effective detector aperture (α).

    下载: 全尺寸图片 幻灯片

    图 2  X射线焦斑大小的半影效应

    Figure 2.  Penumbra effect of X-ray focal spot size.

    下载: 全尺寸图片 幻灯片

    图 3  Micro-CBCT系统, 该系统由一个结合PFXRL的微聚焦X射线源、一个旋转样品台和一个非晶硅平板探测器组成　(a) Micro-CBCT原理图; (b) Micro-CBCT实物图; (c) 采用的PFXRL实物图

    Figure 3.  Micro-CBCT system. The system consists of a microfocus X-ray source combined with a PFXRL, a rotating sample stage and an amorphous silicon-based FPD: (a) Micro-CBCT schematic diagram; (b) desktop micro-CBCT system; (c) the PFXRL.

    下载: 全尺寸图片 幻灯片

    图 4  PFXRL的出口焦斑尺寸和收敛度随X射线能量的依赖关系

    Figure 4.  Energy dependence of the output focal spot size and convergence for the PFXRL, respectively.

    下载: 全尺寸图片 幻灯片

    图 5  PFXRL的传输效率和能量密度增益随X射线能量的依赖关系

    Figure 5.  Energy dependence of transmission efficiency and gain in power density of the PFXRL, respectively.

    下载: 全尺寸图片 幻灯片

    图 6  使用和不使用PFXRL两种条件下的Micro-CBCT图像和对应的MTF　(a)使用PFXRL; (b)不使用PFXRL; (c) MTF曲线

    Figure 6.  Measured images and corresponding MTF of the micro-CBCT system with and without PFXRL: (a) With PFXRL; (b) without PFXRL; (c) MTF curves.

    下载: 全尺寸图片 幻灯片

    图 7  对于测量的水模, 在使用和不使用PFXRL下的Micro-CBCT系统对比度分辨率与管电压的关系

    Figure 7.  Measured contrast resolution of Micro-CBCT system as a function of tube voltages for water phantom.

    下载: 全尺寸图片 幻灯片

    图 8  水模的均匀性响应　(a) 不使用PFXRL, 重建的水模中平横断面图像及绿线对应的CT值; (b) 使用PFXRL, 重建的水模中平横断面图像及绿线对应的CT值; (c) 采用0.5 mm厚的铝片附加滤过, 重建的水模中平横断面图像及绿线对应的CT值. 空气和水的CT值分别归一化为0和50

    Figure 8.  Uniformity response of the water phantom: (a) Reconstructed transaxial image of the uniformity phantom without using PFXRL and radial signal profile taken from the green line; (b) reconstructed transaxial image of the uniformity phantom with using PFXRL and radial signal profile taken from the green line; (c) reconstructed transaxial image of the uniformity phantom with a 0.5 mm thick aluminum sheet as filter and radial signal profile taken from the green line. The CT values of air and water are normalized to 0 and 50, respectively.

    下载: 全尺寸图片 幻灯片

    图 9  水模的噪声响应　(a)不使用PFXRL (附加滤过: 0.5 mm Al); (b)使用PFXRL

    Figure 9.  Noise response of the water phantom: (a) Without PFXRL (additional filtration: 0.5 mm Al); (b) with PFXRL.

    下载: 全尺寸图片 幻灯片

    图 10  小鼠Micro-CBCT成像　(a)使用PFXRL; (b)不使用PFXRL

    Figure 10.  CT images of rats: (a) With PFXRL; (b) without PFXRL.

    下载: 全尺寸图片 幻灯片

    图 11  麻醉小鼠横断面图像比较　(a)使用PFXRL; (b)不使用PFXRL

    Figure 11.  Comparison of axial images of the anesthetized mice: (a) With PFXRL; (b) without PFXRL.

    下载: 全尺寸图片 幻灯片

    图 12  麻醉小鼠肺、肾和下脊柱区域的基于PFXRL的Micro-CBCT图像. 每个图像的窗宽窗位设置不同以呈现出感兴趣的结构. 横断面((a), (d), (g))、冠状面((b), (e), (h))和矢状面((c), (f), (i))切片展示了各向同性的空间分辨率和好的软组织对比度. 肺内支气管结构、肾脏与周围的肌肉和脂肪、椎骨和椎间隙都清晰可见. 图像中的垂直比例尺显示1 cm的间距

    Figure 12.  PFXRL-based Micro-CBCT images of the lung, kidney and lower spine of the anesthetized mice. The window and level settings are varied in each image to allow visualization of the structures of interest. Axial ((a), (d), (g)) , coronal ((b), (e), (h)) , and sagittal ((c), (f), (i)) slices qualitatively demonstrate isotropic spatial resolution, with excellent soft-tissue contrast in each case. Bronchial structure within the lungs is clearly identifiable, the kidney is well delineated from surrounding muscle and fat, and fine detail in the vertebrae and intervertebral spaces is demonstrated. The vertical scale in the images shows 1 cm spacing.

    下载: 全尺寸图片 幻灯片

    表 1  PFXRL的基本参数

    Table 1.  Parameters of the PFXRL.

    PFXRL
    Length /mm69.4
    Input focal distance /mm76.3
    Output focal distance /mm20.1
    Diameter of IFS at 17.4 keV /${\mu }\mathrm{m}$169.2
    Diameter of OFS at 17.4 keV /${\mu }\mathrm{m}$46.7
    Channel inner diameter of capillary at input/output /μm10.4
    下载: 导出CSV
  • [1] Gregory S G, Sekhon M, Schein J, et al. 2002 Nature 418 743 doi: 10.1038/nature00957
    [2] Ntziachristos V, Ripoll J, Wang L V, Weissleder R 2005 Nat. Biotechnol. 23 313 doi: 10.1038/nbt1074
    [3] Guerra A D, Belcari N 2007 Nucl. Instrum. Meth. Phys. Res. A 583 119 doi: 10.1016/j.nima.2007.08.187
    [4] Badea C T, Drangova M, Holdsworth D W, Johnson G A 2008 Phys. Med. Biol. 53 319 doi: 10.1088/0031-9155/53/2/001
    [5] Jan M L, Ni Y C, Chen K W, Ching H 2006 Nucl. Instrum. Meth. Phys. Res. A 569 314 doi: 10.1016/j.nima.2006.08.106
    [6] Biederer J, Mirsadraee S, Beer M, Molinari F, Puderbach M 2012 Insights Into Imaging 3 373 doi: 10.1007/s13244-011-0142-z
    [7] Hoyer C, Gass N, Fahr W W, Sartorius A 2014 Neuropsychobiology 69 187 doi: 10.1159/000360859
    [8] Kunjachan S, Ehling J, Storm G, Kiessling F, Lammers T 2015 Chem. Rev. 115 10907 doi: 10.1021/cr500314d
    [9] Eghtedari M, Oraevsky A, Copland J A, Kotov N A, Conjusteau A, Motamedi M 2007 Nano Lett. 7 1914 doi: 10.1021/nl070557d
    [10] Taruttis A, Ntziachristos V 2015 Nat. Photonics 9 219 doi: 10.1038/nphoton.2015.29
    [11] Paulus M J, Gleason S S, Kennel S J, Hunsicker P R, Johnson D K 2000 Neoplasia 2 62 doi: 10.1038/sj.neo.7900069
    [12] 罗召洋, 杨孝全, 孟远征, 邓勇 2010 物理学报 58 8237 doi: 10.7498/aps.59.8237 Luo Z Y, Yang X Q, Meng Y Z, Deng Y 2010 Acta Phys. Sin. 58 8237 doi: 10.7498/aps.59.8237
    [13] 魏星, 闫镔, 张峰, 李永丽, 席晓琦, 李磊 2014 物理学报 63 058702 doi: 10.7498/aps.63.058702 Wei X, Yan B, Zhang F, Li Y L, Xi X Q, Li L 2014 Acta Phys. Sin. 63 058702 doi: 10.7498/aps.63.058702
    [14] Mazel V, Reiche I, Busignies V, Walter P, Tchoreloff P 2011 Talanta 85 556 doi: 10.1016/j.talanta.2011.04.027
    [15] Sun T, Liu Z, Li Y, Lin X, Wang G, Zhu G, Xu Q, Luo P, Pan Q, Liu H 2010 Nucl. Instrum. Meth. Phys. Res. A 622 295 doi: 10.1016/j.nima.2010.06.098
    [16] Macdonald C A, Gibson W M 2003 X-Ray Spectrom. 32 258 doi: 10.1002/xrs.648
    [17] Albertini V R, Paci B, Generosi A, Dabagov S B, Kumakhov M A 2007 Spectrochim. Acta B 62 1203 doi: 10.1016/j.sab.2007.07.012
    [18] Huang R, Bilderback D H 2006 J. Synchrotron Radiat. 13 74 doi: 10.1107/S0909049505038562
    [19] Balaic D X, Barnea Z, Nugent K A, Garrett R F, Wilkins S W 1996 J. Synchrotron Radiat. 3 289 doi: 10.1107/S0909049596009351
    [20] MacDonald C A, Owens S M, Gibson W M 1999 J. Appl. Crystallogr. 32 160 doi: 10.1107/S0021889898008826
    [21] Bjeoumikhov A, Bjeoumikhova S, Langhoff N, Wedell R 2005 Appl. Phys. Lett. 86 144102 doi: 10.1063/1.1897441
    [22] Sun T, Liu Z, Ding X 2007 Nucl. Instrum. Meth. Phys. Res. B 262 153 doi: 10.1016/j.nimb.2007.05.012
    [23] Sun T, Peng S, Liu Z, Sun W, Ma Y, Ding X 2013 J. Appl. Crystallogr. 46 1880 doi: 10.1107/S0021889813027088
    [24] Sun T, Macdonald C A 2013 J. Appl. Phys. 113 053104 doi: 10.1063/1.4789799
    [25] Lamb J S, Bilderback D H, Pollack L, Kwok L, Smilgies D M 2007 J. Appl. Crystallogr. 40 193 doi: 10.1107/S0021889806044505
    [26] Barrea R A, Huang R, Cornaby S, Bilderback D H, Irving T C 2009 J. Synchrotron Radiat. 16 76 doi: 10.1107/S0909049508039782
    [27] Zeng X, Duewer F, Feser M, Huang C, Lyon A, Tkachuk A, Yun W 2008 Appl. Opt. 47 2376 doi: 10.1364/AO.47.002376
    [28] Li F, Liu Z, Sun T, Jiang B, Zhu Y 2016 J. Chem. Phys. 144 104201 doi: 10.1063/1.4943268
    [29] Li F, Liu Z, Sun T 2016 J. Appl. Crystallogr. 49 627 doi: 10.1107/S1600576716000376
    [30] Li F, Liu Z, Sun T 2016 Rev. Sci. Instrum. 87 093106 doi: 10.1063/1.4962709
    [31] Li F, Liu Z, Sun T 2016 Food Chem. 210 435 doi: 10.1016/j.foodchem.2016.05.012
    [32] Li F, Liu Z, Sun T, Ma Y, Ding X 2015 Food Control 54 120 doi: 10.1016/j.foodcont.2015.01.043
    [33] Abreu C C, Kruger D G, MacDonald C A, Mistretta C A, Peppler W W, Xiao Q F 1995 Med. Phys. 22 1793 doi: 10.1118/1.597635
    [34] Goertzen A L, Nagarkar V, Street R A, Paulus M J, Boone J M, Cherry S R 2004 Phys. Med. Biol. 49 5251 doi: 10.1088/0031-9155/49/23/004
    [35] Kim H K, Min K C, Achterkirchen T, Lee W 2009 IEEE Trans. Nucl. Sci. 56 1179 doi: 10.1109/TNS.2009.2014232
    [36] Feldkamp L A, Davis L C, Kress J W 1984 J. Opt. Soc. Am. A 1 612 doi: 10.1364/JOSAA.1.000612
    [37] Flannery B P, Deckman H W, Roberge W G, D'Amico K L 1987 Science 237 1439 doi: 10.1126/science.237.4821.1439
    [38] Sun T, Ding X 2005 J. Appl. Phys. 97 124904 doi: 10.1063/1.1935134
    [39] Kai Y, Kwan A L C, Miller D W F, Boone J M 2006 Med. Phys. 33 1695 doi: 10.1118/1.2198187
    [40] Kwan A L C, Boone J M, Yang K, Huang S Y 2007 Med. Phys. 34 275
    [41] 余晓锷, 占杰, 李萍, 李婵娟 2006 第四军医大学学报 27 978 doi: 10.3321/j.issn:1000-2790.2006.11.006 Yu X E, Zhan J, Li P, Li C J 2006 J. Fourth Mil. Med. Univ. 27 978 doi: 10.3321/j.issn:1000-2790.2006.11.006
  • 加载中
图( 12) 表( 1)
计量
  • 文章访问数:  381
  • HTML全文浏览数:  381
  • PDF下载数:  1
  • 施引文献:  0
出版历程
  • 收稿日期:  2021-11-28
  • 刊出日期:  2022-05-05

一种基于毛细管X光透镜的微型锥束CT扫描仪

    通讯作者: E-mail: lfz880920@163.com; 
  • 1. 赣南医学院医学信息工程学院, 赣州　341000
  • 2. 赣南医学院, 组织工程材料与生物制造江西省重点实验室, 赣州　341000
  • 3. 北京师范大学核科学与技术学院, 射线束技术教育部重点实验室, 北京　100875
  • 收稿日期:  2021-11-28

摘要: 活体小动物成像系统是疾病研究、新药开发的重要组成部分. 其中, X射线微型锥束计算机断层成像(X-ray micro cone-beam computed tomography, Micro-CBCT)能实现数十至数百微米空间分辨率的解剖结构成像研究. Micro-CBCT成像仪的一个关键挑战是其空间和对比度分辨率主要取决于X射线源焦斑大小、探测器分辨率和系统几何结构等因素. 为提高Micro-CBCT的空间分辨率、对比度分辨率和成像均一性, 本文基于毛细管X光透镜研制了一款能够调控照射X射线束斑孔径的Micro-CBCT扫描仪, 用于小动物成像研究. 此系统由微焦斑X射线源、非晶硅平板探测器、旋转工作台和控制电脑组成, 并采用Feldkamp-Daivs-Kress算法重建投影图像. 对该系统的性能进行了测试, 结果表明, 系统在10%调制传递函数下的空间分辨率为9.1 lp/mm, 提高了1.35倍. 同时, 由于毛细管X光透镜对低能X射线的吸收和散射抑制作用, 实现了2倍以上的对比度增强, 减轻了多色X射线束硬化效应引起的图像均一性恶化问题. 应用该Micro-CBCT系统对麻醉小鼠进行了活体成像, 验证了该系统在小动物成像研究中的实用性.

English Abstract

    全文HTML

    1.   引　言

    • 基因组学研究表明, 实验室中常用的小动物, 如小鼠, 与人类基因具有高度的遗传同源性[1]. 这为在小动物身上研究与人类相关的生物医学问题, 如疾病机制的研究和新药的开发奠定了理论基础. 小动物研究在现代分子生物学实验中越来越重要. 专用的高分辨率小动物成像系统已成为生物医学研究的重要工具, 主要包括微型X射线计算机断层成像(X-ray micro-computed tomography, Micro-CT)、微型磁共振成像(micro-magnetic resonance imaging, Micro-MRI)、微型正电子发射断层成像(micro-positron emission tomography, Micro-PET)和微型单光子发射计算机断层成像(micro-single photon emission computed tomography, Micro-SPECT). 这些成像系统允许研究人员非侵入性地筛选出用于突变或病理研究的动物, 并在体外和/或体内监测疾病进展和对治疗的反应情况[2-10]. Micro-CT是其中一种专用于小动物结构成像的方法. Micro-CT系统可提供高分辨率图像、快速数据采集、对骨组织和软组织的良好敏感性, 尤其是在使用造影剂的情况下. Micro-CT成像空间分辨率主要由X射线焦斑大小、探测器分辨率和系统几何结构决定, 而对比度分辨率主要由X射线通量和探测器分辨率决定[11]. 同时, Micro-CT图像往往会出现模糊伪影, 一般是由于物理效应(如X光束硬化、康普顿散射等)和仪器效应造成的[12,13].

      如何进一步实现Micro-CT的高空间分辨率和对比度分辨率以及消除或降低模糊效应是目前需要解决的问题. X光透镜为解决这些问题提供了一种可能的方法. X光透镜聚焦X射线束可提高其强度增益和调控焦斑大小, 已广泛应用于基于实验室或同步辐射光源的原位分析技术领域[14-32]. 特别是, X光透镜在射线照相仪中可提供光束整形和散射抑制作用, 这对小动物成像是特别有益的[33].

      目前, 用于Micro-CBCT的平板探测器主要有3种类型: 电荷耦合器件(charge coupled device, CCD)、互补型金属氧化物半导体(complementary metal oxide semiconductor, CMOS)和非晶硅平板探测器[34]. CCD平板探测器的像素大小从几微米到几百微米, CMOS平板探测器像素大小通常为50—200 μm, 而非晶硅平板探测器的像素大小通常为100—400 μm. 纯粹从像素大小的角度来看, CMOS和CCD平板探测器更适合高分辨率成像, 而非晶硅平板探测器需要在高放大系数下才能实现高分辨率成像. CCD平板探测器会带来几何畸变和横向畸变问题, 难以实现大面积高质量成像. CMOS平板探测器具有高填充系数, 几乎没有图像延迟, 但由于其对X射线的敏感性而引起的老化问题将大大降低图像质量, 比如噪声水平、分辨率和均匀性[34,35]. 非晶态硅平板探测器对X射线敏感性低, 帧率高, 更适合快速成像, 这对小动物纵向研究具有重要意义.

      本文设计了一种基于X光透镜的Micro-CBCT扫描仪. Micro-CBCT系统由带有X光透镜的微焦斑X射线管、非晶硅平板探测器、旋转工作台和控制计算机组成. 电动旋转台在程序的控制下, 带动小动物实现360°一步一步旋转. 采集投影数据后, 采用锥束Feldkamp-Daivs-Kress (FDK)算法重建三维图像[36]. 对Micro-CBCT系统的放大系数、空间分辨率、对比度分辨率和成像均匀性进行了评价. 最后, 以麻醉小鼠为样品进行活体成像, 验证基于X光透镜的低剂量Micro-CBCT系统的有效性和实用性.

    2.   实验装置和方法

      2.1.   理论基础

    • 影响CT扫描仪空间分辨率的因素主要有两类: 一类是几何因素, 另一类是算法因素[11]. 其中几何因素方面, 有效探测器孔径大小是制约空间分辨率的基本因素. 有效探测器孔径越小, 空间分辨率越高. 有效探测器孔径大小与X射线源焦斑大小、探测器像素单元大小、X射线源焦斑到物体的距离(distance from X-ray focal spot to object, SOD)和X射线源焦斑到平板探测器的距离(distance from X-ray focal spot to detector, SDD)有关. 这些参数通过影响有效探测器孔径大小来影响空间分辨率的大小, 如图1所示. 其中, X射线源焦斑大小与有效探测器孔径大小成正比(如图1(a)所示), X射线源焦斑尺寸越小, 空间分辨率越高. 此外, SOD/SDD越小(即物体离X射线源焦点越近), 空间分辨率越好(如图1(b)所示).

      此外, X射线源焦斑的大小直接影响成像的清晰度. X射线源焦斑越大, 半影面积增大, 图像清晰度将降低(如图2所示). 减小X射线源焦斑尺寸可以提高图像的清晰度, 但较小的焦斑尺寸会导致X射线强度的降低. 要达到相同的照射剂量, 就必须增加照射时间或增加X射线强度增益.

      总体上, 系统分辨率由X射线源焦斑大小、探测器分辨率、SOD和SDD共同决定. 空间分辨率因探测器像素单元大小而产生的展宽可表示为标准偏差为$ {\sigma }_{\mathrm{D}\mathrm{e}\mathrm{t}\mathrm{e}\mathrm{c}\mathrm{t}\mathrm{o}\mathrm{r}} $的高斯函数[11]:

      式中, $ {X}_{\mathrm{D}\mathrm{e}\mathrm{t}\mathrm{e}\mathrm{c}\mathrm{t}\mathrm{o}\mathrm{r}} $为探测器像素单元大小. 同样, 由X射线源焦斑大小导致的空间分辨率展宽可表示为标准偏差为$ {\sigma }_{\mathrm{S}\mathrm{o}\mathrm{u}\mathrm{r}\mathrm{c}\mathrm{e}} $的高斯函数[11]:

      式中, $ {X}_{\mathrm{f}, \mathrm{F}\mathrm{W}\mathrm{H}\mathrm{M}} $是X射线源焦斑的半高宽.

      由探测器像素单元大小和X射线源焦斑大小引起的空间分辨率展宽分别是独立的函数, 而这两个独立函数引起的系统空间分辨率展宽可正交相加得到[11]

      因此, 设计Micro-CBCT扫描仪需要选择合适的X射线源微焦斑大小和合适的探测器像素单元大小, 同时优化SOD和SDD, 以实现理想的空间分辨率、视场大小和数据采集时间分辨率.

      微焦斑X射线源出射的X射线通量有限, 且探测器像素单元尺寸较小, Micro-CBCT光子计数通常满足泊松统计, 对比度分辨率可以表示为[11]

      式中, ${\sigma }_{\text{CT Number}}$是测量的CT值的标准偏差, D是样品的直径, $ \mu $是与能量相关的X射线吸收衰减系数, N是视场内总的X光子计数, $ \Delta x $是探测器像素大小. 对比度分辨率在X射线能量满足μ = 2/D情况下达到最优值[11,37].

      多色X射线能谱导致一个重要的问题, 即光束硬化伪影. X射线吸收衰减系数与X射线能量大小具有很强的相关性, 尤其是小动物成像研究所用的低能量X射线. 当一束X射线透射过样品时, 低能量的X射线在样品表面附近被大量吸收, 导致测量到的CT值在样品边缘附近较高. 通过对X射线束进行预过滤, 可以增加X射线束的单色性, 从而降低硬化伪影. 然而, X射线束的硬化效应很难完全消除.

      多毛细管聚焦X光透镜(polycapillary focusing X-ray lens, PFXRL)可将位于PFXRL入口焦点上的微焦斑X射线源发散的X射线束聚焦到出口焦点上, 获得更小的聚焦尺寸和高的强度增益[16,22-32], 从而实现更小的有效探测器孔径大小, 减少半影模糊. 此外, 由于PFXRL对低能量X射线的吸收, 可以起到低能滤波的作用, 以减轻X射线束硬化伪影效应[24,33].

      • 2.2.   系统设计

      2.2.1.   系统结构

    • Micro-CBCT扫描仪的几何结构见图3. 该系统包括一个带有PFXRL的微焦斑X射线源, 一个非晶硅平板探测器, 一个电动旋转平移台和一台用于信号采集、控制和图像重建的计算机. 微焦斑X射线源为功率50 W的钨阳极靶(美国牛津仪器公司), 最大管电压为50 kV, 最大管电流为1 mA, 焦斑尺寸为70 μm. 非晶硅平板探测器采用1536 × 1920阵列, 每个阵列的像素大小为127 μm, 闪烁屏面板大小为25 cm × 20 cm (2520DX, 美国Varex成像公司). 图3所示几何结构中, SOD为PFXRL出口焦斑到旋转中心的距离, SDD为PFXRL出口焦斑到非晶硅平板探测器成像面之间的距离. SDD与SOD的比值为成像系统的放大系数. 当非晶硅平板探测器成像面积确定时, 通过改变系统的放大系数可以改变系统的有效视场大小和成像分辨率大小. 将样品放置在旋转台上(SURUGA SEIKI, KRW06360, 日本), 由主机控制驱动逐步旋转360°, 投影图像的采集步进角度为1.8°. 由主机控制微焦斑X射线源开关, 设定所需的管电压和管电流大小.

      对于PFXRL, 能量依赖的入口/出口焦斑大小、传输效率、发散度/收敛度和强度增益在我们之前的工作中有详细的介绍[38]. 表1列出了PFXRL的基本参数. 图4给出了PFXRL的出口焦斑尺寸和收敛度随X射线能量的变化关系. 图5给出了PFXRL的传输效率和能量密度增益随X射线能量的变化关系.

      在本系统中, 由于样品是小动物, 选择放大系数在2.0左右. 系统经过几何校准后[39], SOD为390.3 mm, SDD为768.8 mm, 因此系统的放大系数为1.97, 有效像素大小为64.5 μm. 此时成像区域的直径约为10 cm, 足以实现小鼠的成像.

      • 2.2.2.   图像重构

    • 采用c++语言开发Windows平台下的Micro-CBCT系统运动控制数据采集与数据重构软件. 采集投影数据后, 对原始投影图像进行校正, 并采用经典的锥束FDK算法重建三维图像. 在重建过程中, 利用RAM-LAK内核函数进行滤波, 以保持图像的分辨率.

      • 2.3.   性能测试

    • 实验条件为: 微焦斑X射线源管电压为50 kV, 管电流为800 μA, 非晶硅平板探测器帧频为3.0 fps (frames per second), 共采集200幅投影图像.

      • 2.3.1.   系统分辨率

    • 调制传递函数(modulation transfer function, MTF)表示成像系统对不同频率输入信号的调制能力. 在空间域, 高频信号分量代表图像的细节信息. 因此, 高频分量的调制越高, 成像系统对样品细节信息的分辨能力越强.

      测量Micro-CBCT系统空间分辨率通常是采用点扩散函数(point spread function, PSF)法. 本文以直径为20 μm的金属丝为试样, 对使用和不使用PFXRL两种情况下的Micro-CBCT系统的PSF进行了测试[40]. 在整个扫描过程中将金属丝保持铅直, 避免PSF出现相位畸变. 在CT扫描和图像重建后, 选取感兴趣区域对金属丝进行过采样获得PSF, 通过计算PSF的傅里叶变换得到Micro-CBCT系统的MTF. 对比度分辨率表示系统所能分辨的对比度差别的能力. 本文采用直径为19 cm的圆柱形塑料套管水模, 通过改变管电压并采用统计学的方法来评估Micro-CBCT系统的对比度分辨率[41].

      • 2.3.2.   图像均匀性和噪声

    • 重建图像的均匀性意味着对于一个均匀体模, 测量的CT值不应随所选感兴趣区域(region of interest, ROI)的位置或体模相对于系统旋转中心的位置而变化. 其反映了Micro-CBCT系统测量的空间分布的一致性, 主要取决于探测器的均一性和X射线束的硬化效应等因素. 另外, FDK重构算法是一种近似算法, 也会带来重构图像的非均一性. 噪声是指均匀物体的影像中CT值在平均值上下的随机涨落. 图像噪声的存在, 会掩盖或降低图像的对比度. 实验选取蒸馏水作为实验对象, 将蒸馏水放入外径为85 mm, 内径为84 mm的塑料瓶中, 对其进行扫描并重建图像, 选取一定尺寸的ROI进行相应的数据分析来评估Micro-CBCT系统的成像均匀性和分析ROI区域CT数的方差来评估其噪声特征. 噪声的相对平稳性通过检测噪声在整个三维数据集上的均匀性来评估.

      • 2.3.3.   小动物实验

    • 将质量为300 g的雌性小鼠用10%氨基甲酸乙酯混合2%水合氯醛(0.9 mL/kg)进行腹腔麻醉, 并固定在旋转平台上进行CT扫描成像. 本文所述动物实验均经赣南医学院动物护理与利用委员会批准.

    3.   结果与讨论

      3.1.   系统分辨率

    • MTF测试结果如图6所示, 可以看出, 使用PFXRL后, 空间分辨率明显提高. 当MTF为10%时, 采用PFXRL的空间分辨率达到9.3 lp/mm, 而不使用PFXRL的空间分辨率为6.9 lp/mm. 根据(1)式—(3)式, 不使用PFXRL的Micro-CBCT系统的空间分辨率更低, 这是由X射线源焦斑大小和金属丝造成的模糊效应以及投影图像较少引起的. 同时, 基于PFXRL的Micro-CBCT系统的空间分辨率更高, 这是由于PFXRL能实现更小的X射线焦斑大小, 并抑制X射线散射, 进而有效地减少模糊效应.

      图7给出了测量的水模的对比度分辨率. 可以发现, 使用PFXRL的系统对比度分辨率比未使用PFXRL提高了数倍. 由(4)式可知, 对于不同尺寸的低对比度结构, 系统的对比度分辨率不仅取决于噪声, 还取决于系统的空间分辨率. PFXRL的使用可以有效提高系统的空间分辨率和照射物体的X射线强度增益(即提高系统的信噪比), 从而提高Micro-CBCT系统的对比度分辨率.

      • 3.2.   图像均匀性和噪声

    • 图8为均匀水模在使用和不使用PFXRL条件下以及对X射线微焦点源采用0.5 mm厚的铝片进行滤过后的Micro-CBCT系统成像扫描结果. 对于给定的重建和采集参数, 重建完成后, 将空气的CT数设置为0 (即, CT# 0), 将水的CT数设置为50 (即, CT# 50). 图8(a)为不使用PFXRL时, 重建的水模中平横断面图像及绿线对应的CT值; 图8(b) 为使用PFXRL时, 重建的水模中平横断面图像及绿线对应的CT值; 图8(c) 为采用0.5 mm厚的铝片滤过后, 重建的水模中平横断面图像及绿线对应的CT值. 显然, 不使用PFXRL时, 图像出现了明显的不均匀性, 存在典型的杯状伪影, 即均匀水模中间的CT值比边缘部分的CT值低, 造成中间黑, 边缘白的情况. 使用PFXRL和/或铝片滤过器后, 图像均显示出良好的均匀性, 但硬化伪影不能完全消除. PFXRL作为滤过器, 能有效吸收低能X射线, 减少硬化伪影, 提高图像的均一性.

      图9给出了均匀水模在使用PFXRL和不使用PFXRL (附加滤过: 0.5 mm Al)条件下噪声的大小和噪声在整个三维数据集上的均匀性. 重建水模的噪声图像中, 每个点代表5 × 5 × 1像素区域内CT数的方差($ {\sigma }^{2} $). 可以看出, 采用PFXRL的Micro-CBCT系统的噪声水平有较明显的降低. 对应标准偏差($ \sigma $)来表示, 使用PFXRL的Micro-CBCT系统达到CT# 5— CT# 8的噪声水平, 而不使用PFXRL(附加滤过: 0.5 mm Al)达到CT# 7—CT# 11的噪声水平. 从影响噪声的因素来看, ${\sigma }^{2}\propto {{f}^{3}}/({zD})$(其中, f为系统的空间分辨率, z为层厚, D为剂量), 采用PFXRL可以有效提高系统的空间分辨率(即f值减小)和照射物体的X射线强度增益(即D值增大), 因此采用PFXRL能有效地降低Micro-CBCT系统的噪声水平.

      • 3.3.   小动物实验

    • 麻醉小鼠成像结果如图10所示. 从图像上可以清楚地看到小鼠的结构, 使用PFXRL的图像的对比度明显高于未使用PFXRL的图像, 显示了精细的空间分辨率并提供了清晰的软组织边界和细节. 图11给出了使用PFXRL和未使用PFXRL的Micro-CBCT系统的麻醉小鼠的横断面切片. 在等效曝光条件下, 产生的图像在空间分辨率和对比灵敏度方面具有可比性. 基于PFXRL的Micro-CBCT图像展示了更高的空间分辨率, 并提供了更清晰的软组织边界轮廓和肠道细节. 这表明X光透镜可以有效地提高Micro-CBCT系统的成像均一性、对比度分辨率和空间分辨率. 图12给出了麻醉小鼠肺(图12(a)—(c))、肾(图12(d)—(f))和下脊柱(图12(g)—(i))区域的基于PFXRL的Micro-CBCT切片图像, 以说明其定量特性和成像均匀性质量. 每个图像的窗宽窗位设置不同, 以呈现出感兴趣的结构. 图12(a)—(c)的窗宽窗位设置是强调大鼠肺的特征, 肺的细节清晰可见, 说明了基于PFXRL的Micro-CBCT具有软组织高对比度灵敏性. 对于图12(d)—(f), 设置窗宽窗位来描绘脂肪和肌肉, 肾脏与周围的肌肉和脂肪清晰可见. 其他结构, 如胃、肠和肝也清晰可见. 每张图片中的十字线指示了大鼠左肾的位置. 图12(g)—(i)选择相同的窗宽窗位显示小鼠骨骼特征, 椎骨间隙和骨盆非皮质骨清晰可见, 说明了基于PFXRL的Micro-CBCT系统具有各向同性的高空间分辨率.

    4.   结　论

    • 本文设计并实现了一种基于PFXRL的Micro-CBCT系统用于小动物成像. 首先对系统进行了主要性能测试, 结果表明, 基于PFXRL的Micro-CBCT系统的空间分辨率、对比度分辨率和成像均一性均有较大提高. 使用PFXRL将所测Micro-CBCT系统的MTF的极限空间频率提高了1.35倍. PFXRL在所有空间频率上产生MTF的增加, 对于诊断运用, 可能比单独极限MTF的增加更有意义. 同时, PFXRL的散射抑制作用, 实现了2倍以上的对比度增强. 此外, 通过将样品放置在PFXRL的出口焦点附近, 增加了照射样品的X射线束流通量, 可将扫描时间缩短而不增加半影模糊. 基于X光透镜的Micro-CBCT系统的设计有助于提高图像分辨率、图像均一性和消除模糊效应, 将在小动物活体成像领域具备广泛的潜在应用前景.

    参考文献 (41)

目录

/

返回文章
返回