白话空间统计二十一：密度分析（六）R语言实现

白话空间统计二十一：密度分析（六）R语言实现

今天主要讲讲核密度分析的R语言实现，原理神马的，看前面的文章。

下面贴代码和代码注释，还是老规矩，需要源代码还数据的，通过邮箱获取。

还是先看看结果（话说R语言默认的颜色渲染我还是觉得挺不错的，想用其他颜色渲染的，后面有空再说）












那么如果把这四张图都放到同一个颜色渲染尺度下面，如下：（可以点开看大图）,可以再次验证昨天的结论：带宽越大，核曲面的高度越低。




下面贴出代码:


里面各个方法的帮助，参考帮助文档。

下面是一些注意事项：

1、上面我用的是csv直接读取的数据，如果本身就是shapefile的point图层，可以直接通过maptools读取。
2、cellsize参数很重要——如果设置得太小，会占用大量的内存——最后甚至可以爆掉你的内内存——比如下面这样，虾神我机器16G内存，就这样爆掉了：




文章来源：http://blog.csdn.net/allenlu2008/article/details/53258240