R Tmap пакет, присваивающий разрывы, соответствующие ячейкам, принадлежащим другим значениям разрывов, 2 растра - PullRequest
0 голосов
/ 30 сентября 2018

Мне нужно визуализировать карту 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.

tmap

Здесь - воспроизводимый пример.

...