在R中的levelplot上覆盖大陆轮廓

我试图在图层上重叠大陆轮廓，并希望用经度和纬度值标记x轴和y轴。目前，我没有得到一个错误，但没有轮廓出现。我想知道这是否是因为我没有正确定义lat和lon轴

library(rgdal)
library(rasterVis)
library(RColorBrewer)
library(ncdf4)
library(maptools)

f<-nc_open("cases5.nc")
cases<-ncvar_get(f,"cholera_cases")

sen<-getData('GADM', country='Senegal', level=1)

colr <- colorRampPalette(brewer.pal(11, 'RdYlBu'))

for (i in seq(1,1)){

name<-sprintf("plot_%03d.png",i)

png(name)
print(name)

p<-levelplot(cases, 
      margin=FALSE,                       # suppress marginal graphics
      colorkey=list(
        space='bottom',                   # plot legend at bottom
        labels=list(at=0:10, font=4)      # legend ticks and labels 
      ),    
      par.settings=list(
        axis.line=list(col='transparent') # suppress axes and legend outline
      ),
      scales=list(draw=FALSE),            # suppress axis labels
      col.regions=colr,                   # colour ramp
      at=seq(0, 10, len=101))
  +layer(sp.polygons(sen, fill="transparent"))

print(p)

dev.off()
}

库（rgdal）
图书馆（拉斯特维斯）
图书馆（RColorBrewer）
图书馆（ncdf4）
图书馆（地图工具）
f我认为您需要将一个
光栅层
对象传递到levelplot：

library(rgdal)
library(rasterVis)
library(RColorBrewer)
library(ncdf4)
library(maptools)
library(raster)

r <- raster("/Users/christopher/Downloads/cases5.nc")
colr <- colorRampPalette(brewer.pal(11, 'RdYlBu'))

p <- levelplot(
  r,
  margin = FALSE,
  colorkey = list(space = 'bottom',
                  labels = list(at = 0:10, font = 4)),
  par.settings = list(axis.line = list(col = 'black')),
  scales = list(draw = FALSE),
  col.regions = colr,
  at = seq(0, 10, len = 101)
)

sen <- getData('GADM', country = 'Senegal', level = 1)
p + layer(sp.polygons(sen))

库（rgdal）
图书馆（拉斯特维斯）
图书馆（RColorBrewer）
图书馆（ncdf4）
图书馆（地图工具）
图书馆（光栅）
R