我在研究区域内有几个特定物种的丰度图。我一直在成功使用图像绘图。但是现在,我只想绘制存在/不存在的图形,因此值从0到1,并且该图形仍然显示了一个连续的图例。有办法修改吗?
我尝试过使用图像,但图例并不那么酷。 这是我的代码:
colnames(zz)<-c("x","y","z")
ras<-rasterFromXYZ(zz)
asc<-asc.from.raster(ras)
image.plot(asc,
xlim=c(-9.5,-1),ylim=c(43,48),zlim=c(min(asc,na.rm=T),max(asc,na.rm=T)),
ylab="Latitude (ºN)",xlab="Longitude (ºW)",
col = viridis(2, option = "D"))