从在土X射线计算机断层扫描数据量和表面分析:提取指标三维根系

摘要

一种用于获得从土壤获得的X射线计算机断层扫描数据的视觉和定量根部结构信息的方法，提出。

摘要

植物根系在植物 - 土壤 - 微生物相互作用发生在根际的重要作用，以及与处理气候变化和作物管理具有重要意义。在他们自然环境对根的定量大小的信息是非常宝贵的学习根系生长，在涉及植物环境过程。 X射线计算机断层扫描（XCT）已被证明是用于原位根扫描和分析的有效工具。我们的目的是开发一种没有成本和高效率的工具，从三维（3D）断层扫描数据，而不管其形状接近根部的表面和体积。一个草原dropseed的根结构（ 鼠尾heterolepis）标本用XCT成像。根被重建，并且使用许可和开源软件的组合从数据萃取主根结构。等值面多边形网格然后为便于分析创建。我们已经开发Ť他独立的应用程序imeshJ，在MATLAB ^1中产生，以计算从网格根体积和表面积。 imeshJ的输出是表面积（以mm ^2）和体积（毫米^3）。的过程中，利用从成像定量根分析的工具的独特的组合，进行说明。证明了XCT的结合，开源软件是一个强大的组合，以非侵入性的图像植物根系样品，段根数据，并从中提取三维数据的定量信息。处理三维数据的这种方法应适用于在有相似的X射线衰减的部件和困难之间具有分割产生连接的其他材料/样品的系统。

引言

根，作为根际^2-5的一部分，代表植物生物学的"看不见"的一部分，因为土壤使得难以图像根非侵入^6,7，但是，研究土壤环境中的根生长和相互作用是理解临界根/植物的生长和养分循环，进而影响造林，粮食安全和气候。 X射线计算机断层扫描（XCT）已被证明是在其本地环境⁸植物根样本的非侵入性成像的有价值的工具。为了测量在不同条件下的根发育和尺寸的变化，并且可以比较从不同的数据集/标本数据，需要提取该断层摄影数据的定量信息。从周围的土壤的根数据的分割，即，基础图像的从一切它周围其他的隔离（包括，例如，一个相邻植物）是ACCU前一个关键步骤率大小分析可以做到的。然而，一个简单的阈值的方法是经常不可行的根数据。在土壤成像植物根部有关的挑战包括在根材料的X射线衰减特性的变化，和引起的水和有机物质的根和土壤之间的衰减值的重叠。这些问题已被寄望Mairhofer 等近期解决。在他们的视觉跟踪工具RooTrak ^7,9。成功的分割后的下一步是根体积和表面积的准确测定。体积可以通过计数的体素的数量并作为⁷之前所示的体素"尺寸的立方相乘来估计。为了更准确地测定根表面积和体积的，所述分段根系的等值面可以由三角形网格使用被称为移动立方体¹⁰的算法表示。开源的ImageJ ¹¹可以采用近似第基于对移动立方体算法E根卷。据我们所知，致力于在厘米范围计算为根基于标本断层体积/表面数据只有开放源代码软件数量有限，以上是目前可用^12。我们看了¹³张开源软件侧重于根系生长，旨在细胞特征在单细胞分辨率，可以定量分析的体积。一些专门为整根系统¹⁴开源软件是极好的基础上近似，他们的造型实际上是管小口径管的根部系统。然而，随着二维图像一些工作，无法处理3D堆栈^14。此外，当表面粗糙和不均匀的形状，如那些树根系统，进行了研究，管状形状近似可能是无效的。另一种方法¹⁵使用创新绕过第二维（2D）旋转图像序列Ë需要昂贵的CT扫描仪。它可以测量，记录和显示根系长度。我们从那些仅可以买到^16-18测试的软件;一个不出现，以便能够处理三维图像栈^16，第二是一个叶面积和根长度测量工具^17，而第三个是基于颜色分析^18。在此基础上的调查，我们建议，从3D断层扫描数据，而不管其形状接近根部的表面和体积的无成本的选择是理想的。

