其他
ESA全球土地覆被数据1992-2018
ESA Global Land Cover products1992-2018
欧空局发布了2016,2017,2018三期300m的全球土地覆被数据,可以和1992-2015年数据衔接起来进行研究。不过新版数据使用的是nc格式,需要进行一些转换。
数据下载和介绍
数据下载
数据介绍可以从参考文献1网址查看,并进入数据下载界面。数据通过C3S Climate Data Store(CDS)分发。
进入CDS之后,点击下载数据,注册,选择数据,同意,即可下载。
数据介绍
可以从Documentation选项卡中下载数据介绍。关于地类的介绍可以直接从Product user guide中查看。
数据格式转换
重点讲一下数据格式的转换。2016年以后的数据为NC格式,一般我们更习惯于使用TIFF格式数据进行分析,这就需要进行一下格式转换。我主要尝试了两种方法,R语言和QGIS,都很方便,但是QGIS更快,推荐QGIS
R语言nc格式转换
非常简单,三行语句搞定,就是慢。
library(raster)
tmpbr <- brick("C3S-LC-L4-LCCS-Map-300m-P1Y-2016-v2.1.1.nc")
writeRaster(tmpbr, filename = "2016.tif", format="GTiff")
QGIS
QGIS可以直接打开nc文件,非常方便。我们主要是用lccs_calss图层,将这个图层直接拖进去,打开。
需要等一会,数据量比较大,略卡。打开后,在图层上右击,导出数据,需要注意的就是默认设置导出的TIFF文件是没有压缩的,很大,建议进行压缩。将Create Options打对勾,Profile选择High Compression,不过这个High Compression默认的压缩算法我没有用过,不清楚情况。
我一般使用LZW压缩,点击帮助,发现可以输入LZW,指定为LZW压缩方法。
参考文献
http://www.esa-landcover-cci.org/?q=node/197 NetCDF(nc)数据读取与格式转换
点击下方阅读原文即可跳转下载数据网站