4

在 Google 地球引擎中,我已将 Featurecollection 作为 JSON 加载,其中包含一些多边形。我想向此 FeatureCollection 添加列,它为每个多边形和图像集合中包含的多个图像中的每一个提供了两个波段的平均值。

这是我到目前为止的代码。

//Polygons

var polygons = ee.FeatureCollection('ft:1_z8-9NMZnJie34pXG6l-3StxlcwSKSTJFfVbrdBA');

Map.addLayer(polygons);

//Date of interest

var start = ee.Date('2008-01-01');
var finish = ee.Date('2010-12-31');

//IMPORT Landsat IMAGEs
var Landsat = ee.ImageCollection('LANDSAT/LT05/C01/T1') //Landsat images
.filterBounds(polygons)
.filterDate(start,finish)
.select('B4','B3');

//Add ImageCollection to Map
Map.addLayer(Landsat);

//Map the function over the collection and display the result
print(Landsat);

// Empty Collection to fill
var ft = ee.FeatureCollection(ee.List([]))

var fill = function(img, ini) {
  // type cast
  var inift = ee.FeatureCollection(ini)

  // gets the values for the points in the current img
  var mean = img.reduceRegions({
    collection:polygons,
    reducer: ee.Reducer.mean(),
 });

 // Print the first feature, to illustrate the result.
print(ee.Feature(mean.first()).select(img.bandNames()));

  // writes the mean in each feature
  var ft2 = polygons.map(function(f){return f.set("mean", mean)})

  // merges the FeatureCollections
  return inift.merge(ft2)

  // gets the date of the img
  var date = img.date().format()

  // writes the date in each feature
  var ft3 = polygons.map(function(f){return f.set("date", date)})

  // merges the FeatureCollections
  return inift.merge(ft3)
}

// Iterates over the ImageCollection
var newft = ee.FeatureCollection(Landsat.iterate(fill, ft))

// Export
Export.table.toDrive(newft,
"anyDescription",
"anyFolder",
"test")

在控制台中,我收到一条错误消息

元素(错误)无法解码 JSON。错误:对象 '{"type":"ArgumentRef","value":null}' 的字段 'value' 缺失或为空。对象:{"type":"ArgumentRef","value":null}。

在我生成的 csv 文件中,我得到一个名为 mean 的新列,但它填充了实际值,没有实际值。

4

1 回答 1

4

没有理由在iterate()这里使用。您可以做的是嵌套的map(). 在多边形上,然后在图像上。您可以展平列表的结果列表以将其转换为单个列表,如下所示:

// compute mean band values by mapping over polygons and then over images
var results = polygons.map(function(f) {
  return images.map(function(i) {
    var mean = i.reduceRegion({
      geometry: f.geometry(),
      reducer: ee.Reducer.mean(),
    });

    return f.setMulti(mean).set({date: i.date()})
  })
})

// flatten
results = results.flatten()

脚本:https ://code.earthengine.google.com/b65a731c78f78a6f9e08300dcf552dff

同样的方法也可以用于reduceRegions()图像上的映射,然后是区域上的映射。但是,您必须映射生成的功能以设置日期。

images.filterBounds(f)如果您的功能覆盖更大的区域,也可以添加。

PS:你的表没有共享

于 2018-01-23T23:06:33.013 回答