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

NASA的.hgt文件是一种二进制文件格式，用于存储地理高程数据。这些文件通常用于数字高程模型（Digital Elevation Model，DEM）。

要阅读NASA .hgt二进制文件，你需要了解文件的结构和数据格式。.hgt文件以16位有符号整数的形式存储高程数据，每个数据点占据2个字节。

以下是一个示例代码，演示如何读取和解析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])

在上述示例中，read_hgt_file函数接受一个.hgt文件的路径作为输入，并使用open函数以二进制模式打开文件。然后，它读取文件的内容，并使用numpy库的frombuffer函数将二进制数据解析为有符号的16位整数数组。

你可以将示例代码中的file_path替换为你自己的.hgt文件的路径，并使用read_hgt_file函数读取文件中的高程数据。然后，你可以根据需要对高程数据进行进一步处理和分析。