推荐订阅源

GEE中不同的landsat系列影像的波段名称、QA云位数有出入，生搬硬套网上现成的QA去云很容易出错。

不同的影像需要设置不同的QA_PIXEL波段名称和比特位数，都在GEE的数据波段说明文档里。当前版本似乎都是3-云、4-云阴影。了解了原理，后面遇到陌生影像也能处理。

def rmCloud(image):
    # for LANDSAT/LE07/C02/T1_L2
    # 选择QA_PIXEL波段
    qa = image.select('QA_PIXEL')
    
    # 创建云掩膜，QA_PIXEL波段的第3位表示云
    cloud_mask = qa.bitwiseAnd(1 << 3).eq(0)
    
    # 创建云影掩膜，QA_PIXEL波段的第4位表示云影
    cloud_shadow_mask = qa.bitwiseAnd(1 << 4).eq(0)
    
    # 结合云和云影掩膜
    mask_all = cloud_mask.And(cloud_shadow_mask)
    
    # 应用掩膜，去除云和云影
    return image.updateMask(mask_all)

landsat7 = ee.ImageCollection("LANDSAT/LE07/C02/T1_L2")
landsat7.filterDate(i_date, f_date).map(rmCloud).median()