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

简介: 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


目录
打赏
0
0
0
0
213
分享
相关文章
如何用Google Earth Engine快速、大量下载遥感影像数据?
【2月更文挑战第9天】本文介绍在谷歌地球引擎(Google Earth Engine,GEE)中,批量下载指定时间范围、空间范围的遥感影像数据(包括Landsat、Sentinel等)的方法~
3280 1
如何用Google Earth Engine快速、大量下载遥感影像数据?
GEE时序——利用sentinel-2(哨兵-2)数据进行地表物候学分析(时间序列平滑法估算和非平滑算法代码)
GEE时序——利用sentinel-2(哨兵-2)数据进行地表物候学分析(时间序列平滑法估算和非平滑算法代码)
1004 3
|
12月前
|
API
GEE案例分析——利用sentinel-3数据计算空气污染指数(Air Pollution Index,简称API)
GEE案例分析——利用sentinel-3数据计算空气污染指数(Air Pollution Index,简称API)
316 0
基于Google Earth Engine云平台构建的多源遥感数据森林地上生物量AGB估算模型含生物量模型应用APP
基于Google Earth Engine云平台构建的多源遥感数据森林地上生物量AGB估算模型含生物量模型应用APP
363 0
|
12月前
GEE——Google dynamic world中在影像导出过程中无法完全导出较大面积影像的解决方案(投影的转换)EPSG:32630和EPSG:4326的区别
GEE——Google dynamic world中在影像导出过程中无法完全导出较大面积影像的解决方案(投影的转换)EPSG:32630和EPSG:4326的区别
199 0
R语言指数平滑法holt-winters分析谷歌Google Analytics博客用户访问时间序列数据
R语言指数平滑法holt-winters分析谷歌Google Analytics博客用户访问时间序列数据
GEE 案例——如何计算sentinel-2中每一个单景影像的波段的DN值并绘制直方图
GEE 案例——如何计算sentinel-2中每一个单景影像的波段的DN值并绘制直方图
177 5
GEE Colab——如何从本地/Google云盘/Google Cloud Storage (GCS)上传和下载
GEE Colab——如何从本地/Google云盘/Google Cloud Storage (GCS)上传和下载
540 4
GEE Colab——初学者福音快速入门 Google Colab(Colaboratory)
GEE Colab——初学者福音快速入门 Google Colab(Colaboratory)
363 3
GEE数据集——哨兵2号Sentinel-2 云概率数据集
GEE数据集——哨兵2号Sentinel-2 云概率数据集
495 2
AI助理

你好,我是AI助理

可以解答问题、推荐解决方案等