地理信息技术在保障性安居工程审计中的应用

项程光一 万卓敏 杜辰华

近年来，随着计算机技术的发展以及信息大数据的广泛应用，审计中运用大数据审计的重要作用日益凸显。在保障性安居工程审计中，我们引入了地理信息技术进行大数据审计，经过反复运用和验证，开发了地理信息大数据审计方法。采用地理信息技术进行审计，不仅可以有效的盘清“库存”，还可以发现保障性安居工程建设项目位置偏远、土地闲置等问题，充分发挥了大数据审计全面性、精准性、高效性以及直观可视化等优势。现结合实际，介绍如何运用地理信息技术进行大数据审计。

一、源数据的采集

地理坐标源数据采集需要统一坐标系统及格式。由于地理坐标信息格式种类繁多，每一种坐标系中又可以通过一定的投影方式形成新的平面坐标系统，不同坐标系之间的坐标点都会存在或大或小的偏移。为更加方便快捷的采集地理坐標信息数据，保障性安居工程审计按照大家熟悉的高德地图经纬度坐标采集，然后通过算法精确转换成“2000国家大地坐标系”做数据分析与处理，同时可以方便的将采集到的地理信息在谷歌地球软件中展现，并查看对应历史影像。

保障性安居工程审计中涉及到的地理坐标信息主要有两大类，一类是以项目为单位的坐标点数据，另一类是以地块为单位的坐标区域数据，这两类坐标数据的处理以及导入方式各不相同，故分别拟定以下标准表：

1.坐标点数据标准表：

注：上表A、B、C列为关键字，下发此表时可将审计组取得的保障房项目名称及地址汇总后一同提供给被审计单位对照填写，以防出现项目名称不一致的情况。此表其余列可根据数据分析内容及审计需求自行拟定，或后期在SQL中根据C列关键字直接关联分析。

2.坐标区域标准表：

坐标区域标准表应按如下格式填写。将地块界址点坐标的经、纬度按标准格式拆分，分别填列至A、B列，同一地块的界址点对应的C列单元格均填充相同地块名称。

注：地块数据由土地管理局地块登记系统后台导出后填列。

二、坐标数据处理及导入

1.点坐标

将采集的坐标点数据汇总后直接导入KML制作工具“kmlgenerator”中生成KML文件，或直接将汇总的EXCEL表格另存为“CSV（逗号分隔符）”格式，以便下一步导入QGIS。

（直接生成KML文件）

（EXCEL另存为CSV格式）

2.地块坐标

（1）将采集的地块坐标数据复制到“界址坐标KML文件生成工具”的“大地坐标”表中（为避免坐标数据意外损失，请选择“黏贴值”），若需要添加字段可自行随后添加列。

（界址坐标KML文件生成工具）

（2）对于需要在QGIS中进行数学计算的字段，在“字段类型”表中设置相应类型。切换到“字段类型”表，点开第二行下拉框，选择“float”。“字段类型”表无需人工填入数据，当在“大地坐标”表中添加或删除列后会自动在此表中更新相应字段及类型，默认类型为“string”，仅当需要更改字段类型时才手动选择。

（选择字段类型）

注：当“大地坐标”表列标题有变动时，“字段类型”表中第二行数据会全部重置为string。故请先确定好“大地坐标”表的列标题，然后再在此表中设置字段类型。

（3）点击“加载项”中的“生成地块坐标KML文件”按钮，程序会自动生成地块的坐标数据并存入本表的存放路径，文件名为“地块坐标.kml”。

（按钮位置）

（生成结果）

3.坐标导入

将制作好的点坐标及地块坐标KML文件导入QGIS软件中进行进一步分析。

（1）首先在QGIS中新建一个项目。依次点击“项目”—“新建”。

（新建项目）

