MATLAB数字图像处理在根系固土中的应用研究

阎梦晴 于司杭 王培亮 王逸飞 李明宝

摘要 在水土防治中采用植物根系固定土壤的措施，引入基于MATLAB的数字图像处理技术，为植物根系固土机制的研究提供根系和土壤等的形态特征参数。对比Sobel、Prewitt、Roberts、LoG、Canny 5种边缘检测器的效果，对比基于梯度的边缘信息改进全局阈值处理方法和Otsu方法进行最佳全局阈值处理的效果，确定Canny边缘检测器和Otsu阈值处理方法相结合的最优算法。该数字图像处理技术可为植物根系固土机制的研究提供可靠的分析手段。

关键词 根系固土；数字图像处理；Canny边缘检测器；Otsu阈值处理

中图分类号 S157.1 文献标识码 A 文章编号 0517-6611（2014）31-11039-04

Study of Soil Reinforcement by Plant Roots Using MATLAB Digital Image Processing

YAN Mengqing， YU Sihang， WANG Peiliang*， LI Mingbao* et al

（School of Civil Engineering， Northeast Forestry University， Harbin， Heilongjiang 150040）

Abstract For soil conservation， adopting the root reinforcement measure and using digital image processing based on MATLAB can provide morphological character parameters of plant roots and soil for the study of mechanism of soil reinforcement by plant roots. After constrasting the effects of five edge detectors including Sobel edge detector， Prewitt edge detector， Roberts edge detector， LoG edge detector and Canny edge detector， and the effects of two methods of threshold processing， including overall threshold processing which improved the edge information and based on gradient and Otsu threshold processing， the Canny edge detector and the Otsu threshold processing are the best algorithm. This Digital Image processing based on MATLAB provides a dependable method in the study of mechanism of soil reinforcement by plant roots.

Key words Soil reinforcement by plant roots； Digital image processing； Canny edge detector； Otsu threshold processing

我國是世界上水土流失最严重的国家之一^[1]。水土流失会破坏土壤的肥力，造成土壤干层，淤积水库、阻塞河道、抬高河床，恶化生态环境。虽然自新中国成立以来，我国的水土流失治理已经逐渐成为国建的生态建设重点工程，但是由于我国国土面积大，自然地理条件特殊，以及长期以来对水土资源的过度利用，我国水土流失仍然面积大、分布广，防治任务艰巨^[2]。防止水土流失，采用科学、环保的植物固土措施，引入数字图像处理技术研究植物根系固土机制，能够为我国农业水土保持增添新的技术手段，促进我国生态文明建设的发展。

水土流失的防治措施有工程措施和植物措施，其中植物措施不仅可以防治水土流失，还可改善局部的生态环境，节约经济成本，具有生态可持续性。近年来，由于生态环境建设的需要及各学科科学技术的发展，植物根系固定土体机理的研究成为根系研究的焦点^[3]。根系与土壤形成复合结构体，起到固定植物地上部分以及固定土壤、防止水土流失的作用。对植物根系固定土体的机理进行综合定量研究，对于搭建优良的人工水土保持生态系统不仅具有广阔的应用前景，而且具有重要的学术意义^[1]。对于根系固土力学机理的研究，已经得到越来越多学者的重视，植物根系防治水土流失已经成为水土保持领域的一个研究热点。

近年来，国内外科研专家运用各种科学技术方法对根系固定土体的力学机理进行了研究，测定植物根系本身的抗拉力^[4-5]，采用原位剪切试验得到的有根土和无根土的抗挤强度值作为衡量根系固土性质的指标^[6]。随着科学技术手段的不断进步，土壤中根系的几何形态和分布情况能够通过各种试验方法得到，在研究植物根系在土壤中的分布以及根系几何形态的提取与仿真方面，也不可避免地需要对根土断层图像、土壤裂缝图像进行数字图像处理，为科研人员提供根系的形态参数数据。

运用数字图象处理技术，对根系固土抗裂性试验中的试验图像进行特征提取，将裂缝的图像特征（如裂缝宽度、深度等）用以衡量根系固土的抗裂有效性，补充完善了宏观的试验指标，通过这些指标可以直观地反映根系固土的程度。通过对图像特征的提取，建立更加全面的参数，这些参数可作为研究抗裂性及其影响因素之间的桥梁。运用数字图像处理，能够使根系固土抗裂性这一性质转化为更加具体且全面的数学指标，通过这些量化的指标可建立抗裂性与其影响因素之间的关系，探究它们的影响机理，这将为研究植物根系固土防治水土流失奠定很好的理论基础。数字图像处理技术的引入，将成为水土保持研究领域的一个重要环节，将为分析根系固土机制提供有效的技术手段。

1 基于MATLAB的数字图像处理

数字图像处理，把一幅图像定义为一个二维函数f（x，y），其中x、y是空间坐标，f是任意坐标（x，y）处的幅度（称为亮度或灰度）。MATLAB（矩阵实验室）是处理矩阵和线性代数的强大的数值计算软件，将图像作为矩阵来处理，选用MATLAB软件能有效地提高图像处理的效率，准确地提取出数字图像的特征。将坐标值（x，y）数字化称为采样，将幅值8数字化称为量化，采样和量化后得到一个实数矩阵。假设，对图像f（x，y）采样后得到M×N的图像，坐标值为离散量。将基于MATLAB的数字图像处理技术应用于根系固土抗裂性试验，旨在得到土壤裂缝的图像特征。

