(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

一、首先创建一个工程目录

二、数据格式转换

所用的数据是IDRISI中的栅格数据，因此需要将gis中的tif数据转换为IDRISI支持的栅格数据格式。

方法：File→Import→Desktop Publishing Formats→GEOTIFF/TIFF

(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

转换后的格式为.rst

三、获取马尔科夫矩阵

方法：Modeling→Environmental/Simulation models→MARKOV

(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

1表示获取转换矩阵的前一期影像，为我们的87年遥感影像；

2表示获取转换矩阵的后一期影像，为我们的96年遥感影像；

3是这个模型中输出条件概率的前缀，表示的是从87到96变化的一些信息（具体是什么，我也不清楚，但是后续的预测会用到这个文件），一般都是我们自己命名，比如说8796；

4表示第一个与第二个影像之间的时间间隔，这里为9年；

5表示我们向前预测的时间周期，这里也设置为9年，即模拟2005年的土地利用情况；

6是比例误差（我看的资料里面一般都设置的是0.15）。

获取的马尔科夫矩阵记录了在下一个时期，从每个土地利用类型转换为其他土地利用类型的概率。

四、实现CA_Markov模型预测

(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

土地利用变化模拟使用的是IDRISI软件中的CA-Markov模型, 位于Modeling→Environmental/Simulation models→CA_Markov。

1表示模拟05影像需要依据的影像，即为我们的96年遥感影像；

2表示马尔科夫转换矩阵面积文件，这里选择的是马尔科夫转换概率矩阵；

3即为转换适宜性图集（我是把从87转换为96年影像中产生的那个8796文件作为适宜性图集，一般都是自己重新做一个这种图集，需要道路、河流、坡度等信息，我之前也做过，但主观性特别强，而且出来的模拟精度很低，所以就舍弃了这个方法）；

4表示输出的土地利用变化数据，命名为05；

5表示元胞自动机循环次数，一般为两个年份之间间隔的整数倍，这里可以取9、18、27等等，但是数字越大，需要的模拟时间越多；

6是我们讨论的CA模型中邻域结构的设定，系统默认的是5*5型，即为我们所说的5*5的冯诺依曼形状，如下所示：

(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

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） 编辑输入过滤模板

(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

第一行表示7行

第二行表示7列

下面的矩阵构成邻域过滤矩阵，1表示有影响作用的邻域，反之用0表示。

3） 将该文件保存为.fil格式文件，命名为moore7.

五、精度评价

GIS Analysis→Database Query→CROSSTAB

(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

1表示第一幅影像数据，这里的05为模拟产生的数据；

2表示第二幅影像数据，这里的05wuhan1为实际的数据。

其实，上述二者的顺序是可以调换的。

模型运行后得出的结果为

(完整word版)IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤

里面有不同的评价指标值。