我正在尝试创建一个降水量图,其中的降水量超过某个阈值。我在61x101网格上有一个名为tp98slice的数据,其值为0或1。如果每个网格点的值等于1,我想将这些点绘制为点。
我尝试了以下代码,
levelplot(tpslice ~ lon * lat, data=grid, at=cutpts, cuts=50, pretty=T,
col.regions=(brewer.pal(9,"Blues")),
xlab="Longitude", ylab="Lattitude") +
layer(sp.points(tp98_pts, pch=ifelse(tp98slice==1, 1, 4), cex=2))
但是图层位(我的观点)不起作用/未显示。 Current plot