1.1 边缘检测算法

从输入的图像中提取属性，进行图像分割。分割算法的基本原理是基于灰度的突变来分割图像，基于检测灰度值的不连续来检测边缘。灰度的不连续性的检测采用一阶导数和二阶导数。图像处理中的二阶导数由于对噪声过于敏感而很少直接用于边缘检测，仅在与其他边缘检测技术结合时使用，在此不作讨论。

图像二维函数f（x，y）的梯度为

1.2.5 Canny边缘检测器。

与LoG检测器相比，Canny边缘检测器指定了高斯滤波器中的标准差σ，以此来平滑图像，减少噪声。在（x，y）处计算局部梯度[g2x+g2y]1/2和边缘方向arctan（gx/gy），通过寻找f（x，y）的梯度局部最大发现边缘。Canny边缘检测器使用了两个阈值来检测强边缘和弱边缘，此方法更有可能检测到真正的弱边缘。

2 结果与分析

为得到土壤裂缝的形态特征参数，要选用适宜的边缘检测器以及阈值处理方法。通过5种不同边缘检测器的处理结果图与原裂缝图的对比，选择出精度最高的边缘检测器。在此基础上，采用两种不同的阈值处理方法，对比所得图像的精确度，确定最适宜的阈值处理方法。

2.1 边缘检测器的选用

土壤的裂缝图像来自于牛国权^[8]所做的根系固土抗裂有效性试验，见图2（a）。

图2为5种不同边缘检测器处理土壤裂缝图像后的结果图。由图2（b）、（c）可知，Prewitt边缘检测器和Sobel边缘检测器的效果均不佳，二者显然无法检测到图像中的弱边缘，无法有效分割出土壤裂缝的图像。Roberts边缘检测器原理如图1所示，与Prewitt和Sobel边缘检测器相比，Roberts边缘检测器模板是非对称的，虽然简单，但用于检测土壤裂缝的图像特征时，如图2（d）所示，其效果不好。由图2（e）可知，LoG检测器结果比前3种检测器的效果好，但不及Canny边缘检测器。由图2（f）可知，Canny边缘检测器功能强大，能检测到细微的图像特征，但是在检测土壤裂缝时，还需调整参数尽可能地减少不相关的细节。

土壤裂缝图像不同于其他图像，为了获取裂缝的形态参数，通过对比研究了5种边缘检测器以及2种阈值处理方法，确定了最优边缘检测器和阈值处理方法，即结合Canny边缘检测器和Otsu阈值处理的方法，能够得到最精确的裂缝图像。图6所示的图像能够精确地反映原图中土壤的裂缝形态。将二维图像数字化，可得到裂缝边缘的二值化图像，在MATLAB中以矩阵形式存储了图像像素，再后期运用裂缝形态特征参数进行根系固土机制研究时，即可直接调用。

3 结论

该研究以根系固土试验中土壤裂缝图像为例，运用MATLAB软件，对图像进行数字图像处理，对比Sobel、Prewitt、Roberts、LoG、Canny这5种边缘检测器的效果，选择最适合用于检测土壤裂缝的Canny边缘检测器。比较基于梯度的边缘信息改进全局阈值处理方法和Otsu方法进行最佳全局阈值处理的效果，确定使用Otsu方法进行阈值处理能得到最优阈值。结合Canny边缘检测器和Otsu阈值处理方法，是最适用于土壤裂缝图像的数字图像处理方法。

为防止水土流失，对植物根系固定土体机制的研究中，引入基于MATLAB的数字图像处理方法，能够得到宏观图像中如根系、土壤裂缝等图像的形态特征参数。处理后的根系、裂缝等边缘像素的特征以及根系图像像素尺寸与实际尺寸存在线性关系，根据各个参数的几何特性可分析出根系的长度、表面积、平均直径、体积以及根系间的夹角，或是裂缝的长度、宽度、面积等。通过该研究结果，可以为关于根系固土相关领域的研究人员提供可靠分析手段以及根系的形态参数数据。

参考文献

[1] 刘震.水土保持 60 年：成就·经验·发展对策[J].中国水土保持科学，2009（4）：1-6.

[2] 田卫堂，胡维银，李军，等.我国水土流失现状和防治对策分析[J].水土保持研究，2008，15（4）：204-209.

[3] 陈丽华，余新晓，宋维峰，等.林木根系固土力学机制[M].北京：科学出版社，2008：87-115.

[4] 张瑜.微根窗根系的图像处理方法研究[D].哈尔滨：东北林业大学，2012.

[5]

史敏华，王棣，李任敏.石灰岩区主要水保灌木根系分布特征与根抗拉力研究初报[J].山西林业科技，1994（1）：17-19.

[6] 赵丽兵，张宝贵，苏志珠.草本植物根系增强土壤抗剪切强度的量化研究[J].中国生态农业学报，2008，16（3）：718-722.

[7] GONZALEZ R C，WOODS R E.Digital Image Processing[M].3rd ed.Prentice Hall，2007.

[8] 牛國权.基于力学特性的根系固土抗裂性有效性研究[D].呼和浩特：内蒙古农业大学，2010.