Landsat 使用QA波段去云
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()
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· 微软正式发布.NET 10 Preview 1:开启下一代开发框架新篇章
· C# 集成 DeepSeek 模型实现 AI 私有化(本地部署与 API 调用教程)
· spring官宣接入deepseek,真的太香了~