免费文献传递   相关文献

Comparison of three regression analysis methods for application to LAI inversion using Hyperion data

三种回归分析方法在Hyperion影像LAI反演中的比较



全 文 :
摇 摇 摇 摇 摇 生 态 学 报
摇 摇 摇 摇 摇 摇 摇 (SHENGTAI XUEBAO)
摇 摇 第 32 卷 第 24 期摇 摇 2012 年 12 月摇 (半月刊)
目摇 摇 次
从文献计量角度分析中国生物多样性研究现状 刘爱原,郭玉清,李世颖,等 (7635)……………………………
CO2 浓度升高和模拟氮沉降对青川箭竹叶营养质量的影响 周先容,汪建华,张摇 红,等 (7644)………………
陕西中部黄土高原地区空气花粉组成及其与气候因子的关系———以洛川县下黑木沟村为例
吕素青,李月从,许清海,等 (7654)
…………………
……………………………………………………………………………
长三角地区 1995—2007 年生态资产时空变化 徐昔保,陈摇 爽,杨桂山 (7667)…………………………………
基于智能体模型的青岛市林地生态格局评价与优化 傅摇 强,毛摇 锋,王天青,等 (7676)………………………
青藏高原高寒草地生态系统服务功能的互作机制 刘兴元,龙瑞军,尚占环 (7688)……………………………
北京城市绿地的蒸腾降温功能及其经济价值评估 张摇 彪,高吉喜,谢高地,等 (7698)…………………………
武汉市造纸行业资源代谢分析 施晓清,李笑诺,赵吝加,等 (7706)………………………………………………
丽江市家庭能耗碳排放特征及影响因素 王丹寅,唐明方,任摇 引,等 (7716)……………………………………
基于分布式水文模型和福利成本法的生态补偿空间选择研究 宋晓谕,刘玉卿,邓晓红,等 (7722)……………
设施塑料大棚风洞试验及风压分布规律 杨再强,张摇 波,薛晓萍,等 (7730)……………………………………
湖南珍稀濒危植物———珙桐种群数量动态 刘海洋,金晓玲,沈守云,等 (7738)…………………………………
云南岩陀及其近缘种质资源群体表型多样性 李萍萍,孟衡玲,陈军文,等 (7747)………………………………
沙埋和种子大小对柠条锦鸡儿种子萌发、出苗和幼苗生长的影响 杨慧玲,梁振雷,朱选伟,等 (7757)………
栗山天牛天敌花绒寄甲在栎林中的种群保持机制 杨忠岐,唐艳龙,姜摇 静,等 (7764)…………………………
基于相邻木排列关系的混交度研究 娄明华,汤孟平,仇建习,等 (7774)…………………………………………
三种回归分析方法在 Hyperion影像 LAI反演中的比较 孙摇 华,鞠洪波,张怀清,等 (7781)…………………
红松和蒙古栎种子萌发及幼苗生长对升温与降水综合作用的响应 赵摇 娟,宋摇 媛,孙摇 涛,等 (7791)………
新疆杨边材贮存水分对单株液流通量的影响 党宏忠,李摇 卫,张友焱,等 (7801)………………………………
火干扰对小兴安岭毛赤杨沼泽温室气体排放动态影响及其影响因素 顾摇 韩,牟长城,张博文 (7808)………
不同潮汐和盐度下红树植物幼苗秋茄的化学计量特征 刘滨尔,廖宝文,方展强 (7818)………………………
腾格里沙漠东南缘沙质草地灌丛化对地表径流及氮流失的影响 李小军,高永平 (7828)………………………
西双版纳人工雨林群落结构及其林下降雨侵蚀力特征 邓摇 云,唐炎林 ,曹摇 敏,等 (7836)…………………
西南高山地区净生态系统生产力时空动态 庞摇 瑞,顾峰雪,张远东,等 (7844)…………………………………
南北样带温带区栎属树种种子化学组成与气候因子的关系 李东胜,史作民,刘世荣,等 (7857)………………
模拟酸雨对龙眼叶片 PS域反应中心和自由基代谢的影响 李永裕,潘腾飞,余摇 东,等 (7866)………………
沈阳市城郊表层土壤有机污染评价 崔摇 健,都基众,马宏伟,等 (7874)…………………………………………
降雨对旱作春玉米农田土壤呼吸动态的影响 高摇 翔,郝卫平,顾峰雪,等 (7883)………………………………
冬季作物种植对双季稻根系酶活性及形态指标的影响 于天一,逄焕成,任天志,等 (7894)……………………
施氮量对小麦 /玉米带田土壤水分及硝态氮的影响 杨蕊菊,柴守玺,马忠明 (7905)……………………………
微山湖鸟类多样性特征及其影响因子 杨月伟,李久恩 (7913)……………………………………………………
新疆北部棉区作物景观多样性对棉铃虫种群的影响 吕昭智,潘卫林,张摇 鑫,等 (7925)………………………
杭州西湖北里湖沉积物氮磷内源静态释放的季节变化及通量估算 刘静静,董春颖,宋英琦,等 (7932)………
基于实码遗传算法的湖泊水质模型参数优化 郭摇 静,陈求稳,张晓晴,等 (7940)………………………………
气候环境因子和捕捞压力对南海北部带鱼渔获量变动的影响 王跃中,孙典荣,陈作志,等 (7948)……………
象山港南沙岛不同养殖类型沉积物酸可挥发性硫化物的时空分布 颜婷茹,焦海峰,毛玉泽,等 (7958)………
专论与综述
提高植物抗寒性的机理研究进展 徐呈祥 (7966)…………………………………………………………………
植被对多年冻土的影响研究进展 常晓丽,金会军,王永平,等 (7981)……………………………………………
凋落物分解主场效应及其土壤生物驱动 査同刚,张志强,孙摇 阁,等 (7991)……………………………………
街尘与城市降雨径流污染的关系综述 赵洪涛,李叙勇,尹澄清 (8001)…………………………………………
期刊基本参数:CN 11鄄2031 / Q*1981*m*16*374*zh*P* ¥ 70郾 00*1510*40*
室室室室室室室室室室室室室室
2012鄄12
封面图说: 永兴岛海滩植被———永兴岛是中国西沙群岛的主岛,也是西沙群岛及南海诸岛中最大的岛屿。 国务院 2012 年 6 月
批准设立的地级三沙市,管辖西沙群岛、中沙群岛、南沙群岛的岛礁及其海域,三沙市人民政府就驻西沙永兴岛。 永
兴岛岛上自然植被密布,野生植物有 148 种,占西沙野生植物总数的 89% ,主要树种有草海桐(羊角树)、麻枫桐、野
枇杷、海棠树和椰树等。 其中草海桐也称为羊角树,是多年生常绿亚灌木植物,它们总是喜欢倚在珊瑚礁岸或是与
其他滨海植物聚生于海岸沙滩,为典型的滨海植物。
彩图提供: 陈建伟教授摇 北京林业大学摇 E鄄mail: cites. chenjw@ 163. com
第 32 卷第 24 期
2012 年 12 月
生 态 学 报
ACTA ECOLOGICA SINICA
Vol. 32,No. 24
Dec. ,2012
http: / / www. ecologica. cn
基金项目:国家“十二五冶863 项目(2012AA102001);林业公益性行业科研专项(201104028)
收稿日期:2011鄄11鄄29; 摇 摇 修订日期:2012鄄10鄄10
*通讯作者 Corresponding author. E鄄mail: ju@ caf. ac. cn
DOI: 10. 5846 / stxb201111291819
孙华, 鞠洪波, 张怀清,林辉, 凌成星.三种回归分析方法在 Hyperion影像 LAI反演中的比较.生态学报,2012,32(24):7781鄄7790.
Sun H, Ju H B, Zhang H Q, Lin H, Ling C X. Comparison of three regression analysis methods for application to LAI inversion using Hyperion data. Acta
Ecologica Sinica,2012,32(24):7781鄄7790.
三种回归分析方法在 Hyperion影像
LAI反演中的比较
孙摇 华1,2, 鞠洪波1,*, 张怀清1,林摇 辉2, 凌成星1
(1. 中国林业科学研究院资源信息所, 北京摇 100091;2. 中南林业科技大学林业遥感信息工程研究中心, 长沙摇 410004 )
摘要:借助 GPS进行地面精确定位,利用 LAI鄄 2000 冠层分析仪在攸县黄丰桥林场开展 130 个样地(60m伊60m)的叶面积指数
(Leaf Area Index,LAI)测量。 采用 FLAASH模块对 Hyperion数据进行大气校正并与地面同步冠层观测数据进行拟合,通过研究
地面实测 LAI与 Hyperion影像波段及其衍生的系列植被指数(NDVI、RVI等)的相关性,筛选出估算叶面积指数的植被指数因
子。 应用曲线估计、逐步回归及偏最小二乘三种回归分析技术分别建立叶面积指数的最优估算模型。 结果表明:参与建模的因
子中,比值植被指数(RVI)与 LAI的相关性最大,敏感性最高,其次是 SARVI0. 1, NDVI705, NDVI,SARVI0. 1, SARVI0. 25;曲线估
计、逐步回归分析和偏最小二乘回归三种分析方法所建的 6 个回归模型中,偏最小二乘回归的拟合效果最好,预测值与实测值
的决定系数 R2 为 0. 84、曲线估计的拟合效果最低,预测值与实测值的决定系数 R2 为 0. 64;建模精度分析表明,选用 5—6个自
变量因子进行 LAI建模是可靠的,以 6 个植被因子建立的偏最小二乘回归模型预测精度最高。
关键词:遥感反演;叶面积指数;偏最小二乘回归;植被指数;黄丰桥林场
Comparison of three regression analysis methods for application to LAI inversion
using Hyperion data
SUN Hua1,2, JU Hongbo1,*, ZHANG Huaiqing1, LIN Hui2, LING Chengxing1
1 Research Institute of Forest Resources Information Technique Chinese Academy of Forestry, Beijing 100091, China
2 Research Center of Forestry Remote Sensing & Information Engineering Central South University & Technology, Changsha 410004, China
Abstract: This paper focuses on Leaf Area index (LAI) inversion, using EO鄄 1 Hyperion data for Huangfengqiao forest
farm, YouXian County, Hunan Province. First, LAI was acquired using a LAI鄄 2000 canopy analyzer at 130 sample plots
(60 m 伊 60 m), with a Global Positioning System (Trimble GPS Geo XT). Second, atmospheric correction was applied to
Hyperion data using the ground鄄synchronous canopy observation data. Third, effective vegetation indexes were selected to
estimate LAI, according to research on correlation between LAI, bands and vegetation indexes derived from Hyperion
imagery. Finally, an optimal estimation model of LAI was built by curve estimation, stepwise regression, and a partial
least鄄squares regression algorithm. Results show that sensitivity of ratio vegetation index (RVI) was highest among all
model factors, followed by SARVI0. 1, NDVI705, NDVI, SARVI0. 1, and SARVI0. 25 . Among all fit models, the effect of the
partial least鄄squares regression was best, with R2 coefficient 0. 84, whereas the curve estimation effect was worst, with R2
coefficient 0. 64. Model precision analysis shows that it is reliable to build the model using 5 to 6 independent variables,
and prediction accuracy of the partial least鄄square regression was the greatest.
Key Words: remote sensing inversion; leaf area index;partial least鄄squares regression; vegetation index;Huangfengqiao
http: / / www. ecologica. cn
forest farm
叶面积指数(Leaf Area Index, LAI)是表征植被冠层特征的主要指标之一[1鄄3],同时也是决定生态系统净
初级生产力的重要因子[1],是植被定量遥感研究的重要参数,对全球变化和生态系统碳循环研究具有重要意
义。 目前小范围内的叶面积指数一般通过仪器进行测量,大范围的叶面积指数通过遥感反演模型模拟获得。
现有的研究主要集中于主要农作物的叶面积指数与生理生化参数反演[4鄄7],森林郁闭度[8]、森林生物量[9鄄10]
的估测等方面。 用于反演的数据一般以 Landsat、MODIS 等中低分辨率的多光谱影像为主[11鄄20],少量机载高
光谱数据(CASI、PHI、HyMap)及部分成像光谱数据如 Hyperion,CHRIS等。
叶面积指数的遥感反演一般通过两种途径实现:物理模型法和经验模型法。 物理模型法主要以
PROSPECT+SAIL模型为主,从物理机理角度反演植被叶面积指数[21鄄24]。 反演过程需要输入叶片的生化参数
和结构参数。 经验模型法主要通过利用遥感技术获取各种植被指数,计算植被指数与 LAI 的回归关系模型,
达到间接反演林分叶面积指数的效果[25鄄26]。 经验模型法虽然存在模型参数随时间或研究区域变化的缺陷,
但具有输入参数少、计算效率高、容易实现等优点,且大量的研究表明,植被指数与 LAI 指数之间具有良好的
定量关系,这是经验模型反演 LAI的理论基础。 归纳现有经验模型的研究成果可知,LAI 反演回归模型以单
因子的线性与非线性回归居多,大部分论文只采用一种回归方式,使用多种方法进行比较的少,论文以
Hyperion影像为研究对象,采用曲线估计,逐步回归,偏最小二乘回归方法在攸县黄丰桥林场开展叶面积指数
(LAI)反演研究,并对 3 种方法所得结果进行对比,旨在分析不同回归分析方法的特点,找出对 LAI 反演敏感
性高的植被指数因子。
图 1摇 研究区位置图
Fig. 1摇 Location of study area
1摇 研究区概况
黄丰桥国有林场呈带状分布横跨株洲市攸县东西
部,介于东经 113毅04忆至 113毅43忆,北纬 27毅06忆至 27毅24忆
之间。 东北部与江西的莲花、萍乡交界,东南与茶陵县
接壤,西北部与株洲、醴陵毗邻。 现有林地总面积
10122郾 6 hm2,森林覆盖率为 86. 24% 。 全场地貌以中
低山为主,境内最高海拔 1270 m,最低海拔 115 m,坡度
介于 20—35毅度之间。 林场地处中亚热带季风湿润气
候区,年均气温 17. 8益;平均无霜期为 292 d;年均降水
量 1410. 8mm。 此次研究的区域主要包括柏市、广黄和
长冲 3 个分场(图 1)。
2摇 材料与方法
2. 1摇 数据处理
2. 1. 1摇 Hyperion影像处理
研究所用星载高光谱影像为 Hyperion 影像的 L1G 数据,通过美国地质调查局网站下载(http: / / glovis.
usgs. gov / ),影像的成像时间为 2010 年 12 月 28 日,原始数据只进行了简单的辐射校正,在 ENVI4. 5 软件中
使用 FLAASH大气校正模块对影像进行大气校正,校正时选择合适的大气模式、水汽含量、气溶胶模型、波谱
分辨率和多散射模型等参数,经过大气校正后图像质量得到明显提高。 利用 SPOT5 的 10m 分辨率多光谱影
像为参考影像对校正后的 Hyperion数据进行几何校正,共选取 24 个明显地物点,系统均方根误差为 0. 41 个
像元,基本满足研究所需的精度。 为了更好的分析大气校正的结果,在研究区采用美国 ASD 公司的 Field
Spec手持式地物光谱仪,同步测量了 11 个杉木样地的光谱反射率地物,并与大气校正后的数据进行了数据
拟合,经过大气校正后的地物光谱曲线和地面观测到的光谱曲线具有较好的一致性( R2 为 0. 94)。
2. 1. 2摇 地面样地 LAI数据获取
样地内的叶面积指数测量使用的是美国 LI鄄COR公司生产的 LAI鄄2000 仪器,调查期内研究区的天气以阴
2877 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
天为主,比较适合测量。 在 60m伊60m的样地范围内,测量 4 个角上 10m伊10m样方及样地中心位置 10m伊10m
样方的 LAI值,取 5 次测量的均值作为样地的 60m伊60m范围的 LAI值,共测量 130 个样地,同时,在样地范围
内开展了林分郁闭度,土壤类型,土壤厚度,年龄,胸径和平均树高等因子的调查,并进行记录。
2. 2摇 植被指数提取
研究在参考前人工作的基础上,共选取归一化植被指数(NDVI),比值植被指数(RVI),土壤修正植被指
数(SAVI),抗土壤和大气的植被指数(SARVI),修正的归一化差值植被指数(MNDVI),红边归一化植被指数
NDVI(705),三角植被指数(TVI)等 7 个植被指数。 其中 SAVI和 SARVI带有修正因子(L取 0. 1,0. 25,0. 35,
0. 5),实际共计算了 13 个因子。 每个植被指数计算方式及在 Hyperion数据中对应的波段选择如表 1 所示。
表 1摇 用于 LAI反演的植被指数
Table 1摇 Vegetation index for LAI inversion
植被指数
Vegetation index
计算公式
Formula
Hyperion对应波段
Band of Hyperion image
参考文献
References
归一化植被指数(NDVI) NDVI =
籽n - 籽r
籽n + 籽r
band48,band32 Tomoaki Miura,2006[28]
比值植被指数(RVI) RVI =
籽n
籽r
band48,band32 Martin Schlerfa,2005[29]
土壤修正植被指数(SAVI) SAVI =
(籽n - 籽r)(1 + L)
籽n + 籽r + L
band48,band32 Peng Gong,2003[30]
抗土壤和大气的植被指数
(SARVI)
SARVI =
(籽n - 籽rb)(1 + L)
籽n + 籽rb + L
籽rb = 籽r - 茁(籽b - 籽r) , 茁取 0. 5 或 1
band48,band32 和 band13 AR Huete,1997[31]
修正的归一化差值植被指数
(MNDVI)
NDVI =
籽n - 籽r
籽n + 籽r
(1 -
籽s - 籽smin
籽smax - 籽smin
) band48,band32 和 band150 AR Huete,1997[31]
红边归一化植被指数
(NDVI705)
NDVI705 =
籽750 - 籽705
籽750 + 籽705
band40,band35 MR Mobasheri,2007[32]
三角植被指数(TVI) TVI =
120(籽n - 籽g) - 200(籽r - 籽g)
2
band48,band32 和 band22 Shunlin Liang, 2004[33]
2. 3摇 叶面积指数反演
对野外调查所得的 130 个样本数据通过计算残差图分布进行数据筛选,剔除了 8 个离群值大的数据。 随
机选取 80 个样本数据进行叶面积指数反演建模,利用剩余 42 组数据作为测试数据对建立的反演模型进行精
度检验。 叶面积指数反演分 3 个阶段进行,首先开展单因子曲线估计建模,即每个植被指数与叶面积指数建
立回归模型,选择 SPSS软件曲线估计所列举的 11 种模型进行分析;其次,通过分析植被指数与叶面积指数的
相关性,选取逐步回归分析方法以及偏最小二乘回归分析技术进行自变量的筛选,对筛选出的因子建立叶面
积指数反演模型;最后对 3 种建模方式所得结果进行精度分析和比较。
2. 4摇 模型精度评价
模型精度是衡量预测方法是否适用于预测对象的一个重要指标。 论文选用残差标准差和模型评价指标
值(Model Evaluation Indicator,MEI) [34]来衡量 3 种回归分析方法的建模精度。 计算公式为:
SSe =移(yi - y^i) 2
滓^2 = SSe
n - 2
滓^2 R =
RSSe
m - q
滓^2 F =
RSSF
n - m
MEI = 0. 5滓^2 R + 0. 5滓^2
ì
î
í
ï
ï
ï
ï
ï
ï
ï
ï
ïï
F
3877摇 24 期 摇 摇 摇 孙华摇 等:三种回归分析方法在 Hyperion影像 LAI反演中的比较 摇
http: / / www. ecologica. cn
式中,SSe表示残差平方和; 滓^表示残差标准差;n表示样地总数; m 表示参与反演的样地数量;q 参与反演的
自变量个数; RSSe 表示 m个参与反演样地的残差平方和;RSSF 表示用于精度检验的 n-m个测试数据样地残
差平方和。 原则上应选择 滓^2R 和 滓^2 F 都达到最小的最优变量选择方法,考虑建立方程与精度评价同等重要,一
般取模型评价指标值最小的方法为最优模型方法。
3摇 结果与分析
3. 1摇 曲线估计模型
研究所选用的 13 种植被指数均参与建模,为了对拟合结果进行比较,各植被指数与 LAI 单独回归拟合,
在拟合方法上采用 SPSS中的 11 种曲线拟合方法,最终拟合结果保留 2 种拟合曲线,即各植被指数与 LAI 的
线性回归模型以及一个最佳的曲线拟合模型。
表 2摇 各植被指数与 LAI最佳曲线估计拟合结果
Table 2摇 The best curve fitting results for LAI and Vegetation index
植被指数
Vegetation index
线性回归方程
Linear regression model R
2 最佳回归方程
Best regression model R
2
MNDVI Y=5. 065x -0. 164 0. 49 Y=0. 827 x1. 984 0. 55
NDVI Y=5. 332x-0. 348 0. 57 Y=exp(2. 023-0. 8 / x ) 0. 60
RVI Y=0. 418x + 0. 177 0. 73 Y=1. 033 exp(0. 155 x ) 0. 74
NDVI705 Y=4. 223x + 0. 427 0. 53 Y=1. 725 + 2. 877 x2 + 0. 517 x3 0. 54
SARVI0. 1 Y=5. 187x + -0. 508 0. 68 Y=0. 759 exp(1. 96 x ) 0. 71
SARVI0. 25 Y=4. 521x + 0. 03 0. 55 Y=exp(1724 +0. 9 / x ) 0. 59
SARVI0. 35 Y=4. 199x + 0. 092 0. 48 Y=exp(1. 900 -0. 569 / x ) 0. 53
SARVI0. 5 Y=3. 428x + 0. 61 0. 35 Y=exp(1. 35 +1. 115 / x ) 0. 40
SAVI0. 1 Y=3. 493x + 0. 721 0. 37 Y=exp(1. 183 +1. 352 / x ) 0. 41
SAVI0. 25 Y=4. 036x + 0. 281 0. 42 Y=exp(1. 559+0. 973 / x ) 0. 49
SAVI0. 35 Y=4. 428x + 0. 071 0. 53 Y=exp(1. 7+0. 924 / x ) 0. 58
SAVI0. 5 Y=4. 032x + 0. 314 0. 45 Y=exp(1. 055 +1. 569 / x ) 0. 51
TVI Y=0. 505x + 2. 068 0. 24 Y=1. 818 +0. 351x+ 0. 492x2 -0. 148x3 0. 30
从表 2 分析可知,从 LAI与各植被指数的回归模型来看,RVI植被指数与 LAI回归模型的决定系数( R2)
最高(线性 0. 73,指数方程 0. 74),NDVI回归模型次之。 TVI决定系数最低(线性 0. 24,三次方曲线模型 0. 3)
SARVI0. 5 决定系数次之(线性 0. 35,S形曲线模型 0. 40)。 从决定系数中可以看出,当决定系数小于 0. 55 时,
LAI与单个植被指数的最佳回归模型以 S形曲线模型居多。 土壤调节因子 L,对 LAI 与 SAVI 及 SARVI 的曲
线回归模型影响较大,对 SAVI来说 L=0. 35,曲线拟合效果较好,对 SARVI而言, L=0. 1,曲线拟合效果最佳。
3. 2摇 逐步回归分析
从回归分析的角度来说,一般认为回归方程包含的自变量越多,回归平方和越大,剩余的平方和越小,剩
余均方也随之较小,预测值的误差也愈小,模拟的效果也就越好。 变量过多,建模的工作量就会增加;而且自
变量之间本身也存在相关性,因此在进行叶面积指数建模时,一是要尽可能选择变量之间相关性较小的因子,
使变量之间具有较强的独立性,二是需要进行自变量的筛选,选择适宜数量的变量。 在这种情况下可用逐步
回归分析,进行自变量因子的筛选。 对于逐步回归分析而言,首先要建立因变量 y 与自变量 x 之间的总回归
方程,再对总的方程及每—个自变量进行假设检验。 当总的方程不显著时,表明该多元回归方程线性关系不
成立;而当某—个自变量对 y影响不显著时,应该把它剔除,重新建立不包含该因子的多元回归方程。 筛选出
有显著影响的因子作为自变量,并建立“最优冶回归方程。 逐步回归分析结果见表 3 和表 4。
4877 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
表 3摇 逐步回归拟合结果
Table 3摇 Results of stepwise regression
模型
Model
引入变量
Variable R
2 修正的 R
2
Adj. R2
估计标准误差
Std. Error of the Estimate F
显著水平
Sig.
1 RVI 0. 73 0. 73 0. 46 187. 82 0. 000
2 RVI, SARVI0. 1 0. 83 0. 82 0. 37 162. 77 0. 000
3 RVI, SARVI0. 1, NDVI705 0. 85 0. 84 0. 35 124. 329 0. 000
4 RVI, SARVI0. 1, NDVI705, NDVI 0. 87 0. 87 0. 33 107. 050 0. 000
表 3 反映了逐步回归分析过程,共进行了 4 次拟合,最先引入回归方程的自变量为比值植被指数 RVI,其
次为 SARVI0. 1,NDVI705,NDVI。 可以看出 4 个被选择的自变量经过逐步回归都进入了回归方程。 随着模型
中变量的增加, R2、修正的 R2 随之增大,变量增至 4 个时,拟合方程的决定系数 R2 为 0. 87,修正的 R2 为
0郾 87,修正的 R2 和 R2 能较好的反映拟合优度;另一个方面随着模型中变量个数的增加,估计标准误差随之减
少,变量个数为 4 个时,估计标准误差为 0. 33。 回归方程包含不同的变量时虽然 F检验的值随着模型中变量
个数的增加而减少,但是其显著性概率值均小于 0. 001,即拒绝回归系数为 0 的假设,回归方程相关非常显
著。 因此最终的回归方程应该包括 RVI, SARVI0. 1,NDVI705 和 NDVI 等 4 个自变量,变量回归系数如表 4
所示。
表 4摇 逐步回归分析最终模型的统计量
Table 4摇 The final model statistic results of stepwise regression
变量
Variable
非标准的回归系数
Unstandardized coefficients
B Std. Error
标准回归系数
Standardized coefficients
Beta
T 显著水平Sig.
常数 Constant -1. 068 0. 219 -4. 879 0
RVI 0. 213 0. 034 0. 436 6. 354 0
SARVI01 1. 331 0. 514 0. 211 2. 592 0. 012
NDVI705 1. 144 0. 378 0. 198 3. 024 0. 004
NDVI 1. 652 0. 474 0. 234 3. 487 0
把表 4 中“非标准化回归系数冶栏目中的“B冶列数据代入多元回归模型得到最优的逐步回归模型。
Y= -1. 068+0. 213RVI+1. 331NDVI+1. 144NDVI705+1. 652SARVI0. 1
由回归方程式可以看出,研究区的叶面积指数 LAI与 RVI, SARVI0. 1, NDVI705, NDVI呈显著正相关。
3. 3摇 偏最小二乘回归分析
偏最小二乘法(Partial Least鄄Squares Regression, PLS)是一种新型的多元统计分析方法,集多元线性回归
分析、典型相关分析和主成分分析的功能和优点于一体[35]。 论文尝试使用偏最小二乘回归分析方法的变量
投影重要性指标(VIP)结合变量解释能力及变量权重来进行变量筛选,通过设置 3 次潜在因子的最大数量来
开展分析,其中模型效应设定为主效应,得到 3 个较好的回归分析结果(表 5)。
表 5摇 偏最小二乘回归分析结果
Table 5摇 Calibration and validation results of PLS regression
模型
Model
引入变量
Variables Entered
X方差累积解释量
Cumulative X Variance
Y方差累积解释量 R2
Cumulative Y Variance R2
Adj. R2
1 MNDVI, NDVI, RVI,NDVI705 0. 881 0. 890 0. 883
2 MNDVI, NDVI, RVI,NDVI705, SARVI0. 1 0. 903 0. 896 0. 887
3 MNDVI, NDVI, RVI, NDVI705, SARVI0. 1, SARVI0. 25 0. 934 0. 897 0. 887
5877摇 24 期 摇 摇 摇 孙华摇 等:三种回归分析方法在 Hyperion影像 LAI反演中的比较 摇
http: / / www. ecologica. cn
摇 摇 在整个回归分析过程中有 7 个因子的变量投影重要性指标小于 1 而被筛除(图 2)。 对剩余 6 个因子进
行回归建模,从表 5 可知,其实在潜在因子量达到 4 个以后,自变量的累积解释量和因变量的累积解释量的增
幅开始减少,最终得到 3 个偏最小回归模型:
Y= -1. 060+0. 484MNDVI+2. 092NDVI+0. 223RVI+1. 551 NDVI705
Y= -1. 113+0. 255MNDVI+1. 593NDVI+0. 207RVI+1. 145NDVI705+1. 297SARVI0. 1
Y= -1. 11+0. 291MNDVI+1. 57NDVI+0. 199RVI+1. 152NDVI705+2. 4SARVI0. 1-1. 072SARVI0. 25
图 2摇 变量投影重要性指标
Fig. 2摇 Variable Importance in the Projection
3. 4摇 反演模型的适用性检验
模型的适用性检验主要从两个方面对叶面积指数
的反演模型进行检验,即利用 42 个测试数据数据样地,
检验模型预测值与实际观测值是否存在较好的线性拟
合关系;其二则是计算了 3 种回归分析方法的 6 个反演
模型的残差、残差标准差 滓^ 、和模型评价值指标MEI,分
析残差图是否落在置信带内[ -2滓^ , 2滓^ ],分析结果如
图 3 和图 4 所示。
图 3摇 LAI预测结果与实测值比较
Fig. 3摇 Observed values and predicted values of LAI
从图 6 预测结果与实测值比较结果来看,3 种回归
分析方法所对应的 6 个反演模型的预测值与实际观测
值之间拟合度较高,存在较好的线性关系;其中以 5—6
个因子建立的偏最小二乘回归分析方法所得结果的决定系数(R2)最高,为 0. 84,逐步回归分析结果的决定系
6877 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
数次之,为 0. 83,以 4 个因子建立的偏最小二乘回归分析方法所得结果的决定系数为 0. 82;曲线估计中的线
性回归的决定系数最低为 0. 64,曲线估计中指数性回归的决定系数次之,为 0. 66,从建模因子的选择上看,比
值植被指数(RVI)出现在每一个反演模型中,是最敏感的建模因子,其次是 SARVI0. 1, NDVI705, NDVI;从回归
方法的选择上看,逐步回归与偏最小二乘回归分析方法要优于曲线估计所得结果,利用偏最小二乘回归开展
叶面积指数反演,自变量因子采用 5—6个因子所建模型预测结果与实测值的差异变化不大。 在自变量个数
相同时,逐步回归与偏最小二乘回归分析两种方法所得的反演模型的精度基本相近。
图 4摇 回归分析的残差图
Fig. 4摇 Residual results of three regression analysis methods
从图 7 中可知数据的残差绝大部分都在置信带内[-2滓^ , 2滓^ ], 4 因子的偏最小二乘回归模型的预测值
与实测值虽然有较好的拟合关系,但残差全部落在置信带正区间内,随机性较差。 残差分析表明,选用 5—6
个自变量因子进行 LAI建模是可靠的,建模方法以偏最小二乘回归所得模型的精度最高。
由表 6 可知,由 5—6个因子建立的偏最小二乘回归方程,对样点预测比较稳定且残差平方和最小,可见
偏最小二乘法建立的估测模型在三种回归分析方法比较中是最理想的。
4摇 结论与讨论
采用曲线估计、逐步回归分析和偏最小二乘回归分析方法对黄丰桥林场进行叶面积指数进行估测,结合
地面调查结果,分析了叶面积指数与植被指数的关系,建立了叶面积指数估算模型。 利用与残差相关的多种
参数和相关图形表达,不仅检验了建模精度,而且进行了模型的适用性检验,得出结论。
7877摇 24 期 摇 摇 摇 孙华摇 等:三种回归分析方法在 Hyperion影像 LAI反演中的比较 摇
http: / / www. ecologica. cn
表 6摇 3 种回归分析方法建模精度比较
Table 6摇 Modeling precision of three regression analysis methods
逐步回归分析
Stepwise regression
曲线估计
Curve estimation
Linear Exponential
偏最小二乘回归
Partial least鄄square regression
PLS4 PLS5 PLS6
滓^2R 0. 11 0. 211 0. 232 0. 118 0. 11 0. 10
滓^2F 0. 131 0. 272 0. 262 0. 133 0. 127 0. 123
MEI 0. 121 0. 242 0. 247 0. 126 0. 119 0. 112
(1)植被指数反映植被结构和覆盖状态信息,是体现绿色植被相对丰度和活性的辐射量值指标,在不同
领域中得到广泛应用,一般通过计算遥感影像不同光谱波段间的线性和非线性组合来实现,本次研究参与分
析的 13 个植被指数中,RVI与 LAI 的相关性最大,敏感性最高,其次是 SARVI0. 1, NDVI705, NDVI,SARVI0. 1,
SARVI0. 25,且 6 个指数都参与了叶面积指数反演建模,所得模型的拟合效果较好,说明利用植被指数来建立研
究区 LAI的遥感监测模型是可行的。
(2)建立了植被指数与叶面积指数的曲线估计、逐步回归及偏最小二乘回归分析方法的 6 个模型,对这
些模型的分析发现,植被指数与叶面积指数之间的最佳回归模型为偏最小二乘回归模型,模型预测值与实测
值拟合的决定系数 R2 的最大值为 0. 84,参与建模的植被指数因子为 6 个,线性回归模型的决定系数最低为
0. 64。 建模因子相同时,逐步回归分析与偏最小二乘回归所得结果基本接近。
(3)模型的适用性检验表明,5—6个因子建立的偏最小二乘回归分析模型和逐步回归模型的残差点分布
随机性较好,以 0 为基准线上下对称分布。 曲线估计方法所得的结果的残差分布图较差,说明单个因子建模
精度较差,表 6 中模型评价指标值的分析也得到了类似的结论。 认为建模因子在 5—6个之间所得的模型精
度较高。
(4)研究所选建模因子基本上为植被指数因子,建立的回归分析模型多为线性模型,有关非线性模型及
定性因子的引入,有待进一步研究。
References:
[ 1 ]摇 Huang M, Ji J J. The spatial鄄temporal distribution of leaf area index in China: a comparison between ecosystem modeling and remote sensing
reversion. Acta Ecologica Sinica, 2010, 30(11): 3057鄄3064.
[ 2 ] 摇 Wang X Q, Ma L Y, Jia Z K, X俟 X Y. Research and application advances in leaf area index. Chinese Journal of Ecology, 2005, 24 (5):
537鄄541.
[ 3 ] 摇 Wu T, Ni S X, Li Y M, Chen J. A comparison on the algorithms for retrieval of LAI based on gap fraction of vegetation canopy. Journal of Nanjing
Normal University: Natural Science, 2006, 29(1): 111鄄115.
[ 4 ] 摇 Tian Y C, Yang J, Yao X, Zhu Y, Cao W X. Quantitative relationships between hyper鄄spectral vegetation indices and leaf area index of rice.
Chinese Journal of Applied Ecology, 2009, 20(7): 1685鄄1690.
[ 5 ] 摇 Song S L, Li P X, Gong W, Wang J P. PLS wavelength selection by hyper鄄spectral remote sensing in rice. Geomatics and Information Science of
Wuhan University, 2010, 35(2): 219鄄223.
[ 6 ] 摇 Liang L, Yang M H, Zang Z. Improvement research of wheat leaf area index measurement using visible / near鄄infrared spectra. Laser and Infrared,
2010, 40(11): 1205鄄1210.
[ 7 ] 摇 Zhao X, Liu S H, Wang P J, Wang J D, Tian Z K. A method for inverting chlorophyll content of wheat using hyperspectral. Geography and Geo鄄
Information Science, 2004, 20(3): 36鄄39.
[ 8 ] 摇 Du X M, Cai T J, Ju C Y. Estimation of forest canopy closure by using partial least square regression. Chinese Journal of Applied Ecology, 2008,
19(2): 273鄄277.
[ 9 ] 摇 Zhang J, Zhang Q, He J M, Sun B. Estimation and analysis of grass biomass in desert by remote sensing. Acta Ecologica Sinica, 2007, 27(6):
2294鄄2302.
[10] 摇 Guo Z H, Peng S L, Wang B S. Estimating forest biomass in Western Guangdong using Landsat TM data. Acta Ecologica Sinica, 2002, 22(11):
8877 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
1832鄄1840.
[11] 摇 Gonsamo A, Pellikka P, King D J. Large鄄scale leaf area index inversion algorithms from high鄄resolution airborne imagery. International Journal of
Remote Sensing, 2011, 32(4): 3897鄄3916.
[12] 摇 Gonsamo A. Normalized sensitivity measures for leaf area index estimation using three鄄band spectral vegetation indices. International Journal of
Remote Sensing, 2011, 32(7): 2069鄄2080.
[13] 摇 Sun P S, Liu S R, Liu J T, Li C We, Lin Y, Jiang H. Derivation and validation of leaf area index maps using NDVI data of different resolution
satellite imageries. Acta Ecologica Sinica, 2006, 26(11): 3826鄄3834.
[14] 摇 Zhang W C, Zhong S, Hu S Y. Spatial scale transferring study on leaf area index derived from remotely sensed data in the Heihe River Basin,
China. Acta Ecologica Sinica, 2008, 28(6): 2495鄄2503.
[15] 摇 Chen X F, Chen J M, An S Q, Liu Y H, Fang X Q, Wang S M. Comparison of different atmospheric correction models in their effects on Landsat
TM estimation of forest leaf area index. Chinese Journal of Ecology, 2006, 25(7): 769鄄773.
[16] 摇 Yang G J, Huang W J, Wang J H, Xing Z R. Inversion of forest leaf area index calculated from multi鄄source and multi鄄angle remote sensing data.
Chinese Bulletin of Botany, 2010, 45(5): 566鄄578.
[17] 摇 Li K L, Jiang J J, Mao R Z, Ni S X. The modeling of vegetation through leaf area index by means of remote sensing. Acta Ecologica Sinica, 2005,
25(6): 1491鄄1496.
[18] 摇 Chen J, Ni S X, Li J J, Wu T. Scaling effect and spatial variability in retrieval of vegetation LAI from remotely sensed data. Acta Ecologica Sinica,
2006, 26(5): 1502鄄1508.
[19] 摇 Turner D P, Cohen W B, Kennedy R E, Fassnacht KS, Briggs JM. Relationships between leaf area index and landsat TM spectral vegetation
indices across Three Temperate Zone Sites. Remote Sensing of Environment, 1999, 70(1): 52鄄68.
[20] 摇 Huete A R, Liu H Q, Batchily K, van Leeuwen W. A comparison of vegetation indices over a global set of TM images for EOS鄄MODIS. Remote
Sensing of Environment, 1997, 59(3): 440鄄451.
[21] 摇 Cai B F, Shao X. Leaf area index retrieval based on remotely sensed data and Prospect+ Sail model. Remote Sensing for Land and Resource,
2007, (2): 39鄄43.
[22] 摇 Jacquemoud S. Inversion of the PROSPECT+SAIL canopy reflectance model from AVIRIS equivalent spectra: theoretical study. Remote Sensing of
Environment, 1993, 44(2 / 3): 281鄄292.
[23] 摇 Liu Z Y, Ma L L, Tang L L. Inversion of LAI based on SAIL model with multi鄄angle and multi鄄spectral remote sensing data. Arid Land
Geography, 2010, 33(1): 93鄄98.
[24] 摇 Chen Y H, Zhang W C, Yong B. Retrieving leaf area index using a neural network based on classification knowledge. Acta Ecologica Sinica,
2007, 27(7): 2785鄄2793.
[25] 摇 Zhu G L, Ju W M, Chen J M, Fang 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.
[26] 摇 Zhao C Y, Shen W H, Peng H H. Methods for determining canopy leaf area index of Picea Crassifolia Forest in QILIAN Mountain, CHINA.
Chinese Journal of Plant Ecology, 2009, 33(5): 860鄄869.
[27] 摇 Wang Z X, Liu C, Huete A. From AVHRR鄄NDVI to MODIS鄄EVI: Advances in vegetation index research. Acta Ecologica Sinica, 2003, 23(5):
979鄄987.
[28] 摇 Miura T, Huete A, Yoshioka H. An empirical investigation of cross鄄sensor relationships of NDVI and red / near鄄infrared reflectance using EO鄄 1
Hyperion data. Remote Sensing of Environment, 2006, 100(2): 223鄄236.
[29] 摇 Schlerfa M, Atzbergerb C, Hilla J. Remote sensing of forest biophysical variables using HyMap imaging spectrometer data. Remote Sensing of
Environment, 2005, 95(2): 177鄄194.
[30] 摇 Gong P, Pu R L, Biging G S, Larrieu M R. Estimation of forest leaf area index using vegetation indices derived from Hyperion hyperspectral data.
IEEE Transactions of Geoscience and Remote Sensing, 2003, 41(6): 1355鄄1362.
[31] 摇 Huete A R, Liu H Q, Batchily K, van Leeuwen W. A comparison of vegetation indices over a global set of TM images for EOS鄄MODIS. Remote
Sensing of Environment, 1997, 59(3): 440鄄451.
[32] 摇 Mobasheri M R, Rezaei Y, Valadan Zoej M J. A method in extracting vegetation quality parameters using hyperion images, with application to
precision farming. World Applied Sciences Journal, 2007, 2(5): 476鄄483.
[33] 摇 Liang S L. Quantitative Remote Sensing of Land Surfaces. New Jersey: John Wiley and Sons Inc, 2004: 273鄄275.
[34] 摇 Li C G, Zhao X W, Li C G. Theory and realization of estimating forest stock volume by remote sensing. Beijing: Science Press,2006:109鄄110.
[35] 摇 Wang H W, Wu Z B, Meng J. Partial Least鄄Squares Regression Linear and Nonlinear Methods. Beijing: National Defense Industry Press,
2006: 97鄄123.
9877摇 24 期 摇 摇 摇 孙华摇 等:三种回归分析方法在 Hyperion影像 LAI反演中的比较 摇
http: / / www. ecologica. cn
参考文献:
[ 1 ]摇 黄玫, 季劲钧. 中国区域植被叶面积指数时空分布———机理模型模拟与遥感反演比较. 生态学报, 2010, 30(11): 3057鄄3064.
[ 2 ] 摇 王希群, 马履一, 贾忠奎, 徐程扬. 叶面积指数的研究与应用进展. 生态学杂志, 2005, 24(5): 537鄄541.
[ 3 ] 摇 吴彤, 倪绍祥, 李云梅, 陈健. 由冠层孔隙度反演植被叶面积指数的算法比较. 南京师大学报: 自然科学版, 2006, 29(1): 111鄄115.
[ 4 ] 摇 田永超, 杨杰, 姚霞, 朱艳, 曹卫星. 高光谱植被指数与水稻叶面积指数的定量关系. 应用生态学报, 2009, 20(7): 1685鄄1690.
[ 5 ] 摇 宋沙磊, 李平湘, 龚威, 汪金平. 基于水稻高光谱遥感数据的 PLS波长选择研究. 武汉大学学报: 信息科学版, 2010, 35(2): 219鄄223.
[ 6 ] 摇 梁亮, 杨敏华, 臧卓. 利用可见 /近红外光谱测定小麦叶面积指数的改进研究. 激光与红外, 2010, 40(11): 1205鄄1210.
[ 7 ] 摇 赵祥, 刘素红, 王培娟, 王锦地, 田振坤. 基于高光谱数据的小麦叶绿素含量反演. 地理与地理信息科学, 2004, 20(3): 36鄄39.
[ 8 ] 摇 杜晓明, 蔡体久, 琚存勇. 采用偏最小二乘回归方法估测森林郁闭度. 应用生态学报, 2008, 19(2): 273鄄277.
[ 9 ] 摇 张杰, 张强, 何金梅, 孙冰. 荒漠草原牧草生物量的遥感估算及空间分布. 生态学报, 2007, 27(6): 2294鄄2302.
[10] 摇 郭志华, 彭少麟, 王伯荪. 利用 TM数据提取粤西地区的森林生物量. 生态学报, 2002, 22(11): 1832鄄1840.
[13] 摇 孙鹏森, 刘世荣, 刘京涛, 李崇巍, 林勇, 江洪. 利用不同分辨率卫星影像的 NDVI数据估算叶面积指数(LAI)———以岷江上游为例. 生
态学报, 2006, 26(11): 3826鄄3834.
[14] 摇 张万昌, 钟山, 胡少英. 黑河流域叶面积指数(LAI)空间尺度转换. 生态学报, 2008, 28(6): 2495鄄2503.
[15] 摇 陈新芳, 陈镜明, 安树青, 刘玉虹, 方秀琴, 王书明. 不同大气校正方法对森林叶面积指数遥感估算影响的比较. 生态学杂志, 2006, 25
(7): 769鄄773.
[16] 摇 杨贵军, 黄文江, 王纪华, 邢著荣. 多源多角度遥感数据反演森林叶面积指数方法. 植物学报, 2010, 45(5): 566鄄578.
[17] 摇 李开丽, 蒋建军, 茅荣正, 倪绍祥. 植被叶面积指数遥感监测模型. 生态学报, 2005, 25(6): 1491鄄1496.
[18] 摇 陈健, 倪绍祥, 李静静, 吴彤. 植被叶面积指数遥感反演的尺度效应及空间变异性. 生态学报, 2006, 26(5): 1502鄄1508.
[21] 摇 蔡博峰, 绍霞. 基于 PROSPECT+SAIL模型的遥感叶面积指数反演. 国土资源遥感, 2007, (2): 39鄄43.
[23] 摇 刘照言, 马灵玲, 唐伶俐. 基于 SAIL模型的多角度多光谱遥感叶面积指数反演. 干旱区地理, 2010, 33(1): 93鄄98.
[24] 摇 陈艳华, 张万昌, 雍斌. 基于分类知识利用神经网络反演叶面积指数. 生态学报, 2007, 27(7): 2785鄄2793.
[26] 摇 赵传燕, 沈卫华, 彭焕华. 祁连山区青海云杉林冠层叶面积指数的反演方法. 植物生态学报, 2009, 33(5): 860鄄869.
[27] 摇 王正兴, 刘闯, Huete A. 植被指数研究进展: 从 AVHRR鄄NDVI到 MODIS鄄EVI. 生态学报, 2003, 23(5): 979鄄987.
[34] 摇 李崇贵, 赵宪文, 李春干. 森林蓄积量遥感估测理论与实践. 北京: 科学出版社, 2006: 109鄄110.
[35] 摇 王惠文, 吴载斌, 孟洁. 偏最小二乘回归的线性与非线性方法. 北京: 国防工业出版社, 2006: 97鄄123.
0977 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
ACTA ECOLOGICA SINICA Vol. 32,No. 24 December,2012(Semimonthly)
CONTENTS
A bibliometric study of biodiversity research in China LIU Aiyuan, GUO Yuqing, LI Shiying,et al (7635)…………………………
Effects of elevated CO2 and nitrogen deposition on leaf nutrient quality of Fargesia rufa Yi
ZHOU Xianrong, WANG Jianhua, ZHANG Hong,et al (7644)
……………………………………………
……………………………………………………………………
Airborne pollen assemblages and their relationships with climate factors in the central Shaanxi Province of the Loess Plateau:
a case in Xiaheimugou, Luochuan County L譈 Suqing, LI Yuecong, XU Qinghai,et al (7654)…………………………………
Spatial and temporal change in ecological assets in the Yangtze River Delta of China 1995—2007
XU Xibao, CHEN Shuang, YANG Guishan (7667)
……………………………………
…………………………………………………………………………………
Evaluation and optimization of woodland ecological patterns for Qingdao based on the agent鄄based model
FU Qiang, MAO Feng, WANG Tianqing,et al (7676)
……………………………
……………………………………………………………………………
Interactive mechanism of service function of alpine rangeland ecosystems in Qinghai鄄Tibetan Plateau
LIU Xingyuan, LONG Ruijun, SHANG Zhanhuan (7688)
…………………………………
…………………………………………………………………………
Preliminary evaluation of air temperature reduction of urban green spaces in Beijing
ZHANG Biao, GAO Jixi, XIE Gaodi,et al (7698)
…………………………………………………
…………………………………………………………………………………
Resources metabolism analysis for the pulp and paper industry in Wuhan, China
SHI Xiaoqing,LI Xiaonuo,ZHAO Linjia,et al (7706)
………………………………………………………
………………………………………………………………………………
The characteristics and influential factors of direct carbon emissions from residential energy consumption: a case study of Lijiang
City, China WANG Danyin, TANG Mingfang, REN Yin, et al (7716)…………………………………………………………
Spatial targeting of payments for ecosystem services Based on SWAT Model and cost鄄benefit analysis
SONG Xiaoyu,LIU Yuqing,DENG Xiaohong,et al (7722)
…………………………………
…………………………………………………………………………
The wind tunnel test of plastic greenhouse and its surface wind pressure patterns
YANG Zaiqiang,ZHANG Bo,XUE Xiaoping,et al (7730)
………………………………………………………
…………………………………………………………………………
Population quantitative characteristics and dynamics of rare and endangered plant Davidia involucrata in Hunan Province
LIU Haiyang, JIN Xiaoling, SHEN Shouyun,et al (7738)
……………
…………………………………………………………………………
Phenotypic diversity in populations of germplasm resources of Rodgersia sambucifolia and related species
LI Pingping, MENG Hengling, CHEN Junwen,et al (7747)
……………………………
………………………………………………………………………
Effects of sand burial and seed size on seed germination, seedling emergence and growth of Caragana korshinskii Kom. (Fabaceae)
YANG Huiling, LIANG Zhenlei,ZHU Xuanwei,et al (7757)

