将矢量世界地图添加到 R 中的光栅动画

我用以下内容创建了一个动画animate的功能raster包并想在其上添加世界矢量地图。

这是我的代码，在对其进行动画处理之前，我根据月平均温度的 3D 数组创建了一个栅格块：

r <- brick(ncols=nb.cols, nrows=nb.rows)
r <- setValues(r, monthly.mean)
animate(r)

光栅动画效果很好，但如何添加世界矢量地图？

非常感谢

谢谢你！这两种方法都很有效！

有兴趣的人可以快速回答：

使用 saveHTML：

require(animation)
require(rgdal)
require(raster)

countries <- readOGR(dsn=".", layer="Countries")  # Shapefile of world countries

r <- brick(ncols=144, nrows=72)
r <- setValues(r, monthly.mean)          # monthly.mean is a 3D array 144x72x420

saveHTML({
  for(m in 1:420) {                      # serie of 420 months
    plot(r[[m]], zlim=c(min, max))       # min, max of the legend
    plot(countries, add=T)
  }
})

充满乐趣：

require(rgdal)
require(raster)

countries <- readOGR(dsn=".", layer="Countries")

r <- brick(ncols=144, nrows=72)
r <- setValues(r, monthly.mean)

fun <- function() {
  plot(countries, add=T)
}

animate(r, zlim=c(min, max), addfun=fun)