平陆县审计局：如何利用ArcGis地理信息系统软件 进行审计数据分析

随着疫情、洪涝灾害等问题不断爆发，粮食安全问题再次引起全国各地高度重视。因土地利用规划及国土资源涉及面广、数据量大，加之国土资源本身的空间属性，使得在开展审计项目时，必须借助信息化手段，运用ArcGis地理信息系统软件、SQLServer数据库、BIGEMAP地图等软件，来发现审计疑点线索、核实确定问题，下面笔者以基本农田“非粮化”及违规发放耕地地力保护补贴为例，简单介绍如何利用ArcGis地理信息系统软件进行审计数据分析。

一、数据采集

ARCGIS数据分析包括矢量数据分析、权属数据分析、栅格数据分析、 络数据分析等。在这里我们经常使用的就是对矢量数据中的.SHP格式的数据进行分析，所以在采集数据时数据格式为*.SHP。

1. 土地利用现状分类

2017年由国家质量监督检验检疫总局和中国国家标准化管理委员会联合发布的《土地利用现状分类》，如：‘01’开头地类编码代表‘耕地’、‘02’开头代表‘园地’、‘03’代表‘林地’、‘0810’公园绿地、‘1104’坑塘水面（养殖坑塘）等。

1. 案例分析

1–基本农田‘非粮化’

根据《基本农田保护条例》第十七条“禁止任何单位和个人占用基本农田发展林果业和挖塘养鱼”的规定，永久基本农田要实行特殊保护，防止“非粮化”。在这里我们要根据规定将【2020年国土变更调查数据库】（三调）中提取林地、园地、坑塘水面与永久基本农田图斑进行叠加相交分析，汇总计算出基本农田现状为园地、林地、坑塘水面的面积。

步骤1.打开ArcGis软件，在右侧目录栏中进行文件夹连接，将采集到的基本农田数据和2020年调整数据连接至ArcGis。

步骤2.点击2020年调整数据，选中DLTB（地理图斑）拖至左侧图层栏目下，点击基本农田数据，选中JBNTBHTB.SHP拖至左侧图层栏目下。

步骤3.右击DLTB，点击属性-定义查询-查询构建器-在下面的对话框中写入筛选条件，方法一：DLBM like ‘02%’ or DLBM like ‘03%’ or DLBM =’1104’；方法二：DLMC in( ‘果园’ , ‘坑塘水面’, ‘其他园地’, ‘竹林地’ ,’乔木林地’, ‘其他林地’, ‘灌木林地’ ,’灌丛沼泽’ )，标点符 均在英文输入法状态下输入，筛选出2020年调整数据中的林地、园地、坑塘水面数据。

步骤4：右击JBNTBHTB，点击属性-定义查询-查询构建器-在下面的对话框中写入筛选条件，DLBZ not like ‘K’，剔除基本农田中DLBZ（地类标志）为K（可调整园地）的面积，此处需要注意，K需要复制原数据的字段，不能直接输入K因为可能会存在无法识别的情况，不能完全剔除。

步骤5：点击菜单栏中地理处理-相交，得出林、园、渔业占用基本农田图斑，右击相交图斑-打开属性列表，左击表选项选择“添加字段”并给字段起名（MJ），类型选择“双精度”，右击MJ字段选择‘计算几何’，单位选择“平方米”。

步骤6：在数据分析时，大家应该对图斑的表格要素进行了解，我们可以发现永久基本农田的面积= 图斑面积扣除田坎面积扣除线状地物面积扣除零星地物面积，在重新计算以后得到的是图斑面积，不是基本农田的实际面积，所以我们还要扣除田坎、线状地物、零星地物的面积，对上一步的结果进行进一步的细化和精确。因为GIS里面进行复杂的计算容易自动关闭，我们在实际操作中可以将其导出SQL_SERVER进行计算处理，得出扣除后的结果,再通过地类名称进行分类汇总。

步骤7：结合BIGEMAP将【结果图斑】导入，通过定位功能查看现场，拍照取证。