（2）若需要查看坐标点实际位置，还需要添加在线地图。在资源浏览器中找到“XYZ Tiles”，展开节点，双击“OpenStreetMap”，将地图添加进项目，此在线地图为WSG84坐标系，即谷歌地图的坐标系，与“2000国家大地坐标系”通用。

（3）若要查看影像地图，可手动添加地图源。

（右击XYZ Tiles，新建连接）

（命名地图名称自定）

（添加后即可在资源浏览器中找到）

（4）根据影像地图可直接查看对应项目坐标点或地块的实际情况。但是由于中国地区的谷歌影像数据更新较慢，目前很多区域仍是数年前的影像，故如必要，可先将坐标点及地块坐标数据使用“GPS批量查询转换工具”，将“2000国家大地坐标系”统一转换为高德坐标系，再导入QGIS中，然后使用高德影像地图查看项目实际情况。

（5）坐标数据的kml文件导入后，最好再将导入的图层另存为shp格式，然后移除原kml图层，以免后期进行数据操作时出错。

（图层另存为）

（格式选择ESRI Shapfile，设置shp文件存储路径，坐标参照系一般默认）

（勾选“将已保存的文件添加到地图中“，勾选欲保存的字段。）

（保存后在图层管理器中会看到两个相同图层，新添加的那个即是shp图层，将原有kml图层移除。）

三、数据分析及结果

1.筛选出地处偏远的项目

（1）在QGIS中分别测量城市规划区经度和纬度方向的跨度，取其平均值作为该区域的直径R。

（2）将汇总的项目坐标信息复制到“GPS批量查询转换工具”，查询项目点至政府直线距离及所在行政区划，同时核实项目点行政区划是否有误，保存结果，形成“表D01”。在表D01中筛选出距政府直线距离大于（3/4）R（具体距离可根据实际情况自行调整）的项目。

（3）在QGIS中筛选出未落在规划区内的项目点，并将筛选结果导出。

（依次点击菜单栏的设置——选项，在弹出的界面中选择最下面一个标签  ，点开“通用”节点，双击“要素过滤器无效”右侧，将选项改为“不过滤”，设置完成后点“OK”。）

（依次点击矢量——地学数据处理工具——差分。）

（按图示选择项目坐标点及城市规划区图层）

（点击“差分”栏右边按钮，选择比对结果保存路径及格式，一般选择xlsx，取消勾选“执行算法后打开输出文件”，点击下方“运行”即可生成结果。结果表中的项目即为未在规划区内的项目）

（4）将（2）、（3）中两部分筛选结果数据汇总，即为“偏远项目表”。

（5）针对筛选出的偏远项目重点核实，是否存在定向分配、违规将棚改政策扩大到一般建制镇等问题。

2.筛选出空置地块

（1）在QGIS中筛选保障房划拨地块内无项目坐标点的地块。

（依次点击矢量——分析工具——计算多边形内点的数目。）

（设置需要比对的地块及点坐标图层）

（在“计数”栏选择比对结果保存路径及格式，点击下方“运行”即可生成结果。）

（2）将分析结果导出，导出表中“NUMPOINTS”列即为对应地块中的坐标点个数，值为“0”的即为空置地块，将此部分地块筛选整理出来，形成“闲置地块表”。

（闲置地块表）

（3）在“闲置地块表”中筛选出距今超过三年的地块，针对此部分已划拨但闲置超过三年的地块重点核查原因。

通过地理信息技术应用于保障性安居工程审计，突显先进技术手段给审计分析过程中提供大量化、多样化和快速化等特点，不仅大量节省了审计成本，还将原来一些难以认定的情形实现具体可操作化，提高了审计质量和效率。在实际审计中，审计人员应同时结合传统审计方式，借助地理信息技术分析出的疑点，进一步落实保障性安居工程项目选址不科学，建设位置偏远导致空置、违规扩大棚改范围以及保障性安居工程地块闲置等问题，巩固完善审计证据链。

（作者单位：景德镇市审计局、江西省审计厅）