【转】Geary's C
Geary's C is a measure of spatial autocorrelation or an attempt to determine if adjacent observations of the same phenomenon are correlated. Geary's C is also known as Geary's contiguity ratio or simply Geary's ratio. While Moran's I measures global spatial autocorrelation, Geary's C is more sensitive to local spatial autocorrelation. Geary's C is related to Moran's I, but it is not identical.
Geary's C is defined as
where xi denotes the observed value at location i, xbar is the mean of the variable x over the n locations and wij are the elements of the spatial weights matrix, defined as 1 if location i is contiguous to location j and 0 otherwise.
Given the null hypothesis is one of no global spatial autocorrelation, the expected value of Geary's C equals 1. As with Moran's I, inference is based on z-scores:
The p-value for the null hypothesis is given as
The value of Geary's C lies between 0 and 2. 1 means no spatial autocorrelation. Values lower than 1 demonstrate increasing positive spatial autocorrelation, whilst values higher than 1 illustrate increasing negative spatial autocorrelation.
Example
以四张NDVI数据为例,代码测试结果如Fig. 1,与SpaceStat返回结果一致。
Fig. 1
References
[1] Moran's I, fromWikipedia, the free encyclopedia.
[2] Geary's C, from Wikipedia, the free encyclopedia.
[3] Global Spatial Autocorrelation Indices - Moran's I, Geary's C and the General Cross-ProductStatistic. (under the normality assumption the variance, the variance under Randomization方差公式存在两种分布假设:正态分布和随机分布)
[4] BioMedware SpaceStat Help - Geary's C Tutorial.
[5] Geary's C – Springer.
[6] Local Spatial Statistics - ENVI.
[7] Geary's C - AURIN.
[8] Global spatial autocorrelation.
相关文章
- 安卓 android studio 报错 WARNING: Configuration 'compile' is obsolete and has been replaced with 'implementation' and
- SVN更新报错问题(Please execute the 'Cleanup' command)
- 《征服 C 指针》摘录1:什么是空指针?区分 NULL、0 和 '0'
- navcat导入mysql.sql出现:2006, 'MySQL server has gone away'
- manager.add_command('runserver', Server('0.0.0.0', port='8093', extra_files=extra_files))
- [Ramda] Handle Branching Logic with Ramda's Conditional Functions
- Error:A problem occurred configuring project ':app'.
- UnicodeDecodeError: 'utf8' codec can't decode byte 0xce in position 47: invalid continuation byte
- ERROR 1044 (42000): Access denied for user ''@'localhost' to database 'mysql'
- 使用 turtle报错 name 'setup' is not defined
- Universal-image-loader Mason 修复版(ImageLoader Image can't be decoded)
- [Ramda] Convert a Promise.all Result to an Object with Ramda's zip and zipObj
- 关于define('DISCUZ_ROOT', substr(dirname(__FILE__), 0, -7));的理解
- [转]Class 'ThinkLog' not found
- Atitit. http 代理原理 atiHttpProxy 大木马
- Atitit。Cas机制 软件开发 编程语言 无锁机制 java c# php
- Atitit.木马病毒强制强行关闭360 360tray.exe的方法
- Atitit.字节数组转字符串 base64 base16 Quoted-printable 编码原理设计 attilax 总结
- 笔记:Why don't you pull up a chair and give this lifestyle a try?
- cannot bind non-const lvalue reference of type 'std::__cxx11::string&
- ▶ 0001 No application 'E:wwwgoblog' found in your GOPATH
- python中为什么需要使用“if __name__ == '__main__'”语句