繁体   English   中英

减少图像集合 Google Earth Engine

[英]Reduction over an image collection Google Earth Engine

我正在使用 CHIRPS 数据集(每日降水量)来推导特定时间范围内的平均、中值、最小和最大降水量。 然后,我想提取点 shapefile 中包含的特定位置的值并将结果保存在表中。 该脚本似乎工作,但输出表只有零 (0) 作为 4 个变量的值。 请看下面的脚本

var lng = 65.64; 
var lat = 34.35;
var point = ee.Geometry.Point(lat, lng); 
//var aoi = point.buffer(100000); // Create an area (1km buffer around point)
var country = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017')
.filter(ee.Filter.eq('country_co', 'AF'));
var aoi = country;

Map.setCenter(lng, lat, 5); // Center the map on this location, zoom level 10

var start = '2018-02-15'; // initial date of the image collection
var end = '2018-07-15'; //final date of the image collection

var p1 = ee.Geometry.Point([69.78086, 34.65411])
var p2 = ee.Geometry.Point([61.82234, 30.66048])
var table = ee.FeatureCollection(ee.List([ee.Feature(p1),ee.Feature(p2)]))

var AddPrMean = function(image) {
  var PrMean = image.reduce(ee.Reducer.mean()).rename('PrMean');
  return image.addBands(PrMean);
};

var AddPrMedian = function(image) {
  var PrMedian = image.reduce(ee.Reducer.median()).rename('PrMedian');
  return image.addBands(PrMedian);
};

var AddPrMin = function(image) {
  var PrMin = image.reduce(ee.Reducer.min()).rename('PrMin');
  return image.addBands(PrMin);
};

var AddPrMax = function(image) {
  var PrMax = image.reduce(ee.Reducer.max()).rename('PrMax');
  return image.addBands(PrMax);
};

var dataset = ee.ImageCollection('UCSB-CHG/CHIRPS/DAILY')
    .filterDate(start, end)
    .filterBounds(aoi)
    .map(AddPrMean)
    .map(AddPrMedian)
    .map(AddPrMin)
    .map(AddPrMax);
    
var composites = dataset.select(['PrMean','PrMedian','PrMin','PrMax']).first();   

var YieldLocations = ee.FeatureCollection(table);

var YPrec = composites.reduceRegions(YieldLocations, ee.Reducer.max(), 1);

print(YPrec); ``` 

我找到了一个似乎有效的解决方案,但是,结果中位数始终为 0。

var lat = 34.35;
var point = ee.Geometry.Point(lat, lng); 
//var aoi = point.buffer(100000); // Create an area (1km buffer around point)
var country = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017')
.filter(ee.Filter.eq('country_co', 'AF'));
var aoi = country;

Map.setCenter(lng, lat, 5); // Center the map on this location, zoom level 10

var start = '2018-02-15'; // initial date of the image collection
var end = '2018-07-15'; //final date of the image collection

var p1 = ee.Geometry.Point([69.78086, 34.65411])
var p2 = ee.Geometry.Point([61.82234, 30.66048])
var table = ee.FeatureCollection(ee.List([ee.Feature(p1),ee.Feature(p2)]))

var dataset = ee.ImageCollection('UCSB-CHG/CHIRPS/DAILY')
    .filterDate(start, end)
    .filterBounds(aoi);
    
var PrMean = dataset.mean().rename('PrMean');

var PrMedian = dataset.median().rename('PrMedian');

var PrMin = dataset.min().rename('PrMin');

var PrMax = dataset.max().rename('PrMax');

var composites = PrMean
      .addBands(PrMedian)
      .addBands(PrMin)
      .addBands(PrMax);
  



var YieldLocations = ee.FeatureCollection(table);

var YPrec = composites.reduceRegions(YieldLocations, ee.Reducer.max(), 1);

print(YPrec);```

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM