一、景观生态风险指数计算
0 准备工作
数据：
1.不同时期的土地利用数据
2.研究区矢量边界
插件：
格网裁剪影像插件tool(下载链接：https://pan.baidu.com/s/1r5Fjp8KeH5w1f-0xe3d5xQ 提取码：7qy9)
软件：arcgis、total commander、Fragstats4.2
1 生成研究区渔网
1.打开创建渔网工具，设置输出路径、模板范围、像元高度（本文设置20km）、行列数（设置为0或者空值可进行自动计算）、几何类型（面）
注：像元高度请结合研究区大小参考相关文献，本文设置20km仅为方便计算。
2.按位置进行选择，对研究区内的渔网进行选择，导出。
2 分割渔网
1.打开属性表，添加字段LID（文本型，长度20）如图2-1；右键新建的表头LID，打开字段计算器，进行如下字段计算如图2-2；得到结果如图2-3。
2-1 2-2 2-3
2.使用分割工具，进行一下设置，分割结果位于新建的数据库中。
注意：目标工作空间是数据库，需要在分割前新建一个文件地理数据库。
3 影像分割
1.将上文中的插件添加到工具箱中，打开工具箱进行设置参数。计算得到若干tif相关文件（时间取决于分割小格的数量）。
注意：文件名一定设置为%Name%.tif
4 景观格局指数计算
1.将文件中的0.tif(当然也可以是其他的.tif)导入到Fragstats4.2中，对其进行设置，点击ok。
注意：左侧窗口选择GDAL Geo TIFF grid(.tif)，右侧背景值为999，文件的路径必须是英文路径，否则会出错。（本文已对路径进行了修改）
2.然后不进行任何操作进行，将其导出；点击左侧export batch，选择路径，生成2000.fbt文件（默认文件格式）；将其后缀修改为2000.csv，双击打开；
3.使用total commander查找存储的影像小格文件的路径，点击扩展名，扩展名为tif的全选，然后在标记中选中复制完整路径及名称。
4.复制在2000.csv中，并将其他列都修改
5.修改2000.csv的后缀名为2000.fbt，然后打开Fragstats4.2，点击import batch，将2000.fbt导入，导入完成在左侧会出现所有的影像斑块，选择合适的斑块级别、景观级别指标进行运算（Fragstats4.2软件使用请在网上查阅资料），将计算结果另存为2000，系统会默认生成2000.class和2000.land文件
6.将生成的2000.class和2000.land文件后缀名分别修改为2000_class.csv和2000_land.csv，打开进行景观指数运算，得到景观生态风险指数并另存为2000ERI.csv。
（参考： [1]侯蕊,李红波,高艳丽.基于景观格局的武汉市江夏区土地利用生态风险评价
[2]侯蕊,李红波,高艳丽.基于景观格局的武汉市江夏区土地利用生态风险评价
[3陈心怡 谢跟踪 张金萍.海口市海岸带近30年土地利用变化的景观生态风险评价
等文献）
二、景观生态风险制图
0添加X、Y字段
1.打开属性表，新建X、Y字段
1 计算质心
1.计算X、Y的质心坐标。分别右键X、Y字段，进行计算几何。
2 将景观生态风险指数与渔网连接
1.打开arcgis，将渔网和生态风险指数加载，连接2000ERI.csv(LID字段相同，以LID字段进行连接)
3 导出并以XY数据显示，另存
将连接好的属性表导出显示，右键导出文件，选择显示XY数据。将生成的点文件另存为2000point.shp。
4 克里金插值
1.打开克里金插值的工具，将数据添加，设置参数，本文使用默认（根据需求进行修改）。
注意：需要将研究区20km加载，在处理范围里的范围和栅格分析里的掩膜进行选择。
5出图
1.加载边界数据，进行掩膜，设置色带显示。
三、写在最后
1.本文仅此为提供方法，很多细节没有注意，望大家在学习借鉴时注意起来！！
2.这也是我第一次写博客，有错误希望大家提出指正！！