其他
全球陆表特征参量数据产品(GLASS产品)检索与预处理方法
封面图：俄罗斯贝加尔湖 贝加尔湖科考
GLASS产品名称如下，你是在这个名称里找不到关于位置的信息的！
01
GLASS产品介绍
—
全球陆表特征参量数据产品(GLASS产品)，由北京师范大学梁顺林教授团队自主研发，数据产品包括叶面积指数（LAI）、反照率（Albedo）、发射率（BBE）、光合有效辐射（PAR）、下行短波辐射（DSR）、净辐射（NR）、光合有效辐射吸收比（FAPAR）、植被覆盖率（FVC）、潜热（ET）和植被总初级生产力（GPP）等多种产品。GLASS产品是基于多源遥感数据和地面实测数据，反演得到的长时间序列、高精度的全球地表遥感产品。这些产品为研究全球环境变化提供了可靠的依据，能够广泛应用于全球、洲际和区域的大气、植被覆盖、水体等方面的动态监测，并与气温、降水等气候变化表征参数结合起来，应用于全球变化分析。
02
数据检索
—
既然文件标题中没有位置信息，那么，位置信息从哪里来呢？
答案是从XML中找！XML中可以找到影像的四至和中心点坐标，我们就利用这个来查找定位影像。
library(XML)#XML文件所在目录：xmlnames <- list.files(path = "D:/R/GLASSto TIFF/xml", pattern = '.xml$')xmldir <- paste("D:/R/GLASSto TIFF/xml/",xmlnames,sep="") xmlfiles <- lapply(xmldir, xmlParse)MyContent <- lapply(xmlfiles, xmlToList)#提取中心点经度subLon <- function(x){ x$GranuleURMetaData$SpatialDomainContainer$HorizontalSpatialDomainContainer$GPolygon$Center$PointLongitude}CenterLon <- sapply(MyContent, subLon)#提取中心点纬度subLat <- function(x){ x$GranuleURMetaData$SpatialDomainContainer$HorizontalSpatialDomainContainer$GPolygon$Center$PointLatitude}CenterLat<- sapply(MyContent, subLat)#XML文件名、经纬度合并为一个数据框XMLCenter <- cbind(xmlnames, CenterLon, CenterLat)#输出XML文件名和对应数据中心点坐标write.csv(XMLCenter, "XMLCenter.csv")将经纬度坐标在ArcMap中转换为点：03
GLASS数据查看与投影
—
点击阅读原文，即可前往GEODATA下载数据~