写一个python读取遥感图像监督分类的代码

以下是一个Python读取遥感图像监督分类的代码示例：

from osgeo import gdal
import numpy as np

# 读取遥感图像
img_path = 'path/to/image.tif'
img_dataset = gdal.Open(img_path, gdal.GA_ReadOnly)

# 获取图像的元数据
img_cols = img_dataset.RasterXSize
img_rows = img_dataset.RasterYSize
img_bands = img_dataset.RasterCount
img_projection = img_dataset.GetProjection()
img_geotransform = img_dataset.GetGeoTransform()

# 将图像转换为numpy数组
img_array = np.zeros((img_rows, img_cols, img_bands))
for b in range(img_bands):
    band = img_dataset.GetRasterBand(b + 1)
    img_array[:, :, b] = band.ReadAsArray()

# 读取分类标签
label_path = 'path/to/label.tif'
label_dataset = gdal.Open(label_path, gdal.GA_ReadOnly)

# 将标签转换为numpy数组
label_array = label_dataset.ReadAsArray()

# 关闭数据集
img_dataset = None
label_dataset = None

以上代码使用了GDAL库进行读取，可以读取多种遥感图像格式。读取后，将图像转换为Numpy数组，方便进行数据处理和分析。同时，也读取了分类标签图像，以便进行监督分类。