案例分析2-应扣未扣补贴面积，违规发放耕地地力保护补贴

根据XX省财政厅、XX省农业厅《XX省财政厅XX省农业厅关于做好农业支持保护补贴资金工作的通知》（X财农〔20XX〕XX ）第一条“补贴范围原则上为种地农民已确权登记颁证到户的拥有承包权的耕地。以下情形不在补贴范围：1、已作为畜牧养殖场使用的耕地……”和第六条“第四项严格监督检查，各级农业、财政部门应加强政策落实的日常工作监管，落实定期检查指导和重点抽查相结合，严格核实补贴对象和面积，严禁对不符合补贴条件的耕地给予补贴……”的规定，运用ArcGis叠加、计算几何、汇总功能、SQL-SERVER及EXCEL中VLOOKUP函数实现。（‘六不补’主要为‘畜牧养殖场用地面积’‘颁发林权证或退耕还林土地面积’‘转为设施农业用地面积’‘非农业征、占用土地’‘常年抛（撂）荒耕地’‘质量达不到耕种条件土地’）此处以‘设施农用地’为例。

1.将【农村承包经营权】地理图斑连接至ArcGis，找到以DK（地块）开头后缀为*.shp的格式文件，拖拽至左侧图层；

2.右击DK.shp，点击属性-查询构建器，剔除非承包地等内容（如村组机动地）这里，以便与农户申请补贴确权面积一致。（‘DKLB’表示‘地块类别’，‘10’代表‘承包地’；

3.将【2020年度国土变更调查数据库（三调）】连接至ArcGis，找到矢量图层中DLTB.shp拖拽至左侧图层，右击属性，查询出‘设施农用地’数据；

4.通过ArcTOOLBOX工具—叠加分析-相交，将【2020年度国土变更调查数据库（三调）】与DK*.shp相交，相交完成后，打开添加字段-计算几何得出设施农用地确权面积。

5.右击相交后的图层‘打开属性表’，查看所有字段，发现缺少乡镇名、村名，只有农户姓名，与被审计单位填 的耕地地力补贴发放调查表无法比对。对权属单位代码表分析后发现，在DK**********.shp中有DKBM（地块编码）字段，DKBM如：123456789asdfghjklqw ，前123456为县级编 、789为乡级编 、asd为村级编 ，fghjklqw为地块编 ，通过SQL-SERVER来连接后获取乡镇名、村名。

6.将ArcGis与SQL-SERVER连接，打开ARCGIS目录点开数据库连接添加数据库连接，实例为计算机名称，身份验证类型选择数据库身份验证，用户名SA 密码为自己密码，确定。右击连接的库导入表选择即可，通过SQL语句处理以后，得到想要的结果。

7.以村为单位将设施农用地确权面积与农业部门提供的耕地地力发放表中已经扣除的‘六不补’面积（特指“设施农用地”面积）进行比对，比对前先对 送的表格用数据透视表进按村进行汇总，利用VLOOKUP函数，查看是否存在应扣除未扣除补贴面积的问题。

8.审计人员利用BIGEMAP地图将问题疑点图斑转为*.KML格式，导入Google Earth Pro中快速确定设施农用地方位，通过地图影像查看分析设施种类，并结合山西农业补贴资金管理平台确定是否对该处设施农用地面积进行了扣除，最后抽调问题疑点图斑面积较大的利用Google Earth Pro定位进行实地核查，现场取证，确定问题，计算出违规金额。

1. 案例应用延伸

1.绿化造林目标任务完成情况；

2.耕地、永久基本农田保护目标任务完成情况；

3.是否存在入河排污口设立在生态红线保护区内的情况；

4.是否存在在畜禽养殖禁养区内开展畜禽养殖的情况；

5.在限制或禁止开发区域违规审批建设用地情况；

6.居住地不属于施教区范围，违规纳入入学情况；

7.收集各单位（如水利项目、土地整理项目等）土地项

目矢量数据，相互对比发现是否存在在同一地块重复实施，浪费财政资金或骗取套取财政资金的行为等。

平陆县审计局

刘莹