Fragstats从数据预处理到运行结果详解过程

Fragstats从数据预处理到运行结果详解过程

（1）Tif土地利用数据预处理：Arcgis10.2

a.使用重分类后将土地利用数据设置为如图形式，每个数值代表一种类别，我的数据共有5类，数值分别为1、2、5、7、11，分别代表水体、林地、耕地、建筑用地、草地，文件名为"Reclass_tif2.tif"。

b.将影像的地理坐标系GCS_WGC_1984，转为投影坐标系，我用的投影坐标系为Krasovsky_1940_Albers坐标系。（这里使用定义投影，不会投影数据，将得到无效的结果）

首先，新建一个空白.mxd文件，随便导入一个空间参考系为Krasovsky_1940_Albers的数据A，此时数据框的坐标系将随着第一个导入的数据（即数据A）变为Krasovsky_1940_Albers投影坐标系；

然后，加载重分类后的影像"Reclass_tif2.tif"，右键导出数据，在空间参考一栏选择为（数据框参考系），是否将NoData添加到栅格数据集中，选择是。

最后，按掩膜提取出感兴趣的区域，得到“roi_LC.tif”, 之后再次进行数据导出，设置NoData为-9999（不设置容易报错）

（2）景观格局指数计算：Fragstats4.2

a. 新建工程，打开影像，选择GeoTIFF grid格式（我的16G内存输入文件超过600M会报错，使用了32G内存1G以下能跑）

b. 加载四个参数文件

文件1：类别描述文件 "class.fcd’，更改影像值及对应类别名称即可，以下为文件内容：

ID,Name,Enabled,IsBackground
1,water,true,false
2,tree,true,false
5,crop,true,false
7,build,true,false
11,grass,true,false

文件2：边缘深度文件 "edgeDepth.fsq’，更改影像值及对应类别名称，更改各个类别之间边缘深度值（0~100），矩阵为5*5对称矩阵，以下为文件内容：

FSQ_TABLE
CLASS_LIST_LITERAL(water,tree,crop,build,grass)
CLASS_LIST_NUMERIC(1,2,5,7,11)
0,30,40,10,60
30,0,50,10,60
40,50,0,20,60
10,10,20,0,10
60,60,60,10,0

文件3：边缘差异度文件 "contrast.fsq’，更改影像值及对应类别名称，更改各个类别之间差异度值（0~1），矩阵为5*5对称矩阵，以下为文件内容：

FSQ_TABLE
CLASS_LIST_LITERAL(water,tree,crop,build,grass)
CLASS_LIST_NUMERIC(1,2,5,7,11)
0,0.6,0.5,0.4,0.6
0.6,0,0.2,0.6,0.3
0.5,0.2,0,0.6,0.5
0.4,0.6,0.6,0,0.6
0.6,0.3,0.5,0.6,0

文件4：相似度文件 "similarity.fsq’，更改影像值及对应类别名称，更改各个类别之间相似度（0~1），矩阵为5*5对称矩阵，以下为文件内容：

FSQ_TABLE
CLASS_LIST_LITERAL(water,tree,crop,build,grass)
CLASS_LIST_NUMERIC(1,2,5,7,11)
1,0.2,0.3,0.2,0.3
0.2,1,0.5,0.3,0.6
0.3,0.5,1,0.3,0.6
0.2,0.3,0.3,1,0.4
0.3,0.6,0.6,0.4,1

c. 设置分析参数和要计算的景观指数，分析参数中设置斑块、类别、景观尺度，根据需要计算分景观指数所属尺度选择分析参数中的对应尺度，如果多选或少选均会报错，这里我计算了斑块尺度的AREA和类别尺度的PD、AI、LSI。

d. 点击运行，生成结果。

GIS搬砖日常，与君共勉，点赞的UU今年必暴富，嘿嘿！