生态安全格局

记录一下自己的实验

---

1.   实验数据准备

1.1 土地利用/土地覆盖数据

1.1.1 1985 to 2022 30m土地覆盖数据

土地利用和土地覆盖有些区别

GLC_FCS30D: the first global 30-m land-cover dynamic monitoring product with fine classification system from 1985 to 2022

1985、1990、1995、2000~2022（每年）

Zhang, X., Liu, L., Chen, X., Gao, Y., Xie, S., Mi, J., 2021. GLC_FCS30: global land-cover product with fine classification system at 30 m using time-series Landsat imagery. Earth Syst. Sci. Data 13, 2753-2776 , https://doi.org/10.5194/essd-13-2753-2021.

数据类别包含35个土地覆盖子类别的精细分类系统，例如：10，雨养农田；11，草本覆盖农田；...；51~92，针叶林阔叶林等；...；150，稀疏植被；152，稀疏灌木，...；裸地，沼泽等

1.1.2 1985 to 2023 30m土地覆盖数据

The 30 m annual land cover datasets and its dynamics in China from 1985 to 2023

1985、1990~2023（每年）

数据包含：1农田，2森林，灌木，4草地，5水，6雪/冰，7裸地，8不透水湿地

Jie Yang, & Xin Huang. (2024). The 30 m annual land cover datasets and its dynamics in China from 1985 to 2023 [Data set]. 收入 Earth System Science Data (1.0.3, 卷 13, 期 1, 页 3907–3925). Zenodo. The 30 m annual land cover datasets and its dynamics in China from 1985 to 2023

1.1.3 其他的土地利用数据

例如：2015~2019的100m空间分辨率Copernicus Global Land Service的数据

（ee.ImageCollection("COPERNICUS/Landcover/100m/Proba-V-C3/Global"))

The European Space Agency (ESA) WorldCover 10 m

（ee.ImageCollection("ESA/WorldCover/v200"))

2.   实验

2.1 数据预处理

(1) 以30m土地利用数据为例。找到自己研究区的经纬度，文件先以精度划分。

 压缩包下载完成之后，再找到自己研究区的纬度。

 .tif文件 是一个具有23波段的影像，每一个波段对应的是每一年份的土地利用数据。例如 band 1对应2000年，band 23 对应2022年。需要在GIS中进行的步骤：影像镶嵌—转投影——重采样（至30m）——掩膜提取——重分类.....

将波段拉至Contents，右键导出为TIF。

左图为处理好的影像，下一步则是按需分类。

(2) 2019的100m空间分辨率Copernicus Global Land Service的数据,几行代码下载

类型划分如下，部分类型如图： 

Copernicus Global Land Service: Land Cover 100m: version 3 Globe 2015-2019: Product User Manual

下载好的影像（已经在GEE中处理过），直接进行重分类：

2.2 景观格局指数计算（fragstats）