地理加权步骤

地理加权步週OYEAR-YICAIEJINGBIAN

1.数据准备

明确因变量和解释变量指标（有机碳密度是因变量，温度、降水量、粘粒含量是解释变量指标），准备好每个样点的因变量及解释变量指标数据，必须是数值。

2.回归分析

将准备好的数据导入SPSS （—般可以直接在SPSS里面打开）做逐步回归，选择进入的指标作为解释变量，温度、降水量、粘粒含量三个指标全部进入。

t Sig.

G标淮谏差试用版

1 （常列

a z ave 1 121

148

072

.015.444

16.458

9.904

.000

.000

2 （常呈）

平均气温1.706

132

.129

015

008

.397

-.239

13.249

8.97B

-541 B

.000

000

000

3

a^_ave 平

沟气温

'm 1 506

125

・057

001

149

015

009

000

.375

-.307

.131

10.124

8.391

・ 6.034

2.621

000

000

009图1

3・地理加权回归

将带有因变量和解释变量数据的点Shp文件导入ArcGIS o打开工具箱■空间统计工具■空间关系建模■地理加权回归（图2）。

（1）输入要素添加shp文件；

（2）因变量选择作为因变量的字段；

（3）解释变量选择作为解释变量的字段；

（4）选择一个放置“加权点”的位置并命名，生成的加权点包含点数据的回归系数、预测值和残差等；

（5）核类型、带宽方式等如果无特殊要求，默认即可；

（6）打开附加参数选项，在系数栅格工作空间中创建一个新的数据库并命名（图3）,用以存储解释变量的系数栅格以及截距栅格；

（7）像元大小选择合适的栅格大小（旱作区选择2000感觉还可以，可自行调试）；

（8）打开环境设置选项（图4）,在范围处理中选择底图图层（否则默认

范围为点shp, 一般比研究范围要小一些）。

图2

图3

操作完成后运行一段时间，时间大小与ArcGIS版本有关，运行完毕得到解释变量的系数栅格以及截距栅格（图5）。____________________

图5

4.裁剪

打开工具箱■数据管理工具■栅格■栅格处理■裁剪（图6） o

（1）输入需要裁剪的栅格；

（2）输入底图文件（这个图就是我们最后想要得到的范ffl）:

（3）勾选“使用输入要素裁剪儿何S —定要勾选；

（4）在输岀栅格数据集处选择存储裁剪文件的位置并命名；

（5）其他全部默认，点击确定。

图6

得到解释变量和截距的系数栅格图，分析其与有机碳密度的相关关系。