………………………………………………………………………
Population鄄keeping mechanism of the parasitoid Dastarcus helophoroides (Coleoptera: Bothrideridae) of Massicus raddei
(Coleoptera: Cerambycidae) in oak forest YANG Zhongqi, TANG Yanlong, JIANG Jing,et al (7764)…………………………
Study of mingling based on neighborhood spatial permutation LOU Minghua, TANG Mengping, QIU Jianxi,et al (7774)……………
Comparison of three regression analysis methods for application to LAI inversion using Hyperion data
SUN Hua, JU Hongbo, ZHANG Huaiqing,et al (7781)
…………………………………
……………………………………………………………………………
Response of seed germination and seedling growth of Pinus koraiensis and Quercus mongolica to comprehensive action of warming
and precipitation ZHAO Juan, SONG Yuan, SUN Tao, et al (7791)……………………………………………………………
Impacts of water stored in sapwood Populus bolleana on its sap flux DANG Hongzhong, LI Wei,ZHANG Youyan,et al (7801)………
Dynamics of greenhouse gases emission and its impact factors by fire disturbance from Alnus sibirica forested wetland in
Xiaoxing忆an Mountains, Northeast China GU Han,MU Changcheng, ZHANG Bowen (7808)……………………………………
Different tide status and salinity alter stoichiometry characteristics of mangrove Kandelia candel seedlings
LIU Biner, LIAO Baowen, FANG Zhanqiang (7818)
……………………………
………………………………………………………………………………
Effects of shrub encroachment in desert grassland on runoff and the induced nitrogen loss in southeast fringe of Tengger Desert
LI Xiaojun, GAO Yongping (7828)
……
…………………………………………………………………………………………………
Community structure and throughfall erosivity characters of artificial rainforest in Xishuangbanna
DENG Yun, TANG Yanlin , CAO Min, et al (7836)
………………………………………
………………………………………………………………………………
Temporal鄄spatial variations of net ecosystem productivity in alpine area of southwestern China
PANG Rui,GU Fengxue,ZHANG Yuandong, et al (7844)
…………………………………………
…………………………………………………………………………
Relationships between chemical compositions of Quercus species seeds and climatic factors in temperate zone of NSTEC
LI Dongsheng, SHI Zuomin, LIU Shirong, et al (7857)
……………
……………………………………………………………………………
Effects of simulated acid rain stress on the PS域 reaction center and free radical metabolism in leaves of longan
LI Yongyu, PAN Tengfei, YU Dong, et al (7866)
……………………
…………………………………………………………………………………
Assessment of organic pollution for surface soil in Shenyang suburbs CUI Jian,DU Jizhong,MA Hongwei,et al (7874)………………
The impact of rainfall on soil respiration in a rain鄄fed maize cropland GAO Xiang, HAO Weiping, GU Fengxue, et al (7883)………
Effects of winter crops on enzyme activity and morphological characteristics of root in subsequent rice crops
YU Tianyi, PANG Huancheng,REN Tianzhi,et al (7894)
…………………………
…………………………………………………………………………
Dynamic changes of soil moisture and nitrate nitrogen in wheat and maize intercropping field under different nitrogen supply
YANG Ruiju, CHAI Shouxi, MA Zhongming (7905)
…………
………………………………………………………………………………
Characteristics of the bird diversity and the impact factors in Weishan Lake YANG Yuewei, LI Jiuen (7913)………………………
The effect of cropping landscapes on the population dynamics of the cotton bollworm Helicoverpa armigera (Lepidoptera,
Noctuidae) in the northern Xinjiang LU Zhaozhi, PAN Weilin, ZHANG Xin, et al (7925)……………………………………
The seasonal variations of nitrogen and phosphorus release and its fluxes from the sediments of the Beili Lake in the Hangzhou
West Lake LIU Jingjing,DONG Chunying,SONG Yingqi,et al (7932)……………………………………………………………
Optimization of lake model salmo based on real鄄coded genetic algorithm
GUO Jing, CHEN Qiuwen, ZHANG Xiaoqing, et al (7940)
………………………………………………………………
………………………………………………………………………
The influence of climatic environmental factors and fishing pressure on changes of hairtail catches in the northern South China
Sea WANG Yuezhong, SUN Dianrong, CHEN Zuozhi, et al (7948)………………………………………………………………
Seasonal and spatial distribution of acid volatile sulfide in sediment under different mariculture types in Nansha Bay, China
YAN Tingru, JIAO Haifeng, MAO Yuze, et al (7958)
…………
……………………………………………………………………………
Review and Monograph
Research progress on the mechanism of improving plant cold hardiness XU Chengxiang (7966)………………………………………
Influences of vegetation on permafrost: a review CHANG Xiaoli,JIN Huijun,WANG Yongping,et al (7981)…………………………
Home鄄field advantage of litter decomposition and its soil biological driving mechanism: a review
ZHA Tonggang, ZHANG Zhiqiang, SUN Ge, et al (7991)
………………………………………
…………………………………………………………………………
Research progress on the relationship of pollutants between road鄄deposited sediments and its washoff
ZHAO Hongtao, LI Xuyong, YIN Chengqing (8001)
…………………………………
………………………………………………………………………………
8008 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
《生态学报》2013 年征订启事
《生态学报》是中国生态学学会主办的生态学专业性高级学术期刊,创刊于 1981 年。 主要报道生态学研
究原始创新性科研成果,特别欢迎能反映现代生态学发展方向的优秀综述性文章;研究简报;生态学新理论、
新方法、新技术介绍;新书评介和学术、科研动态及开放实验室介绍等。
《生态学报》为半月刊,大 16 开本,300 页,国内定价 90 元 /册,全年定价 2160 元。
国内邮发代号:82鄄7,国外邮发代号:M670
标准刊号:ISSN 1000鄄0933摇 摇 CN 11鄄2031 / Q
全国各地邮局均可订阅,也可直接与编辑部联系购买。 欢迎广大科技工作者、科研单位、高等院校、图书
馆等订阅。
通讯地址: 100085 北京海淀区双清路 18 号摇 电摇 摇 话: (010)62941099; 62843362
E鄄mail: shengtaixuebao@ rcees. ac. cn摇 网摇 摇 址: www. ecologica. cn
摇 摇 编辑部主任摇 孔红梅摇 摇 摇 执行编辑摇 刘天星摇 段摇 靖
生摇 态摇 学摇 报
(SHENGTAI摇 XUEBAO)
(半月刊摇 1981 年 3 月创刊)
第 32 卷摇 第 24 期摇 (2012 年 12 月)
ACTA ECOLOGICA SINICA

