一种气候变化影响下沼泽湿地植被覆盖度估算方法技术

一种气候变化影响下沼泽湿地植被覆盖度估算方法，本发明专利技术涉及一种气候变化影响下沼泽湿地植被覆盖度估算方法。本发明专利技术的目的是为了解决现有方法缺少对大的区域尺度模拟及未来预测的问题。方法为：获取及预处理遥感影像数据和地表观测气象数据；提取时间段内未变化沼泽湿地分布作为研究区范围；对气象数据进行插值，并对结果进行重采样；利用像元二分模型计算逐像元植被覆盖度；提取研究区所有像元对应的植被覆盖度及气象要素数值；利用多元逐步回归构建气候变化影响下沼泽湿地植被覆盖度估算模型。本发明专利技术利用了遥感数据宏观性及长期监测的优点，结合地表观测数据模拟了气候变化影响下沼泽湿地植被覆盖度。本发明专利技术应用于沼泽湿地植被覆盖度估算领域。

【技术实现步骤摘要】
一种气候变化影响下沼泽湿地植被覆盖度估算方法
本专利技术涉及一种气候变化影响下沼泽湿地植被覆盖度估算方法。
技术介绍
湿地是世界上最重要的自然生态系统之一，湿地生态系统在保护生物多样性，控制环境质量，服务社会经济可持续发展等方面起到重要作用。沼泽湿地作为一种非常重要的湿地类型，不仅可以促进水资源储存和净化，调节区域气候，还能为野生动物的繁殖提供栖息地。植被作为沼泽湿地生态系统的核心组成部分，能够将大气圈、水圈和土壤圈紧密联系在一起。沼泽湿地植被在沼泽湿地生态系统中的能量交换、生物地球化学循环和水文循环中起着重要的作用。植被覆盖度是指植被(包括叶、茎、枝)在地面的垂直投影面积占统计区域总面积的百分比，是衡量地表植被状况一个重要的指标，并已成为监测植被动态的重要生物物理参数。当前，有关沼泽湿地植被覆盖度的获取主要有两种方法，一种是地表实测法，主要借助植被覆盖度测量仪来开展地面实地观测；第二种方法是传统的模型估算法，主要借助植被归一化指数等遥感影像数据，利用模型方法计算地表植被覆盖度。第一种地表实测法，由于受实地观测的时间和空间尺度的限制，仅能在相对较小的区域或样点开展，且很难进行长时间监测。第二种传统模型估算法，由于是利用的历史植被指数遥感影像数据，很难对未来区域植被覆盖度变化做出预测。地表植被覆盖度主要受人类活动和气候变化的影响。在人类活动影响下，明确气候变化对沼泽湿地植被覆盖度的影响对预测未来区域沼泽湿地植被覆盖度具有重要的意义。当前有关探究气候变化对沼泽湿地植被覆盖度影响的研究较少，同时缺少在像元尺度上对未来气候变化影响下的区域植被覆盖度做出估算。专利技术内容本专利技术的目的是为了解决现有关于沼泽湿地植被覆盖度观测尺度较小，且无法对未来做出预测的问题，提出一种气候变化影响下沼泽湿地植被覆盖度估算方法。本专利技术一种气候变化影响下沼泽湿地植被覆盖度估算方法包括以下步骤：步骤一、获取覆盖研究区域的中等分辨率成像光谱仪NDVI数据集、地面气象观测数据集、以及研究时间段起止年份的土地利用/覆被数据，并进行数据预处理；步骤二、根据研究时间段起止年份的两期土地利用/覆被数据，提取研究时间段内未发生土地利用/覆被变化的沼泽湿地分布，并将其作为研究区范围；步骤三、采用克里格插值法，将地面气象观测数据集进行插值，得到研究时间段内各气象要素数值空间分布数据集；步骤四、分别将各气象要素数值空间分布数据集重采样为与NDVI数据集相同分辨率的长时间序列气象要素数据集；步骤五、从步骤一NDVI数据集中提取研究区内所有像元逐旬NDVI数值，分别计算各像元累积频率5％和95％的NDVI数值，作为NDVI阈值；步骤六、采用像元二分模型，根据步骤五确定的NDVI阈值，利用NDVI数据集计算逐像元植被覆盖度，得到长时间序列植被覆盖度数据集；步骤七、根据步骤二的未发生土地利用/覆被变化的沼泽湿地分布，从步骤六得到的长时间序列植被覆盖度数据集和步骤四得到的长时间序列气象要素数据集中提取研究区范围内所有像元逐一对应的植被覆盖度值及各气象要素值；步骤八、根据提取出的所有像元逐一对应的植被覆盖度值及各气象要素值，通过多元逐步回归，构建气候变化影响下逐像元植被覆盖度估算模型，利用该模型估算气候变化影响下沼泽湿地植被覆盖度。本专利技术的有益效果是：在传统的植被覆盖度模型估算方法的基础上，提出了一种构建像元尺度上的气候变化与沼泽湿地植被覆盖度定量关系模型，来估算气候变化影响下的沼泽湿地植被覆盖度的方法。该方法改善了传统模型估算法不能对未来气候变化影响下沼泽湿地植被覆盖做出预测的缺点。本专利技术既能够在高的空间分辨率上，对区域沼泽湿地植被覆盖度做出估算，又能够预测未来气候变化影响下的沼泽湿地植被覆盖度。本专利技术对于利用长时间序列的遥感植被数据及地面观测气象数据，在区域尺度上，预测未来沼泽湿地植被覆盖提供了新的方法。附图说明图1为本专利技术的流程图；图2为实施例一中未变化沼泽湿地空间分布图；图3为实施例一中沼泽湿地植被覆盖度估算值及实地观测点图；其中星形图标为实地观测点；图4为实施例一中沼泽湿地植被覆盖度实测值与估算值比较图。具体实施方式本专利技术技术方案不局限于以下所列举的具体实施方式，还包括各具体实施方式之间的任意组合。具体实施方式一：本实施方式一种气候变化影响下沼泽湿地植被覆盖度估算方法包括以下步骤：步骤一、获取覆盖研究区域的中等分辨率成像光谱仪NDVI数据集、地面气象观测数据集、以及研究时间段起止年份的土地利用/覆被数据，并进行数据预处理；步骤二、根据研究时间段起止年份的两期土地利用/覆被数据，提取研究时间段内未发生土地利用/覆被变化的沼泽湿地分布，并将其作为研究区范围；步骤三、采用克里格插值法，将地面气象观测数据集进行插值，得到研究时间段内各气象要素数值空间分布数据集；步骤四、分别将各气象要素数值空间分布数据集重采样为与NDVI数据集相同分辨率的长时间序列气象要素数据集；步骤五、从步骤一NDVI数据集中提取研究区内所有像元逐旬NDVI数值，分别计算各像元累积频率5％和95％的NDVI数值，作为NDVI阈值；步骤六、采用像元二分模型，根据步骤五确定的NDVI阈值，利用NDVI数据集计算逐像元植被覆盖度，得到长时间序列植被覆盖度数据集；步骤七、根据步骤二的未发生土地利用/覆被变化的沼泽湿地分布，从步骤六得到的长时间序列植被覆盖度数据集和步骤四得到的长时间序列气象要素数据集中提取研究区范围内所有像元逐一对应的植被覆盖度值及各气象要素值；步骤八、根据提取出的所有像元逐一对应的植被覆盖度值及各气象要素值，通过多元逐步回归，构建气候变化影响下逐像元植被覆盖度估算模型，利用该模型估算气候变化影响下沼泽湿地植被覆盖度。本实施方式的有益效果是：在传统的植被覆盖度模型估算方法的基础上，提出了一种构建像元尺度上的气候变化与沼泽湿地植被覆盖度定量关系模型，来估算气候变化影响下的沼泽湿地植被覆盖度的方法。该方法改善了传统模型估算法不能对未来气候变化影响下沼泽湿地植被覆盖做出预测的缺点。本实施方式既能够在高的空间分辨率上，对区域沼泽湿地植被覆盖度做出估算，又能够预测未来气候变化影响下的沼泽湿地植被覆盖度。本实施方式对于利用长时间序列的遥感植被数据及地面观测气象数据，在区域尺度上，预测未来沼泽湿地植被覆盖提供。具体实施方式二：本实施方式与具体实施方式一不同的是：步骤一中数据预处理过程为：对NDVI植被指数及土地利用/覆被数据集进行投影转换，将其统一到同一坐标系和投影之下，再运用最大值合成法将旬NDVI数据合成为月NDVI值。其它与具体实施方式一相同。具体实施方式三：本实施方式与具体实施方式一或二不同的是：步骤六中利用像元二分模型，计算逐像元植被覆盖度公式如下：其中FVC为植被覆盖度计算数值，NDVI为研究时间段的生长季内NDVI平均值，NDVIsoil和NDVIveg分别为累积频率为5％和95％的NDVI数值。其它与具体实施方式一或二相同。具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：步骤七中的数据提取方法为：(1)、利用未变化沼泽湿地分布图截取植被覆盖度影像数据；(2)、将截取后的任一植被覆盖度栅格数据转换为点图层；(3)、利用点图层提取插值并重采样后的气象数据集中各气象要素本文档来自技高网...

