一种基于计算机图像处理的草地盖度测算方法与流程

本发明涉及图像处理技术领域，特别涉及一种基于计算机图像处理的草地盖度测算方法。

背景技术：

传统的测算包括目估法、采样法和照相法。目估法简单易行，但是主观随意性较大，测算结果跟目估对象的实际覆盖度大小及测算人的经验密切相关。采样法是通过设定样方，将一根根样针在草地中垂直投下，记录击中草业的样针的百分比作为盖度值，操作复杂且十分耗时。照相法是数码相机普及后发展出来的方法，从数码照片中识别草叶像素，具有简单、经济、现场作业耗时少等优势，目前已经可以针对特定照片设定特定的识别规则的方式自动测算草地盖度。但是，每一张照片拍摄时的光照环境存在差异，识别规则的人为设定又带入了主观因素，造成测算结果因人而异，而且无法自适应的识别规则就不具备批处理海量照片的能力。

技术实现要素：

本发明的目的是提供一种能够实现根据垂直向下拍摄的草地照片测算出草地盖度的基于计算机图像处理的草地盖度测算方法。

为此，本发明技术方案如下：

一种基于计算机图像处理的草地盖度测算方法，包括如下步骤：

s1、利用图像采集设备采用垂直向下拍摄的方式获取能够分辨出植物与背景以及植物的边缘轮廓的草地彩色照片，然后将所述彩色照片转换为exg灰度图，并获取该灰度图中exg最小值和exg最大值；

其中，图像采集设备应保证其拍摄获得的彩色照片的分辨率能够保证从照片中清晰地分辨出叶片的内部和叶片的边缘，相机尽可能不使用短焦距的广角镜头；另外，相机曝光适度，过度则太阳光直接照射的叶片会呈亮白色，不足则叶片呈暗绿色，与背景对比差别不够明显；当直射太阳光强烈而造成遮荫和未遮荫叶片之间的形成强烈反差，可用雨伞等工具为拍摄区域遮荫以屏蔽直射太阳光，或者拍摄时开启闪光灯压制直射太阳光。

s2、依次设定exg最小值至exg最大值之间包含的n个整数作为exg阈值，将exg灰度图转换为n张表征植物和背景的黑白二值图，并计算出每张黑白二值图中植物所占百分比；

s3、逐一判断出在每个设定阈值下获取的黑白二值图中的全部边界像素，并计算出每个边界像素的exg梯度强度值，然后求取边界像素的exg梯度强度平均值作为该阈值对应的分割强度；

s4、绘制以exg阈值为横坐标、分割强度为纵坐标的exg阈值-分割强度曲线；

勾绘出exg阈值-分割强度曲线的下包络线然后通过exg阈值-分割强度曲线与包络线之间的差值运算绘制出以exg阈值为横坐标、分割强度与下包络线之间的间隙间距为纵坐标的间隙曲线，该间隙曲线上的间隙最大值对应的exg阈值即为最佳exg阈值，该最佳exg阈值对应的植物所占百分比确定为草地盖度测算结果。

进一步地，步骤s1的具体步骤为：对该彩色数码照片的全部像素进行提取，并利用过绿指数计算公式：exg＝2×g-b-r对每一个像素的过绿指数进行计算，获得该彩色数码照片的exg灰度图。

其中，由于在步骤s1中拍摄得到的彩色数码照片的图像颜色也基于rgb三原色成像原理，即自然界中所有的颜色都可以用红(r)、绿(g)、蓝(b)这三种颜色波长的不同强度组合得到；因此每张彩色数码照片中包含的所有像素都包含g、r、b三种颜色的数值；进而利用过绿指数(exg)计算公式可以将每个像素的g、r、b值带入得到对应每个像素的exg值，不同的exg值在exg灰度图表现出不同的亮度，相应地，exg值越大，像素越亮，因此通过将rgb三原色成像变换为exg值成像，就将彩色照片转换为了exg灰度图；对应每个像素均还具有一个相应的exg值。

进一步地，步骤s2的具体步骤为：依次设定exg最小值至exg最大值之间包含的n个整数为exg阈值；根据每次设定的exg阈值对exg灰度图中所有像素进行重新赋值：当像素的exg值≥exg阈值时，该像素赋值为1；当像素的exg值＜exg阈值时，该像素赋值为0；然后根据所有像素的新的赋值，将赋值为1的像素区域设为白色，将赋值为0的像素区域设为黑色，进而得到能够表征植物和背景的黑白二值图；其中，每张黑白二值图中，赋值为0的像素代表背景像素，赋值为1的像素代表植物像素，因此可进一步计算出照片中植物所占百分比，具体计算公式为：植物％＝(赋值为1的像素总数)/(赋值为1的像素总数+赋值为0的像素总数)*100％进行计算。

