gdal将三个矩阵组合显示

如果您有三个矩阵，您可以使用numpy库将它们组合成一个多波段数组，然后使用gdal库将其转换为可视化格式（如TIFF、PNG、JPEG等）。

以下是一个示例代码：

import numpy as np
from osgeo import gdal

# 读取三个矩阵
matrix1 = np.loadtxt('matrix1.txt')
matrix2 = np.loadtxt('matrix2.txt')
matrix3 = np.loadtxt('matrix3.txt')

# 将三个矩阵组合成一个三维数组
stacked = np.stack([matrix1, matrix2, matrix3], axis=2)

# 将数组保存为TIFF格式的文件
driver = gdal.GetDriverByName('GTiff')
dataset = driver.Create('output.tif', stacked.shape[1], stacked.shape[0], stacked.shape[2], gdal.GDT_Float32)
for i in range(stacked.shape[2]):
    dataset.GetRasterBand(i+1).WriteArray(stacked[:,:,i])
dataset.FlushCache()

该代码将三个矩阵组合成一个三维数组，然后使用gdal库将其保存为TIFF格式的文件。您可以根据需要修改代码以适应不同的输入和输出格式。