2022年IDRISI软件之CA-Markov模块实现土地利用变化模拟方法及步骤

1、IDRISI 软件之 CA_Markov 模块实现土地利用变化模拟方法及步骤 一、第一创建一个工程目录二、数据格式转换所用的数据是IDRISI 中的栅格数据, 因此需要将gis 中的 tif数据转换为IDRISI 支持的栅格数据格式；方法： File Import Desktop Publishing FormatsGEOTIFF/TIFF 转换后的格式为 .rst 三、猎取马尔科夫矩阵方法： Modeling Environmental/Simulation modelsMARKOV 1 表示猎取转换矩阵的前一期影像,为我们的 87 年遥感影像；2 表示猎取转换矩阵的后一期影像,为我们的 9

2、6 年遥感影像；3 是这个模型中输出条件概率的前缀,表示的是从87 到 96 变化的一些信息（具体是什么,我也不清楚, 但是后续的推测会用到这个文件）,一般都是我们自己命名,比如说 8796；4 表示第一个与其次个影像之间的时间间隔,这里为 9 年；5 表示我们向前推测的时间周期,这里也设置为 9 年,即模拟 2022 年的土地利用情形；6 是比例误差（我看的资料里面一般都设置的是 0.15 ）；猎取的马尔科夫矩阵记录了在下一个时期,型的概率；从每个土地利用类型转换为其他土地利用类四、实现 CA_Markov模型推测CA-Markov模型 , 位于Modeling 土地利用变化模拟使用的是ID

3、RISI软件中的Environmental/Simulation modelsCA_Markov；96 年遥感影像；1 表示模拟 05 影像需要依据的影像,即为我们的2 表示马尔科夫转换矩阵面积文件,这里选择的是马尔科夫转换概率矩阵；3 即为转换适宜性图集（我是把从87 转换为 96 年影像中产生的那个8796 文件作为适宜性图集,一般都是自己重新做一个这种图集,需要道路、河流、坡度等信息,我之前也做过,但主观性特别强,而且出来的模拟精度很低,所以就舍弃了这个方法）；9、18、4 表示输出的土地利用变化数据,命名为05；5 表示元胞自动机循环次数,一般为两个年份之间间隔的整数倍,这里可以取27

4、 等等,但是数字越大,需要的模拟时间越多；6 是我们争辩的CA模型中邻域结构的设定,系统默认的是5*5 型,即为我们所说的5*5的冯诺依曼形状,如下所示：0 0 1 0 0 0 1 1 1 0 1 1 1 1 1 0 1 1 1 0 0 0 1 0 0 假如需要自己设置的话,可以选择上述对话框中的User-defined filter,但是里面的文件需要我们自己制作,方法如下（以7*7 摩尔结构为例） ：1） 点击工具栏上的Edit按钮：2） 编辑输入过滤模板第一行表示 7 行其次行表示 7 列下面的矩阵构成邻域过滤矩阵,1 表示有影响作用的邻域,反之用 0 表示；3） 将该文件储存为 .fil 格式文件,命名为 moore7. 五、精度评判GIS AnalysisDatabase Query CROSSTAB 1 表示第一幅影像数据,这里的 05 为模拟产生的数据；2 表