PLUS模型

准备工作

注意：所有土地利用数据行列号需一致，坐标系需一致，最好使用WGS-84坐标系，影响因子行列号对齐不做要求。文件名称和路径不能出现中文！！！
操作方法：
1、提取一期研究区土地利用数据，例如：2020年的。
2、提取2010年和2000年土地利用数据。按掩膜提取—环境设置—处理范围—捕捉栅格，这里选择2020年土地利用数据。

PLUS模型运行

1. 数据转换 ，Data Processing–Comvert LULCs to…；
   当出现Finish说明转化成功，数据会带有“_uc”的后缀。之后的在PLUS过程中都应该运用该数据

2. 提取土地利用扩张模块，PLUS——E L E
   输入土地利用数据，land expansion：输出路径和文件名，Convert：运行
   

3. 获取各用地发展概率及贡献率
   

4. 限制区制作。限制区一般是自然保护区或者开放性湖泊。
   1）研究区自然保护区和湖泊数据剪裁
   2）删掉保护区和湖泊多余字段，两个数据各添加一个字段，命名一致，赋值一致。这里本人添加的字段名叫“A”，都赋值为“0”
   3）研究区矢量数据也添加一个字段“A”，这里赋值为“1”，这里需要融合处理，得到研究区整体边界数据，内部的行政区划边界不要。
   

3）保护区数据和湖泊数据合并。地理处理——合并，输入合并数据，设置输出路径，其他不变。
4）合并后的数据包括很多个面，这里进行融合，变成1个面。地理处理——融合，融合字段选择“A”。
5）擦除研究区矢量数据中的自然保护区和湖泊数据。直接搜索擦除，输入要素选择研究区矢量数据，擦除数据选择第4步生成的数据。
7）转为栅格，字段选“A”，就得到限制数据，这里注意行列号要和LULC数据一致。
8）数据转换，参考“准备工作”。

5. CARS；
   Land use pattern ：当前LULC；Development potential：第3步生成的结果；Conversion constraint：限制要素，第4步生成结果。
   

Land Demands：需要输入未来LULC的栅格个数（F A 1），S A不需要输入

Transition Matrix：转移矩阵。1表示可以转移，0表示不能转移。根据具体情境设置。

Neighborhood Weights：邻域权重，可参考相关文献，也可以根据各用地类型扩张面积计算。

也可以这样计算权重

Color设置颜色。

6. 精度验证
   G T：2020年真实土地利用数据
   S r：第5步模拟的2020年土地利用数据
   Sampling rate：采样率
   I M：2010年真实土地利用数据
   Caculate 计算。
   
   Kappa Coefficient>0.75就可以了
   Overall Accuracy：整体精度，>0.75
   
   FoM精度
   
   注意：参数需要不断调试，没有固定的值。