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