尹晓爱, 杨秀超, 侯 瑞,2†, 赵龙山,2,, 张金鑫
(1.贵州大学林学院,550025,贵阳;2.贵州大学土壤侵蚀与生态修复研究中心,550025,贵阳;3.中国林科院林业研究所,100091,北京)
叶面积指数(leaf area index, LAI)是指覆盖单位地表面积与所生长的植物叶面积总和之比。LAI既可反映作物群体大小,又是表征作物冠层结构的重要指标[1]。在植被对气候变化响应、遥感作物估产及植被水土保持作用等研究中LAI被广泛采用[2]。为此,关于LAI测量与计算方法一直受到人们的关注。
植物冠层叶面积的测量方法主要有直接法与间接法2种[2]。直接测量法包括手工原位测量、落叶法和系数法等;间接测量方法主要基于光学原理,通过获取植物冠层截获光能量的多少来估算或通过测量入射光在冠层中的衰减量进行估算[1-2]。二者相比,直接法较耗时,而且对冠层结构会造成干扰,影响测量的可重复性和结果的准确性,而间接法需要借助先进的仪器设备完成。虽然间接测量仪器的应用在一定程度上提高了作物冠层参数的获取效率,但该类方法计算流程复杂。如当前使用较多的是植被冠层分析仪,通过测量冠层内透过的太阳辐射量来估算冠层LAI[2-3],但该仪器价格昂贵,野外携带也不方便。
随着智能手机的普及,图像分辨率、内存容量和人机交互功能等不断提高,在野外调查中利用智能手机收集图像、视频资料已较为普遍,故开展基于智能手机的LAI测量和计算方法研究也受到人们关注。如郭文川等[4]基于智能手机开发叶面积测量软件,无需购置硬件设备,也不需要给手机增加额外的硬件,具有很强的便携性。高大帅等[5]利用智能手机方位传感器确定拍摄角度,建立误差校正模型并提出基于智能手机方位传感器的叶面积测量校正公式。
尽管利用智能手机测定冠层参数有一些研究,但对LAI的测量和计算方法仍在进行不断探索和改进;鉴此,本研究探讨一种基于智能手机的玉米LAI快速测量方法,在试验条件下通过对不同生育期玉米植株冠层LAI的模型计算结果与实测结果比较,验证该方法的有效性,为进一步完善利用智能手机测定冠层参数方法和开发有效冠层LAI测定手机APP提供科技条件。
以坡耕地种植玉米为研究对象,选择5个重要生长期(拔节期、小喇叭期、大喇叭期、抽雄期、成熟期)的玉米植株各30株,进行单株冠层结构特征的测量,将测量数据用于制作玉米植株模型(表1)。玉米叶片采用仿真叶模拟,其材质为聚乙烯,透光率约为10%,茎秆采用PVC管代替,将仿真叶片粘贴在PVC管,叶片倾角为30°,茎秆底部固定在木板上。根据田间实际测量数据模拟不同生育期完整植株,模拟玉米株行距分别设置为30 cm×50 cm,株数分别为1、2、3、6和9株,以代表不同采样面积。
表1 不同生育期单株玉米模型数量特征Tab.1 Quantitative characteristics of maize model per plant at different growth stages
试验于2018年3—11月在贵州大学土壤侵蚀模拟实验室进行。实验室为暗室条件,用白炽灯模拟太阳光照,模拟太阳高度角分别为30°、60°和90°,试验场景见图1。先将模拟玉米植株按设定株行距置于坐标纸上,用智能手机进行玉米冠层投影拍摄(使用设备为iPhone 7 Plus MNRM2CH/A型智能手机,手机内置iOS 14.0操作系统,手机屏幕大小5.5 inch,1 200万像素双摄像头,运行内存32 GB)。同时,在坐标纸上勾绘玉米植株实际投影边界,以便对手机拍摄投影值进行建模修正,设置3次重复试验。
图1 试验场景示意图(∠1、∠2和∠3分别为30°、60° 和90°)Fig.1 Schematic diagram of test scene (measure 1, measure 2 and measure 3 are 30°, 60° and 90°, respectively)
玉米冠层实际投影在坐标纸勾绘后,利用扫描仪对坐标纸进行扫描存入计算机,采用Adobe Photoshop CS6软件提取阴影,为便于区分对图像中阴影区域填充为绿色(RGB:0,255,0)的图像;对手机拍摄的冠层阴影图像也采用同样方法提取图像中阴影区域,并填充为红色(RGB:255,0,0)。利用像元比例关系分别计算出投影面积,分别记为SC(坐标纸勾绘的冠层投影面积)和SM(手机拍摄的冠层投影面积),后通过回归分析建立SC与SM之间的函数关系,对手机拍摄的冠层投影面积进行校对,进而可直接通过手机拍摄的冠层投影面积计算冠层LAI。
玉米冠层LAI计算公式为
(1)
式中:LAI为玉米冠层叶表面积指数,量纲为1;Sleaf为玉米冠层叶面积,cm2;SA为土地面积,cm2;本研究中其值为坐标纸面积。
本研究中模拟玉米植株冠层叶面积真值(Sleaf)通过仿真叶片实际面积计算,由于各生长期模拟玉米植株以实际玉米植株参数制作(表1),故模拟的不同生长期玉米叶面积与实际坡耕地玉米接近。5个生育期单个模拟玉米植株叶面积Sleaf分别为177.03、799.62、2 061.25、5 868.40和6 581.03 cm2。土地面积(SA)即为试验中坐标纸面积(9万112.00 cm2)。相应地,1、2、3、6和9株模拟玉米冠层LAI真值(LAIT)见表2。
表2 不同生育期模拟玉米冠层LAI真值 (LAIT)Tab.2 LAI actual value (LAIT) of maize per plant atdifferent growth stages
采用Excel 2016软件绘制图表;采用SPSS 26.0对玉米冠层LAI进行多因素方差分析;通过回归方法对基于手机拍摄图像、坐标纸勾绘的投影面积与LAI真值进行模型分析。
由表3可知,随着太阳高度角的增大,坐标纸勾绘的玉米冠层投影面积与实际面积间的比值(以下简称勾绘面积比)逐渐减小,太阳高度角为30°时勾绘的面积比最接近1,60°下的面积比次之,90°下勾绘的面积比值偏离1的程度较大。不同玉米生长时期与种植株数勾绘面积比变化不大,在太阳高度角为30°时,相比玉米的其他4个生长期,小喇叭期的勾绘面积比最接近1。
表3 坐标纸勾绘的冠层投影面积与实际面积之比Tab.3 Ratio of the projected area of the canopy to theactual area plotted on a coordinate paper
以玉米模型参数(生长时期、种植株数和太阳高度角)为因子,坐标纸勾绘的玉米冠层投影面积为变量进行多因素方差分析,结果见表4。玉米生长期、种植株数分别与坐标勾绘的投影面积之间具有显著相关关系(P<0.05);太阳高度角与坐标勾绘的投影面积之间关系不显著(P>0.05);玉米生长期与种植株数间的相互效应对坐标勾绘的投影面积具有显著的影响(P<0.05);其余2个交互效应(生长期×太阳高度角,太阳高度角×种植株数)对坐标勾绘投影面积的影响不显著(P>0.05)。
表4 各模型参数下坐标勾绘的冠层投影面积相关性分析Tab.4 Analysis of canopy projection area plotted withcoordinates under various model parameters
总体上,太阳高度角为60°时,手机拍摄的玉米冠层投影面积与实际面积间的比值(以下简称拍摄面积比)较接近1,在太阳高度角为90°时次之,太阳高度角为30°时拍摄面积比值偏离1的程度最大(表5)。不同玉米生长期与种植株数拍摄面积比变化波动不大,在太阳高度角60°时,相比玉米的其他4个生长期,大喇叭期下的拍摄面积比最接近1。
表5 手机拍摄的冠层投影面积与实际面积之比
以玉米模型参数(生长时期、种植株数和太阳高度角)为因子,手机拍摄的玉米冠层投影面积为变量进行多因素方差分析,结果见表6。由表6可知,玉米生长期、种植株数分别与手机拍摄的冠层投影面积之间具有显著相关关系(P<0.05);太阳高度角与手机拍摄的投影面积之间关系不显著(P>0.05);玉米生长期与种植株数间的相互效应对手机拍摄的投影面积具有显著的影响(P<0.05);其余2个交互效应(生长期×太阳高度角,太阳高度角×种植株数)对投影面积的影响不显著(P>0.05)。
表6 各玉米模型参数下手机拍摄的冠层投影面积相关性分析
图2 以坐标纸勾绘的叶面积来修正手机拍摄的冠层投影面积Fig.2 Use the leaf area drawn on coordinate paper to correct the canopy projection area taken by mobile phone
图3 坐标纸勾绘的冠层投影面积与实际面积间的回归分析Fig.3 Regression analysis between the projected area and the actual area of the canopy plotted on coordinate paper
LAIE=
(2)
图4 LAI估算值与真值间的关系Fig.4 Relationship between estimated and true LAI values
1)玉米生长期、种植株数及二者的交互效应对玉米冠层叶面积测定具有显著的影响(P<0.05),太阳高度角的影响不显著(P>0.05)。
本研究提出了一种基于智能手机的冠层LAI测定方法,通过模拟试验表明该方法可操作性强,测量结果具有一定的准确性和有效性,可为植物冠层LAI测定提供新的途径。