在R中创建地图上的颜色网格

时间:2012-07-26 09:22:29

标签: r spatial-interpolation

我一直试图以如下所示的方式在我的世界地图上创建一个网格 example of gridded map

我希望能够根据文件中的rate2列为网格着色。我所获得的大多数解决方案只根据国家/地区为世界着色。请问我该如何解决这个问题。我附上了我的csv。

heatmap.csv

1 个答案:

答案 0 :(得分:2)

看一下R中的一些空间插值包,Spatial Task View for R提供了一个很好的列表,例如: automap,gstat,Fields,geoR等。要获得地图,您需要:

  • 定义预测网格
  • 如果需要,请拟合数据的某些数学属性,例如:在克里金法的情况下的变异函数。
  • 将值预测到此新网格上。

我认为对于这些数据,您的第一个最佳选择是一些反距离加权或基于样条的插值(例如,字段中的Tps)。

一些注意事项:

  • 您的数据仅在陆地上,我建议仅在陆地上定义预测网格。
  • 您在同一空间位置有多个观测值。我确信在automap / gstat中的克里金实现不会喜欢这个。扔掉,或聚合(例如平均)重复的观察。有一些方法可以使用重复数据来估计测量不确定度,并在插值过程中考虑到这一点,但这并非易于实现。

有关更多信息,请参阅我之前关于该主题的两篇帖子: