如何从具有不同范围和分辨率的1个栅格中减去3个栅格

from osgeo import gdal

def resample_raster(src_raster, dst_raster, target_resolution):
    gdal.ReprojectImage(src_raster, dst_raster, src_raster.GetProjection(), dst_raster.GetProjection(), gdal.GRA_Bilinear)

def align_rasters(raster_list):
    # 获取所有栅格的最小边界框
    min_x = min([r.GetExtent()[0] for r in raster_list])
    max_y = max([r.GetExtent()[3] for r in raster_list])
    max_x = max([r.GetExtent()[1] for r in raster_list])
    min_y = min([r.GetExtent()[2] for r in raster_list])
    
    # 裁剪或扩展所有栅格到这个边界框
    for r in raster_list:
        r.SetProjection(...)  # 设置相同的投影
        r.SetGeoTransform(...)  # 设置相同的地理变换参数

import numpy as np

def subtract_rasters(raster1, raster2, raster3, raster4):
    data1 = raster1.ReadAsArray()
    data2 = raster2.ReadAsArray()
    data3 = raster3.ReadAsArray()
    data4 = raster4.ReadAsArray()
    
    result = data1 - data2 - data3 - data4
    return result