IDRISI软件之CAMarkov模块实现土地利用变化模拟方法和步骤.docVIP

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摩尔结构为例）： 点击工具栏上的Edit按钮： 编辑输入过滤模板 第一行表示7行 第二行表示7列 下面的矩阵构成邻域过滤矩阵，1表示有影响作用的邻域，反之用0表示。 将该文件保存为.fil格式文件，命名为moore7. 五、精度评价 GIS Analysis→Database Query→CROSSTAB 1表示第一幅影像数据，这里的05为模拟产生的数据； 2表示第二幅影像数据，这里的05wuhan1为实际的数据。 其实，上述二者的顺序是可以调换的。 模型运行后得出的结果为 里面有不同的评价指标值。