空间自相关分析在地学中的应用

  • 格式:pdf
  • 大小:2.74 MB
  • 文档页数:50
2007年科技成果交流
空间自相关分析在地学中的应用
刘永生
汇报内容

空间自相关概念

应用实例 结论

“任何事物间都存在着关系,但距离越近关 系越大” ——Tobler’s First Law (TFL) (Tobler 1970)。 地理学第一定律

时间上和空间上的相关性是自然界存在秩序、 格局和多样性的根本原因之一 (Goodchild, 1996)
在对地学数据的处理研究中,必须面对的一个很重
要的问题就是数据之间广泛存在空间自相关。由于空 间信息的自相关特性,经典的线性回归模型在此无能 为力。 在空间分析中,空间自相关分析是一个非常重要的 领域。
计算空间自相关的方法有很多种。用的最广的主要 有:Moran’s I、Geary’C、Getis、Join count等等,但 这些方法作用不是很相同。一般来说,这些方法从功 用上来说大致可以分为两大类:一为全域型,另一为 区域型。全域型方法是探索整个空间上属性数据的分 布模式以及空间结构情况,而局域型空间分析方法可 以进一步探索属性数据在空间上的聚集行为、聚集位
0
2
4
6
1
1 1
1
1 1
1
1 1
1
1 1
1
1 1
1
1 1
1
1 1
1
1 1
1
1 1
1
1 1
1
1
1
1
1
1
1
1
1
100
1
100
1
1
1
1
1
1
1
1
1
1 1 1
1
1 1 1
1
1 1 1
1
1 1 1
100
1 1 1
100
1 1 1
1
1 1 1
1
1 1 1
1
1 1 1
1
1 1 1
1
1
1
1
1
1
1
1
1
1
0.0
2 4 6 8 距离 10 12 14
实例二:
在元素空间强的点分布 (Al2O3)图 (297,1公里尺度)(浙江 上虞数据)
P值<0.05
空间自相关性强的点 分布(Al2O3)图 (67, 1公里尺度) (浙江上虞数据)
P值<0.05
空间自相关性强的点分布 (Al2O3)图 (379, 2公里尺度) (浙江上虞数据)
P值<0.05
空间自相关性强的点分 布(Al2O3)图 (260, 2公里尺度) (浙江上虞数据)
P值<0.05
1公里尺度
P值<0.05
Al 67
B
Ca
Ce 117
MgO 515
N 293
Na2O OrgC 441 261
S
Se
SiO2 K2O La 61 80
409 475
192 254 196
以Rook方式为例:
a a b c d
9 8 7 6 5 4 3 2 1 0 1 0 1 0 0 0
a d
b e
c f
b
1 0 1 0 1 0 0
c
0 1 0 0 0 1 0
d
1 0 0 0 1 0 1
e
0 1 0 1 0 1 0
f
0 0 1 0 1 0 0
g
0 0 0 1 0 0 0
h
0 0 0 0 1 0 1
n
n
j
x)
ij
(x x )
i 1 i
n
2
1 100 1
100 1 100
1 100 1
100 1 100
1 100 1
100 1 100
1 100 1
100 1 100
1 100 1
100 1 100
100
1
1
100
100
1
1
100
100
1
1
100
100
1
1
100
100
1
1
100
100
1 100 1 100
1
100 1 100 1
100
1 100 1 100
1
100 1 100 1
100
1 100 1 100
1
100 1 100 1
100
1 100 1 100
1
100 1 100 1
100
1 100 1 100
1
100 1 100 1
-1.0
2 4 6 8 距离 10 12 14
Moran'I概 率 密 度 曲 线
1.5 Density 0.0
-1.0
0.5
1.0
-0.5
0.0
0.5
N = 10000 Bandwidth = 0.0353
2.5%
2.5%
-1.96
0
1.96
显著性检验
I (d )
n
W
i 1 j i
n
n

W ( x x )( x
i 1 j i ij i
-0.8
-0.6
moran'I -0.4 -0.2 0.0 0.2
Z(I) -10
2 4 6 8 距离 10 12 14
-5
0
5
1
1 1
1
1 1
1
1 1
1
1 1
1
1 1
100
100 100
100
100 100
100
100 100
100
100 100
100
100 100
1
1
1
1
1
1
1
1
1
1
100
100
1
1 1 1 1
1
100 100 1 1
1
100 100 1 1
1
1 1 1 1
1
1 1 1 1
1
100 100 1 1
1
100 100 1 1
1
1 1 1 1
1
1 1 1 1
0.0
2 4 6 8 距离 10 12 14
0.1
moran'I 0.2 0.3
0.4
Z(I) -2
2 4 6 8 距离 10 12 14
i
0 0 0 0 0 1 0
g
h
i
e f g
h
i
0
0
0
0
0
0
0
0
1
0
0
1
1
0
0
1
1
0
MoranI=0.5
统计意义?
在开源自由统计软件R 2.5.1环境下进行编程处理。 R 是一种为统计计算和图形显示而设计的语言及环境, 它提供了一系列统计和图形显示工具。空间自相关利 用软件包spdep-0.4-6进行计算。
(1040,lag=1km)
(1040,lag=1.5km)
(1040,lag=2km)
(1040,lag=3km)
I (d )
n
Wij
i 1 j i
n
n

W ( x x )( x
i 1 j i ij i 2 ( x x ) i i 1 n
n
n
j
x)
置等空间信息。在许多研究中,Moran’I和Getis是最被
经常使用的方法。
空间统计参数 (The General Cross-Product Statistic)
Wij Cij
i j
(Hubert et al. 1981; Upton and Fingleton, 1985)
空间权重(Wij)构造方式
重金属 地质统计法
As
Cd
Cr
Cu
Hg
Pb
Ni
Zn
15.2 7.7 37.4 13.0 24.8 26.7 23.6 17.4
空间分析法
15
8
24
12
16
26
24
13
(1)从两种途径得到的空间变程来看,空间自相关分
析法得到的空间变程最为稳定,也最为直观,而地
质统计法得到的空间变程随着拟合模型的不同其结
0.1
moran'I 0.2 0.3 0.4
Z(I) -2
2 4 6 8 距离 10 12 14
0
2
4
6
应用实例
实例一:
空间最大变异距离研究
刘永生,李瑞敏.不同方法下土壤重金属元素空间最大变异 距离研究.岩矿测试,2007,26(4)
重金属元素的半变异函数图
土壤重金属元素半变异函数参数
重金属元素 理论模型 块金值 As Cd Cr Cu Hg 指数模型 指数模型 指数模型 指数模型 指数模型 0.080 0.072 0.050 0.086 0.220 基台值 0.280 0.126 0.131 0.182 0.907 变程/km 15.188 7.714 37.374 13.021 24.760 块金系数 拟合度 0.285 0.569 0.379 0.471 0.242 0.896 0.763 0.820 0.894 0.910
100
100
100
100
100
100
100
100
1
1 1 1