Google Earth Engine（GEE）——计算特定区域得每年特点月份得1985-2021年得ndvi影像值

我们需要统计长时间序列的ndvi值得话应该如何做一个时序分析：

这里我们传统得思路就是把这个Landsat 5 7 8和9分别进行时间边界筛选后再进行分析，但是我们如果先合并后再利用到循环中得时候，就会发现我们无法获取其中得波段进行分析，结果可以出来但都是0值。

//随便加载一个研究区
var roi = 
    /* color: #d63000 */
    /* displayProperties: [
      {
        "type": "rectangle"
      }
    ] */
    ee.Geometry.Polygon(
        [[[-118.403125, 39.70193562223249],
          [-118.403125, 39.55809120167372],
          [-118.194384765625, 39.55809120167372],
          [-118.194384765625, 39.70193562223249]]], null, false);
//Landsat5/7/8 SR数据去云函数，随后要用map()遍历函数  
function fmask(image) {  
  var cloudBitMask = (1 << 3);  
  var CloudShadowBitMask = (1 << 4);  
  var SnowBitMask = (1 << 5);  
  var qa = image.select("pixel_qa");  
  // var qa = image.select("pixel_qa");  
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0)  
                 .and(qa.bitwiseAnd(CloudS