基本农田变更方法设计与实现DOC

1、基本农田变更方法设计与实现邓光林1，卢建青2，彭涛2 （1.义乌市国土资源局,义乌322000 2.浙江省测绘大队，杭州 310030）Method of Design and Realization of the Basic Farmland ChangeDENG Guanglin,LU Jianqing, PENG Tao 摘要：基本农田保护片块年度变更是土地利用年度变更的重要工作之一.本文通过综合应用ArcGIS中GP工具，设计了一整套基本农田变更方法和关键步骤的实现算法，利用AO开发实现了变更方法软件化,并成功的应用在义乌市的基本农田保护图斑数据更新工作中。关键词：基本农田；变更；图斑

2、；线状地物；净面积计算；AO；GP工具；插件一、引言基本农田保护片块年度变更的主要内容是对基本农田保护图斑（JBNTBHTB）进行净面积计算，将JBNTBHTB图斑按照同一坐落（ZLDWDM）融合成基本农田保护片块（JBNTBHPK），对JBNTBHPK层中片块按照ZLDWDM进行编号，JBNTBHTB层中图斑按照所在PK号进行编号。从JBNTBHTB层中统计生成基本农田情况统计表、基本农田调查汇总表（二级类）、基本农田调查汇总表等。JBNTBHTB基础数据获取主要有两种方式，一种是从本年度地类图斑层中提取地类类型属于耕地的地类图斑，再根据规划红线进行分割得到JBNTBHTB基础数据，一种是利

3、用上一年度的JBNTBHTB图层扣除本年度新增的建设用地要素得到JBNTBHTB基础数据。基本农田保护片块年度变更数据处理的思路简单，过程复杂，方法多样，但基于不同的GIS 软件平台，基本农田保护片块年度变更的实现途径又有差异。浙江省大部分县市在全国第二次土地调查及土地数据建库采用的都是ArcGIS 平台，本文针对该平台对基本农田保护片块年度变更方法进行了分析与探讨，设计了净面积计算、椭球面积计算、片块编号等关键技术实现方法，利用AO开发了基本农田变更软件，该软件在义乌市的基本农田保护片块图斑数据更新工作中得到了应用。二、处理流程整个工作流程图如图一所示：提取JBNTBHTB基础数据利用本年度

4、DLTB、XZDW、LXDW对JBNTBHTB进行净面积计算生成JBNTBHPKJBNTBHPK编号,JBNTBHTB编号统计结果出表图1 计算流程图基本农田年度变更的整个流程如图一所示：数据准备（面状DLTB，线状XZDW,点状LXDW, JBNTBHTB）第一步：提取JBNTBHTB(基本农田保护图斑)本年度地类图斑层中提取地类类型属于耕地的地类图斑，再根据规划红线进行分割得到JBNTBHTB基础数据，或利用上一年度的JBNTBHTB图层扣除本年度新增的建设用地要素得到JBNTBHTB基础数据。第二步：JBNTBHTB净面积计算JBNTBHTB净面积重新计算，包括提取与该图斑相关联的线状地

5、物、零星地物重新提取XZDWMJ和LXDWMJ，再计算图斑地类面积整个计算过程如下：利用JBNTBHTB层分割XZDW、LXDW层得到JBNTBHTB范围内的XZDW、LXDW。处理分割后的XZDW的KCBL(扣除比例)。 利用空间分析建立分割后的XZDW与JBNTBHTB的对应关系计算JBNTBHTB中XZDWMJ、TKMJ、LXDWMJ.第三步：生成JBNTBHPK从JBNTBHTB中按照ZLDWDM(坐落单位代码),将统一ZLDWDM空间相连的JBNTBHTB图斑进行合并生成JBNTBHPK。第四步：JBNTBHPK编号、JBNTBHTB编号JBNTBHTB编号以坐落在同一单位内的图斑为

6、组合，按照从上到下，从左到右编码。JBNTBHPK编号以坐落在同一单位内的片块为组合，按照从上到下，从左到右编码。第五步：统计结果出表JBNTBHTB中净面积计算好后，生成基本农田情况统计表、第二次全国土地调查基本农田调查汇总表（二级类）、第二次全国土地调查基本农田调查汇总表等表。三、关键技术处理1. 净面积计算基本农田保护图斑是在地类图斑的基础上提取水田、园地、林地等相关地类构成基本农田保护图斑，其中保护图斑范围按照相关的规划红线范围进行分割。因此需要对提取后的保护图斑重新计算净面积。 1.1净面积计算公式地类图斑净面积的计算公式如下:图斑净面积= 图斑毛面积零星地物面积线状地物面积田坎面积

