GeoDa空间计量（一）——空间权重矩阵的生成

一、基于地理邻近的空间权重矩阵

（一）导入shp文件，文件打开后展示为地图形状；

（二）创建权重矩阵。

1、设置各省市的ID变量。在这里需要注意，如果shp文件中包含POLYID变量，即每个省市都有一个ID，则不需要创建POLYID变量；否则，需要创建POLYID变量（点击“Add ID Variable”,在弹出的窗口中会自动出现“POLY_ID”这个变量名，点击‘Add Variable’，创建成功）。这里我们已经包含了POLYID变量，故不需要创建。

2、设置相邻的类型，这里采用“Queen contiguity”。注意，“Queen contiguity”是两个省市有公共边界或公共点则相邻，为1；否则为0。“Rook contiguity”是两个省市有公共边界则相邻，为1；否则为0。“Order of contiguity”是生成的权重矩阵为一阶或多阶；一阶指只将直接相邻的省市记为1；二阶指不仅将直接相邻的省市记为1，相邻省市的相邻省市也记为1；一次类推。一般情况下，我们均采用一阶空间权重矩阵；且因中国不存在省市之间以点相邻的情况，故采用“Queen contiguity”和“Rook contiguity”生成的权重矩阵相同。

3、保存权重矩阵，这里重命名为“rook1.gal”


**（三）权重矩阵的解读。**这里使用‘Notepad++’打开‘rook1.gal’文件。下同。

第一行为各省市的ID号和shp文件名称；第二行为ID为1的省市和它的邻居个数；第三行为ID号为1的省市的邻居的ID号；第四行为ID号为2的省市和他的邻居个数；第五行为ID号为2的省市的邻居的ID号；以此类推。

二、基于地理距离的空间权重矩阵

半径距离空间权重矩阵

前面的操作和基于地理邻接的空间权重矩阵操作相同，不同的是在（二）下面的2——参数设置。本文的设置为：



（三）权重矩阵的解读。

第一行同‘rook1.gal’解释相同；第二行为ID为1的省市与ID为28的省市相邻，距离为758370.096；第三行为ID为1的省市和ID为4的省市相邻，距离为 434322.347；第四行为ID为2的省市和ID为10的省市相邻，距离为875071.467；以此类推。注意：这里只展示了相邻的情况，不相邻的情况没有展示。

具体参数解释：
1、点击“Threshold distance”，默认如下。

2、点击‘Distance metric’右边的向下箭头，出现如下三个选项。其中，‘Euclidean Distance’是在地图为投影的情况下选取的测量距离的方式；其余两种是在地图为地理坐标的情况下选取测量距离的方式，‘mi’为英里（mile）,‘km’为千米（kilometer）。

3、X和Y坐标系变量默认为X轴的核心和Y轴的核心，即某省市的重心。软件自带，不需要设置，了解即可。

4、‘Threshold distance’是最小门槛距离，可通过鼠标滑动改变大小。

k个邻近空间权重矩阵

前面的操作和基于地理邻接的空间权重矩阵操作相同，不同的是在（二）下面的2——参数设置。本文的设置为：

‘Number of neighbors’为邻居个数，这里默认为4，也可根据情况自行修改。


（三）权重矩阵的解读

具体解读同半距离空间权重矩阵相同，不同的是K个邻近空间权重矩阵每个省市的邻居都相同，半距离空间权重矩阵各个省市的邻居数不同。