回归和空间聚类分析
- 格式:pdf
- 大小:2.80 MB
- 文档页数:30
CASE 7-2
Case 7-2 用空间聚类法分析 中国南部地名
首都师范大学,2008年12月8日-16日
Capital Normal University, 8-16 December 2008
• 这个实验研究的是中国南部泰语地名的分布状态, 是正在进行中的“研究中国南部泰语地名历史起源” 大型项目的一部分。少数民族的汉化在中国的历史 中是一个进行了长久并且正在进行的历史过程。随 着时间的变迁,这种历史进程也在地理地名上留下 了标记。由于许多旧的泰语地名是根据泰语语系中 地物特征命名的,因此他们仍然可以被识别出来。 而与此同时,许多曾经的泰语地名也在汉化的过程 中被删改了。这个项目的目的就是重建所有早期泰 语地名的分布,以期能发现在汉族来到南方之前, 泰民族在中国南部的最初居住范围。
2008-12-15
CASE 7-2
17
右键点击
2008-12-15
CASE 7-2
18
9
Mon. 15 December 2008
CASE 7-2
4. 将qztai图层属性表导出为 dBase文件 qztai.dbf
• • • • (1)右键点击属性表下方的Options (2)选择Export (3)在弹出的窗口中,设置导出文件名 称为qztai.dbf (4)点击OK
2008-12-15
CASE 7-2
11
Right_click
2008-12-15
CASE 7-2
12
6
Mon. 15 December 2008
CASE 7-2
3. 准备坐标文件
• (1)选择ArcToolBox • (2)选择Data Management Tools -> Features -> Feature to Point • (3)在弹出窗口中设置,输入文件为qztai 下的point文件,将输入文件命名为qztai • (4)点击OK
2008-12-15
CASE 7-2
3
步骤一: 用ArcGIS准备SaTScan软件需要的数据
在SaTScan软件平台下,用伯努利(Bernoulli) 模型执行基于点的空间聚类分析需要定义三个数据 文件: •事件文件:包括区位ID和每个区位的事件数 •非事件文件:包括区位ID和每个区位的非事件数 •坐标文件:包括区位ID 和对应的笛卡尔坐标或经纬 度坐标
2008-12-15
CASE 7-2
19
右键点击
2008-12-15
CASE 7-2
20
10
Mon. 15 December 2008
CASE 7-2
步骤二: 用SaTScan执行空间聚类分析
• 1. 运行SaTScan软件,选择Create New Session
2008-12-15
CASE 7-2
2008-12-15
CASE 7-2
46
23
Mon. 15 December 2008
CASE 7-2
2008-12-15
CASE 7-2
47
2008-12-15
CASE 7-2
48
24
Mon. 15 December 2008
2008-12-15 CASE 7-2 2
1
Mon. 15 December 2008
CASE 7-2
• 实验所用数据有: 1.钦州市所有乡镇的点数据 qztai (coverage), 其中的TAI属性标识出此地名是否为泰语地 名,是泰语地名(=1),不是泰语地名(=0)。 2.钦州市六个县的面数据 qzcnty (shapefile)。 • 实验所需软件: 1. ArcGIS; 2. SaTScan。
2008-12-15 CASE 7-2 4
2
Mon. 15 December 2008
CASE 7-2
1. 将e00文件转为coverage文件 (准备事件文件)
(1)打开ArcCatalog -> 选择菜单栏中的View -> Toolbars -> ArcView 8x Tools -> Import from Interchange File (2)输入文件(Input File): …\qztai.e00 输出文件(Output File): …\qztai qztai文件属性表中的变量TAI已经定义了每个 区位的事件,其值为1正好表示该区位事件数为 1,为0表示该区位事件数为0。
2008-12-15
CASE 7-2
15
2008-12-15
CASE 7-2
16
8
Mon. 15 December 2008
CASE 7-2
• (9)右键点击qztai.shp文件,选择Open Attribute Table • (10)在属性表中右键单击POINT_X字段 • (11)在弹出窗口中输入计算: POINT_X = POINT_X / 100000 • (12)右键单击POINT_Y字段,输入的计算式: POINT_Y = POINT_Y / 100000
2008-12-15
CASE 7-2
30
15
Mon. 15 December 2008
CASE 7-2
2008-12-15
CASE 7-2
31
• (8)在Output标签下,输出文件名( Result File)定为Taicluster,勾选dBase 下的五个复选框,确定五项均为输出文件 的内容。
2008-12-15
CASE 7-2
44
22
Mon. 15 December 2008
CASE 7-2
2008-12-15
CASE 7-2
45
标识集聚圈范围
• (4)在Drawing工具条里选择New Circle -> 将一级集聚区域用圆标识出来 -> 双击圆,修改其颜色等属性 -> 用同样方法将次级集聚圈标识出来
2008-12-15
CASE 7-2
25
注意: 此时,定义的是非事件文件, 因此非事件数选择的是字段 NONTAI。
2008-12-15
CASE 7-2
26
13
Mon. 15 December 2008
CASE 7-2
• (4)点击 Coordinates File 之后的 按钮,选 择qztai.dbf作为输入文件 • (5)在SaTScan Import Wizard对话框中设置 如下: Location ID: QZTAI_ID Latitude (y-axis): POINT_Y Longitude (x-axis): POINT_x • (6)点击Next后选择输出文件路径(选择默认 即可)
2008-12-15
CASE 7-2
9
2008-12-15
CASE 7-2
10
5
Mon. 15 December 2008
CASE 7-2
• (5) 在属性表中右键点击NONTAI字段 • (6) 选择Field Calculator • (7) 在弹出的对话框中输入如下: NONTAI = 1-TAI • (8) 点击OK
2008-12-15
CASE 7-2
13
qztai
2008-12-15
CASE 7-2
14
7
Mon. 15 December 2008
CASE 7-2
• (5)选择ArcToolBox下的Index标签 • (6)输入关键字Add XY Coordinates, 双击弹出对话框 • (7)选择qztai.shp文件作为输入文件 • (8)点击OK
2008-12-15
CASE 7-2
27
ቤተ መጻሕፍቲ ባይዱ
2008-12-15
CASE 7-2
28
14
Mon. 15 December 2008
CASE 7-2
如图,即设 置好了进行 伯努利空间 聚类分析的 三个文件。
2008-12-15
CASE 7-2
29
• (7)点击Analysis标签,设置如下: 分析类别:纯空间( Purely Spatial ) 概率模型:伯努利(Bernoulli) 区域扫描速度:高速(High Rates)
2008-12-15 CASE 7-2 34
17
Mon. 15 December 2008
CASE 7-2
2008-12-15
CASE 7-2
35
步骤三: 分析结果的制图
2008-12-15
CASE 7-2
36
18
Mon. 15 December 2008
CASE 7-2
• (1)在ArcMap中添加图层qztai.shp, qzcnty.shp和Taicluster.gis • (2)右键点击qztai.shp,选择join, 将Taicluster.gis与qztai.shp连接,公 共字段为QZTAI_ID (qztai.shp)和 LOC_ID (Taicluster.gis)。
2008-12-15 CASE 7-2 5
2008-12-15
CASE 7-2
6
3
Mon. 15 December 2008
CASE 7-2
2. 准备控制文件
• (1) 打开ArcMap • (2) 打开菜单栏中的File -> Add Data -> 选择qztai -> 点击Add添加图层
2008-12-15
21
2.使用Import Wizard来定义事件文件
• (1)在第一个标签(Input)下,点击 Case File 之 后的 按钮 • (2)选择qztai.dbf作为输入文件 • (3)在弹出的SaTScan Import Wizard对话框 中设置如下: Display SaTScan Variables for: Bernoulli Location ID: QZTAI_ID Number of Cases: TAI • (4)点击Next后选择输出文件路径