dbf文件怎么创建_如何创建空间权重矩阵

点击上方蓝字，关注公众号

更多空间计量推文：

超好用的空间数据分析软件GeoDa(点击阅读)

ArcGIS地图制图(点击阅读)

ArcGIS矢量数据空间分析(点击阅读)

ArcGIS栅格数据空间分析(点击阅读)

ArcGIS空间统计分析(点击阅读)

Stata空间计量全面教程(点击阅读)

Stata空间计量问答精选(点击阅读)

MATLAB空间计量实战(点击阅读)

SPSS地图制图(点击阅读)

Excel地图制图(点击阅读)

R地图制图(点击阅读)

空间

世界上万千事物的状态都可以由一个三维的空间坐标系与一个一维的时间坐标系来唯一刻画。时间或空间上距离相近的两个事物的状态是相互关联的，即不能被认为是相互独立的，且两事物越是接近，它们状态的相关性越强。当两点距离为零，实则是同一个体，它们将完全相关。越是相距遥远的事物相关性越弱，当两事物之间距离为无穷远，可近似地认为两者完全不相关。

空间计量即旨在为处理截面数据或面板数据中的空间效应、空间相关性与空间异质性发展专门的建模、估计与统计检验方法。对于空间计量模型而言，权重矩阵的元素是非随机的、外生的。创建空间权重矩阵，是进行空间计量分析的重要一步。

确定空间权重矩阵，常用的方法有：

1. 邻接：假定空间相互作用的强度，取决于地区间的相邻关系。例如，0-1邻接矩阵，当两地区具有公共边界，则定义权重为1，否则为0；但是0-1邻接矩阵是一个对称阵，不符合某些现实情况(模仿效应)，且0-1元素的设置无法区分各邻居空间作用的强弱。可选择的改进方法是，以“两地区边界相同部分的长度”与“该地区与其他相邻接的个体边界的总长”之比确定权重。

2. 距离：假定空间相互作用的强度，取决于地区间的质心距离或者区域行政中心所在地之间的距离。例如，当两地区在定义的距离之内，即两地区相邻，则定义权重为1，否则为0；在这种情况下，不同的权重指标随距离的定义而变化，其取值取决于选定的函数形式，如距离的倒数或倒数的平方、欧氏距离等。当然，还需要定义一个门槛距离，超过了某给定的门槛距离则地区间的相互作用可以忽略不计。

3. 经济社会流量/社会网络结构：除了使用真实的地理坐标计算地理距离外，还有包括经济和社会因素的更加复杂的权重矩阵设定方法。例如，根据地区间交通运输流、通讯量、GDP总额、贸易流动、资本流动、人口迁移、劳动力流等确定空间权重，计算各个地区任何两个变量之间的距离。

关于各种权重矩阵的选择，没有现成的理论根据，一般可考虑空间计量模型对各种空间权重矩阵的适用程度，检验估计结果对权重矩阵的敏感性，最终的依据实际上就是结果的客观性和科学性。

一般而言，尽管二进制邻近矩阵并非适用于所有的空间计量模型，但是由于多数情况下的实用性，在构建空间计量模型时，首选就是从二进制的邻近矩阵开始，逐步选择、确定空间权重矩阵。

以下分别介绍使用GeoDa、ArcGIS、MATLAB创建空间权重矩阵的方法，以我国大陆地级市地理信息作为演示数据：

GeoDa

在GeoDa中，主要有3种方式来产生空间权重：

一是有公共边界或顶点即为相邻，分“Queen contiguity(后相邻)”和“Rook contiguity(车相邻)”两种情况，“后”、“车”指国际象棋走步法；

二是在一定距离阈值范围内为相邻；

三是将距离最近的k个单元设为相邻。这3种方式下生成的空间权重矩阵均为0-1矩阵，即两单元相邻则空间权重矩阵的元素赋值为1，否则为0；

通过“工具-空间权重管理”菜单，或点击“W”图标，均可创建空间权重矩阵，并生成.gal格式(邻接)文件或.gwt格式(距离)文件：

GeoDa 1.12操作界面

ArcGIS

在ArcGIS中，点击工具栏中的ArcToolbox图标，选择工具箱中的“空间统计工具-空间关系建模-生成空间权重矩阵”，依次设置；

注意，“空间关系的概念化”下拉菜单中，选项从上到下依次为：距离倒数、固定距离、k个最近邻、仅有边界相邻、边界和顶点相邻；“距离法(可选)” 下拉菜单中，选项默认为欧几里得距离。

通过该工具，可创建.swm格式文件，该格式文件仅能在ArcGIS中使用，如需要表格文件，可使用ArcGIS将其转换为.dbf格式文件：

ArcMap Desktop 10.2操作界面

MATLAB

在MATLAB中，首先要加载分析所须的distance_weights包：该程序包由Donald J. Lacombe编写。

找到MATLAB安装路径中的toolbox文件夹，例如我的电脑是C:\Program Files\MATLAB\R2012a\toolbox，把distance_weights文件夹复制进去，然后在MATLAB中点击file-set path-add with subfolders，选择distance_weights以设置路径。

执行以下命令，然后打开workspace中的result，然后找到dw，点击dw后可以得到基于距离的空间权重矩阵；若创建基于反距离的空间权重矩阵，则须使用invdistance_wm命令：

clear;clc;

% 读取数据

[data,text]=xlsread('中国地级市经纬度数据.xlsx')

xc=data(:,1)

yc=data(:,2)

% 创建基于距离的空间权重矩阵

result=distance_wm(yc,xc)

result.dw

MATLAB R2012a操作界面

阿虎定量笔记

数据分析方法

软件使用技巧