用GDAL读取GTIFF中的无效高程值

可以使用GDAL库中的RasterIO函数读取GTIFF格式的高程数据，并且可以通过设置读取数据时的选项来处理无效值。

以下是一个基本的示例代码，它使用了GDAL库中的RasterIO函数来读取GTIFF格式的高程数据，并且使用了GDAL库中的GDALRasterBand类来获取无效值的元数据信息，然后将无效值替换为0：

import gdal

# 打开GTIFF文件
dataset = gdal.Open('your_file.tif', gdal.GA_ReadOnly)

# 获取第一个波段
band = dataset.GetRasterBand(1)

# 获取无效值元数据信息
nodata = band.GetNoDataValue()

# 读取高程数据
elevation = band.ReadAsArray()

# 将无效值替换为0
elevation[elevation == nodata] = 0

# 关闭数据集
dataset = None

在上面的代码中，我们首先打开GTIFF文件，然后获取第一个波段。接下来，我们获取无效值的元数据信息，然后使用ReadAsArray函数读取高程数据。最后，我们使用numpy库将无效值替换为0，并且在完成数据处理后关闭数据集。