实例 | 使用QGIS分析Esri2020土地覆盖数据

​ 01 前言

2021年6月24日，Esri发布了空间分辨率为10米的2020年全球陆地覆盖数据（Esri 2020 Land Cover），由728个GeoTIFF文件组成，总数据量接近69GB。关于数据的基本情况，可以参考：

www.esri.com/about/newsroom/announcements/esri-releases-new-2020-global-land-cover-map

由于具有覆盖全球、时效性好（2020年）、分辨率高（10米）、开源（基于知识共享协议）的特点，该数据的发布引发较大反响，成为近期GIS行业的热门话题。

热点易逝，而数据永恒。作为专注于QGIS技术的公众号，本文在简要介绍此数据集的基础上，重点举例说明了如何通过QGIS分析和应用此数据集，完成特定区域土地覆盖地图的制作和各种土地覆盖类型的面积统计。

 02 数据下载方式

可以通过以下地址直接下载全部文件，文件为zip压缩格式，压缩后的大小为60GB左右。由于数据量巨大，假定下载速度为5MB/秒，需要将近3个半小时时间。下载链接如下：

https://ai4edataeuwest.blob.core.windows.net/io-lulc/io-lulc-model-001-v01-composite-v03-supercell-v02-clip-v01.zip

也可以通过arcgis.com网站逐个下载需要的文件，Esri通过以下页面提供了按照单个GeoTiff文件进行下载的方式，点击地图上的某个方格，即可下载该方格所在地理区域对应的数据文件：

https://www.arcgis.com/apps/instant/media/index.html?appid=fc92d38533d440078f17678ebc20e8e2

 03 数据投影与文件命名分析

将下载的数据文件解压缩后，共有728个GeoTIFF文件，文件名从01C_20200101-20210101.tif起，至60W_20200101-20210101.tif止。所有文件名都以_20200101-20210101结尾，代表数据是2020年时间段的。而文件名开头的两位数字加一位字母编号代表什么呢？这要从数据采用的投影（Projection）——通用横轴墨卡托投影（Universal Transverse Mercator，简称UTM投影）说起。

UTM投影基于WGS84坐标系，将地球分为60个南北走向的分带，每个分带所跨经度范围为6度。分带起点位于西经180度，终点位于东经180度，自西向东编号，从1编至60。在QGIS中，UTM投影每个经度分带都有两个独立的CRS，并以N和S为后缀区分北半球地区和南半球地区，如EP