python gdal批量读入遥感影像

使用 GDAL 库可以方便地批量读入遥感影像。下面是一个简单的 Python 代码示例：

import os
from osgeo import gdal

# 遥感影像所在文件夹路径
dir_path = '/path/to/remote/sensing/images'

# 获取文件夹下所有影像文件路径
image_paths = [os.path.join(dir_path, f) for f in os.listdir(dir_path) if f.endswith('.tif')]

# 循环遍历每个影像文件并读取
for image_path in image_paths:
    # 打开影像文件
    dataset = gdal.Open(image_path)

    # 获取影像信息
    cols = dataset.RasterXSize
    rows = dataset.RasterYSize
    bands = dataset.RasterCount

    # 读取影像数据
    data = []
    for i in range(1, bands + 1):
        band = dataset.GetRasterBand(i)
        band_data = band.ReadAsArray(0, 0, cols, rows)
        data.append(band_data)

    # 关闭数据集
    dataset = None

    # 处理影像数据
    # ...

以上代码首先获取遥感影像所在文件夹路径，并通过 `os.listdir` 函数获取文件夹下所有 `.tif` 文件的路径。然后循环遍历每个影像文件，打开文件并获取影像的信息（列数、行数和波段数），接着读取每个波段的数据，并将其存储在 `data` 列表中。最后关闭数据集，并对数据进行处理。