搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

高能电子三维成像技术实验研究

李双双 赵全堂 曹树春 冉朝晖 申晓康 赵书俊 张子民

李双双, 赵全堂, 曹树春, 冉朝晖, 申晓康, 赵书俊, 张子民. 高能电子三维成像技术实验研究. 物理学报, 2021, 70(18): 184204. doi: 10.7498/aps.70.20210686
引用本文: 李双双, 赵全堂, 曹树春, 冉朝晖, 申晓康, 赵书俊, 张子民. 高能电子三维成像技术实验研究. 物理学报, 2021, 70(18): 184204. doi: 10.7498/aps.70.20210686
Li Shuang-Shuang, Zhao Quan-Tang, Cao Shu-Chun, Ran Zhao-Hui, Shen Xiao-Kang, Zhao Shu-Jun, Zhang Zi-Min. Experimental demonstration of three-dimensional high energy electron radiography. Acta Phys. Sin., 2021, 70(18): 184204. doi: 10.7498/aps.70.20210686
Citation: Li Shuang-Shuang, Zhao Quan-Tang, Cao Shu-Chun, Ran Zhao-Hui, Shen Xiao-Kang, Zhao Shu-Jun, Zhang Zi-Min. Experimental demonstration of three-dimensional high energy electron radiography. Acta Phys. Sin., 2021, 70(18): 184204. doi: 10.7498/aps.70.20210686

高能电子三维成像技术实验研究

李双双, 赵全堂, 曹树春, 冉朝晖, 申晓康, 赵书俊, 张子民

Experimental demonstration of three-dimensional high energy electron radiography

