How to use GetPixel() function in GDAL Python for GeoTiff Images?

Question

We need to get the Pixel value for the GeoTIFF images to calculate the mean.

In Python we use GetPixel() to know the each pixel values but which function should we use in GDAL python.

Help me Out!

Please, search the Web more thoroughly next time http://gis.stackexchange.com/questions/225370/get-individual-pixel-values-on-a-raster-image-using-gdal-and-python — mloskot, Feb 16 '17 at 13:10

score 1 · Answer 1 · answered Feb 06 '17 at 16:56

In the simplest case you can simply use ReadAsArray() and numpy's average function:

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)

This statistic also likely exists in the metadata, but may not be accurate, and you may wish to calculate it differently (including/excluding no data values, background values, etc.).

score 1 · Answer 2 · answered Jun 28 '18 at 11:40

1

If you want to get a limited number of pixels use gdallocationinfo

gdallocationinfo input.tif 250 150

answered Jun 28 '18 at 11:40

m.zohary

34
2

How to use GetPixel() function in GDAL Python for GeoTiff Images?

2 Answers2

gdallocationinfo input.tif 250 150