Google Earth Engine——Landsat8谐波模型（NDVI）的时序分析案例

代码：

// 加载LandsatTOA大气层顶反射数据
var landsatCollection = ee.ImageCollection('LANDSAT/LC08/C01/T1_TOA');

// 设置感兴趣的区域
var roi = ee.Geometry.Point([-121.14, 37.98]);

// 我们要建模的因变量。
var dependent = 'NDVI';

// 建立模型的每年的周期数。
var harmonics = 1;

// 编制一个谐波频率的清单来建立模型。
// 这些也可以作为波段名称的后缀。
var harmonicFrequencies = ee.List.sequence(1, harmonics);

// 获取谐波术语的带状名称序列的函数。
var constructBandNames = function(base, list) {
  return ee.List(list).map(function(i) {
    return ee.String(base).cat(ee.Number(i).int());
  });
};


var cosNames = constructBandNames('cos_', harmonicFrequencies);
var sinNames = constructBandNames('sin_', harmonicFrequencies);

// 独立变量。
var independents = ee.List(['constant', 't'])
  .cat(cosNames).cat(sinNames);

// Landsat8 简单去云函数
var maskClouds = function(image) {
  var score = ee.Algorithms.Landsat.simpleCloudScore(image).select(