步骤s3中，边界像素的判断方法为：对每张黑白二值图中的全部像素逐一进行判断，当某一像素与周边8个相邻像素的赋值存在至少一个差异时，即可认定该像素为边界像素；边界像素exg梯度强度值为该边界像素周边8个相邻像素的exg最小值与最大值之差的绝对值。

步骤s4中，由于根据经验，一般草地盖度的范围在0.5％～99.5％之间，因此对绘制出的exg阈值-分割强度曲线中植物所占百分比＜0.5％和＞99.5％的曲线部分进行剔除处理。

该草地盖度测算方法基于计算机图像处理技术无需考虑其所拍摄照片时不同光照环境存在的测算结果差异，识别规则的人为设定带入的主观因素差异，是一种能够实现批量处理照片进行草地盖度测算的方法。

附图说明

图1为本发明的基于计算机图像处理的草地盖度测算方法的流程示意图；

图2为本发明的实施例1获取的数码照片；

图3为本发明的实施例1通过数码照片转换得到的exg灰度图；

图4为本发明的实施例1的exg阈值-分割强度曲线图；

图5为本发明的实施例1的间隙曲线图；

图6为本发明的实施例1的最佳exg阈值为19对应的黑白二值图；

图7为本发明的实施例2获取的数码照片；

图8为本发明的实施例2的勾勒有下包络线的exg阈值-分割强度曲线图；

图9为本发明的实施例2的最佳exg阈值对应的黑白二值图；

图10为本发明的实施例3获取的数码照片；

图11为本发明的实施例3的勾勒有下包络线的exg阈值-分割强度曲线图；

图12为本发明的实施例3的最佳exg阈值对应的黑白二值图；

图13为本发明的实施例4获取的数码照片；

图14为本发明的实施例4的勾勒有下包络线的exg阈值-分割强度曲线图；

图15为本发明的实施例4的最佳exg阈值对应的黑白二值图；

图16为本发明的实施例5获取的数码照片；

图17为本发明的实施例5的勾勒有下包络线的exg阈值-分割强度曲线图；

图18为本发明的实施例5的最佳exg阈值对应的黑白二值图；

图19为本发明的实施例6获取的数码照片；

图20为本发明的实施例6的勾勒有下包络线的exg阈值-分割强度曲线图；

图21为本发明的实施例6的最佳exg阈值对应的黑白二值图；

图22为本发明的实施例7获取的数码照片；

图23为本发明的实施例7的勾勒有下包络线的exg阈值-分割强度曲线图；

图24为本发明的实施例7的最佳exg阈值对应的黑白二值图；。

具体实施方式

下面结合附图及具体实施例对本发明做进一步的说明，但下述实施例绝非对本发明有任何限制。

实施例1

以某地区的一张垂直向下拍摄的草地照片为例，采用该基于计算机图像处理的草地盖度测算方法对该其拍摄的草地的盖度进行自动测算。

如图1所示，具体测算步骤如下：

步骤一、获取指定地区的草地彩色图像，并将其转换为exg灰度图：

手持sonydsc-w90卡片机并保持平举，在距离草地上方1.5m左右的高处垂直向下拍摄照片指定地区的草地的彩色数码照片，如图2所示；由于当天天气条件较好，光线充足，相机设置成自动档；具体地，光圈值f/8、快门速度1/1600秒、is0-400。从该长彩色数码照片可以看出其所在区域的草地属于荒漠草原，背景为沙子和砾石，植物为簇生耐干旱物种；

将该彩色数码照片扫描存储至计算机，然后利用图像软件对的全部像素进行提取，进而利用过绿指数(exg)的计算公式：exg＝2×g-b-r对每一个像素的exg进行计算，获得该彩色数码照片的exg灰度图，以及该灰度图中exg最小值和exg最大值；具体地，经图2转换得到的exg灰度图如图3所示，其中的exg最小值为-81和exg最大值为158；

步骤二、依次设定exg最小值至exg最大值之间包含的全部整数作为exg阈值，将exg灰度图转换为多张表征植物和背景的黑白二值图，并计算每张黑白二值图的植物所占百分比；

具体地，依次将从-81到158之间的240个整数作为exg阈值，根据exg阈值设定对exg灰度图中所有像素进行重新赋值。当像素的exg值≥exg阈值时，该像素赋值为1；当像素的exg值＜exg阈值时，该像素赋值为0；然后根据所有像素的新的赋值，将赋值为1的像素区域设为白色，将赋值为0的像素区域设为黑色，即得到能够表征植物和背景的黑白二值图。据此，通过依次将exg阈值设定为-81至158，最终获得240张表征植物和背景的黑白二值图；

