10米分辨率土地利用图LUCC的下载、出图

目录

一、准备工作：

1.数据准备：

2.软件：

二、具体操作

1.加载土地利用数据，加载研究区行政规划shp文件。

2.裁剪土地利用数据至研究区。

3.合并上一步裁剪后得到的4景研究区内栅格图层。

4.调整坐标系为 WGS-1984

5.导出图层样式：

6.计算不同土地利用面积

7.上传geoserver

8.图例

---

• 一、准备工作：

1.数据准备： 

        ①10米分辨率研究区土地利用数据

        下载地址：Esri | Sentinel-2 Land Cover Explorer网站地址

Esri | Sentinel-2 Land Cover ExplorerThis application displays the land use/land cover (LULC) timeseries layer derived from ESA Sentinel-2 imagery at 10m resolution.https://livingatlas.arcgis.com/landcoverexplorer/#mapCenter=102.595%2C39.45%2C5&mode=step&timeExtent=2022%2C2022&year=2022&renderingRule=0&downloadMode=true

        点击、下载研究区所需土地利用数据

        ②研究区行政规划shp文件

2.其他：

        软件：QGIS

        网站：Geoserver

• 二、具体操作

        加载影像、shp文件→4景土地利用图像分别剪裁至shp→转为同一坐标系（EPSG:32649 - WGS 84 / UTM zone 49N）→调整为同一样式渲染（导出渲染样式，后期发布Geoserver样式）→合并→另存为WGS - 1984 坐标系→计算不同土地利用方式所占面积。

1.加载土地利用数据，加载研究区行政规划shp文件。

2.裁剪土地利用数据至研究区。

        将4景10米土地覆盖数据分别和研究区行政规划shp图层裁剪

        栅格→提取→按掩膜图层裁剪栅格→选择保存文件名、创建选项：配置：高压缩，保存

3.合并上一步裁剪后得到的4景研究区内栅格图层。

        4景图像坐标系不同，需调整至统一坐标系：

        统一样式：调整坐标系后栅格图像显示为默认的灰色，右键原图像→样式→复制样式，右键灰色图像→样式→粘贴样式。

        将4景统一坐标系的栅格图层合并为一张图像：

栅格→杂项→合并→选择保存文件名、CRS（坐标系）、创建选项：配置：高压缩，OK

4.调整坐标系为 WGS-1984

导出→另存为→选择保存文件名、CRS（坐标系）、创建选项：配置：高压缩，OK

5.导出图层样式：

        删除多余无用的分类值（原始土地利用数据样式有255个分类，删除13-255分类值），样式→保存样式→选择保存文件、名称，保存为.qml文件。

6.计算不同土地利用面积

        工具箱→栅格分析→栅格图层唯一值报告→唯一值表格：选择保存路径、名称，格式选为：CSV文件，运行。

        得到的CSV文件为EXCEL表格格式，表格中value、count、（~）2字段分别代表像元取值、像元数量（个）、面积统计（单位：平方米）。（0值为背景值）

        将单位平方米换算为平方公里。

        换算公式：1平方公里=1,000,000平方米

7.上传geoserver

        上传图像数据，上传图层样式，发布图层。

8.图例

土地利用数据图例
土地利用 类型	颜色	颜色值
水		#419bdf
树		#397d49
被淹没植被		#7a87c6
农田		#e49635
建筑		#c4281b
裸露地面		#a59b8f
冰、雪		#a8ebff
云		#616161
草地		#e3e2c3