0

以下示例直接来自 rasterVis::levelplot 文档。请注意,当从 RAT 表中选择的属性是时,图例中的名称不会更改。我自己的数据也是这种情况(我必须绘制一个具有不同土地类别(因子)的分类栅格,我希望这些栅格具有相关的颜色)。有人有解决方法吗?

library(raster)
library(rasterVis)

r <- raster(nrow=10, ncol=10)
r[] = 1
r[51:100] = 3
r[3:6, 1:5] = 5
r <- ratify(r)
     
rat <- levels(r)[[1]]
rat$landcover <- c('Pine', 'Oak', 'Meadow')
rat$class <- c('A1', 'B2', 'C3')
levels(r) <- rat
r

levelplot(r, col.regions=c('palegreen', 'midnightblue', 'indianred1'))

## with 'att' you can choose another variable from the RAT
levelplot(r, att=2, col.regions=c('palegreen', 'midnightblue', 'indianred1'))
levelplot(r, att='class', col.regions=c('palegreen', 'midnightblue', 'indianred1'))

4

0 回答 0