python gdal.Warp tif裁剪、拼接、重采样_gdal.warp进行影像拼接

参考https://www.jianshu.com/p/821c741ff169

拼接

def rasterwarp(base_file,inpath,outpath):
    if not os.path.exists(outpath):
        os.mkdir(outpath)
    out_file = os.path.join(outpath,'mcd12q1.tif')
    inputrasfile1 = gdal.Open(base_file, gdal.GA_ReadOnly)  # 第一幅影像
    inputProj1 = inputrasfile1.GetProjection()
    options = gdal.WarpOptions(srcSRS=inputProj1, dstSRS=inputProj1, format='GTiff', resampleAlg=gdalconst.GRA_Max)
    rasterlist = []
    filelist = pathlib.Path(inpath).glob('*.tif')
    for rasterfile in filelist:
        rasterfile = str(rasterfile)
        inputrasfile2 = gdal.Open(rasterfile, gdal.GA_ReadOnly)  # 第二幅影像
        rasterl