IDRISI软件之CA_Markov模块实现土地利用变化模拟方法及步骤
一、首先创建一个工程目录
二、数据格式转换
所用的数据是IDRISI中的栅格数据，因此需要将gis中的tif数据转换为IDRISI支持的栅格数据格式。
方法：File→Import→Desktop Publishing Formats→GEOTIFF/TIFF
转换后的格式为.rst
三、获取马尔科夫矩阵
方法：Modeling→Environmental/Simulation models→MARKOV
1表示获取转换矩阵的前一期影像，为我们的87年遥感影像；
2表示获取转换矩阵的后一期影像，为我们的96年遥感影像；
3是这个模型中输出条件概率的前缀，表示的是从87到96变化的一些信息（具体是什么，我也不清楚，但是后续的预测会用到这个文件），一般都是我们自己命名，比如说8796；
4表示第一个与第二个影像之间的时间间隔，这里为9年；
5表示我们向前预测的时间周期，这里也设置为9年，即模拟2005年的土地利用情况；
6是比例误差（我看的资料里面一般都设置的是0.15）。
获取的马尔科夫矩阵记录了在下一个时期，从每个土地利用类型转换为其他土地利用类型的概率。
四、实现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
的冯诺依曼形状，如下所示：
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 Analysis→Database Query→CROSSTAB
1表示第一幅影像数据，这里的05为模拟产生的数据；
2表示第二幅影像数据，这里的05wuhan1为实际的数据。
其实，上述二者的顺序是可以调换的。
模型运行后得出的结果为
里面有不同的评价指标值。