全球陆表特征参量数据产品(GLASS产品)检索与预处理方法

“ 全球陆表特征参量数据产品(GLASS产品)用途广泛，就是检索起来有些困难，数据产品的名称没有行列号，和空间位置没有直接的对应关系，需要通过XML文件确定影像位置，进而检索所需的影像。本文对GLASS产品检索和预处理的方法进行一个介绍，欢迎大家前往Geodata网站下载数据，点击阅读原文即可跳转。”

封面图：俄罗斯贝加尔湖 贝加尔湖科考

GLASS产品名称如下，你是在这个名称里找不到关于位置的信息的！

01

GLASS产品介绍

—

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

02

数据检索

—

既然文件标题中没有位置信息，那么，位置信息从哪里来呢？

答案是从XML中找！XML中可以找到影像的四至和中心点坐标，我们就利用这个来查找定位影像。

例如，我要检索中国范围内的长波辐射数据，那么我可以先拿到所有的XML文件，将所有的XML文件放到一个目录下，写一段R语言代码，提取所有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中转换为点：直接添加XY数据生成的点可能在ArcMap中显示不全，建议导出数据到数据库中，形成一个新的点要素文件，如下所示：同样的数据，直接生成的点只显示了上图的一半区域，而下面地理数据库中的点能够显示完整，包含了中国区域：每个点都是一个影像的中心点，代表一幅GLASS数据产品，这个时候，我们在想挑选指定区域的数据就很容易了。直接框选，或者用下图所示的按位置选择工具，都可以。查看属性表里面的xmlnames字段，即可找到所需的GLASS产品名称：

03

GLASS数据查看与投影

—

GLASS数据不同的产品采用了不同的自定义的坐标系，ENVI中不能直接打开，需要手动指定一下坐标情况。如下图所示，将GLASS产品中Latitude和Longitude两个图层指定为Geographic Information，下面的Lwup、Lwdn、Lwnt作为可以显示的波段，即可打开产品。
接下来在ENVI中使用Reproject Raster工具做一下投影变换，同时转换为TIFF格式，就可以将导出的TIFF文件在ArcGIS中使用了。

点击阅读原文，即可前往GEODATA下载数据~