格网数据交换格式NSDTM

格网数据交换格式NSDTM-DEM如何转换到SuperMap的数据源中
前天有个好朋友，从他客户那里拷了一些DEM数据（文件名后缀为.dem），据说他的用户希望能把这个DEM数据转换为EOO，但是不知道怎么转。

我拿到这个数据后，先简单地通过SuperMap Deskpro新建数据源，导入数据集（SuperMap中直接支持DEM数据文件的导入）,可是导入失败。

直觉告诉我，应该是此dem文件中记录的是非SuperMap支持的那种dem格式信息。

我试着用记事本打开这个dem文件，还好，这个文件是明码格式的，而且存储比较规范，大概格式如下
NSDTF-DEM
1.0
M
0.000000
0.000000
39512435.000000
2743120.000000
5.000000
5.000000
985
1328
100
我以前也没看到过这种数据，不过见过ArcInfo的Grid交换数据文件，经验告诉我39512435.000000，27431 20.000000信息为栅格数据的左上角（或者左下角）坐标；5.0为分辨率；985和1328应该为这个DEM数据的行列数，上网搜索了一下NSDTF-DEM，找到了包含NSDTF-DEM的相关介绍，原来这是格网数据交换格式文件，而且是国际标准的，果然跟我上面的判断一样，而且明确了具体的含义，这里我列出来有兴趣的可以详细看一下：
DataMark--------中国地球空间数据交换格式-格网数据交换格式(CNSDTF-RAS或CNSDTF-DEM)的标志。

基本部分，不可缺省。

Version--------该空间数据交换格式的版本号,如1.0。

基本部分，不可缺省。

Unit--------坐标单位,K表示公里,M表示米,D表示以度为单位的经纬度,S表示以度分秒表示的经纬度(此时坐标格式为DDDMM SS.SSSS, DDD为度, MM为分, SS.SSSS为秒)。

基本部分，不可缺省。

Alpha--------方向角。

基本部分，不可缺省。

Compress--------压缩方法。

0表示不压缩,1表示游程编码。

基本部分，不可缺省。

Xo--------左上角原点X坐标。

基本部分，不可缺省。

Yo--------左上角原点Y坐标。

基本部分，不可缺省。

DX--------X方向的间距。

基本部分，不可缺省。

DY--------Y方向的间距。

基本部分，不可缺省。

Row--------行数。

基本部分，不可缺省。

Col--------列数。

基本部分，不可缺省。

HZoom--------高程放大倍率。

基本部分，不可缺省。

设置高程的放大倍率,使高程数据可以整数存贮,如高程精度精确到厘米,高程的放大倍率为100。

如果不是DEM则HZoom为1。

上面这些都是DEM数据的头文件信息，接下来每行记录的就是相应的DEM栅格单元值了。

比如：
-99999 20738 19502 18176 16504 14833 1368
5 12564 11311 11307
11302 11298 11785 12149 12145 12140 1201
3 11763 11635 11630
11625 11621 11617 11612 11608 11603 1135
3 11102 11097 11093
...
这些数据应该都是按照前面记录的行列数来进行组织的，先不管这些。

从机器上找了一个SuperMap的示范Grid数据，在Deskpro中打开保存为ArcInfo的Grid交换格式，用记事本打开后，对比看了一下，其实文件后部分栅格单元数据部分差不多，差别就是ArInfo Grid交换格式文件头主要按下面方式记录了相关信息：
ncols 319----------------数据列数
nrows 198----------------数据行数
xllcorner 131985----------------数据左上角的X值
yllcorner 19020----------------数据左上角的Y值
cellsize 5----------------数据分辨率（栅格单元的宽高）
NODATA_value -99999----------------无值数据标志
对比出这样的不同，转换起来就方便了，依葫芦画瓢，把原先DEM数据的文件头信息改为了ArcInfo Grid的文件头格式（如下），
ncols 1328
nrows 985
xllcorner 39512435
yllcorner 2743120
cellsize 5
NODATA_value -99999
并把dem文件另存为grd后缀的文件，再通过Deskpro导入这个另存后的grd文件后，能导入成功了，而且看到了比较好的效果。

不过出现了新的问题，我在把上面文件临接的文件感觉没拼在一起。

后来感觉应该是yllcorner的值没计算好，试着把yllcorner值2743120减去行数985*分辨率5 后再转换就没有问题（这个过程非常重要，有类似文件需要处理的朋友一定要注意）。

格网数据交换格式文件中还记录了HZoom（高程放大倍率），我们只要在Deskpro的分析菜单，栅格分析子菜单中找到栅格代数运算，把转换进来的GRID数据结合这个高程放大倍率进行一个运算保存为新的数据集就可以了。

前面提到需要把格网数据交换格式DEM文件转换为DEM。

在经过前面格网数据交换格式如何转换到SuperMap的数据源中介绍的方法转换后，DEM数据转换为SuperMap 中的一个Grid数据集了。

SuperMap也提供了比较高级的分析功能，可以从Grid数据中提取等值线，或者提取指定等值线，我们就可以利用这个功能（可以在Deskpro的分析菜单中找到）从转入的Grid数据，提取出等值线，有了等值线数据集，就可以将其转换为E00数据了（导出线数据集时选择转出文件类型为ArcInfo 交换格式）。

BIL, BIP, and BSQ raster files
Band interleaved by line (BIL), band interleaved by pixel (BIP), and band sequenti al (BSQ) are three common methods of organizing image data for multiband images. B IL, BIP, and BSQ are not in themselves image formats but are schemes for storing t he actual pixel values of an image in a file. These files support the display of s ingle and multiband images and handle black-and-white, grayscale, pseudo color, tr ue color, and multispectral image data.
The BIL, BIP, and BSQ files are binary files, and they must have an associated ASC II file header to be interpreted properly by ArcGIS. This header file contains anc illary data about the image such as the number of rows and columns in the image, i f there is a color map, and latitude and longitude.
Band interleaved by line data stores pixel information band by band for each line, or row, of the image. For example, given a three-band image, all three bands of d ata are written for row 1, all three bands of data are written for row 2, and so o n, until the total number of rows in the image is reached.
BIL、 BIP 、或BSQ文件可以提供有三个描述文件：一个头文件（ HDR ）：描述图像像素数据的布局，必须提供；彩色文件（clr）：描述的图像彩色地图；和描述每一个波段图像统计资料的统计文件（ STX ）。

这些ASCII文本文件都可以由一个文本编辑器生成使用这些你知道的图像信息。

coverage
一个coverage是georelational数据模型的数据存储载体，它同时包含空间（地点）和属性（描述性的）数据的地理特征。

coverages使用一套特征类代表的地理特征。

每个特征类存储一套点，线（弧），多边形，或批注（文字）。

coverages可以拓扑结构，它决定特征之间的关系。