可能和很多朋友一樣,我們只能拿到紙質(zhì)地形圖,還得自己掃描,之后再在ArcMap里面配準。對于剛開始接觸的人來說,真是一頭霧水。
我以前的誤區(qū):"剛接觸ArcGIS,感覺挺簡單的,不就先掃描,然后在ArcMap中加載數(shù)據(jù),再Georeferencing配準。而實地采集的GPS數(shù)據(jù)就是通過excel轉(zhuǎn)換到Access數(shù)據(jù)庫后加載到arcmap里面,再display data。這樣,點就投到圖上了。"
呵呵,剛開始,還真以為是這樣做呢。不過到Layout view界面,Insert Scale bar后,非常困惑,才發(fā)現(xiàn)怎么插入的比例尺和原圖中的對不上啊。經(jīng)過多方學習,才知道其中還有許多深奧的知識,任何一點沒注意都會出差錯。下面就與大家分享一下我的學習總結(jié)。
首先要知道什么是地理坐標和大地坐標,簡單的說地理坐標就是球面坐標,而大地坐標是平面坐標。要深入了解這兩者的區(qū)別,大家可以到網(wǎng)上搜相關(guān)資料。
地理坐標到大地坐標的過程叫做投影。
用1:50000地形圖來舉個例子,是北京1954坐標系,6度分帶的高斯克呂格投影。
上面有兩種坐標,一種是常說的方里網(wǎng),就是間隔為1公里的網(wǎng)格。還有就是四角上的經(jīng)緯度。
我們常用的GPS采集的數(shù)據(jù)都是WGS84坐標系的,是常見的度分秒格式。這里要特別注意,地形圖上的經(jīng)緯度不是WGS-84的經(jīng)緯度,而是北京1954基準面下的經(jīng)緯度。而且,這個經(jīng)緯度也不是球面坐標,而是經(jīng)過高斯克呂格投影的平面坐標。論壇上常提及的wgs-84到beijing 1954的轉(zhuǎn)換就是這個問題,兩者的換算涉及到圖幅所在分帶保密點的問題,轉(zhuǎn)換參數(shù)并未公開。在我國,這兩個坐標系同一點經(jīng)緯度有幾十米的差別。
了解了上面的知識后,就開始在arcmap中進行圖像配準了,這里假設(shè)大家都知道基本操作步驟,因此我只注重于工作內(nèi)容的闡述。
1, 首先把圖復制到工作目錄下
2, 在ArcCatalog中找到需要配準的圖,給他定義一個投影系統(tǒng)(注意和投影變換的區(qū)別)。這里我們用的是1:50000的地形圖,它是基于北京1954坐標系,6度分帶的高斯克呂格投影。在地形圖方里網(wǎng)上可以看出本圖幅位于哪個分度帶,這里假設(shè)是19。因此我們要選擇的是Beijing 1954 GK Zone 19.prj。同時目錄里面還有一個Beijing 1954 GK Zone 19N.prj,這個是用于沒有分度帶號的。而我們的圖幅是包括分度帶號。這個可以在網(wǎng)上搜索 "arcgis 坐標文件" 來進一步了解。
3, 給配準的圖定義了一個投影系統(tǒng)后,把圖加載到ArcMap里面,打開GeoReferencing工具,直接利用方里網(wǎng)交點進行配準。這里要注意的問題是,地形圖上的方里網(wǎng)坐標為公里,而需要輸入的應(yīng)該是米。所以要在方里網(wǎng)對應(yīng)坐標后面加000。如地形圖上讀出一個交點為(19387, 3420),19387的19為分帶號,也要一并輸入,那么這個點應(yīng)該輸入(19387000, 3420000)。其余配準過程和網(wǎng)上教程一致。
4, 配準后可以在Lay Properity里面把顯示單位改成度分秒,于是地圖就以經(jīng)緯度格式顯示了。此時的經(jīng)緯度是基于北京1954基準面的,可以把光標指向四角的經(jīng)緯度標記,以核對配準的精確度。
到此為止,地形圖已經(jīng)精確配準了。接下來往里面加載GPS數(shù)據(jù)。
1, 用你能想到的方式把GPS接收機上的數(shù)據(jù)導出來,除了一個一個手工輸入,呵呵。
2, 盡量建立起數(shù)據(jù)庫吧,這樣方便些
3, ArcGIS中加載所得到的數(shù)據(jù)庫,用Defined Query來選擇自己所需要顯示的數(shù)據(jù)
4, 如果需要精確地在地形圖上顯示GPS數(shù)據(jù),那么還需要一個WGS-84到北京1954坐標的轉(zhuǎn)換。有轉(zhuǎn)換參數(shù)最好,如果沒有,但的確需要,就去網(wǎng)上查一下 手持GPS接收機轉(zhuǎn)換參數(shù)估算 這篇文章。我們這里是假設(shè)不需要精確顯示,直接把WGS-84數(shù)據(jù)當成北京1954來用。前面也說了,有一定誤差,但是不大。因此直接跳到5。
5, Display X,Y Datas, 定義為WGS-84坐標系統(tǒng)。
(其實此處我有疑問,因為現(xiàn)在的地形圖是投影后的圖,按我的理解,也應(yīng)該把是球面坐標的WGS-84數(shù)據(jù)進行投影變換后才能用。但是查了很多資料,都是在這里直接定義數(shù)據(jù)來源是WGS-84就可以了,也許是ArcGIS在投點的時候,可以自動進行變換。)
6, OK,現(xiàn)在點也投上去了,ArcGIS工作中最基礎(chǔ)的工作也就完成了。:)
再轉(zhuǎn)到Layout view 界面,插入一個比例尺,發(fā)現(xiàn)和地形圖上的比例尺非常吻合:-),把顯示單位改成kilometer,用measure tool 就可以量地形圖上任意點的距離了,試了一下公里網(wǎng)格,非常精準。