GDAL(Geospatial Data Abstraction Library)是一个用于栅格数据操作的库，是开源地理空间基金会（Open Source Geospatial Foundation，OSGeo）的一个项目。
GDAL可以操作各种栅格地理数据格式的库，包括图像的读取、写入、转换、处理各种栅格数据格式。GDAL目前支持约100种格式的栅格数据读取，包括ERDAS Imagine、ENVI、GRASS、GeoTIFF、HDF4、HDF5、TIFF、JPEG、JPEG2000、PNG、GIF、BMP等。
GDAL的Open(filename)函数用于读栅格数据，函数返回Dataset对象。 通过Dataset对象，可以栅格数据的基本信息，如行列数、波段数、坐标转换参数等。通过Dataset对象可以返回每个波段数据（Band对象），每个波段数据的信息。Dataset对象和Band对象都可以转换成数组，通常情况下，栅格数据是基于数组进行操作。
1. Obtener información básica de los datos ráster.
Las propiedades RasterYSize, RasterXSize y RasterCount del objeto Dataset devuelven el número de filas, columnas y bandas de los datos ráster respectivamente.
de osgeo importar gdal
ds = gdal.Open("D:/img/GF2.tif")
filas = ds.RasterYSize
cols = ds.RasterXSize
bandas = ds.RasterCount
print("filas ",filas)
print("cols ", cols )
print("bands ","bands ")
2. Obtenga la referencia espacial de datos ráster
filepath="D:/img/GF2.tif"
ds = gdal.Open(filepath)
geotransform = ds .GetGeoTransform()
El método GetGeoTransform() del objeto ds devuelve los parámetros de transformación de coordenadas de los datos ráster, es decir, los parámetros de conversión de las coordenadas de fila y columna y las coordenadas espaciales.
El valor devuelto es una tupla con un total de 6 elementos. El primer y cuarto elemento son las coordenadas xey del píxel de la esquina superior izquierda, el segundo y sexto elemento son la escala en las direcciones xey, y el tercero y El quinto elemento es el ángulo de rotación en las direcciones xey.
proyección= ds .GetProjection()
El método GetProjection() del objeto ds devuelve la información del sistema de referencia espacial de los datos ráster