九种NPP\GPP数据集介绍和下载

数据简介

中国要在2030年实现碳达峰,2060年实现碳中和。如何度量生态系统碳源汇,并检测其时空的动态分布?科学家选择了陆地碳循环中的关键部分 NPP,将其视作生态系统能量流动和生态系统功能的指示器。

总初级生产力GPP(Gross Primary Production)是在一定时间内生态系统生产者获得化学能,并将其储存于生物量中的速率。

被植物所固定的能量中的部分被用于自身生长R,剩余部分的生产量就是净初级生产力NPP(Net Primary Production
) 。

目前中文互联网关于NPP的数据集的总结博客几乎没有,我这里简单总结了九种常见的NPP及GPP数据集。

全球数据集

主要是介绍几种覆盖全球的NPP/GPP数据集。

MODIS GPP\NPP (500m)

MODIS GPP\NPP数据是目前使用最广泛的数据集,该数据集已完成产品与地面参考点的比较验证,并提供14个数据类别可供用户下载。

名称 类型 分辨率 合成方式 来源
MOD17A2H v006 GPP 500m 8天 Terra
MOD17A2H v061 GPP 500m 8天 Terra
MOD17A2HGF v006 GPP 500m 8天 Terra
MOD17A2HGF v061 GPP 500m 8天 Terra
MOD17A3H v006 NPP 500m 1年 Terra
MOD17A3HGF v006 NPP 500m 1年 Terra
MOD17A3HGF v061 NPP 500m 1年 Terra
MYD17A2H v061 GPP 500m 8天 Aqua
MYD17A2H v006 GPP 500m 8天 Aqua
MYD17A2HGF v006 GPP 500m 8天 Aqua
MYD17A2HGF v061 GPP 500m 8天 Aqua
MYD17A3H v006 NPP 500m 1年 Aqua
MYD17A3HGF v061 NPP 500m 1年 Aqua
MYD17A3HGF v006 NPP 500m 1年 Aqua

其中,Terra卫星的观测数据即MOD数据,aqua卫星观测数据为MYD数据,Terra是于当地的上午10:30过境,auqa是于当地下午1:30过境。A2是GPP数据,每8天进行数据合成。A3是NPP数据,每一年进行合成。

如果需要下载该系列数据,首先进入USGS的资源搜索栏(https://lpdaac.usgs.gov/product_search/),搜索Primary Produ

在R语言中下载NASA MODIS GPP ( Gross Primary Production ) NPP ( Net Primary Production) 数据通常需要借助一些专门处理遥感数据的包,如`raster`、`MODIStsp` 或 `modisdata`。以下是大致步骤: 1. **安装必要的包**: 首先,确保已经安装了`raster`, 如果没有,可以使用`install.packages("raster")`安装。 2. **访问MODIS TERRAAQUA网站**: NASA的MODIS数据可以从其官方网站(https://ladsweb.modaps.eosdis.nasa.gov/)获取。但是,直接下载大文件可能会很慢,因此推荐使用专门提供MODIS数据的API,如`MODIStsp`。 3. **安装MODIStsp**: 使用`remotes::install_github("ropensci/MODIStsp")`安装`MODIStsp`包。然后,通过`library(MODIStsp)`加载它。 4. **定义时间范围产品ID**: 对于GPPNPP,产品ID分别为MCD17A3H(GPPMCD17A4H(NPP)。确定你感兴趣的起止日期,例如: ```R start_date <- "YYYY-MM-DD" end_date <- "YYYY-MM-DD" ``` 5. **查询并下载数据**: 使用`get.MODIS()`函数从`MODIStsp`下载数据,指定产品ID、起始日期、结束日期以及输出目录(如果还未创建的话,需要先创建): ```R gpp_data <- get.MODIS(product = "MCD17A3H", startDate = start_date, endDate = end_date, path = "path/to/output/directory") npp_data <- get.MODIS(product = "MCD17A4H", ..., path = ...) ``` 其中省略的部分与`gpp_data`部分类似。 6. **处理下载的数据**: 下载的文件可能是Tiff格式,使用`raster`包打开操作它们: ```R gpp_raster <- raster(gpp_data$filename) npp_raster <- raster(npp_data$filename) ``` 注意:这个过程可能会受限于数据权限、网络状况以及MODIS服务器的响应速度。此外,实际操作时可能需要根据具体的时间地理位置调整参数。
评论 7
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

叫我锐多宝

请我喝杯啤酒吧

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

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

打赏作者

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

抵扣说明:

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

余额充值