crop_mask_resample_reproject.py 文件源码

def get_mask(mask_file):
    mask_ds = gdal.Open(mask_file, gdal.GA_ReadOnly)
    mask_data = mask_ds.GetRasterBand(1).ReadAsArray()
    mask = mask_data != MASK_VALUES_TO_KEEP
    mask_ds = None  # close dataset
    return mask