USAGE OF GDAL RASTER API
Import gdal
from osgeo import gdal
Open the file
Dataset=gdal.Open(filename)
Getting dataset information
Xsize=dataset.RasterXsize
Ysize=dataset.RasterYSize
Count=dataset.RasterCount
Geotransform=dataset.GetGeoTransform()
Pro=dataset.GetProjection()
Fetching a raster band
band=dataset.GetRasterBand(1)
type=gdal.GetDataTypeName(band.DataType)
min=band.GetMinimum()
max=band.GetMaximum()
if not min or not max:
(min,max)=band.ComputeRasterMinMax(True)
Reading raster data
scanline=band.ReadRaster(xoff=0,yoff=0,xsize=band.XSize,ysize=1,buf_xsize=band.XSize,buf_ysize=1,buf_type=gdal.GDT_Float32)
import struct
tuple_of_floats=struct.unpack(‘f’*band.XSize,scanline)
CreatCopy
driver=src_ds.GetDriver