DEM高程数据处理教程.doc

DEM高程数据DEM高程数据包括两个部分：ASTER GDEM30米分辨率高程数据和SRTM90米分辨率高程数据。ASTER GDEM数据来源于NASA，数据覆盖范围为北纬83到南纬83之间的所有陆地区域，时间范围为2000年前后；SRTM数据来源于CIAT，数据覆盖范围为北纬60至南纬60之间的所有陆地区域，时间范围为2000年前后。ASTER GDEM 30米分辨率高程数据本数据集利用ASTER GDEM第一版本（V1）的数据进行加工得来，是全球空间分辨率为30米的数字高程数据产品。由于云覆盖，边界堆叠产生的直线，坑，隆起，大坝或其他异常等的影响，ASTER GDEM第一版本原始数据局部地区存在异常，所以由ASTER GDEMV1加工的数字高程数据产品存在个别区域的数据异常现象，可以和全球90米分辨率数字高程数据产品互相补充使用。ASTER GDEM数据采用UTM/WGS84投影，数据格式为IMG栅格影像，数据的值域范围为-152-8806米之间，比例尺为1:25万，其垂直精度20米，水平精度30米。数据命名规则：ASTER GDEM基本的单元按1度X1度分片。每个GDEM数据包有两个文件，一个数据高程文件和一个质量评估（QA）文件。每个文件的命名是根据影像几何中心左下角的经纬度产生。例如，ASTGTM_N29E091代表左下角坐标是北纬29度，东经91度。ASTGTM_N29E091_dem和ASTGTM_N29E091_num对应的分别是高程数据和质量控制数据。SRTM 90米分辨率高程数据SRTM（ShuttleRadarTopographyMission）90米分辨率高程数据由美国太空总署（NASA）和国防部国家测绘局（NIMA）联合测量。2000年2月11日，美国发射的“奋进”号航天飞机上搭载SRTM系统，共计进行了222小时23分钟的数据采集工作，获取北纬60度至南纬60度之间总面积超过1.19亿平方公里的雷达影像数据，覆盖地球80%以上的陆地表面。SRTM系统获取的雷达影像的数据量约为9.8万亿字节，经过两年多的数据处理，制成了数字地形高程模型（DEM），即现在的SRTM地形产品数据。SRTM因插值算法不同，存在不同版本，本平台发布数据为V4.1版本。SRTM数据采用WGS84椭球投影，使用16位的数值表示高程数值的（+-32767米），空数据用-32726表示。数据格式为IMG栅格影像，数据的值域范围为-12000-9000米之间，水平精度20米，高程精度16米。SRTM的数据组织方式为：每5度经纬度方格划分一个文件，共分为24行（-60至60度）和72列（-180至180度）。文件命名规则为srtm_XX_YY.zip，XX表示列数（01-72），YY表示行数（01-24）。示意图如下：高程数据处理方法第一次使用DEM高程数据的朋友常常遇到这个问题，IMG是压缩包么？怎么不能解压呢？为什么我打开之后数据是灰色的呢？明明是平原地区，为什么显示的高程范围却在-32767-32767之间呢？为什么展示图里是五颜六色的，而我打开的却是灰色影像呢？首先IMG不是压缩包，“.img”作为一种栅格影像格式，可以直接在ArcMap、ENVI、ERDAS等遥感软件中打开使用，无需解压。其次，怎么去除高程影像中的空值（如+-32767），让它在一个正常的范围内显示呢？小编这里以TIF格式的DEM高程影像为例（IMG的处理方式同样），一步步带大家来操作。1. 在ARCMAP里打开一幅DEM高程数据（ADD DATA），可以从左边看到其显示的数据范围是-32767-32726，右侧为灰色影像。2. 在ArcMap里打开Spatial analyze工具，选择Raster Calculator，设置DEM高程数据值为0并进行计算（点击Evaluate按钮），页面如下：3. 计算之后结果显示为Calculation，页面如下图：4. 再次打开Spatial analyze工具，选择Reclassify，选中chang missing values to Nodata打钩, 点击OK：5. 然后打开Raster Calculator面板，将上一步计算出来的Reclass of Calculation与原始的DEM 高程数据相乘, 点击Evaluate，界面如下：6. 计算结果，如图：7. 鼠标右键点击Calculation2，打开Properties, 选择Symbology选项，在左侧选中Stretched, 点击确定。如下图：8. 在这里，我们看到，经过处理，高程影像的值域范围已经变为正常的海拔范围（105-4336）了：9. 我们还可以给它设一个漂亮的显示颜色，简单的双击左侧的灰度条即可，到这一步，DEM数据的预处理就做完了。DEM高程数据