Detection Method Based on X-Ray Micro-Computed Laminography for Slab-shaped Fossils

研究背景

随着现代生物学技术在古生物学中的广泛应用，获取更多的生物学信息对古生物学研究至关重要。一些化石的生物信息，如脑颅、神经血管系统、内耳以及一些其他细微或者内部结构，在扩展我们对化石的生物性认识、以及这些结构在生物功能形态和演化研究上有很大的潜力，能为研究生命起源与演化提供关键的证据。古生物研究最初使用的经典磨片法 (Chang, 1982)，以破坏标本为代价获得其内部结构信息。后来，高精度CT技术作为一种无损检测手段在古生物学领域得到了广泛应用 (盖志琨等., 2013; 王燕芳等., 2019)，能够对包埋在岩石内部或者古生物化石内部组织结构、器官形态等进行无损观察研究。常规高精度CT通常将化石固定在一个小载物台上，将其进行360° 旋转扫描采集并重建其三维数据做进一步研究。不过，针对特殊的大型薄板状化石，想要无损的获得其内部结构，常规高精度CT有一定的局限性且存在化石损坏的风险，如在长轴方向穿透样品困难甚至无法穿透，若将化石立起来穿透短轴，又容易损坏薄板化石，再则无法将感兴趣的区域足够靠近光源以实现高分辨成像。
    近年来，在国家重大科学仪器设备开发专项支持下，中国科学院高能物理研究所(The Institute of High Energy Physics of the Chinese Academy of Sciences, IHEP) 研发的国内首台新型工业X射线显微层析成像系统160kV-Micro-CL (Liu et al., 2015, Wei et al., 2017, 尹鹏飞等., 2016)，只需将化石直接水平固定在大载物台上，通过有限角倾斜入射的方法实现高分辨三维成像｡本文主要介绍使用Micro-CL对板状化石进行三维数据采集及重建的方法，并应用到不同门类的大型板状化石研究中，该方法相较于传统高精度CT，能在保障样品安全的前提下，提高空间分辨率和图像衬度。

仪器设备

1. 新型工业X射线显微层析成像系统 (中国科学院高能物理研究所，IHEP-IVPP, 160kV-Micro-CL) ，如图1所示。
2. 225 kV-μCT (中国科学院高能物理研究所, IHEP-IVPP, 225kV-μCT) 
3. GE – Phoenix v|tome|x m成像系统 (General Electric Measurement and Control Solutions, GE, Phoenix v|tome|x m)
   
   
   图1. 160KV-Micro-CL设备图
   A. Micro-CL设备外观图； B. Micro-CL设备内部图。
   

软件

