CN111798509A - 一种基于半球图像法测量叶面积指数的方法 - Google Patents
一种基于半球图像法测量叶面积指数的方法 Download PDFInfo
- Publication number
- CN111798509A CN111798509A CN202010572717.9A CN202010572717A CN111798509A CN 111798509 A CN111798509 A CN 111798509A CN 202010572717 A CN202010572717 A CN 202010572717A CN 111798509 A CN111798509 A CN 111798509A
- Authority
- CN
- China
- Prior art keywords
- image
- pixel
- point
- theta
- canopy
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000003384 imaging method Methods 0.000 claims abstract description 16
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 230000011218 segmentation Effects 0.000 claims description 14
- 241000566145 Otus Species 0.000 claims description 8
- 238000012417 linear regression Methods 0.000 claims description 5
- 230000007797 corrosion Effects 0.000 claims description 3
- 238000005260 corrosion Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000010801 machine learning Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 5
- 230000000903 blocking effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000005068 transpiration Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于半球图像法测量叶面积指数的方法,先通过鱼眼摄像头拍摄植被冠层图像并进行预处理,再确定出植被冠层图像的有效成像区域的几何中心和圆形有效成像区域的半径,从而对鱼眼镜头进行投影方程标定;然后采集待测点的植被冠层图像并处理,再根据投影方程提取出多个角度的冠层孔隙度,进而反演出叶面积指数。
Description
技术领域
本发明属于生态学研究技术领域,更为具体地讲,涉及一种基于半球图像法测量叶面积指数的方法。
背景技术
叶面积指数是一个无量纲量,可认为叶面积与地面积之间的比值。它反映了叶片的数量和分布情况,是研究冠层中光线穿透情况,冠层生产力,冠层下土壤水分蒸发、蒸腾损失总量,冠层截留,及土壤温度的基础因子。
目前叶面积指数的测量方法有直接测量法和间接测量法。直接测量法通常来说包括收集样本和叶面积测量两个步骤,虽然精确度高,但是具有很强的破坏性和工作量大的缺点,仅限于小冠层测量,对于树林、森林大冠层基本上是不可能的。间接测量法通过光线传播模型能相对简单地测量出光线的穿透情况,并基于Lambert-Beer定律对叶面积指数做出精确的估计。间接测量方法主要有斜点样方法、基于辐射测定法和孔隙度反演法,具有快速、简单等优势。其中斜点样方的问题是需要大量的专用探针且操作相当繁琐,因而不适应冠层较大的植物的参数获取;基于辐射测定法常常需要借助激光雷达(根据搭载平台分为地基、车载、机载、星载等类型)获取目标对象的空间信息,处理数据复杂,成本高昂。而孔隙度反演法的代表方法——半球图像法,通过鱼眼镜头获取接近180°冠层图像信息,经过图像处理获取冠层孔隙度,进而反演出叶面积指数。自Even等将其用于植被的冠层的辐射监测以来,半球图像法以其超大视角范围和获取鱼眼图像设备价廉的优势得到了广泛的研究和应用。近年来,随着数码相机和数字图像处理技术的完善和进步,半球图像法的功能也日益完善,目前已成为各领域的重要研究方法。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于半球图像法测量叶面积指数的方法,在受时空约束条件下,快速、准确测量出叶面积指数。
为实现上述发明目的,本发明一种基于半球图像法测量叶面积指数的方法,其特征在于,包括以下步骤:
(1)、图像采集及预处理;
采用鱼眼摄像头拍摄植被冠层图像,然后对植被冠层图像先进行灰度化处理,再使用固定阈值法进行二值化处理,最后进行先膨胀后腐蚀的处理,从而得到轮廓清晰的二值化图像;
(2)、确定植被冠层图像的有效成像区域的几何中心O1(u0,v0)和圆形有效成像区域的半径R;
(2.1)、在二值化图像中,以左上角第一个像素点作为像素坐标系的原点O,建立像素坐标系xOy;
(2.2)、设置有效数据集,初始化为空;
(2.3)、从左往右逐行逐列扫描二值化图像,当扫描到某一像素点(x,y)的像素值为255时,进一步判断像素点(x+1,y)的像素值是否也为255,如果这两个像素点的像素值都为255,则将像素点(x,y)识别为轮廓边界点,并记录像素点(x,y)坐标保存至有效数据集,否则继续扫描,直到最后一个像素点扫描结束;
从右往左逐行逐列扫描二值化图像,当扫描到某一点(x,y)的像素值为255时,进一步判断像素点(x-1,y)的像素值是否也为255,如果这两个像素点的像素值都为255,则将像素点(x,y)识别为轮廓边界点,并记录像素点(x,y)坐标保存至有效数据集,否则继续扫描,直到最后一个像素点扫描结束;
(2.4)、在像素坐标系xOy中,构建圆方程x2+y2=Ax+By+C;
(2.5)、以x2+y2作为标签,x,y,1作为特征值,利用有效数据集训练基于机器学习的多元变量线性回归模型,得到参数A,B,C;
(3)、对鱼眼镜头进行投影方程标定;
(3.1)、选取8×11的黑白棋盘板作为标准面,设置上、下四个相邻黑白棋格交线的交点为标定点,即黑白棋格图像的角点;在与标准面垂直的轴线上固定好镜头位置,使得透镜光轴垂直于标定板平面,拍摄标定板图像;
(3.2)、用激光笔标识好镜头中心在标准面上的投射点O2;
(3.5)、调整标定板与鱼眼镜头中心的距离,在不同距离下重复(3.1)-(3.4)的操作,从而记录n组(r,θ)数据;
(4)、采集待测点的植被冠层图像并处理;
(4.1)、采用鱼眼摄像头拍摄待测点的植被冠层图像,并进行灰度化处理;
(4.2)、对植被冠层的灰度图像进行二值化处理:
(4.2.1)、以图像圆心O1为起点,将灰度图像中的圆形有效成像区域按照拍摄的方位角和入射角分别进行分块,其中,灰度图像按拍摄的方位角一共分为m1个区域块,每个区域块又按入射角分为m2个区域块,共分为m1×m2个扇形块;
(4.2.2)、利用改进的otus算法获取每一个扇形块的二值化阈值T,并按照二值化阈值T进行二值化分割,得到待测点的植被冠层二值化图像;
(5)、提取多个角度的冠层孔隙度;
(5.1)、根据投影方程r=fθ,在待测点的植被冠层二值化图像中按照(0~θ0)入射角范围取n个环,每个环对应宽度dθ=θ°/n的入射角区间,由内到外,各环的半径为r1=f·θ0/n,r2=2f·θ0/n,r3=3f·θ0/n......rn=f·θ0,各环的入射角依次用θ1=1/(2n)θ0,θ2=3/(2n),θ0θ3=5/(2n)θ0......θn=(2n-1)/(2n)θ0表示;
(5.2)、计算各入射角下的冠层孔隙度T(θi);
其中,T(θi)表示入射角θi下的冠层孔隙度,Nleaf(θi)表示在入射角θi下背景像素点个数,Nsky(θi)表示在入射角θi下叶片像素点个数;
(6)、利用反演公式计算叶面积指数;
叶面积指数的反演公式为:
本发明的发明目的是这样实现的:
本发明一种基于半球图像法测量叶面积指数的方法,先通过鱼眼摄像头拍摄植被冠层图像并进行预处理,再确定出植被冠层图像的有效成像区域的几何中心和圆形有效成像区域的半径,从而对鱼眼镜头进行投影方程标定;然后采集待测点的植被冠层图像并处理,再根据投影方程提取出多个角度的冠层孔隙度,进而反演出叶面积指数。
同时,本发明一种基于半球图像法测量叶面积指数的方法还具有以下有益效果:
(1)、一般的LAI测量仪器拍摄条件最好满足万里无云,且避免太阳直射条件下进行;本发明的阈值分割算法具有自适应的优点,对光照条件的要求没有特别的限制,因此对测量环境要求较低,适合在各种天气条件下进行测量;
(2)、基于辐射测量的方法需要同时获取冠层上方和冠层底部的辐射量,当冠层较为高大时获取冠层上方的辐射不易实现,这种情况往往需要在附近空旷地测得的辐射量替代冠层上方的辐射量;且需要搭配激光雷达来获取冠层的三维空间信息,成本较高;与基于辐射测定的方法相比,本发明计算叶面积指数的方法只需要在冠层下方拍摄一次,对采集设备要求较低,只需要一个装配鱼眼镜头即可完成采集,具有价格低廉,且便于携带的优点;
(3)、本发明测量叶面积指数的方法使用冠层孔隙度来替代冠层间隙率,避免了基于辐射测定法测量冠层上方、下方辐射量的不便和对外部设备的依赖。只需通过数字图像处理技术,就可以很轻松地获取冠层孔隙度信息;
(4)、本发明采用的改进Otus算法对图像分块进行二值化处理,减少因光照、鱼眼相机曝光度等因素对图像二值化品质的影响,提高了冠层孔隙度提取的精确度,进而提高叶面积指数计算的准确性。
附图说明
图1是本发明基于半球图像法测量叶面积指数的方法流程图;
图2是为了提取鱼眼图像轮廓信息的经过图像处理后的图像;
图3是使用圆方程拟合拟合鱼眼图像有效区域的效果图;
图4是建立投影方程图;
图5是某一次测量采集到的冠层图像;
图6是叶片对光线的反射率和透射比与波长的关系图;
图7是冠层图像分块示意图;
图8是各种阈值分割效果对比图;
图9是本发明算法测量结果与LAI2200C测量结果的对比图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图1是本发明基于半球图像法测量叶面积指数的方法流程图。
在本实施例中,如图1所示,本发明一种基于半球图像法测量叶面积指数的方法,包括以下步骤:
S1、图像采集及预处理;
采用鱼眼摄像头拍摄植被冠层图像,然后对植被冠层图像先进行灰度化处理,再使用固定阈值法(固定阈值设为30)进行二值化处理,最后进行先膨胀后腐蚀的处理,如图2所示,从而得到轮廓清晰的二值化图像;
S2、确定植被冠层图像的有效成像区域的几何中心O1(u0,v0)和圆形有效成像区域的半径R;
S2.1、在二值化图像中,以左上角第一个像素点作为像素坐标系的原点O,建立像素坐标系xOy;
S2.2、设置有效数据集,初始化为空;
S2.3、从左往右逐行逐列扫描二值化图像,当扫描到某一像素点(x,y)的像素值为255时,进一步判断像素点(x+1,y)的像素值是否也为255,如果这两个像素点的像素值都为255,则将像素点(x,y)识别为轮廓边界点,并记录像素点(x,y)坐标保存至有效数据集,否则继续扫描,直到最后一个像素点扫描结束;
从右往左逐行逐列扫描二值化图像,当扫描到某一点(x,y)的像素值为255时,进一步判断像素点(x-1,y)的像素值是否也为255,如果这两个像素点的像素值都为255,则将像素点(x,y)识别为轮廓边界点,并记录像素点(x,y)坐标保存至有效数据集,否则继续扫描,直到最后一个像素点扫描结束;
S2.4、在像素坐标系xOy中,构建圆方程x2+y2=Ax+By+C;在本实施例中,如图2所示,可以看出由于光线分布不均的影响,有的地方轮廓信息提取不完整,所以需要上述步骤进一步筛选有效的轮廓信息来进行圆方程拟合,最后通过圆方程拟合鱼眼图像有效区域的效果图,如图3所示;
S2.5、以x2+y2作为标签,x,y,1作为特征值,利用有效数据集训练基于机器学习的多元变量线性回归模型,得到参数A,B,C;
S3、对鱼眼镜头进行投影方程标定;
研究表明鱼眼镜头成像时遵循的模型可以近似为单位球面投影模型,可以将鱼眼相机的成像过程分解成两步:第一步,三维空间点线性地投影到一个球面上,它是一个虚拟的单位球面,它的球心与相机坐标系的原点重合;第二步,单位球面上的点投影到图像平面上。鱼眼镜头大多采用等距投影模型设计,这种投影方式最大的特点是来自入射角θ的一点m在成像平面上的像点(m')到图像中心O的距离r与入射角θ成线性关系,用方程表示为:r=f·θ;那么,对鱼眼镜头进行投影方程标定的具体过程为:
S3.1、选取8×11的黑白棋盘板作为标准面,设置上下相邻四个黑白棋格交线的交点为标定点,即黑白棋格图像的角点;在与标准面垂直的轴线上固定好镜头位置,使得透镜光轴垂直于标定板平面,拍摄标定板图像;
S3.2、用激光笔标识好镜头中心在标准面上的投射点O2;
S3.5、调整标定板与鱼眼镜头中心的距离,在不同距离下重复(3.1)-(3.4)的操作,从而记录n=352组(r,θ)数据;
S4、采集待测点的植被冠层图像并处理;
S4.1、采用鱼眼摄像头拍摄待测点的植被冠层图像,如图5所示,共计4幅,并进行灰度化处理;在本实施例中,叶片对光线的反射率和透射比与波长的关系如图6所示,因为在小于490nm区域的光线受叶片反射和透射作用最小,这使得叶簇在天空背景下是黑色的,从而更容易区分天空和植物。蓝光波段在400-450nm,满足我们应用需求。为了更好地分割天空像元和树叶像元,因此选择RGB冠层图像的蓝色分量进行灰度化处理。
S4.2、对植被冠层的灰度图像进行二值化处理:
S4.2.1、以图像圆心O1为起点,将灰度图像中的圆形有效成像区域按照拍摄的方位角和入射角分别进行分块,其中,灰度图像按拍摄的方位角一共分为m1个区域块,每个区域块又按入射角分为m2个区域块,共分为m1×m2个扇形块;
在本实施例中,如图7所示,按拍摄的方位角一共分为8个区域块,每个区域块又按入射角分为10个区域块,共分为80个扇形块;这样分块的目的好处在于某一分块区域灰度分布直方图都存在明显的波峰波谷,从而可以更好地应用Otus算法进行后续阈值分割;
S4.2.2、利用改进的otus算法获取每一个扇形块的二值化阈值T,然后将每一个扇形块中像素值大于等于阈值T的灰度值设置为255,并归为天空背景;将像素值小于阈值T的灰度值设置为0,归为叶片,得到待测点的植被冠层图像的二值化图像。
具体过程为:
其中,w1为天空像素占比,Ni为某一扇形块内灰度值为i的像素点个数,Sum为某一扇形块内总像素点个数;w2为叶面像素占比,w2=1-w1;μ1为天空像素的平均灰度值,pi为某一扇形块内灰度值为i的像素占比,μ2为叶面像素的平均灰度值,
在本实施例中,图8给出了各种阈值分割效果对比图,其中,图(a)RGB原图;(b)固定阈值法分割,阈值为100;在圈出的两个地方,错误地将树叶像元识别为天空像元,阈值的选择对固定阈值分割影响较大,而阈值的选择往往只有通过人眼看到视觉效果人为调节,不具备实际应用的要求;(c)otus算法,因为无效区域大量黑色像素的影响,使得otus算法计算出的阈值往往偏低,往往容易将灰度较暗的天空像元识别为树叶像元;(d)基于高斯加权和的自适应阈值分割,邻域大小为81×81;该算法对于树干只提出了轮廓部分,很多细节丢失。这是因为高斯加权和的自适应阈值分割算法要求处理的图片,灰度不能太复杂,尽量符合混合高斯模型;(e)本文提出的阈值算法,在原有的otus算法上进行了改进,未将有效区域外的像素值参与计算,消除了无效区域大量黑色像素导致阈值偏低的影响,同时结合鱼眼镜头光线分布与入射角和方位角的关系,对图像进行分块处理,具有自适应性,阈值分割效果较好,满足本文对阈值分割的需求。
S5、提取多个角度的冠层孔隙度;
S5.1、根据投影方程r=fθ,在待测点的植被冠层二值化图像中按照(0~θ0)入射角范围取5个环,每个环对应宽度dθ=75°/5=15°的入射角区间,由内到外各环的半径为r1=f*15,r2=f*30,r3=f*45,r4=f*60,r5=f*75,各环的入射角依次用θ1=7°,θ2=23°,θ3=38°,θ4=53°,θ5=68°表示。
S5.2、在假设叶片是黑色不透光的前提下,冠层孔隙度是光线穿透冠层不被捕获的概率,具体计算各入射角下的冠层孔隙度T(θi)为:
其中,T(θi)表示入射角θi下的冠层孔隙度,Nleaf(θi)表示在入射角θi下背景像素点个数,Nsky(θi)表示在入射角θi下叶片像素点个数;
S6、利用反演公式计算叶面积指数;
叶面积指数反演的理论公式为:
其中,p(θ)为入射角θ下的冠层间隙率,在利用半球图像法计算叶面积指数时通常使用冠层孔隙度T(θ)代替冠层间隙率p(θ)。
在实际计算中,通过利用5个环(θ1=7°,θ2=23°,θ3=38°,θ4=53°,θ5=68°,dθ=15°)的孔隙度信息对上面的公式离散化处理,得到实际的叶面积指数反演公式:
如表1所示,各环的权重因子如下:
表1
在本实施例中,采集到学校四个样本地不同植被类型近50组冠层图像,并同时用LAI-2200C测量每个实验点的叶面积指数。对采集到的冠层图像,采用本发明的算法进行叶面积指数的估算。然后对本发明的算法结果和LAI-2200C的测量结果进行线性回归、相关性等分析,如图9所示。从全体数据分析结果(图9(a))来看,本文基于半球图像法LAI的计算结果与LAI-2200C的LAI测量结果之间相关系数达到了0.90244。证明本论文的算法与LAI-2200C的测量结果之间具有线性关系,说明本算法的精确度较高(约0.9)。
从更深层次上分析,两者的较高相关性从源头上是半球图像法采用的孔隙度变量与LAI-2200C的间隙率变量之间较高的相关性。在将冠层看作泊松分布的假设下,两者从不同角度去量化树叶对光束的拦截概率。基于这点出发,可以用线性模型回归训练的参数对本文算法结果进行修正,即xr=kx+b,其中xr为本文算法测量LAI修正后的结果,k,b分别为图9(a)线性回归模型的斜率参数和截距参数。修正后的结果具有较高的精确度和准确度。
通过与LAI2200C测量结果进行比较,二者相关系数达到0.9,肯定了本发明提供的基于半球图像法测量叶面积指数方法的可行性和可信性。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (2)
1.一种基于半球图像法测量叶面积指数的方法,其特征在于,包括以下步骤:
(1)、图像采集及预处理;
采用鱼眼摄像头拍摄植被冠层图像,然后对植被冠层图像先进行灰度化处理,再使用固定阈值法进行二值化处理,最后进行先膨胀后腐蚀的处理,从而得到轮廓清晰的二值化图像;
(2)、确定植被冠层图像的有效成像区域的几何中心O1(u0,v0)和圆形有效成像区域的半径R;
(2.1)、在二二值化图像中,以左上角第一个像素点作为像素坐标系的原点O,建立像素坐标系xOy;
(2.2)、设置有效数据集,初始化为空;
(2.3)、从左往右逐行逐列扫描二值化图像,当扫描到某一像素点(x,y)的像素值为255时,进一步判断像素点(x+1,y)的像素值是否也为255,如果这两个像素点的像素值都为255,则将像素点(x,y)识别为轮廓边界点,并记录像素点(x,y)坐标保存至有效数据集,否则继续扫描,直到最后一个像素点扫描结束;
从右往左逐行逐列扫描二值化图像,当扫描到某一点(x,y)的像素值为255时,进一步判断像素点(x-1,y)的像素值是否也为255,如果这两个像素点的像素值都为255,则将像素点(x,y)识别为轮廓边界点,并记录像素点(x,y)坐标保存至有效数据集,否则继续扫描,直到最后一个像素点扫描结束;
(2.4)、在像素坐标系xOy中,构建圆方程x2+y2=Ax+By+C;
(2.5)、以x2+y2作为标签,x,y,1作为特征值,利用有效数据集训练基于机器学习的多元变量线性回归模型,得到参数A,B,C;
(3)、对鱼眼镜头进行投影方程标定;
(3.1)、选取8×11的黑白棋盘板作为标准面,设置上、下四个相邻黑白棋格交线的交点为标定点,即黑白棋格图像的角点;在与标准面垂直的轴线上固定好镜头位置,使得透镜光轴垂直于标定板平面,拍摄标定板图像;
(3.2)、用激光笔标识好镜头中心在标准面上的投射点O2;
(3.5)、调整标定板与鱼眼镜头中心的距离,在不同距离下重复(3.1)-(3.4)的操作,从而记录n组(r,θ)数据;
(4)、采集待测点的植被冠层图像并处理;
(4.1)、采用鱼眼摄像头拍摄待测点的植被冠层图像,并进行灰度化处理;
(4.2)、对植被冠层的灰度图像进行二值化处理:
(4.2.1)、以图像圆心O1为起点,将灰度图像中的圆形有效成像区域按照拍摄的方位角和入射角分别进行分块,其中,灰度图像按拍摄的方位角一共分为m1个区域块,每个区域块又按入射角分为m2个区域块,共分为m1×m2个扇形块;
(4.2.2)、利用改进的otus算法获取每一个扇形块的二值化阈值T,并按照二值化阈值T进行二值化分割,得到待测点的植被冠层二值化图像;
(5)、提取多个角度的冠层孔隙度;
(5.1)、根据投影方程r=fθ,在待测点的植被冠层二值化图像中按照(0~θ0)入射角范围取n个环,每个环对应宽度dθ=θ°/n的入射角区间,由内到外,各环的半径为r1=f·θ0/n,r2=2f·θ0/n,r3=3f·θ0/n......rn=f·θ0,各环的入射角依次用θ1=1/(2n)θ0,θ2=3/(2n),θ0θ3=5/(2n)θ0......θn=(2n-1)/(2n)θ0表示;
(5.2)、计算各入射角下的冠层孔隙度T(θi);
其中,T(θi)表示入射角θi下的冠层孔隙度,Nleaf(θi)表示在入射角θi下背景像素点个数,Nsky(θi)表示在入射角θi下叶片像素点个数;
(6)、利用反演公式计算叶面积指数;
叶面积指数的反演公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010572717.9A CN111798509B (zh) | 2020-06-22 | 2020-06-22 | 一种基于半球图像法测量叶面积指数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010572717.9A CN111798509B (zh) | 2020-06-22 | 2020-06-22 | 一种基于半球图像法测量叶面积指数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111798509A true CN111798509A (zh) | 2020-10-20 |
CN111798509B CN111798509B (zh) | 2021-12-10 |
Family
ID=72803660
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010572717.9A Active CN111798509B (zh) | 2020-06-22 | 2020-06-22 | 一种基于半球图像法测量叶面积指数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111798509B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112381863A (zh) * | 2020-11-12 | 2021-02-19 | 中国电建集团江西省电力设计院有限公司 | 森林场景的地面激光点云方法 |
CN114821325A (zh) * | 2022-04-24 | 2022-07-29 | 南京林业大学 | 基于模型提取植被冠层孔隙度的方法、装置、设备及介质 |
CN117929467A (zh) * | 2024-03-25 | 2024-04-26 | 上海农林职业技术学院 | 一种基于热成像技术的叶片水分数据检测系统及方法 |
CN118397075A (zh) * | 2024-06-24 | 2024-07-26 | 合肥工业大学 | 基于鱼眼相机山地森林有效叶面积指数的计算方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101916438A (zh) * | 2010-07-20 | 2010-12-15 | 浙江大学 | 半球摄影法获取水稻冠层叶面积指数和平均叶倾角的方法 |
CN102855625A (zh) * | 2012-08-02 | 2013-01-02 | 浙江大学 | 一种半球摄影法获取水稻冠层丛生指数的方法 |
CN102997871A (zh) * | 2012-11-23 | 2013-03-27 | 南京大学 | 一种利用几何投影和激光雷达反演有效叶面积指数的方法 |
CN104200457A (zh) * | 2014-07-21 | 2014-12-10 | 华南农业大学 | 基于广角摄像的离散型冠层叶面积指数检测系统和方法 |
CN104330058A (zh) * | 2014-11-25 | 2015-02-04 | 东北林业大学 | 一种利用植物冠层分析仪测定叶面积指数的校正方法 |
CN104330057A (zh) * | 2014-11-25 | 2015-02-04 | 东北林业大学 | 一种利用半球摄影法测定叶面积指数的校正方法 |
CN104778693A (zh) * | 2015-04-08 | 2015-07-15 | 云挺 | 基于投影算法和活动轮廓模型的叶面积指数计算 |
JP2016131494A (ja) * | 2015-01-15 | 2016-07-25 | 株式会社日立ソリューションズ | 生育管理装置 |
CN106770289A (zh) * | 2016-12-09 | 2017-05-31 | 东华大学 | 一种用于分布式叶面积指数半球图像采集的节点 |
CN107610066A (zh) * | 2017-09-21 | 2018-01-19 | 电子科技大学 | 一种叶面积指数测量方法 |
CN109859231A (zh) * | 2019-01-17 | 2019-06-07 | 电子科技大学 | 一种基于光学图像的叶面积指数提取阈值分割方法 |
-
2020
- 2020-06-22 CN CN202010572717.9A patent/CN111798509B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101916438A (zh) * | 2010-07-20 | 2010-12-15 | 浙江大学 | 半球摄影法获取水稻冠层叶面积指数和平均叶倾角的方法 |
CN102855625A (zh) * | 2012-08-02 | 2013-01-02 | 浙江大学 | 一种半球摄影法获取水稻冠层丛生指数的方法 |
CN102997871A (zh) * | 2012-11-23 | 2013-03-27 | 南京大学 | 一种利用几何投影和激光雷达反演有效叶面积指数的方法 |
CN104200457A (zh) * | 2014-07-21 | 2014-12-10 | 华南农业大学 | 基于广角摄像的离散型冠层叶面积指数检测系统和方法 |
CN104330058A (zh) * | 2014-11-25 | 2015-02-04 | 东北林业大学 | 一种利用植物冠层分析仪测定叶面积指数的校正方法 |
CN104330057A (zh) * | 2014-11-25 | 2015-02-04 | 东北林业大学 | 一种利用半球摄影法测定叶面积指数的校正方法 |
JP2016131494A (ja) * | 2015-01-15 | 2016-07-25 | 株式会社日立ソリューションズ | 生育管理装置 |
CN104778693A (zh) * | 2015-04-08 | 2015-07-15 | 云挺 | 基于投影算法和活动轮廓模型的叶面积指数计算 |
CN106770289A (zh) * | 2016-12-09 | 2017-05-31 | 东华大学 | 一种用于分布式叶面积指数半球图像采集的节点 |
CN107610066A (zh) * | 2017-09-21 | 2018-01-19 | 电子科技大学 | 一种叶面积指数测量方法 |
CN109859231A (zh) * | 2019-01-17 | 2019-06-07 | 电子科技大学 | 一种基于光学图像的叶面积指数提取阈值分割方法 |
Non-Patent Citations (2)
Title |
---|
ZHOU L ET AL: "Methodology comparison for effective LAI retrieving based on digital hemispherical photograph in rice Canopy", 《SPRING BERLIN HEIDELBERG》 * |
吴伟斌等: "基于半球图像法的叶面积指数反演研究进展", 《科学技术与工程》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112381863A (zh) * | 2020-11-12 | 2021-02-19 | 中国电建集团江西省电力设计院有限公司 | 森林场景的地面激光点云方法 |
CN114821325A (zh) * | 2022-04-24 | 2022-07-29 | 南京林业大学 | 基于模型提取植被冠层孔隙度的方法、装置、设备及介质 |
CN114821325B (zh) * | 2022-04-24 | 2024-04-02 | 南京林业大学 | 基于模型提取植被冠层孔隙度的方法、装置、设备及介质 |
CN117929467A (zh) * | 2024-03-25 | 2024-04-26 | 上海农林职业技术学院 | 一种基于热成像技术的叶片水分数据检测系统及方法 |
CN118397075A (zh) * | 2024-06-24 | 2024-07-26 | 合肥工业大学 | 基于鱼眼相机山地森林有效叶面积指数的计算方法 |
CN118397075B (zh) * | 2024-06-24 | 2024-09-06 | 合肥工业大学 | 基于鱼眼相机山地森林有效叶面积指数的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111798509B (zh) | 2021-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111798509B (zh) | 一种基于半球图像法测量叶面积指数的方法 | |
Zhu et al. | Improving leaf area index (LAI) estimation by correcting for clumping and woody effects using terrestrial laser scanning | |
Woodgate et al. | Understanding the variability in ground-based methods for retrieving canopy openness, gap fraction, and leaf area index in diverse forest systems | |
Simonse et al. | Automatic determination of forest inventory parameters using terrestrial laser scanning | |
Seidel et al. | Analyzing forest canopies with ground-based laser scanning: A comparison with hemispherical photography | |
CN106772417B (zh) | 一种动目标多维度多尺度红外光谱特征测量方法及系统 | |
Grotti et al. | An intensity, image-based method to estimate gap fraction, canopy openness and effective leaf area index from phase-shift terrestrial laser scanning | |
Hassaan et al. | Precision forestry: trees counting in urban areas using visible imagery based on an unmanned aerial vehicle | |
Tuominen et al. | Local radiometric correction of digital aerial photographs for multi source forest inventory | |
Larsen | Crown modelling to find tree top positions in aerial photographs | |
CN112396019A (zh) | 一种基于无人机的植被分布识别方法、系统及可读存储介质 | |
CN112285710A (zh) | 一种多源遥感水库蓄水量估算方法与装置 | |
CN116295285A (zh) | 基于区域自适应的浅海水深遥感反演方法 | |
CN114332085B (zh) | 一种光学卫星遥感影像检测方法 | |
CN108318458B (zh) | 一种适用于不同天气条件下的室外典型地物pBRDF的测量方法 | |
Larsen | Finding an optimal match window for Spruce top detection based on an optical tree model | |
Wang et al. | The impact of variable illumination on vegetation indices and evaluation of illumination correction methods on chlorophyll content estimation using UAV imagery | |
Schwalbe et al. | Hemispheric image modeling and analysis techniques for solar radiation determination in forest ecosystems | |
US20230377163A1 (en) | Image-based method of defining a scanning area | |
CN108416760A (zh) | 一种震后高分辨率sar图像单体建筑物倒塌损毁检测方法 | |
CN111175231A (zh) | 冠层植被指数的反演方法、装置及服务器 | |
Del Pozo et al. | Multispectral imaging: Fundamentals, principles and methods of damage assessment in constructions | |
CN113469104A (zh) | 基于深度学习的雷达遥感影像地表水体变化检测方法和设备 | |
Yun et al. | Study of subtropical forestry index retrieval using terrestrial laser scanning and hemispherical photography | |
Kedzierski et al. | Integration of point clouds and images acquired from a low-cost NIR camera sensor for cultural heritage purposes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |