基于植被指数的秦岭森林健康评价方法研究
2022-04-06
来源:好走旅游网
西北林学院学报2013,28(6):145~150 Journal of Northwest Forestry University doi:10.3969/j.issn.1001—7461.2013.06.29 基于植被指数的秦岭森林健康评价方法研究 麻坤,赵鹏祥 季菲菲,翟晓江,孙 凯 (西北农林科技大学林学院,陕西杨陵712lOO) 摘要:利用HJ一1/HIS L2级高光谱数据提取敏感反应森林健康状况的植被指数,用敏感性分析法 筛选指标,建立评价森林健康的指标体系与评价模型。基于层次分析法获得评价模型中各指标的权 重。采用自然间断点分级法将森林健康指数分为4个等级,分别为不健康、亚健康、健康和优质健康, 绘制健康状况分布图。从小班和样地二水平验证该森林健康评价模型的精度,结果表明,评价结果在 二水平上均与传统调查方法获得的结果显著相关。决定系数分别为0.707和0.665。该模型所得森 林健康指数小于传统调查结果,但趋势一致。 关键词:植被指数;高光谱影像;敏感性分析;层次分析法 中图分类号:¥757.1 文献标志码:A 文章编号:1001—7461(2013)06—0145—06 Forest Health Assessment Method in Qinling Mountains Based on the Vegetation Indices MA Kan,ZHA0 Peng-xiang ,jI Fei_fei,ZHAI Xiao-jiang,SUN Kai Northwest A&F University,gangLing,Shaanxi 712100,China) (College of Forestry, Abstract:Based on the HJ一1 A L2 hyperspectral data,this paper extracted vegetation indexes which were sensitive to the forest health and screened the reasonable evaluation indexes by means of sensitivity analy— sis.The model to evaluate the health status of the sample area was established.We worked out the weights of the indexes by analytic hierarchy process and reclassified the FHI into four levels(unhealthy, sub—healthy,healthy,good—healthy)based on the natural discontinuity classification method.The distri— bution map of the forest health leve1 was drawn;the precision of the model on the stand leve1 and the sam— ple plot leve1 was verified.The results demonstrated that there was significant correlation between the re— sult we obtained from the model and the traditional investigation method both on the two levels.The de— termination coefficients reached 0.707 and 0.655,respectively.The FHI of the model was below which of the traditional investigation method.Nevertheless,the trend was the same. Key words:vegetation index;hyperspectral imaging;sensitivity analysis;analytic hierarchy process 植被对电磁波的响应,即植被的光谱反射或发 射特性是由其化学和形态学特征决定的,这种特征 与植被的发育,健康状况以及生长条件密切相关 J。 高光谱遥感是通过对不同类型植被的生物物理化学 据的将敏感反应林木健康信息的植被指数作为指标 体系,来评价森林健康的模型和软件。而我国目前 对于森林健康的研究仅局限于传统的样地调查方 法,评价周期长、尺度小、工作量大,且对森林健康状 况的变化无法及时有效地反映。另一方面,我国高 光谱遥感在农业和林业领域的应用已取得了长足的 进步,主要集中在叶面积指数反演l_4]、植被的生理生 成分含量的估算,可获得较为精细的植被生态学信 息[2],尤其对红边位移的监测可获得有关植被受环 境胁迫的信息[3]。因此,敏感表示此类信息的植被 指数可以从林木的生理生化水平评价和监测森林健 康状况。目前,国外已经开发出基于高光谱遥感数 收稿日期:2013—04一10 修回日期:2013-04—28 基金项目:林业公益性行业科研专项(200804022B)。 化参量反演 和植被精细分类_6 等研究领域,在森 林健康评价方面的研究相对较少。本研究选取与林 作者简介:麻坤,男,在读硕士。研究方向:3s技术在资源与环境中的应用。E—mail:conanmakun2036@163.corn 通信作者:赵鹏祥,男,副教授,硕士生导师。研究方向:林业遥感、地理信息系统应用。E—mail:zpx666@263.net 146 西北林学院学报 28卷 木健康状况密切相关的植被指数作为指标体系,采 用层次分析法确定指标权重以建立评价森林健康的 模型。并基于高光谱影像,通过波段运算得到森林 健康状况分布图,以确立一种大尺度评价森林健康 状况的方法。 1 研究区概况与数据处理 1.1研究区概况 研究区域位于秦岭南坡中段国家生态观测研究 站火地塘林区内,该地区属于中国北亚热带和暖温 带的过渡地带。年平均气温8~10℃,年降水量 1 000~1 200 mm,海拔1 420~2 470 m。中山地带 在20世纪60~7O年代进行过全面的森林主伐,现 森林植被恢复较好,全区天然次生林面积较大,分布 集中,在部分地段有一定面积人工林,森林总覆盖率 达93.8 。主要的植被类型包括寒温带针叶林、针 阔叶混交林和落叶阔叶林。成林树种以华山松 (Pinus armandii)、油松(Pinus tabulaeformis)、锐 齿栎(Quercus aliena var.acuteserrata)、红桦(Bet— ula albo—sinensis)等为主 ]。 1.2数据获取及预处理 采用国产环境与灾害监测预报小卫星上装载的 超光谱成像仪于2011年4月13日拍摄的~幅HJ一 1/HIS L2级高光谱数据。影像数据覆盖范围为 107。57 3O ~1O8。37 OO E,33。17 。0 ~33。49 ()() N 的区域。包含了火地塘林区在内的秦岭中段森林。 影像空间分辨率为i00 m,重访周期为96 h,具有 115个波段,光谱范围0.45~0.95 m,平均光谱分 辨率为4.32 nm。L2级HSI数据经过了辐射校正 和几何粗校正,完成了由影像DN值向绝对辐射亮 度值转换和坏线去除等前期工作。为了达到使用要 求,应做的后续处理包括降噪 ]、大气校正和几何精 校正。基本流程如图1所示。 1.3植被指数提取 提取能够敏感反映植被生长健康状况的植被指 数包含如下三类:绿度指数、光利用率指数和叶色素 指数(表1)。 1.3.1 绿度指数 主要包含宽带绿度指数和窄带 绿度指数,前者主要用来简单度量绿色植被的数量 和分布情况。它对植物的叶绿素含量和叶子表面冠 层结构比较敏感。后者主要基于绿色植物的“红边” 效应(健康植被的光谱响应在红外和近红外波段区 域陡然增加)描述植物的健康及色素状态。 1.3.2 光利用率指数 主要通过度量植被在光合 作用中对入射光的利用效率来评价碳的吸收效率、 植被的生长速度和光合有效辐射(FAPAR)。据此 可从森林的生长速度方面评估森林的健康状况。 1.3.3 叶色素指数 主要是通过度量植被中与胁 迫性相关的色素来监测森林的衰减状况,如类胡萝 卜素和花青素,这类色素大量存在于衰减植被中。 该类色素含量多,证明植物目前处于衰老、病虫害以 及干旱等影响之下。因此,该类指数和森林健康状 况负相关,也可作为评价森林生长健康状况的重要 指标。 数据格式转换 HDF5 成G 1 条纹去除 全局去 1 r 大气校正 FLAA 1 r 几何精校亚 图1 HSI影像处理流程图 Fig.1 Flow diagram of the HIS data processing 2 评价方法 2.1指标筛选 初选的植被指数可以很好地反映森林的健康状 况,但植被指数在不同立地条件下对植被生长状况 的敏感性不同,为避免信息冗余,采用许明祥Ll 等 2005年提出的敏感性分析法,剔除在该地区对植被 生长状况反应不明显的植被指数(表2)。 敏感性分析法的原理主要是计算影像内各个植 被指数的变异系数,变异系数越大,表示该植被指数 在研究区域内可以敏感反应植被生长状况;变异系 数越小,表示该植被指数对植被生长情况的反应不 敏感,应剔除。计算后的各植被指数的变异系数情 况如表2,根据敏感性划分依据(表3),发现Vo— gelmann红边指数1和红边位置指数对该地区森林 生长情况的反应不敏感,故剔除此2指标。 2.2评价模型的建立 根据筛选出的11个植被指数指标,从表示植被 分布情况的绿度指数(g),表示植被生长状况的光 利用率指数(z)和表示植被衰减和胁迫情况的叶色 素指数(P)三方面建立评价森林生长健康状况 (FHD的模型(公式1)。 第6期 麻坤等:基于植被指数的秦岭森林健康评价方法研究 147 植被指数 最小值 最大值 平均值 标准差 变异系数 排序 注:为了剔除噪声对实验结果的干扰,v了…取VI原值在统计直方图上累计统计个数达o.05 处的值, J…取整幅影像VJ原值统计直方图上 累计统计个数在99.95 处的值。 FHI一(∑"OiW )×W +(∑ 叫,)×W + (∑ )×W (1) alytic herarchy process,AHP)来确定。该方法先把 问题层次化,根据问题的性质和目标,将其分解为不 同组成因素,按照因素间的相互关联影响以及隶属 式中, , ,, ,分别表示三方面涉及的植被指 数,w , ,w 分别表示三方面的权重,而叫 ,训 , 叫 表示每个方面涉及的植被指数所占的权重。且 各权重值应满足公式2。 关系将因素不同层次聚集组合,形成一个多层次的 分析结构模型。本研究中,将森林健康状况问题分 为3个层次(图2)。 指标权重的确定对于评价结果的客观性与合理 性具有重要影响。本研究采用层次分析法并结合德 尔菲法和该区域内植被指数敏感性分析结果,邀请 (2) J 1 k=I + + ===l∑ 一1∑ 一1∑ 一1 z 1 了数十位林学、生态学、植物生理学等领域的专家对 表3感性分析表 Table 3 Sensitivity analysis table 各项指标进行一一对比评判与量化运算。确定了 1~5及其倒数10个判断标度,最终得到了各层次指 标权重的判断矩阵(表4)。 计算指标的权重可归纳为求判断矩阵的特征根 2.3指标权重的确定 各植被指数指标的权重采用美国运筹学家T. L.Saaty 20世纪70年代提出的层次分析法(an 和特征向量问题。对于判断矩阵A,如果满足AW一 MAxw,AMAX为矩阵的最大特征值,w为对应于 M 的 特征向量,则w的各个分量即各指标的权重值。 148 西北林学院学报 28卷 准 } 子 准 则 层 2 2 图2森林健康评价指标层次 Fig.2 Indicator system of the FHI evaluation model 2.4一致性检验 由于判断矩阵是采用专家打分法和敏感性分析 CI=0时,矩阵具有完全一致性;cI≠0时,求其与 同阶的随机一致性指标RI之比CR,如CR<O.10, 法结合而得出,存在人为因素干扰,因此,必须对结 果进行一致性检验以判断矩阵是否合理。 2.4.1 层次单排序 层次单排序是指根据判断矩 则认为判断矩阵具有很好的一致性。检验(表4), 可见所有判断矩阵均具有很好的一致性,结果合理。 2.4.2层次总排序 利用同一层次中所有层次单 阵计算对于上一层某因素而言,本层次内与之有联 系的因素的重要性次序。按照本层次各组分内因素 相对上一层次而言的重要性进行排序的基础。为了 检验矩阵的一致性,需计算一致性指标cI。 一排序的结果,就可以计算针对上一层次而言,本层次 所有组分中因素重要性的权值,这就是层次 9,al ̄序。 其组合权重就是评价模型中各指标的权重。cI和 RI通过求上层指标的加权平均值获得(表5)。结 果显示,指标权重合理。 ∑—(AW—)I,CI一 Table 4 (3) 表4层次单排序计算结果 The results of single—level sequencing 判断矩阵 权重 一致性检验 CI RI B4 2 B5 3 1/3 1/5 1/2 1 2 1/3 1/2 1 B11 1/2 2 1/2 1 2.5数据标准化 标准化方法对数据进行处理(公式4)。对于与森林 由于各植被指数指标计算公式的差异,其结果 通常不在同一数量级上,因此本研究采用min—max 健康状况负相关的植被指数(主要指叶色素指数), 取其结果的相反数为实际输出结果。 第6期 麻坤等:基于植被指数的秦岭森林健康评价方法研究 l49 w一± 暑 3 结果与分析 3,1森林健康分级结果 脚 粥矾 邯 邯鼬眦 0.72之间。本研究采用Jenks提出的自然间断 ~㈩ 示,整幅遥感影像的森林健康指数(FH工)范围在0 点分级法,将FHJ值分为不健康、亚健康、健康和 优质健康4个等级。该方法主要根据统计学中的最 优化观点得出各等级的分界点,能够使各等级FH 值的内部方差之和最小。分类结果和等级划分依据 如图3。 据上述森林健康评价模型,基于ENVI软件中 的波段运算模块,对影像逐象元进行处理。结果显 表5层次总排序计算结果 Table 5 The results of general sorting 准则层A 子准则层B 组合权重 排序 一致性检验 RI一兰a RI,一0.582 5 J一1 CR一丽CI一0090 2<.1 果进行对比(图4)。结果表明,基于植被指数的评 幂●tJt 王.t▲0-0.勰 0.18-0.26 价模型得出的森林健康指数略低于传统的基于样地 和小班数据的评价结果。首先,因为高光谱影像受 大气校正、尺度转换、土壤背景等因素影响较大,提 取的植被指数有不同程度的衰减;其次,植被指数评 价模型中选取的部分指标标准化后取值在(一1,1) 区间,而传统模型指标标准化后取值在(0,1)区间。 冀It 甓曩40.0.26-0.13 033-0. ̄ | l ・, N 但曲线趋势一致,说明该模型可以理想地反映森林 的健康程度。因此,对健康等级的划分结果是合理 的。线性回归结果表明两种水平的评价结果和本研 究得出的评价结果显著相关。判定系数R2分别为 0.707和0.665, <0.000 1。 健康等扳 一■■■■a ̄.1q.,llt uR-Ml1. I健曩.-优质健纛 一…,,i I8千米 ’ 4 结语 随着高光谱遥感的迅速发展,可以从高光谱数 据中提取的植被信息越来越详尽,本研究从林木生 .. l 图3森林健康等级分布图 Fig.3 The distribution map of forest health level 3.2与传统调查方法的比较 为了验证评价模型的合理性,选取影像区域内 20个小班数据和2O个样地数据从小班和样地2水 理生化水平选取了敏感反应森林健康状况的11个 植被指数,建立评价模型,对实现秦岭地区森林健康 状况快速有效的评估具有意义。可为大尺度的森林 健康评价和监测研究提供一种新的思路。但由于国 平与传统的样地调查结果进行比较。在保积存口 等(2010)提出的评价方法基础上进行优化得到基于 产高光谱数据的一些缺陷,本模型依然存在问题,如 空间分辨率较低,对于小尺度林分变化反映不灵敏; 影像数据无法有效反映林分结构和生物多样性,无 样地数据和小班数据评价森林健康程度的模型。由 于基于植被指数对森林健康的评价主要考虑林木健 康,而非森林生态系统健康,故只选择样地和小班数 据中能反应林木健康的指标进行评价,并对评价结 法从生态系统尺度去实现森林生态系统健康状况的 有效评价等。这类问题将随着高光谱遥感科学的进 一步发展和植被指数研究的完善会得到有效解决。 15O 两北林学院学报 28卷 O.90 0.9O 0.8O 0.7O 0.40 O.20 O.1O 0 0 5 1O 15 20 0 15 0.2O 0.25 0-30 0.35 0.40 0.45 验证小班号 1.O 实测值 验证点一回归曲线 叫~植被指数法一小班水平 靼0.7 髅8音 O 5 10 15 20 0.1 0.15 0.20 0.25 0.30 0.35 0.40 0.45 验证小班号 实测值 +植被指数法一样地水平 ・验证点——回归曲线 图4两种水平上对模型的验证 Fig.4 The model verification on two levels EO 1 Hyperion hyperspectral data[J].Remote Sensing Infor 参考文献: Eli BOOCHS F,KUPFER G,DOCKTER K,e£a1.shaDe。f0OO0OOOOO the red edge as vitality indicator fnr plants r¨. International 蛐加∞如∞如加 O Journal of Remote Sensing,1 990,11(10):1741 1 753. 器聪 mation, oo (0):36 41.(i“chi“。 ) [9]SELLERS P J・Canopy reflectance,photosynthesis and tran— 印i a i。n_r・I“ter“ationa1 Journal of Re mf】te Sensing・1985, 6(8):1335 1372・ r2] JOHNS()N I F,HI AVKA C A,PETERSON D I .Multiva riate analysis 0f AVRIS data f。r can。Dv bi。chemical estima— [1O] HUETE A R,I.IU H Q,BATCHII Y K, £al・A compari— son of vegetation indices over a global set of TM images for ti。n a1ong the Oregon transectEJ].Remote Sensing 0f Envi r0nment,l994,47(2):21 6 232. EOS—MODIS[J]・Remote Sensing of EnVironment,1 997,59 L 1:44u一4bl‘ [33 MII I ER J R,HARE E W,wu J.Quantitative characteriza一 ti0n of the vegetation red edge reflectance 1. an inverted gaussian reflectance model r J]. International Journal of Re一 [11] L()BEI I D B,G P ASNER・Hyperion studies of crop stress in Mexico[C]・12 JP[ AVIRIS Geoscience Workshop.Pas— d n ‘LA ZOO3. m0te Sensing,1990,11(1O):1755—1773. [12]SIMS D A,GAMON J A.Relationships between leaf pig [4]赵丽芳,谭炳香,杨华,等.高光谱遥感森林叶面积指数估测研 究现状[J].世界林业研究,2007,20(2):50—54. ZHAO I F,TAN B X,YANG H,etd .Review ofthe estimating m content and。P ⅡaI retlectance across wide “g ot species,leaf structures and developmental stages[J]・Remote Sensing of Environment,2002,81(2—3):337—354. forest 1eaf area index with hyperspectra1 remote sensing data[J]. World Forestrv Research,2007,20(2):50—54.(in Chinese) [13]GAMON J A,SERRANO L,SURFUS J S.The photochem— cal reilectance index:an optical indicator ot photosynthetic [5]刘秀英,熊建利,臧卓,等.基于植被指数的马尾松叶绿素含量 估算模型[J].西北林学院学报,2012,27(3):44—47. I Iu x Y,XIONG J I ,ZANG Z,etⅡ . Estimation models f。r chlorophyll c。ntent。f P “5 mn s0nin d based on vegeta— ble indexes rI]. Journal of Northwest Forestrv Universitv, radiation use fl“ Y across P ,luncttonal tYP and nutrient levels[J]・Oecologia 1997,112(4):492—501. [1 4] GAMON J A,SURFUS J S.Assessing leaf pigment content and activity with a reflectometer[J].New Phytologist,1999, 143(1):105—117. 201 2,27(3):44—47.(in Chinese) [6] 吴见,彭道黎.基于空间信息的高光谱遥感植被分类技术[J]. 农业工程学报,2012,28(5):1 50 1 53. WU J,PENG D I . Vegetati。n classification techno1。gy。f [151 GITEI SON A A,ZUR Y,CH1VKUNOVA O B,et a1.As sesslng carotenoid content ln plant leaves with retleetance spectroscopy[J].Photochemistry and Photobiology,2002,75 (3):272—281. hyperspectral remote sensing based on spatial information[J]. Transactions of the Chinese Society of Agricultural Engineer— [16]GITEI SON A A,MERZI YAK M N・CHIVKUNOVA O B.Optical properties and nondestructive estimation of antho ing,2012,28(5):150 153.(in Chinese) [7]保积存,王得祥,陈飞,等.基于灰色关联的火地塘林区锐齿栎 林天然更新研究[J].西北林学院学报,2011,26(5):121 126. BAO I C,WANG D X,CHEN F,et nz.Naturalregeneration of Quercu ̄aliena var.acuteserrata forests in Huoditang cyanin content in plant leaves[J].Photochemistry and Photo一 biology,2001,74(1):38 45. [17]许明祥,刘国彬,赵允格.黄土丘陵区土壤质量评价指标研究 [J].应用生态学报,2005,16(1o):1843 1848. XU M X,LIU G B,ZHAO Y( .Assessment indicators of based on the grey relational analysis method[J].Journal of Northwest Forestry University,2011,26(5):121—126.(in soil quality in hilly Loess Plateau[J].Chinese Journal of Ap— plied Ecology,2005,1 6(10):1 843—1848.(in Chinese) Chinese) [8]谭炳香,李增元,陈尔学,EO-1Hyperion高光谱数据的预处理 [J].遥感信息,2005(6):36—41. TAN B X,g1 Z Y,CHEN E X,PANG Y.Preprocessing of [18] 托马斯・萨迪.领导者:面临挑战与选择 层次分析法在决 策中的应用[M].张录,译北京:中国经济出版社,1993. [19]保积存.秦岭火地塘林区锐齿栎林健康状况及其动态分析 (D).陕西杨陵:西北农林科技大学.2O11.