CN103983217B - 一种叶面积指数测量方法及系统 - Google Patents

一种叶面积指数测量方法及系统 Download PDF

Info

Publication number
CN103983217B
CN103983217B CN201410229812.3A CN201410229812A CN103983217B CN 103983217 B CN103983217 B CN 103983217B CN 201410229812 A CN201410229812 A CN 201410229812A CN 103983217 B CN103983217 B CN 103983217B
Authority
CN
China
Prior art keywords
blade
area
image
sample prescription
remaining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201410229812.3A
Other languages
English (en)
Other versions
CN103983217A (zh
Inventor
张兵
焦全军
刘良云
胡勇
刘新杰
唐凤莉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Remote Sensing and Digital Earth of CAS
Original Assignee
Institute of Remote Sensing and Digital Earth of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Remote Sensing and Digital Earth of CAS filed Critical Institute of Remote Sensing and Digital Earth of CAS
Priority to CN201410229812.3A priority Critical patent/CN103983217B/zh
Publication of CN103983217A publication Critical patent/CN103983217A/zh
Application granted granted Critical
Publication of CN103983217B publication Critical patent/CN103983217B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供一种叶面积指数测量方法及系统,其中叶面积指数测量方法,包括:获取样方内采集叶片的图像;以采集叶片的图像为依据,获取样方内所采集的全部叶片的面积;获取样方内残余叶片的图像;基于残余叶片的图像,获取样方内残余叶片的面积;将采集的全部叶片的面积和残余叶片的面积加和,得到样方内的叶片总面积A,这样将采集的全部叶片的面积和残余叶片的面积加和后得到的叶片总面积基本等于样方内所有叶片的面积,再依据公式:LAI=A/S得到样方的叶面积指数LAI的准确度相对提高,其中S为样方面积。

Description

一种叶面积指数测量方法及系统
技术领域
本发明涉及生态系统技术领域,特别涉及一种叶面积指数测量方法及系统。
背景技术
叶面积指数(LeafAreaIdex,LAI)又叫叶面积系数,是指单位土地面积上植物叶片总面积占土地面积的倍数,即叶面积指数=叶片总面积/土地面积是指一块地上植物叶片的总面积与占地面积的比值。在生态系统中,叶面积指数是生态系统的一个重要结构参数,用来反映植物叶面数量、冠层结构变化、植物群落生命活力及其环境效应,为植物冠层表面物质和能量交换的描述提供结构化的定量信息,并在生态系统碳积累、植被生产力和土壤、植物、大气间相互作用的能量平衡、植被定量遥感等研究中是十分重要的参数。
目前,测量叶面积指数的方法有两种:冠层水平的非接触光学测量反演的方法和地面采样-扫描叶片方法。其中冠层水平的非接触光学测量反演的方法则是利用普通相机或鱼眼镜头获得叶片的冠层孔隙度,进而通过孔隙率-叶面积指数的理论关系,反演冠层叶面积指数。
地面采样-扫描叶片方法的工作原理是:将叶片层置于激光器和光电变送器之间;在任意叶片上任意选取若干点作为拟合测试点进行测试。光电变送器将接收到的激光信号转换为电压信号,电压信号经数据采集卡采集至计算机;经若干次测试后,由计算机拟合得到电压值Y和重叠叶片层数X之间的关系函数:Y=A×exp(-B×X);其中A和B为拟合得到的正数A>0,B>0;根据测试得到的电压值,利用关系函数Y=A×exp(-B×X)计算测试点的重叠叶片层数X;计算测试点叶片面积=测试点叶片层数×光电变送器的接收面积;对各测试点叶片面积求和,和值为叶片层的总叶片面积;计算叶面积指数:叶面积指数=总叶片面积/叶片层的占地面积。
但是上述两种方法在对叶片采样时是对从样方内采集的全部叶片进行扫描,而采集的全部叶片并不代表样方内的全部叶片,使得得到的叶面积指数的准确度降低。
发明内容
有鉴于此,本发明实施例提供一种叶面积指数测量方法及系统,用于解决现有叶面积指数测量方法得到的叶面积指数的准确度降低的问题,技术方案如下:
本发明实施例提供一种叶面积指数测量方法,包括:
获取样方内采集叶片的图像,其中所述采集叶片的图像为从所述样方内采集的叶片不重叠平铺在采样板上形成的图像,所述样方为用于调查植物群路状况而随机抽取的取样地块;
以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积;
获取所述样方内残余叶片的图像,其中所述残余叶片的图像为从所述样方内采集叶片后剩余的叶片形成的样方图像;
基于所述残余叶片的图像,获取所述样方内残余叶片的面积;
将所述采集的全部叶片的面积和所述残余叶片的面积加和,得到所述样方内的叶片总面积;
依据公式:LAI=A/S得到所述样方的叶面积指数LAI,其中A为叶片总面积,S为样方面积。
优选地,所述获取样方内采集叶片的图像,包括:获取由商用扫描仪扫描到的所述样方内采集叶片的图像。
优选地,以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积,包括:
通过称重装置获取形成所述采样板上所述叶片的重量;
对所述采集叶片的图像中进行颜色分类;
利用颜色分类结果,提取所述采集叶片的图像中叶片分布区域,统计所述采集叶片的图像中的叶片面积;
利用公式AA=AS*WA/WS得到所述样方内所采集的全部叶片的面积AA,其中AS为所述叶片面积,WS为所述采样板上所述叶片的重量,WA为从所述样方内采集的全部叶片的重量。
优选地,统计所述采集叶片的图像中的叶片面积,包括:
获取所述颜色分类结果中叶片元素的像素数NUM_LEAF以及叶片元素所占面积A_PIXEL,其中所述叶片元素所占面积A_PIXEL为所述采样板面积A_PAPER除以所述颜色分类结果中所有元素的像素数NUM_PAPER;
将所述叶片元素的像素数NUM_LEAF和叶片元素所占面积A_PIXEL相乘,得到所述采集叶片的图像中的叶片面积。
优选地,所述获取样方内残余叶片的图像,包括:获取由照相机拍摄到的所述样方内残余叶片的图像。
优选地,以所述残余叶片的图像为依据,获取所述样方内所残余的全部叶片的面积,包括:
对所述残余叶片的图像按照样方范围进行剪裁,并获取所述剪裁后的图像中叶片的覆盖比例;
依据所述叶片的覆盖比例乘以样方面积,获取所述采集叶片的面积。
本发明实施例还提供一种叶面积指数测量系统,包括:
图像获取设备,用于获取样方内采集叶片的图像和样方内残余叶片的图像,其中所述采集叶片的图像为从所述样方内采集的叶片不重叠平铺在采样板上形成的图像,所述样方为用于调查植物群路状况而随机抽取的取样地块,所述残余叶片的图像为从所述样方内采集叶片后剩余的叶片形成的样方图像;
电子设备,用于以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积,基于所述残余叶片的图像,获取所述样方内残余叶片的面积,并将所述采集的全部叶片的面积和所述残余叶片的面积加和,得到所述样方内的叶片总面积,依据公式:LAI=A/S得到所述样方的叶面积指数LAI,其中A为叶片总面积,S为样方面积。
优选地,所述图像获取设备包括:商用扫描仪和照相机;
所述商用扫描仪扫描,用于扫描不重叠平铺在采样板上的采集叶片,得到所述样方内采集叶片的图像;
所述照相机,用于得到所述残余叶片的图像。
优选地,所述图像获取设备为照相机。
优选地,所述电子设备以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积,包括:
通过称重装置获取形成所述采样板上所述叶片的重量,对所述采集叶片的图像中进行颜色分类,利用颜色分类结果,提取所述采集叶片的图像中叶片分布区域,统计所述采集叶片的图像中的叶片面积,利用公式AA=AS*WA/WS得到所述样方内所采集的全部叶片的面积AA,其中AS为所述叶片面积,WS为所述采样板上所述叶片的重量,WA为从所述样方内采集的全部叶片的重量。
优选地,所述电子设备统计所述采集叶片的图像中的叶片面积,包括:获取所述颜色分类结果中叶片元素的像素数NUM_LEAF以及叶片元素所占面积A_PIXEL,其中所述叶片元素所占面积A_PIXEL为所述采样板面积A_PAPER除以所述颜色分类结果中所有元素的像素数NUM_PAPER;
将所述叶片元素的像素数NUM_LEAF和叶片元素所占面积A_PIXEL相乘,得到所述采集叶片的图像中的叶片面积。
与现有技术相比,本发明包括以下优点:
通过本发明实施例提供的叶面积指数测量方法,在获取样方内所采集的全部叶片的面积之外,还可以获取样方内残余叶片的面积,这样将采集的全部叶片的面积和残余叶片的面积加和后得到的叶片总面积基本等于样方内所有叶片的面积,相对于现有单纯以采集叶片的面积为叶片总面积来说,提高叶片总面积的准确度,基于此得到的叶面积指数的准确度相对提高。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的叶面积指数测量方法的流程图;
图2是图1中步骤102的流程图;
图3是基于图1所示方法测量到的叶面积指数的示意图;
图4是基于图1所示方法的过程示意图;
图5是实验结果的对比示意图;
图6是本发明实施例提供的叶面积指数测量系统的示意图。
具体实施方式
为了使本领域技术人员更好地理解本发明实施例,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,其示出了本发明实施例提供的叶面积指数测量方法的一种流程图,可以包括以下步骤:
101:获取样方内采集叶片的图像。
其中采集叶片的图像为从样方内采集的叶片不重叠平铺在采样板上形成的图像,样方为用于调查植物群路数量而随机抽取的取样地块。
在采样板上平铺采集的叶片时,可以将从样方内采集的全部叶片平铺到采样板上,也可以将样方内采集的全部叶片中抽取一部分叶片,然后将抽取的叶片平铺在采样板上。并且在平铺叶片时可以选择两片透明硬塑料纸片,其中一片作为衬底板、一片作为覆盖板,将叶片不重叠的铺贴在衬底板上并用覆盖板覆盖。
为了让叶片能够较牢固的固定在衬底板上,防止在移动衬底板时将平铺的叶片重叠,本发明实施例在平铺叶片时采用粘贴剂将叶片粘贴在衬底板上。
102:以采集叶片的图像为依据,获取样方内所采集的全部叶片的面积。在本发明实施例中,采集叶片的图像是从样方内采集的全部叶片的图像时,可以在获取到采集叶片的图像后,直接获取到全部叶片的面积。而当采集叶片的图像是从样方内采集的全部叶片中,抽取部分叶片后得到的图像时,则需要首先获取采集叶片的图像中叶片面积,然后再获取采集的全部叶片的面积。
当获取采集叶片的图像方式不同时,获取采集的全部叶片的面积也不同,其中获取采集叶片的图像的一种可行方式是:获取由商用扫描仪扫描到的样方内采集叶片的图像。由于商用扫描仪的通用性较高,所以在实际操作过程中,使用的商用扫描仪故障后,可以较容易找到另一个商用扫描仪代替,提高本发明实施例提供的叶面积指数测量方法的通用性。
相应的在由商用扫描仪获取采集叶片的图像时,其获取样方内采集的全部叶片的面积的过程可以参阅图2所示,可以包括以下步骤:
1021:通过称重装置获取形成采样板上叶片的重量。其中称重装置可以是天平或者电子称,具体可以依据实际操作过程中可以找到的称重装置获取形成采样板上叶片的重量。
1022:对采集叶片的图像中进行颜色分类。目前对采集叶片的图像进行颜色分类的方式可以包括非监督分类和监督分类,其中非监督分类(UnsupervisedClassification)是以不同影像物在特征空间中类别特征的差别为依据的一种无先验(已知)类别标准的图像分类,是以集群为理论基础,通过计算机对图像进行集聚统计分析的方法。监督分类(supervisedclassification)则是以建立统计识别函数为理论基础,依据典型样本训练方法进行分类的方法。
在本发明实施例可以采用非监督分类,例如利用ENVI软件(TheEnvironmentforVisualizingImages,遥感影像处理)的非监督分类中的k-means(k均值)模块对采集叶片的图像进行分类,鉴于采集叶片的图像中绿白(或灰白)两种颜色分布,设置非监督分类使用的k-means模块的类别为2类。
在本发明实施例中也可以采用监督分类方法,例如利用ENVI软件的监督分类中的MLC(MaximumLikelihoodClassification,最大似然分类)模块进行分类,鉴于采集叶片的图像中绿白(或灰白)两种颜色分布,选择叶片扫描区域和非叶片扫描区域作为兴趣区,参与MLC模块分类并获得叶片扫描区域和非叶片扫描区域的分类结果,叶片扫面区域和非叶片扫面区域以颜色区分,因此该分类结果也可以称为颜色分类结果。
1023:利用颜色分类结果,提取采集叶片的图像中叶片分布区域,统计采集叶片的图像中的叶片面积。
由于一幅图像中颜色分布有两种,所以可以依据该分类模块对采集叶片的图像进行分类时,以将图像中的叶片区域和背景区域区分。在本发明实施例中,统计采集叶片的图像中的叶片面积,包括:
首先,获取颜色分类结果中叶片元素的像素数NUM_LEAF以及叶片元素所占面积A_PIXEL,其中叶片元素所占面积A_PIXEL为采样板面积A_PAPER除以颜色分类结果中所有元素的像素数NUM_PAPER;
然后将叶片元素的像素数NUM_LEAF和叶片元素所占面积A_PIXEL相乘,得到采集叶片的图像中的叶片面积。
1024:利用公式AA=AS*WA/WS得到样方内所采集的全部叶片的面积AA,其中AS为叶片面积,WS为采样板上叶片的重量,WA为从样方内采集的全部叶片的重量,全部叶片的重量WA也可以由称重装置获取。
103:获取样方内残余叶片的图像,其中残余叶片的图像为从样方内采集叶片后剩余的叶片形成的样方图像。残余叶片的图像可以由照相机拍摄得到,并且在拍摄的时候要对样方进行垂直拍摄。
104:基于残余叶片的图像,获取样方内残余叶片的面积。在本发明实施例中,获取残余叶片的面积可以是利用ENVIMLC模块,选择图像中的叶片区域和背景区域作为兴趣区,这样可以从兴趣区中获取叶片区域,以区分开样方内残余叶片和背景;然后统计出图像中叶片区域的面积,即得到残余叶片的面积。
在本发明实施例中,获取残余叶片的面积还可以通过方式:对所述残余叶片的图像按照样方范围进行剪裁,并获取所述剪裁后的图像中叶片的覆盖比例;依据所述叶片的覆盖比例乘以样方面积,获取所述采集叶片的面积
105:将采集的全部叶片的面积和残余叶片的面积加和,得到样方内的叶片总面积。
106:依据公式:LAI=A/S得到样方的叶面积指数LAI,其中A为叶片总面积,S为样方面积。在本发明实施例中样方的形状可以随意设置,如正方形、圆形或三角形,样方面积S的计算可以是由人工计算得到。
通过上述技术方案,在获取样方内所采集的全部叶片的面积之外,还可以获取样方内残余叶片的面积,这样将采集的全部叶片的面积和残余叶片的面积加和后得到的叶片总面积基本等于样方内所有叶片的面积,基于此得到的叶面积指数的准确度相对提高。
以2011年在西藏低矮草地采集的样本N01-0716为例,其草平均高度低于5CM,样本N01-0716的LAI测定结果,总LAI计算结果为0.62,其中采样叶片的LAI为0.4,残余叶片照片提取的LAI为0.22,如图3所示。因为稀疏低矮草叶片难以采集,残余叶片量的比例非常高,本发例中残余叶片LAI占总样方LAI的35.5%,因此本发明实施例提供的叶面积指数测量方法将残余叶片计算在内,对于低矮草地叶面积指数的测定精度提高是明显的。
上述阐述的叶面积指数测量方法的过程可以如图4所示,首先人工对样方内的叶片进行采集,其次通过商用扫描仪对采集叶片进行扫描,得到采集叶片的图像,然后利用ENVI软件及权重法计算采集叶片的面积,然后获取残余叶片的图像得到残余叶片的面积,最后可以得到叶片总面积并利用公式LAI=A/S得到样方的叶面积指数LAI(图中最后一幅图片)。
在本发明实施例中,获取采集叶片的图像的另一种可行方式是:获取由照相机拍摄到的样方内采集叶片的图像。相应的获取全部叶片的面积的过程可以是:首先对采集叶片的图像进行剪裁,并获取剪裁后的图像中叶片的覆盖比例;然后依据叶片的覆盖比例和采样板的面积,获得采集的全部叶片的面积其中采样板为放置采集叶片的工具,如A4纸。
但是采用照相机获取采集叶片的图像会由于拍摄高度和角度不稳定性导致图像畸形,使得叶面积指数的准确度降低。下面将列举实验数据进行说明。
下面为一个叶面积实际测量的案例,以某光机扫描仪(如型号为LI-3000的叶面积扫描仪)、商用扫描仪和照相机(垂直、倾斜30度、倾斜40度)测量叶面积。并对光机扫描仪测量结果进行精密校验,校验得到的叶面积值(一般认为接近真值)为120.19cm-2,进而计算不同方法测量的相对精度。测量及精度计算结果如图5所示,发现实际使用的某光机扫描仪测量经度偏低,相对精度为-9.60%,而其他测量偏高,其中商用扫描仪的方法测量精度最高,相对精度为3.81%,照相法随着倾斜角度加大畸变效应增加,误差在加大。
与上述方法实施例相对应,本发明实施例还提供一种叶面积指数测量系统,其示意图如图6所示,可以包括:图像获取设备11和电子设备12,其中图像获取设备11和电子设备12之间通过数据线进行信息传输。
图像获取设备11,用于获取样方内采集叶片的图像和样方内残余叶片的图像,其中采集叶片的图像为从样方内采集的叶片不重叠平铺在采样板上形成的图像,样方为用于调查植物群路状况而随机抽取的取样地块,残余叶片的图像为从样方内采集叶片后剩余的叶片形成的样方图像。
在采样板上平铺采集的叶片时,可以将从样方内采集的全部叶片平铺到采样板上,也可以将样方内采集的全部叶片中抽取一部分叶片,然后将抽取的叶片平铺在采样板上。并且在平铺叶片时可以选择两片透明硬塑料纸片,其中一片作为衬底板、一片作为覆盖板,将叶片不重叠的铺贴在衬底板上并用覆盖板覆盖。
为了让叶片能够较牢固的固定在衬底板上,防止在移动衬底板时将平铺的叶片重叠,本发明实施例在平铺叶片时采用粘贴剂将叶片粘贴在衬底板上。
在本发明实施例中,图像获取设备11可以包括:商用扫描仪和照相机。商用扫描仪扫描,用于扫描不重叠平铺在采样板上的采集叶片,得到样方内采集叶片的图像。照相机,用于得到残余叶片的图像。
商用扫描仪和照相机得到的图像可以通过USB(UniversalSerialBus,通用串行总线)数据线由于商用扫描仪的通用性较高,所以在实际操作过程中,使用的商用扫描仪故障后,可以较容易找到另一个商用扫描仪代替,提高本发明实施例提供的叶面积指数测量方法的通用性,并且商用扫描仪的精度高于照相机的精度。当然在本发明实施例中,图像获取设备11可以为照相机来获取采集叶片的图像和残余叶片的图像。
电子设备12,用于以采集叶片的图像为依据,获取样方内所采集的全部叶片的面积,基于残余叶片的图像,获取样方内残余叶片的面积,并将采集的全部叶片的面积和残余叶片的面积加和,得到样方内的叶片总面积,依据公式:LAI=A/S得到样方的叶面积指数LAI,其中A为叶片总面积,S为样方面积。
其中,电子设备12以采集叶片的图像为依据,获取样方内所采集的全部叶片的面积的过程如下:
通过称重装置获取形成采样板上叶片的重量,对采集叶片的图像中进行颜色分类,利用颜色分类结果,提取采集叶片的图像中叶片分布区域,统计采集叶片的图像中的叶片面积,利用公式AA=AS*WA/WS得到样方内所采集的全部叶片的面积AA,其中AS为叶片面积,WS为采样板上叶片的重量,WA为从样方内采集的全部叶片的重量。
统计采集叶片的图像中的叶片面积的过程为:获取颜色分类结果中叶片元素的像素数NUM_LEAF以及叶片元素所占面积A_PIXEL,其中叶片元素所占面积A_PIXEL为采样板面积A_PAPER除以颜色分类结果中所有元素的像素数NUM_PAPER。再将叶片元素的像素数NUM_LEAF和叶片元素所占面积A_PIXEL相乘,得到采集叶片的图像中的叶片面积。
在实际操作过程中,电子设备12中可以运行图像分类软件(如ENVI)来获取颜色分类结果。例如利用ENVI软件(TheEnvironmentforVisualizingImages,遥感影像处理)的非监督分类中的k-means(k均值)模块对采集叶片的图像进行分类,鉴于采集叶片的图像中绿白(或灰白)两种颜色分布,设置非监督分类使用的k-means模块的类别为2类。
或者又例如利用ENVI软件的监督分类中的MLC(MaximumLikelihoodClassification,最大似然分类)模块进行分类,鉴于采集叶片的图像中绿白(或灰白)两种颜色分布,选择叶片扫描区域和非叶片扫描区域作为兴趣区,参与MLC模块分类并获得叶片扫描区域和非叶片扫描区域的分类结果,叶片扫面区域和非叶片扫面区域以颜色区分,因此该分类结果也可以称为颜色分类结果。
需要说明的是,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。对于装置类实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
最后,还需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
以上对本发明所提供的一种叶面积指数测量方法及系统进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (7)

1.一种叶面积指数测量方法,其特征在于,包括:
获取样方内采集叶片的图像,其中所述采集叶片的图像为从所述样方内采集的叶片不重叠平铺在采样板上形成的图像,所述样方为用于调查植物群路状况而随机抽取的取样地块;
以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积;
获取所述样方内残余叶片的图像,其中所述残余叶片的图像为从所述样方内采集叶片后剩余的叶片形成的样方图像;
基于所述残余叶片的图像,获取所述样方内残余叶片的面积;
将所述采集的全部叶片的面积和所述残余叶片的面积加和,得到所述样方内的叶片总面积;
依据公式:LAI=A/S得到所述样方的叶面积指数LAI,其中A为叶片总面积,S为样方面积;
所述获取样方内采集叶片的图像,包括:获取由商用扫描仪扫描到的所述样方内采集叶片的图像;
以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积,包括:
通过称重装置获取形成所述采样板上所述叶片的重量;
对所述采集叶片的图像中进行颜色分类;
利用颜色分类结果,提取所述采集叶片的图像中叶片分布区域,统计所述采集叶片的图像中的叶片面积;
利用公式AA=AS*WA/WS得到所述样方内所采集的全部叶片的面积AA,其中AS为所述叶片面积,WS为所述采样板上所述叶片的重量,WA为从所述样方内采集的全部叶片的重量。
2.根据权利要求1所述的方法,其特征在于,统计所述采集叶片的图像中的叶片面积,包括:
获取所述颜色分类结果中叶片元素的像素数NUM_LEAF以及叶片元素所占面积A_PIXEL,其中所述叶片元素所占面积A_PIXEL为所述采样板面积A_PAPER除以所述颜色分类结果中所有元素的像素数NUM_PAPER;
将所述叶片元素的像素数NUM_LEAF和叶片元素所占面积A_PIXEL相乘,得到所述采集叶片的图像中的叶片面积。
3.根据权利要求1所述的方法,其特征在于,所述获取样方内残余叶片的图像,包括:获取由照相机拍摄到的所述样方内残余叶片的图像。
4.根据权利要求3所述的方法,其特征在于,以所述残余叶片的图像为依据,获取所述样方内所采集的全部叶片的面积,包括:
对所述残余叶片的图像进行按照样方范围进行剪裁,并获取所述剪裁后的图像中残余叶片的覆盖比例;
依据所述残余叶片的覆盖比例乘以样方面积,获取所述残余叶片的面积。
5.一种叶面积指数测量系统,其特征在于,包括:
图像获取设备,用于获取样方内采集叶片的图像和样方内残余叶片的图像,其中所述采集叶片的图像为从所述样方内采集的叶片不重叠平铺在采样板上形成的图像,所述样方为用于调查植物群路状况而随机抽取的取样地块,所述残余叶片的图像为从所述样方内采集叶片后剩余的叶片形成的样方图像;
电子设备,用于以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积,基于所述残余叶片的图像,获取所述样方内残余叶片的面积,并将所述采集的全部叶片的面积和所述残余叶片的面积加和,得到所述样方内的叶片总面积,依据公式:LAI=A/S得到所述样方的叶面积指数LAI,其中A为叶片总面积,S为样方面积;
所述图像获取设备包括:商用扫描仪和照相机;
所述商用扫描仪扫描,用于扫描不重叠平铺在采样板上的采集叶片,得到所述样方内采集叶片的图像;
所述电子设备以所述采集叶片的图像为依据,获取所述样方内所采集的全部叶片的面积,包括:
通过称重装置获取形成所述采样板上所述叶片的重量,对所述采集叶片的图像中进行颜色分类,利用颜色分类结果,提取所述采集叶片的图像中叶片分布区域,统计所述采集叶片的图像中的叶片面积,利用公式AA=AS*WA/WS得到所述样方内所采集的全部叶片的面积AA,其中AS为所述叶片面积,WS为所述采样板上所述叶片的重量,WA为从所述样方内采集的全部叶片的重量。
6.根据权利要求5所述的系统,其特征在于,所述图像获取设备为照相机。
7.根据权利要求5或6所述的系统,其特征在于,所述电子设备统计所述采集叶片的图像中的叶片面积,包括:获取所述颜色分类结果中叶片元素的像素数NUM_LEAF以及叶片元素所占面积A_PIXEL,其中所述叶片元素所占面积A_PIXEL为所述采样板面积A_PAPER除以所述颜色分类结果中所有元素的像素数NUM_PAPER;
将所述叶片元素的像素数NUM_LEAF和叶片元素所占面积A_PIXEL相乘,得到所述采集叶片的图像中的叶片面积。
CN201410229812.3A 2014-05-28 2014-05-28 一种叶面积指数测量方法及系统 Expired - Fee Related CN103983217B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410229812.3A CN103983217B (zh) 2014-05-28 2014-05-28 一种叶面积指数测量方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410229812.3A CN103983217B (zh) 2014-05-28 2014-05-28 一种叶面积指数测量方法及系统

Publications (2)

Publication Number Publication Date
CN103983217A CN103983217A (zh) 2014-08-13
CN103983217B true CN103983217B (zh) 2016-07-06

Family

ID=51275283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410229812.3A Expired - Fee Related CN103983217B (zh) 2014-05-28 2014-05-28 一种叶面积指数测量方法及系统

Country Status (1)

Country Link
CN (1) CN103983217B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106018284B (zh) * 2016-05-12 2018-10-30 中国科学院遥感与数字地球研究所 一种作物长势定量遥感监测方法及系统
CN106468536B (zh) * 2016-11-22 2019-08-09 天津市气候中心 一种叶面积测量方法
CN110068299B (zh) * 2019-05-05 2020-10-13 华北水利水电大学 一种温室作物叶面积指数的计算方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102243069A (zh) * 2011-06-22 2011-11-16 华南农业大学 一种叶面积指数测定方法及装置
CN102331244A (zh) * 2011-08-10 2012-01-25 浙江大学 利用普通数码相机获取水稻冠层叶面积指数的方法
CN102538717A (zh) * 2010-12-30 2012-07-04 北京师范大学 叶面积指数自动观测系统及其方法
CN102778212A (zh) * 2012-08-02 2012-11-14 浙江大学 一种利用数码相机获取水稻冠层叶面积指数的方法及装置
CN102997871A (zh) * 2012-11-23 2013-03-27 南京大学 一种利用几何投影和激光雷达反演有效叶面积指数的方法
US20130182911A1 (en) * 2010-12-02 2013-07-18 Nec Corporation Leaf area index measurement system, device, method, and program

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130182911A1 (en) * 2010-12-02 2013-07-18 Nec Corporation Leaf area index measurement system, device, method, and program
CN102538717A (zh) * 2010-12-30 2012-07-04 北京师范大学 叶面积指数自动观测系统及其方法
CN102243069A (zh) * 2011-06-22 2011-11-16 华南农业大学 一种叶面积指数测定方法及装置
CN102331244A (zh) * 2011-08-10 2012-01-25 浙江大学 利用普通数码相机获取水稻冠层叶面积指数的方法
CN102778212A (zh) * 2012-08-02 2012-11-14 浙江大学 一种利用数码相机获取水稻冠层叶面积指数的方法及装置
CN102997871A (zh) * 2012-11-23 2013-03-27 南京大学 一种利用几何投影和激光雷达反演有效叶面积指数的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
叶面积指数的主要测定方法;谭一波等;《林业调查规划》;20080630;第33卷(第3期);第1节-第3节 *
基于数字照片特征的小麦覆盖度自动提取研究;李存军等;《浙江大学学报》;20040630;第30卷(第6期);全文 *

Also Published As

Publication number Publication date
CN103983217A (zh) 2014-08-13

Similar Documents

Publication Publication Date Title
CN109489637B (zh) 水量变化监测方法、装置、计算机设备和存储介质
Malambo et al. Automated detection and measurement of individual sorghum panicles using density-based clustering of terrestrial lidar data
Tang et al. Three-dimensional surface reconstruction of tree canopy from lidar point clouds using a region-based level set method
AU2014256411B2 (en) Yield measurement and base cutter height control systems for a harvester
McCombs et al. Influence of fusing lidar and multispectral imagery on remotely sensed estimates of stand density and mean tree height in a managed loblolly pine plantation
CN110988909B (zh) 基于tls进行高寒脆弱区沙地植被的植被盖度测定方法
Carr et al. Individual tree segmentation from a leaf-off photogrammetric point cloud
Latella et al. A density-based algorithm for the detection of individual trees from LiDAR data
CN109815833A (zh) 一种基于ccd与三维激光传感信息融合特征提取的茶尖识别方法
CN103983217B (zh) 一种叶面积指数测量方法及系统
CN109325431A (zh) 草原放牧绵羊采食路径中植被覆盖度的检测方法及其装置
Brocks et al. Toward an automated low-cost three-dimensional crop surface monitoring system using oblique stereo imagery from consumer-grade smart cameras
Kathuria et al. Development of an automated individual tree detection model using point cloud LiDAR data for accurate tree counts in a Pinus radiata plantation
CN111397709A (zh) 一种小麦千粒重快速测量方法
Tao et al. Recognition of ecological vegetation fairy circles in intertidal salt marshes from UAV LiDAR point clouds
Jónsson RGB and Multispectral UAV image classification of agricultural fields using a machine learning algorithm
Smigaj et al. Capturing hedgerow structure and flowering abundance with UAV remote sensing
Sánchez-Lopera et al. Classification of lidar bare-earth points, buildings, vegetation, and small objects based on region growing and angular classifier
Mason et al. Measurement of habitat predictor variables for organism-habitat models using remote sensing and image segmentation
CN112017057A (zh) 一种保险理赔处理方法及装置
CN117671498A (zh) 一种雷电防护用无人机自动检测系统
CN105116407B (zh) 一种利用手持型激光测距仪测量植被覆盖度的方法
CN116721344A (zh) 基于航拍设备的植被检测方法、装置和设备
CN116258844A (zh) 一种棉花叶片表型性状快速精准鉴定方法
Tang et al. A new method for extracting trees and buildings from sparse LiDAR data in urban areas

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160706

Termination date: 20190528