r - 使用多个属性为shapefile着色

library(raster)
admin <- getData('GADM', country='FRA', level=2)

set.seed(123)

id <- data.frame(ID_2 = admin@data$ID_2, day1 = sample(1:20,96,replace = T),
                                   day2 = sample(50:80,96,replace = T),
                                   day3 = sample(120:140,96,replace = T),
                                   day4 = sample(200:230,96,replace = T))

 admin.shp <- merge(admin,id)

如果我想使用day1，day2，day3或day4为绘图着色，则可以执行以下操作:

plot(admin.shp, col = admin.shp@data$day1)

我正在尝试做的是生成从第一天到第二百三十天的某种.gif文件或动画:

  range <- c(min(id$day1),max(id$day4))

如果日期与day1中的值匹配，则多边形应变为绿色，
如果在第2天匹配，多边形应变为蓝色，
如果在第3天匹配，多边形应变为橙色，
如果在第4天匹配，则多边形应变为红色

如果我必须对单个列执行此操作(例如对于day1)，则可以执行以下操作:

library(magick)
c(min(id$day1),max(id$day1)) # 1, 20

for(i in 1:20){

   breaks.pl <- c(0, i, 21)
   col.pl <- c("green4","white")
   cuts.pl <- cut(data.frame(admin.shp)[, "day1"],breaks = breaks.pl)

   png(paste0(i,".png"), width = 1000, height = 600)
   plot(admin.shp, col = col.pl[cuts.pl], border = 'transparent', main = paste0("day:",i))
   plot(admin.shp, add = T)
   dev.off()
 }

这将产生一系列的png文件，然后我可以生成这20个png文件的.gif文件，以创建day1的动画

list.png <- list()

for(i in 1:20){
png.file <- image_read(path = paste0(i,".png"))
list.png[[i]] <- png.file
}

png.stack1 <- list.png[[1]]

for(i in 1:20){
    png.stack <- list.png[[i]]
    png.stack1 <- c(png.stack1,png.stack)
 }

 png.img <- image_scale(png.stack1)
 png.ani <- image_animate(png.img, fps = 1, dispose = "previous")

 image_write(png.ani, "my.animation.gif")

但是，我想扩展此范围，以便将所有四列day1都使用day4

最佳答案

根据您的代码改编，这应该可以解决问题:

colors = c("white", "green4", "blue", "orange", "red")

library(animation)
ani.options(interval=.05)

i = 0
saveGIF({
  for(k in 2:5){
    while (i < max(id[,k])) {
      print(i)
      i = i + 1
      breaks.pl <- c(0, i, max(id$day4)+1)
      col.pl <- c(colors[k], colors[k-1])
      cuts.pl <- cut(data.frame(id)[, k], breaks = breaks.pl)

      plot(admin.shp, col = col.pl[cuts.pl], main = paste0("day:",i))
    }
  }
})

演示:

关于r - 使用多个属性为shapefile着色，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/49830708/