Google Earth Engine(GEE)—geemap无云影像合成案例分析Cloud free composite

举报
此星光明 发表于 2022/04/16 01:19:07 2022/04/16
【摘要】 import ee import geemap Map = geemap.Map() Map states = ee.FeatureCollection('TIGER/2018/States') TN = states.filter(ee.Filter.eq("NAME", "Tennessee")) Map.addLayer(TN...
import ee
import geemap

Map = geemap.Map()
Map

states = ee.FeatureCollection('TIGER/2018/States')
TN = states.filter(ee.Filter.eq("NAME", "Tennessee"))
Map.addLayer(TN, {}, "Tennessee")
years = ee.List.sequence(2013, 2020)
years.getInfo()

def yearly_image(year):
    
    start_date = ee.Date.fromYMD(year, 1, 1) 
    end_date = start_date.advance(1, "year")
    
    collection = ee.ImageCollection('LANDSAT/LC08/C01/T1') \
        .filterDate(start_date, end_date) \
        .filterBounds(TN) 
    #这里用到了一个简单合成和裁剪矢量集合的操作,因为都是对集合所以函数有些变化
    image =  ee.Algorithms.Landsat.simpleComposite(collection).clipToCollection(TN)

    return image

images = years.map(yearly_image)
vis_params = {'bands': ['B5',  'B4',  'B3'], 'max': 128}
#上面是进行了一年的操作,下面这个就是分别计算每一遍的,通过时间遍历i每一年的影像然后加载到图层!
for index in range(0, 8):
    image = ee.Image(images.get(index))
    layer_name = "Image " + str(index + 2013)
    Map.addLayer(image, vis_params, layer_name)

文章来源: blog.csdn.net,作者:此星光明2021年博客之星云计算Top3,版权归原作者所有,如需转载,请联系作者。

原文链接:blog.csdn.net/qq_31988139/article/details/119983950

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。