gdal获取栅格最大最小值

使用GDAL库可以方便地获取栅格数据的最大和最小值。下面是一个示例代码：

import gdal

# 打开栅格文件
dataset = gdal.Open('path/to/raster.tif')

# 获取栅格数据的波段数量
band_count = dataset.RasterCount

# 循环遍历每个波段并获取最大最小值
for i in range(1, band_count + 1):
    band = dataset.GetRasterBand(i)
    band_min, band_max, _, _ = band.ComputeStatistics(False)
    print(f'Band {i}: Min={band_min}, Max={band_max}')

# 关闭栅格文件
dataset = None

在这个示例中，我们首先打开要处理的栅格文件。然后，我们使用`RasterCount`属性获取波段数量，并使用循环遍历每个波段。对于每个波段，我们使用`GetRasterBand`方法获取波段对象，并使用`ComputeStatistics`方法计算最小值和最大值。最后，我们输出结果并关闭文件。

需要注意的是，在使用`ComputeStatistics`方法时，我们将参数`approx_ok`设置为`False`，以确保得到准确的结果。