[发明专利]一种制备逐年土地利用数据的方法

权利要求书

1.一种制备逐年土地利用数据的方法，其特征在于：步骤如下：

1)数据获取

发明采用的人工解译土地利用数据A-LUCC，数据格式为矢量，时间为不连续的年份，时间间隔为5年，分辨率较精细，为30m；采用的逐年土地利用数据M-LUCC，时间间隔为1年，分辨率较粗，为500m；

2)数据预处理

首先对多期A-LUCC数据进行几何精度校正，几何精度校正是通过图像上若干已知点，即地面控制点，构造一个图像坐标与地理坐标之间的关系式，假定图像坐标为(X，Y)，地理坐标为(U,V)，则地理坐标与图像坐标之间存在下列关系：

U＝a+bX+cY+dX²+eXY+fY²+…+pX^αY^β+…+zY^γ

V＝a′+b′X+c′Y+d′X²+e′XY+f′Y²+…+p′X^αY^β+…+z′Y^γ

为保证精度和效率，采用9个地面控制点的坐标进行拟合，拟合方程为：

U＝a+bX+cY+dX²+eXY+fY²+gX³+hX²Y+kXY²+lY³

V＝a′+b′X+c′Y+d′X²+e′XY+f′Y²+g′X³+h′X²Y+k′XY²+l′Y³

然后，对多期A-LUCC数据进行矢量转栅格，并基于最近邻算法进行栅格数据重采样，即将空间分辨率重采样为500m，以保证和M-LUCC的空间分辨率一致；然后对A-LUCC和M-LUCC的土地利用类型，采用大类的分类标准，分为农田、林地、草地、湿地、建设用地和其他未利用地六大类；

3)预设分析窗口

本发明通过均匀网格化(Createfishnet)功能来生成分析格网，并以此判定最小的分析窗口；即，把整个研究区分成q×q，(q+1)×(q+1),…规则亚区作为分析窗口，在亚区内逐一判读2个研究时段内均有的土地利用转类，最小的网格亚区即为满足条件的最优窗口；其中q为≥1的自然数；

4)时间变化速率判定

对不同土地利用类型的时间变化速率，以A-LUCC的时间段各土地利用类型的变化值为限定值，参考M-LUCC年际变化值，采用对应数据差值的方法来计算，此处称之为逐年类归一化线性拟合处理，公式表示如下：

A-LUCC数据中某土地利用转类i的j+n年中j＝1,2,…，且n≥1；变化总量ATR_{ni_sum}用如下公式表示：

ATR_{ni_sum}＝ATR_j+n-ATR_j

M-LUCC数据中某土地利用转类i的逐年变化量之和MTR_{ni_sum}等于ATR_{ni_sum}或不等于ATR_{ni_sum}，MTR_{ni_sum}不等于ATR_{ni_sum}表达式为：

MTR_{ni_sum}＝∑(MTR_j+n-MTR_j+n-1，…，MTR_j+1-MTR_j)

式中，MTR为M-LUCC对应转类i的对应年份j+n的面积；(MTR_j+n-MTR_j+n-1)、(MTR_j+1-MTR_j)、…为逐年际间各土地利用转类的面积变化差值；则，如果MTR_{ni_sum}≠ATR_{ni_sum}时，需要采用如下方法处理：

使M-LUCC数据各土地利用类型变化的量等于两期A-LUCC数据各土地利用类型变化的量，即：

式中，j+m，j+n和j分别为研究时间段的任意年份、末年和首年，其中1<m<n；则有此式推导出逐年LUCC类归一化的土地利用类型i对应的面积值为：

式中，TR_j+m为LUCC数据中某土地利用转类i的面积；ATR和MTR分别为A-LUCC和M-LUCC中某一转类i的面积；

5)空间变化的定位

以A-LUCC数据空间分布为准，把对应的前后两期数据变化位置圈定不变；有两种情况：一是如果斑块的动态分布范围与在A-LUCC连续两期研究时段内t和t+1的同类斑块(Grid)动态分布范围一致，则为正确的动态转类空间定位，即任意栅格位置对应满足：Grid_t-Grid_t+1＝0；二是如果空间范围不一致时，则以M-LUCC的同转类栅格单元分布位置作为参考判定空间转换位置，通过观察各个变化单元的中心点所处位置，先填充相同方向上的A-LUCC变化空间，如果该空间填充不完，再向周边扩展；

6)精度检验

基于以上对时间变化速率判定和空间变化定位判定的处理，融合A-LUCC数据和M-LUCC数据，获取得到研究区2000-2010年的500m空间分辨率的逐年土地利用数据；逐年土地利用数据用研究时段内覆盖部分研究区的不连续遥感数据进行验证；或选用时间变化速率逐年累加得到的值作为标准的理论值，空间变化定位逐年视域判读，进行精度自检验。

2.根据权利要求1所述的一种制备逐年土地利用数据的方法，其特征在于步骤1)所述的数据源均为土地利用现状图层。