ArcGIS在土地报批中的应用*

2023-01-18于凤兵谢克文祖鹏蕾

地矿测绘 2022年4期

于凤兵，谢克文，祖鹏蕾

(几核(浙江)科技有限公司，浙江湖州 313000)

0 引言

土地报批是自然资源管理的日常工作之一。ArcGIS的功能较强，可以对各种空间数据进行分析和处理[1]，能为土地报批工作提供技术支持，大大地提高数据统计、叠加分析等工作的效率[2]。

土地报批中土地勘测定界、功能区划分及指标计算、耕地占补平衡、永久基本农田占补平衡、违法用地查处等多个流程都需要为自然资源行政主管部门提供科学准确的数据资料[3]，其中单独选址项目更是工作量大，时间紧。ArcGIS的应用，可以有效避免重复工作，保证数据的准确性。

1 ArcGIS在土地报批中的应用

1.1 在土地勘测定界中的应用

土地勘测定界，是根据土地征收、征用、划拨、出让、农转用等工作需要，实地界定土地使用范围、测定界址位置、调绘土地利用现状，计算用地面积，为自然资源行政主管部门用地审批等提供科学、准确的基础资料而进行的技术性服务，是土地报批的基础工作。

考虑到用地红线多为设计单位提供，设计软件基于CAD，其注重线形以及图形质量，而实体间的拓扑关系则不太注重[4]，运用ArcGIS的拓扑检查功能对用地红线进行拓扑检查，拓扑规则为不能重叠和不能有空隙，如图1所示。

图1 拓扑规则设置图Fig.1 Topology rule setup diagram

土地勘测定界主要是地类分类面积统计和土地勘测定界图的制作。主要内容都是计算土地利用现状图的地类图斑面积，地类图斑是指被单一线状地物以及权属和行政界线分割的单一地类地块，计算公式如下：

图斑面积=图斑毛面积-田坎面积-零星地物面积-线状地物面积

其中，线状地物面积=线状地物长度×线状地物宽度×扣除比例，扣除比例由线状地物与图斑的关系来确定，如果线状地物与地类图斑重合，扣除比例为0.5，若线状地物全部位于图斑内，则扣除比例为1[5]。计算公式原理看似简单，但由于所有数据均涉及到空间关系的查询和判别，加上图斑量大等因素，人工计算面积工作量大，效率低，容易出错。基于ArcGIS二次开发编写的小插件，实现点、线、面图层分别与用地红线的叠加分析，获取图斑面积，结果如图2所示。

图2 图斑面积自动获取数据成果图Fig.2 Automatic acquisition of map spot area

土地勘测定界图需要在图中按从上到下、从左到右的顺序标注出每一个图斑的地块号、地类代码及地块面积，在CAD或CASS中标注工作量大，地块总面积汇总极易出错，结合ArcGIS的排序、标注等命令可以快速、准确地完成土地勘测定界图的图斑标注。

地块号结合ArcGIS的排序命令实现，排序方法选择“ASCENDING”，按照值从低到高的顺序对记录进行排序，空间排序方法按图形选择UR、UL、LR、LL、PEANO进行选择，如图3所示。

图3 图斑地块号排序图Fig.3 Map plot number ordination

添加字段“地块号”，类型为文本型，结合Python顺序编码插件进行赋值。

为了实现地块号位数一致，需要进行补零，然后字段计算器执行命令：right("000"& [顺序号]，X)，若图斑总数为两位，X为2，图斑总数为三位，X为3。

地块号、地类编码和面积完成赋值后，新建字段“标注”，类型为文本型，字段计算器执行命令：[顺序号]&vbnewline&"——"& [MJ] &"m2"&vbnewline& [DLBM]，完成图斑标注。另外。也可以采用VBScript解析程序完成标注，如图4所示。

图4 标注表达式图Fig.4 An expression

具体表达式如下：

FUNCTION strlen(str)

dim p＿len

p＿len=0

strlen=0

p＿len=len(str)

FOR xx=1 to p＿len

IF asc(mid(str，xx，1))<0 then

strlen=int(strlen)+2

ELSE

strlen=int(strlen)+1

END if

END function

FUNCTION myFind(cunname，DJH，SHAPE＿Area)

dim str

str=SHAPE＿Area

肠道菌群与呼吸道抗感染能力：肠道菌群能够通过TLRs促进相应的免疫细胞的活动，促进IFN-α、ENF-α等炎性细胞因子的分泌，调节呼吸道严重反应和对病毒的清除反应。

dim d

d=strlen(str)

dim d1

dim d2

d1=strlen(cunname)/2

if d1<1 then

d1=1

end if

d2=strlen(DJH)/2

d2=1

end if

if d2>d1 then

d1=d2

end if

myFind=cunname & space(d)& vbnewline & string(d1，"—")& str &""& vbnewline & DJH & space(d)

END Function

FunctionFindLabel ([DKH]，[DLBM]，[SHAPE＿Area])

FindLabel=

myFind([DKH]，[DLBM]，[SHAPE＿Area])

End Function

将标注转换为注记、要素转CAD后完成土地勘测定界图图斑标注。

1.2 在功能区划分及指标计算中的应用

