Landsat 使用QA波段去云

GEE中不同的landsat系列影像的波段名称、QA云位数有出入,生搬硬套网上现成的QA去云很容易出错。
image
不同的影像需要设置不同的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()

image

posted @   GeoRepublic  阅读(108)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· 微软正式发布.NET 10 Preview 1:开启下一代开发框架新篇章
· C# 集成 DeepSeek 模型实现 AI 私有化(本地部署与 API 调用教程)
· spring官宣接入deepseek,真的太香了~
点击右上角即可分享
微信分享提示