全球陆表特征参量数据产品GLASS数据批量下载、裁剪、格式转换

GLASS数据批量下载、裁剪、格式转换

全球陆表特征参量数据产品(GLASS产品)，由北京师范大学梁顺林教授团队自主研发，数据产品包括叶面积指数（LAI）、反照率（Albedo）、发射率（BBE）、光合有效辐射（PAR）、下行短波辐射（DSR）、净辐射（NR）、光合有效辐射吸收比（FAPAR）、植被覆盖率（FVC）、潜热（ET）和植被总初级生产力（GPP）等多种产品。GLASS产品是基于多源遥感数据和地面实测数据，反演得到的长时间序列、高精度的全球地表遥感产品。这些产品为研究全球环境变化提供了可靠的依据，能够广泛应用于全球、洲际和区域的大气、植被覆盖、水体等方面的动态监测，并与气温、降水等气候变化表征参数结合起来，应用于全球变化分析。

GLASS数据在国内由国家地球系统科学数据中心分发，提供免费下载、预处理服务

• 服务人员联系方式：

• 电话：[phone]；
• 邮箱：[e-mail]；
• QQ：3298863160

• 服务时间：

  如何找到靠谱的科研数据？免费的国家数据中心为你服务！

• 周一至周五
• 早8点-晚5点，中间12：00-1：30午休

下面给大家提供一个自己动手的方法：

数据下载

GLASS数据在美国马里兰大学提供免费的下载，随便下，无限制，自己动手

• 数据简介
• http://www.glass.umd.edu/Overview.html

Download界面提供数据下载

在这里我下载了2015年的发射率数据，以这个数据作为例子讲讲怎么下载和裁剪

• 下载地址
• http://www.glass.umd.edu/BBE/MODIS/0.05D/2015/

进入到2015年的发射率Broadband Emissivity(BBE)数据界面

• 可以看到，一年有很多数据，一个个点多麻烦，肯定要批量下载
• 批量下载方法：使用wget
• wget详细介绍见下面图片链接推文

wget批量下载遥感数据

• 打开wget后，先输入下载影像存储目录
• 然后用下面的代码启动下载，并查看下载进度

wget -m -np -nH --cut-dir=1 --wait=1  http://www.glass.umd.edu/BBE/MODIS/0.05D/2015/  -c -o MCD12Q1-data.log -b
tail -f 1624  MCD12Q1-data.log

数据转换

• 原始的GLASS数据是HDF格式，我们一般常用TIFF格式做分析
• 原始的数据是全球的，要用自己研究区的需要裁剪
• 在这我使用的是8天合成的数据，为了方便，我把图层名用时间命名
• 为了实现以上效果，我使用R语言terra包完成以上工作

library(terra)
gllist = list.files(path = "F:/RSimg/GLASS/MODIS/0.05D/2015", pattern = ".hdf$")   #读取HDF列表
gldir = paste0("F:/RSimg/GLASS/MODIS/0.05D/2015/", gllist)
gl1 = rast(gldir)   #读取所有的HDF
gldate = seq.Date(from = as.Date("2015-01-01"), to = as.Date("2015-12-31"), by=8)   #生成2015年8天时间序列
names(gl1) = gldate   #图层用日期命名
cmr = vect("./SHP/CMR.shp")    #读取研究区矢量
glcmr = trim(mask(gl1, cmr))        #裁剪
writeRaster(glcmr, filename="./GLASS/BBE2015.tif", names= glcmr@ptr[["names"]])   #输出TIFF

2015年前四期发射率数据

参考文献

1. http://www.geodata.cn/thematicView/GLASS.html

点击阅读原文前往马里兰大学下载GLASS数据