毕业论文需要用到空间自相关,所以摸索摸索了好久,终于弄出了大概的流程了。
情景1
如果你没有shp格式的文件数据,
那么我建议你下载geoda095i这个版本(因为最新版本的我不太会操作)
明确问题:假如我们要对广东省各市2005人均GDP进行空间统计分析。
1.将数据输入为txt文本。
如上,第一行的意义是总共有21个观测点,共5列;第二行是标题行,意义分别是城市、唯一的ID号、维度、经度、人均gdp;剩下的21行为观测值的数据。
2.打开geoda095i软件(新版的不行哦),选择tools>shape>points from ASCII;
3.点击Inputfile行后的文件夹,在弹出窗口中通过路径选择gdp.txt文件。点击Outputfile行后的图标输入文件名:gdp。在X一coords栏中选择LON,在y-coords栏中选择LAT,点击Create完成。
这样就创建好一个shp格式的文件。
4.构建空间权重文件
我比较笨,所以这块儿构建空间权重矩阵用最新版的geoda了。
打开新版geoda,首先用file<-new project from
打开weights file creation框,点击add id variable...在existing
variables中选择ID项(也可以选其他的,但是其中包含的数值要唯一),如下图,
然后,点击add variable;点击 threshold distance,将distance metric中选择arc distance(km),其他为默认,从而创建出空间权重矩阵。如下图,
4.全局空间自相关分析和局部空间自相关分析全局,在新版geoda软件space<-univariate moran's I,在弹出的variable settings中选择first
variable(x)为GDP,点击ok.在弹出的select weight中的select from file选择刚才创建的空间权重矩阵gdp.gwt,点击ok.结果如下,
局部,在新版geoda软件 space<-univariate local moran’s I,在弹出的variable settings中选择GDP,点击ok.接下来,添加空间权重矩阵,在弹出来的what windows to open ?全选那几个框,点击OK.
结果如下,
转载自:https://www.zhihu.com/question/25810650