GEE开发之降雨(CHIRPS)数据获取和分析

news/2024/12/2 16:32:38/

GEE开发之降雨CHIRPS数据获取和分析

  • 1.数据介绍
  • 2.初识CHIRPS
    • 2.1 代码一
    • 2.2 代码二
  • 3.逐日数据分析和获取
  • 4.逐月数据分析和获取
    • 4.1 代码一
    • 4.2 代码二(简洁)
  • 5.逐年数据分析和获取
    • 5.1 代码一
    • 5.2 代码二(简洁)

前言:主要获取和分析UCSB-CHG/CHIRPS/DAILY的日数据、月数据和年数据。


1.数据介绍

  • 数据集: CHIRPS Daily: Climate Hazards Group InfraRed Precipitation With Station Data (Version 2.0 Final)
  • 数据说明: Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) 是一个记录了1981年到当前全球降雨量的数据集。CHIRPS将0.05°分辨率的卫星图像与in-situ站点数据结合,创建网格化的降雨时间序列,用于趋势分析和季节性干旱监测。

2.初识CHIRPS

2.1 代码一

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
var dataset = ee.ImageCollection('UCSB-CHG/CHIRPS/DAILY')
.filterDate('2020-01-01', '2020-12-31')
.filterBounds(geometry)
.select('precipitation');
print(dataset)

在这里插入图片描述

2.2 代码二

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,7);
var dataset = ee.ImageCollection('UCSB-CHG/CHIRPS/DAILY')
.filterDate('2020-01-01', '2020-12-31')
.filterBounds(geometry)
.select('precipitation');
var precipitationVis = {min: 1.0,max: 17.0,palette: ['001137', '0aab1e', 'e7eb05', 'ff4a2d', 'e90000'],
};
Map.addLayer(dataset.mean().clip(geometry), precipitationVis, 'pre');

在这里插入图片描述

