一种草本沼泽植被物候信息提取方法

1.本发明涉及一种草本沼泽植被物候信息提取方法，特别是涉及一种基于遥感数据提取区域草本沼泽植被物候信息的方法。


背景技术：

2.湿地和森林、海洋是世界三大生态系统，湿地在全球碳循环过程中扮演着极其重要的角色。草本沼泽作为一种重要的湿地类型，在调节区域气候、净化环境等方面具有重要作用。植被作为草本沼泽生态系统的重要组成部分，在沼泽生态系统的能量转化和循环中发挥重要作用。植被物候是植物随着气候的季节性变化而发生萌芽、抽枝、展叶、开花、结果及落叶等规律性变化的现象。植被物候是植物长期适应季节变化而形成的生长发育节律，是指示植被生长发育和生理活动最常用的一个指标。在全球变化背景下，草本沼泽植被物候发生了明显变化，准确获取草本沼泽植被物候信息，明确草本沼泽植被物候时空变化是评估草本沼泽固碳潜力的关键。
3.当前，有关草本沼泽植被物候信息的获取主要通过实地观测法，通过定点观测，获取草本沼泽植被的物候期。这种观测方法虽然在小范围内较为精确，但会耗费大量的人力和物力，无法进行连续观测，且无法获取整个区域的草本沼泽植被物候信息。


技术实现要素：

