第29卷第2期2010年02月地理科学进展PROGRESSINGEOGRAPHYVol.29,No.2Feb.,2010收稿日期：2009-08；修订日期：2009-11.基金项目：国家自然科学基金项目(70873118、70821140353)、国家科技支撑项目(2006BAC08B03、2006BAC08B06、2008BAK47B02)。作者简介：邓祥征(1971-)，男，博士，副研究员，主要从事区域环境变化、土地系统变化及效应研究。E-mail:dengxz.ccap@igsnrr.ac.cn179~185页基于栅格面积成分数据的土地利用格局解释模型稳健估计邓祥征1，林英志1，战金艳2，何书金1(1.中国科学院地理科学与资源研究所，北京100101；2.北京师范大学环境学院水环境模拟国家重点实验室，北京100875)摘要：针对最小二乘估计不能应用于栅格尺度以面积成分表征的土地利用格局驱动机理分析的难题，本文提出了一种利用偏最小二乘回归法稳健估计该类土地利用格局解释模型的方法。利用该方法能在解释变量间存在多重共线性的情况下，获得基于栅格面积成分数据的土地利用格局解释模型的稳健估计。本文推导了应用偏最小二乘回归分析的数据处理和建模估计过程，并运用该方法开展了针对黄淮海地区耕地、建设用地分布格局及其驱动因子的建模分析，得到了拟合优度高的估计结果。研究表明，偏最小二乘回归分析方法在开展栅格尺度以面积成分表征的土地利用格局驱动机理分析时具备高效与稳健的特征，适宜在类似研究中推广应用。关键词：栅格成分数据；偏最小二乘回归；土地利用格局；驱动力；黄淮海地区1引言随着GIS空间数据分析技术的不断发展与大量空间数据集的积累，栅格数据逐渐成为地理学相关研究领域最主要、最基本的数据格式之一[1-3]。而基于栅格数据发展而来的栅格面积成分数据，在土地利用格局演替、土地利用结构变化模拟等领域发挥着日益重要的作用。虽然广受地学分支学科及其研究领域的重视，迄今为止，栅格面积成分数据的建模及估计方法仍未取得实际性进展[4-7]。当前，为了避免土地利用格局驱动机理分析中的共线性问题，常用的土地利用格局驱动机理回归分析方法仍主要依靠通过增减自变量提升模型估计结果的准确性。利用栅格面积成分数据估计土地利用分布解释模型统计建模也面临着两个技术难题：一是回归分析得到的栅格内各种土地利用类型面积之和常常不等于栅格总面积；二是变量间的多重共线性问题导致模型估计结果存在偏误[8]。针对上述问题，本文提出并推导了栅格尺度以面积成分表征的土地利用格局解释模型的偏最小二乘回归分析方法，并在黄淮海地区开展了案例研究，建立起土地利用格局解释模型，揭示了2000年各驱动因子对黄淮海地区土地利用格局的驱动机理，为土地利用格局演替、土地利用变化驱动机理分析等研究提供可资借鉴的方法。本文第二部分介绍栅格成分数据的概念和基于栅格面积成分数据进行偏最小二乘回归分析的步骤。第三部分针对黄淮海地区开展案例研究，验证栅格成分偏最小二乘回归分析方法的有效性与推广应用价值。最后一部分是结论与讨论。2模型与方法2.1栅格成分数据成分数据的概念最早由Ferrers于1866年提出[9]，经一个多世纪的发展，现已成为统计分析中的一种重要数据类型。通常，成分数据可用如下向量集的形式表达：S=(s1，s2，…，sm)T∈Rm|∑i=1msi=1,2
