上回说到,对于空间统计学里面的分析,主要面对的数据一般都是面状数据。嗯,说到这里,大家回想(想不到的,可以去回看一下),整个儿新版白话空间统计,到现在为止第18期,基本上讲的分析例子,用的都是面状要素。(那么,插一个问题:以后会不会有点要素的分析呢?那当然是有的,对于点数据分析,是有专门的分析方法的,我们挖坑待填)要说空间关系的直观体验,做得最好自然是GeoDa,(关于GeoDa的介绍略)关于GeoDa软件的下载,可以在公众号发送消息"geoda",获取X度云地址……(吐槽一句,geoda的更新太慢了,现在依然是2019年8月发布的1.14版本……另外,源码最后也是更新到2019年11月,不知道是不是因为疫情导致的。)话说,作为专门的分析软件,而不是一个大而全的平台,GeoDa还是简单好用——几乎不用预先学习,也不用看啥操作手册,只要懂一点空间统计学的同学,能看懂他菜单上的名词的,就可以直接上手操作了(难道这就是GeoDa到现在为止,也没有一个比较正式的的操作手册的原因?)具体的GeoDa的操作,大家自己简单摸索一下,就能学会,灰常简单,作为一个专业的空间统计工具软件,不存在码农和产品经理那么多的脑洞,设计者的思想很明确:分析工具,作为工具,一个按钮就要有一个按钮的作用,憋整那么多关联控制……(手动狗头包括ArcGIS在内的一堆软件,一个按钮一个菜单的点击下去,ArcGIS甚至可以下沉到19级界面功能卡片……)GeoDa一打开,就要求你先加载一份数据,我们这里就直接选择Shape file了:从这里就看出分析工具和软件平台的不同了,GeoDa一上来就要求你得打开一个数据,因为它认为你要做分析,不打开数据算神马?好吧,你说了算。简明扼要的W,点看之后,会发现也是一个非常简明扼要的空间权重管理选项卡:然后需要选择一个ID(唯一值字段),如果没有的话,需要“添加ID变量”。老版本只能用ID,也就是只能用整数,在新版本里面,可以用字符型了,但是只能用字母,不能用中文字符,所以我们这里只能选择CODE:面要素默认的空间邻接模式,就是Queen,所谓的共点共边即相邻,然后直接点击创建,会弹出一个保存框,让你输入你要创建的空间权重矩阵的文件:而在中国地图,使用queen的话,会出现独立要素(海南、台湾),会弹出一个警告,直接点确定就完:接下去,就是比ArcGIS好用的地方了——GeoDa自带空间权重的可视化功能,大家注意下面三个按钮:这里再给出一个好评的是GeoDa可以直接实时联动,比如我选择最多的这个粉红色的柱子:表示周边有四-五个邻接的,在地图上就会直接显示出来,如下:2、连通性地图,可以直接查看那些要素相互连接,比如我们选择湖北,上面会非常直观的列出湖北省的邻接:那么看到这里,肯定有同学就会问,台湾没有数据,没法加入计算我能忍,but海南不能忍啊,作为中国疆域面积最大的省(海南省,包括海南岛和南海海域,加起来是中国疆域(非陆地面积)最大的省),不加入计算,这个不能忍啊,有办法破解么?下面我们就来看看GeoDa的空间权重矩阵文件,这里有个好消息,就是Geoda的空间权重矩阵文件是ASCII码,我们可以直接打开、读取和修改,为了保险,我先拷贝了一个副本(这个好习惯,希望大家都能养成):然后我们用文本编辑器打开,就可以看见这个gal文件,采用的是稀疏矩阵的格式存储的,第一行是元数据,然后从第二行开始,两行为一组:表示,编号为110000的要素(北京),有两个邻接,分别是120000(天津)和130000(河北)那么我们来做个简单的修改,比如我让北京和山东也有关系,只需要改成这样:然后在空间权重管理里面,重新去加载这个文件。删除掉原来的空间权重参数,然后点击打开:就发现北京和山东变成邻接了,但是选山东,却没有北京,这是因为我们只改了北京,没有修改山东,所以在空间权重矩阵里面,非特殊情况,都是无向图,所以要成对的去修改。改北京和山东,没啥意义啊,我们下面来修改一下海南,我们设定海南与广东广西两省作为邻居,原始情况是这样的:460000的海南,邻居是0,我们来修改一下,把0改成2,下面加上广东和广西的ID:然后同样的,要回去把广东广西修改一下,挂接上海南:注意,修改的时候,要先改后面的邻接数量,再添加信息,否则就会不成功保持之后,重新加载,选择海南的时候,就有广东广西两个省了下一章我们用Geoda来讲讲,不同空间关系对莫兰指数的影响。
|