SRTM数据提取（使用Python语言）

馨宝猪的驰大帅

已于 2022-08-17 10:51:36 修改

阅读量742

点赞数 1

分类专栏： GIS相关文章文章标签： python

于 2022-08-17 10:46:50 首次发布

本文链接：https://blog.csdn.net/weixin_42538763/article/details/126379389

版权

本文介绍如何利用Python从SRTM的.hgt文件中提取高程值。内容涉及从SRTM数据中获取地理信息的研究背景，重点是提供一段Python代码，该代码根据输入的经纬度自动查找.hgt文件并读取对应点的海拔高度。

摘要由CSDN通过智能技术生成

SRTM数据处理（使用Python语言）

背景

笔者在做GIS相关研究的期间，需要从SRTM数据中提取高程值（数据后缀为.hgt），在经过漫长的网络检索过程后，从https://librenepal.com/article/reading-srtm-data-with-python/处获取到灵感，略微修改，得到下面的代码，希望可以帮助到大家。代码的主要功能是：通过入参传入点的经度和纬度，自动寻找.hgt文件的路径，并通过路径提取该点的高程值。

import os
import json
import numpy as np

SAMPLES = 1201  # Change this to 3601 for SRTM1
# HGTDIR = 'hgt'  # All 'hgt' files will be kept here uncompressed


def get_elevation(lon, lat):
    hgt_file = get_file_name(lon, lat)
    if hgt_file:
        return read_elevation_from_file(hgt_file, lon, lat)
    # Treat it as data void as in SRTM documentation
    # if file is absent
    return -32768


def read_elevation_from_file(hgt_file, lon, lat):
    with open(hgt_file, 'rb') as hgt_data:
        # HGT is 16bit signed integer(i2) - big endian(>)
        elevati