福建省196地质大队 福建 漳州 363000
第三次全国国土调查(以下简称“三调”)的主要目标是在第二次全国土地调查成果基础上,全面细化和完善土地利用基础数据,掌握翔实准确的土地利用现状,满足生态文明建设、空间规划编制、自然资源管理体制改革和统一确权登记、国土空间规划体系建设等各项工作的需要。其中耕地作为专项调查内容之一,不仅要计算耕地图斑地类面积,还要计算耕地净面积。坡度是计算耕地净面积的重要参考因素,如何准确快速的实现坡度赋值对三调工作具有重要意义。
本文以云霄县第三次全国国土调查为例,介绍利用ArcGIS对耕地坡度赋值的方法,极大提高工作效率。
1.1 三调地类图斑 地类图斑属性结构表见表1。
表1 三调地类地类图斑属性结构表(部分)
6 图斑地类面积 TBDLMJ Float 15 M 7 耕地坡度级别 GDPDJB Char 2 C………………
1.2 坡度数据 利用下发的云霄县坡度数据作为耕地坡度赋值基础数据。
坡度数据属性结构表见表2。
表2 坡度数据属性结构表
坡度赋值应当遵循一定的原则,根据第二次全国土地调查中《利用DEM确定耕地坡度分级技术规定(试行)》要求:
(1)原则上不能打破调查的耕地图斑界线,每个耕地图斑确定一个坡度级;
(2)当调查的耕地图斑涉及两个以上坡度级时,面积最大的坡度级为该耕地图斑的坡度级。
(3)当耕地图斑面积较大(如从山顶到山底为一个图斑)、含有两个以上坡度级时,且各坡度级耕地面积相当时,可参照坡度分级界线,依据调查底图(DOM)上明显地物界线,可将该耕地图斑划分为两个以上不同坡度级的图斑。
(4)对于破碎耕地,其整体视为一个图斑,按上述要求确定坡度分级。
(5)当DEM存在缺陷时,应通过其他手段补充、完善,确定耕地图斑坡度级。
由于云霄县(3)(5)情况较少,故本次赋值的技术方法只针对(1)(2)(4)适用。
利用ArcGIS对耕地坡度赋值的技术流程:
(1)对云霄县三调数据库的地类图斑标识码进行赋值,使每一个地类图斑具有唯一的标识码,便于后续坡度数据准确挂接。
(2)在三调数据库DLTB属性表中按属性[DLBM]<='0103'选择耕地图斑,共计9636个。
(3)利用ArcGIS工具箱中的分析工具将耕地数据和坡度数据做相交处理,保留两个数据的所有属性连接。
(4)由于一个耕地图斑内可能存在多个不同级别的坡度数据,也可能一个耕地图斑内具有多个相同坡度级别的数据分布在不同的区域;在这种情况下将耕地数据和坡度数据做相交处理会出现同一编码的耕地图斑对应多个坡度级别,且每个坡度级别对应的图斑面积是独立的,即使是同一个坡度级别的数据也未进行汇总,因此也无法判断哪个坡度级别占比是最大的。
(5)利用ArcGIS工具箱分析工具对相交的数据进行汇总统计数据,统计字段为SHAPE_Area,统计类型为SUM,案例分组字段为BSM、PDJB。
(6)为了快速获取同一标识码内面积占比最大的坡度级别,首先需要将上一步获取数据进行高级排序,标识码升序,SUM_SHAPE_Area为次排序降序。再利用ArcGIS工具箱分析工具对上一步的数据进行汇总统计。统计字段SUM_SHAPE_Area,统计类型 MAX;统计字段PDJB,统计类型FIRST;案例分组字段BSM。
(7)在三调数据库地类图斑要素中利用标识码挂接(6)汇总统计的坡度数据,再利用字段计算器将PDJB数据赋值给GDPDJB,最终完成三调耕地坡度赋值。
本文介绍了利用ArcGIS各项工具完成三调耕地坡度赋值的方法。该方法只需采用ArcGIS一种软件便可完成,无需编写程序,操作简单,可一次性完成整个县域的坡度赋值,极大提高了内页数据处理的效率和准确性,为云霄县三调工作助力提速。