一种人工林森林年龄空间制图的新方法与流程

本发明涉及人工林森林年龄空间制图领域，特别涉及一种人工林森林年龄空间制图的新方法。

背景技术：

联合国2013年发布的第五次《气候变化评估报告决策者摘要》再次确认了气候变暖的事实，并预测到本世纪末地表温度可能比1850年至1900年增长1.5至2摄氏度。联合国政府间气候变化委员会(ipcc)第五次评估报告更进一步警示世人气候变化的严重性和气候治理的紧迫性。在应对气候变化的研究中，生态系统碳储量和碳循环模式被认为对气候变化具有决定性影响，从而成为各国政府和学术界重点关注的最紧迫的环境问题。生态系统中积聚的碳储量的消长分别对应着大气中co2的浓度高低变化，这直接影响全球气温变化的速率和幅度。地球上碳库存在的5大生态系统(森林、湿地、灌木林/草地、农田和水生生态系统)中，森林是陆地生态系统中最大的碳库，其单位面积的碳储量是农田的20至100倍。现今，发展可靠的方法来充分理解森林干扰及恢复模式变化对碳蓄积能力、循环模式及气候变化的影响已成为众多科学计划的优先领域，并且这些计划的长期目标是降低全球和区域尺度碳核算的不确定性以及改进对大气碳源、碳汇动态特征的理解。

然而由于森林结构参数、森林干扰模式和强度存在广泛的时空异质性，以及估算模型和方法内在的缺陷，当今的森林碳估算面临诸多挑战和不确定性。利用森林清查数据记录的森林种植或干扰时间或年轮数据可推算林龄，而取森林清查数据或者年轮数据通常较为困难，费时费力，其对于天然林更为困难。

技术实现要素：

本发明所要解决的技术问题是提供一种人工林森林年龄空间制图的新方法，以解决现有技术中导致的上述多项缺陷。

为实现上述目的，本发明提供以下的技术方案：一种人工林森林年龄空间制图的新方法，包括如下步骤：

(1)遥感数据源选择

landsat影像按照下载地址http://glovis.usgs.gov/提示安装java，然后输入目标地的经纬度坐标或者条带号，并按年份逐步筛选影像，选择productlevel1下载数据；

(2)数据标准化处理方法

接着将ledaps数据预处理的运行环境设置为linuxubuntu，按照数据标准化处理方法将数据分为ledaps可用的标准数据处理和非标准数据格式处理；

(3)ledaps大气校正方法

然后对ledaps进行校正，所述ledaps校正方法含四个具体操作命令，分别为：lndpm、lndcal、lndcsm和lndsr；

(4)高程数据

目标地的dem文件，下载地址：

ftp://ftp.glcf.umd.edu/glcf/srtm/ft/ft0012.filled_finished-03sec-wrs2-global/glcf.tsm.b4-004.00.filled_finished-03sec-wrs2-global/wrs2_tiles/p000/srtm_ff03_p000r000/，并在envi软件中保存为envi标准文件格式，并命名为dem000000_utm。(000需替换为目标地条带号)；

(5)土地覆盖数据

建立目标地的土地覆盖数据(landcover)文件，从下载的源影像中挑选近期且质量较好的一幅影像进行有监分类，分类方法选择svm方法，分类类别为4大类：城市、农田、森林、水体；将分类好的影像在erdas或arcgis软件中修改类码，类码修改参照土地覆盖分类标准，城市类码调整为21，农田调整为82，森林为41，水体为11，将调整好类码的影像重命名为landcover01000000_utm(01，代表土地覆盖数据库的发展年份；前三个000为pathnumber，后三个000为rownumber，其余保持不变)；

(6)森林变化制图

将统一命名后的ledaps所得地表反射率图像，dem数据，landcover数据在基于单一影像掩膜与归一化、时间序列分析和森林干扰产品测算原理的基础上编译形成vct自动化模型，得到的森林干扰年份和干扰量结果产品，建立森林干扰和恢复数据集。根据已建立的森林干扰和恢复数据集，依据不同年龄森林的生物量累积模式或特征差异，将森林划分为3个龄组：幼龄林、中龄林和成熟龄，最后将每年的成熟林、中龄林和幼龄林合并得到林龄分布图；