项目用地面积和功能分区需符合建设用地标准控制等节约集约用地要求，符合工程初步设计批复要求。建设用地审批必须严格执行用地标准，对超出用地标准的，要核减用地面积。

功能区划分及指标计算即依据设计提供的用地图进行功能区划分，然后按照建设用地指标进行复核。以新建客运专线铁路为例，功能区分为：路基、桥梁、隧道、牵引变电所、分区所、车辆基地等。使用CAD分割图形工作量大，易产生拓扑图斑，将CAD格式的用地红线转为SHP文件，应用ArcGIS的裁剪面工具和追踪工具，可以快速准确地实现功能区划分。功能区划分完成后，可以运用ArcGIS的拓扑检查功能对用地红线进行拓扑检查，避免功能分区间存在重叠或缝隙，确保用地总面积不变。新建字段“功能分区”，依据设计用地图对功能分区进行赋值。其中，桥梁、路基等用地规模计算时需以长度为基础，可以结合ArcGIS的提取中心线工具获取中心线。

先把用地红线进行融合，将融合后的数据转成线数据，然后对线数据提取中心线，为了防止只提取出外围轮廓线，需将转换后闭合的线数据进行打断处理，在道路任意一个末端使用打断工具进行打断处理，后使用ArcGIS自带的提取中心线工具提取道路中心线。中心线提取成功就可获得整条道路的长度，将道路中心线与用地红线的功能区进行相交，即得到每个功能分区的用地长度。路基、桥梁等功能分区就可通过指标数*长度=用地规模来复核建设用地的用地面积是否符合建设用地规模。

对国家和地方尚未颁布土地使用标准和建设标准的建设项目，以及确需突破土地使用标准确定的规模和功能分区的建设项目，应当开展节地评价。

1.3 在耕地占补平衡中的应用

耕地占补平衡是指建设占用多少耕地，各地人民政府就应补充划入多少数量和质量相当的耕地行为。耕地占补平衡在数量质量上严格立足“占一补一、占优补优，占水田补水田，确保耕地占补平衡数量质量双到位”。

首先要获取耕地数量，定义查询DLTB为水田和旱地，与用地红线进行叠加分析，提取用地红线范围内耕地图斑。耕地质量等级是耕地图斑与耕地质量等别年度更新评价成果中的“县级分等单元”图层进行叠加，提取图斑内耕地质量的“国家利用等”属性信息。考虑到耕地质量等别年度更新评价成果与耕地图斑不一定完全重叠，应使用ArcGIS的标识功能，将耕地质量等别字段标识到耕地图斑属性表中，这适用于耕地图斑与耕地质量分等单元有重叠部分；与耕地质量分等单元不重叠的耕地图斑，需参照邻近周边的耕地质量等别信息逐个进行赋值。获得每个耕地图斑的面积和耕地质量等别信息后，使用加权平均的方法来获取耕地的平均质量等级。

补充耕地数量、水田规模和粮食产能均需符合要求，做到耕地占补平衡数量、质量双到位。

1.4 在永久基本农田占补平衡中的应用

重大建设项目占用永久基本农田的，按照“数量不减、质量不降、布局稳定”的要求进行补划。利用ArcGIS将用地红线与永久基本农田数据进行叠加分析，得到项目占用的永久基本农田，按照永久基本农田数据的TBMJ、JBNTMJ来计算占用永久基本农田的净面积。添加字段“净面积”，用字段计算器的计算几何得到图形面积，用净面积计算公式得到占用永久基本农田的净面积，计算公式如下：

占用永久基本农田面积=图形面积/TBMJ*JBNTMJ

依据永久基本农田每块图斑的耕地质量等级，用加权平均的方式获取平均耕地质量等级。按照永久基本农田补划原则“占用城市周边永久基本农田的，原则上在城市周边范围内补划”。

需在长期稳定利用的耕地上落实永久基本农田补划任务，补划地块选取时，按照补划地块在第三次国土调查中为优质耕地(稳定耕地)和实地为耕地的补划原则进行补划地块筛选，结合与原永久基本农田连片，耕地质量等级和城镇周边等因素选取出符合条件的补划地块，做到永久基本农田数量不减少，质量不降低。

永久基本农田示范区补划中补划地块需满足补划数量不减少、质量不降低、布局基本稳定，确保城镇周边永久基本农田示范区占补平衡，达到了永久基本农田(示范区)连片度不降低的要求。

1.5 在违法用地查处的应用

违法用地是依托主管部门提取疑似动工图斑，结合自然资源部下发的即时卫片，综合利用时序遥感影像数据，通过人工识别地表目标的空间分布状况、变化及其变化情况，分析项目用地内疑似动工迹象的位置、范围以及其他变化情况等信息。

利用ArcGIS将下发的疑似动工图斑与农转、供地、临时用地等审批资料进行叠加分析，结合现场实地踏勘，对疑似动工图斑逐个分析，明确是否存在违法用地情况，如涉及违法用地，应明确违法用地的动工时间，查处确认情况；将疑似动工图斑与生态保护红线进行叠加分析，涉及生态保护红线范围内违法用地的，应说明主观故意情况。