登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
前言:近期因研究需要,对区域地形高程数据进行处理,将SRTM地形高程数据多次转换处理成micaps第四类格式,并把资料地址、处理经验和转换所写的小程序以及所得的结果在这里和大家共享,欢迎提出不足。另外,由于转换用的是csharp程序所以就把此帖子发在了C#/ASP.NET板块里了,增加增加C#的人气!
言归正传,先看几张结果图
0.05°分辨率地形高程数据micaps显示
0.05°分辨率地形高程数据
0.25°分辨率地形高程数据显示
地形高程数据数值显示
0.25°分辨率地形高程数据、数值
格栅数据转micaps程序
转换程序
处理说明=================================
所用工具Global Mapper、Csharp、Micaps3.2版本
一、数据下载 数据下载地址: 或者:(2017.03 use this address) 简单注册即可下载,数据路径: 数据资源->免费数据->DEM数字高程数据->GDEMV2 30M分辨率数字高程数据
下载所需经纬度区域所有数据块1°×1°,放同一目录下。 数据信息 For 1x1 文件名中包含了该单元左下角位置的信息;比如文件名类似N37W105 的数据文件,其单元左下角的坐标为 (37N,105W); 二、数据Global Mapper预处理 1、打开文件,[File->OpenData Files],等待读取加载所有文件…….. 2、导出地形高度数据 [File->Export->Export Elevation GridFormat…] 导出数据 根据输出数据需要设置X-axis/Y-axis方向分辨率。 得到栅格高程数据结果文件 ncols 97 nrows 81 xllcorner 104.8750000000 yllcorner 29.8750000000 cellsize 0.25 nodata_value -9999.0 1101 865 775 643 822 708 1032 1180 793 634 638 733 975 724 1312 1617 1236 1411 1560 1308 1359 1419 1384 1658 1731 1434 1537 1285 1177 1349 1244 851 1158 925 1049 866 780 764 780 802 697 644 584 581 695 634 848 700 723 770 771 707 675 741 658 532 520 521 546 713 756 770 856 859 883 1003 802 891 792 776 959 555 917 495 427 424 434 383 383 426 262 279 370 390 533 551 546 441 355 311 130 125 162 174 176 184 166 712 848 966 801 670 786 892 685 792 900 685 815 1015 1051 1242 1565 1259 1373 1447 1464 2144 1474 1300 1477 1621 1386 1498 1159 1227 1147 1056 803 1088 929 1002 940 866 859 773 691 667 631 647 633 651 709 684 666 839 …………………… 三、编程再处理 读取格栅数据文件,把数据坐标方向作调整,输出为micaps-diamond4类数据格式文件;缺测值-9999数据处理(一般为海洋区域,海拔设为0即可)。
四、高程数据介绍(网络上搜集哒) SRTM数据主要是由美国太空总署(NASA)和国防部国家测绘局(NIMA)联合测量的,SRTM的全称是ShuttleRadar Topography Mission,即航天飞机雷达地形测绘使命。获取的雷达影像数据经过两年多的处理,制成了数字地形高程模型,该测量数据覆盖中国全境。SRTM数据每经纬度方格提供一个文件,精度有1arc-second和3arc-seconds两种,称作SRTM1和SRTM3,或者称作30M和90M数据,SRTM1的文件里面包含3601*3601个采样点的高度数据,SRTM3的文件里面包含1201*1201个采样点的高度数据。目前能够免费获取中国境内的RTM3文件,是90米的数据,每个90米的数据点是由 9个30米的数据点算术平均得来的。 SRTM数据下载的地址是: ftp://e0mss21u.ecs.nasa.gov/srtm/。我国境内的数据在Eurasia目录下,每经纬度方格一个文件,文件命名方法是X1X2X3X4.hgt.zip,X1是N或S表示南北,X2是下方纬度数,X3是E或W表示东西,X4是左方经度数。 我们在进行地学研究时有时会用到这种免费下载但分辨率适中的数据,可是SRTM数字高程模型的格式是hgt,我们怎么把其转换成ArcGIS支持的格式以供我们查看呢?美国地质调查局的工程师们用AML和地图代数语言混合编程写了一个小程序,专门用来将SRTM数据由hgt转换成GRID格式。把下列代码拷贝到记事本上,另存为srtmgrid.aml。
最后重要的:-------------------------------------------- |