上可免费获得的RooTrak和ImageJ的基础上，我们已经开发出一种方案，名为imeshJ（参见补充编码文件），其处理来自分段根数据产生的等值面目（表面立体文件），并通过计算该根的体积和表面积做网格上的三角形索引数据简单的几何计算。在这里，我们报告说，综合运用XCT成像的方法，数据重建和可视化（软件CT临3D和VG工作室），从3D数据（开源软件ImageJ的和RooTrak）土壤样品的根分割，并且表面和体积信息提取从三角形网格（ImageJ的和计算机代码imeshJ）。

方案

注意:X射线断层摄影扫描仪的操作需要两个一般辐射训练，和特定仪器的辐射安全训练。相关实验者的实验室所有相应的程序应遵循。

1.根成像

注意:此步骤描述了在管状塑料盆原文土壤保持的草标本的摄像（塑料管，其直径为40毫米，210毫米的高度，以及约2mm壁厚）。

1. 放在仪器的样品操纵器盆栽植物在期望目标倍率的距离。在一个2英寸直径的支架的植物中，样品到源的距离应为约3英寸（7厘米）。
2. 调整所述X射线扫描设置，以达到最佳的颜色（灰度级）的检测器图像的对比度。注意:这些设置中使用的仪器控制软件可用。
   1. 一套X射线电源设置; 85千伏和190μA在这次考试中使用PLE。
   2. 设置曝光时间。这里，被用于更好信噪比1秒的相对长的曝光时间。
   3. 集数每个投影预测和帧;每个投影4帧，共计3,142预测建议为良好的数据统计。
   4. 使用运行通过选择"阴影校正"选项卡，并单击"创建"上面设置的测量条件阴影校正。
      注意:当具有恒定集合X射线通量的照亮变化在成像装置的像素的响应的阴影校正补偿。这个过程需要空白图像（与光束路径移除样品）与X射线束接通，并且与梁关闭。这种校正被施加到收集到的所有图像。
   5. 选择"最小化环假象"选项（也称为"穿梭模式"）;样品将在角的步骤，而在投影图像获取旋转。这导致的DAT一个采集以较慢的速度，但有助于消除环假象。
   6. 通过点击采集标签下的"采集"按钮开始扫描（与上面列出的设置，图像采集将需要大约4小时）。

2.数据重构

注意:本节介绍了从原始图像（从CT扫描X光片）三维数据的重建。

1. 加载的原始数据到程序中。
2. 比较第一和最后一个图像，以确保样品不移动或扫描设置数据采集过程中没有改变（作为最后图象被检体的^360°旋转后取它们应该是几乎相同）。
3. 通过选择"旋转中心"选项卡，并单击"开始"计算旋转（COR）的中心;使用选项"自动"COR寻找具有"高品质"的精度，而"双反"（上，下）片小号选举COR计算。
4. 选择样本量进行重构:选择"音量"选项卡，并编辑使用缩略图音量选择窗口。
5. 执行重建通过单击"开始"，以创建一个包含3D数据的卷文件。

3.数据处理/分割

注:这部分描述了要采取的步骤，以在程序RooTrak跟踪根准备用于进一步处理重建的数据，因为他们从任何周围材料分支通过土壤，并隔离根部产生的仅仅是二进制图像的堆叠根本身。

1. 在ImageJ的数据量的处理准备一个RooTrak处理的图像堆栈:
   1. 加载卷文件到ImageJ的。
   2. 通过调整亮度和对比度设置优化根系与土壤之间的图像对比（点击图像/调整/亮度/对比度）。当的图像内感兴趣的区域是可见的和clearly区分，设置将被视为进行了优化。
   3. 保存为JPEG，BMP，或PNG格式的图像堆栈。
2. 处理RooTrak分割的根源:
   1. 加载图像进栈RooTrak（进入"工具"选项卡，然后按"跟踪"）。
   2. 内根设置种子点:点击内的各相关根段卷中的数据的顶视图切片可见的几个点。
   3. 设置跟踪参数"平滑度"和"相似性"分别为0.3和0.8。
   4. 运行跟踪功能。这会跟随根从顶部图像切片一路底部片。
   5. 查看体数据后，选择根据该可用数据卷片的数目;在这种情况下，跟踪被停止了在200切片，相当于6.2毫米的深度，其中，根的边界变得不明确（被跟踪的根的图像开始混合到该土壤）。
      注意:产生的图像堆栈将被自动保存到哪里输出目录已创建。

