Python使用技巧:使用Python处理SRTM(.hgt)文件
【摘要】
数据类型
航天飞机雷达地形任务(SRTM),顾名思义,是一个研究任务,产生一个通用的免费数字高程模型。SRTM以.hgt为结尾,文件的名字解释了hgt文件的范围。比如,文件名字为N30E11,表示范围时30°N11°E到31°N12°E的正方形网格范围。
分为两种类型,SRTM1,SRTM3,分别对应的网格(1201*1201或36...
数据类型
航天飞机雷达地形任务(SRTM),顾名思义,是一个研究任务,产生一个通用的免费数字高程模型。SRTM以.hgt
为结尾,文件的名字解释了hgt文件的范围。比如,文件名字为N30E11
,表示范围时30°N11°E到31°N12°E的正方形网格范围。
分为两种类型,SRTM1,SRTM3,分别对应的网格(1201*1201或3601*3601)。
资料格式
可以简单的理解为hgt文件将每一个经纬度分为了一个1201*1201(3601*3601)的网格,位置与范围如数据类型所示。
因此使用Numpy打开hgt文件的代码下面所示:
import numpy as np
SAMPLES = 1201 # Change this to 3601 for SRTM1
def read_hgt(f_name, lat, lon):
with open(f_name, 'rb') as hgt_data:
elevations = np.fromfile(hgt_data, np.dtype('>i2'), SAMPLES * SAMPLES) \
.reshape((SAMPLES, SAMPLES))
lat_range = np.arange(lat, lat + 1 / 1200 + 1, 1 / 1200)
lon_range = np.arange(lon, lon + 1 + 1 / 1200, 1 / 1200)
retur
文章来源: wenyusuran.blog.csdn.net,作者:文宇肃然,版权归原作者所有,如需转载,请联系作者。
原文链接:wenyusuran.blog.csdn.net/article/details/123334243
【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱:
cloudbbs@huaweicloud.com
- 点赞
- 收藏
- 关注作者
评论(0)