马尔科夫模型预测景观动态变化.docVIP

基于RS、GIS的马尔克夫模型在景观格局变化分析 实验原理?? 景观演变是一个十分复杂的过程，自然过程、人类活动一直在影响着湿地景观格局的发展变化。特别是在人口不断增长、社会经济结构在不断调整变化的今天，人们对景观不同的利用方式、不同的认识都在改变景观格局。马尔科夫过程较多地用于模拟和预测自然植被的变化，，因此我们将学习马尔科夫模型模拟预测景观格局的动态变化。 ????? 马尔科夫过程是一种特殊的随机运动过程。一个运动系统在t+1时刻的状态和t时刻的状态有关，而和以前的状态无关，利用这一点进行湿地景观格局变化的预测是适合的，运用马尔科夫模型的关键在于景观转移矩阵和转移概率的确定。 土地利用变化的转移矩阵，是一种表格，可以表示不同时期各种土地利用类型的面积变化, 流入率、流出率等,定量分析土地利用变化的强度、范围和方向。转移矩阵说明不同时间段内同一区域内土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。 比如某一类别的土地有多少面积分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的。 时期1 时期2 类1 类2 类3 合计 类1 43 5 2 50 类2 2 45 3 50 类3 0 1 49 50 合计 45 51 44 150 ?????二 实验方法 本次实验利用两期遥感数据分类，获得分类结果，然后以年为单位，把景观格局变化分成一系列离散的演化状态，通过各时间段内景观要素的年平均转化率获得从一个状态到另一个状态的转化速率（即转移率，步长为1a）。各类景观要素转移率确定后就可以构成一个转移率矩阵，其数学表达式为： ????? Pik为景观要素i转化为j的转移概率，其特点有：①0≤Pij≤1；②，即每行元素之和为1， n为景观要素的数目。 根据马尔科夫模型模拟过程的基本方程： ?????????? n：状态转移时段或时期（本文中表示要预测的某年与起始年差值）；????? Pij(n)：转移率矩阵元素，表示事件在初始（k＝0）状态为已知的条件下，经过n次状态转移后，在第n个时刻（时期）由状态Pi向状态Pj的状态转移率（本文中指景观要素经过n年后的状态转移率）；????? k=0，1，2，3，……n-1；（表示1991年到要预测的某年间）。? ?代入初始状态矩阵和状态转移概率矩阵即可得到一个马尔柯夫链模型,求得起始年间以后任一年各土地利用类型的转移概率j ,进而计算各土地利用类型的面积百分比及数量需求变化。 三、实验步骤 1利用卫星遥感影像图进行矢量化，获得1991年和1999年各景观要素的分类图，计算各类别总面积； 时间 1991 1999 2 计算初始状态矩阵 初始状态矩阵就是1991年各种类别的面积矩阵： 类别 COUNT area（平方公里） 1 88215 79.3935 2 231606 208445400 3 268666 241799400 4 12446 5 532573 479315700 6 885807 797226300 7 1681998 1513798200 8 355281 319752900 9 433167 389850300 10 26123 11 675441 607896900 ?3????计算状态转移概率矩阵： 先在ARCGIS中计算1991~1999年间的景观类别的面积转移矩阵，利用栅格数据分析工具中的区域分析/面积制表工具，制作二维表 。 ?平均每年的转移概率矩阵：?由于1991~1999年间相差8年，所以可由各景观要素面积的转化数值除以8，求出湿地各景观要素面积的年平均转化值（km2/a），再用湿地各景观要素的年平均转化值除以1996年各景观要素总面积，得到各湿地景观要素的状态转移概率矩阵（步长为1a）（表2-3）。 4 预测景观动态变化 利用excel表模拟：根据马尔科夫模型的模拟，?代入初始状态矩阵和状态转移概率矩阵即可得到一个马尔柯夫链模型,求得起始年间以后任一年各土地利用类型的转移概率j ,进而计算各土地利用类型的面积百分比及数量需求变化。从而预测2010、2020年各景观要素的面积。 （例）表2-4 2010年-2050年湿地各景观要素的面积 四、结果分析与讨论 要求列出结果并分析： 初始状态矩阵：（表格） 景观平均年变化概率矩阵： 2015、2020年各景观面积： ????? 利用地理信息系统软件对卫星遥感影像图分类结果进行分析，获得不同时段的景观要素面积数据，运用马尔科夫模型以不同时期各景