gdal_raster_open_read_projection_geotransform_numpy_array.py

1from osgeo import gdal
2import numpy as np
3
4# Open the dataset
5dataset = gdal.Open("example.tif")
6
7# Get basic metadata
8print(f"Driver: {dataset.GetDriver().ShortName}/{dataset.GetDriver().LongName}")
9print(f"Size is {dataset.RasterXSize} x {dataset.RasterYSize} x {dataset.RasterCount}")
10print(f"Projection is {dataset.GetProjection()}")
11
12# Get geotransform (Top Left X, W-E pixel resolution, rotation, Top Left Y, rotation, N-S pixel resolution)
13geotransform = dataset.GetGeoTransform()
14if geotransform:
15    print(f"Origin = ({geotransform[0]}, {geotransform[3]})")
16    print(f"Pixel Size = ({geotransform[1]}, {geotransform[5]})")
17
18# Read the first band
19band = dataset.GetRasterBand(1)
20print(f"Band Type={gdal.GetDataTypeName(band.DataType)}")
21
22# Read data as a numpy array
23data = band.ReadAsArray(0, 0, dataset.RasterXSize, dataset.RasterYSize)
24
25# Clean up
26dataset = None