(Semimonthly,Started in 1981)

Vol郾 32摇 No郾 24 (December, 2012)
编摇 摇 辑摇 《生态学报》编辑部
地址:北京海淀区双清路 18 号
邮政编码:100085
电话:(010)62941099
www. ecologica. cn
shengtaixuebao@ rcees. ac. cn
主摇 摇 编摇 冯宗炜
主摇 摇 管摇 中国科学技术协会
主摇 摇 办摇 中国生态学学会
中国科学院生态环境研究中心
地址:北京海淀区双清路 18 号
邮政编码:100085
出摇 摇 版摇
摇 摇 摇 摇 摇 地址:北京东黄城根北街 16 号
邮政编码:100717
印摇 摇 刷摇 北京北林印刷厂
发 行摇
地址:东黄城根北街 16 号
邮政编码:100717
电话:(010)64034563
E鄄mail:journal@ cspg. net
订摇 摇 购摇 全国各地邮局
国外发行摇 中国国际图书贸易总公司
地址:北京 399 信箱
邮政编码:100044
广告经营
许 可 证摇 京海工商广字第 8013 号
Edited by摇 Editorial board of
ACTA ECOLOGICA SINICA
Add:18,Shuangqing Street,Haidian,Beijing 100085,China
Tel:(010)62941099
www. ecologica. cn
Shengtaixuebao@ rcees. ac. cn
Editor鄄in鄄chief摇 FENG Zong鄄Wei
Supervised by摇 China Association for Science and Technology
Sponsored by摇 Ecological Society of China
Research Center for Eco鄄environmental Sciences, CAS
Add:18,Shuangqing Street,Haidian,Beijing 100085,China
Published by摇 Science Press
Add:16 Donghuangchenggen North Street,
Beijing摇 100717,China
Printed by摇 Beijing Bei Lin Printing House,
Beijing 100083,China
Distributed by摇 Science Press
Add:16 Donghuangchenggen North
Street,Beijing 100717,China
Tel:(010)64034563
E鄄mail:journal@ cspg. net
Domestic 摇 摇 All Local Post Offices in China
Foreign 摇 摇 China International Book Trading
Corporation
Add:P. O. Box 399 Beijing 100044,China
摇 ISSN 1000鄄0933
CN 11鄄2031 / Q
国内外公开发行 国内邮发代号 82鄄7 国外发行代号 M670 定价 70郾 00 元摇