绘制手机基站泰森多边形及建立邻接矩阵
最近拿到了手机信令(CDR)数据,准备用于做一些科研研究,记录一下怎么用ArcGIS绘制泰森多边形,并建立邻接矩阵。
excel数据转为点要素数据
我们拿到的数据往往是文本数据或者csv数据,以excel数据为例。如果需要在ArcGIS中使用上述数据,首先需要将其转化为shp格式的数据。
- 在ArcMap中导入excel数据
先点击左上角导入按钮,后根据提示选择需要导入的excel数据。注意可能有的版本不支持.xlsx格式,需要转换为.xls格式。
- 将导入的数据转换为shp文件
右击左侧图层中导入的sheet数据,选择“显示XY数据”,弹出如下对话框
根据具体内容选择XYZ字段(X为经度,Y为纬度,Z为高程),点击确定,即可将点打到地图上
完成后,右击左侧图层中的sheet,选择数据–导出,编辑好名称点击确定,即可导出shp数据。
根据点元素数据建立泰森多边形
将excel数据转化为shp数据后,使用ArcGIS中的工具箱,输入元素选择刚才获得的基站shp数据,即可快速完成泰森多边形的建立。
完成后的效果如下:
刚完成图形建立的时候,整个图形会是矩形,需要使用裁剪工具,将边缘不需要的点裁剪删除掉。
建立邻接矩阵
完成泰森多边形的建立后,考虑获取泰森多边形的邻接矩阵,使用ArcGIS工具箱的命令完成。
使用上图中的空间连接,双击后弹出以下对话框
目标要素和连接要素都填刚才得到的泰森多边形,同时点击上图右下角的+号按钮,并编辑如下
字段名称可自行修改,属性中的长度可根据实际情况更改,分隔符为英文逗号,完成后确定。
完成后连接要素的字段映射中出现刚才新建的字段名称(此处为邻接),右击该字段选择添加输入字段,之后选择第一个尾缀为.Input_FID的字段(该字段为该泰森多边形的中心点的字段,以该字段建立映射,方便后续完成与基站点相关的操作),点击确定,即可得到邻接矩阵。
右击左侧图层中的邻接矩阵文件,选择打开属性表,如上图,可发现有一列列名为邻接,该列即为该面的邻接面的Input_FID;
选择需要数据,导出至excel后,使用python进行简单处理,即可得到邻接矩阵。