R在 map 上绘制网格值

Question

我有几个格式如下的文件:lat,lon,value。我想为每个纬度、经度点绘制一个具有此类值的彩色网格图，并将其重叠在欧盟 map 上。

感谢支持。

Answer 1

大致有两种选择:一种使用lattice和sp包，另一种使用ggplot2包:

格子/sp

首先，您必须将数据转换为 SpatialPixelsDataFrame(由 sp-package 提供)，假设您的数据名为 df:

require(sp)
gridded(df) = c("lat","lon")

然后绘图:

spplot(df, "value")

可以使用 sp.layout 参数覆盖国家边界等附加信息:

spplot(df, "value", sp.layout = list("sp.lines", cntry_boundaries))

其中 cntry_boundaries 是 SpatialLines(DataFrame) 或 SpatialPolygons(DataFrame)。可以使用 rgdal 包或 maptools 包 (readShapeLines) 将多边形集读入 R。

ggplot2

我个人更喜欢使用 ggplot2 而不是 spplot。 ggplot2 更灵活，语法更清晰。请注意，ggplot2 适用于普通的 data.frame，而不适用于 sp-package 的空间对象。

一个最小的例子看起来像这样:

ggplot(aes(x = lat, y = lon), data = df) + geom_tile(aes(fill = value)) + 
    geom_path(data = cntry_boundaries)

有关更多信息，请参阅我的这些早期答案:

• Plotting interpolated data on map
• Average values of a point dataset to a grid dataset