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

最新推荐文章于 2025-10-12 11:15:37 发布

原创最新推荐文章于 2025-10-12 11:15:37 发布 · 7.2k 阅读

44 ·

CC 4.0 BY-SA版权

QGIS基本功同时被 3 个专栏收录

31 篇文章

订阅专栏

QGIS教程

10 篇文章

订阅专栏

qgis

7 篇文章

订阅专栏

本文介绍了如何通过QGIS利用Esri 2020全球陆地覆盖数据进行区域分析，包括数据下载、投影分析、数据裁剪合并以及土地覆盖类型面积统计，为GIS从业者提供了实用的分析方法。

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为后缀区分北半球地区和南半球地区，如EPSG:32650 - WGS 84 / UTM zone 50N和EPSG:32750 - WGS 84 / UTM zone 50S，分别代表50号分带在北半球和南半球区域的CRS。

UTM投影包含的纬度范围从南纬80度到北纬84度，采用8度分带，编号从C到X。通过经度分带带号与纬度分带带号的组合（如50T）即可确定一块地理区域范围。如下图所示：

2020年全球陆地覆盖数据GeoTIFF文件名的前三位（两位数字加一位字母）代表的正是UTM投影经度带号和纬度带号的组合编号。因此，根据GeoTIFF文件名，结合上图，即可判断出文件对应的大致地理区域。

我国位于UTM投影43-53经度分带内，将以43-53之间的数字为开头的GeoTIFF文件加载到QGIS中，将地图窗口CRS设为WGS84，全幅显示和按照栅格图层原始分辨率（100%）显示的效果分别如下：

04 数据分析案例

下面以河南省为研究区域，制作全省范围的土地覆盖图，并统计各种覆盖类型的土地面积。

示范数据：

1:100万全国行政境界图层（BOUA.shp），来自天地图网站（www.tianditu.gov.cn）数据API服务。
与河南省境界范围有交集的2020年10米分辨率土地覆盖数据，包括三个GeoTIFF文件：49S_20200101-20210101.tif、50R_20200101-20210101.tif、50S_20200101-20210101.tif，来自Esri网站（www.arcgis.com）。

示范数据打包下载地址：

https://pan.baidu.com/s/1B6ruGVG1MiPiCUQR4I85lA

提取码：mezq

特别说明：进行后续详细操作步骤之前，请先将解压得到的名为Henan的文件夹放置在不包含中文、特殊字符及空格的文件路径中，否则后续操作步骤可能无法正确执行。

具体操作步骤如下：

1.启动QGIS，加载BOUA.shp图层，将图层“填充样式”设为“无填充”，便于与栅格图层叠加显示。

2.在“图层”面板中的BOUA图层名上点击鼠标右键，在弹出菜单中选择“过滤”，在弹出对话框下方的“提供过滤表达式”中输入："PAC" between 410102 and 419001，然后点击“OK”按钮，BOUA图层将仅显示河南省范围。

备注：PAC取值在410102和419001之间的区域为河南省县级行政单元。

3.将三个栅格图层加载到地图窗口中，然后将BOUA图层移动到“图层”面板最上方。

4.点击菜单“栅格-->提取-->按掩膜图层裁剪栅格”。

5.弹出“按掩膜图层裁剪栅格”对话框，在“输入图层”下拉列表中选择50R_20200101-20210101，在“掩膜图层”下拉列表中选择BOUA，确保勾选“将剪切栅格的范围扩展到掩膜图层范围”，在“高级参数”面板的“配置”下拉列表中选择“高压缩”，该选项有助于显著压缩裁剪生成的栅格文件大小。设置完毕后，点击“运行”按钮。

6.视计算机配置不同，裁剪操作通常需要十几分钟乃至更长时间。裁剪生成的新栅格图层默认名称为“已裁剪（掩膜）”，在该图层名上点击鼠标右键，在弹出菜单中选择“重命名图层”，将图层名称修改为CLIP_50R。

7.针对图层50S_20200101-20210101重复上述步骤4-6生成新栅格图层CLIP_50S，针对图层49S_20200101-20210101 重复上述步骤4-6生成新栅格图层CLIP_49S，在“图层”面板中移除三个原始栅格图层，仅保留BOUA图层和三个新栅格图层。可以看出，这时地图窗口仅显示河南省范围内的土地覆盖数据。

8.下面需要将新生成的三个栅格图层合并为一个栅格图层，以便于进行后续的数据统计。QGIS的栅格合并操作要求输入的栅格图层必须具有相同的CRS，CLIP_50R和CLIP_50S图层的CRS都是EPSG:32650 - WGS 84 / UTM zone 50N，而CLIP_49S图层的CRS则是EPSG:32649 - WGS 84 / UTM zone 49N，需要先进行统一处理。在CLIP_49S图层名上点击右键菜单“导出-->另存为”，在弹出的“栅格地图另存为”对话框中，将输出的栅格文件命名为“CLIP_49S_REPRJ”，将输出图层的CRS设为EPSG:32650 - WGS 84 / UTM zone 50N（通过点击CRS下拉列表右侧的“选择CRS”按钮），在“创建选项”面板的“配置”下拉列表中选择“高压缩”，点击“OK”按钮。

9.导出的CLIP_49S_REPRJ图层将会自动添加到地图窗口中，默认显示为灰色。在CLIP_49S图层上点击鼠标右键“样式-->复制样式”。

10.在CLIP_49S_REPRJ图层上点击鼠标右键“样式-->粘贴样式”。

11.CLIP_49S_REPRJ图层将和CLIP_50S、CLIP_50R图层一样显示为彩色，且三个图层具有完全相同的CRS：EPSG:32650 - WGS 84 / UTM zone 50N。移除CLIP_49S图层，保留BOUA、CLIP_50S、CLIP_50R、CLIP_49S_REPRJ四个图层。

12.点击菜单“栅格-->杂项-->合并”。

13.在弹出的“合并”对话框中，点击“输入图层”右侧的省略号按钮进而选中全部三个图层，在“输出类型”下拉列表中选择“Byte”（与原始栅格文件保持一致），在“高级参数”面板中的“配置”下拉列表中选择“高压缩”，点击“合并输出”文本框右侧的省略号按钮，指定合并生成的栅格文件名称（本文将其命名为Merge.tif）和存放路径，最后点击“运行”按钮。