7、其中:线状地物面积=长度×宽度田坎面积= ( 图斑毛面积零星地物面积线状地物面积)×田坎系数如果线状地物同时也是地类图斑的边界则半扣，即:图斑净面积=图斑毛面积零星地物面积线状地物面积×0.5田坎面积若线状地物位于地类图斑之内则全部扣除。根据上述思路图斑净面积计算的难点是确定图斑所关联的线状地物的扣除系数。1.2线状地物扣除系数处理扣除系数的确定原则是: 线状地物位于地类图斑之内，则扣除系数为1.0，若位于地类图斑边界，则为0.5，可以根据线状地物与地类边界的空间关系来判断。但不适用于特殊情况，位于地类边界但扣在一个图斑内，扣除比例为1。线状地物扣除系数处理的基本

8、思路如下：（1）重构线状地物与基本农田保护图斑拓扑结构，拓扑结构是指线状地物与地类图斑的空间关系，线状地物跨越地类图斑边界时，要断开。利用ArcGIS中Identify工具生成线状地物与保护图斑叠加层，层中线跨面边界自动断开，同时记录了线层中左右图斑面信息。通过从层中提取左右面中信息中包含保护图斑信息的线要素。构建保护图斑与线状地物的拓扑结构。（2）重新计算提取后的线状地物的长度、线状地物面积。（3）重新计算提取后的线状地物的扣除比例：a.先将提取后线状地物层中左右扣除图斑信息为同一图斑的设置扣除系数为1。其他扣除系数为0.5的保留。b.原扣除系数为1，但提取后线状地物层中左右扣除图斑信息不一

9、致处理。这部分系数确定不能直接按照扣除系数系统确定的原则来确定，因为如下图所示：图2 线包含在图斑内分割图图3 线与两图斑相交分割图如图2和图3所示，原线状地物扣除比例都为1,图2原线状地物包含在图斑内，图3原来线与两图斑相交。图2处理只需要通过空间分析判断原线是否包含在图斑内，若包含在图斑内则扣除系数改为0.5。图3处理确定分割后扣除系数扣在哪一个图斑上，首先通过空间分析获取分割后的面状图斑，利用该图斑与分割前图斑做空间比较，根据以下条件：条件1：分割后图斑包含在分割前图斑中条件2：分割前扣除的图斑信息是否与分割后图斑一致。通过上述条件即可确定扣除系数扣在哪一个图斑上。2片块编号基本农田保护

10、图斑层片块编号的基本思路是取坐落在同一个行政村范围内的图斑，获取每个图斑的几何中心点坐标构造数组，图斑按照从左到右，从上到下的原则进行排序。确定每个图斑的排列序号。该算法充分利用ArcGIS面状图形接口IArea获取中心点IPoint，根据中心点X，Y坐标进行排序。图4 片块编号算法四、功能实现功能实现了自动计算净面积、自动计算椭球面积、自动片块编号、自动生成基本农田保护片块、根据自定义报表格式自动输出成果报表。并对计算结果自动核对。整个软件功能在ArcObject的基础上开发，通过设计功能插件，将功能以插件命令的方式嵌入到ArcMAP界面中。功能结构图如图5所示。生成报表JNBTBHTB编号

11、赋值生成JBNTBHPK计算图斑面积图5 功能结构图 （1）计算图斑面积功能：通过导入地类图斑、线状地物、零星地物、基本农田保护图斑。重新计算保护图斑椭球体面积，提取保护图斑范围内线状地物并构造图斑与线状地物拓扑结构，重新计算保护图斑净面积。（2）生成基本农田保护片块层：根据基本农田保护图斑层，按照坐落在同一行政区内且空间上相连接的面状图斑合并为一个片块。生成基本农田保护片块层并计算片块编号。（3）基本农田保护图斑编号: 根据基本农田保护图斑层，按照坐落在同一行政区内的图斑自动编号。系统操作界面如下图所示:图6 系统界面五、结束语本文利用ArcGIS软件数据处理工具设计了一整套基本农田数据库变

12、更的工作流程，通过利用ArcGIS Object开发实现了操作过程自动化。软件解决了地球椭球面积自动计算、片块编号等关键技术。尤其是对于线与两图斑相交但扣除系数为1的关键技术的处理解决了以往人工核对效率慢的问题；同时软件在数据处理过程中对图斑、片块数据的拓扑及其属性逻辑错误自动检查。保证了数据处理的效率以及数据处理的正确性。参考文献1徐志红.谈新形势下的土地变更调查J.浙江国土资源,2011(4):36-382邹建胜,屠海东,金卓妮.县级年度农村土地变更调查实践J.浙江测绘,2011(3):44-453史大成,温瑞智,任叶飞,周宝峰.基于GIS的场地分类方法研究J.地理信息世界,2011(1):23-27.4兰小机,刘德儿,魏瑞娟. 基于ArcObjects与C#.NET的GIS应用开发M.北京,冶金工业出版社,20115刘洪江,曹玉香.基于ArcGIS 实现地类图斑净面积的计算J.城市勘测,2012(5):115-1176潘雅辉,冯杭