Мне нужно визуализировать карту 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.
Здесь - воспроизводимый пример.