3.逐日数据分析和获取

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,7);
var dataset = ee.ImageCollection('UCSB-CHG/CHIRPS/DAILY')
.filterDate('2020-05-01', '2020-05-31')
.filterBounds(geometry)
.select('precipitation');
print(ui.Chart.image.series(dataset, geometry, ee.Reducer.mean(), 1000));
function exportImageCollection(imgCol) {var indexList = imgCol.reduceColumns(ee.Reducer.toList(), ["system:index"]).get("list");indexList.evaluate(function(indexs) {for (var i=0; i<indexs.length; i++) {var image = imgCol.filter(ee.Filter.eq("system:index", indexs[i])).first();image = image.clip(geometry);//tif数据下载Export.image.toDrive({image: image,description: 'pre_'+indexs[i],fileNamePrefix: 'pre_'+indexs[i],folder: 'pre',region: geometry,scale: 1000,crs: "EPSG:4326",maxPixels: 1e13});}});
}
exportImageCollection(dataset);

在这里插入图片描述
在这里插入图片描述

4.逐月数据分析和获取

4.1 代码一

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
var dataset = ee.ImageCollection('UCSB-CHG/CHIRPS/DAILY').filterBounds(geometry).select('precipitation');
//时间范围
var years = ee.List.sequence(2020, 2020);
var months = ee.List.sequence(1, 12);//将逐日数据生成月平均数据
var monthlysum =  ee.ImageCollection.fromImages(years.map(function (y) {return months.map(function(m) {return dataset.filter(ee.Filter.calendarRange(y,y, 'year')).filter(ee.Filter.calendarRange(m, m, 'month')).sum().set('year', y).set('month', m).set('system:time_start', ee.Date.fromYMD(y, m, 1));});}).flatten());print(monthlysum);
print(ui.Chart.image.series(monthlysum, geometry, ee.Reducer.mean(), 1000));function exportImageCollection(imgCol) {var indexList = imgCol.reduceColumns(ee.Reducer.toList(), ["system:index"]).get("list");indexList.evaluate(function(indexs) {for (var i=0; i<indexs.length; i++) {var image = imgCol.filter(ee.Filter.eq("system:index", indexs[i])).first();image = image.clip(geometry);//tif数据下载Export.image.toDrive({image: image,description: 'pre_2020_'+indexs[i],fileNamePrefix: 'pre_2020_'+indexs[i],folder: 'pre',region: geometry,scale: 1000,crs: "EPSG:4326",maxPixels: 1e13});}});
}
exportImageCollection(monthlysum);

在这里插入图片描述
在这里插入图片描述

4.2 代码二(简洁)

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset =ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY");
for(var i=2018;i<=2020;i++){for(var j=1;j<=12;j++){var data_collection = null;switch(j){case 1:case 3:case 5:case 7:case 8:case 10:case 12:data_collection = dataset.filterDate(i+'-'+j+'-01',i+'-'+j+'-31').select('precipitation');break;case 4:case 6:case 9:case 11:data_collection = dataset.filterDate(i+'-'+j+'-01',i+'-'+j+'-30').select('precipitation');break;case 2:data_collection = dataset.filterDate(i+'-'+j+'-01',i+'-'+j+'-28').select('precipitation');break;}var YR_collection = data_collection.sum().clip(geometry);Export.image.toDrive({image: YR_collection,description: i+'-'+j,fileNamePrefix: i+'-'+j,scale: 1000,region: geometry,maxPixels: 1e13,folder: 'pre'})}
}

在这里插入图片描述

5.逐年数据分析和获取

5.1 代码一

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
// 选择数据集并进行波段比例换算
var collection = ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY").filterDate('2000-01-01', '2020-12-31').select("precipitation");//进行年平均值的计算
var years = ee.List.sequence(2000, 2020);
var collectYear = ee.ImageCollection(years.map(function(y) {var start = ee.Date.fromYMD(y, 1, 1);var end = start.advance(12, 'month');return collection.filterDate(start, end).reduce(ee.Reducer.sum()).float().set('system:time_start',y).set('year',y);
}));
print(collectYear);
//年均值的时间序列展示
var Yearly_chart = ui.Chart.image.series({imageCollection: collectYear.select('precipitation_sum'),region: geometry,reducer: ee.Reducer.mean(),scale: 500,xProperty: 'year',}).setOptions({interpolateNulls: true,lineWidth: 2,title: 'pre Yearly Seires',vAxis: {title: 'pre'},hAxis: {title: 'Date'},});
print(Yearly_chart);function exportImageCollection(imgCol) {var indexList = imgCol.reduceColumns(ee.Reducer.toList(), ["system:index"]).get("list");indexList.evaluate(function(indexs) {for (var i=0; i<indexs.length; i++) {var image = imgCol.filter(ee.Filter.eq("system:index", indexs[i])).first();image = image.clip(geometry);//tif数据下载Export.image.toDrive({image: image,description: 'pre_'+indexs[i],fileNamePrefix: 'pre_'+indexs[i],folder: 'pre',region: geometry,scale: 1000,crs: "EPSG:4326",maxPixels: 1e13});}});
}
exportImageCollection(collectYear);

在这里插入图片描述
在这里插入图片描述

5.2 代码二(简洁)

var geometry = ee.FeatureCollection('users/www1573979951/luyixian');
Map.centerObject(geometry,6);
var dataset =ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY");for(var i=2018;i<=2020;i++){var data_collection = dataset.filterDate(i+'-01-01',i+'-12-31').select('precipitation');print(data_collection)var YR_collection = data_collection.sum().clip(geometry);Export.image.toDrive({image: YR_collection,description: i,fileNamePrefix: i,scale: 1000,region: geometry,maxPixels: 1e13,folder: 'pre'});
}

在这里插入图片描述


http://www.ppmy.cn/news/29909.html

相关文章

Elasticsearch的搜索命令

Elasticsearch的搜索命令 文章目录Elasticsearch的搜索命令数据准备URI Searchq&#xff08;查询字符串&#xff09;analyzer&#xff08;指定查询字符串时使用的分析器&#xff09;df&#xff08;指定查询字段&#xff09;_source&#xff08;指定返回文档的字段&#xff09;s…

第十四届蓝桥杯三月真题刷题训练——第 1 天

目录 题目1&#xff1a;数列求值 代码&#xff1a; 题目2&#xff1a;质数 代码&#xff1a; 题目3&#xff1a;饮料换购 代码&#xff1a; 题目1&#xff1a;数列求值 题目描述 本题为填空题&#xff0c;只需要算出结果后&#xff0c;在代码中使用输出语句将所填结果输出…

4. 字符设备驱动高级--- 下篇

文章目录一、字符设备驱动高级1.1 注册字符设备驱动新接口1.1.1 新接口与旧接口1.1.2 cdev介绍1.1.3 设备号1.1.4 编程实践1.1.5 alloc_chrdev_region自动分配设备号1.1.6 中途出错的倒影式错误处理方法二、字符设备驱动注册代码分析2.1 旧接口register_chrdev2.2 新接口regist…

Docker学习(十七)save 和 export 命令的区别

Docker 中有两个命令可以将镜像导出为本地文件系统中的 tar 文件&#xff1a;docker save 和 docker export。尽管它们的作用类似&#xff0c;但它们之间有一个重要的区别。 1.使用方式的不同&#xff1a; docker save 的使用示例&#xff1a; docker save -o test.tar image…

【C++】内存管理

C/C内存分布C语言中动态内存管理方式C内存管理方式new/delete操作内置类型new和delete操作自定义类型operator new与operator delete函数new和delete的实现原理定位new表达式(placement-new)C/C内存分布 在C/C中关于内存的知识是非常重要的&#xff0c;只有清楚的了解内存的相…

将微信小程序页面转为图片

最近做项目遇到一个需求,那就是要将某个页面转为图片然后传给后端,我仔细找了一圈,发现官方那个Api也就是wx.canvasToTempFilePath生成的图片很有可能为空,太坑了,于是我放弃用它了,选择了用wxml2canvas。 安装wxml2canvas npm init npm install wxml2canvas --save --…

MDK软件使用技巧

本文主要汇总MDK软件使用技巧 一、字体大小及颜色修改 第一步点击工具栏的这个小扳手图标 进去后显示如下&#xff0c;先设置 Encoding 为&#xff1a;Chinese GB2312(Simplified)&#xff0c;然后设置 Tab size 为&#xff1a;4 以更好的支持简体中文&#xff0c;否则&…

kafka

安装jdk;下载zookeepe并启动&#xff1b;下载地址&#xff1a;https://archive.apache.org/dist/zookeeper/zookeeper-3.5.8/apache-zookeeper-3.5.8-bin.tar.gz将conf目录下的zoo_sample.cfg复制一份重命名为zoo.cfg。启动&#xff1a;下载kafka并启动&#xff1b;下载地址&am…