CN116108312B - 一种基于单图斑成果数据的土地变化面积计算方法及装置 - Google Patents
一种基于单图斑成果数据的土地变化面积计算方法及装置 Download PDFInfo
- Publication number
- CN116108312B CN116108312B CN202211384671.3A CN202211384671A CN116108312B CN 116108312 B CN116108312 B CN 116108312B CN 202211384671 A CN202211384671 A CN 202211384671A CN 116108312 B CN116108312 B CN 116108312B
- Authority
- CN
- China
- Prior art keywords
- area
- land
- change
- data
- spot
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
Abstract
本发明公开了一种基于单图斑成果数据的土地变化面积计算方法及装置,所述方法包括:分别获取关于待调查土地的基础数据和单图斑成果数据;对所述基础数据与所述单图斑成果数据进行矢量联合叠加得到叠加数据,对所述叠加数据进行赋值计算得到待调查土地的变更面积;采用所述变更面积计算土地变化面积。本发明可以在对土地进行调查时,可以利用单图斑成果数据与土地的各类基础数据预测计算得到土地变化面积,可以在年度变更调查数据库成果出来之前,能够以单个图斑为对象,提前得到单图斑的土地利用现状数据变化及流向,实时掌握全省各区县的即时地类流向变化情况,满足对土地利用现状的提前预判和监测的预警。
Description
技术领域
本发明涉及预估土地变化面积技术领域,尤其涉及一种基于单图斑成果数据的土地变化面积计算方法及装置。
背景技术
图斑指在一定比例的调查底图上,将土地利用现状类型相同的地块分为一类,以其为基础调查单元,然后将单元勾绘成的图,图斑成果数据则是勾绘的图的相关数据。由于图斑成果数据包含调查图上关于土地的土地现状、土地权属、行政区划变化等相关数据,可利用各个数据确定土地的面积变化,供相关部门或单位进行土地分析、统计、规划和管理。
目前常用的土地面积变化的计算方法是:先确定待调查的土地区域,在调查图上将土地区域划分成多个图斑并由调查人员实地采集每个图斑的地类变化,土地权属、行政区划数据等数据。待调查人员将外业调查图斑积累到一定程度后,再根据调查时间分批内业建库,集中进行数据库质量检查最后统一增量数据库成果生成,将新建数据库与在先的数据库进行比较,以确定土地的变化面积。
但目前常用的方法有如下技术问题:调查与采集数据的耗时较长,当调查的图斑数量越多所需的时间越长,最后形成增量数据库成果的周期也就越长,计算土地变化面积的效率低。并且整个处理过程属于事后行为,在调查过程中若土地出现新的变化或其应用出现调整,需要重新进行数据采集与调查,不但增加了处理的工作量,而且时效性较低,难以满足相关部门的土地管理需求。
发明内容
本发明提出一种基于单图斑成果数据的土地变化面积计算方法及装置,所述方法在对土地进行调查时,分别采集土地的各类基础数据和一个图斑成果数据,将两种数据叠加并根据土地的类型对叠加后的数据进行赋值计算,得到土地变化面积,从而减少处理的工作量,缩短计算时长,提高计算的效率。
本发明实施例的第一方面提供了一种基于单图斑成果数据的土地变化面积计算方法,所述方法包括:
分别获取关于待调查土地的基础数据和单图斑成果数据,所述基础数据包括:土地利用现状数据、土地范围图层、土地坡度图、田坎系数;
对所述基础数据与所述单图斑成果数据进行矢量联合叠加得到叠加数据,对所述叠加数据进行赋值计算得到待调查土地的变更面积;
采用所述变更面积计算土地变化面积。
在第一方面的一种可能的实现方式中,所述对所述叠加数据进行赋值计算得到待调查土地的变更面积,包括:
确定所述叠加数据对应的土地属性,基于所述土地属性在预设的赋值规则表中查找对应的数据赋值;
采用所述数据赋值计算变更参数,所述变更参数包括变更前地类编码、变更后地类编码、变更前扣除地类编码、变更后扣除地类编码、变更后地类图斑面积、变更前地类图斑面积、变更后扣除地类面积和变更前扣除地类面积;
获取当前面积参数和预设的调整面积参数,确定所述当前面积参数与所述预设的调整面积参数的参数比较结果,根据所述参数比较结果采用所述变更参数计算变更面积。
在第一方面的一种可能的实现方式中,所述当前面积参数,包括:当前地类图斑面积和当前田坎面积,所述预设的调整面积参数,包括:变更后地类图斑面积和变更后田坎面积;
所述根据所述参数比较结果采用所述变更参数计算变更面积,包括:
若所述当前地类图斑面积≥所述变更后地类图斑面积,且所述当前田坎面积≤所述变更后田坎面积时,确定待调查土地的第一地类流向;
若所述第一地类流向为变更前地类流向变更后地类,则所述变更面积为所述变更后地类图斑面积;
若所述第一地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为所述变更前扣除地类面积;
若所述第一地类流向为变更前地类流向变更后扣除地类,则所述变更面积为所述变更前地类图斑面积减去所述变更后地类图斑面积。
在第一方面的一种可能的实现方式中,所述预设的调整面积参数,还包括:变更后田坎面积;
所述根据所述参数比较结果采用所述变更参数计算变更面积,还包括:
若所述当前地类图斑面积<所述变更后地类图斑面积,且所述当前田坎面积≥所述变更后田坎面积时,确定待调查土地的第二地类流向;
所述第二地类流向为变更前地类流向变更后地类,则所述变更面积为变更前地类图斑面积;
所述第二地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为变更后扣除地类面积;
所述第二地类流向为变更前扣除地类流向变更后地类,则所述变更面积为所述变更前扣除地类面积减去所述变更后扣除地类面积。
在第一方面的一种可能的实现方式中,所述数据赋值包括:后扣除地类系数、前扣除地类系数;
计算所述变更后扣除地类面积,具体为:
变更后扣除地类面积=图斑面积*后扣除地类系数;
计算所述变更前扣除地类面积,具体为:
所以前扣除地类面积=图斑面积*前扣除地类系数;
计算所述变更后地类图斑面积,具体为:
变更后地类图斑面积=图斑面积-变更后扣除地类面积;
计算所述变更前地类图斑面积,具体为:
变更前地类图斑面积=图斑面积-变更前扣除地类面积。
在第一方面的一种可能的实现方式中,计算所述图斑面积,具体为:
计算待调查土地对应图斑的原椭球面积,并采集待调查土地对应图斑的历史图斑面积;
叠加并裁剪所述原椭球面积与所述历史图斑面,得到处理椭球面积;
将所述处理椭球面积除以所述原椭球面积与所述历史图斑面积的乘积,得到图斑面积。
在第一方面的一种可能的实现方式中,计算所述土地变化面积具体为:
上式中,Kj为土地利用类型j的土地变化面积,Kj,a为土地利用类型j的监测图斑面积、Kj,b为土地利用类型j的变更面积,Kj,A为土地利用类型j的监测总图斑面积。
在第一方面的一种可能的实现方式中,在所述分别获取关于待调查土地的基础数据和单图斑成果数据的步骤后,所述方法还包括:
对所述基础数据和所述单图斑成果数据进行预处理,所述预处理包括:对所述土地利用现状数据和所述土地范围图层进行联合叠加,得到图层命名为基础库,以及根据所述土地坡度图和所述田坎系数对进行所述单图斑成果数据赋值处理。
在第一方面的一种可能的实现方式中,所述方法还包括:
在所述获取关于待调查土地的基础数据和单图斑成果数据的步骤后、在所述对所述基础数据和所述单图斑成果数据进行预处理的步骤后、在所述对所述叠加数据进行赋值计算得到待调查土地的变更面积的步骤后或者在所述采用所述变更面积计算土地变化面积的步骤后,进行校验分析,以确定数据的一致性。
本发明实施例的第二方面提供了一种基于单图斑成果数据的土地变化面积计算装置,所述装置包括:
获取模块,用于分别获取关于待调查土地的基础数据和单图斑成果数据,所述基础数据包括:土地利用现状数据、土地范围图层、土地坡度图、田坎系数;
叠加赋值模块,用于对所述基础数据与所述单图斑成果数据进行矢量联合叠加得到叠加数据,对所述叠加数据进行赋值计算得到待调查土地的变更面积;
计算模块,用于采用所述变更面积计算土地变化面积。
相比于现有技术,本发明实施例提供的一种基于单图斑成果数据的土地变化面积计算方法及装置,其有益效果在于:本发明可以在对土地进行调查时,分别采集土地的各类基础数据和一个图斑成果数据,将两种数据叠加并根据土地的类型对叠加后的数据进行赋值计算得到土地的变更面积,采用变更面积预测计算得到土地变化面积。可以在年度变更调查数据库成果出来之前,能够以单个图斑为对象,提前得到单图斑的土地利用现状数据变化及流向,实时掌握全省各区县的即时地类流向变化情况,满足对土地利用现状的提前预判和监测的预警。
附图说明
图1是本发明一实施例提供的一种基于单图斑成果数据的土地变化面积计算方法的流程示意图;
图2是本发明一实施例提供的叠加数据的示意图;
图3是本发明一实施例提供的一种基于单图斑成果数据的土地变化面积计算方法的操作流程图;
图4是本发明一实施例提供的一种基于单图斑成果数据的土地变化面积计算装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前常用的土地面积变化的计算方法是:先确定待调查的土地区域,在调查图上将土地区域划分成多个图斑并由调查人员实地采集每个图斑的地类变化,土地权属、行政区划数据等数据。待调查人员将外业调查图斑积累到一定程度后,再根据调查时间分批内业建库,集中进行数据库质量检查最后统一增量数据库成果生成,将新建数据库与在先的数据库进行比较,以确定土地的变化面积。
但目前常用的方法有如下技术问题:调查与采集数据的耗时较长,当调查的图斑数量越多所需的时间越长,最后形成增量数据库成果的周期也就越长,计算土地变化面积的效率低。并且整个处理过程属于事后行为,在调查过程中若土地出现新的变化或其应用出现调整,需要重新进行数据采集与调查,不但增加了处理的工作量,而且时效性较低,难以满足相关部门的土地管理需求。
为了解决上述问题,下面将通过以下具体的实施例对本申请实施例提供的一种基于单图斑成果数据的土地变化面积计算方法进行详细介绍和说明。
参照图1,示出了本发明一实施例提供的一种基于单图斑成果数据的土地变化面积计算方法的流程示意图。
其中,作为示例的,所述基于单图斑成果数据的土地变化面积计算方法,可以包括:
S11、分别获取关于待调查土地的基础数据和单图斑成果数据,所述基础数据包括:土地利用现状数据、土地范围图层、土地坡度图、田坎系数。
土地利用现状数据可以是土地利用现状数据库内的相关数据,其中,土地利用现状数据库是国土调查相关部门每10年进行一次全国土地调查;根据土地管理工作的需要,每年进行土地变更调查构成的数据库。其中可以包括土地利用现状及变化情况,包括地类、位置、面积、分布等状况,土地权属及变化情况和土地条件等内容。土地利用现状数据库记录每个地块的土地利用现状情况。
在一可选的实施例中,土地利用现状的分类可以涉及12个一级类73个二级类,地类区分的精细度越高,对数据要求和精确度越高。
土地范围图层可以是城镇村范围图层,具体可以包括根据国土要求所划定的城市、建制镇、村庄范围和边界。
土地坡度图(又简称坡度图)可以是表示地面倾斜率的地图。具体可以是根据国土调查坡度图生产技术规定,用DEM生产坡度栅格图,进行坡度分级,图形综合等处理形成全省坡度图成果。
田坎系数可以是田坎系数表,具体可以是计算耕地范围内的田坎面积,田坎系数由省(区、市)统一组织测算、测算方案以及测算结果。
单图斑成果数据可以是以单个图斑为对象的相关成果数据,具体可以是通过外业调查举证、单图斑建库、单图斑自检、单图斑成果县级审核、单图斑省级审核,待所有单图斑成果通过审核后由省级统一下发地方,地方整合建库经省级核查通过后进行上报的成果数据。单图斑成果数据主要包含外业图斑是否变更、变化范围、变更后地类、未变更原因等外业图斑变更行为,以及需要变更的图斑绘制变更后的图形范围。
由于上述各个数据可能格式或者内容暂未符合后续的计算需求,可以对其进行预处理,使其满足后续的计算要求。
其中,作为示例的,在步骤S11后,所述方法还可以包括以下步骤:
S21、对所述基础数据和所述单图斑成果数据进行预处理,所述预处理包括:对所述土地利用现状数据和所述土地范围图层进行联合叠加,得到图层命名为基础库,以及根据所述土地坡度图和所述田坎系数对进行所述单图斑成果数据赋值处理。
在一预处理的实现方式中,可以将土地利用现状数据库叠加到城镇村范围层,具体可以利用基础库模型算法。
具体地,将土地利用现状数据库中的地类图斑层与城镇村范围层进行联合叠加,得到图层命名为基础库(JCK)。
在一预处理的实现方式中,可以根据单图斑成果数据对坡度进行赋值和对田坎系数进行赋值处理。
由于单图斑成果数据中缺少图斑坡度属性和田坎系数属性,需要根据坡度图对单图斑成果进行坡度赋值,根据单图斑成果中的坡度和耕地类型配对田坎系数表进行田坎系数赋值。
首先,单图斑成果坡度赋值计算图斑内优势坡度等级,按面积比例占比最大的坡度等级作为该图斑的坡度等级,坡度赋值属性字段名为单图斑坡度(DTBPD)。具体可以采用现有的坡度赋值模型算法。
然后,根据单图斑坡度和耕地类型,结合田坎系数表对耕地图斑进行赋值田坎系数,田坎系数赋值属性字段名为单图斑田坎系数(DTBKCXS)。具体可以采用现有的田坎系数模型算法。
通过预处理,可以让基础数据和单图斑成果数据符合后续的计算要求,以提高计算的准确率。
S12、对所述基础数据与所述单图斑成果数据进行矢量联合叠加得到叠加数据,对所述叠加数据进行赋值计算得到待调查土地的变更面积。
参照图2,示出了本发明一实施例提供的叠加数据的示意图。
在一实施例中,可以将数据进行叠加分析,叠加后如图2所示。
具体地,可以将在上述预处理后的成果将基础库与单图斑成果进行联合叠加分析,得到叠加数据,该叠加数据具体可以为图层命名为年末库(nmk)。
在完成叠加后,可以对叠加数据进行赋值与计算,从而得到待调查土地的变更面积。
在一可选的实施例中,步骤S12可以包括以下子步骤:
S121、确定所述叠加数据对应的土地属性,基于所述土地属性在预设的赋值规则表中查找对应的数据赋值。
在一实施例中,土地属性可以是其字段名称,可以获取叠加数据内每个字段的字段名称,然后按照其字段名称在赋值规则表中查找对应的字段代码,在根据字段代码查找对应的阈值。阈值具体可以根据实际需求进行调整。
在一实施方式中,预设的赋值规则表可以如下表所示:
在一可选的实施例中,具体字段赋值逻辑如下:
前地类编码(QDLBM):使用基础库中的地类编码。
前耕地类型(QGDLX):使用基础库中的耕地类型。
前扣除地类编码(QCKDLBM):使用基础库中的扣除地类编码。
前扣除地类系数(QKCXS):使用基础库中的扣除地类系数。
前耕地坡度级别(QGDPDJB):使用基础库中的耕地坡度级别。
前城镇村属性码(QCZCSXM):使用基础库叠加城镇村范围层的城镇村类型代码。
前种植属性码(QZZSXDM):使用基础库中种植属性代码。
前图斑细化代码(QTBXHDM):使用基础库中图斑细化代码。
后地类编码(HDLBM):使用单图斑层中地类编码更新QDLBM。
后耕地类型(HGDLX):后地类编码为耕地的图斑,用单图层中的耕地类型更新QGDLLX。
后扣除地类编码(HCKDLBM):后耕地类型为梯田的,将后扣除地类编码赋值为1203即田坎。
后扣除地类系数(HKCXS):使用单图斑田坎系数更新基础库中的扣除地类系数。
后耕地坡度级别(HGDPDJB):使用单图斑更新前耕地坡度级别。
后种植属性码(HZZSXDM):使用单图斑中的种植属性码更新前种植属性码。
后图斑细化代码(HTBXHDM):使用单图斑中的图斑细化代码更新前图斑细化代码。
其中,后城镇村属性码(QCZCSXM):根据变更调查规则具体赋值逻辑如下:
第一:当单图斑图层中标注城市(201、201A)、建制镇(202、202A)及村庄(203、203A)的地类图斑,标注工矿用地(204)和风景名胜特殊用地(205)的建设用地图斑,用单图斑的城镇村代码更新前城镇村属性码。
第二:当后城镇村属性码为工矿用地(204)和风景名胜特殊用地(205)的后地类编码为非建设用地部分,后城镇村属性码赋值为空。
第三:当后地类编码为铁路用地、轨道交通用地、公路用地、机场用地、港口码头用地、管地运输用地、水工建筑用地的图斑,后城镇村属性码赋值为空。
S122、采用所述数据赋值计算变更参数,所述变更参数包括:变更前地类编码、变更后地类编码、变更前扣除地类编码、变更后扣除地类编码、变更后地类图斑面积、变更前地类图斑面积、变更后扣除地类面积和变更前扣除地类面积。
在一实现方式中,所述数据赋值包括:后扣除地类系数、前扣除地类系数;
具体地,计算所述变更后扣除地类面积,具体为:
变更后扣除地类面积=图斑面积*后扣除地类系数。
具体地,计算所述变更前扣除地类面积,具体为:
所以前扣除地类面积=图斑面积*前扣除地类系数。
具体地,计算所述变更后地类图斑面积,具体为:
变更后地类图斑面积=图斑面积-变更后扣除地类面积。
具体地,计算所述变更前地类图斑面积,具体为:
变更前地类图斑面积=图斑面积-变更前扣除地类面积。
其中,作为示例的,计算所述图斑面积,具体为:
S31、计算待调查土地对应图斑的原椭球面积,并采集待调查土地对应图斑的历史图斑面积。
S32、叠加并裁剪所述原椭球面积与所述历史图斑面,得到处理椭球面积。
S33、将所述处理椭球面积除以所述原椭球面积与所述历史图斑面积的乘积,得到图斑面积。
在一实现方式中,图斑面积可以按比例分割面积的方法计算得到。
在具体实现时,可以先计算该待调查土地对应的图斑的原椭球面积,可以新建字段为原椭球面积Y,其中,新增字段不分割面积,作为图形属性,记录图斑分割前的图形面积。(需要说明的是,原椭球面积是矢量数据,具体可以根据《全国国土调查技术规程附录D图幅理论面积与图斑椭球面积计算公式及要求》进行计算,也可以借助常用的空间数据库软件进行计算,如arcgis软件)。
接着,可以计算历史图斑面积,该历史图斑面积可以是待调查土地原来所对应的图斑的面积,可以记为YTBMJ。
然后,叠加并裁剪所述原椭球面积与所述历史图斑面积,再重新计算被裁切后的新图斑椭球面积,得到处理椭球面积H。
最后,可以计算图斑面积,具体可以如下式所示:
图斑面积=新椭球面积H/原椭球面积Y*原图斑面积YTBMJ。
在计算时,可以将每块图斑重新进行图斑面积计算。需要说明的是,在计算图斑面积的过程中,如裁切后形成的多个要素的“图斑面积”之后不等于原始的“图斑面积”,以将差值平到面积最大的要素中。具体就是将差异值添加到最大面积中,为了保证图形分割前后面积保持一致。
S123、获取当前面积参数和预设的调整面积参数,确定所述当前面积参数与所述预设的调整面积参数的参数比较结果,根据所述参数比较结果采用所述变更参数计算变更面积。
在一实施例中,在计算变更参数后,可以获取待调查土地的面积参数,包括当前面积参数和预设的调整面积参数,将两个面积参数进行比较,再基于其比较结果确定土地的变更情况,进而能根据请变更情况采用变更参数进行变更面积的计算,确定具体的变更面积。
在一实施例中,所述当前面积参数,包括:当前地类图斑面积和当前田坎面积,所述预设的调整面积参数,包括:变更后地类图斑面积和变更后田坎面积。
需要说明的是,当前地类图斑面积是待调查土地实时检测到的地类图斑面积,当前田坎面积可以是根据待调查土地的当前田坎系数计算耕地范围内的田坎面积,具体可以通过技术人员实地调查后计算得到。变更后地类图斑面积则为上述步骤计算的变更后地类图斑面积,变更后田坎面积是根据待调查土地在变更后的田坎系数计算耕地范围内的田坎面积。
参照下表,变更前为A地,变更后为B地和C地。
土地在变更的过程中,可能变更成不同的类型,使得土地的变更流向也不同,其变更的面积也不同,因此,需要根据其变更流向确定变更面积。
其中,作为示例的,步骤S123可以包括以下子步骤:
S1231、若所述当前地类图斑面积≥所述变更后地类图斑面积,且所述当前田坎面积≤所述变更后田坎面积时,确定待调查土地的第一地类流向。
S1232、若所述第一地类流向为变更前地类流向变更后地类,则所述变更面积为所述变更后地类图斑面积。
S1233、若所述第一地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为所述变更前扣除地类面积。
S1234、若所述第一地类流向为变更前地类流向变更后扣除地类,则所述变更面积为所述变更前地类图斑面积减去所述变更后地类图斑面积。
在一实施例中,所述预设的调整面积参数,还包括:变更后田坎面积。
其中,变更后田坎面积可以是根据待调查土地的调整后的田坎系数计算耕地范围内的田坎面积。
其中,作为示例的,步骤S123可以包括以下子步骤:
S1235、若所述当前地类图斑面积<所述变更后地类图斑面积,且所述当前田坎面积≥所述变更后田坎面积时,确定待调查土地的第二地类流向。
S1236、所述第二地类流向为变更前地类流向变更后地类,则所述变更面积为变更前地类图斑面积。
S1237、所述第二地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为变更后扣除地类面积。
S1238、所述第二地类流向为变更前扣除地类流向变更后地类,则所述变更面积为所述变更前扣除地类面积减去所述变更后扣除地类面积。
其变更情况如下表所示:
S13、采用所述变更面积计算土地变化面积。
在上表基础上,根据具体分析需要,可以汇总统计不同地类流向数据。以县级为单位开展数据汇总工作,汇总年度内每块土地的利用变化情况,并按照地(市)级、省级依次统计汇总方式,形成年度各级土地利用变化汇总成果。如年度国土变更调查流向重要地类统计表。
根据地方上报单图斑成果可知国家下发监测图斑范围内实际变更地类和面积,通过下发监测图斑面积和实际变更面积的比例得到监测图斑变更率,用监测图斑变更率乘以总监测图斑面积得到某类土地利用变化预测面积,从而实现对土地利用变化提前预估,由于单图斑成果上报实时动态更新,土地利用变化也实现了提前预估与实时预测。
在其中一种的实施例中,计算所述土地变化面积具体为:
/>
上式中,Kj为土地利用类型j的土地变化面积,Kj,a为土地利用类型j的监测图斑面积、Kj,b为土地利用类型j的变更面积,Kj,A为土地利用类型j的监测总图斑面积。
为了进一步加强对各个数据的监管质量,其中,作为示例的,在步骤S11、S12、S13或S21后,所述方法还可以包括以下步骤:
S14、进行校验分析,以确定数据的一致性。
具体地,在分别获取关于待调查土地的基础数据和单图斑成果数据后,可以对基础数据和单图斑成果数据进行校验分析,以确定基础数据和单图斑成果数据是否符合数据要求,若符合,则执行后续的矢量联合叠加操作,反之,若不符合要求,则重新采集基础数据和单图斑成果数据。
同理,在完成预处理后,可以对完成预处理的数据进行数据校验分析,确定预处理数据是否符合数据要求,若符合要求,则执行矢量联合叠加操作,若不符合,则重新采集基础数据和单图斑成果数据。
同理,在完成矢量联合叠加操作并计算得到待调查土地的变更面积后,也可以对变更面积进行校验分析,确定计算的变更面积是否符合数据要求,若符合,则采用变更面积计算土地变化面积,反之,若不符合要求,则重新采集基础数据和单图斑成果数据。
同理,在计算得到土地变化面积后,也可以对土地变化面积进行校验分析,确定计算的土地变化面积是否符合数据要求,若符合,则向用户展示土地变化面积,反之,若不符合要求,则重新采集基础数据和单图斑成果数据。
各个数据要求可以根据实际需要进行调整,例如,变更面积可以设置一个阈值,判断变更面积是否大于对应的阈值,若是,则符合要求,反之,则不符合要求。
本发明在完成每一个步骤后,均对步骤计算得到数据进行校验,校验不通过需查找原因并返回上一环节进行重新分析。校验主要方法是数据分析叠加前后面积保持一致,数据统计汇总前后总面积保持一致。
参照图3,示出了本发明一实施例提供的一种基于单图斑成果数据的土地变化面积计算方法的操作流程图。
在一具体的操作方式中,可以先采集相关数据,包括单图斑成果数据、土地利用现状数据、土地范围图层、土地坡度图、田坎系数。接着,对上述各个数据进行预处理,得到预处理数据。然后,对预处理数据进行分析,分析操作包括:联合叠加、算法运行(具体包括:赋值计算)以及构建表。最后,对表内的数据进行汇总与计算,得到预测的土地变化面积。
本发明可以利用高效准确的土地利用现状情况,结合土地利用地类变化及流向信息,对土地利用情况实现提前预估、对土地利用地类流量实现提前预判、对异常变更地类实现提前预警。及时进行国土变更调查变更过程实时监测,进一步加强国土变更调查成果质量全流程监管。
而且结合土地调查数据库、坡度图、田坎系数表等数据,利用单图斑成果数据进行土地利用现状变化及流向预估主要分为四个阶段,分别是数据收集、数据预处理、数据分析和成果输出等阶段。为了保证数据分析成果质量,在每一环节结束后都需要进行校验,保证数据的正确性。
本发明可以在年度变更调查数据库成果出来之前,能够以单个图斑为对象,提前得到单图斑的土地利用现状数据变化及流向,实时掌握全省各区县的即时地类流向变化情况,满足对土地利用现状的提前预判和监测的预警。可以在年度变更调查数据库成果出来之前,能够以单个图斑为对象,提前得到单图斑的土地利用现状数据变化及流向,实时掌握全省各区县的即时地类流向变化情况,满足对土地利用现状的提前预判和监测的预警。
同时本发明使用的土地利用现状,其分类共涉及12个一级类73个二级类,地类区分的精细度越高,对数据要求和精确度越高。结合国土调查数据库更新变更规则,能全面梳理了不同地类流向和流量统计逻辑,规范不同地类流向的计算方法,在基于对单图斑成果数据进行土地利用变化及流向预估的前提下,保障了数据分析成果的计算精度,以提高计算的准确率。
在本实施例中,本发明实施例提供了一种基于单图斑成果数据的土地变化面积计算方法,其有益效果在于:本发明可以在对土地进行调查时,分别采集土地的各类基础数据和一个图斑成果数据,将两种数据叠加并根据土地的类型对叠加后的数据进行赋值计算得到土地的变更面积,采用变更面积预测计算得到土地变化面积。整个过程无需采集多个图斑数据,也可以避免因调整土地应用而重新进行数据采集的工作,从而减少处理的工作量,大大缩短计算时长,提高预测土地变化面积的处理效率。
本发明实施例还提供了一种基于单图斑成果数据的土地变化面积计算装置,参见图4,示出了本发明一实施例提供的一种基于单图斑成果数据的土地变化面积计算装置的结构示意图。
其中,作为示例的,所述基于单图斑成果数据的土地变化面积计算装置可以包括:
获取模块401,用于分别获取关于待调查土地的基础数据和单图斑成果数据,所述基础数据包括:土地利用现状数据、土地范围图层、土地坡度图、田坎系数;
叠加赋值模块402,用于对所述基础数据与所述单图斑成果数据进行矢量联合叠加得到叠加数据,对所述叠加数据进行赋值计算得到待调查土地的变更面积;
计算模块403,用于采用所述变更面积计算土地变化面积。
可选地,所述叠加赋值模块,还用于:
确定所述叠加数据对应的土地属性,基于所述土地属性在预设的赋值规则表中查找对应的数据赋值;
采用所述数据赋值计算变更参数,所述变更参数包括:变更前地类编码、变更后地类编码、变更前扣除地类编码、变更后扣除地类编码、变更后地类图斑面积、变更前地类图斑面积、变更后扣除地类面积和变更前扣除地类面积;
获取当前面积参数和预设的调整面积参数,确定所述当前面积参数与所述预设的调整面积参数的参数比较结果,根据所述参数比较结果采用所述变更参数计算变更面积。
可选地,所述当前面积参数,包括:当前地类图斑面积和当前田坎面积,所述预设的调整面积参数,包括:变更后地类图斑面积和变更后田坎面积;
所述叠加赋值模块,还用于:
若所述当前地类图斑面积≥所述变更后地类图斑面积,且所述当前田坎面积≤所述变更后田坎面积时,确定待调查土地的第一地类流向;
若所述第一地类流向为变更前地类流向变更后地类,则所述变更面积为所述变更后地类图斑面积;
若所述第一地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为所述变更前扣除地类面积;
若所述第一地类流向为变更前地类流向变更后扣除地类,则所述变更面积为所述变更前地类图斑面积减去所述变更后地类图斑面积。
可选地,所述预设的调整面积参数,还包括:变更后田坎面积;
所述叠加赋值模块,还用于:
若所述当前地类图斑面积<所述变更后地类图斑面积,且所述当前田坎面积≥所述变更后田坎面积时,确定待调查土地的第二地类流向;
所述第二地类流向为变更前地类流向变更后地类,则所述变更面积为变更前地类图斑面积;
所述第二地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为变更后扣除地类面积;
所述第二地类流向为变更前扣除地类流向变更后地类,则所述变更面积为所述变更前扣除地类面积减去所述变更后扣除地类面积。
可选地,所述数据赋值包括:后扣除地类系数、前扣除地类系数;
计算所述变更后扣除地类面积,具体为:
变更后扣除地类面积=图斑面积*后扣除地类系数;
计算所述变更前扣除地类面积,具体为:
所以前扣除地类面积=图斑面积*前扣除地类系数;
计算所述变更后地类图斑面积,具体为:
变更后地类图斑面积=图斑面积-变更后扣除地类面积;
计算所述变更前地类图斑面积,具体为:
变更前地类图斑面积=图斑面积-变更前扣除地类面积。
可选地,计算所述图斑面积,具体为:
计算待调查土地对应图斑的原椭球面积,并采集待调查土地对应图斑的历史图斑面积;
叠加并裁剪所述原椭球面积与所述历史图斑面,得到处理椭球面积;
将所述处理椭球面积除以所述原椭球面积与所述历史图斑面积的乘积,得到图斑面积。
可选地,计算所述土地变化面积具体为:
上式中,Kj为土地利用类型j的土地变化面积,Kj,a为土地利用类型j的监测图斑面积、Kj,b为土地利用类型j的变更面积,Kj,A为土地利用类型j的监测总图斑面积。
可选地,所述装置还包括:
预处理模块,用于对所述基础数据和所述单图斑成果数据进行预处理,所述预处理包括:对所述土地利用现状数据和所述土地范围图层进行联合叠加,得到图层命名为基础库,以及根据所述土地坡度图和所述田坎系数对进行所述单图斑成果数据赋值处理。
可选地,所述装置还包括:
校验分析模块,用于在所述获取关于待调查土地的基础数据和单图斑成果数据的步骤后、在所述对所述基础数据和所述单图斑成果数据进行预处理的步骤后、在所述对所述叠加数据进行赋值计算得到待调查土地的变更面积的步骤后或者在所述采用所述变更面积计算土地变化面积的步骤后,进行校验分析,以确定数据的一致性。
所属技术领域的技术人员可以清楚地了解到,为方便的描述和简洁,上述描述的装置的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
进一步的,本申请实施例还提供了一种电子设备,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述实施例所述的基于单图斑成果数据的土地变化面积计算方法。
进一步的,本申请实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行程序,所述计算机可执行程序用于使计算机执行如上述实施例所述的基于单图斑成果数据的土地变化面积计算方法。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (7)
1.一种基于单图斑成果数据的土地变化面积计算方法,其特征在于,所述方法包括:
分别获取关于待调查土地的基础数据和单图斑成果数据,所述基础数据包括:土地利用现状数据、土地范围图层、土地坡度图、田坎系数;
对所述基础数据与所述单图斑成果数据进行矢量联合叠加得到叠加数据,对所述叠加数据进行赋值计算得到待调查土地的变更面积;
采用所述变更面积计算土地变化面积;
所述对所述叠加数据进行赋值计算得到待调查土地的变更面积,包括:
确定所述叠加数据对应的土地属性,基于所述土地属性在预设的赋值规则表中查找对应的数据赋值;
采用所述数据赋值计算变更参数,所述变更参数包括变更前地类编码、变更后地类编码、变更前扣除地类编码、变更后扣除地类编码、变更后地类图斑面积、变更前地类图斑面积、变更后扣除地类面积和变更前扣除地类面积;
获取当前面积参数和预设的调整面积参数,确定所述当前面积参数与所述预设的调整面积参数的参数比较结果,根据所述参数比较结果采用所述变更参数计算变更面积;
所述当前面积参数,包括:当前地类图斑面积和当前田坎面积,所述预设的调整面积参数,包括:变更后地类图斑面积和变更后田坎面积;
所述根据所述参数比较结果采用所述变更参数计算变更面积,包括:
若所述当前地类图斑面积≥所述变更后地类图斑面积,且所述当前田坎面积≤所述变更后田坎面积时,确定待调查土地的第一地类流向;
若所述第一地类流向为变更前地类流向变更后地类,则所述变更面积为所述变更后地类图斑面积;
若所述第一地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为所述变更前扣除地类面积;
若所述第一地类流向为变更前地类流向变更后扣除地类,则所述变更面积为所述变更前地类图斑面积减去所述变更后地类图斑面积;
所述预设的调整面积参数,还包括:变更后田坎面积;
所述根据所述参数比较结果采用所述变更参数计算变更面积,还包括:
若所述当前地类图斑面积<所述变更后地类图斑面积,且所述当前田坎面积≥所述变更后田坎面积时,确定待调查土地的第二地类流向;
所述第二地类流向为变更前地类流向变更后地类,则所述变更面积为变更前地类图斑面积;
所述第二地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为变更后扣除地类面积;
所述第二地类流向为变更前扣除地类流向变更后地类,则所述变更面积为所述变更前扣除地类面积减去所述变更后扣除地类面积。
2.根据权利要求1所述的基于单图斑成果数据的土地变化面积计算方法,其特征在于,所述数据赋值包括:后扣除地类系数、前扣除地类系数;
计算所述变更后扣除地类面积,具体为:
变更后扣除地类面积=图斑面积*后扣除地类系数;
计算所述变更前扣除地类面积,具体为:
所以前扣除地类面积=图斑面积*前扣除地类系数;
计算所述变更后地类图斑面积,具体为:
变更后地类图斑面积=图斑面积-变更后扣除地类面积;
计算所述变更前地类图斑面积,具体为:
变更前地类图斑面积=图斑面积-变更前扣除地类面积。
3.根据权利要求2所述的基于单图斑成果数据的土地变化面积计算方法,其特征在于,计算所述图斑面积,具体为:
计算待调查土地对应图斑的原椭球面积,并采集待调查土地对应图斑的历史图斑面积;
叠加并裁剪所述原椭球面积与所述历史图斑面,得到处理椭球面积;
将所述处理椭球面积除以所述原椭球面积与所述历史图斑面积的乘积,得到图斑面积。
4.根据权利要求1所述的基于单图斑成果数据的土地变化面积计算方法,其特征在于,计算所述土地变化面积具体为:
上式中,Kj为土地利用类型j的土地变化面积,Kj,a为土地利用类型j的监测图斑面积、Kj,b为土地利用类型j的变更面积,Kj,A为土地利用类型j的监测总图斑面积。
5.根据权利要求1-4任意一项所述的基于单图斑成果数据的土地变化面积计算方法,其特征在于,在所述分别获取关于待调查土地的基础数据和单图斑成果数据的步骤后,所述方法还包括:
对所述基础数据和所述单图斑成果数据进行预处理,所述预处理包括:对所述土地利用现状数据和所述土地范围图层进行联合叠加,得到图层命名为基础库,以及根据所述土地坡度图和所述田坎系数对进行所述单图斑成果数据赋值处理。
6.根据权利要求5所述的基于单图斑成果数据的土地变化面积计算方法,其特征在于,所述方法还包括:
在所述获取关于待调查土地的基础数据和单图斑成果数据的步骤后、在所述对所述基础数据和所述单图斑成果数据进行预处理的步骤后、在所述对所述叠加数据进行赋值计算得到待调查土地的变更面积的步骤后或者在所述采用所述变更面积计算土地变化面积的步骤后,进行校验分析,以确定数据的一致性。
7.一种基于单图斑成果数据的土地变化面积计算装置,其特征在于,所述装置包括:
获取模块,用于分别获取关于待调查土地的基础数据和单图斑成果数据,所述基础数据包括:土地利用现状数据、土地范围图层、土地坡度图、田坎系数;
叠加赋值模块,用于对所述基础数据与所述单图斑成果数据进行矢量联合叠加得到叠加数据,对所述叠加数据进行赋值计算得到待调查土地的变更面积;
计算模块,用于采用所述变更面积计算土地变化面积;
所述对所述叠加数据进行赋值计算得到待调查土地的变更面积,包括:
确定所述叠加数据对应的土地属性,基于所述土地属性在预设的赋值规则表中查找对应的数据赋值;
采用所述数据赋值计算变更参数,所述变更参数包括变更前地类编码、变更后地类编码、变更前扣除地类编码、变更后扣除地类编码、变更后地类图斑面积、变更前地类图斑面积、变更后扣除地类面积和变更前扣除地类面积;
获取当前面积参数和预设的调整面积参数,确定所述当前面积参数与所述预设的调整面积参数的参数比较结果,根据所述参数比较结果采用所述变更参数计算变更面积;
所述当前面积参数,包括:当前地类图斑面积和当前田坎面积,所述预设的调整面积参数,包括:变更后地类图斑面积和变更后田坎面积;
所述根据所述参数比较结果采用所述变更参数计算变更面积,包括:
若所述当前地类图斑面积≥所述变更后地类图斑面积,且所述当前田坎面积≤所述变更后田坎面积时,确定待调查土地的第一地类流向;
若所述第一地类流向为变更前地类流向变更后地类,则所述变更面积为所述变更后地类图斑面积;
若所述第一地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为所述变更前扣除地类面积;
若所述第一地类流向为变更前地类流向变更后扣除地类,则所述变更面积为所述变更前地类图斑面积减去所述变更后地类图斑面积;
所述预设的调整面积参数,还包括:变更后田坎面积;
所述根据所述参数比较结果采用所述变更参数计算变更面积,还包括:
若所述当前地类图斑面积<所述变更后地类图斑面积,且所述当前田坎面积≥所述变更后田坎面积时,确定待调查土地的第二地类流向;
所述第二地类流向为变更前地类流向变更后地类,则所述变更面积为变更前地类图斑面积;
所述第二地类流向为变更前扣除地类流向变更后扣除地类,则所述变更面积为变更后扣除地类面积;
所述第二地类流向为变更前扣除地类流向变更后地类,则所述变更面积为所述变更前扣除地类面积减去所述变更后扣除地类面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211384671.3A CN116108312B (zh) | 2022-11-07 | 2022-11-07 | 一种基于单图斑成果数据的土地变化面积计算方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211384671.3A CN116108312B (zh) | 2022-11-07 | 2022-11-07 | 一种基于单图斑成果数据的土地变化面积计算方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116108312A CN116108312A (zh) | 2023-05-12 |
CN116108312B true CN116108312B (zh) | 2023-09-19 |
Family
ID=86258610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211384671.3A Active CN116108312B (zh) | 2022-11-07 | 2022-11-07 | 一种基于单图斑成果数据的土地变化面积计算方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116108312B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107516282A (zh) * | 2017-08-01 | 2017-12-26 | 青岛市房地产信息与交易资金监管中心 | 三查合一的土地变更调查信息管理系统 |
CN109739943A (zh) * | 2018-12-14 | 2019-05-10 | 中国测绘科学研究院 | 面向自然资源矢量地表覆盖变化统计处理方法 |
CN112948759A (zh) * | 2021-03-04 | 2021-06-11 | 中国科学院生态环境研究中心 | 一种生态环境质量变化情况的判断方法及系统 |
-
2022
- 2022-11-07 CN CN202211384671.3A patent/CN116108312B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107516282A (zh) * | 2017-08-01 | 2017-12-26 | 青岛市房地产信息与交易资金监管中心 | 三查合一的土地变更调查信息管理系统 |
CN109739943A (zh) * | 2018-12-14 | 2019-05-10 | 中国测绘科学研究院 | 面向自然资源矢量地表覆盖变化统计处理方法 |
CN112948759A (zh) * | 2021-03-04 | 2021-06-11 | 中国科学院生态环境研究中心 | 一种生态环境质量变化情况的判断方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN116108312A (zh) | 2023-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Aburas et al. | Land suitability analysis of urban growth in Seremban Malaysia, using GIS based analytical hierarchy process | |
CN102434210B (zh) | 地下工程画像信息与监测信息安全监控的方法与系统 | |
CN102122395B (zh) | 一种保持地形特征的自适应尺度dem建模方法 | |
CN106547853A (zh) | 基于一张图的林业大数据建设方法 | |
CN105630988A (zh) | 一种快速检测空间数据变化并更新的方法及系统 | |
CN104778369A (zh) | 一种基于地面沉降监测的决策与预警方法及其系统 | |
US20160003792A1 (en) | Functional soil maps | |
CN113487123B (zh) | 高光谱监测与gis耦合山洪灾害动态风险评估方法 | |
CN109948737A (zh) | 基于大数据和机器学习的贫困空间分类识别方法和装置 | |
CN102184423B (zh) | 一种全自动的区域不透水面遥感信息精确提取方法 | |
CN111539904B (zh) | 一种基于降雨量的灾害易发性预测方法 | |
CN107119657A (zh) | 一种基于视觉测量基坑监测方法 | |
Guo et al. | Monitoring and simulation of dynamic spatiotemporal land use/cover changes | |
Belyakov et al. | Model Of Intellectual Visualization Of Geoinformation Service. | |
CN113340286A (zh) | 一种土地规划勘察勘测项目测绘信息数据分析方法、设备及计算机存储介质 | |
CN113487105A (zh) | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 | |
Zhao et al. | Combining ICESat-2 photons and Google Earth Satellite images for building height extraction | |
CN112923900B (zh) | 一种基于rs技术的建筑物不均匀沉降的判断方法及装置 | |
CN116108312B (zh) | 一种基于单图斑成果数据的土地变化面积计算方法及装置 | |
Rastogi et al. | Comparison of different spatial interpolation techniques to thematic mapping of socio-economic causes of crime against women | |
Oduro Appiah et al. | Measuring forest change patterns from oil and gas land use dynamics in northeastern British Columbia, 1975 to 2017 | |
KR20200052500A (ko) | 고해상도 드론 영상을 활용한 산림지역 변화 탐지 및 변화지역 투수/불투수도 갱신 방법 | |
CN109685895B (zh) | 一种适用于智能电网建设面向对象的地图地形均值方法 | |
CN113807724B (zh) | 一种基于综合风险评价的弃碴场选址方法 | |
Cai et al. | A schema of ecological environment sensitivity evaluation based on GIS |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |