r - 仅根据邮政编码在 R 中绘制热图

Question

我想在 R 中绘制一个热图，但我的数据文件是这样的

Lat, Long, Zip, Zvalue

我基本上需要在 lat 和 long 值之间进行插值，并根据 zvalue 绘制颜色

我怎么能在 R 中做到这一点

我最终想要得到这样的东西

Answer 1

包裹spatstat是你的 friend !

http://www.spatstat.org/spatstat/

如果您的数据是一组事件(例如许多 (x,y) 的“犯罪发生在 (x,y)”)，那么您可以使用核密度估计来生成热图。你可以在这里看到一个例子:

https://github.com/drewconway/WikiLeaks_Analysis/blob/master/density.r

特别是第 72 行。

如果您的 z 值是真实值(阅读您的问题可能就是这种情况)，那么您可以使用 spatstat 的 smooth函数使用高斯核来执行插值并取回一组通过插值数据生成的像素。

或者，您可以使用 akima包(由 spatstat 作者推荐)用于在您指定的位置进行插值。它使用线性或基于样条的插值，看起来非常简单(尽管我没有经验!)。

更一般地说，您尝试执行的操作通常称为“克里金法”，因此如果您为此使用谷歌搜索，您会得到很多结果。见例如gstat和 geoR包。

最后(和 FTW)你可以使用高斯过程来做同样的事情。根据您的数据，这实际上会给您一个可能的插值分布。 R包kernlab有一个实现，但我不知道如何使用它。