香格里拉县森林生物量遥感估测研究
【摘要】:近年来,在全球气候变暖的背景下,森林碳储量研究成为了人们关注的焦点。而森林生物量的估算是进行陆地生态系统碳循环和碳储量及变化分析的基础,也是生态学和全球变化研究的重要内容之一。香格里拉县地处横断山脉,位于世界面积最大的滇西北“三江并流”世界自然遗产保护区腹地。研究这一区域森林植被的生物量/碳储量,对进一步开展高原陆地生态系统碳循环研究,科学评价自然遗产保护区的生态服务功能及其在全球气候变化中的作用,提高全社会对保护区生态价值的认识都具有重要意义。随着3S技术的出现和发展,为快速、无破坏地进行较大尺度森林生物量和碳储量的研究提供了可能。
本文以香格里拉县为研究区,以2008、2009年陆地卫星数据TM为信息源,结合2008的森林资源二类调查数据和本文研究中调查的186个样地数据,在国家基金项目“滇西北三江并流区森林生态系统碳储量遥感定量研究”(No:40861009)支持下开展了香格里拉森林生物量和碳储量的遥感估测研究,建立了香格里拉县4种主要森林类型的遥感生物量估测模型。利用样地调查方法按4个森林类型估算了林下灌木、草本和枯枝落叶层的生物量,与遥感模型估测的乔木层生物量汇总,完成了香格里拉县主要森林生态系统的森林生物量和碳储量的估测。由于研究地区地处高山峡谷区,地形深切割,高差大(从1503米-5545米),导致遥感数据的预处理和标准化难度大。针对此情况,论文重点在遥感数据的地形阴影去除、森林生物量遥感特征提取、森林生物量遥感建模方法等方面重点开展研究。主要研究结论如下:
(1)高山峡谷区遥感数据受地形影响较大,而本文采用的坡度匹配地形辐射校正方法较之常用的方法能更好地减小地形对遥感数据的影响,改善了遥感数据与森林蓄积量之间的相关性。
研究表明:仅采用辐射定标处理的各波段数据与森林蓄积量的相关性都很低,TM的6个波段与蓄积量相关都未通过显著性检验;经过大气校正后相关性有了部分改善,但仅有TM4的相关达到p=0.05显著水平;在大气校正的基础上进行的地形校正明显地改善了各波段与样地生物量之间的相关关系。通过0.05显著性检验的波段达到了3个波段,其中TM4的相关系数通过了0.01水平的显著性检验。
(2)与传统的遥感特征因子相比,本文提出的通过像元分解获取的森林阴影丰度特征与森林生物量的具有较好的相关性,并在3个树种的逐步回归建模中入选回归方程。
研究表明,像元分解后的阴影丰度图像与森林生物量的相关性较之单波段的相关性都有不同程度的提高,除云冷杉外,相关性都通过了0.01水平的显著性检验。表明像元分解得出的阴影丰度特征可以作为森林生物量遥感估测的重要因子。
(3)三种建模方法的结果表明,线性回归模型精度最低,神经网络模型精度次之,支持向量机回归模型(SVM)精度最高。
利用相同的的建模样本数据和检验数据,采用多元线性回归、神经网络和支持向量机方法进行了森林生物量遥感估测模型的建模实验。检验结果表明,多元线性回归的检验精度最低,神经网络次之,支持向量机精度最高,除云冷杉外,模型预测值的平均相对误差小于20%。
多元回归由于要求数据分布满足一些假设条件(如线性、正态性、等方差、独立性等),使模型中得以保留下的变量很少(2-3个),造成了遥感信息利用不充分,导致模型的估测效果较差。
神经网络方法可以得到小的训练误差,但对于建模外的新数据,其泛化能力较差,存在过学习问题,从而导致模型的估测精度不高。
支持向量机不同于神经网络和线性回归等传统方法以训练误差最小化作为优化目标,而是以训练误差作为优化问题的约束条件,以置信范围值最小化作为优化目标。因此,SVM的泛化能力要明显优越于神经网络等传统学习方法。另外,SVM的求解最后转化成二次规划问题的求解,因此,SVM的解是唯一的、也是全局最优的。本文的研究结果中支持向量机模型精度最高也充分证明了SVM理论的正确和有效性。
(4)研究区主要森林生态系统的生物量和碳储量
利用支持向量机回归模型估算了香格里拉县主要森林生态系统的乔木层生物量,利用样地调查法估算了林下层生物量和土壤碳,得到香格里拉县主要森林生态系统的总碳储量为302.984 TgC,其中,乔木层、灌木层、草本层、枯落物层和土壤层的碳储量分别为60.196 TgC、5.433 TgC、1.080 TgC、3.582 TgC、232.692 TgC,分别占总碳储量的19.87%、1.79%、0.36%、1.18%和76.80%,各层按照碳储量大小排序为:土壤层乔木层灌木层枯落物层草本层。
不同的森林类型生态系统的碳储量差异较大。研究区主要的森林生态系统按照碳储量大小排序为:云冷杉栎类高山松云南松。
(5)研究区主要森林生态系统碳密度
香格里拉县主要森林生态系统的平均碳密度为403.480 t/hm2。其中,云冷杉林的碳密度最大,为576.889 t/hm2,其次是栎类,为326.947 t/hm2,高山松和云南松分别为279.993 t/hm2、255.792 t/hm2。香格里拉县的平均碳密度403.480 t/hm2要高于临近地区估算的四川森林生态系统平均碳密度为232.81 t/hm2。这说明,位于世界最大的自然遗产保护区的香格里拉除了具有重要的自然、人文和生物多样性的保护价值外,还是该地区的一个重要碳库,对整个区域生态系统的稳定和平衡发挥着重要的作用。
【关键词】:森林生物量 TM 支持向量机 遥感估测模型 香格里拉县
【学位授予单位】:北京林业大学
【学位级别】:博士
【学位授予年份】:2012
【分类号】:S718.5
【目录】:
【学位授予单位】:北京林业大学
【学位级别】:博士
【学位授予年份】:2012
【分类号】:S718.5
【目录】:
- 摘要3-6
- ABSTRACT6-15
- 图目录15-17
- 表目录17-19
- 1 绪论19-36
- 1.1 森林生物量/碳储量研究的意义19-20
- 1.2 森林生物量/碳储量国内外研究概况20-33
- 1.2.1 生物量研究以大尺度和小尺度为主20-21
- 1.2.2 森林生物量/碳储量研究方法研究概述21-32
- 1.2.2.1 基于森林调查的森林生物量研究21-22
- 1.2.2.2 基于生理生态模型法的生物量研究22-25
- 1.2.2.3 基于遥感技术的森林生物量研究25-30
- 1.2.2.4 森林碳储量与碳循环研究30-32
- 1.2.3 存在的问题与展望32-33
- 1.3 研究目标与研究内容33-36
- 1.3.1 研究目标33
- 1.3.2 研究内容和研究方法33-34
- 1.3.4 技术路线34-36
- 2 研究区概况与试验数据36-44
- 2.1 地理位置36-37
- 2.2 自然条件37-39
- 2.2.1 地质地貌37
- 2.2.2 气候37
- 2.2.4 森林土壤37-38
- 2.2.5 森林植被类型38-39
- 2.3 社会经济状况39-40
- 2.4 数据源40-44
- 2.4.1 遥感数据源40
- 2.4.2 DEM数据40
- 2.4.3 森林资源二类调查数据40
- 2.4.4 样地的乔木生物量实测数据40-43
- 2.4.5 林下层植被生物量样地数据43
- 2.4.6 土壤碳储量样地数据43-44
- 3 研究区遥感图像预处理44-58
- 3.1 LANDSAT TM数据几何校正44-47
- 3.1.1 几何校正方法44-45
- 3.1.2 校正结果及精度评价45-46
- 3.1.3 遥感图像的剪切与拼接46-47
- 3.2 LANDSAT TM数据辐射校正47-57
- 3.2.1 Landsat TM辐射大气校正47-51
- 3.2.1.1 大气校正方法48-49
- 3.2.1.3 TM遥感图像大气校正结果及评价49-51
- 3.2.2 Landst TM地形辐射校正51-57
- 3.2.2.1 地形校正方法52-53
- 3.2.2.2 TM图像地形校正结果53-56
- 3.2.2.3 图像辐射校正对森林蓄积量与遥感因子相关性的影响56-57
- 3.3 本章小结57-58
- 4. 森林生物量遥感特征因子提取及筛选58-87
- 4.1 森林生物量遥感特征因子的提取58-78
- 4.1.1 树木阴影特征因子的提取58-71
- 4.1.1.1 树木阴影的形成58-60
- 4.1.1.2 混合像元分解模型60-61
- 4.1.1.3 基于混合像元分解的树木阴影信息提取61-71
- 4.1.2 植被指数因子提取71-74
- 4.1.2.1 简单植被指数72
- 4.1.2.2 多波段线性组合72-73
- 4.1.2.3 归一化植被指数73
- 4.1.2.4 复杂植被指数73-74
- 4.1.3 K-T因子提取74-75
- 4.1.4 K-L因子提取75-76
- 4.1.5 纹理特征因子的提取76-78
- 4.2 森林生物量遥感估测模型的变量优选78-86
- 4.2.1 备选变量的确定78-81
- 4.2.2 变量的优选81-86
- 4.2.2.1 变量优选方法81
- 4.2.2.2 变量优选结果81-86
- 4.7 本章小结86-87
- 5 森林生物量遥感模型建立87-133
- 5.1 研究思路87-89
- 5.2 多元线型回归生物量遥感模型89-92
- 5.2.1 建模样本和采用的遥感因子89-90
- 5.2.2 森林生物量回归建模结果90-91
- 5.2.3 模型有效性检验91-92
- 5.2.4 结论92
- 5.3 基于人工神经网络的森林生物量遥感模型92-100
- 5.3.1 BP人工神经网络原理与方法92-96
- 5.3.2 基于MATLAB的人工神经网络模型实现方法96
- 5.3.3 香格里拉县森林生物量估测神经网络模型的建立96-100
- 5.3.3.1 BP神经网络建模96-99
- 5.3.3.2 BP神经网络模型精度评价99-100
- 5.4 支持向量机森林生物量估测模型的建立100-129
- 5.4.1 支持向量机(Support Vector Machine,SVM)基本原理100-114
- 5.4.1.1 线性支持向量分类100-102
- 5.4.1.2 非线性支持向量分类102-105
- 5.4.1.3 ε-支持向量回归机105-112
- 5.4.1.4 ν-支持向量回归机112-113
- 5.4.1.5 支持向量机回归算法步骤113
- 5.4.1.6 支持向量机回归的特点113-114
- 5.4.3 支持向量机森林生物量遥感回归模型的建立114-129
- 5.4.3.1 Libsvm软件介绍115-116
- 5.4.3.2 建模特征116
- 5.4.3.3 建立样本训练集和检验集116
- 5.4.3.4 数据格式与归一化处理116-117
- 5.4.3.5 核函数选取117
- 5.4.3.6 最优参数选取及模型确定117-129
- 5.5 三种建模方法精度评价与模型选优129-131
- 5.6 本章小结131-133
- 6 森林生物量/碳储量估测结果及空间分异规律133-150
- 6.1 乔木生物量估测133-135
- 6.1.1 基于支持向量机的乔木层生物量估测133-134
- 6.1.2 乔木层生物量估测结果与分析134-135
- 6.2 林下植物生物量的估测135-136
- 6.2.1 灌木生物量测定方法135-136
- 6.2.2 草本及枯落物生物量测定方法136
- 6.2.3 林下灌草及枯落物生物量的估算136
- 6.3 森林生物量估测136-138
- 6.4 森林生态系统碳储量估测138-146
- 6.4.1 4个优势树种(组)含碳率测定138-140
- 6.4.1.1 样品的采集与处理138-139
- 6.4.1.2 测定方法139
- 6.4.1.3 含碳率测定结果139-140
- 6.4.2 研究区主要森林生态系统乔木层碳储量140-142
- 6.4.3 研究区主要森林生态系统灌草层和枯落物层碳储量142-143
- 6.4.4 研究区主要森林生态系统土壤层碳储量143-144
- 6.4.5 研究区主要森林生态系统总碳储量144-146
- 6.5 研究区主要森林生态系统碳储量空间格局146-148
- 6.5.1 研究区主要森林生态系统碳储量垂直分布特征146-147
- 6.5.2 研究区主要森林生态系统碳储量坡度分布特征147
- 6.5.3 研究区主要森林生态系统碳储量坡向分布特征147-148
- 6.6 本章小结148-150
- 7 结论与讨论150-155
- 7.1 结论150-152
- 7.2 讨论152-153
- 7.2.1 关于模型精度152-153
- 7.2.2 关于建模方法153
- 7.3 创新与展望153-155
- 参考文献155-166
- 导师简介166-167
- 个人简介167-168
- 读博期间发表文章168-169
- 致谢169
| 【引证文献】 | ||
|
|||
|
|||
| 【参考文献】 | ||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||
|
|||||
| 【共引文献】 | ||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
| 【同被引文献】 | ||
|
|||||||||||||
|
|||||||||||||
| 【二级参考文献】 | ||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||
|
|||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
| 【相似文献】 | ||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||
|
|||||||||
|
|||
|
|||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||
|
|||||||||||||||||||||



