如何使用for循环从多个棱镜.bil栅格创建栅格堆栈？

import os
from osgeo import gdal

# 设置输入文件夹路径
input_folder = "path/to/prism_bil_files"

# 获取文件夹中的所有棱镜.bil栅格文件
bil_files = [file for file in os.listdir(input_folder) if file.endswith(".bil")]

# 创建一个空的栅格堆栈
stack = []

# 遍历每个棱镜.bil栅格文件
for bil_file in bil_files:
    # 构建栅格文件的完整路径
    bil_path = os.path.join(input_folder, bil_file)
    
    # 打开栅格文件
    dataset = gdal.Open(bil_path)
    
    # 读取栅格数据
    raster_data = dataset.ReadAsArray()
    
    # 将栅格数据添加到栅格堆栈中
    stack.append(raster_data)
    
    # 关闭栅格文件
    dataset = None

# 在这里可以对栅格堆栈进行进一步处理或保存为新的栅格文件