from osgeo import gdal
gdal.UseExceptions()
ds = gdal.Open('test_data.tif')
band = ds.GetRasterBand(1)
elevation = band.ReadAsArray()
print(elevation.shape)
print(elevation)
Read elevation using gdal python from geotiff