Python实现遥感影像批量裁剪和重采样

from osgeo import gdal, ogr
import glob
import  os
from tqdm import tqdm

vector_path = r'F:\Solar_Radiance\BSRN_SSR\Re_clip_Data\A2022311_0200\common_bounds.shp'

#输入路径
src_root = r'F:\Solar_Radiance\BSRN_SSR\MODIS\MCD43C3\TIFF\*.tif'
#输出路径
dst_root = r'F:\Solar_Radiance\BSRN_SSR\MODIS\MCD43C3\clip'
src_list = glob.glob(src_root)
# print(len(src_list))
# print(src_list,'src_list')

for src_dir in tqdm(src_list):
    dst_dir = os.path.join(dst_root,os.path.basename(src_dir))
    print('src_dir',src_dir)
    print('dst_dir',dst_dir)
    x_res = 5000
    y_res = 5000
    resample_alg = gdal.GRA_Bilinear

    # 进行重采样
    
    gdal.Warp(dst_dir, src_dir, xRes=x_res, yRes=y_res, resampleAlg=resample_alg,cutlineDSName = vector_path,cropToCutline=True)