第二篇GEE中使用降水数据集求年平均降水（1）

Gee Explorer

已于 2024-12-03 16:42:23 修改

阅读量792

点赞数 13

分类专栏： Google Earth Engine 文章标签： javascript

于 2023-12-30 11:58:40 首次发布

本文链接：https://blog.youkuaiyun.com/weixin_64898072/article/details/135303339

版权

Google Earth Engine 专栏收录该内容

16 篇文章

订阅专栏

一、主要信息

（一）数据提供者和数据集

JAXA Earth Observation Research Center、"JAXA/GPM_L3/GSMaP/v6/operational"

（二）可用时间

2014年03月01日00:00:00–至今

（三）空间分辨率

11132米

（四）时间分辨率

3小时

（五）所含主要波段

如下表

名字	单位	最小值	最大值	描述
hourlyPrecipRate	毫米/小时	0	204.88	逐小时降水率快照
hourlyPrecipRateGC	毫米/小时	0	200.36	根据雨量计调整的每小时降水率快照
observationTimeFlag	小时	-124.72	16.06	从文件开始时间开始的相对时间到微波辐射计（成像仪/测深仪）观测的时间

二、数据集的使用

（一）具体代码

使用此数据集计算一个区域（roi）的年平均降水（注意：1、这里需要自己定义一个感兴趣区域roi。2、使用此数据集不可超出它的可用时间范围。）

//1.JAXA/GPM_L3/GSMaP/v6/operational

//创建了一个变量dataset，用于存储从JAXA/GPM_L3/GSMaP/v6/operational数据集中获取的图像集合。

//通过.filter方法筛选了日期范围为'2014-03-01'至'2023-12-25'之间的图像。

var dataset = ee.ImageCollection('JAXA/GPM_L3/GSMaP/v6/operational')

                  .filter(ee.Filter.date('2014-03-01', '2023-12-25'));

                 

//创建一个变量precipitation，用于存储从dataset中选择的名为'hourlyPrecipRate'的降水率图层

var precipitation = dataset.select('hourlyPrecipRate');



//将名为roi的图层添加到地图中，roi表示感兴趣区域

Map.addLayer(roi);



//以roi为中心，将地图视图设置为缩放级别8

Map.centerObject(roi,8);



// 定义空间和时间范围

var startyear = 2015;

var endyear = 2023;

var startmonth = 1;

var endmonth = 1;



//使用ee.Date.fromYMD方法创建了一个起始日期，其中年、月和日分别由startyear、startmonth和1指定

var startdate = ee.Date.fromYMD(startyear, startmonth, 1);



//使用ee.Date.fromYMD方法创建了一个结束日期，其中年、月和日分别由endyear+1（加1表示下一年）、endmonth和1指定。

var enddate = ee.Date.fromYMD(endyear + 1, endmonth, 1);



// 创建了一个包含起始年份到结束年份范围的年份列表

var years = ee.List.sequence(startyear, endyear);



// 使用filterDate方法根据起始日期和结束日期筛选出在该时间范围内的降水数据，

//并使用sort方法按照'系统时间开始'属性进行逆序排序。

var P = precipitation.filterDate(startdate, enddate)

  .sort('system:time_start', false);

 

// 计算每年的降水量，创建了一个名为yearlyP的图像集合，其中每个图像代表每年的降水量。

//通过map方法遍历年份列表，对每一年的降水数据进行求和并设置相关属性。

var yearlyP = ee.ImageCollection.fromImages(

  years.map(function(y) {

    var w = P.filter(ee.Filter.calendarRange(y, y, 'year')).sum();

    return w.set('year', y)

      .set('date', ee.Date.fromYMD(y, 1, 1))

      .set('system:time_start', ee.Date.fromYMD(y, 1, 1));

  })

);



// 预定义了可视化参数，包括最小值、最大值和颜色调色板。

var p_viz = {

  min: 0.0,

  max: 30.0,

  palette:

      ['1621a2', 'ffffff', '03ffff', '13ff03', 'efff00', 'ffb103', 'ff2300'],

};



// 预定义了图表的标题、横轴标题和纵轴标题。

var title = {

  title: 'Yearly precipitation',

  hAxis: { title: 'Time' },

  vAxis: { title: 'Precipitation (mm)' },

};



//创建了一个图表，通过ui.Chart.image.seriesByRegion方法以区域为单位绘制图像集合的时间序列。

var chart = ui.Chart.image.seriesByRegion({

    imageCollection: yearlyP,

    regions: roi,

    reducer: ee.Reducer.mean(),

    scale: 2500,

    xProperty: 'system:time_start',

    seriesProperty: 'SITE'

  })

  .setOptions(title)

  .setChartType('ColumnChart');



// 打印图表，将其输出显示在控制台上

print(chart);