GeoDa莫兰指数、局部莫兰指数的生成及制图

莫兰指数具体分为全局莫兰指数(Global Moran's I)和局部莫兰指数(Local Moran's I)两种,是用来测算空间关系的一种相关系数值。实际研究过程中,通常先计算一个区域内的全局莫兰指数,当全局莫兰指数通过了显著性检验,即认为存在相关性时,紧接着做局部莫兰指数,来具体查看在哪些区域存在集聚效用。Moran's I大于0时,表示数据在区域内呈现出正相关性;Moran's I小于0时,表示数据在区域内呈现出负相关性;Moran's I等于0时,表示数据在区域内呈现随机分布的特性。一般而言,可使用Arcgis、Stata、Python和GeoDa来对区域的莫兰指数进行计算。对比其他软件,GeoDa在计算时更为方便快捷,软件界面为点按式操作,对于小白来说更容易上手。


GeoDa界面

这里本文以山东半岛区域为例子,使用市的GDP数据来展示如何计算莫兰指数。首先,我们需要导入相关底图,并选择合适的方法生成对应的空间邻接矩阵。选择Weights Manager来生成对应空间权重矩阵。


底图导入后界面

GeoDa中支持多种空间权重矩阵的生成,可以分为连通性(contiguity)和距离(distance)两种大的权重分类方式。连通性权重具体可以分为Queen和Rook两种,这两种均为0-1权重矩阵。Rook要求两个区域必须有公共邻边才算相邻,Queen除了边相邻外,只存在一点相邻也会被判定为邻接区域,判定方式更为宽松。Order of contiguity表示的是相邻阶数,假设阶数为k时,表示两个区域之间有k个相邻时,两区域也可以算作相邻区域。


注(Select ID Variable处必须选择无重复的ID才可以)

距离选择为欧式距离(euclidean distance),即表示为不考虑现实环境、地面因素的直线距离。Distance Band表示为一个范围内的被认为是邻接区域,以外的则认为不邻接。K阶近邻则表示为距离A区域最近的K个地区作为A的邻接区域。自适应核距离矩阵(adaptive kernel)则是基于自适应的关联区域损失函数来确定权重。

选择适合你研究的方法生成空间权重后,可对应生产相应的邻接直方图和邻接地图来查看所选择的空间权重是否合适。一般而言,可以同时选择连通性和距离两种权重矩阵相互验证空间自相关的结果。


邻接区域数量直方图

基于所生成的权重矩阵,可以在菜单栏Space-univariate Moran's I中计算单变量全局莫兰指数,选择相对应的变量和权重矩阵后可得下图。即对应的全局莫兰指数值。


进一步,为了获取相对应的P值和z值,我们可以右键单击该图,选择Randomization-999(对应999次随机试验次数,也可以自行选择其他)后可以得到下图的p值和z值。这里对应的p值为0.408,表示未通过显著性检验,不存在全局自相关。

为了继续演示局部莫兰指数的计算,我们先忽略未通过显著性检验的全局莫兰指数。同样,在菜单栏Space-univariate Local Moran's I中计算局部莫兰指数。勾选Sig map和Cluster map后选择ok。可以发现,有三个地区通过了局部莫兰指数的显著性检验。其中,淄博,泰安表现为高高聚集,聊城表现为低-高聚集。

左侧为显著性地图,右侧为集聚地图

最后,我们可以右键单击该图将结果保存到对应的shp文件的数据表中,然后选择在Arcgis中继续完善局部莫兰指数的绘制。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 212,383评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,522评论 3 385
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 157,852评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,621评论 1 284
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 65,741评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 49,929评论 1 290
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,076评论 3 410
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,803评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,265评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,582评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,716评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,395评论 4 333
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,039评论 3 316
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,798评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,027评论 1 266
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,488评论 2 361
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,612评论 2 350

推荐阅读更多精彩内容