地理加权回归（GWR）步骤：ArcGIS+GWR4+GeoDa实现

0.数据准备

01-21平均NPP栅格、01-21平均年降水栅格、01-21平均年均温栅格、01-21测区内（省级）各地市GDP shp、01-21各地市人口密度shp

1.创建渔网

将研究区域矢量数据划分为网格（渔网）数据（要选择创建标注点），将测区（省级）划分为10 km ×10 km 的网格，并按需要对测区域进行提取。

先看（整个流程）：

https://blog.csdn.net/weixin_47779247/article/details/127473294

再看（细节处理）：ArcGIS创建渔网步骤及一些问题-CSDN博客

2.GDP、人口密度shp数据插值并转为矢量

（1）将面要素的行政区shp图层通过要素转点转换成点要素图层，生成的点为各行政区单元的几何中心，并且生成的点要素均带有面要素图层的属性。

（2）将点要素图层运用反距离权重法进行插值（根据采样点值创建连续表面，反距离权重法是一种确定性插值方法，其原理是将根据已知周围采样点的测量值和用于确定所生成表面平滑度的指定数学公式将值指定给其他所有位置，使数值铺满整个表面），生成栅格表面。

具体操作见视频mp.weixin.qq.com/s/cZVF9RvSGM2HNr9m8ftIaw

3.将要分析的NPP、降水、温度、GDP（插值后的栅格）、人口密度（插值后的栅格）等栅格影像多值提取至点

具体操作：

基于ArcGIS和SPSS对两幅或者多幅遥感影像进行回归分析建模_spss可以处理遥感数据吗-CSDN博客

4.将提取至点的矢量数据进行数据筛选

（1）打开属性表，剔除异常数据（或者将正常数据导出）

我用的是将正常数据导出，分析工具-提取分析-筛选，表达式如下（属性表中我发现异常数据均是-9999）

（2）为每个点创建唯一ID（FID会随数据变动，所以需要再创建一个字段），并计算每个点经纬度

首先给属性表创建三个字段：

创建ID：打开属性表-添加字段-长整型；

longitude：打开属性表-添加字段-双精度；

latitude：打开属性表-添加字段-双精度。

其次，对于ID，右键字段里的ID，打开字段计算器，输入以下即可：

然后，对于X/Y，右键字段，计算几何：

单位十进制度是经纬度，米是地理坐标系。

（3）表转Excel

5.用ArcGIS做OLS普通线性回归（操作看我的图，记得输出报表）

OLS和GWR的结果解读看这里：

https://mp.weixin.qq.com/s/BtFqGsc4ogEeGuh6iNoBIQ

6.下载GWR4软件，进行GWR

完全按这个来：

GWR4软件怎么用+结果解读+结果在ArcGIS中可视化-CSDN博客

7.将点数据（我的是Huigui1.shp）导入GeoDa中进行双变量Moran’s I分析，得到I值、p值、z值

步骤：

（1）创建权重文件

（2）双变量Moran’s I

步骤：空间分析-双变量局部Moran’s I

 

 做完上面把上面的图存起来，最上面两张LISA显著性地图、LISA聚类地图可以加背景地图影像。

Moran’s I 结果解读：https://mp.weixin.qq.com/s/LlOEGbe-iK7lKbI2xDMbzA

欢迎联系讨论！

Best wishes！

                                                                                                                                    Gwy

                                                                                                                             2024年3月15日