(7)基于光谱波段以及第一主成分波段的相关性分析

然后通过对辐射定标、大气校正以及地形校正后的遥感影像(特别是森林样方位置上landsat地表反射率时间序列值)进行相应波段以及特征值的提取从而进行相关性分析获得其与由地面清查数据得到的森林年龄值之间的关联性；

(8)森林年龄建模

首先，通过随机抽样的方法将样本中70％的数据作为训练样本采用随机森林算法进行森林年龄估算建模，利用剩下30％的样本进行模型的精度验证，通过建模可以获得森林年龄与各因子之间的定量化关联，然后将构建的模型应用于完整的研究区域完成研究区的森林年龄制图。

优选的，所述步骤(2)中linuxubuntu运行环境设置包括：

a)ubuntu系统home下工作空间设置；

b)ubuntu下设置环境变量，gcc，gfortran编译器安装，配置源163.com，4.6.3，4.6，以找到模型编译所需下载站点；

c)安装ledaps校正模型。

优选的，所述步骤(6)中单一影像掩膜与归一化的原理为在landsat时间序列堆栈中的影像需要进行水体、云、云阴影掩膜来识别森林类别，识别的森林类别被用于计算类似森林的光谱指数。

优选的，所述步骤(6)中时间序列分析的原理为基于森林指数和掩膜结果形成时间序列轨迹用于生成森林变化产品，主要的两个步骤概述为时间序列中的每个影像用于创建掩膜并计算森林光谱指数，获得的掩膜结果和指数用于时间序列的干扰制图。

优选的，所述步骤(6)中森林干扰产品测算为每个监测到的干扰计算两个要素，即森林干扰年份和干扰量。

优选的，所述步骤(6)中幼龄林、中龄林和成熟龄的划分方法为核查森林干扰产品和森林恢复产品，在landsat数据覆盖的30年间，持续保持森林状态的persistingforest类别所在空间像元位置被识别为本研究中“成熟林”；在30年间，被识别为“disturbedforest”类别所在的空间像元位置被标定为本研究中的“幼龄林”；早于我们设定的30年间隔的起始年份发生干扰事件，若干年后恢复了其森林特征光谱的那些空间像元被识别为“中龄林”。

采用以上技术方案的有益效果是：本发明提供的一种人工林森林年龄空间制图的新方法，具有以下优点：(1)基于森林干扰和恢复产品，提出了基于三层结构的森林年龄和森林生物量估算理论和方法体系，(2)针对导致当前大多数碳计量模型估计值变动不一的两大误差源：森林干扰和恢复历史的信息不全和缺乏生物量密度数据，本研究给出了完整解决方案，这使得本研究的结果相对于现存的区域尺度的研究结果而言，时空分辨率上有明显进步，市场前景广阔。

附图说明

图1是人工林龄组划分流程图；

图2是林龄划分原理图。

具体实施方式

下面详细说明本发明的优选实施方式。

实施例：

一种人工林森林年龄空间制图的新方法，包括如下步骤：

(1)遥感数据源选择

landsat影像按照下载地址http://glovis.usgs.gov/提示安装java，然后输入目标地的经纬度坐标或者条带号，并按年份逐步筛选影像，选择productlevel1下载数据；

(2)数据标准化处理方法

接着将ledaps数据预处理的运行环境设置为linuxubuntu，按照数据标准化处理方法将数据分为ledaps可用的标准数据处理和非标准数据格式处理，所述linuxubuntu运行环境设置包括：

a)ubuntu系统home下工作空间设置；

b)ubuntu下设置环境变量，gcc，gfortran编译器安装，配置源163.com，4.6.3，4.6，以找到模型编译所需下载站点；

c)安装ledaps校正模型；

(3)ledaps大气校正方法

然后对ledaps进行校正，所述ledaps校正方法含四个具体操作命令，分别为：lndpm、lndcal、lndcsm和lndsr；

