Geoda空间自相关分析—局部Moran‘I指数（运用GeoDa制作Lisa图）

Geoda空间自相关分析—局部Moran’I指数（运用GeoDa制作Lisa图）

首先介绍一下什么是空间自相关（引用：[https://wenku.baidu.com/view/4ab39c93650e52ea54189885.html]）
各位要是觉得我写的不好可以去看看虾神（这是一个大神）的文章。PS：虾神主页：[https://blog.csdn.net/allenlu2008]

根据 Tobler 的地理学第一定律，即地理分析基本定律：“地理事物之间都存在着相应的关系，其中相距较近的事物比相距较远的关系更密切”。空间自相关作为一种空间统计方法，其全局与局部空间自相关都能够较好地描述地理事物间的关系，衡量出事物空间要素属性间的聚合或离散的程度。

全局空间自相关：

** 全局空间自相关用于描述空间要素属性值在全区域内的空间特征，并反映其邻域属性值相似度，对于某一特定空间属性在全区域分布特征分析描述具有重要意义，常用Geary’s C、 Moran I 指数表示。Moran I 指数的值域为[-1,1]之间，小于0负相关，等于0不相关，大于0正相关。
**

局部空间自相关：

**相较于全局空间自相关，局部空间自相关则是计算分析区域内各个空间对象与其邻域对象间的空间相关程度，计算分析空间对象分布中所存在的局部特征差异，反映局部区域内的空间异质性与不稳定性。同样有Geary’s C、 Moran I 指数，但常用Local Moran I衡量，Anselin也被称为Lisa。
1.首先打开GeoDa加载数据

ps：路径不可有中文名称
加载完成
2.新建空间权重矩阵（我理解的是在进行相关性分析之前要建权重矩阵）
PS：这里要选择唯一值ID或者其他的唯一值字段
权重类型根据自己的需要选择
Queen contiguity：后相邻
Queen contiguity：车相邻

创建成功！
3.打开属性表
PS：属性表不可编辑，但可根据需要进行合并和计算

4.进行局部空间自相关分析—以Univariate Local Moran’s I（单变量莫兰指数）为例
因为我进行侵蚀强度分析所以选择侵蚀编号（就是自己要实验的数据）

PS：最好三个都勾上

莫兰指数为0.085通过显著性水平a=0.05的检验，说明区域分布呈正相关和有一定聚集性。
其中莫兰指数散点图的那些点表示数据数量
第一象限为HH
第二象限为LH
第三象限为LL
第四象限为HL

5.统计第一象限数据的面积（没办法，脑子不好使，只有这种办法。望各位看官不吝赐教。）
首先选中第一象限的点

点击“添加选择”

这样就出现新的字段
点击“聚合”
选中“求和”保存为“.csv”格式的文件

然后在删除变量

其他数据依次类推
最后在excel打开文件进行统计。
PS：虽然办法蠢，但有效。
6.导出Lisa图
目前的GeoDa不支持制图输出，所以只能导出图片，然后在PS中进行处理。

在PS中处理的步骤，各位自行去摸索吧。相信各位PS技术一定比我好。