数据类型
航天飞机雷达地形任务(SRTM),顾名思义,是一个研究任务,产生一个通用的免费数字高程模型。SRTM以 .hgt
为结尾,文件的名字解释了hgt文件的范围。比如,文件名字为 N30E11
,表示范围时30°N11°E到31°N12°E的正方形网格范围。
分为两种类型,SRTM1,SRTM3,分别对应的网格(12011201或36013601)。
资料格式
可以简单的理解为hgt文件将每一个经纬度分为了一个12011201(36013601)的网格,位置与范围如数据类型所示。
因此使用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
Original: https://blog.csdn.net/wenyusuran/article/details/123334243
Author: 文宇肃然
Title: Python使用技巧:使用Python处理SRTM(.hgt)文件
原创文章受到原创版权保护。转载请注明出处:https://www.johngo689.com/762173/
转载文章受原作者版权保护。转载请注明原作者出处!