本发明专利技术公开了一种基于CT扫描技术量化土壤团聚体孔隙空间结构方法，包括首先通过土壤筛选取不同粒级的土壤团聚体样品；将样品放置于CT扫描设备进行扫描；界定团聚体边界；划分出团聚体中心部位和边缘部位；选取孔隙阈值；通过Interactive Thresholding工具，将选取的孔隙阈值，应用于掩膜Mask1、Mask2和Mask3下的团聚体样品，获得团聚体及其中心和边缘部位孔隙的二值图；通过Avizo软件的体积渲染工具进行图像三维可视化；再通过Avizo软件的计算模块得到三维孔隙结构特征参数。本发明专利技术通过将CT扫描仪器与专业三维分析软件相结合，实现对团聚体边界的提取以及对团聚体中心部位和边界部位的划分，采用最合适的阈值分割方法，实现对土壤团聚体孔隙空间结构的三维量化分析。对土壤团聚体孔隙空间结构的三维量化分析。对土壤团聚体孔隙空间结构的三维量化分析。
【技术实现步骤摘要】
一种基于CT扫描技术量化土壤团聚体孔隙空间结构方法
[0001]本专利技术涉及农业工程
，具体来说，涉及一种基于CT扫描技术量化土壤团聚体孔隙空间结构方法。
技术介绍
[0002]迄今为止，在团聚体尺度上对土壤孔隙度进行完整的测量是实验无法达到的。传统的技术，如树脂浸渍土块的薄切片非常耗时，因此只能提供有限的(相当厚的)截面。x射线计算机断层(CT)已被提出作为一种可行的方法来测量完整的土壤孔隙系统，它是快速和无损的，利用CT图像可以直接观察孔隙度和孔隙形态。目前，国际上一些科研工作者已经将CT扫描法作为土壤空间网络结构研究的前沿方法，但是由于扫描成本，这一方法并未广泛应用，在国内土壤学中的应用更是少之又少。在国际上，虽然这一方法已经得到应用，但是团聚体边界的确定以及内外部划分的方面尚存未解决问题。另外，在国内虽然有学者开始利用CT扫描研究土壤大孔隙结构，但主要是对土柱尺度的分析，很少对土壤团聚体内部结构进行三维量化分析；在分析软件方面主要借助于CT扫描仪的自带软件或者Image J等开源软件，Avizo等专业软件的使用较少。
技术实现思路
[0003]针对相关技术中的上述技术问题，本专利技术提出一种基于CT扫描技术量化土壤团聚体孔隙空间结构方法，能够克服现有技术方法的上述不足。
[0004]为实现上述技术目的，本专利技术的技术方案是这样实现的：一种基于CT扫描技术量化土壤团聚体孔隙空间结构方法，包括以下步骤：S1: 首先选取样品，通过土壤筛选取不同粒级选取土壤团聚体；S2: 将选取好的样品放置于CT扫描设备进行扫描；S3: 扫描完成后，界定团聚体边界；S4: 然后划分出团聚体中心部位和边缘部位；S5: 选取孔隙阈值；S6: 通过工具(Interactive Thresholding)，将选取的孔隙阈值，应用于掩膜Mask1、Mask2和Mask3下的团聚体样品，获得团聚体及其中心和边缘部位孔隙的二值图；S7: 通过Avizo软件的体积渲染工具，实现团聚体及其中心和边缘部位孔隙的三维可视化；S8: 再通过Avizo软件的计算模块得到三维孔隙结构特征参数。
[0005]进一步地，步骤S3中，界定边界包括以下步骤：S31：将扫描好的团聚体样品的图片导入Avizo软件中；S32：然后用Volume Edit工具将团聚体外的杂点去掉；S33：在segmentation模块进行掩膜提取，使用Avizo软件中的魔棒工具选中土壤基质，再用Fill工具将其中的空隙填满，获得团聚体掩膜Mask1/土壤边界。
[0006]进一步地，步骤S4中，划分出团聚体中心部位和边界部位包括以下步骤：S41：使用shrink工具，将团聚体边界向内收缩一部分体素；S42：再通过膨胀收缩3
‑
5个体素，得到团聚体的中心部位及其掩膜Mask2；S43：将步骤S33获得的掩膜Mask1去除步骤S42获取的掩膜Mask2，即得到团聚体边界部位的掩膜Mask3；S44：在掩膜Mask2的基础上，通过渲染步骤S33后的团聚体图像，即得到整个团聚体三维图像。
[0007]S45：在掩膜Mask2的基础上，通过渲染步骤S32后的团聚体图像，即得到团聚体中心部位三维图像。
[0008]S46：在掩膜Mask3的基础上，通过渲染步骤S32后的团聚体图像，即得到团聚体边缘部位三维图像。
[0009]进一步地，步骤S5中，选取孔隙阈值包括以下步骤：S51：使用Histogram工具，得到步骤S31后的团聚体图像的灰度直方图，途中有明显的双峰，分别代表土壤基质和孔隙部分；S52：基本确定出Min部分为孔隙，Max部分为土壤基质，Mix部分为阈值混合区；S53：针对阈值混合区Mix，提取这部分阈值的图像，通过大律算法，进一步计算孔隙Min
’
和土壤基质Max
’
的阈值划分；S54：最终团聚体孔隙阈值为0
‑ꢀ
Min
’
。
[0010]进一步地，步骤S6中，孔隙阈值分割包括以下步骤：S61：在掩膜Mask1的基础上，通过工具(Interactive Thresholding)，将S44后团聚体进行孔隙阈值分割，利用S54得到的孔隙阈值0
‑
Min
’
，获得团聚体孔隙的二值图；S62：在掩膜Mask2的基础上，通过工具(Interactive Thresholding)，将S45后团聚体进行孔隙阈值分割，利用S54得到的孔隙阈值0
‑
Min
’
，获得团聚体中心部位孔隙的二值图；S63：在掩膜Mask3的基础上，通过工具(Interactive Thresholding)，将S46后团聚体进行孔隙阈值分割，利用S54得到的孔隙阈值0
‑
Min
’
，获得团聚体边缘部位孔隙的二值图。
[0011]进一步地，三维孔隙结构特征参数包括孔隙的体积、孔隙数量、孔隙等效直径、孔隙表面积和孔隙长度。
[0012]进一步地，将三维孔隙结构特征参数中的孔隙度、孔隙数量密度、孔隙表面积密度和孔隙长度密度分别定义为孔隙体积、孔隙数量、孔隙总表面积和孔隙总长度与Mask1体积的比值。
[0013]本专利技术的有益效果：通过将CT扫描仪器与专业三维分析软件相结合，实现对团聚体边界的提取以及对团聚体中心部位和边界部位的划分，采用最合适的阈值分割方法，实现对土壤团聚体孔隙空间结构的三维量化分析。
附图说明
[0014]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施
例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0015]图1是根据本专利技术实施例所述的基于CT扫描技术量化土壤团聚体孔隙空间结构方法的土壤团聚体解译流程框图。
[0016]图2是根据本专利技术实施例所述的基于CT扫描技术量化土壤团聚体孔隙空间结构方法的土壤团聚体内
‑
外部孔隙解译流程框图。
[0017]图3是根据本专利技术实施例所述的基于CT扫描技术量化土壤团聚体孔隙空间结构方法的利用灰度直方图选择土壤孔隙阈值示意图。
具体实施方式
[0018]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本专利技术保护的范围，为了方便理解本专利技术的上述技术方案，以下通过具体使用方式上对本专利技术的上述技术方案进行详细说明。
[0019]如图1
‑
2所示，根据本专利技术实施例所述的基于CT扫描技术量化土壤团聚体孔隙空间结构方法，步骤包括：首先选取样品A，通过土壤筛划分不同粒级选取土壤团聚体样品A，通常选取外形比较规则、砾石含本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种基于CT扫描技术量化土壤团聚体孔隙空间结构方法，其特征在于，包括以下步骤：S1: 首先选取样品，通过土壤筛选取不同粒级选取土壤团聚体；S2: 将选取好的样品放置于CT扫描设备进行扫描；S3: 扫描完成后，界定团聚体边界；S4: 然后划分出团聚体中心部位和边缘部位；S5: 选取孔隙阈值；S6: 通过工具(Interactive Thresholding)，将选取的孔隙阈值，应用于掩膜Mask1、Mask2和Mask3下的团聚体样品，获得团聚体及其中心和边缘部位孔隙的二值图；S7: 通过Avizo软件的体积渲染工具，实现团聚体及其中心和边缘部位孔隙的三维可视化；S8: 再通过Avizo软件的计算模块得到三维孔隙结构特征参数。2.根据权利要求1所述的基于CT扫描技术量化土壤团聚体孔隙空间结构方法，其特征在于，步骤S3中，界定边界包括以下步骤：S31：将扫描好的团聚体样品的图片导入Avizo软件中；S32：然后用Volume Edit工具将团聚体外的杂点去掉；S33：在segmentation模块进行掩膜提取，使用Avizo软件中的魔棒工具选中土壤基质，再用Fill工具将其中的空隙填满，获得团聚体掩膜Mask1/土壤边界。3.根据权利要求1所述的基于CT扫描技术量化土壤团聚体孔隙空间结构方法，其特征在于，步骤S4中，划分出团聚体中心部位和边界部位包括以下步骤：S41：使用shrink工具，将团聚体边界向内收缩一部分体素；S42：再通过膨胀收缩3
‑
5个体素，得到团聚体的中心部位及其掩膜Mask2；S43：将步骤S33获得的掩膜Mask1去除步骤S42获取的掩膜Mask2，即得到团聚体边界部位的掩膜Mask3；S44：在掩膜Mask2的基础上，通过渲染步骤S33后的团聚体图像，即得到整个团聚体三维图像。4.S45：在掩膜Mask2的基础上，通过渲染步骤S32后的团聚体图像，即得到团聚体中心部位三维图像。5.S46：在掩膜Mask3的基础上，通过渲染步骤S32后的团聚体图像，即得到团聚体边缘部位三维图像。6.根据权利要求1所述的基于CT扫描技术量化土壤团聚体孔隙空间结构方法，其特征在...
【专利技术属性】
技术研发人员：胡霞，赵云朵，李小雁，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。