目录

一、准备工作:

1.数据准备:

2.软件:

二、具体操作

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

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

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

4.调整坐标系为 WGS-1984

5.导出图层样式:

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

7.上传geoserver

8.图例


一、准备工作:

1.数据准备: 

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

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

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

        ②研究区行政规划shp文件

2.其他:

        软件:QGIS

        网站:Geoserver

二、具体操作

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

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

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

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

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

土地利用和弦图 R语言 土地利用图怎么获取_栅格

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

        4景图像坐标系不同,需调整至统一坐标系:

土地利用和弦图 R语言 土地利用图怎么获取_土地利用和弦图 R语言_02

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

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

栅格→杂项→合并→选择保存文件名、CRS(坐标系)、创建选项:配置:高压缩,OK

4.调整坐标系为 WGS-1984

导出→另存为→选择保存文件名、CRS(坐标系)、创建选项:配置:高压缩,OK

土地利用和弦图 R语言 土地利用图怎么获取_arcgis_03

5.导出图层样式:

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

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

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

土地利用和弦图 R语言 土地利用图怎么获取_图层_04

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

土地利用和弦图 R语言 土地利用图怎么获取_数据_05

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

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

7.上传geoserver

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

8.图例

土地利用数据图例

土地利用

类型

颜色

颜色值


#419bdf


#397d49

被淹没植被


#7a87c6

农田


#e49635

建筑


#c4281b

裸露地面


#a59b8f

冰、雪


#a8ebff


#616161

草地


#e3e2c3