(4)高程数据

目标地的dem文件，下载地址：

ftp://ftp.glcf.umd.edu/glcf/srtm/ft/ft0012.filled_finished-03sec-wrs2-global/glcf.tsm.b4-004.00.filled_finished-03sec-wrs2-global/wrs2_tiles/p000/srtm_ff03_p000r000/，并在envi软件中保存为envi标准文件格式，并命名为dem000000_utm。(000需替换为目标地条带号)；

(5)土地覆盖数据

建立目标地的土地覆盖数据(landcover)文件，从下载的源影像中挑选近期且质量较好的一幅影像进行有监分类，分类方法选择svm方法，分类类别为4大类：城市、农田、森林、水体；将分类好的影像在erdas或arcgis软件中修改类码，类码修改参照土地覆盖分类标准，城市类码调整为21，农田调整为82，森林为41，水体为11，将调整好类码的影像重命名为landcover01000000_utm(01，代表土地覆盖数据库的发展年份；前三个000为pathnumber，后三个000为rownumber，其余保持不变)；

(6)森林变化制图

将统一命名后的ledaps所得地表反射率图像，dem数据，landcover数据在基于单一影像掩膜与归一化、时间序列分析和森林干扰产品测算原理的基础上编译形成vct自动化模型，得到的森林干扰年份和干扰量结果产品，建立森林干扰和恢复数据集。根据已建立的森林干扰和恢复数据集，依据不同年龄森林的生物量累积模式或特征差异，将森林划分为3个龄组：幼龄林、中龄林和成熟龄，最后将每年的成熟林、中龄林和幼龄林合并得到林龄分布图，所述单一影像掩膜与归一化的原理为在landsat时间序列堆栈中的影像需要进行水体、云、云阴影掩膜来识别森林类别，识别的森林类别被用于计算类似森林的光谱指数，时间序列分析的原理为基于森林指数和掩膜结果形成时间序列轨迹用于生成森林变化产品，主要的两个步骤概述为时间序列中的每个影像用于创建掩膜并计算森林光谱指数，获得的掩膜结果和指数用于时间序列的干扰制图，所述森林干扰产品测算为每个监测到的干扰计算两个要素，即森林干扰年份和干扰量，所述幼龄林、中龄林和成熟龄的划分方法为核查森林干扰产品和森林恢复产品，在landsat数据覆盖的30年间，持续保持森林状态的persistingforest类别所在空间像元位置被识别为本研究中“成熟林”；在30年间，被识别为“disturbedforest”类别所在的空间像元位置被标定为本研究中的“幼龄林”；早于我们设定的30年间隔的起始年份发生干扰事件，若干年后恢复了其森林特征光谱的那些空间像元被识别为“中龄林”；

(7)基于光谱波段以及第一主成分波段的相关性分析

然后通过对辐射定标、大气校正以及地形校正后的遥感影像(特别是森林样方位置上landsat地表反射率时间序列值)进行相应波段以及特征值的提取从而进行相关性分析获得其与由地面清查数据得到的森林年龄值之间的关联性；

(8)森林年龄建模

首先，通过随机抽样的方法将样本中70％的数据作为训练样本采用随机森林算法进行森林年龄估算建模，利用剩下30％的样本进行模型的精度验证，通过建模可以获得森林年龄与各因子之间的定量化关联，然后将构建的模型应用于完整的研究区域完成研究区的森林年龄制图。

本发明提供了一种人工林森林年龄空间制图的新方法，具有以下优点：(1)基于森林干扰和恢复产品，提出了基于三层结构的森林年龄和森林生物量估算理论和方法体系，(2)针对导致当前大多数碳计量模型估计值变动不一的两大误差源：森林干扰和恢复历史的信息不全和缺乏生物量密度数据，本研究给出了完整解决方案，这使得本研究的结果相对于现存的区域尺度的研究结果而言，时空分辨率上有明显进步，市场前景广阔。

以上所述的仅是本发明的优选实施方式，应当指出，对于本领域的普通技术人员来说，在不脱离本发明创造构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。