我需要可视化一个ET映射,其中包含来自另一个ED映射上相同的中断单元格/位置的中断值。
这是我到目前为止所拥有的以及预期的地图。
library(tmap)
library(sp)
library(rgdal)
library(raster)
ED1 <- tm_shape(ED)+
tm_raster(breaks=c(0,500,1000,2000,3000,
4000,5000,cellStats(ED,"max",na.rm=TRUE)),
style = "fixed")+
tm_legend(outside = TRUE, outside.position = c("right"))+
tm_layout(legend.title.size = 0.5,
legend.text.size = 0.5,
legend.bg.color = "white")
ED1
ET1 <- tm_shape(ET)+
tm_raster()+ # here I need the trick. I need to visualise the values in the cells belongin to the breaks of ED1
tm_legend(outside = TRUE, outside.position = c("right"))+
tm_layout(legend.title.size = 0.5,
legend.text.size = 0.5,
legend.bg.color = "white")
ET1
tmap_arrange(ED1, ET1)
到目前为止,ET中的中断是默认情况下的中断,而ED中的中断是由我自定义的。但是,我希望ET的中断处的值对应于属于ED中断值的单元格。
Here是可复制的示例。