2

我一直在 RStudio 中运行以下代码来测试 R 中的向量图功能：

library(raster)
library(rasterVis)
library(lattice)

proj <- CRS('+proj=longlat +datum=WGS84')

df <- expand.grid(x=seq(-2, 2, .01), y=seq(-2, 2, .01))
df$z <- with(df, (3*x^2 + y)*exp(-x^2-y^2))
r1 <- rasterFromXYZ(df, crs=proj)
df$z <- with(df, x*exp(-x^2-y^2))
r2 <- rasterFromXYZ(df, crs=proj)
df$z <- with(df, y*exp(-x^2-y^2))
r3 <- rasterFromXYZ(df, crs=proj)
s <- stack(r1, r2, r3)
names(s) <- c('R1', 'R2', 'R3')

vectorplot(r1)

上面的代码给了我以下错误：

subset.default(sa, 1) 中的错误：“子集”必须是逻辑的

该代码在 Mac OS 上的某些先前版本的包（raster、rastervis）上运行良好。您能否建议如何解决错误并使向量图再次起作用？

目前，我正在使用 Windows 10 和以下版本的 R 和 RStudio：

  • R: 4.1.2 (2021-11-01)
  • RStudio：2021.09.1+372“幽灵兰花”发布 (2021-11-08) 适用于 Windows
4

0 回答 0