4.音量和表面分析

注:此步骤描述了从RooTrak创建的图像堆栈等值面网格生成。

1. 转换RooTrak图像堆栈成ImageJ的二进制图像格式。选择"进程"，然后是"二进制"，然后"让二进制"。
2. 使用开放源代码的ImageJ插件，BoneJ，创造了三角形网格;在ImageJ的选择"插件"，然后在"BoneJ"，然后"等值面"。
3. 设置"重采样"和"阈值"为6，120，分别为（默认设置）。勾选"显示面"，然后按"确定"按钮。
4. 在"3D查看器"，单击文件选项卡，然后在"导出面"，然后保存为"STL（二进制）"。
5. 打开imeshJ，选择STL文件，并在微米像素输入大小。点击"Calculate表面积"在²毫米收购总样本根表面积。同样，点击"计算容积"，取得³毫米的总样本根卷。

结果

由两个试样茎天然草大草原dropseed（ 鼠尾heterolepis）和原土围绕它是从一个住宅区取和放置在如图1所示的小管形夹持器的该重构的数据的体素尺寸为约31微米×31微米×31微米。重构的卷文件被用于使用开源图像处理程序的ImageJ 1.6 ¹¹从一个选定的方向（顶视图）创建图像的堆叠。该卷的数据也亮堂这个计划，以增加根系与土壤值之间的​​对比。从重建的数据，很明显，根，土壤，最有可能的有机物的部分组件，具有导致小到图像中没 ​​有灰度对比度非常相似的X射线衰减系数（ 图2）。

- [RooTrak，用于分割方案，是在该中心在诺丁汉⁷大学开发为一体的植物生物学的一个开源项目。它是专门设计，因为他们分支通过土壤，并从周围的材料，以产生二进制图像的堆叠隔离根部跟踪根。 RooTrak已经显示，以产生分割比简单的阈值...

讨论

X射线计算机断层扫描和证明一些开放源代码程序的组合是一个强大的组合，以非侵入性的图像植物根系样品，段根数据，并从中提取三维数据的定量信息（表面积和体积）。我们的可视化和测量功能的能力总是由扫描分辨率，以及由RooTrak软件的限制的限制。然而，扫描分辨率足以捕捉大多数的本研究中的样品的特性，并RooTrak能够成功段根的显​​著一部分。在这项工作中所用的RooTrak的版本没有跟踪向上行进根段（向上分支不存在于一个显著程度在试样这里所研究）;该方案的更新版本解决了这个问题^9。

正如引言中暗示，该文献的调查表明，一个无成本的软件选件来计算根volum从厘米级及以上标本的3D断层扫描数据E /表面积是非常可取的。这里所报告的方法的意义在于，根被分割后，它的工作原理与由广泛使用的程序的ImageJ创建的3D数据格式。分析的心脏是从其图像根由三角网格近似的体积和表面积的计算。在.STL格式由ImageJ的生成的三角网格，由构成该表面近似三角形的指数来表示。imeshJ提取从二进制文件格式指数与能够计算从指数的表面积（这是矢量之间的接合的叉积的简单幅...

材料

本文使用的材料清单

姓名	公司	目录编号	评论
X-Tek/Metris XTH 320/225 kV	尼康计量		X 射线断层扫描仪
Inspect X	尼康计量		仪器控制软件
CT Pro 3D	尼康计量		重建软件，版本 XT 2.2
VG Studio MAX	Visual Graphics GmbH		用于 3D 体积的可视化软件，版本 2.1.5
ImageJ	开源		图像处理和分析软件，版本 1.6
RooTrak	开源		根分割软件，版本 0.3.1-b1 beta
imeshJ	EMSL		MATLAB 脚本由作者开发
草原水滴草样本			从住宅区地面获得的样本