NSDTF-DEM格式显示转换方法--兔八哥的Blog

页面载入中....

NSDTF-DEM格式显示转换方法　

By 兔八哥发表于 2007-11-7 8:36:00

去年从浙江测绘局拿的DEM数据说是GeoStar格式的，因为没有这个软件一直没有用。昨天打开一看原来是明码格式，文件头如下：

NSDTF-DEM
1.0
M
0.0
0.0
493997.500000
3380167.500000
5.000000
5.000000
927
1203
100

看起来和ESRI的GRID ASC明码格式相仿，网上搜一下NSDTF-DEM是中华人民共和国国家标准地球空间数据交换格式意思，估计国产GIS软件应该都支持。不过ArcGIS就不支持了。网上有关于NSDTF DEM转GRID的详细方法  http://www.v2sky.com/u/v2sky/archives/2007/36.html

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。

而ESRI GRID交换格式文件头如下：

ncols         数据列数
nrows         数据行数
xllcorner     数据左上角的X值
yllcorner     数据左上角的Y值
cellsize      数据分辨率
NODATA_value  无值数据标志

按照v2sky上的思路，将NSDTF-DEM文件头修改成ESRI GRID交换格式文件头，果然ArcGIS就可以识别了，不过新的数据中还需要除一个HZOOM值。我是在在ArcMap里面用Raster Calculator计算一下，切记不能取整数，不然会丢很多信息的。

v2sky这篇博客上海提醒需要yllcorner不能简单的等于DY，正确的应该是yllcorner=DY-Row*减去行数985*cellsize，不然会拼接时候会出问题，这个没有试验，不过相信作者说的没有错。

阅读全文 | 回复(2) | 引用通告 | 编辑

标签：NSDTF DEM