Google Earth Engine(GEE)——容易犯的错误7(JOIN 和Fiter)

2024-02-02 08:13:08 浏览数 (2)

JOIN vs 地图过滤器

假设您想根据时间、位置或某些元数据属性加入集合。通常,这是通过连接最有效地完成的。以下示例在 Landsat 8 和 Sentinel-2 集合之间进行时空连接:

代码语言:javascript复制
var s2 = ee.ImageCollection('COPERNICUS/S2')
    .filterBounds(ee.Geometry.Point([-2.0205, 48.647]));

var l8 = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR');

//两者进行合并保存链接
var joined = ee.Join.saveAll('landsat').apply({
  primary: s2,
  secondary: l8,
  condition: ee.Filter.and(
    ee.Filter.maxDifference({//这里影像合并的容许误差在一天内以毫秒计算
      difference: 1000 * 60 * 60 * 24, // One day in milliseconds
      leftField: 'system:time_start',//链接的方式通过这个属性波段
      rightField: 'system:time_start',
    }),
    ee.Filter.intersects({//两者相交的时候的名称
      leftField: '.geo',
      rightField: '.geo',
    })
  )
});
print(joined);

尽管您应该首先尝试连接(Export如果需要),但有时 a filter()内的 amap()也可能有效,尤其是对于非常大的集合。

代码语言:javascript复制
var s2 = ee.ImageCollection('COPERNICUS/S2')
    .filterBounds(ee.Geometry.Point([-2.0205, 48.647]));

var l8 = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR');

var mappedFilter = s2.map(function(image) {//这个map比较男鞋,里面主要包含时间,然后时间的筛选和误差,最后返回值的范围和要筛选的返回值大小等等
  var date = image.date();
  var landsat = l8
      .filterBounds(image.geometry())
      .filterDate(date.advance(-1, 'day'), date.advance(1, 'day'));//时间的差异
  // 在属性中返回具有匹配场景的输入图像。
  return image.set({
    landsat: landsat,
    size: landsat.size()
  });
}).filter(ee.Filter.gt('size', 0));
print(mappedFilter);

0 人点赞