土地利用覆被变化时空信息分析方法及应用47486.pdfVIP

第 4 章 土地利用/覆被变化模型模拟方法 第4 章 4 土地利用/覆被变化模型模拟方法 4.1 马尔可夫模型 4.1.1 模型介绍 马尔可夫（Markov ）法是以俄国数学家A.A.Markov 名字命名的一种方法，他是基于马 尔可夫过程，通过对事物不同状态的初步概率和状态之间转移概率的研究，确定状态变化趋 势，以预测事物的未来( 吕效国，2009) 。 马尔可夫过程是一种具有“无后效性” 的特殊随机过程。所谓“无后效性”是指过程（或系 统）在时刻 t0 所处的状态为已知的条件下，过程在时刻 t t0 所处状态的条件分布与过程在 时刻 t0 之前所处的状态无关，即过程“将来” 的情况与“过去” 的情况是无关的。这对于研究土 地利用的动态变化较为适宜，因为在一定条件下，土地利用的动态演变具有马尔可夫过程 性质： （1）在一定区域内，不同土地利用类型之间具有相互可转化性。 （2 ）土地利用类型之间的相互转化过程包含着较多尚难用函数关系准确描述的事件。 运用马尔可夫模型关键在于确定土地利用类型之间相互转化的初始转移概率矩阵 P ， 是由一步转移概率pij 所组成的矩阵，状态空间I {1,2,}，称 p11 p12  p1n  P p21 p22  p2n          且满足以下两个条件 （1）pij  0 ，i,j I （2 ） pij 1，i I j I 式中 n 为研究区土地利用类型的数目，pij 为初始到末期由类型 I 转移为类型 J 的概率。 根据马尔可夫模型，系统所研究的事务在任何时刻的状态概率向量(n)，可由它的初始 状态矩阵(n  1)和转移概率矩阵pij 来确定，可得公式（张军等，2006 ）： (n) (n  1)ij 马尔可夫模型已经在土地利用变化的研究中得到了广泛的应用，比较典型的案例有，张 春华等以三江源地区 1995 年，2000 和 2005 年 3 个时间点的土地利用数据为基础，运用马 -157- 土地利用/覆被变化时空信息分析方法及应用 尔可夫法预测未来 10 年内该区土地变化趋势，为该区土地管理决策、生态环境保护提供参 考（张春华等，2009 ）。朱懋等依据德化县2002~2005 年土地利用变更调查数据，采用马尔 可夫模型对其土地利用结构变化进行了模拟与检验，用此方法对 2010 年和 2020 年德化县土 地利用结构的演变趋势进行预测，并根据预测结果对今后土地利用和管理提出建议（朱懋等， 2009 ）。迟斌等根据2004~2009 年重庆市土地利用数据，按年份进行加权处理，得出转移概 率矩阵，运用加权马尔可夫法预测重庆土地利用资源的变化，对重庆市未来土地利用规划起 到了指导作用（迟斌等，2011 ）。孟庆香等以黄土高原 1986~1997 年土地利用变化为研究对 象，并运用马尔可夫法对未来土地利用进行预测，为其后进行的生态效应研究提供数据依据 （孟庆香等，2009 ）。 4.1.2 模型应用 （1）数据准备与处理 本应用以京津冀为例，开展基于马尔可夫模型的土地利用模拟预测。京津冀即京津冀都 市圈，是国家“十一五”规划中的一个重要的区域规划，