GEE数据集:欧美1950-2022年扩展春季指数(SI-x)

目录

高分辨率扩展春季指数数据库

简介

数据集说明

空间信息

代码

代码链接

APP链接

结果

引用

许可

网址推荐

0代码在线构建地图应用

机器学习


高分辨率扩展春季指数数据库

简介

扩展春季指数(SI-x)为研究春季开始的时间及其与气候变化的关系提供了一个全面的数据集。 这些由每日最低和最高气温得出的模型可追踪主要植物物种的初叶和初花事件。 通过将温度数据转化为一致的指数,SI-x 可以计算冻害指数。 该数据集对北美(1980-2022 年)和欧洲(1950-2020 年)的春季物候进行了多年代、高分辨率(1 公里)分析。

数据集说明

空间信息

目前,1 公里 SI-x 产品可用于北美和中美两个研究区域(位于北纬 14°02'31.3 "N 至 55°37'04.1 "N 和西经 56°05'50.7 "W 至 126°22'06.1 "W 之间)。 该数据集采用了 1980 年至 2022 年的 Daymet 第 4 版。 每日最高和最低气温以及昼长可在 GEE 中获得。 这些 SI-x 产品是 Izquierdo-Verdiguier 等人的研究成果的更新版。

欧洲(位于北纬 35°55'48.7 "N 和 73°32'47.1 "N 之间,西经 10°36'29.5 "W 和 44°50'29.5 "E 之间)。 日最高气温和日最低气温来自 1950 年至 2020 年欧洲观测资料缩减版(E-OBS)第 3 版,可在此处下载。 昼长是在将数据导入 GEE 后计算(模拟)得出的(导入是因为 GEE 无法直接获得该数据)。

代码

var geometry = /* color: #d63000 */ee.Geometry.MultiPoint();
function linspace(a,b,n) {
    if(typeof n === "undefined") n = Math.max(Math.round(b-a)+1,1);
    if(n<2) { return n===1?[a]:[]; }
    var i,ret = Array(n);
    n--;
    for(i=n;i>=0;i--) { ret[i] = (i*b+(n-i)*a)/n; }
    return ret;
}

var lf_eu = ee.ImageCollection('projects/sat-io/open-datasets/SIx_products/LeafEuropev3')
.map(function(img){return img.mask(img.neq(0)).select([3],['leaf']).copyProperties(img)});
var bl_eu = ee.ImageCollection('projects/sat-io/open-datasets/SIx_products/BloomEuropev3')
.map(function(img){return img.mask(img.neq(0)).select([3],['bloom']).copyProperties(img)});
var di_am = ee.ImageCollection('projects/sat-io/open-datasets/SIx_products/DI_Daymetv4')
.map(function(img){return img.select([3],['damage']).copyProperties(img)});
var di_er = ee.ImageCollection('projects/sat-io/open-datasets/SIx_products/DI_Europev3')
.map(function(img){return img.select([3],['damage']).copyProperties(img)});

var yrs = linspace(1950,2022,73);

function compositeAreas(yr){
var image = ee.Image(ee.ImageCollection('projects/sat-io/open-datasets/SIx_products/LeafDaymetv4')
  .filter(ee.Filter.eq('system:index',yr.toString())).select('leaf').first());

var image2 = ee.Image(lf_eu.filter(ee.Filter.eq('system:index',yr.toString())).first());

var image3 = ee.Image(ee.ImageCollection('projects/sat-io/open-datasets/SIx_products/BloomDaymetv4')
  .filter(ee.Filter.eq('system:index',yr.toString())).select('bloom').first());

var image4 = ee.Image(bl_eu.filter(ee.Filter.eq('system:index',yr.toString())).first());
  
var image5 = ee.Image(di_am.filter(ee.Filter.eq('system:index',yr.toString())).first());

var image6 = ee.Image(di_er.filter(ee.Filter.eq('system:index',yr.toString())).first());

var final = ee.ImageCollection([image,image2]).mosaic()
.addBands(ee.ImageCollection([image3,image4]).mosaic())
.addBands(ee.ImageCollection([image5,image6]).mosaic());
return final.set('year',yr);
}

var pheno = ee.ImageCollection(yrs.map(compositeAreas));
print(pheno);

var vis_vi1 = {bands:['leaf'],min:0 , max: 250, palette:["af0000","eb1e00","ff6400","ffb300","ffeb00","9beb4a","33db80","00b4ff","0064ff","000096"]};
Map.addLayer(pheno,vis_vi1,'Leaf');

var vis_vi2 = {bands:['bloom'],min:0, max: 250, palette:["af0000","eb1e00","ff6400","ffb300","ffeb00","9beb4a","33db80","00b4ff","0064ff","000096"]};
Map.addLayer(pheno,vis_vi2,'Bloom');

var vis_vi3 = {bands:['damage'],min:-150 , max: 150, palette:["ff3636","ff7b71","ffafaf","ffe8de","ffffff","dae4ff","b1c5ff","6d93ff","3e70ff"]};
Map.addLayer(pheno,vis_vi3,'DI');

Map.setCenter(-32.639,40.407,3)

代码链接

https://code.earthengine.google.com/?scriptPath=users/sat-io/awesome-gee-catalog-examples:analysis-ready-data/EXTENDED-SPRING-INDICES

APP链接

 https://emma.users.earthengine.app/view/spring-onset

结果

 

引用

Izquierdo-Verdiguier, Emma, Raúl Zurita-Milla, Toby R. Ault, and Mark D. Schwartz. "Development and analysis of spring plant phenology products: 36
years of 1-km grids over the conterminous US." Agricultural and forest meteorology 262 (2018): 34-41.

许可

This work is licensed under a CC BY-NC 4.0 license.

Created by: Izquierdo-Verdiguier. 2024

Curated in GEE by : Samapriya Roy

Keyworks: spring onset, phenology, climate change

Last updated in GEE: 2024-08-29

网址推荐

0代码在线构建地图应用

https://www.mapmost.com/#/?source_inviter=CnVrwIQs 

机器学习

https://www.cbedai.net/xg 

  • 13
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是使用Python在GEE中计算2000-2022月降水数据和降水状况指数PCI的代码: ```python import ee ee.Initialize() # 定义区域 region = ee.Geometry.Polygon([[70, 25], [70, 40], [90, 40], [90, 25]]) # 加载数据集 chirps = ee.ImageCollection('UCSB-CHG/CHIRPS/DAILY').filterDate('2000-01-01', '2022-12-31').select('precipitation').filterBounds(region) # 计算月降水 monthly_precipitation = chirps.reduce(ee.Reducer.mean().combine(ee.Reducer.stdDev(), sharedInputs=True)).resample('monthly').set('region', region) # 计算PCI def compute_pci(image): pci = image.subtract(image.reduce(ee.Reducer.mean())).divide(image.reduce(ee.Reducer.stdDev())) return pci.set('system:time_start', image.get('system:time_start')).set('region', region) pci = chirps.map(compute_pci) # 导出数据 task = ee.batch.Export.table.toDrive(collection=monthly_precipitation, description='monthly_precipitation', fileFormat='CSV', selectors=['system:time_start', 'mean', 'stdDev', 'region']) task.start() task = ee.batch.Export.table.toDrive(collection=pci, description='pci', fileFormat='CSV', selectors=['system:time_start', 'precipitation', 'region']) task.start() ``` 这段代码会将GEE中UCSB-CHG/CHIRPS/DAILY数据集中2000-2022间的降水数据加载进来,然后计算每个月的平均降水和标准差,并将结果导出为CSV文件。同时,还会计算每个时间点的PCI,并将结果导出为CSV文件。 在这段代码中,我们使用了GEE提供的`ee.Reducer`类来计算平均值和标准差。我们还定义了一个`compute_pci`函数来计算PCI。PCI是每个时间点的降水值减去该时间点的平均降水值,然后除以该时间点的标准差。最后,我们使用`ee.batch.Export.table.toDrive`函数将结果导出为CSV文件。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

此星光明

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

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

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

打赏作者

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

抵扣说明:

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

余额充值