前两天,以前的一位师兄问我关于核密度分析的问题,好久没拾起ArcGIS了,顿时眼前一片空白,于是打开软件,研究了一下,在这里做一下总结,以便对这个感兴趣的人可以参考一下。(PS:偶只是个学了七到八年GIS的小小白)。。。
回归正题,师兄的问题是:手里有一份全国矢量图,一份所有省会(也包括其他几个重要城市)的多项指标数据,那么,该如何得到每项指标,在全国范围内的分布特征呢?
在这里,首先对问题进行简单化,第一,全国矢量图可以认为是空间域,第二,先考虑特定的一个指标,且看作是权重,那么,问题变为某个空间域内有几个点,每个点都有一个对应的指标权重,该如何得到该指标的空间分布特征。
好,到这里问题说清楚了。那么,何谓核密度呢?标准解释为“使用核函数根据点或折线要素计算每单位面积的量值以将各个点或折线拟合为光滑锥状表面。”(这里我们只考虑点要素)如图所示。
输入的是点数据,输出的则是栅格数据,即可得知上述定义里面的单位面积即为栅格大小,量值即一个栅格内的所有点,根据各自的权重来计算出一个综合的数值,光滑锥状表面就是来形容栅格数据里每个深蓝色的部分是由内向外颜色逐渐变浅的(很形象啊。。。一开始我还琢磨了半天)。
关键的问题来了,也是难点,这个所谓的量值到底是用什么公式计算出来的呢?对于每一个点,把一个碗倒扣在这个点上,碗底中心对应栅格的中心点,这样就形成了一个半球体,该半球体的体积