我试图在图层图上重叠大陆轮廓,并想用经度和纬度值标记 x 和 y 轴。目前我没有收到错误,但没有出现轮廓。我想知道这是否是因为我没有正确定义纬度和经度轴?
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()
}
此处提供了一个示例 netcdf 数据文件:
http://clima-dods.ictp.it/Users/tompkins/stackoverflow/cases5.nc