如何在GDAL Python中使用GetPixel（）函数来处理GeoTiff图像？

Question

我们需要获取GeoTIFF图像的Pixel值来计算平均值。

在Python中，我们使用GetPixel（）来了解每个像素值，但我们应该在GDAL python中使用哪个函数。

帮帮我吧！

Answer 1

在最简单的情况下，你可以简单地使用ReadAsArray（）和numpy的平均函数：

import numpy
from osgeo import gdal

ds = gdal.Open(tiffilepath)
data = ds.ReadAsArray().astype(numpy.float32) # Type depends on the type you want for the average.
avg = numpy.average(data)
print(avg)

此统计信息也可能存在于元数据中，但可能不准确，您可能希望以不同方式计算（包括/排除无数据值，背景值等）。

Answer 2

如果要获得有限数量的像素，请使用gdallocationinfo

gdallocationinfo input.tif 250150