如何阅读NASA .hgt二进制文件

import numpy as np

def read_hgt_file(file_path):
    with open(file_path, 'rb') as file:
        data = file.read()
    
    # 每个数据点占据2个字节，因此数据点数量为文件大小除以2
    num_data_points = len(data) // 2
    
    # 使用numpy将二进制数据解析为有符号的16位整数数组
    elevations = np.frombuffer(data, dtype='>i2', count=num_data_points)
    
    return elevations

# 示例使用
file_path = 'path/to/your/file.hgt'
elevations = read_hgt_file(file_path)

# 打印前10个高程数据点
print(elevations[:10])