一种基于多源数据的逐年土地利用空间数据的制备方法
【技术领域】
[0001] 本发明涉及土地利用变化逐年数据的制备。主要运用非同源土地利用数据的时 间和空间值的精度互补的思想，以A-LUCC空间分布为准，以同分析窗口下M-LUCC的时间变 化为准，通过非同源数据的有效结合实现逐年的土地利用变化数据的制备，并且可以在一 定的分析窗口和尺度单元下保证精度，从而为土地资源规划管理和发展评估提供科学的参 考资料和决策依据。属于地球科学领域。
【背景技术】
[0002] 土地利用/ 土地覆盖变化（Land Use / Cover Change,简称LUCC)时空格局变化 可延伸为时间变化和空间变化这两个方面。时间变化研究的主线为：不同时期土地利用变 化类型的面积动态指数和转移矩阵，以其为评价指标来量化和分级土地利用的变化速度、 变化幅度和转移比例等。空间变化的研究主线为：以地理信息科学技术为研究工具，依据土 地利用动态信息，划分区域土地利用类型的空间分布格局，进而明确土地利用类型随时间 演化的动态转化过程。
[0003] 对LUCC空间和时间过程变化的研究主要集中在以下两种类型：第一种是模型模 拟研究，De Almeida等将元胞自动机模型应用到巴西城市Bauru的土地利用变化和扩展模 拟。黎夏、邓祥征等人分别提出的城市CA模型和区域土地利用类型变化模型；第二种是年 际变化研究，如刘纪远等以1980年代末、1995年、2000年和2005为时间截面，以5年为一 研究阶段，揭示了 20世纪80年代末以来中国土地利用的时空变化过程。
[0004] 基于时间序列的土地利用变化宏观研究已经日臻成熟，受限于数据，对LUCC"年际 变化"的长时间序列的过程研究却少之又少。而从时间尺度和空间尺度同时入手，以统计学 理论和空间分析技术为工具，探索土地利用时空演化规律，克服土地利用变化研究中，高精 度的土地利用数据无法逐年提供的问题，并在一定尺度范围内保证精度，已逐步成为土地 利用变化研究的前沿思维。
【发明内容】
[0005] 本发明旨在提出一种逐年土地利用空间数据的制备方法，以便可以实现长时间序 列的逐年土地利用变化分析。该方法采用不同源的两种土地利用数据，结合其各自提供的 有效信息，通过数据融合方式，实现逐年土地利用空间数据的制备。该方法可以在一定程度 上替代模型模拟，是一种新的、简便易行的、可快速用于实践研究的方法。
[0006] 为实现上述目的，本发明提供一种逐年土地利用空间数据的制备方法，其步骤依 次为： 1)数据获取： 两种不同源数据：一是利用较高分辨率对地观测卫星影像30m空间分辨率的Landsat TM/ETM+，通过人工解译制备的土地利用变化数据集（A-LUCC);二是500 m空间分辨率 MODIS的土地覆盖数据（M-LUCC)。
[0007] 2)数据预处理：首先对多期A-LUCC数据进行地理空间分析算法中的矢量转栅格， 并基于最近邻算法进行栅格数据重采样，即将空间分辨率重采样为500 m，以保证和M-LUCC 的空间分辨率一致。然后对A-LUCC和M-LUCC的土地利用类型，采用大类的分类标准，分 为农田、林地、草地、湿地(含水域)、建设用地和其他未利用地等六大类(简称农、林、草、湿、 城、荒)。
[0008] 3)预设分析窗口 本发明通过均匀网格化来判定最小的分析窗口，即把整个研究区分成1X1，4X4，… 等亚区作为分析窗口，在亚区内逐一判读2个研究时段内均有的土地利用转类，而满足此 条件的最小网格亚区即为所需的最优分析窗口，分析窗口的上限是整个研究区。在取值区 间内分析窗口越小精度越高。
[0009] 4)时间变化的判定 对不同土地利用类型的时间变化速率，以A-LUCC的时间段变化值为限定值，参考 M-LUCC的年际变化值来计算，这里称之为逐年类归一化线性拟合处理，具体推导过程如 下： A-LUCC数据中某土地利用转类i的J 年（J = 7,之…,· /?夕7)变化总量 _可用如下公式表示：
A-LUCC数据中某土地利用转类i的逐年变化量之和#7?, _并非一定等于，此 值的表达式为：
式中，#7?为M-LUCC对应转类i的对应年份J 的变化速率；（#7仏" - (#7^, - #%)、···为逐年际间各土地利用转类的面积变化差值。则，如果#7?, _时，需要米用如下方法处理： 使M-LUCC数据各土地利用类型变化的量等于两期A-LUCC数据各土地利用类型变化的 量，即：
式中，J - λ和J分别为研究时间段的任意年份、末年和首年对应的土地利用类 型i变化值，其中7 则有此式可以推导出逐年LUCC类归一化的土地利用类型Y 对应的面积值为：
式中，Mj7J^LUCC数据中某土地利用转类i的面积；J7^/7 #7?分别为A-LUCC和M-LUCC 中某一转类i的面积。
[0010] 5)空间变化的判定 以A-LUCC空间分布为准，同分析窗口下M-LUCC逐年的同转类斑块分布位置作为参考， 在变化速率值的限定下，利用条件语句判定空间转换位置。即： 如果该斑块的动态分布范围与在A-LUCC连续两期研究时段内（?和?-7)的同类斑 块（Grid)动态分布范围一致，则为正确的动态转类空间定位，即任意栅格位置对应满足： = <9;否则，空间范围不一致时，则定位该转类斑块到同窗口下A-LUCC的最 邻近同转类斑块的所属空间，超出A-LUCC对应转类斑块空间的依照原M-LUCC的分布空间， 进行同分析窗口下的质点所在空间趋势位置分为八个方向（东、南、西、北、东北、西北、东 南、西南），再随机分布定位：即，当质点偏向某一方向，就在某方向所在位置随机布置空间 转类栅格类型。
[0011] 其中，质心是描述地理目标空间分布最有用的单一量算量，它是保持目标均匀分 布的平衡点。这里引入质心是旨在通过几何中心的基础上的加权运算来追踪土地类型的变 化。任意多边形的质点求算可以通过公式表示： 任意边λ边形顶点為jv*，G = Λ之…，W依照逆时针排列，其面积为：
贝ij，对应任意边λ边形的质点可表达为：
6)精度验证和长时间序列时空分析 利用研究时段内的A-LUCC数据，随机选取研究区中某一年的数据，对应检验各类型的 空间分布和面积是否一致，进而给出精度。或者更进一步，随机选取某年或某区域内人机交 互解译的A-LUCC数据，以此对比检验同一年份的新制备数据，给出对应精度。
[0012] 土地利用长时间序列时空分析是指根据制备出来的数据结果，看其空间变化情 况，主要统计该研究时间段的各土地利用类型的数量。
[0013] 本发明方法的优势和特点： 1)可间接获取连续时间序列的土地利用数据，对不连续的遥感数据进行精度分析和质 量评价，从而提高土地利用变化模拟的效果。
[0014] 2)将规则网格方法运用到土地利用转类的判读中，使分析对象更加地精细化，也 使逐年信息的扩展更加地合理和准确。
[0015] 3)通过中等空间分辨率和较高分辨率的遥感数据之间建立定量关系，互取其优 点，为其他多源数据的融合处理提供了技术参照。
[0016]
【附图说明】
[0017] 图1是制备数据的分析窗口判读示意图 图2是制备逐年土地利用空间数据的方法流程图
【具体实施方式】
[0018] 下面结合
【发明内容】
和附图1和2对本发明的原理及【具体实施方式】作进一步阐述。 包括如下步骤： 1.数据获取 利用较高分辨率对地观测卫星影像30 m空间分辨率的Landsat TM/ETM+，通过人工 解译制备的土地利用变化数据集（A-LUCC);利用500 m空间分辨率MODIS的土地覆盖数据 (M-LUCC)0
[0019] 2.数据预处理 (I)A-LUCC人机交互解译。首先对多期A-LUCC数据进行地理空间分析算法中的矢量 转栅格，并基于最近邻算法进行栅格数据重采样，即将空间分辨率重采样为500 m，以保证 和M-LUCC的空间分辨率一致。利用屏幕数字化方法以研究区内对相同的视图比例人机交 互解译放大的较高空间分辨率影像，同时开展野外考察，确定分类标准，并更正分类结果， 以满足本研究的需要。
[0020] (2)对M-LUCC 土地覆盖数据对应A-LUCC的解译分类标准，最终统一将两套数据的 土地利用类型分为农、林、草、湿、城、荒六类。A-LUCC的解译结果即为农、林、草、湿、城、荒六 类，因此不涉及栅格属性的调整；而对M-LUCC，其土地覆盖类型(表1)需要，根据表2统一 属性为：农、林、草、湿、城、荒六类。
[0021] 表1五种土地覆盖分类方案 CN 105184048 A 说明书 5/8 页
注：标准I: IGBP，国际地圈生物圈计划-全球植被分类方案；标准2: UMD，马里兰大 学方案；标准3: LAI/fPAR，基于叶面积指数/光合有效辐射吸收比例的MODIS方案；标准 4: NPP，基于净初级生产量的MODIS方案；标准5: PFT，植物功能类型方案。
[0022] 表2对LUCC和MODIS进行重分类规则表
3.规则网格 网格尺度即分析窗口大小的选择非常重要，它直接影响了分析窗口的大小和信息提取 的精度。最佳的网格尺度是既能避免对象多边形太破碎，又能保证A-LUCC和M-LUCC在多 个时间段均有相同的土地利用转换类型（图1)。本发明通过每次设定不同的网格宽度和高 度值，如（1X1)，（4X4)，…等来把整个研究区划分成iXi，i+IXi+l，…等规则亚区， 如果第n+1个亚区内两