[发明专利]基于像元尺度的盐碱地盐碱化程度时空变化的表达方法

说明书

基于像元尺度的盐碱地盐碱化程度时空变化的表达方法，涉及一种基于像元尺度的盐碱地盐碱化程度表达方法。本发明是为了解决现有方法难以准确地模拟盐碱地内部变化状况的问题，方法：一、构建盐碱地像元尺度上植被状况变化趋势模型；二、构建盐碱地像元尺度上植被状况变化波动模型；三、对构建模型的精度验证。本发明通过构建试验区盐碱地像元尺度上的植被状况NDVI的趋势和波动模型，从空间和时间尺度上提高盐碱地图斑内植被状况的变化，从植被的时空精细化变化反应盐碱化程度，充分体现出利用卫星遥感数据表达盐碱化程度的科学性。本发明属于盐碱地植被状况时序变化表达领域。

技术领域

本发明涉及一种基于像元尺度的盐碱地盐碱化程度表达方法。

背景技术

植被状况是盐碱地盐碱化程度的一个重要表征指标。关于区域尺度上盐碱化程度的时间序列变化特征研究往往采用对区域内盐碱地植被状况变化构建时间序列模型，进而表征该区域内盐碱地变化的总体情况。由于盐碱化程度受盐碱土壤、地形、水文等条件的综合作用而存在显著的空间差异性，区域内盐碱地植被状况变化趋势呈现显著不同。现有模型主要有二分模型、线性光谱混合模型等，但缺少对盐碱地植被状况变化趋势的相关像元模型，也就是说现有方法难以准确地模拟盐碱地内部变化状况。

发明内容

本发明的目的是为了解决现有方法难以准确地模拟盐碱地内部变化状况的技术问题，提供了一种基于像元尺度的盐碱地盐碱化程度时空变化的表达方法。

基于像元尺度的盐碱地盐碱化程度时空变化的表达方法如下：

一、构建盐碱地像元尺度上植被状况变化趋势模型

在获取研究区不同年份国家土地利用数据库中相同月份时间断面盐碱地类别属性未发生改变的盐碱地斑块基础上，进一步获得盐碱地图斑范围内表征盐碱地植被状况的NDVI指数和气象数据，以表征盐碱地植被状况的NDVI为自变量，分别对研究区内的盐碱地图斑建立像元尺度上的趋势模型，并以方差和均方误差确立最优趋势模型；

二、构建盐碱地像元尺度上植被状况变化波动模型

对比研究区各盐碱地像元时间序列趋势模型的NDVI模拟数值与实际数值，对每个像元的差值部分构建有关降水量和蒸发量的二元线性回归波动模型，得到研究区盐碱地像元植被状况时空变化精细化模型；

三、对构建模型的精度验证

确定已构建的研究区盐碱地像元植被状况时空变化精细化模型的精度，以评判建立模型的优劣和实用性。

步骤一获取盐碱地类别属性未发生改变的盐碱地斑块，其数据精度为像元大小30m。

步骤一所述相同月份时间为每年5-9月份。

步骤一所述获得盐碱地图斑范围内表征盐碱地植被状况的NDVI指数和气象数据，其数据精度为像元大小1km。

步骤一所述趋势模型包括线性函数趋势模型、幂函数趋势模型、对数函数趋势模型、二次多项式函数趋势模型。

本发明在像元尺度上构建东北地区松嫩平原盐碱地植被状况时序变化趋势模型，模拟每一盐碱地像元植被状况变化趋势。本发明通过对粗空间分辨率长时间序列精细时间分辨率的植被覆盖数据建立趋势和波动参数，来修正高空间分辨率、粗时间分辨率的盐碱地数据，达到提高对盐碱地盐碱化程度时空精细化表达的目的。

本发明通过构建试验区盐碱地像元尺度上的植被状况NDVI的趋势和波动模型，从空间和时间尺度上提高盐碱地图斑内植被状况的变化，从植被的时空精细化变化反应盐碱化程度，充分体现出利用卫星遥感数据表达盐碱化程度的科学性。

附图说明

图1是实验一中研究区盐碱地像元NDVI线性函数趋势模型a为0.0049的参数分布图；