生成矩阵_空间权重矩阵生成与使用

1 空间权重矩阵生成

1.1 准备工作

• ArcMap

• STATE_poly.shp

• Stata 15

1.2 软件操作

step 1 新建文档：

新建 —— 空白地图 —— 确定

step 2 加载地图数据：

[1] 图层 —— 右击 —— 添加数据 —— {名称：路径-shp文件} —— 添加

[2] 或者 文件夹链接 —— 所在文件夹 —— 文件 —— 拖进来

step 3 生成id字段：

STATE_poly.shp —— 右击打开属性表 —— 表选项 —— 添加字段 —— { 名称： id} —— 确定

系统工具箱 —— Data Management Tools.tbx —— 字段 —— 计算字段 —— { 输入表：STATE_poly —— 字段名：id —— 表达式：[FID] } —— 确定

step 4 生成空间权重矩阵：

系统工具箱 —— Spatial Statistics Tools.tbx —— 空间关系建模 —— 生成空间权重矩阵 —— {输入要素：选择STATE_poly.shp文件 —— 唯一ID段：id —— 输出：myfile.swm —— 空间关系：CONTIGUITY_EDGES_CORNERS —— 行标准化：【可选(本演示未选)】} —— 确定

【因为gif图有点大没有上传成】

step 5 转换空间权重矩阵为表：

 系统工具箱 —— Spatial Statistics Tools.tbx —— 工具 —— 将空间权重矩阵转换为表 —— {输入空间权重矩阵文件：myfile.swm —— 输出表：路径 + myfile } —— 确定

【因为gif图有点大没有上传成】

step 6 导入stata：

import dbase using "C:\Users\wuxiong\Desktop\myfile.dbf",clearformat %2.0f WEIGHTreshape wide WEIGHT ,i(ID ) j(NID )
mvencode WEIGHT* , mv(.=0)

2 空间权重矩阵使用

spmat dta W WEIGHT*, normali