Li Shuang-Shuang, Zhao Quan-Tang, Cao Shu-Chun, Ran Zhao-Hui, Shen Xiao-Kang, Zhao Shu-Jun, Zhang Zi-Min
PDF
HTML
导出引用
  • 高能电子成像技术被首次提出作为温稠密物质和惯性约束聚变实验研究的高时空分辨诊断工具之一, 现已通过前期实验证明其对中尺度科学诊断的可行性. 为了进一步提高高能电子成像技术诊断样品的能力, 来获取样品内部信息, 将高能电子成像技术和三维重建算法结合, 提出了高能电子三维成像技术. 本文主要通过实验研究了高能电子三维成像技术的可行性. 不同三维重建算法重建样品的结果首次证实了高能电子三维成像技术的可行性, 使用的三维重建算法包括滤波反投影算法、迭代算法-代数重建技术和联立代数重建技术, 最终重建的x, y, z方向上的不同重建切片图像清楚地显示了样品的详细结构. 实验证实的高能电子三维成像技术将有利于拓展高能电子成像技术的应用领域, 尤其是在中尺度科学领域.
    High energy electron radiography (HEER) proposed first for real-time high spatial and temporal resolution diagnosis of warm dense matter (WDM) and inertial confinement fusion (ICF) has proved experimentally feasible for mesoscale sciences diagnosis. Until now, the spatial resolution of the images close to 1 μm has been reached experimentally which is better than that of X-rays and neutron radiography. However, traditional HEER obtains two-dimensional images which cannot accurately present the three-dimensional structure of the sample. To further improve the capability of HEER to diagnose and obtain the internal information of samples, three-dimensional high energy electron radiography (TDHEER) was put forward by combining HEER with three-dimensional (3D) reconstruction tomography technology. The validity and usage of the TDHEER method have been confirmed through simulation of the fully 3D diagnostic of static mesoscale sample. This paper focuses mainly on the experimental demonstration of the 3D high energy electron radiography. The feasibility of TDHEER is for the first time confirmed by the results achieved with different 3D reconstruction algorithms. The 3D reconstruction algorithms, analytical algorithm-filtered back projection (FBP), iterative algorithms-algebraic reconstruction technique (ART), and simultaneous algebraic reconstruction technique (SART) are used here. In this experiment, the less projected data are used, so it takes the less time to obtain two-dimensional (2D) HEER images and the reconstruction. In order to spend the time as little as possible and obtain the satisfactory quality of reconstruction result, there are three groups of projected image sets, 180, 36 and 18, acquired in our experiment. When all three algorithms are adopted in 180 projected images, the reconstructed images show that all three algorithms FBP, ART and SART are feasible for TDHEER. The different reconstructed slice images of the sample in X-, Y-, and Z- direction clearly show the detailed structure of the sample. The images reconstructed by ART and SART algorithm are equivalent. Comparing with ART and SART, the reconstruction results by FBP can show more details, but there are some artifacts. Because the 36 2D HEER images fail to satisfy the Nyquist sampling theory, the analytic algorithm FBP is not used. Taking the result of FBP reconstructed by 180 images as a standard reference to compare the result of ART with the results of SART, the images reconstructed by the SART algorithm are closer to the original images. Testing 18 images, the results of the ART and SART both have lots of artifacts but the SART algorithm spends less time in reconstruction. As fewer projected images are used, more artifacts are found in the reconstructed images. Therefore, it is advantageous to combine the SART algorithm with 36 HEER projected images, which obtains high-quality reconstruction images and spends less time. The feasibility of TDHEER is confirmed experimentally for the first time and all three dimensions of the sample structures are obtained. Of the three different 3D reconstruction algorithms, the SART algorithm is the most suitable for reconstructing the few-view images. The TDHEER technology will extend HEER’s application fields, especially for mesoscale sciences.
      PACS:
      42.30.Wb(Image reconstruction; tomography)
      41.85.Ja(Particle beam transport)
      41.75.-i(Charged-particle beams)
      07.78.+s(Electron, positron, and ion microscopes; electron diffractometers)
      通信作者: 赵全堂, zhaoquantang@impcas.ac.cn ; 赵书俊, zhaosj@zzu.edu.cn
    • 基金项目: 国家重点研发计划(批准号: 2019YFA0404900)、国家自然科学基金(批准号: 11875303)和中国科学院国际合作局对外合作重点项目(批准号: 113462KYSB20160036)资助的课题
      Corresponding author: Zhao Quan-Tang, zhaoquantang@impcas.ac.cn ; Zhao Shu-Jun, zhaosj@zzu.edu.cn
    • Funds: Project supported by the National Key R&D Program of China (Grant No. 2019YFA0404900), the National Natural Science Foundation of China (Grant No. 11875303), and the Key Program of the International Partnership of Bureau of International Cooperation Chinese Academy of Sciences (Grant No. 113462KYSB20160036)

    近年来, 有关高能密度物理(high energy density physics, HEDP)和惯性约束聚变(inertial confinement fusion, ICF)的研究日益受到重视并得到发展, 特别是其诊断方法是面临的一大挑战[1]. 为了对HEDP和ICF进行深入研究, 提出了一种具有高时空分辨的高能电子成像技术(high energy electron radiography, HEER)[2]. 目前高能电子成像技术在仿真和实验中都得到了很好的研究[3-8]. 由于高能电子束穿透力强, HEER可以用于较厚样品(几十微米到毫米)的诊断, 现阶段通过实验表明HEER的空间分辨率接近1 μm[9].

    目前许多射线照相诊断工具都基于X射线、中子和质子束. 以软X射线为基础的X射线显微镜是研究生物材料的理想选择. 硬X射线(光子能量大于10 keV)由于其超强的穿透力被用于临床医学诊断疾病、工程和工业检测缺陷. X射线点投影背光成像已被用于高能量密度物质(high energy density materials, HEDM)实验[10], 其空间分辨率约为10 μm[11]. 在SG-II升级激光设备上, 用短脉冲激光驱动的硬X射线进行了双壳内爆实验[12], 但通过X射线获得的图像空间分辨率不及HEER. 短脉冲中子被应用于了解温稠密物质的动态过程[13]. 中子射线照相的优势在于对低原子序数元素的材料具有极高的敏感性[14]. HEDM的材料一般由重金属组成, 中子射线对其不敏感, 且中子与金属之间的反应截面太小, 无法有效地对HEDM成像. 中子成像诊断是NIF进行惯性聚变研究的工具, 但其成像的空间分辨率大于10 μm[15], 大于HEER的空间分辨率. 高能质子射线照相术在诊断HEDM时显示了其优异的性能[16]. 然而在实验室中不容易产生超短脉冲(ps量级)的高能质子束, 而高能电子成像具有高时空分辨率的优势. 由于质子与物质之间的相互作用相对单一, 因此很难反映物质的不同性质. 此外对比高能质子加速器, 高能电子加速器和相应的成像系统造价更低廉.

    高能电子三维重建(three-dimensional high energy electron radiography, TDHEER)技术将HEER与三维(three dimensions, 3D)重建算法相结合, TDHEER技术的有效性和可行性已通过静态中尺度样品的模拟实验研究得到证实[17]. 为了通过实验验证该技术, 在本文研究中进行了毫米级样本的重建. 根据Nyquist采样理论, 重建使用的投影数据越少, 则同种算法重建所花费的时间就越少, 但会增加伪影降低图像质量. 重建的目标是在满足重建结果质量的同时, 花费尽可能少的时间. 因此在实验中分别使用了三组不同数目的投影图像集: 180, 36和18, 以及三种不同的重建算法, 包括解析和迭代算法来重建样本. 结果发现, 使用不同的算法及投影集重建的图像质量和重建时间是不同的. 为提高重建图像的质量并且减少重建时间, 有必要选择合适的算法.

    TDHEER的主要思想是使用3D重建算法重建获取的2D HEER数据以获得样本的切片图像. HEER的原理如图1所示: 由电子直线加速器产生的高能电子束(50 MeV)穿过样品, 电子束与样品相互作用后将富含样品的信息, 再经由具有一定放大倍数的四极磁透镜系统在荧光屏上实现点对点成像, 最后通过CCD相机采集高能电子图像. 电子束与物质的相互作用主要包括电离、激发、轫致辐射, 这三种相互作用都会导致电子能量的损失[18]. 由于实验采用的是高能电子束, 轫致辐射会更加明显, 并且由其带来的能量损失占据主导地位, 其能量损失可表示为

    图 1  HEER的布局原理图
    Fig. 1.  Layout and principle of HEER.
    $$f(E) = \frac{1}{{{E_0}}}\Bigg\{ {{{\left[ {\ln \left( {\frac{{{E_0}}}{E}} \right)} \right]}^{\tfrac{{4\delta x}}{{3{X_0}}} - 1}}\bigg/\varGamma \left( {\frac{{4\delta x}}{{3{X_0}}}} \right)} \Bigg\}, $$ (1)
    $$ \frac{1}{{{X_0}}} = 4\alpha \frac{{{N_{\text{a}}}}}{A}Z(Z + 1)r_{\text{e}}^2\ln \left(183{Z^{ - \tfrac{1}{3}}}\right), $$ (2)

    其中${E_0}$, E分别为电子的初始能量和穿过物质后的能量; $\delta $是物质材料的密度; x是物质的厚度; $\alpha {\text{ = }}1/137$为精细结构常数; ${r_e} = {{ {e}}^2}/{m_{\text{e}}}{c^2}$是古典电子半径; ${N_{\text{a}}}$, ZA分别是阿伏伽德罗常数, 物质原子的原子序数和质量数; ${X_0}$是电子辐射长度, 其表示为当物体厚度为${X_0}$时, 电子在物质中因辐射损失而使得能量衰减至原初始能量的$1/{ {\rm e}}$.

    电子束与物质相互作用后的角度发散主要由于电子的多重散射, 可由Molière模型描述多库仑散射作用下出射电子的角度分布:

    $$ f(\phi ) = \frac{1}{{{\phi _0}(t)\sqrt {2{\text{π}}} }}{{\text{e}}^{ - \tfrac{1}{2}{{\left( {\tfrac{\phi }{{{\phi _0}(t)}}} \right)}^2}}} , $$ (3)
    $$ {\phi _0}(t) = \frac{{13.6\;{\text{MeV}}}}{{\beta cp}}\sqrt t (1 + 0.038\ln (t)) ,$$ (4)

    其中${\phi _0}(t)$$\phi $分别为电子的初始角度和散射后的角度, p是电子的动量, $\beta c$是电子的速度, t是物质的厚度x与电子辐射长度${X_0}$之比. 实验中电子束的能散较低(约0.1%), 主要根据透射电子的散射角分布对物质进行成像.

    成像透镜系统中存在傅里叶平面, 用于筛除由多库仑散射导致散角过大的电子, 从而提高空间分辨率. 荧光屏是Ce:YAG屏幕, 实现将电信号转换为光信号. 当CCD相机直接对准荧光屏成像时, 由于受到强辐射采集到的图像会有大量噪声. 因此在荧光屏后添加一面反光镜, 变换光路, 使得光线偏转90º后再由相机采集. 3D重建算法包括非迭代和迭代算法, 已广泛应用于医学图像重建[19]、工业CT[20]和其他3D射线照相中. TDHEER类似于电子断层扫描技术, 因此对TDHEER采用电子断层扫描中使用的算法, 即滤波反投影(filtering back projection, FBP)[21]、迭代算法-代数重建技术(algebraic reconstruction technique, ART)[22]和联立代数重建技术(simultaneous algebraic reconstruction technique, SART)[23]来重建二维HEER数据.

    采集2D HEER数据的实验原理如图2所示, 沿Z轴传播的高能电子束穿过目标样品的X -Y平面, 经由成像透镜系统聚焦透射电子实现点对点成像, 最终获得X -Y平面2D HEER图像. 将样品放置在旋转平台上, 该平台可以绕Y轴(± 90°)旋转, 每次旋转的角度间隔为1°, 每旋转一次采集一幅2D HEER图像, 通过CCD相机收集了从0°到179°的共计180张图像. 为便于样品制备, 实验样品是手表的1个零件如图3所示, 其尺寸如下: 长度 × 厚度 × 高度为3.34 mm × 0.6 mm × 2.18 mm. 图3(a)图3(b)分别是样品的正视图和左视图, 以其长度为X轴, 高度为Y轴, 厚度为Z轴, 建立3D坐标系如图3(c)所示. 该实验在兰州高能电子成像实验平台上进行, 实验中使用的电子束参数如下: 束团能量为50 MeV, 能散为0.1%, 归一化发射度约为8.0 mm·mrad, 束团电荷量约为100 pC, 均方根束长约为1.5 ps[24].

    图 2  获取2D HEER数据的实验原理
    Fig. 2.  Experimental principle of acquiring 2D HEER data.
    图 3  实验样品 (a) 正视图; (b) 左视图; (c) 3D坐标系
    Fig. 3.  Experimental sample: (a) The front view; (b) the left view; (c) the 3D coordinate diagram.

    在实验中, 使用优化后的四极成像透镜系统来获得2D HEER图像, 其放大倍数为–2.88[5]. 图4展示了样品分别在0°, 45°, 135°和180°旋转角度下拍摄的2D HEER图像. 分析图像发现含有噪声, 其主要来源是CCD的固有噪声以及一些由X射线引起的噪声. 由于在图像采集过程中这些噪声不可避免地引入并影响最终的重建质量, 因此首先分析噪声类别并选择合适的滤波器将其滤除. 经分析得出噪声的主要成分是椒盐噪声, 分别用3种滤波器: 均值滤波器[25]、中值滤波器[26]和小波滤波器[27]进行滤波降噪, 同一图像分别经过3个滤波器处理后的信噪比(signal to noise ratio, SNR)结果分别为19.6165, 19.6411和16.6871. 图像的SNR越高, 图像质量越好, 因此选择中值滤波器来处理图像. 由于重建算法要求投影图像具有相同的长度和宽度, 因此将原始大小为2448 × 2048像素的图像裁剪到仅包含样本的1670 × 1670像素的区域. 此外当图像分辨率相同时, 图像尺寸越小, 重建时间越短, 将尺寸1670 × 1670压缩为167 × 167像素, 以加快重建速度. 最终获得用于重建的180幅2D HEER图像, 其尺寸为167 × 167像素.

    图 4  2D HEER图像 (a)−(d) 旋转角度分别为0º, 45º, 135º和 180º
    Fig. 4.  2D HEER images: (a)−(d) Rotation angle is 0º, 45º, 135º, 180º, respectively.

    解析重建算法FBP能够克服常规反投影的局限性, 当FBP被应用于180幅2D HEER图像的全角度投影重建时, 其重建结果要优于稀疏角度投影重建[20]. 因此采用由HEER平台收集的样本180幅HEER图像(从0º到179º, 间隔1º)用于三维重建. 最基础的迭代重建算法是ART, 当使用稀疏投影数据时, 诸如ART之类的迭代算法的重建结果要优于FBP[28], 但ART要比FBP花费更多的重建时间. 当使用与ART相同的参数时, 采用SART可获得与ART大致相同的图像质量但其重建时长在FBP和ART之间[29]. 在重建过程中, 调整ART和SART的松弛因子、迭代次数和迭代初始值, 以获得良好的重建结果. FBP, ART和SART算法重建的样品分成X, YZ三个方向, 每个方向中共有236个重建切片, 切片大小为6.24 cm × 6.24 cm. 在X方向上重建的从X –到X +的第95, 134和145层的切片如图5所示. 在Y方向上重建的从Y –到Y +的第89, 124和150层的切片如图6所示. 在Z方向上重建的从Z –到Z +的第52, 120和155层的切片如图7所示.

    图 5  X–到X+在第95, 134和145层使用不同算法重建的切片的结果 (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART
    Fig. 5.  Results of reconstructed slices with different algorithms, at the 95th, 134th and 145th layers from X– to X+: (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART.
    图 6  Y–到Y+的第89, 124和150层使用不同算法重建的切片的结果 (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART
    Fig. 6.  Results of reconstructed slices with different algorithms, at the 89th, 124th and 150th layers from Y– to Y+: (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART.
    图 7  Z–到Z+在第52, 122和155层使用不同算法重建的切片的结果 (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART
    Fig. 7.  Results of reconstructed slices with different algorithms, at the 52nd, 122th and 155th layers from Z– to Z+: (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART.

    根据图57中的重建图像可知, 三种算法FBP, ART和SART都可用于高能电子三维成像. 与ART和SART相比, FBP的重建结果可以显示更多细节, 但存在一些伪影. 分析FBP, ART和SART重建的Y方向上的第124层切片中第151列的像素灰度分布图, 如图8(b)所示, 即图8(a)中红线位置的像素灰度分布. 根据曲线走势分布表明, ART和SART算法重建的图像是等效的.

    图 8  (a) FBP重建的从Y–到Y +第124个切片图像; (b) 是(a)中红线的像素灰度分布图
    Fig. 8.  (a) The 124th slice image from Y– to Y+ reconstructed by FBP; (b) the pixels grayscale of the red line of panel (a).

    为了减少噪声影响、2D HEER数据采集时间及加快重建速度, 使用从0°到179°间隔5º的共36张2D HEER图像数据集来重建. 由于36张2D HEER图像无法满足Nyquist采样理论, 使用解析算法FBP会产生影响图像质量的伪影. 与解析类算法(假定采集位置连续)相比, 基于ART的方法都是非统计性的, 能更好地对采集过程的几何模型进行建模. 因此基于ART的方法可以更好地处理稀疏数据和采样位置不规则的采样. 为了定量分析ART和SART的结果, 以180张图像重建的FBP结果作为基准, 使用两个评价指标作为参考.

    归一化均方距离准则d :

    $$ d = \sqrt {\frac{{\displaystyle\sum\limits_{u = 1}^N {\sum\limits_{v = 1}^N {{{\left( {{t_{u,v}} - {r_{u,v}}} \right)}^2}} } }}{{\displaystyle\sum\limits_{u = 1}^N {\sum\limits_{v = 1}^N {{{\left( {{t_{u,v}} - \overline t } \right)}^2}} } }}} . $$ (5)

    归一化平均绝对距离准则r :

    $$ r = \left({{\displaystyle\sum\limits_{u = 1}^N {\sum\limits_{v = 1}^N {\left| {{t_{u,v}} - {r_{u,v}}} \right|} } }}\right)\bigg/\left({\displaystyle\sum\limits_{u = 1}^N {\sum\limits_{v = 1}^N {\left| {{t_{u,v}}} \right|} } }\right). $$ (6)

    在(5)式和(6)式中, ${t_{u, v}}$${r_{u, v}}$分别表示参考图像和重建图像的第u行和第v列的像素密度; N是每行或每列像素总数; $\overline t $是参考图像的像素平均值, 图像的像素数为236 × 236. dr的值反映了重建图像的质量, 值越小则重建图像与原始图像越接近. ART和SART重建的Y方向上从Y –到Y +的第100层和第127层的切片如图9所示.

    图 9  Y–到Y+的第100层和第127层使用不同的算法得出的结果 (a), (b) ART; (c), (d) SART
    Fig. 9.  Reconstructed slices results with different algorithms, at the 100th and 127th layers from Y– to Y+: (a), (b) ART; (c), (d) SART.

    分析表1数据, 当重建同一张图片时, SART算法重建的图像的dr值小于ART重建的图像的dr值, 即SART算法重建的图像更接近原始图像. 将图6图9进行比较, 由于使用较少的2D HEER图像数量, 所以图9中的伪影更加严重.

    表 1  ART, SART重建切片的dr
    Table 1.  Value of d and r of the image reconstructed by ART and SART.
    切片位置重建算法dr
    Y 100thART11.07360.8820
    SART10.90600.8672
    Y 127thART9.09470.9076
    SART8.95200.8933
    下载: 导出CSV 
    | 显示表格

    为了在获得高质量图像的同时尽可能缩短整个重建过程, 因此测试了旋转间隔为10°的2D HEER数据集. 不同算法重建的Y轴的第127个切片的结果如图10所示, 图10(a)图10(b)分别是ART和SART算法的结果, 图10(a)图10(b)相似, 尽管使用不同的重建算法两者都具有许多伪影, 但此时SART算法的重建时间更短. 图10(c)对比了同种算法重建不同投影子集时图像的像素灰度分布, 与5°间隔的投影数据集相比, 10°间隔的投影数据集重建的图像具有更多的伪影. 这与重建中使用的2D HEER数据越少, 重建图像质量越差的规则相符. 因此, 将SART算法与间隔5°的2D HEER数据相结合更有利于本实验数据的重建, 在获得高质量的重建图像的同时花费更少的时间.

    图 10  Y–到Y +在第127层重建切片的结果 (a) ART; (b) SART; (c) 是(b)中红线的像素灰度分布
    Fig. 10.  Results of reconstructed slices at the 127th layer from Y– to Y+: (a) ART; (b) SART; (c) the pixels grayscale of the red line of panel (b).

    通过本实验研究首次证明了高能电子三维成像技术的可行性, 获得了样品的3个维度切片. 通过将不同的3D重建算法应用于TDHEER, 表明了大间隔采样角度下重建的可行性, 并对比了TDHEER中3种不同算法的优缺点. 结果表明, SART算法更适合重建本实验平台获取的高能电子图像. 但由于实验设置的原因, 例如束流不稳定性、成像采集过程和样品的安装位置以及可旋转平台的精度限制了最终的重建图像质量. 为了进一步优化TDHEER, 需要研究束流稳定性, 样品准直和旋转平台精度. 此外通过将SART与基于总变分的算法结合使用, 可以改进SART算法以减少伪影[30-32], 进而改善图像质量. 相信高能电子三维成像技术的实验验证能拓展HEER技术的应用领域, 尤其在中尺度科学领域中的应用.

    感谢兰州大学李公平教授和张催博士在三维重建算法方面提供的帮助.

    [1]

    Wei G, Qiu J, Jing C 2014 Pro. SPIE 9211 921104

    [2]

    Zhao Y T, Zhang Z, Gai W, et al. 2016 Laser Part. Beams 34 338Google Scholar

    [3]

    Merrill F, Harmon F, Hunt A, Mariam F, Morley K, Morris C, Saunders A, Schwartz C 2007 Nucl. Instrum. Methods Phys. Res., Sect. B 261 382Google Scholar

    [4]

    Zhao Q T, Cao S C, Cheng R, Shen X K, Zhang Z M, Zhao Y T, Gai W, Du Y C 2014 Proceedings of the LINAC2014 Geneva, Switzerland, August 31–September 5, 2014 p76

    [5]

    Zhao Q T, Cao S C, Liu M, et al. 2016 Nucl. Instrum. Methods Phys. Res., Sect. A 832 144Google Scholar

    [6]

    Zhou Z, Du Y C, Cao S C, et al. 2018 Phys. Rev. Accel. Beams 21 074701Google Scholar

    [7]

    Zhao Q T, Cao S. C, Cheng R, et al. 2018 Laser Part. Beams 36 313Google Scholar

    [8]

    Zhao Q T, Cao S C, Shen X K, Wang Y R, Zong Y, Xiao J H, Zhu Y L, Zhou Y W, Liu M, Cheng R, Zhao Y T, Zhang Z M, Gai W 2017 Laser Part. Beams 35 579Google Scholar

    [9]

    Zhou Z, Fang Y, Chen H, Wu Y P, Du Y C, Yan L X, Tang C X, Huang W H 2019 Phys. Rev. Appl. 11 034068Google Scholar

    [10]

    Maddox B R, Park H S, Remington B A, et al. 2011 Phys. Plasmas 18 168Google Scholar

    [11]

    Park H S, Maddox B R, Giraldez E, et al. 2008 Phys. Plasmas 15 3048Google Scholar

    [12]

    Tian C, Yu M H, Shan L Q, Wu Y C, Zhang T K, Bi B, Zhang F, Zhang Q Q, Liu D X, Wang W W, Yuan Z Q, Yang S Q, Yang L, Zhou W M, Gu Y Q, Zhang B H 2019 Nucl. Fusion 59 046012Google Scholar

    [13]

    Higginson D P, Vassura L, Gugiu M, et al. 2015 Phys. Rev. Lett. 115 054802Google Scholar

    [14]

    Strobl M, Manke I, Kardjilov N, Hilger A, Dawson M, Banhart J 2009 J. Phys. D: Appl. Phys. 42 243001Google Scholar

    [15]

    Merrill F E, Bower D, Buckles R, et al. 2012 Rev. Sci. Instrum. 83 051003Google Scholar

    [16]

    King N S P, Ables E, Adams K, et al. 1999 Nucl. Instrum. Methods Phys. Res., Sect. A 424 84Google Scholar

    [17]

    Zhao Q T, Ma Y Y, Xiao J H, Cao S C, Zhang Z M 2019 Appl. Sci. 9 3764Google Scholar

    [18]

    谢一冈, 陈昌, 王曼, 吕军光, 孟祥承, 王锋, 顾树棣, 过雅南 2003 粒子探测器与数据获取 (北京: 科学出版社) 第5−16页

    Xie Y G, Cheng C, Wang M, Lv J G, Meng X G, Wang F, Gu S D, Guo Y N 2003 Particle Detector and Data Acquisition (Beijing: Science Press) pp5−16 (in Chinese)

    [19]

    Padole A, Khawaja Rd Ali, Kalra M K, Singh S 2015 Am. J. Roentgenol. 204 384Google Scholar

    [20]

    Chen B X, Yang M, Zhang Z, Bian J G, Han X, Sidky E, Pan X C 2014 Biochim. Biophys. Acta 1581 1856Google Scholar

    [21]

    Schofield R, King L, Tayal U, Castellano I, Stirrup J, Pontana F, Earls J, Nicol E 2020 J. Cardiovasc Comput Tomogr. 14 219Google Scholar

    [22]

    Gordon R, Bender R, Herman G T 1970 J. Theor. Biol. 29 471Google Scholar

    [23]

    Trampert J, Leveque J J 1990 J. Geophys. Res. B:Solid Earth 95 12553Google Scholar

    [24]

    Zhu Y L, Yuan P, Cao S C, et al. 2018 Nucl. Instrum. Methods Phys. Res., Sect. A 911 74Google Scholar

    [25]

    Blackledget J M 2006 Digital Signal Processing (2nd Ed.) (Cambridge: Woodhead Press) pp522−540

    [26]

    拉斐尔C, 理查E, 史蒂文L 著 (阮秋琦 译) 2014 数字图像处理(MATLAB版) (第二版) (北京: 电子工业出版社) 第45−48页

    Rafael C, Richard E, Steven L (translated by Ruan Q Q) 2014 Digital Image Processing Using MATLAB (2nd Ed.) (Beijing: Electronics Industry Press) pp45−48 (in Chinese)

    [27]

    Vetterli M, Herley C 1992 IEEE Trans. Acoust., Speech, Signal Process 40 2207Google Scholar

    [28]

    Zhang Y, Zhang W H, Lei Y J, Zhou J L 2014 J. Opt. Soc. Am. A 31 981Google Scholar

    [29]

    Soleimani M., Pengpen T 2015 Philos Trans. R. Dov. London, Ser. A 373 20140399Google Scholar

    [30]

    Sidky E Y, Kao C M, Pan X C 2009 J. X-Ray Sci. Technol. 14 119Google Scholar

    [31]

    Liu Y, Ma J H, Fan Y, Liang Z R 2012 Phys. Med. Biol. 57 7923Google Scholar

    [32]

    Deng L Z, Mi D L, He P, Feng P, Yu P W, Chen M Y, Li Z C, Wang J, Wei B 2015 Bio-med. Mater. Eng. 26 1685Google Scholar

    期刊类型引用(0)

    其他类型引用(1)

  • 图 1  HEER的布局原理图

    Fig. 1.  Layout and principle of HEER.

    图 2  获取2D HEER数据的实验原理

    Fig. 2.  Experimental principle of acquiring 2D HEER data.

    图 3  实验样品 (a) 正视图; (b) 左视图; (c) 3D坐标系

    Fig. 3.  Experimental sample: (a) The front view; (b) the left view; (c) the 3D coordinate diagram.

    图 4  2D HEER图像 (a)−(d) 旋转角度分别为0º, 45º, 135º和 180º

    Fig. 4.  2D HEER images: (a)−(d) Rotation angle is 0º, 45º, 135º, 180º, respectively.

    图 5  X–到X+在第95, 134和145层使用不同算法重建的切片的结果 (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART

    Fig. 5.  Results of reconstructed slices with different algorithms, at the 95th, 134th and 145th layers from X– to X+: (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART.

    图 6  Y–到Y+的第89, 124和150层使用不同算法重建的切片的结果 (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART

    Fig. 6.  Results of reconstructed slices with different algorithms, at the 89th, 124th and 150th layers from Y– to Y+: (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART.

    图 7  Z–到Z+在第52, 122和155层使用不同算法重建的切片的结果 (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART

    Fig. 7.  Results of reconstructed slices with different algorithms, at the 52nd, 122th and 155th layers from Z– to Z+: (a)−(c) FBP; (d)−(f) ART; (g)−(i) SART.

    图 8  (a) FBP重建的从Y–到Y +第124个切片图像; (b) 是(a)中红线的像素灰度分布图

    Fig. 8.  (a) The 124th slice image from Y– to Y+ reconstructed by FBP; (b) the pixels grayscale of the red line of panel (a).

    图 9  Y–到Y+的第100层和第127层使用不同的算法得出的结果 (a), (b) ART; (c), (d) SART

    Fig. 9.  Reconstructed slices results with different algorithms, at the 100th and 127th layers from Y– to Y+: (a), (b) ART; (c), (d) SART.

    图 10  Y–到Y +在第127层重建切片的结果 (a) ART; (b) SART; (c) 是(b)中红线的像素灰度分布

    Fig. 10.  Results of reconstructed slices at the 127th layer from Y– to Y+: (a) ART; (b) SART; (c) the pixels grayscale of the red line of panel (b).

    表 1  ART, SART重建切片的dr

    Table 1.  Value of d and r of the image reconstructed by ART and SART.

    切片位置重建算法dr
    Y 100thART11.07360.8820
    SART10.90600.8672
    Y 127thART9.09470.9076
    SART8.95200.8933
    下载: 导出CSV
  • [1]

    Wei G, Qiu J, Jing C 2014 Pro. SPIE 9211 921104

    [2]

    Zhao Y T, Zhang Z, Gai W, et al. 2016 Laser Part. Beams 34 338Google Scholar

    [3]

    Merrill F, Harmon F, Hunt A, Mariam F, Morley K, Morris C, Saunders A, Schwartz C 2007 Nucl. Instrum. Methods Phys. Res., Sect. B 261 382Google Scholar

    [4]

    Zhao Q T, Cao S C, Cheng R, Shen X K, Zhang Z M, Zhao Y T, Gai W, Du Y C 2014 Proceedings of the LINAC2014 Geneva, Switzerland, August 31–September 5, 2014 p76

    [5]

    Zhao Q T, Cao S C, Liu M, et al. 2016 Nucl. Instrum. Methods Phys. Res., Sect. A 832 144Google Scholar

    [6]

    Zhou Z, Du Y C, Cao S C, et al. 2018 Phys. Rev. Accel. Beams 21 074701Google Scholar

    [7]

    Zhao Q T, Cao S. C, Cheng R, et al. 2018 Laser Part. Beams 36 313Google Scholar

    [8]

    Zhao Q T, Cao S C, Shen X K, Wang Y R, Zong Y, Xiao J H, Zhu Y L, Zhou Y W, Liu M, Cheng R, Zhao Y T, Zhang Z M, Gai W 2017 Laser Part. Beams 35 579Google Scholar

    [9]

    Zhou Z, Fang Y, Chen H, Wu Y P, Du Y C, Yan L X, Tang C X, Huang W H 2019 Phys. Rev. Appl. 11 034068Google Scholar

    [10]

    Maddox B R, Park H S, Remington B A, et al. 2011 Phys. Plasmas 18 168Google Scholar

    [11]

    Park H S, Maddox B R, Giraldez E, et al. 2008 Phys. Plasmas 15 3048Google Scholar

    [12]

    Tian C, Yu M H, Shan L Q, Wu Y C, Zhang T K, Bi B, Zhang F, Zhang Q Q, Liu D X, Wang W W, Yuan Z Q, Yang S Q, Yang L, Zhou W M, Gu Y Q, Zhang B H 2019 Nucl. Fusion 59 046012Google Scholar

    [13]

    Higginson D P, Vassura L, Gugiu M, et al. 2015 Phys. Rev. Lett. 115 054802Google Scholar

    [14]

    Strobl M, Manke I, Kardjilov N, Hilger A, Dawson M, Banhart J 2009 J. Phys. D: Appl. Phys. 42 243001Google Scholar

    [15]

    Merrill F E, Bower D, Buckles R, et al. 2012 Rev. Sci. Instrum. 83 051003Google Scholar

    [16]

    King N S P, Ables E, Adams K, et al. 1999 Nucl. Instrum. Methods Phys. Res., Sect. A 424 84Google Scholar

    [17]

    Zhao Q T, Ma Y Y, Xiao J H, Cao S C, Zhang Z M 2019 Appl. Sci. 9 3764Google Scholar

    [18]

    谢一冈, 陈昌, 王曼, 吕军光, 孟祥承, 王锋, 顾树棣, 过雅南 2003 粒子探测器与数据获取 (北京: 科学出版社) 第5−16页

    Xie Y G, Cheng C, Wang M, Lv J G, Meng X G, Wang F, Gu S D, Guo Y N 2003 Particle Detector and Data Acquisition (Beijing: Science Press) pp5−16 (in Chinese)

    [19]

    Padole A, Khawaja Rd Ali, Kalra M K, Singh S 2015 Am. J. Roentgenol. 204 384Google Scholar

    [20]

    Chen B X, Yang M, Zhang Z, Bian J G, Han X, Sidky E, Pan X C 2014 Biochim. Biophys. Acta 1581 1856Google Scholar

    [21]

    Schofield R, King L, Tayal U, Castellano I, Stirrup J, Pontana F, Earls J, Nicol E 2020 J. Cardiovasc Comput Tomogr. 14 219Google Scholar

    [22]

    Gordon R, Bender R, Herman G T 1970 J. Theor. Biol. 29 471Google Scholar

    [23]

    Trampert J, Leveque J J 1990 J. Geophys. Res. B:Solid Earth 95 12553Google Scholar

    [24]

    Zhu Y L, Yuan P, Cao S C, et al. 2018 Nucl. Instrum. Methods Phys. Res., Sect. A 911 74Google Scholar

    [25]

    Blackledget J M 2006 Digital Signal Processing (2nd Ed.) (Cambridge: Woodhead Press) pp522−540

    [26]

    拉斐尔C, 理查E, 史蒂文L 著 (阮秋琦 译) 2014 数字图像处理(MATLAB版) (第二版) (北京: 电子工业出版社) 第45−48页

    Rafael C, Richard E, Steven L (translated by Ruan Q Q) 2014 Digital Image Processing Using MATLAB (2nd Ed.) (Beijing: Electronics Industry Press) pp45−48 (in Chinese)

    [27]

    Vetterli M, Herley C 1992 IEEE Trans. Acoust., Speech, Signal Process 40 2207Google Scholar

    [28]

    Zhang Y, Zhang W H, Lei Y J, Zhou J L 2014 J. Opt. Soc. Am. A 31 981Google Scholar

    [29]

    Soleimani M., Pengpen T 2015 Philos Trans. R. Dov. London, Ser. A 373 20140399Google Scholar

    [30]

    Sidky E Y, Kao C M, Pan X C 2009 J. X-Ray Sci. Technol. 14 119Google Scholar

    [31]

    Liu Y, Ma J H, Fan Y, Liang Z R 2012 Phys. Med. Biol. 57 7923Google Scholar

    [32]

    Deng L Z, Mi D L, He P, Feng P, Yu P W, Chen M Y, Li Z C, Wang J, Wei B 2015 Bio-med. Mater. Eng. 26 1685Google Scholar

  • [1] 付亚鹏, 孙乾东, 李博艺, 他得安, 许凯亮. 基于RCA阵列三维超快超声血流成像方法仿真研究. 物理学报, 2023, 72(7): 074302. doi: 10.7498/aps.72.20222106
    [2] 陈松懋, 苏秀琴, 郝伟, 张振扬, 汪书潮, 朱文华, 王杰. 基于光子计数激光雷达的自适应门控抑噪及三维重建算法. 物理学报, 2022, 71(10): 104202. doi: 10.7498/aps.71.20211697
    [3] 钟虓䶮, 李卓. 原子尺度材料三维结构、磁性及动态演变的透射电子显微学表征. 物理学报, 2021, 70(6): 066801. doi: 10.7498/aps.70.20202072
    [4] 张天奎, 于明海, 董克攻, 吴玉迟, 杨靖, 陈佳, 卢峰, 李纲, 朱斌, 谭放, 王少义, 闫永宏, 谷渝秋. 激光高能X射线成像中探测器表征与电子影响研究. 物理学报, 2017, 66(24): 245201. doi: 10.7498/aps.66.245201
    [5] 代冰, 王朋, 周宇, 游承武, 胡江胜, 杨振刚, 王可嘉, 刘劲松. 小波变换在太赫兹三维成像探测内部缺陷中的应用. 物理学报, 2017, 66(8): 088701. doi: 10.7498/aps.66.088701
    [6] 潘安, 张晓菲, 王彬, 赵青, 史祎诗. 厚样品三维叠层衍射成像的实验研究. 物理学报, 2016, 65(1): 014204. doi: 10.7498/aps.65.014204
    [7] 张辉, 杨洋, 李志青. 三维a-IGZO薄膜中的电子-电子散射. 物理学报, 2016, 65(16): 167301. doi: 10.7498/aps.65.167301
    [8] 张宇, 唐志列, 吴泳波, 束刚. 基于声透镜的三维光声成像技术. 物理学报, 2015, 64(24): 240701. doi: 10.7498/aps.64.240701
    [9] 陈鹤, 于斌, 陈丹妮, 李恒, 牛憨笨. 超衍射成像中双螺旋点扩展函数的三维定位精度. 物理学报, 2013, 62(14): 144201. doi: 10.7498/aps.62.144201
    [10] 王芳, 赵星, 杨勇, 方志良, 袁小聪. 基于人眼视觉的集成成像三维显示分辨率的比较. 物理学报, 2012, 61(8): 084212. doi: 10.7498/aps.61.084212
    [11] 杨超, 刘大刚, 刘腊群, 夏蒙重, 王辉辉, 王小敏. 负氢离子源中电子能量沉积三维数值模拟研究. 物理学报, 2012, 61(15): 155205. doi: 10.7498/aps.61.155205
    [12] 周光照, 王玉丹, 任玉琦, 陈灿, 叶琳琳, 肖体乔. 相干X射线衍射成像三维重建的数字模拟研究. 物理学报, 2012, 61(1): 018701. doi: 10.7498/aps.61.018701
    [13] 刘广东, 张业荣. 乳腺癌检测的三维微波热声成像技术. 物理学报, 2011, 60(7): 074303. doi: 10.7498/aps.60.074303
    [14] 彭凯, 刘大刚, 廖臣, 刘盛纲. 三维电子回旋脉塞的数值模拟研究. 物理学报, 2011, 60(9): 091301. doi: 10.7498/aps.60.091301
    [15] 刘运全, 张 杰, 武慧春, 盛政明. 超短电子脉冲在聚焦强激光场中的三维有质动力散射. 物理学报, 2006, 55(3): 1176-1180. doi: 10.7498/aps.55.1176
    [16] 王少宏, B.Ferguson, 张存林, 张希成. Terahertz波计算机辅助三维层析成像技术. 物理学报, 2003, 52(1): 120-124. doi: 10.7498/aps.52.120
    [17] 向际鹰, 吴 震, 曾绍群, 骆清铭, 张 平, 黄德修. 弱相干扫描层析成像系统的三维传递函数分析. 物理学报, 1999, 48(10): 1831-1838. doi: 10.7498/aps.48.1831
    [18] 解文方. 三维量子点的三电子系统. 物理学报, 1997, 46(3): 563-567. doi: 10.7498/aps.46.563
    [19] 王晓, 蔡建华. 三维紧束缚电子气的等离激元理论. 物理学报, 1993, 42(7): 1149-1156. doi: 10.7498/aps.42.1149
    [20] 张世昌, 王文耀. 电磁波荡器及导引磁场中相对论电子三维运动的线性与非线性分析. 物理学报, 1991, 40(5): 748-755. doi: 10.7498/aps.40.748
  • 期刊类型引用(0)

    其他类型引用(1)

计量
  • 文章访问数:  4562
  • PDF下载量:  73
  • 被引次数: 1
出版历程
  • 收稿日期:  2021-04-12
  • 修回日期:  2021-05-06
  • 上网日期:  2021-06-07
  • 刊出日期:  2021-09-20

/

返回文章
返回