項(xiàng)程光一 萬卓敏 杜辰華
近年來,隨著計(jì)算機(jī)技術(shù)的發(fā)展以及信息大數(shù)據(jù)的廣泛應(yīng)用,審計(jì)中運(yùn)用大數(shù)據(jù)審計(jì)的重要作用日益凸顯。在保障性安居工程審計(jì)中,我們引入了地理信息技術(shù)進(jìn)行大數(shù)據(jù)審計(jì),經(jīng)過反復(fù)運(yùn)用和驗(yàn)證,開發(fā)了地理信息大數(shù)據(jù)審計(jì)方法。采用地理信息技術(shù)進(jìn)行審計(jì),不僅可以有效的盤清“庫存”,還可以發(fā)現(xiàn)保障性安居工程建設(shè)項(xiàng)目位置偏遠(yuǎn)、土地閑置等問題,充分發(fā)揮了大數(shù)據(jù)審計(jì)全面性、精準(zhǔn)性、高效性以及直觀可視化等優(yōu)勢?,F(xiàn)結(jié)合實(shí)際,介紹如何運(yùn)用地理信息技術(shù)進(jìn)行大數(shù)據(jù)審計(jì)。
一、源數(shù)據(jù)的采集
地理坐標(biāo)源數(shù)據(jù)采集需要統(tǒng)一坐標(biāo)系統(tǒng)及格式。由于地理坐標(biāo)信息格式種類繁多,每一種坐標(biāo)系中又可以通過一定的投影方式形成新的平面坐標(biāo)系統(tǒng),不同坐標(biāo)系之間的坐標(biāo)點(diǎn)都會存在或大或小的偏移。為更加方便快捷的采集地理坐標(biāo)信息數(shù)據(jù),保障性安居工程審計(jì)按照大家熟悉的高德地圖經(jīng)緯度坐標(biāo)采集,然后通過算法精確轉(zhuǎn)換成“2000國家大地坐標(biāo)系”做數(shù)據(jù)分析與處理,同時(shí)可以方便的將采集到的地理信息在谷歌地球軟件中展現(xiàn),并查看對應(yīng)歷史影像。
保障性安居工程審計(jì)中涉及到的地理坐標(biāo)信息主要有兩大類,一類是以項(xiàng)目為單位的坐標(biāo)點(diǎn)數(shù)據(jù),另一類是以地塊為單位的坐標(biāo)區(qū)域數(shù)據(jù),這兩類坐標(biāo)數(shù)據(jù)的處理以及導(dǎo)入方式各不相同,故分別擬定以下標(biāo)準(zhǔn)表:
1.坐標(biāo)點(diǎn)數(shù)據(jù)標(biāo)準(zhǔn)表:
注:上表A、B、C列為關(guān)鍵字,下發(fā)此表時(shí)可將審計(jì)組取得的保障房項(xiàng)目名稱及地址匯總后一同提供給被審計(jì)單位對照填寫,以防出現(xiàn)項(xiàng)目名稱不一致的情況。此表其余列可根據(jù)數(shù)據(jù)分析內(nèi)容及審計(jì)需求自行擬定,或后期在SQL中根據(jù)C列關(guān)鍵字直接關(guān)聯(lián)分析。
2.坐標(biāo)區(qū)域標(biāo)準(zhǔn)表:
坐標(biāo)區(qū)域標(biāo)準(zhǔn)表應(yīng)按如下格式填寫。將地塊界址點(diǎn)坐標(biāo)的經(jīng)、緯度按標(biāo)準(zhǔn)格式拆分,分別填列至A、B列,同一地塊的界址點(diǎn)對應(yīng)的C列單元格均填充相同地塊名稱。
注:地塊數(shù)據(jù)由土地管理局地塊登記系統(tǒng)后臺導(dǎo)出后填列。
二、坐標(biāo)數(shù)據(jù)處理及導(dǎo)入
1.點(diǎn)坐標(biāo)
將采集的坐標(biāo)點(diǎn)數(shù)據(jù)匯總后直接導(dǎo)入KML制作工具“kmlgenerator”中生成KML文件,或直接將匯總的EXCEL表格另存為“CSV(逗號分隔符)”格式,以便下一步導(dǎo)入QGIS。
(直接生成KML文件)
(EXCEL另存為CSV格式)
2.地塊坐標(biāo)
(1)將采集的地塊坐標(biāo)數(shù)據(jù)復(fù)制到“界址坐標(biāo)KML文件生成工具”的“大地坐標(biāo)”表中(為避免坐標(biāo)數(shù)據(jù)意外損失,請選擇“黏貼值”),若需要添加字段可自行隨后添加列。
(界址坐標(biāo)KML文件生成工具)
(2)對于需要在QGIS中進(jìn)行數(shù)學(xué)計(jì)算的字段,在“字段類型”表中設(shè)置相應(yīng)類型。切換到“字段類型”表,點(diǎn)開第二行下拉框,選擇“float”?!白侄晤愋汀北頍o需人工填入數(shù)據(jù),當(dāng)在“大地坐標(biāo)”表中添加或刪除列后會自動在此表中更新相應(yīng)字段及類型,默認(rèn)類型為“string”,僅當(dāng)需要更改字段類型時(shí)才手動選擇。
(選擇字段類型)
注:當(dāng)“大地坐標(biāo)”表列標(biāo)題有變動時(shí),“字段類型”表中第二行數(shù)據(jù)會全部重置為string。故請先確定好“大地坐標(biāo)”表的列標(biāo)題,然后再在此表中設(shè)置字段類型。
(3)點(diǎn)擊“加載項(xiàng)”中的“生成地塊坐標(biāo)KML文件”按鈕,程序會自動生成地塊的坐標(biāo)數(shù)據(jù)并存入本表的存放路徑,文件名為“地塊坐標(biāo).kml”。
(按鈕位置)
(生成結(jié)果)
3.坐標(biāo)導(dǎo)入
將制作好的點(diǎn)坐標(biāo)及地塊坐標(biāo)KML文件導(dǎo)入QGIS軟件中進(jìn)行進(jìn)一步分析。
(1)首先在QGIS中新建一個(gè)項(xiàng)目。依次點(diǎn)擊“項(xiàng)目”—“新建”。
(新建項(xiàng)目)
(2)若需要查看坐標(biāo)點(diǎn)實(shí)際位置,還需要添加在線地圖。在資源瀏覽器中找到“XYZ Tiles”,展開節(jié)點(diǎn),雙擊“OpenStreetMap”,將地圖添加進(jìn)項(xiàng)目,此在線地圖為WSG84坐標(biāo)系,即谷歌地圖的坐標(biāo)系,與“2000國家大地坐標(biāo)系”通用。
(3)若要查看影像地圖,可手動添加地圖源。
(右擊XYZ Tiles,新建連接)
(命名地圖名稱自定)
(添加后即可在資源瀏覽器中找到)
(4)根據(jù)影像地圖可直接查看對應(yīng)項(xiàng)目坐標(biāo)點(diǎn)或地塊的實(shí)際情況。但是由于中國地區(qū)的谷歌影像數(shù)據(jù)更新較慢,目前很多區(qū)域仍是數(shù)年前的影像,故如必要,可先將坐標(biāo)點(diǎn)及地塊坐標(biāo)數(shù)據(jù)使用“GPS批量查詢轉(zhuǎn)換工具”,將“2000國家大地坐標(biāo)系”統(tǒng)一轉(zhuǎn)換為高德坐標(biāo)系,再導(dǎo)入QGIS中,然后使用高德影像地圖查看項(xiàng)目實(shí)際情況。
(5)坐標(biāo)數(shù)據(jù)的kml文件導(dǎo)入后,最好再將導(dǎo)入的圖層另存為shp格式,然后移除原kml圖層,以免后期進(jìn)行數(shù)據(jù)操作時(shí)出錯。
(圖層另存為)
(格式選擇ESRI Shapfile,設(shè)置shp文件存儲路徑,坐標(biāo)參照系一般默認(rèn))
(勾選“將已保存的文件添加到地圖中“,勾選欲保存的字段。)
(保存后在圖層管理器中會看到兩個(gè)相同圖層,新添加的那個(gè)即是shp圖層,將原有kml圖層移除。)
三、數(shù)據(jù)分析及結(jié)果
1.篩選出地處偏遠(yuǎn)的項(xiàng)目
(1)在QGIS中分別測量城市規(guī)劃區(qū)經(jīng)度和緯度方向的跨度,取其平均值作為該區(qū)域的直徑R。
(2)將匯總的項(xiàng)目坐標(biāo)信息復(fù)制到“GPS批量查詢轉(zhuǎn)換工具”,查詢項(xiàng)目點(diǎn)至政府直線距離及所在行政區(qū)劃,同時(shí)核實(shí)項(xiàng)目點(diǎn)行政區(qū)劃是否有誤,保存結(jié)果,形成“表D01”。在表D01中篩選出距政府直線距離大于(3/4)R(具體距離可根據(jù)實(shí)際情況自行調(diào)整)的項(xiàng)目。
(3)在QGIS中篩選出未落在規(guī)劃區(qū)內(nèi)的項(xiàng)目點(diǎn),并將篩選結(jié)果導(dǎo)出。
(依次點(diǎn)擊菜單欄的設(shè)置——選項(xiàng),在彈出的界面中選擇最下面一個(gè)標(biāo)簽? ,點(diǎn)開“通用”節(jié)點(diǎn),雙擊“要素過濾器無效”右側(cè),將選項(xiàng)改為“不過濾”,設(shè)置完成后點(diǎn)“OK”。)
(依次點(diǎn)擊矢量——地學(xué)數(shù)據(jù)處理工具——差分。)
(按圖示選擇項(xiàng)目坐標(biāo)點(diǎn)及城市規(guī)劃區(qū)圖層)
(點(diǎn)擊“差分”欄右邊按鈕,選擇比對結(jié)果保存路徑及格式,一般選擇xlsx,取消勾選“執(zhí)行算法后打開輸出文件”,點(diǎn)擊下方“運(yùn)行”即可生成結(jié)果。結(jié)果表中的項(xiàng)目即為未在規(guī)劃區(qū)內(nèi)的項(xiàng)目)
(4)將(2)、(3)中兩部分篩選結(jié)果數(shù)據(jù)匯總,即為“偏遠(yuǎn)項(xiàng)目表”。
(5)針對篩選出的偏遠(yuǎn)項(xiàng)目重點(diǎn)核實(shí),是否存在定向分配、違規(guī)將棚改政策擴(kuò)大到一般建制鎮(zhèn)等問題。
2.篩選出空置地塊
(1)在QGIS中篩選保障房劃撥地塊內(nèi)無項(xiàng)目坐標(biāo)點(diǎn)的地塊。
(依次點(diǎn)擊矢量——分析工具——計(jì)算多邊形內(nèi)點(diǎn)的數(shù)目。)
(設(shè)置需要比對的地塊及點(diǎn)坐標(biāo)圖層)
(在“計(jì)數(shù)”欄選擇比對結(jié)果保存路徑及格式,點(diǎn)擊下方“運(yùn)行”即可生成結(jié)果。)
(2)將分析結(jié)果導(dǎo)出,導(dǎo)出表中“NUMPOINTS”列即為對應(yīng)地塊中的坐標(biāo)點(diǎn)個(gè)數(shù),值為“0”的即為空置地塊,將此部分地塊篩選整理出來,形成“閑置地塊表”。
(閑置地塊表)
(3)在“閑置地塊表”中篩選出距今超過三年的地塊,針對此部分已劃撥但閑置超過三年的地塊重點(diǎn)核查原因。
通過地理信息技術(shù)應(yīng)用于保障性安居工程審計(jì),突顯先進(jìn)技術(shù)手段給審計(jì)分析過程中提供大量化、多樣化和快速化等特點(diǎn),不僅大量節(jié)省了審計(jì)成本,還將原來一些難以認(rèn)定的情形實(shí)現(xiàn)具體可操作化,提高了審計(jì)質(zhì)量和效率。在實(shí)際審計(jì)中,審計(jì)人員應(yīng)同時(shí)結(jié)合傳統(tǒng)審計(jì)方式,借助地理信息技術(shù)分析出的疑點(diǎn),進(jìn)一步落實(shí)保障性安居工程項(xiàng)目選址不科學(xué),建設(shè)位置偏遠(yuǎn)導(dǎo)致空置、違規(guī)擴(kuò)大棚改范圍以及保障性安居工程地塊閑置等問題,鞏固完善審計(jì)證據(jù)鏈。
(作者單位:景德鎮(zhèn)市審計(jì)局、江西省審計(jì)廳)