白话空间统计二十四:地理加权回归(七)ArcGIS的GWR工具扩展参数说明
白话空间统计二十四:地理加权回归(七)ArcGIS的GWR工具扩展参数说明
近期无论是开发者大会,还是个人工作,相当的忙,所以停了一段时间……不过地理加权回归写到第七章,自我感觉也差不多了,无论是基础理论还是来历,包括基础参数的意义,都应该介绍得比较清楚了,当然,后面可能是大家更关心的内容,也就是在ArcGIS(或者其他软件里面),怎么去执行地理加权回归,执行完成之后,生成的结果又如何解读?
下面我们正式进入使用ArcGIS进行地理加权回归分析的步骤。
至于工具的位置和参数说明,请大家看以前的文章:
白话空间统计二十四:地理加权回归(六)ArcGIS的GWR工具参数说明一
当然,前面还缺了一部分,就是GWR的扩展参数部分,作为一只有始有终的虾,准备先把这遗留下来的部分说一说:
从上面的提示上可以看见,整个扩展部分,都是可选,也就是说,有它行,没他也不影响结果。这些参数的作用如下所示:
1、Coefficient raster workspace :系数栅格的工作空间
所谓的工作空间,是ArcGIS工具一个专有术语,在使用过arcpy的同学看来,应该是非常容易理解的的一个概念,设置了工作空间之后,所有的分析结果,默认都会优先写入到这个工作空间里面。
在这里一旦设定了工作空间,那么在回归过程中的截距和各个系数,都会被生成为一个栅格文件,存储在这个工作空间中。
这个栅格有啥用呢?下面可以看看我这里给出来的例子来看看。
数据还是一样,山东省分区县的各类数据:
把五个维度的数据,都做成专题图,进行对比,然后进行GWR,把这五个维度的系数进行栅格化,如下:
从左边矢量的专题图进行对比,其实是很难对数据的差异性进行直观的感受的,但是进行了栅格化之后,这种感受就很直观很明显了。
系数的栅格化,体现出的是该维度数据(各自变量)在不同区域位置对因变量的作用的强弱,很明显的就可以看出在改回归模型中,不同区域的变化强度。
可以通过栅格化,来体现空间的异质性。
第一个财政收入作为因变量,那么没有系数的话,出现的是截距的栅格表面。剩下的自变量,都是各自系数的栅格表面。
根据官方说法,系数栅格化有如下作用:
- 通过统计学上表示较小局部变化的较大全局变量可确定全局策略。
- 通过统计学上表示较强局部变化的较大全局变量可确定局部策略。
- 某些变量可能并不是在全局范围内各区域中均比较显著,因为在某些区域中,它们是正相关的关系,而在其他区域中它们则是负相关的关系。
所谓的局部变化,也就是在某个区域的变化是否剧烈,用地形学的说法,就是等高线(等值线)在一定距离内的线的密度情况,越密,就表示变化越剧烈。所以要找局部变化,最简单的方法就提取等值线,这个不属于我们今天的内容了,有空继续再说。
2、输出单元格大小
这个参数是关于生成栅格的粒度,设置得越小,表示越清晰,平滑,但是所用的计算时间和存储空间越大。
3——5号参数是关于用于填充预测值的。
在空间分析里面,进行未知区域的值预测最常用的方法也就是插值了,但是插值生成栅格之后,再进行指定位置区域的数据提取,中间还有一个过程——而且是插值对面和线要素支持得并不是太好(当然ArcGIS现在已经支持面插值了,但是线还是不行),对多字段数值,支持得也不是太好,特别是多值间有各种关系的,插值都支持不了。这种情况,使用GWR的区域预测就有意义了。
要进行区域预测,首先要准备好一份结构完全一样的数据,从空间参考到需要预测(填充)的字段属性都要一模一样,然后按照前面回归模型的设置,设置好要素类和属性字段,以及输出就可以了。
到此,GWR工具所有的参数设置就讲完了,后面可能还有几个小章节,讲讲结果和应用。
待续未完
转载自:https://blog.csdn.net/allenlu2008/article/details/72861013