CN107818519A - 一种分时序省级土地利用数据处理方法及系统 - Google Patents
一种分时序省级土地利用数据处理方法及系统 Download PDFInfo
- Publication number
- CN107818519A CN107818519A CN201711036172.4A CN201711036172A CN107818519A CN 107818519 A CN107818519 A CN 107818519A CN 201711036172 A CN201711036172 A CN 201711036172A CN 107818519 A CN107818519 A CN 107818519A
- Authority
- CN
- China
- Prior art keywords
- land use
- msub
- mrow
- provincial
- land
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/284—Relational databases
- G06F16/285—Clustering or classification
Abstract
本发明公开了一种分时序省级土地利用数据处理方法及系统，通过对土地利用数据的处理，构建土地利用变化影像，可快速清晰的获取土地利用变化的空间特征，进行土地利用变化的定性分析；通过计算土地利用结构、单一土地利用动态度和空间信息分析模型数据三个参数，对省域级别的土地利用数据进行精确的计算，对土地利用变化研究提供数据支撑，具有良好的实用性。
Description
技术领域
本发明涉及土地利用研究领域，具体涉及一种分时序省级土地利用数据处理方法及系统。
背景技术
自“土地利用/覆被变化研究计划(LUCC)”和“全球土地计划(GLP)”启动后，学者们对于土地利用/土地覆盖变化(LUCC)的研究主要围绕土地利用变化时空特征及其驱动力机制研究、土地利用变化预测模拟和土地利用变化对环境的产生的效应等各方面。
目前在土地利用研究领域，基础数据大部分来自于遥感地图影像，如何对遥感地图影像数据进行加工处理，使其能快速应用到土地利用研究中，是该领域内的一个研究方向。
发明内容
为了弥补现有土地利用研究领域的空白，本发明提供一种分时序省级土地利用数据处理方法及系统，在多个连续的子时序上，在地图上对土地利用发生变化的类型进行统计并在地图上以不同颜色标出，可较为清晰的观察土地利用变化的空间分布特征，从而进行土地利用变化的定性分析；通过土地利用结构、单一土地利用动态度和空间信息分析模型三个参数，可对土地利用变化进行定量分析，为全国省域级别的土地利用研究提供数据支持。
相应的，本发明提供了一种分时序省级土地利用数据处理方法，包括以下步骤：
确定分时序省级分析目标和分析时间段；
在所述分析时间段中添加分隔节点时间，所述分隔节点时间将所述分析时间段划分为一个以上的子时序；
从全国土地利用数据库中提取所述分析时间段起始时间、结束时间和所述分隔节点时间的所述分时序省级分析目标的土地利用数据；
构建所述土地利用数据的土地利用类型分类标准；
在所述分时序省级分析目标的地图上，对在子时序结束时间相对于起始时间发生土地利用类型变化的区域进行标示。
按所述土地利用类型分类标准计算在所述分时序省级分析目标中，每个土地利用类型在所述分析时间段起始时间、结束时间和所述分隔节点时间的土地利用结构、单一土地利用动态度和空间信息分析模型；
将所述土地利用结构、单一土地利用动态度、空间信息分析模型信息汇总。
所述土地利用数据为通过分辨率为30米的LandsatTM数字影像解译获取的分辨率为1km的土地利用分类栅格数据。
所述土地利用类型分类标准包括6个一级类型和25个二级类型。
所述方法还包括，将所述分时序省级分析目标的土地利用数据中的二级类型数据向所述一级类型进行转换。
在所述分时序省级分析目标的地图上，所述土地利用变化类型相同的区域用同一颜色进行标示。
所述土地利用结构用于表示第i种土地利用类型面积所占的比重，公式为Ai表示第i种土地利用类型的面积；n表示土地利用类型的类别数量。
所述土地利用度用于描述第i种土地利用类型在[t1,t2]时间段内的变化速度，用单一土地利用动态度定量描述，公式为其中，LA(i,t1)和LA(i,t2)为研究期起始和结束时第i种土地利用类型的面积；t1,t2研究期起始时间和结束时间。
所述空间信息分析模型包括转移速率TRLi、新增速率IRLi和变化速率CCLi，计算公式为
其中，转移速率TRLi用于表示第i种土地利用类型在t1至t2期间转移面积与期初面积比重，IRLi用于表示第i种土地利用类型在t1至t2期间新增面积与期初面积比重的新增速率，CCLi第i种土地利用类型在t1至t2期间转移面积和新增面积之和与期初面积比重的变化速率；LA(i,t1)表示t1时间第i种土地利用类型的面积，LA(i,t2)表示t2时间第i种土地利用类型的面积，ULAi表示第i种土地利用类型在时序开始与结束没发生改变的面积。
相应的，有一种分时序省级土地利用变化数据处理系统，包括
分析目标模块：用于选择分析省份及分析时间段；
子时序划分模块：用于在所述分析时间段中添加分隔节点时间，将所述分析时间段划分为一个以上的子时序；
土地利用数据提取模块：用于从全国土地利用数据库中提取所述分时序省级分析目标的土地利用数据；
土地利用变化计算模块：用于计算所述土地利用类型在所述子时序结束时间相对于起始时间的土地利用变化类型，并对不同的土地利用变化类型分类；
土地利用变化影像构建模块：用于在分时序省级分析目标的地图，用不同颜色对不同类型的土地利用变化进行标示；
土地利用类型分类标准模块：用于根据预设的土地利用类型分类标准，将所述分时序省级分析目标的土地利用数据分类；
土地利用结构计算模块：用于计算土地利用结构；
单一土地利用动态度计算模块：用于计算单一土地利用动态度；
空间信息分析模型计算模块：用于计算空间信息分析模型；
制表储存模块：用于将所述土地利用数据，土地利用结构、单一土地利用动态度、空间信息分析模型信息汇总并制表储存。
所述系统还包括：
转换模块：用于将所述分时序省级目标的土地利用数据中的二级类型的土地利用转换为一级类型的土地利用。
本发明提供一种分时序省级土地利用数据处理方法及系统，在多个连续的子时序上，在地图上对土地利用发生变化的类型进行统计并在地图上以不同颜色标出，可较为清晰的观察土地利用变化的空间分布特征，从而进行土地利用变化的定性分析；通过土地利用结构、单一土地利用动态度和空间信息分析模型三个参数，可对土地利用变化进行定量分析。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1为本发明实施例的分时序省级土地利用数据处理方法流程图；
图2为土地利用类型分类标准示意图。
图3为1990年－2000年土地利用变化示意图。
图4为2000年－2010年土地利用变化示意图。
图5为2010年－2015年土地利用变化示意图。
图6为1990年－2000年土地利用结构、土地利用度、空间信息分析模型数据表。
图7为2000年－2010年土地利用结构、土地利用度、空间信息分析模型数据表。
图8为2010年－2015年土地利用结构、土地利用度、空间信息分析模型数据表。
图9为本发明实施例的分时序省级土地利用数据处理系统框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
图1示出了本发明实施例的分时序省级土地利用数据处理方法的流程图。本发明实施例提供一种分时序省级土地利用数据处理方法及系统，包括以下步骤：
S101：确定分时序省级分析目标和分析时间段，本发明实施例选择的分析目标为广东省，分析时间段为1990年－2015年。
S102：划分子时序；在所述时间段中添加分隔节点时间，该分隔节点时间将分析时间段划分为一个以上的子时序；本发明实施例的分析时间段起始时间为1990年，结束时间为2015年，选取的时序分隔节点时间为2000年和2010年。具体实施中，书序分隔节点时间的选取避免过长或过短，以免时序过长造成丢失部分土地利用变化特征或时序过短造成土地利用变化不明显，影响分析结果；亦可采用将时间段均分成时序的方法，使划分时序更为合理；
S103：提取分时序省级分析目标的土地利用数据；本发明实施例从全国土地利用数据库中提取在所述分析时间段起始时间、结束时间以及所述时序分隔节点时间上的所述分析目标的土地利用数据。本发明实施例的土地利用数据来源于中国科学院资源环境数据中心全国土地利用数据库，该数据通过采用土地利用变化遥感信息人机交互快速提取方法，解译1990年、2000年、2010年、2015年的分辨率为30米的Landsat TM数字影像，截取其中的广东省土地利用数据，进行后续分析研究；具体实施中，Landsat TM数字影像的数据包中，包括有单位距离为1km的坐标数据以及坐标数据所对应的遥感图像信息，通过计算机软件，如Arcmap，可生成该分时序省级分析目标的遥感信息图；
S104：构建土地利用类型分类标准；具体实施中，在我国土地资源遥感调查与监测技术规程中，刘纪远等提出了基于遥感TM影像的土地利用/覆盖分类系统设计原则，并在此基础上设计了完整的土地利用覆盖分类系统；本文采用该土地分类体系，将研究区土地利用类型划分为耕地、林地、草地、水域、城乡工矿居民用地和未利用地6个一级类型及25个二级类型，其中一级类土地利用类型解译精度达到94.3％以上，具有较理想的精度；
S105：类型转换；本发明实施例对土地利用一级类型进行研究，因此，需要将二级类型取整转换为一级类型，转换规则参照图2所示的土地利用类型分类标准表，将二级类型转换为其所属的一级类型；具体实施中，可根据数据分析的目的决定分类标准，以满足研究需求；
S106：构建土地利用变化影像；在所述分时序省级分析目标的地图上，对在子时序结束时间相对于起始时间发生土地利用类型变化的区域进行标示。图3～图5示出了1990、2000、2010、2015年土地利用变化示意图；本发明实施例提供的分时序省级土地利数据处理方法，主要用于研究社会经济数据、宏观政策、地方政策三个驱动力因子对土地利用变化的影响，而土地利用变化大部分是从非人工用地向不同的非人工用地或人工用地转化，因此，图3～图5主要标出了耕地、林地、草地、水域向其余土地利用类型的转化区域；具体实施中，可按照研究的目标选择标示的土地利用变化类型；比如研究耕地复耕情况，应着重于非耕地土地利用类型向耕地土地利用类型的装换；
S107：计算土地利用结构；所述土地利用结构用于表示第i种土地利用类型面积所占的比重，公式为Ai表示第i种土地利用类型的面积；n表示土地利用类型类别数；耕地、林地、草地、水域、城乡工矿居民用地和未利用地依次分别为P1、P2、P3、P4、P5、P6；
S108：计算土地利用度；所述土地利用度用于描述第i种土地利用类型在[t1,t2]时间段内的变化速度，用单一土地利用动态度定量描述，公式为其中，LA(i,t1)和LA(i,t2)为研究期起始和结束时第i种土地利用类型的面积；t1,t2研究期起始时间和结束时间；耕地、林地、草地、水域、城乡工矿居民用地和未利用地依次分别为K1、K2、K3、K4、K5、K6；
S109：计算所述空间信息分析模型；所述空间信息分析模型包括转移速率TRLi、新增速率IRLi和变化速率CCLi，计算公式为
其中，转移速率TRLi用于表示第i种土地利用类型在t1至t2期间转移面积与期初面积比重，IRLi用于表示第i种土地利用类型在t1至t2期间新增面积与期初面积比重的新增速率，CCLi第i种土地利用类型在t1至t2期间转移面积和新增面积之和与期初面积比重的变化速率
LA(i,t1)表示t1时间第i种土地利用类型的面积，LA(i,t2)表示t2时间第i种土地利用类型的面积，ULAi表示第i种土地利用类型在时序开始与结束没发生改变的面积，IRLi用于表示第i种土地利用类型在t1至t2期间新增面积与期初面积比重的新增速率；
S110：信息汇总；本发明实施例选用土地利用结构、单一土地利用动态度、空间信息分析模型三项指标来对广东省1990－2015年过程中各时序土地利用变化的空间格局进行定量的对比、分析和描述，具体实施中，可选用多种指标综合分析；图6～图8汇总了本发明实施例的1990－2000年、2000－2010年、2010－2015年土地利用类型面积变化及土地利用结构、单一土地利用动态度、空间信息分析模型计算值数据图。
以上为本发明实施例提供的分时序省级土地利用数据处理方法的实施步骤，通过该数据处理方法，可以对分时序省级土地利用数据进行处理。
本发明实施例通过对1990－2000年、2000－2010年、2010－2015年三个子时序的土地利用数据进行分析，将广东省划分为珠三角平原区、粤东沿海区、粤西沿海区、粤西北山区四个区域，通过对图3～图5的观察，其土地利用类型变化方面的空间特征如下：
1990－2000年，全省土地利用变化规模较小，主要集中于珠三角平原区，主要表现为耕地数量的减少和水域、城乡工矿居住用地的增加。2000－2010年，全省土地利用变化更大，珠三角地区、粤东沿海区尤为剧烈，总体趋势为各类用地向城乡工矿居住用地的转变，其中以耕地及林地向城乡工矿居住用地的转变最为突出。2010－2015年，全省土地利用向城乡工矿居住用地转变的趋势持续，但在空间上趋于分散，各区主要的土地变化以耕地、林地向城乡工矿居住用地转变及林地向草地转变为主。
通过图6～图8分析每种土地利用类型的土地利用结构、土地利用度和所述空间信息分析模型并结合图3～图5的时空特征，可分析出：
就全省而言，土地利用类型总体的变化趋势表现为耕地、林地的持续减少及城乡工矿居住用地的持续增加，但土地利用类型变化的速度有所放缓。珠三角地区为全省土地利用类型变化规模最大的区域，以耕地、水域转换为城乡工矿居住用地为主；粤东、粤西北两区的土地利用变化相对珠三角区域较小，但通过对三个时序的比较可看出，区域内草地、城乡工矿居住用地的规模有不同程度的增加，此外，粤西北的东侧在2010－2015年间出现了林地向耕地的变化；粤西地区土地利用变化不大，以耕地的减少，林地、城乡工矿居住用地的增加为主。
就用地类型而言，1990年耕地和林地作为最主要的用地类型其面积合计占整个广东省土地面积85％以上，随时间推移，耕地面积减少明显，至三个时间节点分别减少18.56万公顷、24.71万公顷、4.35万公顷，而单一土地利用动态度Ki则经历了由加快至转慢的过程，占全省用地的比重亦从1990年的27.53％降低至2015年的24.87％；
林地面积逐步缩减，至三个时间节点分别减少0.37万公顷、6.14万公顷、6.2万公顷，而单一土地利用动态度Ki逐步加快，占全省用地的比重变化则从1990年的60.77％降低至2015年的60.41％；
草地的面积和比重经历了先减后增的过程，从1990年的面积83.41万公顷减少至2010年的81.43万公顷再增加至2015年的83.31万公顷，单一土地利用动态度Ki则逐步加快。
城乡工矿居住用地的变化在各用地类型中最为明显，面积由1990年的46.43万公顷增加至2015年的103.77万公顷，比重则由2.59％快速上升至5.78％，单一土地利用动态度Ki保持较高速度；水域面积总体上虽从1990年的78.31万公顷增加至81.74万公顷，比重亦从4.36％上升至4.55％，但实际上则经历了先快速增加后逐步减少的过程，单一土地利用动态度Ki也由1990年的1.11％逐步加快减少为2015年的﹣0.27％。
未利用地的面积、比重变化均不大，整体呈逐年减少的趋势。
其中，针对分析时间段的子时序进行研究分析，可得
耕地土地应用类型变化基本特征：
耕地总量持续减少。耕地总量在1990年、2000年、2010年、2015年分别为494.29万公顷、475.73万公顷、451.02万公顷、446.67万公顷，而2000年、2010年、2015年新增耕地面积仅分别为2万公顷、0.21万公顷、0.38万公顷，增加的数量远远低于减少量。
耕地减少的速度由快减慢。通过变化速率CCLi可更为直观地看出耕地减少的速度经历了先加速后减速的过程：1990－2010年为加速阶段，耕地的变化速率从1990－2000年的0.46％扩大到2000－2010年的0.53％；2010－2015年为减速阶段，耕地的变化速率CCLi下降至0.23％。
耕地向其他地类的转变在空间上逐步扩散。通过图3～图5广东省土地利用变化空间分区，1990－2000年，耕地的转移高度集中于珠三角区；2000－2010年，耕地的转移集中于珠三角区但其他区域的耕地的转移亦有所增加；至2010－2015年，耕地的转移在广东省内分布已较为均匀。
城乡工矿居住用地土地利用类型变化基本特征：
城乡工矿居住用地总量持续增加，城乡工矿居住用地总量在1990年、2000年、2010年、2015年分别为46.43万公顷、61.66万公顷、94万公顷、103.77万公顷，而2000年、2010年、2015年城乡工矿居住用地转移面积仅分别为0.06万公顷、0.04万公顷、0.22万公顷，转移量远远低于新增量。
城乡工矿居住用地增加的速度由快减慢，结合数据汇总表可知，城乡工矿居住用地增加的速度经历了先加速后减速的过程：1990－2010年为加速阶段，城乡工矿居住用地的变化速率CCLi从1990－2000年的3.3％快速提高到2000－2010年的5.26％；2010－2015年为减速阶段，城乡工矿居住用地的变化速率CCLi回落至2.17％。
其他用地向城乡工矿居住用地的转移在空间上逐步扩散，1990－2000年，其他用地向城乡工矿居住用地的转移高度集中于珠三角区；2000－2010年，其他用地向城乡工矿居住用地的转移集中于珠三角区，但其他区域的其他用地向城乡工矿居住用地的转移亦有所增加；至2010－2015年，其他用地向城乡工矿居住用地的转移在广东省内分布已较为均匀。
非人工土地利用类型变化基本特征：
水域与各用地类型的转换集中发生在珠三角区域。水域面积总体增加，但经历了大量增加再逐渐减少的过程：1990－2000年，由耕地转水域增加面积5.58万公顷；2000－2010年，水域主要转为城乡工矿居住用地，耕地转水域的面积小于水域面积的转换，面积减少1.04万公顷；2010－2015年，水域主要转为城乡工矿居住用地，耕地转水域的面积小于水域面积的转换，面积减少1.1万公顷。水域的变化速率CCLi逐渐减缓，从1990－2000年的1.41％减缓至2000－2010年的0.94％、2010－2015年的0.36％林地总量持续减少，2000年、2010年、2015年林地的减少面积分别为0.37万公顷、6.14万公顷、6.2万公顷；减少的速度有增加趋势，变化速率CCLi从1990－2000年、2000－2010年的0.08％增加至2010－2015年的0.18％。林地向其他地类的转变主要发生在珠三角、粤西及粤西北区域，而其他地类向林地的转变1990－2010年主要发生在粤西区域，2010－2015年主要发生在珠三角和粤西北区域。
草地总量先减少后增加，总量基本维持不变。1990年、2000年、2010年草地的面积分别为83.41万公顷、81.51万公顷、81.43万公顷而2015年则增加至83.31万公顷，1990年与2015年相比仅相差0.1万公顷。草地面积的变化速率CCLi经历了减少的速度放缓到快速的面积增加的过程：1990－2000年的变化速率为0.33％放缓至2000－2010年的0.03％，2010－2015年提高至1.56％。草地的转移与新增在1990－2010年间主要发生在珠三角、粤西及粤西北区域，2010－2015年则于全省较为均匀分布。
未利用地土地利用类型变化基本特征
未利用地总量基本维持不变，2000年、2010年、2015年未利用地的面积增减分别为0.03万公顷、－0.38万公顷、0.01万公顷；变化的速度先增加后放缓，变化速率CCLi从1990－2000年的0.16％增加至2000－2010年的1.94％再放缓至2010－2015年的0.71％。未利用地与各用地类型的转换于全省内零星分布。
相应的，图9示出了本发明实施例中的分时序省级土地利用数据处理系统，包括
分析目标模块：用于选择分析省份及分析时间段；
子时序划分模块：用于在所述分析时间段中添加分隔节点时间，将所述分析时间段划分为一个以上的子时序；
土地利用数据提取模块：用于从全国土地利用数据库中提取所述分时序省级分析目标的土地利用数据；
土地利用类型分类标准模块：用于根据预设的土地利用类型分类标准，将所述分时序省级分析目标的土地利用数据分类；
转换模块：用于将所述分时序省级目标的土地利用数据中的二级类型的土地利用转换为一级类型的土地利用；
土地利用变化计算模块：用于计算所述土地利用类型在所述子时序结束时间相对于起始时间的土地利用变化类型，并对不同的土地利用变化类型分类；
土地利用变化影像构建模块：用于在分时序省级分析目标的地图，用不同颜色对不同类型的土地利用变化进行标示；
土地利用结构计算模块：用于计算土地利用结构；
单一土地利用动态度计算模块：用于计算单一土地利用动态度；
空间信息分析模型计算模块：用于计算空间信息分析模型；
制表储存模块：用于将所述土地利用数据，土地利用结构、单一土地利用动态度、空间信息分析模型信息汇总并制表储存；
本发明实施例提供的分时序省级土地利用数据处理方法及系统主要用于对省域土地利用数据进行处理，可较为直观的得出在目标分析时间段中，各子时序的土地利用变化数据，有利于快速获取土地利用变化类型的空间分布特征和土地利用变化的土地利用结构、单一土地利用动态度、空间信息分析模型三个主要参数的信息，加快了土地利用研究的数据获取速度。
以上对本发明实施例所提供的一种分时序省级土地利用数据处理方法及系统进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种分时序省级土地利用数据处理方法，其特征在于，包括以下步骤：
确定分时序省级分析目标和分析时间段；
在所述分析时间段中添加分隔节点时间，所述分隔节点时间将所述分析时间段划分为一个以上的子时序；
从全国土地利用数据库中提取所述分析时间段起始时间、结束时间和所述分隔节点时间的所述分时序省级分析目标的土地利用数据；
构建所述土地利用数据的土地利用类型分类标准；
在所述分时序省级分析目标的地图上，对在子时序结束时间相对于起始时间发生土地利用类型变化的区域进行标示；
按所述土地利用类型分类标准计算在所述分时序省级分析目标中，每个土地利用类型在所述分析时间段起始时间、结束时间和所述分隔节点时间的土地利用结构、单一土地利用动态度和空间信息分析模型；
将所述土地利用结构、单一土地利用动态度、空间信息分析模型信息汇总。
2.如权利要求1所述的分时序省级土地利用数据处理方法，其特征在于，所述土地利用数据为通过分辨率为30米的Landsat TM数字影像解译获取的分辨率为1km的土地利用分类栅格数据。
3.如权利要求2所述的分时序省级土地利用数据处理方法，其特征在于，所述土地利用类型分类标准包括6个一级类型和25个二级类型。
4.如权利要求3所述的分时序省级土地利用数据处理方法，其特征在于，所述方法还包括，将所述分时序省级分析目标的土地利用数据中的二级类型数据向所述一级类型进行转换。
5.如权利要求1所述的分时序省级土地利用数据处理方法，其特征在于，在所述分时序省级分析目标的地图上，所述土地利用变化类型相同的区域用同一颜色进行标示。
6.如权利要求1所述的分时序省级土地利用数据处理方法，其特征在于，所述土地利用结构用于表示第i种土地利用类型面积所占的比重，公式为Ai表示第i种土地利用类型的面积；n表示土地利用类型的类别数量。
7.如权利要求1所述的分时序省级土地利用数据处理方法，其特征在于，所述土地利用度用于描述第i种土地利用类型在[t1,t2]时间段内的变化速度，用单一土地利用动态度定量描述，公式为其中，LA(i,t1)和LA(i,t2)为研究期起始和结束时第i种土地利用类型的面积；t1,t2研究期起始时间和结束时间。
8.如权利要求1所述的分时序省级土地利用数据处理方法，其特征在于，所述空间信息分析模型包括转移速率TRLi、新增速率IRLi和变化速率CCLi，计算公式为
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>TRL</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>LA</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>ULA</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>LA</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>IRL</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>LA</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>ULA</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>LA</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>CCL</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>LA</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>ULA</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<msub>
<mi>LA</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>ULA</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>LA</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
<mo>=</mo>
<msub>
<mi>TRL</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>IRL</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，转移速率TRLi用于表示第i种土地利用类型在t1至t2期间转移面积与期初面积比重，IRLi用于表示第i种土地利用类型在t1至t2期间新增面积与期初面积比重的新增速率，CCLi第i种土地利用类型在t1至t2期间转移面积和新增面积之和与期初面积比重的变化速率；表示t1时间第i种土地利用类型的面积，表示t2时间第i种土地利用类型的面积，ULAi表示第i种土地利用类型在时序开始与结束没发生改变的面积。
9.一种分时序省级土地利用数据处理系统，其特征在于，包括
分析目标模块：用于选择分析目标省份及分析时间段；
子时序划分模块：用于在所述分析时间段中添加分隔节点时间，将所述分析时间段划分为一个以上的子时序；
土地利用数据提取模块：用于从全国土地利用数据库中提取所述分时序省级分析目标的土地利用数据；
土地利用类型分类标准模块：用于根据预设的土地利用类型分类标准，将所述分时序省级分析目标的土地利用数据分类；
土地利用变化计算模块：用于计算所述土地利用类型在所述子时序结束时间相对于起始时间的土地利用变化类型，并对不同的土地利用变化类型分类；
土地利用变化影像构建模块：用于在分时序省级分析目标的地图，用不同颜色对不同类型的土地利用变化进行标示；
土地利用结构计算模块：用于计算土地利用结构；
单一土地利用动态度计算模块：用于计算单一土地利用动态度；
空间信息分析模型计算模块：用于计算空间信息分析模型；
制表储存模块：用于将所述土地利用数据，土地利用结构、单一土地利用动态度、空间信息分析模型信息汇总并制表储存。
10.如权利要求9所述的分时序省级土地利用数据处理系统，其特征在于，所述系统还包括：
转换模块：用于将所述分时序省级目标的土地利用数据中的二级类型的土地利用转换为一级类型的土地利用。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711036172.4A CN107818519A (zh) | 2017-10-30 | 2017-10-30 | 一种分时序省级土地利用数据处理方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711036172.4A CN107818519A (zh) | 2017-10-30 | 2017-10-30 | 一种分时序省级土地利用数据处理方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107818519A true CN107818519A (zh) | 2018-03-20 |
Family
ID=61603446
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711036172.4A Pending CN107818519A (zh) | 2017-10-30 | 2017-10-30 | 一种分时序省级土地利用数据处理方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107818519A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108460372A (zh) * | 2018-04-25 | 2018-08-28 | 天津城建大学 | 一种生态用地的流失风险识别方法及系统 |
| CN110580539A (zh) * | 2018-06-11 | 2019-12-17 | 深圳市数字城市工程研究中心 | 基于土地时空模型反馈土地利用变化的方法及系统 |
| CN110737745A (zh) * | 2019-10-15 | 2020-01-31 | 中国地质大学（北京） | 一种土地利用动态变化地图的构建方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105184048A (zh) * | 2015-07-23 | 2015-12-23 | 河南大学 | 一种基于多源数据的逐年土地利用空间数据的制备方法 |
| JP2016051388A (ja) * | 2014-09-01 | 2016-04-11 | 日本放送協会 | 集計装置、集計方法及び集計プログラム |
-
2017
- 2017-10-30 CN CN201711036172.4A patent/CN107818519A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016051388A (ja) * | 2014-09-01 | 2016-04-11 | 日本放送協会 | 集計装置、集計方法及び集計プログラム |
| CN105184048A (zh) * | 2015-07-23 | 2015-12-23 | 河南大学 | 一种基于多源数据的逐年土地利用空间数据的制备方法 |
Non-Patent Citations (1)
| Title |
|---|
| 冯异星: "遥感与GIS支持下的土地利用动态变化研究-以鄂尔多斯市为例", 《中国优秀硕士学位论文全文数据库 经济与管理科学辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108460372A (zh) * | 2018-04-25 | 2018-08-28 | 天津城建大学 | 一种生态用地的流失风险识别方法及系统 |
| CN110580539A (zh) * | 2018-06-11 | 2019-12-17 | 深圳市数字城市工程研究中心 | 基于土地时空模型反馈土地利用变化的方法及系统 |
| CN110580539B (zh) * | 2018-06-11 | 2021-07-30 | 深圳市数字城市工程研究中心 | 基于土地时空模型反馈土地利用变化的方法及系统 |
| CN110737745A (zh) * | 2019-10-15 | 2020-01-31 | 中国地质大学（北京） | 一种土地利用动态变化地图的构建方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fan et al. | Infrastructure and regional economic development in rural China | |
| CN107818519A (zh) | 一种分时序省级土地利用数据处理方法及系统 | |
| CN105260390B (zh) | 一种面向群组的基于联合概率矩阵分解的项目推荐方法 | |
| Yang et al. | Impact of LUCC on landscape pattern in the Yangtze River Basin during 2001–2019 | |
| CN108549858B (zh) | 一种城市热岛效应的定量评价方法 | |
| CN105787488B (zh) | 由全局向局部传递的图像特征提取方法及装置 | |
| Boix et al. | The Roots of the Industrial Revolution: Political Institutions or (Socially Embedded) Know-How? | |
| Kaiyong et al. | The research on impact factors and characteristic of cultivated land resources use efficiency---take henan province, China as a case study | |
| Kuc-Czarnecka | Is the regional divergence a price for the international convergence? The case of the Visegrad group | |
| Garbin et al. | Spatial segregation of subordinate species is not controlled by the dominant species in a tropical coastal plant community | |
| CN107390187A (zh) | 基于混合高斯分布的杂波抑制方法 | |
| Xu et al. | Differentiation and progress of urban regionalization in China: Perspectives of land use and geography | |
| Funnell | Urban-rural linkages: Research themes and directions | |
| Zhang | Case dependence of multiscale interactions between multisource perturbations for convection-permitting ensemble forecasting during SCMREX | |
| CN102208027A (zh) | 基于间隙度维数的土地利用空间格局评价方法 | |
| Bandara et al. | An empirical analysis of Sri Lanka's Manufacturing Productivity slow-down | |
| Ghosh | Estimating the demographic dividend: Evidence from Indian states | |
| Hildebrandt | CSR und Sportmanagement | |
| Li et al. | Improved local accumulate histogram-based Thangka Image Retrieval | |
| CN111353628A (zh) | 一种中国主导树种标准生长指数研建方法 | |
| CN102708412A (zh) | 一种城乡生态规划的生态类因子分析方法 | |
| CN102236898A (zh) | 基于无限成分数的t混合模型的图像分割方法 | |
| CN117011698B (zh) | 多维度和多模型的地表全覆盖解译样本集评价方法 | |
| Yuan et al. | Researching the comprehensive use of fuzzy comprehensive evaluation and Markov chain | |
| CN108241593B (zh) | 一种空间抽样目标定位方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180320 |