爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 22649|回复: 15

[原创]高分辨率高程地形数据处理成micaps第四类

[复制链接]

新浪微博达人勋

发表于 2017-3-17 09:35:33 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
前言:近期因研究需要,对区域地形高程数据进行处理,将SRTM地形高程数据多次转换处理成micaps第四类格式,并把资料地址、处理经验和转换所写的小程序以及所得的结果在这里和大家共享,欢迎提出不足。另外,由于转换用的是csharp程序所以就把此帖子发在了C#/ASP.NET板块里了,增加增加C#的人气


言归正传,先看几张结果图

0.05°分辨率地形高程数据micaps显示

0.05°分辨率地形高程数据micaps显示

0.05°分辨率地形高程数据

0.25°分辨率地形高程数据显示

0.25°分辨率地形高程数据显示

地形高程数据数值显示

地形高程数据数值显示

0.25°分辨率地形高程数据、数值

格栅数据转micaps程序

格栅数据转micaps程序

转换程序


处理说明=================================
所用工具Global MapperCsharpMicaps3.2版本
一、数据下载
数据下载地址:
或者:(2017.03 use this address)
http://www.gscloud.cn/(测试正常)
简单注册即可下载,数据路径:
数据资源->免费数据->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。

最后重要的:--------------------------------------------
dem_125.asc (117.42 KB, 下载次数: 19)

评分

参与人数 1金钱 +20 贡献 +5 收起 理由
mofangbao + 20 + 5 赞一个!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-17 10:17:34 | 显示全部楼层
补充一下:附件提供的数据区域范围为30N-50N,105E-129E,其它区域未下载
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-25 13:20:36 | 显示全部楼层
共享一下源程序
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-29 10:16:24 | 显示全部楼层
方法很好呀
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-29 10:50:01 | 显示全部楼层
谢谢楼主分享,楼主辛苦了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-31 10:07:22 | 显示全部楼层
{:eb502:}{:eb502:}{:eb502:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-11-15 14:23:52 | 显示全部楼层
{:eb500{:eb505:}{:eb512:}{:eb500{:eb505:}{:eb512:}{:eb500
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-11-15 14:23:58 | 显示全部楼层
{:eb500{:eb505:}{:eb512:}{:eb500{:eb505:}{:eb512:}{:eb500
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-11-15 14:24:03 | 显示全部楼层
{:eb500{:eb505:}{:eb512:}{:eb500{:eb505:}{:eb512:}{:eb500
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-5-30 15:33:53 | 显示全部楼层
谢谢分享,十分给力!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表