Google Earth Engine(GEE)——sentinel-5p数据每三日SO2数据平均值

问题:

我刚接触GEE,关于您发的Google Earth Engine——Sentinel-5P (Sentinel-5P OFFL SO2)二氧化硫的使用和下载(中国区域案例分析和下载)这篇文章,我还想请教您个问题。您分享了中国区域SO2展示和下载的代码,但是这个代码里是11天的数据取平均后,镶嵌在一张图上?不知道我理解的是不是对的。我现在如果要下载三个月的数据,每三天的数据平均镶嵌在一张tif图上,(比如21年6月到8月的数据,0601-0603的数据均值在一张图上,0604-0606的数据均值在另一张图上)。请问我应该怎么修改代码呢?

 

这个问题其实很简单:只要修改日期就好了,然后把握们统计的量编程平均值即可,然后导出下载,但是这个问题有一个很严重的问题就是,因为时间短,所以很多地方都会有空值,所以建议还是逐日的数据要根据研究区内是否会产生空值来决定采取什么样的方式进行修补。

 这里代码:




// 分别定义两年的影响数据筛选
//用时间节点获取你想要的影像
var y2019 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-01","2019-06-3");
var y20191 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-04","2019-06-6");
var y20192 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-07","2019-06-9");
var y20193 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-10","2019-06-12");
var y20194 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-12","2019-06-15");
var y20195 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-16","2019-06-18");
var y20196 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-19","2019-06-20");
var y20196 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-19","2019-06-20");
var y20197 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-20","2019-06-22");
var y20198 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-23","2019-06-25");
var y20199 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-26","2019-06-28");
var y201910 = ee.ImageCollection("COPERNICUS/S5P/OFFL/L3_SO2").filterDate("2019-06-29","2019-06-30");

//获取中国边界
var countries = ee.FeatureCollection("USDOS/LSIB_SIMPLE/2017").filter(ee.Filter.eq("country_co", "CH"))

// 上色
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(countries),band_viz," 2019");
Map.addLayer(y20191.mean().select("SO2_column_number_density").clip(countries),band_viz," 20191");
Map.addLayer(y20192.mean().select("SO2_column_number_density").clip(countries),band_viz," 20192");
Map.addLayer(y20193.mean().select("SO2_column_number_density").clip(countries),band_viz," 20193");
Map.addLayer(y20194.mean().select("SO2_column_number_density").clip(countries),band_viz," 20194");
Map.addLayer(y20195.mean().select("SO2_column_number_density").clip(countries),band_viz," 20195");
Map.addLayer(y20196.mean().select("SO2_column_number_density").clip(countries),band_viz," 20196");
Map.addLayer(y20197.mean().select("SO2_column_number_density").clip(countries),band_viz," 20197");
Map.addLayer(y20198.mean().select("SO2_column_number_density").clip(countries),band_viz," 20198");
Map.addLayer(y20199.mean().select("SO2_column_number_density").clip(countries),band_viz," 20199");
Map.addLayer(y201910.mean().select("SO2_column_number_density").clip(countries),band_viz," 201910");

至于下载的化:

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

具体代码链接参照:

https://code.earthengine.google.com/5d10259b1601045adf9a7054d70c4996

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

此星光明

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

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

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

打赏作者

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

抵扣说明:

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

余额充值