4.本发明的目的是为了解决现有草本沼泽植被物候实地观测范围小、周期长、无法快速准确地提取区域草本沼泽植被物候信息的问题，提供了一种草本沼泽植被物候信息提取方法。
5.一种草本沼泽植被物候信息提取方法按照以下步骤进行：
6.步骤一、获取研究时间段内覆盖研究区域的逐旬归一化植被指数ndvi数据集、以及两期沼泽湿地分布数据，并进行数据预处理；
7.步骤二、利用最大值合成法将逐旬ndvi数据集合成为逐月ndvi数据集，并利用算术平均法计算得到多年平均生长季ndvi；
8.步骤三、根据两期沼泽湿地分布数据，提取研究时间段内未变化的草本沼泽分布；
9.步骤四、在未变化的草本沼泽分布范围内，提取多年平均生长季ndvi≥0.1的像元分布，得到未变化的草本沼泽植被分布，并将此作为研究区；
10.步骤五、利用研究区裁切逐旬ndvi数据集，提取出逐像元的每年逐旬ndvi值，并利用算术平均法计算得到逐像元的多年平均逐旬ndvi值；
11.步骤六、计算逐像元的多年平均逐旬ndvi值的变化率，并确定最大变化率对应的旬数和最小变化率对应的旬数；
12.步骤七、将最大变化率对应的旬数作为物候开始期对应的旬数，并将旬数换算为物候开始期对应的儒略日数；将最小变化率对应的下一旬数作为物候结束期对应的旬数，并将旬数换算为物候结束期对应的儒略日数；
13.步骤八、根据步骤五得到的逐像元的每年逐旬ndvi值及多年平均逐旬ndvi值，利用多元拟合最大值法构建一元六次多项式函数，拟合出逐像元的每年逐日ndvi拟合曲线及逐像元的多年平均逐日ndvi拟合曲线；
14.步骤九、将步骤七得到的物候开始期对应的儒略日数及结束期对应的儒略日数代入到逐像元的多年平均逐日ndvi拟合曲线中，得到物候开始期及结束期对应的ndvi阈值；
15.步骤十、分别将物候开始期及结束期对应的ndvi阈值带入到步骤八得出的逐像元的每年逐日ndvi拟合曲线中，得到ndvi阈值对应的天数，即研究时间段内每年的草本沼泽植被物候开始期及结束期。
16.步骤一所述数据预处理过程为：将逐旬归一化植被指数ndvi数据集以及两期沼泽湿地分布数据统一为相同的投影和坐标系。
17.步骤三所述未变化的草本沼泽分布是两期沼泽湿地分布数据中均为草本沼泽的像元分布。
18.步骤六所述多年平均逐旬ndvi值的变化率计算公式如下：
[0019][0020]
其中，t为旬数，ndvi
ratio
(t)是多年平均第t旬ndvi值的变化率，ndvi(t)是多年平均第t旬ndvi值，ndvi(t+1)是多年平均第t+1旬ndvi值。
[0021]
步骤七所述儒略日数是在一年中，从1月1日开始，到某一天结束所经历的日数。
[0022]
步骤八所述利用多元拟合最大值法构建一元六次多项式函数，拟合出逐像元的每年逐日ndvi拟合曲线及逐像元的多年平均逐日ndvi拟合曲线，公式如下：
[0023]
ndvi(x)＝a+a1x1+a2x2+a3x3+a4x4+a5x5+a6x6ꢀꢀꢀꢀ
公式(2)
[0024]
其中：x为儒略日数，ndvi(x)为第x儒略日对应的ndvi拟合值，a1，a2，a3…
a6为对应的系数。
[0025]
本发明的优点在于：
[0026]
与实地观测数据相比，遥感数据具有覆盖范围广，连续性强，数据易获取的优点。归一化植被指数(ndvi)能够在大的空间尺度上准确地反映植被的生长状态。本发明方法改善了传统物候观测范围小、周期长、观测物种有限，无法快速准确获取区域草本沼泽植被物候信息的缺点。本发明对于利用遥感植被数据及沼泽湿地分布数据，在区域尺度上提取草本沼泽植被物候信息提供了新的方法。本发明方法能利用遥感数据准确获取草本沼泽植被物候信息，保证了时效性和宏观性的同时，能够快速准确地提取整个区域草本沼泽植被物候信息。
附图说明
[0027]
图1是本发明的流程图；
[0028]
图2是实验一中2001-2015年三江平原未变化草本沼泽植被分布图；
[0029]
图3是实验一中2015年三江平原草本沼泽植被开始期空间分布图；
[0030]
图4是实验一中2015年三江平原草本沼泽植被结束期空间分布图。
具体实施方式
[0031]
本发明技术方案不局限于以下所列举具体实施方式，还包括各具体实施方式间的任意组合。
[0032]
具体实施方式一：本实施方式一种草本沼泽植被物候信息提取方法按照以下步骤进行：
[0033]
步骤一、获取研究时间段内覆盖研究区域的逐旬归一化植被指数ndvi数据集、以及两期沼泽湿地分布数据，并进行数据预处理；
[0034]
步骤二、利用最大值合成法将逐旬ndvi数据集合成为逐月ndvi数据集，并利用算术平均法计算得到多年平均生长季ndvi；
[0035]
步骤三、根据两期沼泽湿地分布数据，提取研究时间段内未变化的草本沼泽分布；
[0036]
步骤四、在未变化的草本沼泽分布范围内，提取多年平均生长季ndvi≥0.1的像元分布，得到未变化的草本沼泽植被分布，并将此作为研究区；
[0037]
步骤五、利用研究区裁切逐旬ndvi数据集，提取出逐像元的每年逐旬ndvi值，并利用算术平均法计算得到逐像元的多年平均逐旬ndvi值；
[0038]
步骤六、计算逐像元的多年平均逐旬ndvi值的变化率，并确定最大变化率对应的旬数和最小变化率对应的旬数；
[0039]
步骤七、将最大变化率对应的旬数作为物候开始期对应的旬数，并将旬数换算为物候开始期对应的儒略日数；将最小变化率对应的下一旬数作为物候结束期对应的旬数，并将旬数换算为物候结束期对应的儒略日数；
[0040]
步骤八、根据步骤五得到的逐像元的每年逐旬ndvi值及多年平均逐旬ndvi值，利用多元拟合最大值法构建一元六次多项式函数，拟合出逐像元的每年逐日ndvi拟合曲线及逐像元的多年平均逐日ndvi拟合曲线；
[0041]
步骤九、将步骤七得到的物候开始期对应的儒略日数及结束期对应的儒略日数代入到逐像元的多年平均逐日ndvi拟合曲线中，得到物候开始期及结束期对应的ndvi阈值；
[0042]
步骤十、分别将物候开始期及结束期对应的ndvi阈值带入到步骤八得出的逐像元的每年逐日ndvi拟合曲线中，得到ndvi阈值对应的天数，即研究时间段内每年的草本沼泽植被物候开始期及结束期。
[0043]
具体实施方式二：本实施方式与具体实施方式一不同的是步骤一所述数据预处理过程为：将逐旬归一化植被指数ndvi数据集以及两期沼泽湿地分布数据统一为相同的投影和坐标系。其他与具体实施方式一相同。
[0044]
具体实施方式三：本实施方式与具体实施方式一或二不同的是步骤三所述未变化的草本沼泽分布是两期沼泽湿地分布数据中均为草本沼泽的像元分布。其他与具体实施方式一或二相同。
[0045]
具体实施方式四：本实施方式与具体实施方式一至三之一不同的是步骤六所述多年平均逐旬ndvi值的变化率计算公式如下：
[0046][0047]
其中，t为旬数，ndvi
ratio
(t)是多年平均第t旬ndvi值的变化率，ndvi(t)是多年平均第t旬ndvi值，ndvi(t+1)是多年平均第t+1旬ndvi值。其他与具体实施方式一至三之一相
同。
[0048]
具体实施方式五：本实施方式与具体实施方式一至四之一不同的是步骤七所述儒略日数是在一年中，从1月1日开始，到某一天结束所经历的日数。其他与具体实施方式一至四之一相同。
[0049]
具体实施方式六：本实施方式与具体实施方式一至五之一不同的是步骤八所述利用多元拟合最大值法构建一元六次多项式函数，拟合出逐像元的每年逐日ndvi拟合曲线及逐像元的多年平均逐日ndvi拟合曲线，公式如下：
[0050]
ndvi(x)＝a+a1x1+a2x2+a3x3+a4x4+a5x5+a6x6ꢀꢀꢀꢀ
公式(2)
[0051]
其中：x为儒略日数，ndvi(x)为第x儒略日对应的ndvi拟合值，a1，a2，a3…
a6为对应的系数。其他与具体实施方式一至五之一相同。
[0052]
采用下述实验验证本发明效果：
[0053]
实验一：
[0054]
选择三江平原为实施区域，一种草本沼泽植被物候信息提取方法按照以下步骤进行：
[0055]
步骤一、获取研究时间段内覆盖研究区域的逐旬归一化植被指数ndvi数据集、以及两期沼泽湿地分布数据，并进行数据预处理；
[0056]
获取覆盖三江平原区域的2001-2015年modis ndvi数据集、2001和2015年两期沼泽湿地分布数据，将逐旬归一化植被指数ndvi数据集以及两期沼泽湿地分布数据统一为相同的投影和坐标系。
[0057]
步骤二、利用最大值合成法将逐旬ndvi数据集合成为逐月ndvi数据集，并利用算术平均法计算得到2001-2015年多年平均生长季(5-9月)ndvi；
[0058]
步骤三、根据两期沼泽湿地分布数据，提取研究时间段内未变化的草本沼泽分布；
[0059]
提取2001和2015年三江平原沼泽湿地分布数据中均为草本沼泽的像元，即为2001-2015年三江平原未变化的草本沼泽分布。
[0060]
步骤四、在未变化的草本沼泽分布范围内，提取多年平均生长季ndvi≥0.1的像元分布，得到未变化的草本沼泽植被分布，并将此作为研究区；
[0061]
图2为2001-2015年三江平原未变化的草本沼泽植被分布图。
[0062]
步骤五、利用研究区裁切逐旬ndvi数据集，提取出逐像元的每年逐旬ndvi值，并利用算术平均法计算得到逐像元的多年平均逐旬ndvi值；
[0063]
步骤六、计算逐像元的多年平均逐旬ndvi值的变化率，并确定最大变化率对应的旬数和最小变化率对应的旬数；
[0064]
步骤七、将最大变化率对应的旬数作为物候开始期对应的旬数，并将旬数换算为物候开始期对应的儒略日数；将最小变化率对应的下一旬数作为物候结束期对应的旬数，并将旬数换算为物候结束期对应的儒略日数；
[0065]
步骤八、根据步骤五得到的逐像元的每年逐旬ndvi值及多年平均逐旬ndvi值，利用多元拟合最大值法构建一元六次多项式函数，拟合出逐像元的每年逐日ndvi拟合曲线及逐像元的多年平均逐日ndvi拟合曲线；
[0066]
步骤九、将步骤七得到的物候开始期对应的儒略日数及结束期对应的儒略日数代入到逐像元的多年平均逐日ndvi拟合曲线中，得到物候开始期及结束期对应的ndvi阈值；
[0067]
步骤十、分别将物候开始期及结束期对应的ndvi阈值带入到步骤八得出的逐像元的每年逐日ndvi拟合曲线中，得到ndvi阈值对应的天数，即2001-2015年每年的草本沼泽植被物候开始期及结束期。
[0068]
以2015年为例，分别将物候开始期及结束期对应的ndvi阈值带入到逐像元的2015年逐日ndvi拟合曲线中，得到2015年三江平原草本沼泽植被开始期(图3)及结束期(图4)的空间分布。
[0069]
步骤三所述未变化的草本沼泽分布是两期沼泽湿地分布数据中均为草本沼泽的像元分布。
[0070]
步骤六所述多年平均逐旬ndvi值的变化率计算公式如下：
[0071][0072]
其中，t为旬数，ndvi
ratio
(t)是多年平均第t旬ndvi值的变化率，ndvi(t)是多年平均第t旬ndvi值，ndvi(t+1)是多年平均第t+1旬ndvi值。
[0073]
步骤七所述儒略日数是在一年中，从1月1日开始，到某一天结束所经历的日数。
[0074]
步骤八所述利用多元拟合最大值法构建一元六次多项式函数，拟合出逐像元的每年逐日ndvi拟合曲线及逐像元的多年平均逐日ndvi拟合曲线，公式如下：
[0075]
ndvi(x)＝a+a1x1+a2x2+a3x3+a4x4+a5x5+a6x6ꢀꢀꢀꢀ
公式(2)
[0076]
其中：x为儒略日数，ndvi(x)为第x儒略日对应的ndvi拟合值，a1，a2，a3…
a6为对应的系数。