格网数据交换格式NSDTM

1、 格网数据交换格式NSDTM-DEM如何转换到SuperMap的数据源中    前天有个好朋友，从他客户那里拷了一些DEM数据（文件名后缀为.dem），据说他的用户希望能把这个DEM数据转换为EOO，但是不知道怎么转。        我拿到这个数据后，先简单地通过SuperMap Deskpro新建数据源，导入数据集（SuperMap中直接支持DEM数据文件的导入）,可是导入失败。直觉告诉我，应该是此dem文件中记录的是非SuperMap支持的那种dem格式信息。我试着用记事本打开这个dem文件

2、，还好，这个文件是明码格式的，而且存储比较规范，大概格式如下NSDTF-DEM 1.0 M 0.000000 0.000000 39512435.0000002743120.0000005.000000 5.000000 985 1328 100         我以前也没看到过这种数据，不过见过ArcInfo的Grid交换数据文件，经验告诉我39512435.000000，2743120.000000信息为栅格数据的左上角（或者左下角）坐标；5.0为分辨率；985和1328应该为这个DEM数据的行列数，上网搜索了一下NSD

3、TF-DEM，找到了包含NSDTF-DEM的相关介绍，原来这是格网数据交换格式文件，而且是国际标准的，果然跟我上面的判断一样，而且明确了具体的含义，这里我列出来有兴趣的可以详细看一下：DataMark-中国地球空间数据交换格式-格网数据交换格式(CNSDTF-RAS或CNSDTF-DEM)的标志。基本部分，不可缺省。Version-该空间数据交换格式的版本号,如1.0。基本部分，不可缺省。Unit-坐标单位,K表示公里,M表示米,D表示以度为单位的经纬度,S表示以度分秒表示的经纬度(此时坐标格式为DDDMMSS.SSSS, DDD为度, MM为分, SS.SSSS为秒)。基本部分，不可缺省。A

4、lpha-方向角。基本部分，不可缺省。Compress-压缩方法。0表示不压缩,1表示游程编码。基本部分，不可缺省。Xo-左上角原点X坐标。基本部分，不可缺省。Yo-左上角原点Y坐标。基本部分，不可缺省。DX-X方向的间距。基本部分，不可缺省。DY-Y方向的间距。基本部分，不可缺省。Row-行数。基本部分，不可缺省。Col-列数。基本部分，不可缺省。HZoom-高程放大倍率。基本部分，不可缺省。设置高程的放大倍率,使高程数据可以整数存贮,如高程精度精确到厘米,高程的放大倍率为100。如果不是DEM则HZoom为1。      

5、0; 上面这些都是DEM数据的头文件信息，接下来每行记录的就是相应的DEM栅格单元值了。比如：    -99999      20738      19502      18176      16504      14833      13685  

6、    12564      11311      11307      11302      11298      11785      12149      12145     

7、; 12140      12013      11763      11635      11630      11625      11621      11617      11612 

8、0;    11608      11603      11353      11102      11097      11093      .        这些数据应该都是按照前面记录的行列数来进行组织的

9、，先不管这些。        从机器上找了一个SuperMap的示范Grid数据，在Deskpro中打开保存为ArcInfo的Grid交换格式，用记事本打开后，对比看了一下，其实文件后部分栅格单元数据部分差不多，差别就是ArInfo Grid交换格式文件头主要按下面方式记录了相关信息：ncols         319-数据列数nrows         198-数据

10、行数xllcorner     131985-数据左上角的X值yllcorner     19020-数据左上角的Y值cellsize      5-数据分辨率（栅格单元的宽高）NODATA_value  -99999-无值数据标志        对比出这样的不同，转换起来就方便了，依葫芦画瓢，把原先DEM数据的文件头信息改为了ArcInfo Grid的文件头格式（如下），ncols&#

11、160;        1328nrows         985xllcorner     39512435yllcorner     2743120cellsize      5NODATA_value  -99999并把dem文件另存为grd后缀的文件，再通过Deskpro导入这个另存后的grd文件后，

12、能导入成功了，而且看到了比较好的效果。        不过出现了新的问题，我在把上面文件临接的文件感觉没拼在一起。后来感觉应该是yllcorner的值没计算好，试着把yllcorner值2743120减去行数985*分辨率5 后再转换就没有问题（这个过程非常重要，有类似文件需要处理的朋友一定要注意）。        格网数据交换格式文件中还记录了HZoom（高程放大倍率），我们只要在Deskpro的分析菜单，栅格分析子菜单中找到栅格代数运算，把转换进来的G

13、RID数据结合这个高程放大倍率进行一个运算保存为新的数据集就可以了。前面提到需要把格网数据交换格式DEM文件转换为DEM。        在经过前面格网数据交换格式如何转换到SuperMap的数据源中介绍的方法转换后，DEM数据转换为SuperMap中的一个Grid数据集了。        SuperMap也提供了比较高级的分析功能，可以从Grid数据中提取等值线，或者提取指定等值线，我们就可以利用这个功能（可以在Deskpro的分析菜单中找到）从转入的G

14、rid数据，提取出等值线，有了等值线数据集，就可以将其转换为E00数据了（导出线数据集时选择转出文件类型为ArcInfo 交换格式）。BIL, BIP, and BSQ raster filesBand interleaved by line (BIL), band interleaved by pixel (BIP), and band sequential (BSQ) are three common methods of organizing image data for multiband images. BIL, BIP, and BSQ are not in themselves

15、image formats but are schemes for storing the actual pixel values of an image in a file. These files support the display of single and multiband images and handle black-and-white, grayscale, pseudo color, true color, and multispectral image data.The BIL, BIP, and BSQ files are binary files, and they

16、 must have an associated ASCII file header to be interpreted properly by ArcGIS. This header file contains ancillary data about the image such as the number of rows and columns in the image, if 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 data are written for row 1, all three bands of data are written for row 2, and so on, until the total number of rows in the image is reached. BIL、 BIP 、或BSQ文件可以提供有三个描述文件：一个