【技术保护点】
1.一种气候变化影响下沼泽湿地植被覆盖度估算方法，其特征在于气候变化影响下沼泽湿地植被覆盖度估算方法包括以下步骤：步骤一、获取覆盖研究区域的中等分辨率成像光谱仪NDVI数据集、地面气象观测数据集、以及研究时间段起止年份的土地利用/覆被数据，并进行数据预处理；步骤二、根据研究时间段起止年份的两期土地利用/覆被数据，提取研究时间段内未发生土地利用/覆被变化的沼泽湿地分布，并将其作为研究区范围；步骤三、采用克里格插值法，将地面气象观测数据集进行插值，得到研究时间段内各气象要素数值空间分布数据集；步骤四、分别将各气象要素数值空间分布数据集重采样为与NDVI数据集相同分辨率的长时间序列气象要素数据集；步骤五、从步骤一NDVI数据集中提取研究区内所有像元逐旬NDVI数值，分别计算各像元累积频率5％和95％的NDVI数值，作为NDVI阈值；步骤六、采用像元二分模型，根据步骤五确定的NDVI阈值，利用NDVI数据集计算逐像元植被覆盖度，得到长时间序列植被覆盖度数据集；步骤七、根据步骤二的未发生土地利用/覆被变化的沼泽湿地分布，从步骤六得到的长时间序列植被覆盖度数据集和步骤四得到的长时间序列气象要素数据集中提取研究区范围内所有像元逐一对应的植被覆盖度值及各气象要素值；步骤八、根据提取出的所有像元逐一对应的植被覆盖度值及各气象要素值，通过多元逐步回归，构建气候变化影响下逐像元植被覆盖度估算模型，利用该模型估算气候变化影响下沼泽湿地植被覆盖度。...