1. VGStudio 3.4 (Volume Graphics GmbH, https://www.volumegraphics.com/).

实验步骤

Micro-CL扫描板状化石的整个流程分为三个步骤：Micro-CL扫描前准备；Micro-CL三维数据采集；Micro-CL三维数据重建。整个扫描流程如图2所示，详细步骤如下：


图2. Micro-CL扫描流程图
A. Micro-CL扫描前准备流程图；B. Micro-CL三维数据采集流程图； C. Micro-CL三维数据重建流程图

一、Micro-CL扫描前准备

1. 扫描系统初始化
   打开系统电源，光机抽真空､预热､聚焦，机械运动系统回零，复位。初始位置：载物台运动到正中心(x=-150, y=-150)，探测器运动到载物台正上方倾角α=0，转角β=0。
2. 摆放固定样品
   借助激光十字定位样品扫描中心，通过升降台或者垫泡沫板等方式调整样品高度。本台设备的高度范围为0-50 cm。根据样品的大小和感兴趣区域的大小，在保证安全的情况下，调整样品高度。越靠近光源，视野越小，分辨率越高。当样品和光源的距离一定时，以载物台默认高度为测量点，记录视野大小和倾角的关系如图3所示，同一高度下，倾角越大，视野越大。
   
   
   图3. 倾角和视野大小的关系图
   
   
3. 调整扫描参数(x,y, α,r,U,I)
   首先，打开X光用DR (Digital Radiography) 实时扫描查看样品位置，分别沿X､Y轴方向移动载物台，使感兴趣区域移动到视野中心，并记录中心点即圆心的坐标 (x0,y0) 以及此时的图像。当扫描对比度较差的样品时，DR实时扫描很难直接清晰的定位感兴趣区域，此时需要借助外部标记物，常用的有铜丝､铁丝等密度大又可以贴近样品表面的物体，在正式扫描前移走。
       其次，将探测器沿C型臂移动到某一角度α(α∈[0,60°])时，再沿X轴方向移动载物台，使感兴趣区域的中心再次出现在视野中心，记录载物台在X轴方向相对运动距离，即为半径r(r∈(0,150))。若两个中心点无法重合，则适当减小倾角或者提高样品高度重新调整参数。
       然后，计算空间分辨率判断是否满足用户需求｡CT空间分辨率与放大比M、探测器像素尺寸和X射线源焦点尺寸等有关, 放大比M= SOD/SDD，SOD（source to object distance）是光源到样品的距离， SDD (source to detector distance) 是光源到探测器距离。本台设备在检测板状化石时， X射线源焦点尺寸的影响可忽略不计，空间分辨率=M*探测器像素尺寸，其中SOD= r/sin(α)，SDD=580.5mm，探测器像素尺寸=0.127mm。得到空间分辨率=r/sin(α) /580.5 *0.127 (mm) （r为旋转半径，α为探测器倾斜角度）。
       最后，根据X光穿透扫描样品后的计数值，调整扫描电压U电流I。将样品移出视野，选择叠加数目，在探测器倾斜的状态下做偏置､增益校正，采集并保存空扫文件，重建时作为背景计数使用。恢复样品位置，并将其固定在载物台上，保证在整个扫描过程中不发生位移。

二、Micro-CL三维数据采集
160kV-Micro-CL三维数据采集主要使用投影数据获取模块，通过X光DR实时扫描检查位置并微调。设置探测器参数包括模式帧频、采集时间、叠加张数，选择探测器数据校正方式包括偏置、增益、坏像素，扫描角度为360度，设置采集参数、叠加帧数、扫描数目和保存路径。开始采集投影文件，当载物台移动到设定位置时，采集数据保存并显示投影图像，然后更新扫描参数文件，载物台和探测器再移动到下一个位置，依次采集，直至扫描结束。采集到的三维数据是图像尺寸为1536*1920的序列图像集。

三、Micro-CL三维数据重建
使用160kV-Micro-CL软件打开采集到的投影文件及参数文件，设置修改重建参数，包括重建起始层、终止层、选择滤波器、调节偏置等，选择保存路径，开始重建。使用修正的FDK算法反投影重建出对应的三维数据，重建后的三维数据是图像尺寸为2048*2048的序列图像集，后期根据需求做不同的图像处理。

使用案例

一、使用Micro-CL扫描化石标本案例
早期哺乳动物化石以板状形态保存的有很多，如燕辽生物群和热河生物群中的哺乳型动物。以晚侏罗世早期燕辽生物群中的阿霍氏树贼兽 (Han et al., 2017) 为例，岩性为泥岩，平均厚度只有0.3 cm，大小尺寸为42 cm *26 cm。对于这块超薄的板状化石标本，如果使用常规高精度CT，需要将其立于载物台上，这样的载样方式对标本有着极高的风险。然而，使用Micro-CL就很安全且容易实现，将标本水平放置于50 cm*50 cm的载物台上，固定后就可以进行CL扫描。这块标本的感兴趣区域分别是该哺乳动物的整体骨架、头部以及耳区。由于样品距离光源越近，视野越小，分辨率越高。因此，用低分辨率展现阿霍氏树贼兽的整体骨架，用高分辨率放大其头部，用超高分辨率清晰精确的再现其耳区精细结构。采用低电压60 kV高电流80 μA的方式增强图像对比度，并针对不同尺寸的感兴趣区域设置不同的参数，具体参数详见表1。相关扫描结果如图4所示，从整体聚焦到精细的耳区，展示了一个全新的哺乳动物中耳类型，为研究中生代哺乳动物中耳结构演化提供了可靠的证据，也扩展了燕辽生物群的多样性。

表1. 阿霍氏树贼兽标本不同区域扫描参数表

图4. 阿霍氏树贼兽标本的Micro-CL断层图
A.整体骨架Micro-CL断层图；B.头部Micro-CL断层图；C.中耳Micro-CL断层图｡

 对于不同的化石类群，根据这些化石标本不同的岩性、大小、厚度、感兴趣区域等调整合适的参数，用Micro-CL针对感兴趣区域进行精细扫描，譬如三叠复兴鱼 (马昕莹等., 2017)、施氏慈母鸟 (Bailleul et al., 2019)、谢氏红山蜥 (Dong et al., 2019) 等脊椎动物，以及南京花 (Fu et al.,2018) 等维管植物，详细参数如表2所示，相关扫描结果如图5所示。通过大量的化石检测， Micro-CL无损获得了相关区域的内部高精度三维结构，为科研人员研究化石内部形态结构提供了新的证据，对于探索脊椎动物和维管植物演化具有重要意义｡

表2. 三叠复兴鱼､施氏慈母鸟､谢氏红山蜥､南京花标本扫描参数表

图5. 不同标本的Micro-CL扫描断层图
A. 三叠复兴鱼Micro-CL断层图；B.施氏慈母鸟Micro-CL断层图；C.谢氏红山蜥Micro-CL断层图；D.南京花的Micro-CL断层图。

二、与常规CT对比
对同样的标本，如果使用225 kV-μCT扫描，根据标本尺寸、感兴趣区域大小位置，距光源距离按照该CT相关参数计算对应的空间分辨率 (225 kV-μCT空间分辨率=　SOD/1237*0.194mm)，结果如表3所示。例如在扫描18 cm大区域时，CL可以一次扫描完成分辨率为84.5 μm，225 kV-μCT单次最大视野为9.64 cm，需要分两次扫描后期拼接，分辨率为62.73 μm；在扫描5.5 cm区域时，CL比常规CT精度略微提高；在扫描更小区域1.8 cm时，由于标本尺寸限制了旋转半径，导致物体无法足够靠近光源，使得常规CT无法获得高分辨率成像，而CL可以很靠近光源，分辨率比常规CT有较大改善。

表3. 扫描同一块化石Micro-CL和225kV-μCT空间分辨率对比表

另一方面，由于较大的板状化石长轴和短轴方向相差一般大于10倍，这就使得用常规高精度CT扫描这类样品时，X光穿过长轴时所剩无几，进一步使密度分辨降低。中科院古脊椎所毛方园在研究翔兽牙齿时，分别使用GE – Phoenix v|tome|x m和Micro-CL扫描了同一块标本的同一区域，实验结果如图6所示，A为GE常规高精度CT扫描的断层图，无法区分牙齿根部和密度相近的围岩；B为Micro-CL扫描的断层图，可以明显区分。


图6. 不同设备扫描翔兽牙齿的断层图
A. GE Phoenix v|tome|x m CT扫描断层图；B. Micro-CL扫描断层图。

讨论

大量的科研实验研究表明，显微层析成像技术(Micro-CL)可以方便安全的将样品水平固定在载物台上，解决了常规高精度CT载样需要将化石倾斜或者垂直放置可能存在的风险；Micro-CL扫描样品时，在一定范围内不受样品本身长宽尺寸限制，可以很靠近光源，获得比常规高精度CT更高的空间分辨率；同时避免了样品长轴方向计数低甚至无法穿透的问题，衬度也得到明显的提升。扫描大区域譬如18 cm阿霍氏树贼兽骨架时，Micro-CL可以一次完成；扫描大尺寸样品局部微小区域譬如1.8 cm阿霍氏树贼兽耳区时，Micro-CL可以达到高分辨率成像，清晰再现其内部的精细结构；扫描低对比度样品如翔兽牙齿时，Micro-CL可以获得更好的密度分辨，可以从密度相近的围岩中分割提取出完整的牙齿结构。

结论

针对不同的板状化石标本，根据标本的岩性､埋藏环境､厚度､尺寸等适当调整Micro-CL的扫描参数，能获取并重建化石清晰的内部形态与三维结构，为古生物学研究提供大量关键的形态学基础数据｡通过大量的实验研究表明，该方法在大型薄板状动植物化石的研究中发挥了非常重要的作用，支撑了一批创新成果的产出，这些成果对于研究脊椎动物和维管植物的起源、演化具有非常重要的意义。

致谢

感谢中国科学院古脊椎动物与古人类研究所毛方园、徐光辉、Bailleul Alida M.、董丽萍，中国科学院南京地质古生物研究所王鑫为本文提供的数据支持，感谢毛方园对本文的修改提供指导意见。本文工作得到了中国科学院战略性先导科技专项 (XDA19050102; XDB26000000)；国家自然科学基金项目 (41530102)；中国科学院前沿科学重点研究项目QYZDJ-SSW-DQC002联合资助｡

利益冲突

我声明，作者没有符合NCSP定义的竞争利益或其他可能被认为影响本文报道的结果和讨论的竞争利益。

参考文献

1. Chang, M. M. (1982). The braincase of Youngolepis, a Lower Devonian crossopterygian from Yunnan, south-western China, Stockholm University.
2. 盖志琨，朱敏. (2013). 显微CT技术在古生代鱼类研究中的应用.生命科学 25(08): 779-786.
3. 王燕芳，魏存峰，阙介民，张文定，孙翠丽，舒岩峰，侯叶茂，张久昶，戎坚，魏龙. (2019). 古生物CT装置的研制及应用. 古脊椎动物学报 57(01): 84-92.
4. Liu, B. D., Wei, Z. H., Wei, C. F., Wang, Y. F., Yuan, L. L., Shu, Y. F., Que, J. M., Sun, C. L., Wang, Y. X., Shao, Y. M., Chai, J. B. and Wei, L. (2015). An Industrial Computed Laminography Imaging System.Proceedings of Digital Industrial Radiology and Computed Tomography. (DIR2015): 22-25.
5. Wei, Z. H., Yuan, L. L., Liu, B. D., Wei, C. F., Sun, C. L., Yin, P. F., & Wei, L. (2017). A micro-CL system and its applications. Review of Scientific Instruments. 88(11): 115107
6. 尹鹏飞，侯叶茂 (2016). Micro-CL成像系统在古生物中的应用.2016年全国射线数字成像与CT新技术研讨会论文集， 2016: 79-83.
7. Han, G., Mao, F. Y., Bi, S. D., Wang, Y. and Meng, J. (2017). A Jurassic gliding euharamiyidan mammal with an ear of five auditory bones. Nature 551(7681): 451-456.
8. 马昕莹，徐光辉. (2017). 云南中三叠世(安尼期)预言鱼目(全骨鱼类:近鲱形类)一新属种，古脊椎动物学报，55(02): 162-176.
9. Bailleul, A. M., O'Connor, J., Zhang, S., Li, Z., Wang, Q., Lamanna, M. C., Zhu, X. and Zhou, Z. (2019). An Early Cretaceous enantiornithine (Aves) preserving an unlaid egg and probable medullary bone. Nat Commun 10(1): 1275.
10. Dong, L., Wang, Y., Mou, L., Zhang, G. and Evans, S. E. (2019). A new Jurassic lizard from China. Geodiversitas 41(16): 623-641, 619.
11. Fu, Q., Diez, J. B., Pole, M., Garcia Avila, M., Liu, Z. J., Chu, H., Hou, Y., Yin, P., Zhang, G. Q., Du, K. and Wang, X. (2018). An unexpected noncarpellate epigynous flower from the Jurassic of China. Elife 7: e38827.