帮我写一个Python代码实现MOD09A1影像的一天所有波段tif格式文件融合成一个tif文件

好的，我可以回答这个问题。以下是Python代码实现MOD09A1影像的一天所有波段tif格式文件融合成一个tif文件的示例代码：

import os
from osgeo import gdal

# 设置工作目录
os.chdir("MOD09A1文件夹路径")

# 获取所有tif文件名
tif_files = [f for f in os.listdir() if f.endswith(".tif")]

# 打开第一个tif文件获取影像信息
ds = gdal.Open(tif_files[])
band = ds.GetRasterBand(1)
xsize = band.XSize
ysize = band.YSize
projection = ds.GetProjection()
geotransform = ds.GetGeoTransform()

# 创建输出影像
driver = gdal.GetDriverByName("GTiff")
out_ds = driver.Create("输出影像路径", xsize, ysize, len(tif_files), band.DataType)
out_ds.SetProjection(projection)
out_ds.SetGeoTransform(geotransform)

# 逐个波段写入输出影像
for i, tif_file in enumerate(tif_files):
    ds = gdal.Open(tif_file)
    band = ds.GetRasterBand(1)
    out_band = out_ds.GetRasterBand(i+1)
    out_band.WriteArray(band.ReadAsArray())
    ds = None

# 关闭输出影像
out_ds = None

希望这个代码能够帮到你。