Google Earth Engine(GEE)——Sentinel-5p数据的正确展示和下载

博客讨论了在使用Sentinel-5P数据时遇到的问题,即使用filterBounds()方法筛选福建省数据无效,导致返回全球数据。问题出在使用中文命名上,解决方案是避免使用中文并简化代码。提供了修正后的代码示例,通过加载中国区域矢量边界,筛选指定日期的 Sentinel-5P 数据,并展示了如何进行数据下载和地图显示。
摘要由CSDN通过智能技术生成

问题:

各位,请问一下sentinel-5p数据筛选的时候,为什么filterBounds()方法不起作用呢?(我设置范围为福建省,但是他还是把全球都是筛选出来了)

来看一下他的原始代码:

var O3 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_O3"),
var China = ee.FeatureCollection("users/Zhouxiaokang/Shape/China_Province_Zone"),
var Fujian = ee.FeatureCollection("users/Zhouxiaokang/Fujian");


var star_date=ee.Date.fromYMD(2019,1,1);
var end_date=ee.Date.fromYMD(2020,12,31);

var Fujian=China.filter(ee.Filter.eq("省","福建省"));
var S5p_Fujian=O3.filter(ee.Filter.bounds(Fujian))
                  .filterDate(star_date,end_date);

print(S5p_Fujian.limit(100));
Map.centerObject(Fujian,6);
Map.addLayer(Fujian);
Map.addLayer(S5p_Fujian.first());


首先来看问题,影像和适量添加的都没有问题,时间时日筛选也没有问题,当然,有更简单的方法,一会看我的代码即可。这里的重点是没有我们筛选的时候不能有中文的出现,所以错误出现在这里,这是第一个,我们一般会以英文进行命名,而且不必要向上面代码这么麻烦,直接一个 filterBounds()来解决即可,然后打印这里如果不是为了获取影像的数量不建议打印,这一步有点多余,不过初学者来说可以坚持这么做,毕竟可以发现问题的所在之处。

filterBounds(geometry)

Shortcut to filter a collection by intersection with geometry. Items in the collection with a footprint that fails to intersect the given geometry will be excluded.

This is equivalent to this.filter(ee.Filter.bounds(...)).

Caution: providing a large or complex collection as input can result in poor performance. Collating the geometry of collections does not scale well, use the smallest collection (or geometry) that is required to achieve the desired outcome.

Returns the filtered collection.

Arguments:

this:collection (Collection):

The Collection instance.

geometry (ComputedObject|FeatureCollection|Geometry):

The geometry, feature or collection to intersect with.

Returns: Collection

改正

今天看到有个朋友说无法正确的使用filterBounds()这个函数,在对于sentinel-5P的数据中,对于中国区域没法正确使用,我很是纳闷,随后我就直接自己赶快从之前的博客中找来我的代码,然后进行了一番改写之后实现了中国区域的sentinel-5P数据的展示和下载。

//加载中国区域的矢量边界
//如果你是省的矢量也直接夹在即可
var table2 = ee.FeatureCollection("users/bqt2000204051/CHINA");

 
// 分别定义两年的影响数据筛选
//用时间节点获取你想要的影像
var y2019 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-01","2019-06-3").filterBounds(table2);

// 上色
var band_viz = {
  min: 0.0,
  max: 0.0005,
  palette: ['black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']
};
//图层加载
Map.addLayer(y2019.mean().select("SO2_column_number_density").clip(table2),band_viz," 2019");

//分别利用已经定义的影像名称来下载即可
Export.image.toDrive({
//这里的mean就代表求均值
    image: y2019.mean().select("SO2_column_number_density").clip(table2),
    region:table2,
    scale:1000,
    description: "CHINA_so2_1Km",
    folder: 'CHINA_so2_1Km',
  });

最后结果

 

  • 4
    点赞
  • 13
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值