地理加权步骤

1.数据准备明确因变量和解释变量指标(有机碳密度是因变量，温度、降水量、粘粒含量是解释变量指标)，准备好每个样点的因变量及解释变量指标数据，必须是数值。2.回归分析将准备好的数据导入spss（一般可以直接在spss里面打开）做逐步回归，选择进入的指标作为解释变量，温度、降水量、粘粒含量三个指标全部进入。图13.地理加权回归将带有因变量和解释变量数据的点shp文件导入ArcGIS。打开工具箱-空间统计工具-空间关系建模-地理加权回归（图2）。（1）输入要素添加shp文件；（2）因变量选择作为因变量的字段；（3）解释变量选择作为解释变量的字段；（4）选择一个放置“加权点”的位置并命名，生成的加权点包含点数据的回归系数、预测值和残差等；（5）核类型、带宽方式等如果无特殊要求，默认即可；（6）打开附加参数选项，在系数栅格工作空间中创建一个新的数据库并命名（图3），用以存储解释变量的系数栅格以及截距栅格；（7）像元大小选择合适的栅格大小（旱作区选择2000感觉还可以，可自行调试）；（8）打开环境设置选项（图4），在范围处理中选择底图图层（否则默认范围为点shp，一般比研究范围要小一些）。图2图3图4操作完成后运行一段时间，时间大小与ArcGIS版本有关，运行完毕得到解释变量的系数栅格以及截距栅格（图5）。图54.裁剪打开工具箱-数据管理工具-栅格-栅格处理-裁剪（图6）。（1）输入需要裁剪的栅格；（2）输入底图文件（这个图就是我们最后想要得到的范围）；（3）勾选“使用输入要素裁剪几何”，一定要勾选；（4）在输出栅格数据集处选择存储裁剪文件的位置并命名;（5）其他全部默认，点击确定。图6得到解释变量和截距的系数栅格图，分析其与有机碳密度的相关关系。