进一步根据赋值为1的像素总数和彩色数码照片包含的像素总数，按照公式：植物％＝(赋值为1的像素总数)/(赋值为1的像素总数+赋值为0的像素总数)*100％计算出每张黑白二值图中植物所占百分比(植物％)；当然，其中“赋值为1的像素总数+赋值为0的像素总数”即为彩色照片包含的像素总数。

步骤三、逐一判断出每个设定阈值下对应的黑白二值图中的全部边界像素，并计算出每个边界像素exg梯度强度值，以求取边界像素的exg梯度强度平均值作为该阈值下的分割强度；

其中，边界像素的判断标准为：逐个像素判断该像素与周边8个相邻像素的赋值是否完全相同，若存在至少一个赋值差异，即可认定该像素为边界像素；边界像素exg梯度强度值即为该边界像素周边8个相邻像素的exg最小值与最大值之差的绝对值；边界像素的exg梯度强度平均值即为同一张黑白二值图中全部边界像素exg梯度强度值的平均值；

以此类推，依次对240张表征植物和背景的黑白二值图进行如上步骤，最终得到240个一一对应240个exg阈值的分割强度；

步骤四、如图4所示，绘制出以exg阈值为横坐标、分割强度为纵坐标的exg阈值-分割强度曲线，并剔除exg阈值-分割强度曲线中植物所占百分比＜0.5％和＞99.5％的曲线部分，以减少噪音干扰；根据步骤二计算出的各设定阈值下的植物所占百分比的值，当exg阈值等于-19时，此时的盖度测算结果为0.5％，而exg阈值等于75时，此时的盖度测算结果为99.5％；因此，对原exg阈值-分割强度曲线的两侧曲线进行剔除处理，仅保留阈值为-19至75之间exg阈值-分割强度曲线；

步骤五、依据经过剔除处理后的exg阈值-分割强度曲线勾绘出下包络线，然后通过exg阈值-分割强度曲线与包络线之间的差值运算得到如图5所示的间隙曲线；根据图5所示的间隙曲线可以获知，exg阈值-分割强度曲线与下包络线之间的间隙最大值为2.65，对应的exg阈值为39；对应从exg阈值-分割强度曲线也可以看出，以exg阈值等于39位置为界，其左侧曲线斜向上，右侧曲线接近平直，该处位置是拐点；因此，以39为阈值得到的如图5所示的表征植物和背景的黑白二值图，如图6所示，应为与彩色数码照片中的植物相符程度最高的图像，其对应的草地盖度测算结果等于15％。

实施例2

如图7～图9所示为对某地区公路边坡绿化初期类型的草地采用该基于计算机图像处理的草地盖度测算方法依次获得彩色数码照片、勾勒有下包络线的exg阈值-分割强度曲线图和最佳exg阈值对应的黑白二值图，其草地盖度测算结果为10％。

实施例3

如图10～图12所示为对某地区沙地绿化类型的草地采用该基于计算机图像处理的草地盖度测算方法依次获得彩色数码照片、勾勒有下包络线的exg阈值-分割强度曲线图和最佳exg阈值对应的黑白二值图，其草地盖度测算结果为10％。

实施例4

如图13～图15所示为对某地区庭院绿化类型的草地采用该基于计算机图像处理的草地盖度测算方法依次获得彩色数码照片、勾勒有下包络线的exg阈值-分割强度曲线图和最佳exg阈值对应的黑白二值图，其草地盖度测算结果为10％。

实施例5

如图16～图18所示为对某地区草坪绿化类型的草地采用该基于计算机图像处理的草地盖度测算方法依次获得彩色数码照片、勾勒有下包络线的exg阈值-分割强度曲线图和最佳exg阈值对应的黑白二值图，其草地盖度测算结果为10％。

实施例6

如图19～图21所示为对某地区三叶草草坪类型的草地采用该基于计算机图像处理的草地盖度测算方法依次获得彩色数码照片、勾勒有下包络线的exg阈值-分割强度曲线图和最佳exg阈值对应的黑白二值图，其草地盖度测算结果为10％。

实施例7

如图22～图24所示为对某地区三叶草草坪类型的草地采用该基于计算机图像处理的草地盖度测算方法依次获得彩色数码照片、勾勒有下包络线的exg阈值-分割强度曲线图和最佳exg阈值对应的黑白二值图，其草地盖度测算结果为10％。

综上所述，该基于计算机图像处理的草地盖度测算方法适用于各类草地的盖度测算，其测算过程中不存在人为主观判断，且不存在对不同条件下相同地点拍摄的彩色照片处理结果的差异，因此相对于现有的测算方法准确率更高，误差更小。