基于OLS和GWR模型的区域土地适宜性模拟研究.docVIP

基于OLS和GWR模型的区域土地适宜性模拟研究 　　【摘 要】本文以合肥市包河区为研究区，考虑驱动因子对土地利用适宜性的影响，用全局最小二乘法OLS Logistic回归模型和GWR Logistic模型来建立土地利用格局与其驱动因子之间的回归模型，并进行土地利用格局模拟的实例研究与模型比较。同时，运用此模型定量地分析影响区域土地利用变化的内在机制，预测土地利用未来变化的趋势，可以为包河区及其类似地区的土地利用规划决策提供科学依据。 　　【关键词】土地适宜性；地理加权回归；模拟 　　土地利用/覆被变化（Land Use/Cover Change，简称LUCC）研究是对土地利用的数量、质量和空间布局等变化的综合分析。目前流行的方式是在进行土地利用变化驱动力建模时，运用基于全局最小二乘法（Ordinary Least Squares，简称OLS）的回归分析方法，对各因子进行全局的参数估计。当考虑某一驱动因子在不同空间位置对土地利用格局的影响程度不同时，常采用地理加权回归（Gergraphically Weighted Regression，简称GWR）模型。该模型是对传统回归分析的扩展，允许局部而不是全局的参数估计，通过附加表达空间对象本身相关性和异质性的变化参数，反映样本对回归方程贡献在空间上的分异，使回归结果更加可信。两种方法适用情况不同，在地理要素的空间分析和建模中都有广泛应用。 　　据此，本文以合肥市包河区为例，考虑驱动因子对土地利用格局的影响，采用OLS和GWR模型来建立土地利用格局与其驱动因子的回归模型，并进行验证和分析。 　　1 研究区与数据准备 　　1.1 研究区概况 　　包河区位于合肥东南迎风口，襟五河（包河、南淝河、十五里河、塘西河、派河）而带“一湖”（全国五大淡水湖之一的巢湖）。全区区域面积340平方公里（其中巢湖水面面积70平方公里），辖7个街道、2个镇、1个省级工业园区，常驻人口81万，是合肥市四个城区中面积第一大区、人口第一大区。 　　包河区属亚热带湿润季风气候年平均气温15～16℃，年日照时数2100H ，多年平均降雨量900～1000mm，且时空分布不均，年内和年际变化很大，最大年降雨量与最小年降雨量相差2～3倍，梅雨显著，夏雨集中。地形、地貌及生态环境状况具典型的江淮分水岭特征，表现为丘陵地貌，地形破碎，干旱缺水，水土流失严重。随着经济发展、人口持续增长和省会城市的辐射影响，城镇和农村建设用地快速扩张，耕地锐减，已严重影响了区域粮食安全和生态系统的稳定。 　　1.2 数据准备 　　本研究所使用的土地利用数据来源于包河区2005年土地利用变更调查数据，参考国际粮农组织（Federal Agriculture Organization， 简称FAQ）的土地覆盖分类系统（Land Cover Classification System， 简称LCCS），考虑包河区土地利用特点，结合TM遥感影像对矢量的土地利用现状图进行重分类，得到农业用地、林业用地和城市用地3种土地利用类型。同时，研究利用包河区数字高程模型、自然环境、经济社会等数据，选取坡度、离水系距离这两个自然环境因子和离交通距离、人口密度这两个社会经济因子，作为包河区土地利用空间格局形成和演变的主要驱动因子，应用GIS空间分析的方法对其进行空间量化。 　　2 研究方法 　　土地利用变化驱动力建模常采用对评价因子采取某种或多种评价方法，从而获取评价结果的方式。对于评价方法，可以有多种多样的选择，如图1。 　　本研究采用线性关系，将全局OLS Logistic回归模型和GWR模型应用于土地利用格局模拟的回归建模中，考虑驱动因子对土地利用格局影响的空间不稳定性，对两者进行分析与比较。具体来说，分以下几步： 　　Step1：根据已有的2005年包河区土地利用类型，并采用重分类将各种不同类型用地分开。 　　Step2：数据预处理的另一部分工作是将驱动因子空间量化，包括距离交通、人口密度、距离水系和坡度四类因子。 　　Step3：将上述土地利用类型图和驱动因子空间量化图栅格化，并将每个栅格图层中的所有栅格单元值转换为文本文件，并导入Excel中。 　　Step4：将所得数据导入SPSS 17.0进行OLS Logistic逐步回归和ROC曲线图。 　　Step5：将土地利用类型图层和驱动因子图层导入ArcMap 9.3中，利用工具分别对其进行OLS和GWR回归分析。 　　3 结果与分析 　　3.1 定性分析结果 　　基于前述理论，首先在ArcMap中应用工具对土地利用类型和驱动因子图层定性的进行OLS Logistic回归和GWR分析，并与2005年现状对比，如图2所示。 　　定性分析中，对2005年各类型土