Python使用技巧：使用Python处理SRTM(.hgt)文件
Posted 文宇肃然
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了Python使用技巧：使用Python处理SRTM(.hgt)文件相关的知识，希望对你有一定的参考价值。
数据类型
航天飞机雷达地形任务（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以上是关于Python使用技巧：使用Python处理SRTM(.hgt)文件的主要内容，如果未能解决你的问题，请参考以下文章