文章信息
- 王修信, 孙涛, 朱启疆, 刘馨, 高凤飞, 胡玉梅, 陈声海
- WANG Xiuxin, SUN Tao, ZHU Qijiang, LIU Xin, GAO Fengfei, HU Yumei, CHEN Shenghai
- 林地叶面积指数遥感估算方法适用分析
- Assessment of different methods for estimating forest leaf area index from remote sensing data
- 生态学报, 2014, 34(16): 4612-4619
- Acta Ecologica Sinica, 2014, 34(16): 4612-4619
- http://dx.doi.org/10.5846/stxb201212181819
-
文章历史
- 收稿日期:2012-12-18
- 网络出版日期:2014-3-4
2. 北京师范大学, 遥感科学国家重点实验室, 北京100875;
3. 广西师范大学, 生命科学学院, 桂林541004
2. State key Laboratory of Remote Sensing Science, Beijing Normal University, Beijing 100875, China;
3. College of life Science, Guangxi Normal University, Guilin 541004, China
近年来随着北京生态城市的建设,城区绿化覆盖率已达到了43%,其中林地占大部分比例,城区现有林地改善城市生态环境的效果急需定量数据以建立评估体系。由于植被覆盖率无法准确地描述林地空间结构的差异,而叶面积指数(LAI)是表征林分冠层结构的最基本参量之一,与植被的蒸腾、光合作用、净初级生产力等密切相关[1]。因此,获取可靠的北京城区林地叶面积指数空间分布,是评估林地生态环境效应的研究内容。
林地叶面积指数的获取可利用地面直接或间接观测法、遥感估算法,地面观测法只能获取点或局地值[2],而遥感估算法是快速和准确估算大范围林地LAI的唯一经济可行的方法[3]。遥感估算法可分为物理模型、统计模型、神经网络方法;物理模型的物理意义明确,具有普适性[4],但参数较多,有些参数较难获取,迭代时有些反函数可能不收敛导致结果的不确定性甚至错误,计算量较大。统计模型通过建立遥感植被指数和地面实测LAI的线性或者非线性关系来估算LAI,输入参数较少、计算简单且效率较高,是估算小区域LAI的常用方法[3, 5, 6, 7, 8, 9],也常用于验证MODIS LAI产品的精度,但缺乏普适性。神经网络具有很强的非线性拟合能力[10],成为估算LAI的一种重要手段,但需要足够数量的训练样本。
尽管森林LAI的遥感估算已有大量研究论文发表,但主要针对于流域、山地等大范围自然地表森林[3, 5, 9],较少针对城市林地。Richardson对西雅图市的LIDAR数据使用4种模型估算LAI,发现基于单参数Beer-Lambert定律的模型精度最高[11];Jensen对印第安纳州Terre Haute市的AISA高光谱数据使用回归方程估算LAI,发现当LAI在中、低范围时的估算精度很高,但在5.0—8.0范围时不够准确[12];吴文友研究发现广州市不同类型森林LAI与NDVI的回归模型均可用逻辑斯蒂方程表达[13];林文鹏研究表明上海城市林地实测LAI与SPOT5数据3种植被指数均具有很好的线性回归关系[14]。然而这些研究都是基于单一植被指数和LAI之间的统计关系,由于LAI与植被指数的关系是非线性的,同一植被指数在不同LAI数值范围的统计模型精度不尽相同,导致使用单一植被指数统计模型遥感估算LAI精度提高受到限制。
考虑高光谱遥感数据和LIDAR数据受实验条件的限制较难获取,MODIS LAI产品分辨率较低不适用于城市,而中分辨率Landsat TM遥感图像易于获取且价格较低。因此,以TRAC仪器测定北京城区67块森林样地的LAI,结合同时相TM图像计算的NDVI、SR、RSR、SAVI植被指数,在此基础上分别建立估算林地LAI的多植被指数和单植被指数线性统计模型、单植被指数非线性统计模型、不同多植被指数组合的改进BP神经网络,探讨估算北京城市林地LAI的遥感适用方法和提高精度的途径。
1 研究方法选择与LAI地面观测值具有较好相关性的遥感植被指数,建立LAI估算的统计回归模型和神经网络。
1.1 林地叶面积指数的地面观测研究区为北京市城区,属暖温带半湿润大陆性季风气候,年平均气温约11.5 ℃,年平均相对湿度约56.8%,夏季降水量占全年70%—80%,年平均降水量470—560 mm。
经实地调查后,LAI观测在城区较大面积的林地中进行,根据树种、林木稀疏程度选取30×30 m的样地67块,基本涵盖研究区的主要林分类型和主要树种,样地位于海淀公园、元大都城垣遗址公园、紫竹院公园、玉渊潭公园、龙潭湖公园、陶然亭公园、天坛公园、太阳宫公园、朝阳公园、奥林匹克森林公园、志新路安翔桥旁、机场路旁、北沙滩桥旁等地点,树种主要为国槐(Sophora japonica)、毛白杨(Populus tomentosa)、刺槐(Robinia pseudoacacia)、旱柳(Salix matsudana)、垂柳(Salix babylonica)、金丝垂柳(Salix aureopendula)、银杏(Ginkgo biloba)、洋白蜡(Fontanesia pennsylvanica)、油松(Pinus tabulaeformis)、华山松(Pinus armandii)、白皮松(Pinus bungeana)、侧柏(Platycladus orientalis)、桧柏(Sabina chinensis)、龙柏(Sabina chinensis cv. Kaizuca)、梧桐(Firmiana simplex)等当地优势树种。
LAI的观测选用TRAC仪器进行,TRAC是基于辐射透过率间接测量森林LAI的光学仪器,其从间隙大小分布可以观测到叶片集聚指数,有效地解决了集聚效应问题,LAI的计算不必假设叶片在空间随机分布,减小了从有效LAI转换到实际LAI的计算误差。
2 006年10月17—29日在所选每个样地内取4条长30 m的测线,在每条测线上用TRAC每隔10 m观测采样,将全部测量数据的平均值作为样地的有效LAI值,测量前在样地附近空旷地上测量参考辐射,并用GPS定位样地位置,Trimber差分GPS的定位精度在1 m之内。同时普查树木种类、总株数,测量树木的胸径与树高等数据。
利用TRAC观测的叶片集聚指数ΩE和有效LAIe,扣除树干和树枝等非树叶部分的影响,获得实际LAI:
式中,α为树干等非树叶部分对LAIe的贡献率,γE为簇内所有针叶总表面积与整簇总表面积的比率,对阔叶林γE=1,α和γE通过实地测量估算并参考经验值确定[15]。
1.2 LAI遥感估算方法参考已有的研究,选取NDVI、SR、RSR、SAVI 4种较常用的遥感植被指数作为输入参数:
式中,ρ3、ρ4、ρ5分别为TM图像的红色、近红外、短波红外波段的反射率,ρ5min、ρ5max分别为波段5的最小值、最大值,L为土壤调整参数,北京城区林地一般为中等植被覆盖度,可取L=0.5。
将样地LAI的观测值与对应的遥感植被指数分为训练集和测试集,BP神经网络估算LAI需足够训练样本,根据现有样本总数经实践得到训练集样本需取50个左右,估算结果才可靠。
1.2.1 统计模型通过统计回归分析所选各植被指数训练集与对应LAI观测值之间的关系,从线性、对数、多项式、乘幂、指数等回归方程选取复相关系数的平方R2最接近1的方程作为最优统计模型,然后利用测试集检验模型精度。
1.2.2 改进BP神经网络方法增加动量项提高网络的训练速度,调整学习率缩短学习时间。将所选各植被指数VI进行归一化处理到[0,1]范围内,消除数量级大小不同的影响:
然后将VInor作为输入,而LAI观测值作为输出,根据实践经验网络设置3层,输入层与隐层之间取tansig传递函数,隐层与输出层之间取purelin线性传递函数,初设学习率、动量系数、训练误差、最大训练次数,使用训练样本训练网络达到设定精度,获得估算LAI的神经网络,网络经测试集检验后应用于整幅图像。
1.3 遥感估算方法的评价指标为比较不同方法估算LAI的适用性,对预留的观测值LAImea与VI测试集,由VI计算获得遥感估算值LAImod。对LAImod和相应的LAImea,计算线性回归方程的R2和均方根差RMSE作为评价指标:
式中,n为测试集样本数,R2越接近1,RMSE越接近0,则LAI的估算值越接近于观测值。
2 实验结果与分析由于2006年秋季覆盖北京市城八区的TM卫星图像只有11月1日是无云层存在、可用的,该日最高气温为19℃,LAI观测的起始日期至该日期间无强降温天气、无明显落叶的迹象。因此选取该日图像,经辐射校正、配准和几何精纠正,误差小于0.5个像元;利用MODTRAN进行大气校正,计算植被指数NDVI、SR、RSR、SAVI。
将林地样方LAI观测值与对应的植被指数随机分为训练集47个样本和测试集20个样本,利用训练集建立LAI遥感估算方法,而使用测试集检验精度。
2.1 线性统计模型的比较对训练集植被指数与LAI观测值,分别建立4植被指数和单植被指数与LAI的线性回归方程,计算R2,然后利用测试集计算RMSE,结果见表 1。经F检验,模型精度均达到0.01极显著水平;4植被指数线性回归模型与单植被指数模型相比,R2明显提高,而RMSE明显降低。原因是利用每一植被指数估算LAI时均存在一定的缺陷,多植被指数组合可实现信息互补,故一般多植被指数模型估算LAI的精度比单植被指数模型高,但哪几个植被指数的组合最优将在后面的神经网络方法中进行比较。
植被指数
Vegetation index | 回归方程
Regression equation |
决定系数
R2 | 均方根差
RMSE |
NDVI归一化植被指数Normalized difference vegetation index; SR比值植被指数Simple ratio index; RSR减化比值植被指数Reduced simple ratio index; SAVI土壤调整植被指数Soil adjusted vegetation index; LAI叶面积指数Leaf area index | |||
NDVI | LAI=12.632 NDVI-4.033 | 0.553 | 0.528 |
SR | LAI=1.427 SR -2.071 | 0.566 | 0.512 |
RSR | LAI=1.672 RSR -1.948 | 0.593 | 0.497 |
SAVI | LAI=11.475 SAVI -0.572 | 0.551 | 0.539 |
NDVI、SR、RSR、SAVI | LAI=-0.079NDVI +0.229SR +0.934RSR+ 4.678SAVI -1.945 | 0.715 | 0.394 |
利用训练集分别建立各单植被指数与LAI的对数、二项式、乘幂、指数等非线性回归方程,结果见表 2。由表 1、表 2可见二项式回归方程的R2均高于线性回归方程和其它非线性回归方程,RSR的二项式回归方程为单植被指数估算LAI的最优统计模型,其R2为0.632,最接近于1。单植被指数与LAI呈非线性关系。
植被指数
Vegetation index | 方程类型
Equation type | 回归方程
Regression equation | 决定系数
R2 |
NDVI | 对数 | LAI=6.7669+ 6.4530 ln NDVI | 0.534 |
二项式 | LAI=59.408 NDVI2-49.469 NDVI+12.070 | 0.583 | |
乘幂 | LAI=12.427 NDVI2.4798 | 0.530 | |
指数 | LAI=0.1970 e4.8351NDVI | 0.543 | |
SR | 对数 | LAI=-2.750+4.5366 ln SR | 0.550 |
二项式 | LAI=0.7278 SR 2-3.3216 SR +5.5667 | 0.588 | |
乘幂 | LAI=0.3191 SR 1.7475 | 0.540 | |
指数 | LAI=0.4175e0.5475SR | 0.551 | |
RSR | 对数 | LAI=-1.7367+4.3607 ln RSR | 0.575 |
二项式 | LAI=1.1089 RSR 2-4.2897 RSR +5.9503 | 0.632 | |
乘幂 | LAI=0.4672 RSR1.6894 | 0.562 | |
指数 | LAI=0.4346 e0.6441RSR | 0.575 | |
SAVI | 对数 | LAI=6.4822+2.9893 ln SAVI | 0.510 |
二项式 | LAI=40.448 SAVI 2-11.273 SAVI +2.5281 | 0.580 | |
乘幂 | LAI=10.655 SAVI 1.1158 | 0.469 | |
指数 | LAI=0.7716 e4.254SAVI | 0.501 |
从前面数据分析可知,多植被指数模型估算LAI的精度比单植被指数模型高,单植被指数与LAI二项式回归方程的R2最高,然而多植被指数与LAI的非线性统计模型随植被指数组合个数的增加变为很复杂,4植被指数的二项式回归方程为:
需要确定的系数有15个,构建模型较难,考虑神经网络具有较强的非线性问题求解能力,表 2数据表明测试的4种遥感植被指数与实测LAI数据之间存在统计相关性,因此后面的实验将采用这些植被指数作为神经网络的输入估算LAI。
建立估算LAI的3层BP神经网络,利用训练集训练网络,输入层节点数与所选多植被指数的个数相同,输出为LAI观测值,根据实践经验设置学习率0.03、动量系数0.8、最大训练误差1×10-5、最大训练次数1000。采用试探法确定隐层节点数,调整学习率和动量因子直至网络训练达到设定精度要求。
然后,利用学习训练确定的神经网络输入测试集多植被指数组合,获得LAI估算值,图 1为LAI估算值与观测值构成的散点图,实线为线性回归方程线,虚线为估算值与观测值一致的1 ∶ 1线,回归方程线越接近1 ∶ 1虚线,则估算值与观测值的误差越小。计算R2和RMSE,结果见表 3。
植被指数组合
Vegetation index combination | 决定系数
R2 | 均方根差
RMSE |
NDVI,SR | 0.753 | 0.218 |
NDVI,RSR | 0.765 | 0.220 |
NDVI,SAVI | 0.648 | 0.271 |
SR,RSR | 0.809 | 0.192 |
SR,SAVI | 0.739 | 0.225 |
RSR,SAVI | 0.754 | 0.221 |
NDVI,SR,RSR | 0.796 | 0.194 |
NDVI,SR,SAVI | 0.666 | 0.249 |
NDVI,RSR,SAVI | 0.827 | 0.189 |
SR,RSR,SAVI | 0.795 | 0.198 |
NDVI,SR,RSR,SAVI | 0.776 | 0.206 |
由图 1可见NDVI、RSR、SAVI组合神经网络的LAI估算值与观测值的回归方程线最接近1 ∶ 1虚线,由表 3可知其R2为0.827,最高,RMSE为0.189,误差最低。因此,多植被指数组合时个数的增加,将使得LAI的估算方法的复杂性和计算量增加,但并不一定能提高LAI的估算精度;NDVI、RSR、SAVI是估算LAI的最优组合,精度最高,NDVI、RSR、SAVI组合的神经网络构成北京城市林地LAI遥感估算的适用方法。
2.4 城市林地叶面积指数的遥感制图对北京城区TM图像,结合混合像元分解与支持向量机方法获取的地物分类图[16],对林地LAI的估算利用上述NDVI、RSR、SAVI组合神经网络,对农业用地LAI的估算由于缺少观测数据使用已有研究文献的统计模型[17],得到估算的LAI空间分布图,见图 2。由图可知,LAI的空间分布呈现在主城区外围数值较大,而在主城区内数值较小,但主城区中的天坛公园、朝阳公园、玉渊潭公园、元大都城垣遗址公园、紫竹院公园、龙潭湖公园、陶然亭公园、太阳宫公园、奥林匹克森林公园、海淀公园等公园林地的LAI数值相对较大。
利用样地LAI的全部观测值检验图 2中估算值,估算值与观测值的均方根差RMSE为0.113,LAI的估算值与观测值之间误差较小,精度较高。
3 结果与讨论北京城区林地中较大面积相同树种的纯林较少,纯阔叶林和针叶林的样本太少导致在建立估算LAI的遥感模型时只能对阔叶林和针叶林混合建模,而阔叶林和针叶林的叶面积指数与遥感植被指数之间的非线性关系一般存在差异,因此,遥感单植被指数统计模型估算LAI的精度不高,单植被指数非线性统计模型估算LAI的精度高于线性统计模型。
神经网络具有较强的描述叶面积指数与多植被指数组合的非线性关系的能力,解决了统计模型中多植被指数组合与叶面积指数的非线性回归方程系数的数目较多导致其确定较复杂的难题。
NDVI、RSR、SAVI组合神经网络构成了北京城市林地叶面积指数遥感估算的适用方法。多植被指数组合作为输入估算LAI的神经网络中,NDVI、RSR、SAVI组合是最优组合,精度最高,估算值与观测值线性回归方程的R2最高,为0.827,RMSE最低,为0.189。
估算LAI的多植被指数组合模型比单植被指数模型的精度高,因为LAI与植被指数的关系是非线性的,每种植被指数在估算LAI时都存在一定的缺陷,不同植被指数之间可以相互补充信息。并非模型输入的植被指数个数越多越好,随植被指数个数的增加,模型的复杂性和计算量增加,但并不一定能提高估算精度,神经网络提供了在实测数据的基础上寻找估算LAI的最优多植被指数组合的方法,计算效率较高。
致谢: 中国科学院地理科学与资源研究所张仁华研究员、孙晓敏研究员、北京园林研究所李延明高工、中国科学院大气物理研究所胡非研究员、北京师范大学遥感科学国家重点实验室杨胜天教授、孙睿教授为“北京城市绿地对水、热、CO2通量调节功能的遥感定量研究”项目研究和实验提供了大量建议与支持,特此致谢。
[1] | Yan H, Wang S Q, Billesbach D, Oechel W, Zhang J H, Meyers T, Martin T A, Matamala R, Baldocchi D, Bohrer G, Dragoni D, Scott R. Global estimation of evapotranspiration using a leaf area index-based surface energy and water balance model. Remote Sensing of Environment, 2012, 124(9): 581-595. |
[2] | Majasalmi T, Rautiainen M, Stenberg P, Rita H. Optimizing the sampling scheme for LAI-2000 measurements in a boreal forest. Agricultural and Forest Meteorology, 2012, 154-155(3): 38-43. |
[3] | Zhu G L, Ju W M, Chen J M, Fan W Y, Zhou Y L, Li X F, Li M Z. Forest canopy leaf area index in Maoershan Mountain: Ground measurement and remote sensing retrieval. Chinese Journal of Applied Ecology, 2010, 21(8): 2117-2124. |
[4] | González-Sanpedro M C, Toan T L, Moreno J, Kergoat L, Rubio E. Seasonal variations of leaf area index of agricultural fields retrieved from Landsat data. Remote Sensing of Environment, 2008, 112(3): 810-824. |
[5] | Madugundu R, Nizalapur V, Jha C S. Estimation of LAI and above-ground biomass in deciduous forests: Western Ghats of Karnataka, India. International Journal of Applied Earth Observation and Geoinformation, 2008, 10(2): 211-219. |
[6] | Liu J G, Pattey E, Jégo G. Assessment of vegetation indices for regional crop green LAI estimation from Landsat images over multiple growing seasons. Remote Sensing of Environment, 2012, 123(8): 347-358. |
[7] | Zarate-Valdez J L, Whiting M L, Lampinen B D, Metcalf S, Ustin S L, Brown P H. Prediction of leaf area index in almonds by vegetation indexes. Computers and Electronics in Agriculture, 2012, 85(7): 24-32. |
[8] | Liu Y B, Ju W M, Zhu G L, Chen J M, Xing B L, Zhu J F, Zhou Y L. Retrieval of leaf area index for different grasslands in Inner Mongolia prairie using remote sensing data. Acta Ecologica Sinica, 2011, 31(18): 5159-5170. |
[9] | He L, Tang S Y, Miao F, Li Y X. Remote sensing model and inversion of leaf area index in the typical area of Minjiang Valley upriver. Research of Soil and Water Conservation, 2010, 17(1): 218-221. |
[10] | Heiskanen J, Rautiainen M, Korhonen L, Mõttus M, Stenberg P. Retrieval of boreal forest LAI using a forest reflectance model and empirical regressions. International Journal of Applied Earth Observation and Geoinformation, 2011, 13(4): 595-606. |
[11] | Richardson J J, Moskal L M, Kim S H. Modeling approaches to estimate effective leaf area index from aerial discrete-return LIDAR. Agricultural and Forest Meteorology, 2009, 149(6/7): 1152-1160. |
[12] | Jensen R R, Hardin P J, Bekker M, Farnes D S, Lulla V, Hardin A. Modeling urban leaf area index with AISA+ hyperspectral data. Applied Geography, 2009, 29(3): 320-332. |
[13] | Wu W Y, Wu Z M, Hu H R, Jia B Q. Estimation on leaf area index of urban forests in Guangzhou by remote sensing technique. Journal of Northeast Forestry University, 2010, 38(3): 38-41. |
[14] | Lin W P, Zhao M, Zhang Y F, Liu Y L, Liu D Y, Gao J. Study on estimation of urban forest LAI models based on SPOT5. Science of Surveying and Mapping, 2008, 33(2): 57-59, 63-63. |
[15] | Chen J M, Govind A, Sonnentag O, Zhang Y Q, Barr A, Amiro B. Leaf area index measurements at Fluxnet-Canada forest sites. Agricultural and Forest Meteorology, 2006, 140(1/4): 257-268. |
[16] | Wang X X, Wu H, Lu X C, Wu X J, Lou L E, Zhu Q J. Extracting urban green space with mixed pixel decomposing and SVM. Computer Engineering and Applications, 2009, 45(33): 216-217, 226-226. |
[17] | Liu D S, Li S M. Statistical relationship between LAI indices and canopy spectral data of winter wheat in BEIJING area. Remote Sensing for Land & Resources, 2008, 20(4): 32-34, 42-42. |
[3] | 朱高龙, 居为民, Chen J M, 范文义, 周艳莲, 李显风, 李明泽. 帽儿山地区森林冠层叶面积指数的地面观测与遥感反演. 应用生态学报, 2010, 21(8): 2117-2124. |
[8] | 柳艺博, 居为民, 朱高龙, 陈镜明, 邢白灵, 朱敬芳, 周艳莲. 内蒙古不同类型草地叶面积指数遥感估算. 生态学报, 2011, 31(18): 5159-5170. |
[9] | 何磊, 唐姝娅, 苗放, 李玉霞. 岷江上游典型流域叶面积指数的遥感模型及反演. 水土保持研究, 2010, 17(1): 218-221. |
[13] | 吴文友, 吴泽民, 胡鸿瑞, 贾保全. 基于遥感技术的广州市城市森林叶面积指数推算. 东北林业大学学报, 2010, 38(3): 38-41. |
[14] | 林文鹏, 赵敏, 张翼飞, 柳云龙, 刘冬燕, 高峻. 基于SPOT5遥感影像的城市森林叶面积指数反演. 测绘科学, 2008, 33(2): 57-59, 63-63. |
[16] | 王修信, 吴昊, 卢小春, 吴学军, 罗兰娥, 朱启疆. 利用混合像元分解结合SVM提取城市绿地. 计算机工程与应用, 2009, 45(33): 216-217, 226-226. |
[17] | 刘东升, 李淑敏. 北京地区冬小麦冠层光谱数据与叶面积指数统计关系研究. 国土资源遥感, 2008, 20(4): 32-34, 42-42. |