-
云南松(Pinus yunnanensis Franch.)是我国西南地区主要树种,云南省主要森林植被类型之一,占林地面积52%,在林业生产中占有举足轻重的地位[1]。松小蠹(Tomicus spp.)是云南松的主要害虫之一,云南至今已有超过600 km2的云南松林死亡[2]。松小蠹蛀梢和蛀干习性使其长时间隐藏在树体内部,不易发现和防治。5—11月为蛀梢为害期,成虫蛀食当年生枝梢;11月至次年4月为蛀干为害期,成虫在树干韧皮部和木质部钻蛀、产卵,蛀梢和蛀干交替出现,历时长,危害严重。松小蠹还具备迁飞习性,扩散快,导致松林成片死亡[3]。
高光谱遥感技术主要用于研究植被受病虫害危害前后变化及受害程度和植被光谱曲线、植被指数变化之间的关系,确定不同树种的危害敏感波段和时期,建立病虫害监测模型。刘畅[4]、Ju[5]等利用非成像高光谱仪测量锐齿栎、马尾松叶片光谱曲线,测定其光合色素含量,建立了光谱和光合的联系。王震[6]、伍南[7]等基于高光谱特征建立了马尾松松材线虫病、赤枯病的病情程度反演模型。地面成像光谱能获取地物光谱数据及实现“图谱合一”,利于混合光谱分析和病虫害监测模型反演。徐爽[8]、Shafri[9]等利用地面成像高光谱仪实现了长枣表面虫伤、灵芝基茎腐病监测。于泉州[10]、宋小宁[11]基于高光谱数据进行森林资源调查,提取森林冠层光谱特征、冠层含水量等信息。刘良云等[12]利用多时相高光谱航空影像,实现了冬小麦锈病监测。
目前,利用遥感技术监测松小蠹危害的研究多基于多光谱遥感进行,而对地面高光谱成像仪监测松小蠹危害情况的关注较少。因此,本研究通过实地调查松小蠹危害时期和发生情况,室内采集云南松枝梢针叶光谱反射率,分析光谱特征,建立松小蠹危害时期监测模型,实现危害早期监测。并结合野外测定的不同危害时期云南松针叶光合参数,建立光谱与光合的联系,利用光谱数据实现光合参数预测,估算云南松的生长健康状况。
-
试验地位于石林彝族自治县西街口镇,海拔1 700~1 950 m,地处103°29′20″ E,24°49′25″ N。亚热带低纬度高原山地季风气候,干湿分明,5—10月为雨季,11月到次年4月为旱季,年均气温14.7 ℃。试验地以云南松为优势种,属针叶林,大多为中幼龄林,约80%为纯林,林分单一,且云南松大多生长在土壤贫瘠、水分条件缺乏,立地条件差的地段,易受病虫害侵扰,具体采样地点如图1所示。
-
试验材料为受松小蠹危害的云南松枝梢,根据松小蠹对云南松的实际危害情况,以单梢为单位,将松小蠹危害时期分为4个时期,健康时期、危害早期、危害中期与危害末期(表1)[13]。通过前期调查,根据外观特征分别选择不同危害时期的云南松枝梢,每个时期各选择30个枝梢。
表 1 松小蠹危害时期划分
Table 1. Hazard time division of Tomicus spp.
危害时期Hazard periods 特征描述Feature description 健康时期Health 云南松枝梢的针叶为绿色,未受到松小蠹的危害,枝梢上没有蛀孔 危害早期Early damage 松小蠹开始蛀食枝梢,枝梢出现蛀孔,但枝梢的针叶仍显绿色或嫩绿色 危害中期Mid-hazard 松小蠹仍停留在枝梢内,即将完成对枝梢的蛀食,枝梢的针叶颜色为嫩黄色或黄绿色 危害末期Terminal stage 松小蠹完成对枝梢的蛀食,针叶逐渐发生变化,枝梢的针叶颜色表现为红褐色,最后呈现灰褐色 -
本研究采用的高光谱成像系统由SOC710VP光谱仪、2个500 W的钨灯、计算机、暗箱等组成。针叶光谱测定在可控光照条件的实验室进行,光谱仪架设在暗箱中,用钨灯进行照明。在图像采集时,设置以下参数:曝光时间35 ms,光圈大小5.6,物镜高度50 cm;钨灯光源入射点与镜头垂直于载物台的焦点重合。将针叶平铺在载物台上,打开钨灯,用标准参考板校正。启动光谱仪,进行光谱采集。
-
2018年10—12月,于晴朗天气下午15:00—17:00,在野外使用LI-6400便携式光合作用测定系统测定健康时期、危害早期、危害中期3个时期的云南松枝梢针叶净光合速率(Photo)、气孔导度(Cond)、胞间CO2浓度(Ci)及蒸腾速率(Trmmol)等参数(危害末期云南松枝梢接近枯死,故不测定其气体交换参数)。在自然光源条件下,选取阳面枝梢针叶进行测定,每组处理中每个枝梢重复测定5次,取平均值。利用Photo和Trmmol计算水分利用效率[14](WUE=Photo/Trmmol);利用Ci和空气中的CO2浓度(Ca)计算气孔限制值[15](Ls=1−Ci/Ca)。
-
利用SRAnal710软件,对高光谱图像数据进行光谱、黑场、空间光谱辐射标定等预处理,将DN值转换为反射率。采用Savitzky-Golay滤波法对针叶光谱反射率曲线进行平滑,以消除原始数据中的噪声干扰[16],在OriginPro8.0中绘制针叶光谱曲线。
采用Pearson相关分析与逐步回归分析筛选出特征变量,用以建模。本研究基于前人研究基础,结合云南松光谱反射率数据定义了8个植被指数、23个高光谱特征参数[17-23](表2)。
表 2 特征变量表
Table 2. Characteristic variable scale
变量类别Variable categories 变量定义和公式Variable definitions and formulas 植被指数Vegetation indexes $NDVI1 = {{({{\bar R}_{(650\sim 670)}} - {{\bar R}_{(760\sim 850)}})} / {({{\bar R}_{(650\sim 670)}} + {{\bar R}_{(760\sim 850)}})}}$ $NDVI2 = {{({R_{838}} - {R_{677}})} / {({R_{838}} + {R_{677}})}}$ $NDV{I_{705}} = {{({R_{750}} - {R_{705}})} / {({R_{750}} + {R_{705}})}}$ $SR1 = {{{{\bar R}_{(650\sim 670}}} / {{{\bar R}_{(760\sim 850)}}}}$ $SR2 = {{{R_{838}}} / {{R_{677}}}}$ $WBI = {{{R_{900}}} / {{R_{970}}}}$,水波段指数 $PRI = {{({R_{531}} - {R_{570}})} / {({R_{531}} + {R_{570}})}}$,光化学植被指数 $HI = {{({R_{739}} - R{}_{399})} / {({R_{739}} + {R_{399}}) - 0.5{R_{404}}}}$,植被健康指数 高光谱特征参数
Hyperspectral Features反射率位置特征参数
Reflectance location parametersRg(绿峰),波长520~560 nm范围内最大的波段反射率 Rr(红谷),波长 640~680 nm范围内最小的波段反射率 $D = 1 - \dfrac{{{R_{670}}}}{{{R_{560}} + \dfrac{{{R_{760}} - {R_{560}}}}{{760 - 560}} \times (670 - 560)}}$ $H = 1 - \dfrac{{{R_{500}} + \dfrac{{{R_{670}} - {R_{500}}}}{{670 - 500}} \times (560 - 500)}}{{{R_{560}}}}$ 一阶微分位置特征参数
First derivative location parametersDb,蓝边(470~520 nm)内一阶微分光谱中的最大值 Dy,黄边(560~590 nm)内一阶微分光谱中的最大 Dr,红边(660~740 nm)内一阶微分光谱中的最大值 Dnir,近红外(760~1036 nm)内一阶微分光谱中的最大值 λr,红边(660~740 nm)内一阶微分光谱中的最大值对应的波段 一阶微分面积特征参数
First derivative area parametersSDb,蓝边内一阶微分值的总和 SDy,黄边内一阶微分值的总和 SDr,红边内一阶微分值的总和 SDnir,近红外范围内一阶微分值的总和 植被指数特征参数
Vegetation indexes parametersRg/Rr D/H SDr/SDb SDr/SDy SDnir/SDb SDnir/SDr (Rg-Rr)/(Rg+Rr) (D-H)/(D+H) (SDr-SDb)/(SDr+SDb) (SDrSDy)/(SDr+SDy) -
将数据量的70%作为建立模型的数据,30%作为精度验证的数据。主要涉及以下模型:
多元线性回归模型:y=β0+β1x1+β2x2+β3x3+...+βnxn
线性函数模型:y=β0+β1x
二次函数模型:y=β0+β1x+β2x2
三次函数模型:y=β0+β1x+β2x2+β3x3
式中:y代表拟合值;x, x1, x2, x3, ..., xn代表主要特征参数;β0, β1, β2, β3, ..., βn代表拟合系数。
本研究使用决定系数R2和均方根误差RMSE进行精度验证。R2值在0~1之间,R2越接近1,说明预测模型参考价值越高。RMSE值越小,说明预测模型精度更高。
-
不同危害时期针叶光谱反射率见图2。针叶光谱反射率变化可反映松小蠹危害时期变化,健康时期、危害早期和危害中期针叶光谱反射率曲线具有一致性:在740~1 036 nm,随着危害时期变化,反射率逐渐降低;在680~740 nm(“红边”),健康时期斜率最大,然后依次为危害早期和危害中期;在510~560 nm(“绿峰”)和640~680 nm(“红谷”),健康时期光谱反射率低于危害中期和危害早期。危害末期光谱反射率曲线“绿峰”、“红谷”消失,大致为一条缓慢上升的曲线。从整体上看,反射率大小依次为健康时期、危害早期、危害中期、危害末期。
-
不同危害时期云南松光谱一阶微分曲线见图3,在500~540 nm,健康时期、危害早期和危害中期存在明显波峰,且值大小依次为危害早期、危害中期、健康时期,而危害末期曲线不存在明显波峰;在660~740 nm,存在明显波峰,且峰值达到最高,其值大小依次为健康时期、危害早期、危害中期。随着危害时期变化,出现“红边”“蓝移 ”现象,峰值向短波方向移动且逐渐降低;在760~1 036 nm,各危害时期的波峰、波谷出现波段范围大致相似。
-
松小蠹危害时期与光谱反射率、光谱一阶微分值相关性分析见图4、图5,在399、409、503~539、553~574、610~677、687~697、708~791、807~838、849~870、881~988、1 004、1 025~1 036 nm,松小蠹危害时期与针叶光谱反射率及光谱一阶微分值极显著相关,可有效反映松小蠹危害时期,利用以上波段构建的高光谱参数建立松小蠹危害时期监测模型是可行的。
-
通过相关分析和逐步回归分析(a)分别对植被指数和高光谱特征参数筛选用于建模的变量,获得显著或极显著变量:NDVI2、PRI;D、Dr、SDb、SDy(表3)。将植被指数和高光谱特征变量共同输入SPSS中,进行逐步回归分析(b),获得显著或极显著变量:D、Dr、NDVI705(表3)。
表 3 植被指数变量和高光谱特征参数变量筛选结果
Table 3. The screening results of vegetation index and hyperspectral characteristic parameter
特征变量Characteristic variables 相关分析Correlation analysis 逐步回归分析Stepwise regression analysis(a) 特征变量Characteristic variables 相关分析Correlation analysis 逐步回归分析Stepwise regression analysis(b) NDVI2 0.918** <0.001** D 0.923** <0.001** PRI 0.740** 0.003** Dr 0.911** <0.001** Dr 0.911** 0.019* NDVI705 0.773** <0.001** D 0.923** 0.031* (SDr-SDb)/(SDr+SDb) 0.267 <0.001** SDb 0.381* 0.001** SDy -0.887** 0.005** 注:*: P<0.05,显著相关;**: P<0.01,极显著相关,下同。
Note: P < 0.05, significant correlation; * *: P < 0.01, extremely significant correlation, the same below.基于筛选的变量分别建立松小蠹危害时期的预测模型,并进行精度验证,根据拟合精度和预测精度的R2、RMSE选择精度最高的预测模型(表4),其中以D、Dr、NDVI705为变量的多元线性回归模型精度最高。植被指数类参数中,以NDVI2为变量的模型精度更高;高光谱特征位置参数中,以D为变量的模型精度更高。几个模型中,多元线性回归模型精度较高。
表 4 危害时期监测模型(n=84)及检验(n=36)
Table 4. Monitoring models(n = 84) and validation(n = 36) of damage time
变量类型
Variable categories监测模型
Monitor models回归方程
Regression equations拟合精度
Fitted precision预测精度
Predicted precisionR2 RMSE R2 RMSE 植被指数
Vegetation indexNDVI2 $y = 0.355 + 4.460{x_1} - 6.974{x_1}^2 + 7.410{x_1}^3$ 0.872 0.415 0.854 0.424 PRI $y = 3.637 + 13.108{x_2} - 31.072{x_2}^2 - 208.144{x_2}^3$ 0.623 0.711 0.636 0.791 NDVI2、PRI $y = 0.622 + 3.662{x_1} + 3.780{x_2}$ 0.878 0.391 0.889 0.357 高光谱特征参数
Hyperspectral characteristic
parametersD $y = 1.139 + 0.811{x_3} - 3.174{x_3}^2 + 6.235{x_3}^3$ 0.900 0.367 0.895 0.498 Dr $y = 0.755 + 83.599{x_4} + 27\;496.449{x_4}^2 - 1\;321\;854.674{x_4}^3$ 0.839 0.465 0.796 0.517 SDb $y = - 0.802 + 935.581{x_5} - 67\;441.244{x_5}^2 + 1\;450\;214.907{x_5}^2$ 0.441 0.867 0.106 1.010 SDy $y = 2.019 - 235.870{x_6} + 6\;602.191{x_6}^2 + 1\;404\;056.992{x_6}^3$ 0.828 0.481 0.804 0.508 D、Dr、SDb、SDy $y = 1.343 + 1.121{x_3} + 97.985{x_4} - 47.574{x_5} - 80.177{x_6}$ 0.916 0.314 0.918 0.331 综合Comprehensive D、Dr、NDVI705 $y = 0.540 + 1.107{x_3} + 108.960{x_4} + 2.545{x_7}$ 0.937 0.276 0.925 0.258 D、Dr、NDVI705、(SDr−SDb)/(SDr+SDb) $y = 0.886 + 1.332{x_3} + 86.860{x_4} + 3.636{x_7} - 0.794{x_8}$ 0.969 0.189 0.946 0.212 注:y代表松小蠹的危害时期预测值,x1、x2、x3、x4、x5、x6、x7、x8分别代表NDVI2、PRI、D、Dr、SDb、SDy、NDVI705、(SDr−SDb)/(SDr+SDb)。
Note: y represents the predicted value of damage period, x1, x2, x3, x4, x5, x6, x7 and x8 represent NDVI2, PRI, D, Dr, SDb, SDy, NDVI705, (SDr − SDb) / (SDr + SDb), respectively.基于逐步回归分析(b)得到的4个参数(D、Dr、NDVI705、(SDr−SDb)/(SDr+SDb))建立危害时期的多元线性回归模型(表4),其模型精度大于以D、Dr、NDVI705为变量的多元线性回归模型。对上述4个参数进行建模重要性分析,重要性依次为0.129、0.114、0.575、0.183,结果表明:(SDr−SDb)/(SDr+SDb)参数虽未与松小蠹的危害时期显著相关,但仍具有一定重要性,将其纳入模型变量中可提高模型精度。
-
主要光合指标见表5,健康时期、危害早期、危害中期针叶Photo分别为4.621、2.089、0.643 μmol·m−2s−1,健康时期是危害中期的7倍,受松小蠹危害后,云南松针叶Photo逐渐下降。方差分析及Duncan’s多重比较显示:松小蠹不同危害时期的Photo、Cond、Ci及Ls存在显著差异(P<0.05),健康时期、危害早期、危害中期的值依次下降;相较于健康时期针叶,危害早期和危害中期Trmmol明显下降,存在显著差异;危害中期针叶WUE显著低于健康时期和危害早期。结合松小蠹不同时期针叶Ci和Ls分析,表明:针叶Photo降低是松小蠹危害云南松导致叶肉细胞中与光合作用相关的酶活性下降引起,气孔因素并非主要因素,其下降是为了减少针叶蒸腾引起的水分损失。
表 5 云南松受松小蠹危害的不同时期的针叶气体交换参数的影响
Table 5. The effects of gas exchange parameters in different periods of Pinus yunnanensis damaged by Tomicus spp.
主要光合指标Principal photosynthetic indexs 云南松受松小蠹危害的不同时期Different periods of harm 健康时期Health 危害早期Early damage 危害中期Mid-hazard 净光合速率Photo/(μmol·m−2s−1) 4.621 ± 0.105a 2.089 ± 0.087b 0.643 ± 0.048c 气孔导度Cond/(mol·m−2s−1) 0.047 ± 0.001a 0.023 ± 0.002b 0.016 ± 0.001c 胞间CO2浓度Ci/(μmol·mol−1) 129.098 ± 19.591a 174.854 ± 48.175b 235.635 ± 26.579c 蒸腾速率Trmmol/(mmol·m−2s−1) 1.938 ± 0.014a 1.053 ± 0.111b 0.889 ± 0.034b 水分利用效率WUE/(μmol·mmol−1) 2.385 ± 0.053a 2.631 ± 0.201a 0.755 ± 0.058b 气孔限制值Ls 0.578 ± 0.011a 0.425 ± 0.026b 0.277 ± 0.013c 注:表中数据为测得数据的平均值±标准误差,同行中不同字母表示不同处理间Duncan’s多重比较差异显著(P<0.05,n =90)。
Note: the data in the table is the mean ± standard error of the measured data. Different letters in the same line indicate that there are significant differences in Duncan's multiple comparisons between different treatments (P < 0.05, n = 90). -
不同危害时期针叶Photo与其他光合参数的相关性发生了不同程度变化(表6)。简单相关分析表明:3个不同危害时期针叶Photo与Cond、Ci、Trmmol均呈不同相关性,Photo与Ci呈显著负相关(r健康时期=−0.989**,r危害早期=−0.476**,r危害中期=−0.508**),而Photo与其他因子在不同危害时期呈不一致相关性和显著性。偏相关分析表明:健康时期、危害早期针叶Photo与Ci均呈极显著负相关(r健康=−0.732**,r危害早期=−0.730**),而危害中期呈正相关(r危害中期=0.069);健康时期针叶Photo与Cond呈显著负相关(r健康=−0.410*),危害早期呈极显著负相关(r危害早期=−0.520**),而危害中期呈正相关(r危害中期=0.205);针叶Photo与Trmmol在3个不同危害时期呈不一致相关性和显著性。说明云南松受到松小蠹危害后针叶净光合速率变化是多个因子共同作用的结果。
表 6 针叶主要光合参数的相关性分析
Table 6. The correlation analysis of main photosynthetic parameters of needles
危害时期Hazard periods 因子Factors 净光合速率Photo 气孔导度Cond 胞间CO2浓度Ci 蒸腾速率Trmmol 健康时期Health 净光合速率Photo 1 0.708** −0.989** 0.383 气孔导度Cond −0.410* 1 0.623** 0.371 胞间CO2浓度Ci −0.732** 0.315 1 0.486** 蒸腾速率Trmmol 0.218 0.967** 0.497** 1 危害早期Early damage 净光合速率Photo 1 0.309 −0.476** −0.320 气孔导度Cond −0.520** 1 −0.278 0.999** 胞间CO2浓度Ci −0.730** 0.528** 1 0.291 蒸腾速率Trmmol −0.543** 0.999** 0.551** 1 危害中期Mid-hazard 净光合速率Photo 1 −0.360 −0.508** 0.702** 气孔导度Cond 0.205 1 −0.519** 0.663** 胞间CO2浓度Ci 0.069 0.124 1 0.691** 蒸腾速率Trmmol 0.579** 0.512** 0.401 1 注:主对角线左下方为简单相关系数,右上方为偏相关系数。 对松小蠹危害时期和针叶主要光合参数进行简单相关分析、偏相关分析及复相关分析(表7)。简单相关分析表明:危害时期与Photo、Cond、Trmmol呈极显著正相关(rPhoto=0.943**,rCond=−0.797**,rTrmmol=0.688**),与Ci呈极显著负正相关(rCi=−0.764**)。偏相关分析表明:危害时期与Photo呈极显著正相关(r=0.588**),与其他光合参数相关性均未达到显著性水平;复相关分析表明危害时期与光合参数的相关性达到极显著水平(r=0.909**)。
表 7 松小蠹危害时期和针叶主要光合参数的相关性分析
Table 7. The correlation between damage time of Tomicus spp.andmain photosynthetic parameters
因子Factors 危害时期Hazard periods 净光合速率Photo 气孔导度Cond 胞间CO2浓度Ci 蒸腾速率Trmmol 危害时期Hazard periods 0.909** 0.943** 0.797** −0.764** 0.688** 净光合速率Photo 0.588** 0.926** 0.758** −0.821** 0.635** 气孔导度Cond 0.187 0.293** 0.979** −0.433** 0.972** 胞间CO2浓度Ci −0.155 −0.484** −0.016 0.772** −0.287** 蒸腾速率Trmmol −0.076 −0.217* 0.960** 0.169 0.971** 注:主对角线上为复相关系数,其左下方为简单相关系数,右上方为偏相关系数。
Note: on the main diagonal is the complex correlation coefficient, the lower left is the simple correlation coefficient, and the upper right is the partial correlation coefficient. -
分别建立危害时期与Photo、Cond、Ci及Trmmol的预测模型,拟合和预测精度的R2、RMSE选择精度最高的预测模型(表8):危害时期与Photo、Cond及Trmmol的三次函数模型,以及与Ci的线性函数模型。其中危害时期与Photo的三次函数模型的精度最高(拟合:R2=0.931,RMSE=0.204;预测:R2=0.961,RMSE=0.240)。
表 8 危害时期监测模型(n=63)及检验(n=27)
Table 8. Monitoring models(n = 63) and validation(n = 27) of damage time
监测模型Monitoring models 回归方程 Regression equations 拟合精度Fitted precision 预测精度Predicted precision R2 RMSE R2 RMSE 危害时期与Photo $y = 1.637 + 0.738{x_1} - 0.030{x_1}^2 - 0.004{x_1}^3$ 0.931 0.204 0.961 0.240 危害时期与Cond $y = 3.735 - 186.107{x_2} + 8\;349.646{x_2}^2 - 91\;980.434{x_2}^3$ 0.646 0.459 0.608 0.411 危害时期与Ci $y = 5.116 - 0.011{x_3}$ 0.609 0.488 0.632 0.560 危害时期与Trmmol $y = 10.319 - 20.166{x_4} + 15.464{x_4}^2 - 3.476{x_4}^3$ 0.854 0.295 0.871 0.361 危害时期与Photo、Cond、Ci、Trmmol $y = 2.358 - 0.259{x_1} - 17.174{x_2} + 0.003{x_3} + 0.092{x_4}$ 0.918 0.236 0.930 0.272 危害时期与Photo、Cond、Ci $y = 2.27 + 0.309{x_1} + 11.608{x_2} - 0.002{x_3}$ 0.916 0.225 0.911 0.261 注:y代表松小蠹的危害时期预测值,x1、x2、x3、x4分别代表Photo、Cond、Ci、Trmmol。
Note: y represents the predicted value of damage period of pine bark beetle, and x1, X2, X3 and X4 represent photo, Cond, Ci and Trmmol respectively.另根据危害时期与Photo、Cond、Ci及Trmmol的复相关分析结果,建立多元线性回归模型(拟合:R2=0.918,RMSE=0.236;预测:R2=0.930,RMSE=0.272),能较好的模拟预测松小蠹的危害时期。此外,还对Photo、Cond、Ci、Trmmol在建模分析中进行了重要性分析,前三个参数重要依次为0.707、0.247、0.047,而Trmmol在模型中的可以忽略,因此建立基于Photo、Cond、Ci的多元线性回归模型(拟合:R2=0.916,RMSE=0.225;预测:R2=0.911,RMSE=0.261)。
-
Hernandez[18]在针叶树木上研究发现PRI与气孔导度显著相关;NDVI与SR均能与Photo建立一定的相关性。结合本研究,对植被指数、高光谱特征参数与Photo、Cond、Ci及Trmmol进行相关性分析和逐步回归分析,筛选显著或极显著变量(表9),分别建立植被指数高光谱特征参数的拟合方程,并进行精度验证,根据拟合和预测精度的R2、RMSE选择精度最高的拟合方程(表10)。其中以参数NDVI705为变量的Photo拟合方程精度最高,拟合和预测R2均大于0.8。对Ci植被指数拟合方程的NDVI705、HI及高光谱特征参数拟合方程的λr、SDb、SDnir/SDb进行重要性分析,结果表明:NDVI705和HI的重要性分别为0.825、0.125;λr、SDb、SDnir/SDb的重要性依次为0.370、0.419、0.211。HI、SDnir/SDb虽与Ci未呈现显著相关,但该变量对其预测具有一定贡献,将该变量纳入方程中,提高方程拟合效果。
表 9 特征变量筛选结果
Table 9. The screening results of characteristic parameter
光合参数Photosynthetic parameters 特征变量Characteristic variables 相关分析Correlation analysis 逐步回归分析Stepwise regression analysis 净光合速率Photo NDVI705 0.899** <0.001** λr 0.832** <0.001** 气孔导度Cond SR2 0.773** <0.001** SDr 0.643** 0.002** SDr/SDb 0.731** 0.019* 胞间CO2浓度Ci NDVI705 -0.815** <0.001** HI -0.241 0.008** λr -0.767** 0.005** SDb 0.478* 0.003** SDnir/SDb 0.232 0.024* 蒸腾速率Trmmol SR1 0.715** <0.001** SDr 0.625** 0.011* SDr/SDb 0.660** 0.026* 表 10 光合参数估算方程(n=63)及验证(n=27)
Table 10. Estimation equation(n = 63) and validation(n = 27) of photosynthetic parameters
监测模型
Monitor models变量
Variable回归方程
Regression equations拟合精度
Fitted precision预测精度
Predicted precisionR2 RMSE R2 RMSE 净光合速率Photo NDVI705 $y = 0.770 - 2.041{x_1} + 66.972{x_1}^2 - 97.740{x_1}^3$ 0.822 0.426 0.863 0.539 λr $y = - 101.808 + 0.149{x_5}$ 0.692 0.704 0.674 0.720 气孔导度Cond SR2 $y = 0.023 - 0.005{x_3} + 0.001{x_3}^2 - 3.103{{\rm E}^{ - 5}} \cdot {x_3}^3$ 0.632 0.018 0.629 0.018 SDr $y = 0.010 + 0.339{x_6}$ 0.414 0.012 0.381 0.012 SDr/SDb $y = 0.012 + 0.003{x_7}$ 0.535 0.011 0.552 0.010 SDr、SDr/SDb $y = 0.005 + 0.202{x_6} + 0.002{x_7}$ 0.648 0.009 0.560 0.009 胞间CO2浓度Ci NDVI705 $y = 211.338 - 450.203{x_1} - 3\;444.328{x_1}^2 + 4\;506.702{x_1}^3$ 0.733 23.636 0.695 25.658 NDVI705、HI $y = 180.597 - 302.673{x_1} + 101.860{x_4}$ 0.766 22.148 0.672 27.904 λr $y = 2\;794.375 - 3.737{x_5}$ 0.588 29.347 0.500 31.132 λr、SDb、SDnir/SDb $y = 1\;731.709 - 2.309{x_5} + 4\;379.721{x_8} + 6.697{x_9}$ 0.745 23.112 0.693 29.808 蒸腾速率Trmmol SR1 $y = 0.784 - 0.061{x_2} + 0.026{x_2}^2 - 0.001{x_2}^3$ 0.544 0.476 0.604 0.467 SDr $y = 0.547 + 13.504{x_6}$ 0.391 0.503 0.385 0.495 SDr/SDb $y = 0.710 + 0.106{x_7}$ 0.435 0.485 0.436 0.450 SDr、SDr/SDb $y = 0.373 + 8.712{x_6} + 0.075{x_7}$ 0.561 0.427 0.564 0.415 注:y代表Photo、Cond、Ci、Trmmol的预测值,x1、x2、x3、x4、x5、x6、x7、x8、x9分别代表NDVI705、SR1、SR2、HI、λr、SDr、SDr/SDb、SDb、SDnir/SDb。
Note: y represents the predicted values of Photo, Cond, Ci and Trmmol, and x1, x2, x3, x4, x5, x6, x7, x8 and x9 represent NDVI 705、 SR1 、 SR2 、 HI 、 λr 、 SDr 、 SDr / SDb 、 SDb 、 SDnir / SDb, respectively.
基于高光谱特征和光合参数监测松小蠹不同危害时期的相关分析
Correlation Analysis of Different Hazard Periods of Tomicus spp. Based on Hyperspectral Features and Photosynthetic Parameters
-
摘要:
目的 测定云南松枝梢针叶高光谱数据及光合数据,建立松小蠹危害时期监测模型及光合参数模拟方程,有效、快速地预测松小蠹危害时期,诊断云南松健康状况,为大面积应用高光谱遥感技术实现森林病虫害监测提供依据。 方法 调查松小蠹危害情况,获取不同危害时期针叶高光谱及光合数据,选取与松小蠹危害时期显著相关的参数,建立危害时期监测模型;对针叶光合与高光谱特征参数进行相关分析,建立光合参数的相关方程。 结果 危害时间越久,在740~1 036 nm,云南松枝梢针叶光谱反射率逐渐降低;光谱一阶微分曲线在660~740 nm,出现“红边”“蓝移”现象,且峰值逐渐降低;基于光谱特征建立的松小蠹危害时期多元线性回归模型拟合效果最优;松小蠹危害时期与净光合速率Photo密切相关,其三次函数模型拟合效果较好;建立光谱指数和光合参数之间的联系,得到光合参数最优拟合方程。 结论 根据云南松光谱指数和光合参数分别建立的模型可有效监测松小蠹危害时期;云南松光谱指数与光合参数存在显著相关性,可建立相关模型预估云南松生长健康状况。 Abstract:Objective To establish monitoring models and photosynthetic parameter simulation equations by measuring the hyperspectral data and photosynthetic data of Pinus yunnanensis, so as to effectively and quickly predict the damage period of Tomicus spp. and diagnose the health status of P. yunnanensis, and provide a reference for large-area application of hyperspectral remote sensing technology to achieve forest pest monitoring. Method The damage caused by Tomicus spp. was investigated and the hyperspectral and photosynthetic data were obtained. The parameters which were significantly correlated with the damage period of Tomicus spp. were selected to establish monitoring models. The correlation between photosynthetic parameters and hyperspectral characteristic parameters was analyzed and the correlation equations of photosynthetic parameters were established. Result As the hazard period prolonged, the spectral reflectance gradually decreases in the range of 740-1036 nm. "Red edge" and "blue shift" appeared in the first-order differential curve of the spectrum at 660-740 nm, and the peak value gradually decreased. The multivariate linear regression model which was established based on spectral characteristics showed the best fitting effect. The damage period of Tomicus spp. was closely correlated with Photo, and its cubic function model fitted well. The correlation between spectral index and photosynthetic parameters was established to obtain the optimal fitting equation of photosynthetic parameters. Conclusion The models based on spectral index and photosynthetic parameters of P. yunnanensis can effectively monitor the hazard period of Tomicus spp. There is a significant correlation between the spectral index and photosynthetic parameters of P. yunnanensis, and a correlation model can be established to estimate the growth and health status of P. yunnanensis. -
Key words:
- hyperspectral
- / spectral characteristic parameter
- / Tomicus spp.
- / photosynthesis
- / correlation
-
表 1 松小蠹危害时期划分
Table 1. Hazard time division of Tomicus spp.
危害时期Hazard periods 特征描述Feature description 健康时期Health 云南松枝梢的针叶为绿色,未受到松小蠹的危害,枝梢上没有蛀孔 危害早期Early damage 松小蠹开始蛀食枝梢,枝梢出现蛀孔,但枝梢的针叶仍显绿色或嫩绿色 危害中期Mid-hazard 松小蠹仍停留在枝梢内,即将完成对枝梢的蛀食,枝梢的针叶颜色为嫩黄色或黄绿色 危害末期Terminal stage 松小蠹完成对枝梢的蛀食,针叶逐渐发生变化,枝梢的针叶颜色表现为红褐色,最后呈现灰褐色 表 2 特征变量表
Table 2. Characteristic variable scale
变量类别Variable categories 变量定义和公式Variable definitions and formulas 植被指数Vegetation indexes $NDVI1 = {{({{\bar R}_{(650\sim 670)}} - {{\bar R}_{(760\sim 850)}})} / {({{\bar R}_{(650\sim 670)}} + {{\bar R}_{(760\sim 850)}})}}$ $NDVI2 = {{({R_{838}} - {R_{677}})} / {({R_{838}} + {R_{677}})}}$ $NDV{I_{705}} = {{({R_{750}} - {R_{705}})} / {({R_{750}} + {R_{705}})}}$ $SR1 = {{{{\bar R}_{(650\sim 670}}} / {{{\bar R}_{(760\sim 850)}}}}$ $SR2 = {{{R_{838}}} / {{R_{677}}}}$ $WBI = {{{R_{900}}} / {{R_{970}}}}$,水波段指数 $PRI = {{({R_{531}} - {R_{570}})} / {({R_{531}} + {R_{570}})}}$,光化学植被指数 $HI = {{({R_{739}} - R{}_{399})} / {({R_{739}} + {R_{399}}) - 0.5{R_{404}}}}$,植被健康指数 高光谱特征参数
Hyperspectral Features反射率位置特征参数
Reflectance location parametersRg(绿峰),波长520~560 nm范围内最大的波段反射率 Rr(红谷),波长 640~680 nm范围内最小的波段反射率 $D = 1 - \dfrac{{{R_{670}}}}{{{R_{560}} + \dfrac{{{R_{760}} - {R_{560}}}}{{760 - 560}} \times (670 - 560)}}$ $H = 1 - \dfrac{{{R_{500}} + \dfrac{{{R_{670}} - {R_{500}}}}{{670 - 500}} \times (560 - 500)}}{{{R_{560}}}}$ 一阶微分位置特征参数
First derivative location parametersDb,蓝边(470~520 nm)内一阶微分光谱中的最大值 Dy,黄边(560~590 nm)内一阶微分光谱中的最大 Dr,红边(660~740 nm)内一阶微分光谱中的最大值 Dnir,近红外(760~1036 nm)内一阶微分光谱中的最大值 λr,红边(660~740 nm)内一阶微分光谱中的最大值对应的波段 一阶微分面积特征参数
First derivative area parametersSDb,蓝边内一阶微分值的总和 SDy,黄边内一阶微分值的总和 SDr,红边内一阶微分值的总和 SDnir,近红外范围内一阶微分值的总和 植被指数特征参数
Vegetation indexes parametersRg/Rr D/H SDr/SDb SDr/SDy SDnir/SDb SDnir/SDr (Rg-Rr)/(Rg+Rr) (D-H)/(D+H) (SDr-SDb)/(SDr+SDb) (SDrSDy)/(SDr+SDy) 表 3 植被指数变量和高光谱特征参数变量筛选结果
Table 3. The screening results of vegetation index and hyperspectral characteristic parameter
特征变量Characteristic variables 相关分析Correlation analysis 逐步回归分析Stepwise regression analysis(a) 特征变量Characteristic variables 相关分析Correlation analysis 逐步回归分析Stepwise regression analysis(b) NDVI2 0.918** <0.001** D 0.923** <0.001** PRI 0.740** 0.003** Dr 0.911** <0.001** Dr 0.911** 0.019* NDVI705 0.773** <0.001** D 0.923** 0.031* (SDr-SDb)/(SDr+SDb) 0.267 <0.001** SDb 0.381* 0.001** SDy -0.887** 0.005** 注:*: P<0.05,显著相关;**: P<0.01,极显著相关,下同。
Note: P < 0.05, significant correlation; * *: P < 0.01, extremely significant correlation, the same below.表 4 危害时期监测模型(n=84)及检验(n=36)
Table 4. Monitoring models(n = 84) and validation(n = 36) of damage time
变量类型
Variable categories监测模型
Monitor models回归方程
Regression equations拟合精度
Fitted precision预测精度
Predicted precisionR2 RMSE R2 RMSE 植被指数
Vegetation indexNDVI2 $y = 0.355 + 4.460{x_1} - 6.974{x_1}^2 + 7.410{x_1}^3$ 0.872 0.415 0.854 0.424 PRI $y = 3.637 + 13.108{x_2} - 31.072{x_2}^2 - 208.144{x_2}^3$ 0.623 0.711 0.636 0.791 NDVI2、PRI $y = 0.622 + 3.662{x_1} + 3.780{x_2}$ 0.878 0.391 0.889 0.357 高光谱特征参数
Hyperspectral characteristic
parametersD $y = 1.139 + 0.811{x_3} - 3.174{x_3}^2 + 6.235{x_3}^3$ 0.900 0.367 0.895 0.498 Dr $y = 0.755 + 83.599{x_4} + 27\;496.449{x_4}^2 - 1\;321\;854.674{x_4}^3$ 0.839 0.465 0.796 0.517 SDb $y = - 0.802 + 935.581{x_5} - 67\;441.244{x_5}^2 + 1\;450\;214.907{x_5}^2$ 0.441 0.867 0.106 1.010 SDy $y = 2.019 - 235.870{x_6} + 6\;602.191{x_6}^2 + 1\;404\;056.992{x_6}^3$ 0.828 0.481 0.804 0.508 D、Dr、SDb、SDy $y = 1.343 + 1.121{x_3} + 97.985{x_4} - 47.574{x_5} - 80.177{x_6}$ 0.916 0.314 0.918 0.331 综合Comprehensive D、Dr、NDVI705 $y = 0.540 + 1.107{x_3} + 108.960{x_4} + 2.545{x_7}$ 0.937 0.276 0.925 0.258 D、Dr、NDVI705、(SDr−SDb)/(SDr+SDb) $y = 0.886 + 1.332{x_3} + 86.860{x_4} + 3.636{x_7} - 0.794{x_8}$ 0.969 0.189 0.946 0.212 注:y代表松小蠹的危害时期预测值,x1、x2、x3、x4、x5、x6、x7、x8分别代表NDVI2、PRI、D、Dr、SDb、SDy、NDVI705、(SDr−SDb)/(SDr+SDb)。
Note: y represents the predicted value of damage period, x1, x2, x3, x4, x5, x6, x7 and x8 represent NDVI2, PRI, D, Dr, SDb, SDy, NDVI705, (SDr − SDb) / (SDr + SDb), respectively.表 5 云南松受松小蠹危害的不同时期的针叶气体交换参数的影响
Table 5. The effects of gas exchange parameters in different periods of Pinus yunnanensis damaged by Tomicus spp.
主要光合指标Principal photosynthetic indexs 云南松受松小蠹危害的不同时期Different periods of harm 健康时期Health 危害早期Early damage 危害中期Mid-hazard 净光合速率Photo/(μmol·m−2s−1) 4.621 ± 0.105a 2.089 ± 0.087b 0.643 ± 0.048c 气孔导度Cond/(mol·m−2s−1) 0.047 ± 0.001a 0.023 ± 0.002b 0.016 ± 0.001c 胞间CO2浓度Ci/(μmol·mol−1) 129.098 ± 19.591a 174.854 ± 48.175b 235.635 ± 26.579c 蒸腾速率Trmmol/(mmol·m−2s−1) 1.938 ± 0.014a 1.053 ± 0.111b 0.889 ± 0.034b 水分利用效率WUE/(μmol·mmol−1) 2.385 ± 0.053a 2.631 ± 0.201a 0.755 ± 0.058b 气孔限制值Ls 0.578 ± 0.011a 0.425 ± 0.026b 0.277 ± 0.013c 注:表中数据为测得数据的平均值±标准误差,同行中不同字母表示不同处理间Duncan’s多重比较差异显著(P<0.05,n =90)。
Note: the data in the table is the mean ± standard error of the measured data. Different letters in the same line indicate that there are significant differences in Duncan's multiple comparisons between different treatments (P < 0.05, n = 90).表 6 针叶主要光合参数的相关性分析
Table 6. The correlation analysis of main photosynthetic parameters of needles
危害时期Hazard periods 因子Factors 净光合速率Photo 气孔导度Cond 胞间CO2浓度Ci 蒸腾速率Trmmol 健康时期Health 净光合速率Photo 1 0.708** −0.989** 0.383 气孔导度Cond −0.410* 1 0.623** 0.371 胞间CO2浓度Ci −0.732** 0.315 1 0.486** 蒸腾速率Trmmol 0.218 0.967** 0.497** 1 危害早期Early damage 净光合速率Photo 1 0.309 −0.476** −0.320 气孔导度Cond −0.520** 1 −0.278 0.999** 胞间CO2浓度Ci −0.730** 0.528** 1 0.291 蒸腾速率Trmmol −0.543** 0.999** 0.551** 1 危害中期Mid-hazard 净光合速率Photo 1 −0.360 −0.508** 0.702** 气孔导度Cond 0.205 1 −0.519** 0.663** 胞间CO2浓度Ci 0.069 0.124 1 0.691** 蒸腾速率Trmmol 0.579** 0.512** 0.401 1 注:主对角线左下方为简单相关系数,右上方为偏相关系数。 表 7 松小蠹危害时期和针叶主要光合参数的相关性分析
Table 7. The correlation between damage time of Tomicus spp.andmain photosynthetic parameters
因子Factors 危害时期Hazard periods 净光合速率Photo 气孔导度Cond 胞间CO2浓度Ci 蒸腾速率Trmmol 危害时期Hazard periods 0.909** 0.943** 0.797** −0.764** 0.688** 净光合速率Photo 0.588** 0.926** 0.758** −0.821** 0.635** 气孔导度Cond 0.187 0.293** 0.979** −0.433** 0.972** 胞间CO2浓度Ci −0.155 −0.484** −0.016 0.772** −0.287** 蒸腾速率Trmmol −0.076 −0.217* 0.960** 0.169 0.971** 注:主对角线上为复相关系数,其左下方为简单相关系数,右上方为偏相关系数。
Note: on the main diagonal is the complex correlation coefficient, the lower left is the simple correlation coefficient, and the upper right is the partial correlation coefficient.表 8 危害时期监测模型(n=63)及检验(n=27)
Table 8. Monitoring models(n = 63) and validation(n = 27) of damage time
监测模型Monitoring models 回归方程 Regression equations 拟合精度Fitted precision 预测精度Predicted precision R2 RMSE R2 RMSE 危害时期与Photo $y = 1.637 + 0.738{x_1} - 0.030{x_1}^2 - 0.004{x_1}^3$ 0.931 0.204 0.961 0.240 危害时期与Cond $y = 3.735 - 186.107{x_2} + 8\;349.646{x_2}^2 - 91\;980.434{x_2}^3$ 0.646 0.459 0.608 0.411 危害时期与Ci $y = 5.116 - 0.011{x_3}$ 0.609 0.488 0.632 0.560 危害时期与Trmmol $y = 10.319 - 20.166{x_4} + 15.464{x_4}^2 - 3.476{x_4}^3$ 0.854 0.295 0.871 0.361 危害时期与Photo、Cond、Ci、Trmmol $y = 2.358 - 0.259{x_1} - 17.174{x_2} + 0.003{x_3} + 0.092{x_4}$ 0.918 0.236 0.930 0.272 危害时期与Photo、Cond、Ci $y = 2.27 + 0.309{x_1} + 11.608{x_2} - 0.002{x_3}$ 0.916 0.225 0.911 0.261 注:y代表松小蠹的危害时期预测值,x1、x2、x3、x4分别代表Photo、Cond、Ci、Trmmol。
Note: y represents the predicted value of damage period of pine bark beetle, and x1, X2, X3 and X4 represent photo, Cond, Ci and Trmmol respectively.表 9 特征变量筛选结果
Table 9. The screening results of characteristic parameter
光合参数Photosynthetic parameters 特征变量Characteristic variables 相关分析Correlation analysis 逐步回归分析Stepwise regression analysis 净光合速率Photo NDVI705 0.899** <0.001** λr 0.832** <0.001** 气孔导度Cond SR2 0.773** <0.001** SDr 0.643** 0.002** SDr/SDb 0.731** 0.019* 胞间CO2浓度Ci NDVI705 -0.815** <0.001** HI -0.241 0.008** λr -0.767** 0.005** SDb 0.478* 0.003** SDnir/SDb 0.232 0.024* 蒸腾速率Trmmol SR1 0.715** <0.001** SDr 0.625** 0.011* SDr/SDb 0.660** 0.026* 表 10 光合参数估算方程(n=63)及验证(n=27)
Table 10. Estimation equation(n = 63) and validation(n = 27) of photosynthetic parameters
监测模型
Monitor models变量
Variable回归方程
Regression equations拟合精度
Fitted precision预测精度
Predicted precisionR2 RMSE R2 RMSE 净光合速率Photo NDVI705 $y = 0.770 - 2.041{x_1} + 66.972{x_1}^2 - 97.740{x_1}^3$ 0.822 0.426 0.863 0.539 λr $y = - 101.808 + 0.149{x_5}$ 0.692 0.704 0.674 0.720 气孔导度Cond SR2 $y = 0.023 - 0.005{x_3} + 0.001{x_3}^2 - 3.103{{\rm E}^{ - 5}} \cdot {x_3}^3$ 0.632 0.018 0.629 0.018 SDr $y = 0.010 + 0.339{x_6}$ 0.414 0.012 0.381 0.012 SDr/SDb $y = 0.012 + 0.003{x_7}$ 0.535 0.011 0.552 0.010 SDr、SDr/SDb $y = 0.005 + 0.202{x_6} + 0.002{x_7}$ 0.648 0.009 0.560 0.009 胞间CO2浓度Ci NDVI705 $y = 211.338 - 450.203{x_1} - 3\;444.328{x_1}^2 + 4\;506.702{x_1}^3$ 0.733 23.636 0.695 25.658 NDVI705、HI $y = 180.597 - 302.673{x_1} + 101.860{x_4}$ 0.766 22.148 0.672 27.904 λr $y = 2\;794.375 - 3.737{x_5}$ 0.588 29.347 0.500 31.132 λr、SDb、SDnir/SDb $y = 1\;731.709 - 2.309{x_5} + 4\;379.721{x_8} + 6.697{x_9}$ 0.745 23.112 0.693 29.808 蒸腾速率Trmmol SR1 $y = 0.784 - 0.061{x_2} + 0.026{x_2}^2 - 0.001{x_2}^3$ 0.544 0.476 0.604 0.467 SDr $y = 0.547 + 13.504{x_6}$ 0.391 0.503 0.385 0.495 SDr/SDb $y = 0.710 + 0.106{x_7}$ 0.435 0.485 0.436 0.450 SDr、SDr/SDb $y = 0.373 + 8.712{x_6} + 0.075{x_7}$ 0.561 0.427 0.564 0.415 注:y代表Photo、Cond、Ci、Trmmol的预测值,x1、x2、x3、x4、x5、x6、x7、x8、x9分别代表NDVI705、SR1、SR2、HI、λr、SDr、SDr/SDb、SDb、SDnir/SDb。
Note: y represents the predicted values of Photo, Cond, Ci and Trmmol, and x1, x2, x3, x4, x5, x6, x7, x8 and x9 represent NDVI 705、 SR1 、 SR2 、 HI 、 λr 、 SDr 、 SDr / SDb 、 SDb 、 SDnir / SDb, respectively. -
[1] 邓喜庆, 皇宝林, 温庆忠, 等. 云南松林在云南的分布研究[J]. 云南大学学报: 自然科学版, 2013, 35(6):843-848. [2] 段兆尧, 杨子祥, 王健敏, 等. 云南切梢小蠹的室内饲养及隐蔽期生活习性的观察方法[J]. 林业科学研究, 2013, 26(3):389-392. doi: 10.3969/j.issn.1001-1498.2013.03.020 [3] 李 燕. 云南玉溪松小蠹虫成灾原因及防治对策[J]. 林业资源管理, 2013(3):40-42. doi: 10.3969/j.issn.1002-6622.2013.03.010 [4] 刘 畅, 孙鹏森, 刘世荣. 基于反射光谱指数的锐齿栎叶片色素含量估算[J]. 林业科学研究, 2017, 30(1):88-98. [5] Pan J, Ju Y, Zhang H, et al. Detection of Bursaphelenchus xylophilus infection in Pinus massoniana from hyperspectral data[J]. Nematology, 2014, 16(10): 1197-1207. doi: 10.1163/15685411-00002846 [6] 王 震, 张晓丽, 安树杰. 松材线虫病危害的马尾松林木光谱特征分析[J]. 遥感技术与应用, 2007, 22(3):367-370. doi: 10.3969/j.issn.1004-0323.2007.03.012 [7] 伍 南, 刘君昂, 闫瑞坤, 等. 马尾松赤枯病冠层光谱特征及严重度反演[J]. 中国农学通报, 2012, 28(4):51-57. doi: 10.3969/j.issn.1000-6850.2012.04.011 [8] 徐 爽, 易 东. 利用高光谱成像技术检测长枣表面虫伤[J]. 电子制作, 2013,(21):47-48. [9] Helmi Zulhaidi M S, Hamdan N. Hyperspectral imagery for mapping disease infection in oil palm plantation using vegetation indices and red edge techniques[J]. American Journal of Applied Sciences, 2009, 6(6): 1031-1035. doi: 10.3844/ajassp.2009.1031.1035 [10] 于泉洲, 周 蕾, 王绍强, 等. 基于EO-1 Hyperion的中国典型森林冠层高光谱特征分析[J]. 云南大学学报: 自然科学版, 2018, 40(5):125-132. [11] 宋小宁, 马建威, 李小涛, 等. 基于Hyperion高光谱数据的植被冠层含水量反演[J]. 光谱学与光谱分析, 2013, 33(10):2833-2837. doi: 10.3964/j.issn.1000-0593(2013)10-2833-05 [12] 刘良云, 黄木易, 黄文江, 等. 利用多时相的高光谱航空图像监测冬小麦条锈病[J]. 遥感学报, 2004, 8(3):275-281. doi: 10.11834/jrs.20040313 [13] 季 梅, 刘宏屏, 董谢琼, 等. 林木受小蠹虫危害的表征及其遥感监测[J]. 西部林业科学, 2016, 35(4):150-153. [14] 马华文, 许翠清, 李海朝. 立地条件对长白落叶松光合特性的影响[J]. 东北林业大学学报, 2008, 36(8):4-7. doi: 10.3969/j.issn.1000-5382.2008.08.002 [15] 刘 琦, 崔世茂, 宋 阳, 等. NaCl胁迫对番茄苗期根形态和光合作用的影响[J]. 北方农业学报, 2018, 46(3):32-36. doi: 10.3969/j.issn.2096-1197.2018.03.07 [16] Savitzky A, Golay M. Smoothing and Differentiation of Data by Simplified Least Squares Procedures[J]. Analytical Chemistry, 1964, 36(8): 1627-1639. doi: 10.1021/ac60214a047 [17] 汪 红, 石 雷, 马云强, 等. 云南切梢小蠹危害云南松监测模型与判定规则[J]. 林业科学研究, 2018, 31(4):53-60. [18] Rocío Hernández-Clemente, Navarro-Cerrillo R M, Lola Suárez, et al. Assessing structural effects on PRI for stress detection in conifer forests[J]. Remote Sensing of Environment, 2011, 115(9): 2360-2375. doi: 10.1016/j.rse.2011.04.036 [19] Gamon J A, Surfus J S. Assessing leaf pigment content and activity with a reflectometer[J]. New Phytologist, 2010, 143(1): 105-117. [20] Gamon J A, Pe Uelas J, Field C B. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency[J]. Remote Sensing of Environment, 1992, 41(1): 35-44. doi: 10.1016/0034-4257(92)90059-S [21] Merzlyak M N, Chivkunova O B, Gitelson A A, et al. Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening[J]. Physiologia Plantarum, 2010, 106(1): 135-141. [22] 王晓庆, 冉 烈, 彭 萍, 等. 炭疽病胁迫下的茶树叶片高光谱特征分析[J]. 植物保护, 2014, 40(6):13-17. doi: 10.3969/j.issn.0529-1542.2014.06.003 [23] 伍 南, 刘君昂, 周国英, 等. 利用高光谱微分指数反演油茶炭疽病病情指数的研究[J]. 中国生态农业学报, 2012, 20(6):777-781.