【技术特征摘要】
1.一种气候变化影响下沼泽湿地植被覆盖度估算方法，其特征在于气候变化影响下沼泽湿地植被覆盖度估算方法包括以下步骤：步骤一、获取覆盖研究区域的中等分辨率成像光谱仪NDVI数据集、地面气象观测数据集、以及研究时间段起止年份的土地利用/覆被数据，并进行数据预处理；步骤二、根据研究时间段起止年份的两期土地利用/覆被数据，提取研究时间段内未发生土地利用/覆被变化的沼泽湿地分布，并将其作为研究区范围；步骤三、采用克里格插值法，将地面气象观测数据集进行插值，得到研究时间段内各气象要素数值空间分布数据集；步骤四、分别将各气象要素数值空间分布数据集重采样为与NDVI数据集相同分辨率的长时间序列气象要素数据集；步骤五、从步骤一NDVI数据集中提取研究区内所有像元逐旬NDVI数值，分别计算各像元累积频率5％和95％的NDVI数值，作为NDVI阈值；步骤六、采用像元二分模型，根据步骤五确定的NDVI阈值，利用NDVI数据集计算逐像元植被覆盖度，得到长时间序列植被覆盖度数据集；步骤七、根据步骤二的未发生土地利用/覆被变化的沼泽湿地分布，从步骤六得到的长时间序列植被覆盖度数据集和步骤四得到的长时间序列气象要素数据集中提取研究区范围内所有像元逐一对应的植被覆盖度值及各气象要素值；步骤八、根据提取出的所有像元逐一对应...

【专利技术属性】
技术研发人员：神祥金，姜明，吕宪国，
申请(专利权)人：中国科学院东北地理与农业生态研究所，
类型：发明
国别省市：吉林,22