GEE 时序插值实现MODIS 8天无云影像代码免费分享

最新推荐文章于 2024-08-16 11:45:00 发布

CDUT_HS_2021

最新推荐文章于 2024-08-16 11:45:00 发布

阅读量2.6k

点赞数 7

分类专栏： GEE 遥感文章标签：计算机视觉人工智能

本文链接：https://blog.csdn.net/m0_51479488/article/details/122252156

版权

本文利用MOD09Q1数据，通过GEE算法实现8天每景无云遥感影像，全年共计46景。算法珍贵，与君共享。文章代码大部分来源于https://code.earthengine.google.com/7eb9f43fd3da7fdfda5bfe432a5316e7?noload=true 算法扩展性强，可以利用于Landsat系列和Sentinel-2系列，代码：

var mod_sr = ee.ImageCollection("MODIS/006/MOD09Q1"),
    point = /* color: #d63000 */ee.Geometry.Point([104.12108820077502, 30.576261504840417]),
    table = ee.FeatureCollection("China_province_2019");

var roi = table.filter(ee.Filter.eq('provinces','sichuan'))
Map.centerObject(roi,6)
function replace_mask(img, newimg, nodata) {
   
    nodata   = nodata || 0;
    var mask = img.mask();
    img = img.unmask(nodata);
    img = img.where(mask.not(), newimg);
    img = img.updateMask(img.neq(nodata));
    return img;
}
/** Interpolation not considering weights */
function addTimeBand(img) {
   
    /** make sure mask is consistent */
    var mask = img.mask();
    var time = img.metadata('system:time_start').rename("time").mask(mask);
    return img.addBands(time);
}
var getQABits = function(image, start, end) {
   
    var pattern