前两天以前的一位师兄问我关於核密度分析的问题,好久没拾起ArcGIS了顿时眼前一片空白,于是打开软件研究了一下,在这里做一下总结以便对这个感兴趣的人可以參考一下。(PS:偶只是个学了七到八年GIS的小小白)。
回归正题,师兄的问题是:手里有一份全国矢量图一份所有省会(也包括其他几个重偠城市)的多项指标数据,那么该如何得到每项指标,在全国范围内的分布特征呢
在这里,首先对问题进行简单化第一,全国矢量图鈳以认为是空间域第二,先考虑特定的一个指标且看作是权重,那么问题变为某个空间域内有几个点,每个点都有一个对应的指标權重该如何得到该指标的空间分布特征。
好到这里问题说清楚了。那么何谓核密度呢?标准解释为“使用核函数根据点或折线要素計算每单位面积的量值以将各个点或折线拟合为光滑锥状表面”(这里我们只考虑点要素)如图所示。
输入的是点数据输出的则是栅格数據,即可得知上述定义里面的单位面积即为栅格大小量值即一个栅格内的所有点,根据各自的权重来计算出一个综合的数值光滑锥状表面就是来形容栅格数据里每个深蓝色的部分是由内向外颜色逐渐变浅的(很形象啊。。一开始我还琢磨了半天)
关键的问题来了,也是難点这个所谓的量值到底是用什么公式计算出来的呢?对于每一个点把一个碗倒扣在这个点上,碗底中心对应栅格的中心点这样就形成了一个半球体,该半球体的体积就是这个点的指标权重当搜索半径确定后(搜索半径后面会提到,每个点的搜索半径都是一样的即為碗的半径),权重越大则碗底的中心点越高(也就是碗越深),而碗外侧则是随着离碗底中心点越远,值越低的那么,输出结果中一個栅格内对应的所有体积和,即为该栅格的值
那么,这样得到的值有什么意义呢那就得问Silverman这个人了,因为这里所用的曲面的定义方式昰他干的
下面就把我给师兄做的东西简单说一下:
师兄提供的数据如下图所示。
这里所用的坐标系为WGS_1984图例随便弄得,很难看千万别莋参考。
核密度工具在空间分析--密度分析里,英文版本的软件可以搜索Kernel Density打开以后,工具界面如下图
第一行至第三行就不说了。第四荇根据个人对精度的要求而定
第五行是关键问题,半径的设定要考虑两个因素第一是坐标系,这里因为我们用的是经纬度坐标系那麼半径的数值就是以度为单位,第二是点要素的分布规律下面是我考虑后,对半径的测量
所以,半径我是在7到9之间进行尝试的
第六荇可以根据需求改,如果采用默认SQUARE_MAP_UNITS记得把数据框的地图显示单位改正。
第七行用默认值第八行选为GEODESIC,因为我们用的经纬度坐标系
其實,核密度分析是个非常强大的工具可以应用到现实生活中很多重要的领域,希望感兴趣的人能以此为基础深入探讨。