如何在ArcGIS中使用後綴.dem的數據

轉自：http://www.gisall.com/html/07/5607-2111.html

拿到一批文件，後綴是.dem，領導要求轉成ArcGIS可以識別的格式。剛開始也是一直沒有頭緒，後來發現了一文章，再加上support.esri.com一搜，終於解決問題啦。總結如下：

首先，用記事本打開.dem，如圖，發現是明碼格式，如圖

發現是NSDTF-DEM，上網一搜，得到一文章，說明NSDTF-DEM和ESRI-GRID的文件頭的轉換關係，將此文引用如下：

==================以下部分引用=====================

NSDTF-DEM格式數據

Posted on 2008-10-09 11:29 cosophy 閱讀(815) 評論(0)  編輯收藏 網摘 所屬分類: GIS

【參考文章】 ：http://www.v2sky.com/u/v2sky/archives/2007/36.html http://www.gissky.cn/blog.asp?ID=146#more

NSDTF-DEM是一種國家標準空間數據的交換格式，後綴一般是.dem，使用記事本打開會發現它是明碼。 如：

頭文件信息：

NSDTF-DEM
1.0
M
0.000000
0.000000
39512435.000000
2743120.000000
5.000000
5.000000
985
1328
100

下面就是具體的柵格數據值：

    -99999      20738      19502      18176      16504      14833      13685      12564      11311      11307
     11302      11298      11785      12149      12145      12140      12013      11763      11635      11630
     11625      11621      11617      11612      11608      11603      11353      11102      11097      11093

頭文件每行數據對應的含義爲：

DataMark--------中國地球空間數據交換格式-格網數據交換格式(CNSDTF-RAS或CNSDTF-DEM)的標誌。基本部分，不可缺省。
Version--------該空間數據交換格式的版本號,如1.0。基本部分，不可缺省。
Unit--------座標單位,K表示公里,M表示米,D表示以度爲單位的經緯度,S表示以度分秒錶示的經緯度(此時座標格式爲DDDMMSS.SSSS, DDD爲度, MM爲分, SS.SSSS爲秒)。基本部分，不可缺省。
Alpha--------方向角。基本部分，不可缺省。
Compress--------壓縮方法。0表示不壓縮,1表示遊程編碼。基本部分，不可缺省。
Xo--------左上角原點X座標。基本部分，不可缺省。
Yo--------左上角原點Y座標。基本部分，不可缺省。
DX--------X方向的間距。基本部分，不可缺省。
DY--------Y方向的間距。基本部分，不可缺省。
Row--------行數。基本部分，不可缺省。
Col--------列數。基本部分，不可缺省。
HZoom--------高程放大倍率。基本部分，不可缺省。設置高程的放大倍率,使高程數據可以整數存貯,如高程精度精確到釐米,高程的放大倍率爲100。如果不是DEM則HZoom爲1。

要使用這種數據就需要轉換成目標軟件支持的數據格式了,如轉換成ArcGIS支持的GRID數據，就可以使用以下方法。

由於柵格單元數據值記錄方式基本一樣，主要是頭文件信息不同。如ArcGIS GRID數據的文件頭按如下方式記錄相關信息：

ncols         319----------------數據列數
nrows         198----------------數據行數
xllcorner     131985----------------數據左上角的X值
yllcorner     19020----------------數據左上角的Y值
cellsize      5----------------數據分辨率（柵格單元的寬高）
NODATA_value  -99999----------------無值數據標誌

所以，對於以上的DEM數據，將文件頭改成下面的形式並保存爲.grd文件就可以在ArcGIS中打開使用了。

ncols         1328
nrows         985
xllcorner     39512435
yllcorner     2743120
cellsize      5
NODATA_value  -99999

注意：按原創作者的意思應把yllcorner的值賦爲Yo-Row*DY,不然無法拼接相鄰圖像。 爲什麼這麼做他沒說，難道是其中某種數據文件頭記錄的是左下角的座標值，而不是左上角？？

此外，HZoom值好像沒有使用到，可以在ArcMap裏面用Raster Calculator計算轉換後的GRID（注意不能填整數，不然會丟失很多信息），使它和原始數據保持一致。

=================以上引用=========================

OK，知道之間的關係之後，就直接改一下文件頭，然後存成.grd格式。但轉完之後，發現.grd格式，ESRI是不認的。後來，想了一個辦法，將後綴轉成txt，然後用Toolbox中的ASCII to Raster工具，終於將.dem的數據轉換出來了。

在ArcScene中加載，成果如下：

OVER了，呵呵~~