Python使用技巧：使用Python处理SRTM(.hgt)文件

文宇肃然

于 2022-03-19 06:00:00 发布

阅读量2.2k

点赞数

分类专栏： python快速学习实战应用系列课程文章标签： python 开发语言后端

本文链接：https://blog.youkuaiyun.com/wenyusuran/article/details/123334243

版权

python快速学习实战应用系列课程专栏收录该内容

139 篇文章 ¥29.90 ¥99.00

订阅专栏

超级会员免费看

本文介绍了如何使用Python处理SRTM (.hgt) 文件，包括数据类型、资料格式的解析，并展示了如何利用Numpy读取文件，以及通过循环处理多个.hgt文件来绘制地形图的应用示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

数据类型

航天飞机雷达地形任务（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)