r - 使用 R 中“rasterVis”包中的 levelplot 函数绘制只有 2 个像素的栅格

Question

我正在尝试使用“rasterVis”包中的 levelplot 函数绘制以下栅格，但这给了我一条直线。我尝试使用“at”命令对图例进行不同的中断，但无济于事。我在文档中尝试了我能做的任何事情，但找不到答案。此外，使用 R 中的“plot()”和“image()”命令可以很好地绘制栅格。请帮助。

library(raster)
library(rasterVis)
sr <- "+proj=utm +zone=12 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0"
r <- raster(resolution =c(24600, 28900), crs = sr, nrows = 1, ncols = 2, ext = extent(c(573800, 623000, 3508713, 3537613)))
r[] <- c(1,2)
levelplot(r, margin = FALSE)    

levelplot 的输出图

编辑 1 当我定义没有范围和投影的栅格时，栅格用 levelplot 绘制得很好。我找不到原因。

r2 <- raster(nrows =1, ncols =2)
r2[] <- c(1, 2)
levelplot(r2, margin = FALSE)

栅格图 r2

Answer 1

如果resolution使用，参数ncols和nrows被忽略。删除resolution并nrow用 绘制它levelplot。

library(raster)
library(rasterVis)
sr <- "+proj=utm +zone=12 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0"
r <- raster(ncol = 2, crs = sr, ext = extent(c(573800, 623000, 3508713, 3537613)))
r[] <- c(1, 2)
levelplot(r, margin = FALSE)