摘要:本文以延吉市DEM為例,探NNSDTF與usGs之間的轉(zhuǎn)換問題。
關(guān)鍵詞:NSDTF-DEM、USGS-DEM、Global Mapper8
1、概述
DEM是數(shù)字高程模型(DigitalElevafion Model的簡(jiǎn)寫),是以數(shù)字形式按一定結(jié)構(gòu)組織在一起的,表示實(shí)際地形特征空間分布的一種實(shí)體地面模型,也是地形地理位置和起伏的數(shù)字描述。由于DEMN述的是地面高程信息,它在測(cè)繪、水文、氣象、地貌、地質(zhì)、土壤、工程建設(shè)、通訊、氣象、軍事等國(guó)民經(jīng)濟(jì)和國(guó)防建設(shè)以及人文和自然科學(xué)領(lǐng)域有著廣泛的應(yīng)用,它也是數(shù)字城市的核心數(shù)據(jù)。延吉市虛擬三維系統(tǒng)是在DEM基礎(chǔ)上建立的,地形的制作、通視分析及淹沒分析等都離不開它,由于DEM大部分是USGS-DEM(USGS是美國(guó)地質(zhì)調(diào)查局(U.S.Geological Survey)的英文縮寫,USGS負(fù)責(zé)管理美國(guó)全國(guó)的數(shù)字地圖數(shù)據(jù)的采集與分發(fā)。USGS DEM數(shù)據(jù)以ASCII碼形式存儲(chǔ),是一種公開格式的DEM數(shù)據(jù)格式標(biāo)準(zhǔn),使用范圍較廣格式的,而我們購(gòu)買的都是NSDTF-DEM格式,NSDTF-DEM是中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)地球空間數(shù)據(jù)交換格式,是屬于格網(wǎng)數(shù)據(jù)交換格式,一般的GIS軟件都不支持這種格式,因此我們將討論利用Global Mapper8軟件,將NSDTF-DEM格式轉(zhuǎn)換為USGS-DEM格式。
2、設(shè)計(jì)思想與數(shù)據(jù)轉(zhuǎn)換
與其他數(shù)據(jù)格式轉(zhuǎn)換一樣,首先需要透徹了解數(shù)據(jù)格式,找到溝通兩者的方法和橋梁,制定設(shè)計(jì)方案。
Global Mapper8是一款地圖繪制軟件,不僅能夠?qū)?shù)據(jù)顯示為光柵地圖、高程地圖、矢量地圖,還可以對(duì)地圖作編輯、轉(zhuǎn)換、打印、記錄GPS及利用數(shù)據(jù)的GIS(地理信息系統(tǒng))功能直接訪問USGS(美國(guó)地質(zhì)勘探局)衛(wèi)星照片TerraServer數(shù)據(jù)庫(kù)和Global Mapper內(nèi)部的地形圖及以真實(shí)的3D方式查看高程地圖的功能。我們知道Global Mapper可以把ARClI轉(zhuǎn)化為USGS-DEM格式的DEM,只要把NSDTF-DEM的格式轉(zhuǎn)為AR C II格式的文件就可以,所以對(duì)NSDTF-DEM進(jìn)行處理一下。
我們先用記事本打開一個(gè)NSDTF-DEM格式的DEM文件
根據(jù)中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)地球空間格網(wǎng)數(shù)據(jù)我們知道,第一行NSDTF-DEM告訴我們這個(gè)DEM是中國(guó)地球空間數(shù)據(jù)交換格式一格網(wǎng)數(shù)據(jù)交換格式(CNSDTF-RAS或CNSDTF-DEM)的標(biāo)志;第二行1.0是該空間數(shù)據(jù)交換格式的版本號(hào);第三行M是坐標(biāo)單位,M表示米,D表示經(jīng)緯度;第四行與第五行是方向角和壓縮方式;第六行是左上角原點(diǎn)x坐標(biāo);第七行左上角原點(diǎn)Y坐標(biāo);第八行與第九行是X、Y方向的間距;第十行與第十一行是行數(shù)與列數(shù);第十二行1000是高程放大倍率。而剩下的數(shù)據(jù)就是高程數(shù)據(jù)。我們知道ARCINFO GRID的頭文件是以下面方式記錄的
ncols 319 數(shù)據(jù)列數(shù)
nrows 198 數(shù)據(jù)行數(shù)
xllcorner 131985 數(shù)據(jù)左上角的x值
yncorner 19020 數(shù)據(jù)左上角的Y值
ceHsize 5 數(shù)據(jù)分辨率
(柵格單元的寬高)

NODATA_value-99999 無(wú)值數(shù)據(jù)標(biāo)志
所以我們把DEM數(shù)據(jù)另存為,txt數(shù)據(jù),并把文件頭改為下面的格式
我們要注意yllcorner的值改為左上角原點(diǎn)Y坐標(biāo)一列數(shù)+數(shù)據(jù)分辯率,如果想要本地坐標(biāo),可以改動(dòng)xllcorner和yllcorner。
打開global mapper8,單擊open yourOWll data files,文件類型選擇allfile(*,*),打開剛才改動(dòng)的TXT文件,彈出一個(gè)窗口,按確定,彈出select projection窗口,根據(jù)本地坐標(biāo)和投影進(jìn)行修改,比如,延吉市是北緯129°,所以在ZONE中選擇52(126°E-132°E-northernhemisphere),由于我們使用的是1954北京坐標(biāo)系,所以在datum中選擇BEI JING1954,由于高程放大率為1000,所以在elevationunits中選擇millimeters,如果很多TXT文件,選擇use selected projectionfor all selected files,修改之后按確定,高程圖顯現(xiàn)在我們眼前,打開FILE-EXPROTRAASTER AND ELEVATION DATA-ExPORT DEM,彈出一個(gè)窗口,選擇“是”,彈出DEM EXPORT OPTIONS窗口,點(diǎn)擊確定,彈出另存為的窗口,起個(gè)名,一個(gè)USGS-DEM格式的DEM就生成了。
結(jié)束語(yǔ)
應(yīng)用這個(gè)方法已經(jīng)完成延吉市80平方公路的DEM數(shù)據(jù)轉(zhuǎn)換工作,希望對(duì)你的工作也有幫助。