gdal加载影像数据的代码示例

以下是使用GDAL加载影像数据的Python代码示例：

import gdal

# 打开影像文件
ds = gdal.Open('path/to/image.tif')

# 获取影像的宽度和高度
width = ds.RasterXSize
height = ds.RasterYSize

# 获取影像的波段数
band_count = ds.RasterCount

# 获取影像的投影信息和地理变换参数
projection = ds.GetProjection()
geotransform = ds.GetGeoTransform()

# 读取影像数据
band_data = ds.GetRasterBand(1).ReadAsArray()

# 关闭影像文件
ds = None

在上面的示例中，`gdal.Open()`函数用于打开影像文件。`RasterXSize`和`RasterYSize`属性用于获取影像的宽度和高度。`RasterCount`属性用于获取影像的波段数。`GetProjection()`和`GetGeoTransform()`方法用于获取影像的投影信息和地理变换参数。`GetRasterBand()`方法用于获取指定波段的数据，这里读取的是第1个波段的数据。最后，`ReadAsArray()`方法用于读取数据并返回一个NumPy数组。读取完数据后，需要将文件关闭以释放资源。