如何在Google Earth Engine中的时间序列图像中计算区域统计？

Question

我从Landsat图像计算了2018年1月至2018年12月的NDVI，给出了21个NDVI波段。我有一个带有三个多边形的shapefile。我想要三个与给定shapefile中的多边形相对应的NDVI时间序列值。

进口（2个条目） var L8：ImageCollection“ USGS Landsat 8 Collection 1 Tier 1和实时数据TOA反射率” var roi：表users / user_name / sample_shapefile

// NDVI function

  function addNDVI(image) {
    users/user_name/sample_shapefilevar ndvi = image.normalizedDifference(['B5', 'B4']);
      return image.addBands(ndvi);`
       }
// Date range 2018 jan to dec

   var filtered = L8.filterDate('2018-01-01', '2018-12-31')       
.filterBounds(roi);` 

var with_ndvi = filtered.map(addNDVI);

    var ndvi = with_ndvi.select('nd');

var dt = ndvi.reduceRegion({
  reducer: ee.Reducer.mean(),
  geometry: roi.geometry(),
  scale: 30,
});
print(dt)

我想为给定shapefile中的每个多边形接收三个时间序列NDVI值。