一、基本情況
近年來,憑借強大的空間分析能力,地理信息系統(tǒng)(GIS)技術在自然資源資產審計工作中的應用越來越廣泛。地理信息系統(tǒng)中常用數(shù)據(jù)包括柵格數(shù)據(jù)和矢量數(shù)據(jù)。常說的遙感影像就屬于柵格數(shù)據(jù)。常用的矢量數(shù)據(jù)一般是指Shapefile格式(*.shp)的包含空間圖形和屬性表格的一組數(shù)據(jù),圖形加屬性表的配置使其可以填充容納大量的數(shù)據(jù)信息。
自然資源資產審計經常會使用土地利用方面的矢量數(shù)據(jù),審計人員將疑點進行現(xiàn)場核實后,可使用GPS工具箱等App拍攝帶經緯度的現(xiàn)場照片,常需要將外業(yè)采集的或審計資料中的坐標點導入地理信息系統(tǒng)常用的ArcGIS軟件,并對以上數(shù)據(jù)進行空間疊加分析。
本文介紹的是使用ArcMap 10.2.2將坐標點落圖的常用方法。常用的第三次國土調查、林草生態(tài)綜合監(jiān)測成果、三區(qū)三線等資環(huán)相關數(shù)據(jù)的XY通常以米為單位,因此可以通過將經緯度格式的XY數(shù)據(jù)轉換為米、跨帶數(shù)據(jù)統(tǒng)一帶號等操作落圖,便于審計人員將各類數(shù)據(jù)疊加使用。
二、零散點快速跳轉
當審計人員只需要臨時查詢某一坐標點位置地類時,可進行臨時跳轉并按需添加注記。
(一)轉到XY
在ArcMap新建的地圖文檔中先添加遙感影像,使用工具欄“轉到XY”工具,展開右邊單位選項欄,選擇XY單位,輸入XY數(shù)值后回車進行快速查找,目標點會在地圖中閃爍并居中。
(二)臨時標注點
審計人員可選擇點擊“添加點”在地圖中顯示該點,點擊“添加標注點”在地圖中顯示該點及其坐標。如果查看地類后需要取消顯示該標注,可在地圖中右鍵選擇“數(shù)據(jù)框屬性”,在“注記組”中取消勾選框。
三、點集批量處理
通常,審計人員會外業(yè)采集較多的坐標點,并期望知道現(xiàn)場核實地塊的大致面積。這時可以選擇批量導入坐標點,并生成采集點所圍成圖形的面積。
(一)十進制度轉為米進行落圖
審計人員使用GPS工具箱等手機App拍攝現(xiàn)場照片的XY坐標通常為十進制的經緯度,需要設置投影坐標系轉換為米。
1.編輯經緯度信息
在Excel中新建*01.x1s文件,新建列X、Y填寫外業(yè)核實地點的經緯度,可根據(jù)需要填寫其他屬性,為防止導入過程中出錯,建議表頭避免使用漢字,可用字母、數(shù)字代替。
表1.新建x1s格式的表格填入十進制度
|
ID
|
X
|
Y
|
BZ
|
|
1
|
124.……
|
49.……
|
節(jié)點1
|
|
2
|
……
|
……
|
節(jié)點2
|
|
……
|
|
|
|
2.顯示XY數(shù)據(jù)
ArcGIS中新建空白地圖文檔,添加*01.xls。
在內容列表中右鍵單擊表Sheet1$選擇“顯示XY數(shù)據(jù)”。
在“指定X、Y和Z坐標字段”中X字段選擇“X”、Y字段選擇“Y”。“輸入坐標的坐標系”中點擊編輯,選擇“地理坐標系→Asia→China Geodetic Coordinate System 2000”,點擊右上角☆可添加至收藏便于下次使用。
3.導出并進行投影
右鍵單擊生成的Sheet1$個事件→數(shù)據(jù)→導出數(shù)據(jù),輸出要素類*02.shp文件,將導出的矢量數(shù)據(jù)添加到地圖文檔中。
打開ArcToolbox工具→數(shù)據(jù)管理工具→投影和變換→要素→投影。“輸入數(shù)據(jù)集或要素”選擇*02.shp。“輸出數(shù)據(jù)集或要素”選擇*03.shp。
展開“輸出坐標系”右側的空間參考屬性卡。選擇“投影坐標系→Gauss Kruger→CGCS2000→CGCS2000_3_Degree_GK_Zone_41(選擇數(shù)據(jù)所在地對應的帶號)”。
4.添加遙感影像
新建空白地圖文檔。先添加遙感影像數(shù)據(jù),再添加*03.shp即可落圖。
5.顯示XY(米)坐標
按照需要,可選擇計算并顯示經緯度坐標點對應的X、Y(米)坐標值:在矢量數(shù)據(jù)“表選項”中添加字段“X1”和“Y1”,右鍵點擊新添加的字段選擇“計算幾何”,屬性分別選擇“點的X坐標”和“點的Y坐標”,單位選擇m(米)。
圖層右鍵勾選標注要素,打開圖層屬性,在標注表達式中輸入[X1]&","&[Y1],即可在地圖中顯示XY的值。
6.計算圖形面積
當需要計算外業(yè)采集坐標點連接形成幾何圖形的面積時。打開ArcToolbox工具→數(shù)據(jù)管理工具→投影和變換→要素→“點集轉線”。再將形成的閉合線數(shù)據(jù)作為輸入要素,選擇“要素轉面”,生成面數(shù)據(jù)。
打開面數(shù)據(jù)屬性表,在“表選項”中添加字段“面積”,右鍵“計算幾何”,選擇所需單位,計算圖形面積。
(二)同帶號XY(米)落圖
當XY坐標單位為米且2000坐標系三度帶的帶號與遙感影像一致時,落圖時基本與前文所述十進制度轉為米的流程相同,差別在于進行“顯示XY數(shù)據(jù)”操作時,可直接選擇投影坐標系,省去在工具箱中投影的操作,落圖方法較為快捷。
(三)跨帶號疊加XY(米)
需要注意的是,當被審計地區(qū)東西方向長度較大、投影坐標系跨多個2000坐標系三度帶時,直接跳轉點的XY坐標會產生跨帶錯位,坐標點落圖時統(tǒng)一帶號方便在同一地圖文檔中疊加使用。可選取當?shù)卣急容^大的帶號為轉換目標。以下操作以XY單位同為米、將帶號40數(shù)據(jù)顯示在帶號41的遙感影像上為例。
1.添加遙感影像
在ArcMap中新建空白地圖文檔,添加投影坐標系為2000坐標系三度帶、帶號41的遙感影像。
2.顯示XY數(shù)據(jù)
在Excel中新建x1s文件,新建列X、Y填寫40帶號X、Y坐標數(shù)據(jù)。將表格添加至ArcMap。單擊表右鍵選擇“顯示XY數(shù)據(jù)”。X字段選擇“X”、Y字段選擇“Y”。“輸入坐標系”選擇“CGCS2000_3_Degree_GK_Zone_40”。
3.導出矢量數(shù)據(jù)
右鍵生成的事件→數(shù)據(jù)→導出數(shù)據(jù),此時40帶號的坐標點已經顯示在影像對應的正確位置,審計人員可對生成的shp數(shù)據(jù)進行編輯或繼續(xù)空間分析。
4.轉換投影
需要計算41帶號坐標值時,可繼續(xù)轉換投影坐標系。打開ArcToolbox工具→數(shù)據(jù)管理工具→投影和變換→要素→投影。輸出坐標系”選擇“CGCS2000_3_Degree_GK_Zone_41”。此時在屬性表中重新計算幾何,X坐標值即為41開頭。
四、小結
以上便是審計人員常用的坐標點落圖方法,依據(jù)被審計單位提供的電子數(shù)據(jù)并借助ArcMap軟件可以直觀查看已落圖坐標點的歷年遙感影像,審計人員通過影響紋理可以初步判斷地類,疊加其他數(shù)據(jù)例如某項目區(qū)范圍可以繼續(xù)進行空間分析。