r - ラスタースタックにポリゴンレイヤーを重ねる qplot

Question

を使用して 7 つの TIFF ファイルを並べてプロットgplotし、それらの上にポリゴン レイヤーをオーバーレイしたいと考えています。を使用してシェープファイルを読み取ってから、またはをst_read使用してプロットしてみました。しかし、うまくいきませんでした。に関するエラーをスローしていました。私は最終的に関数を使用し、それをプロットするために使用しました。オーバーレイされますが、それは私が望むものではありません。ポリゴンの境界をオーバーレイしたいだけです。geom_sfgeom_polygonaesfortifygeom_path

PS再現可能なラスターおよびベクターレイヤーを作成する方法がわかりません。そのため、必要に応じてサンプル ファイルをアップロードできます。

shp = fortify(shapefile("E:/Vect/BLOCKS.shp"))
rfiles = list.files(path = "E:/Rast", pattern = "*.tif", full.names = TRUE)
pfiles = stack(rfiles[c(7:14)])

plot1 = gplot(pfiles1) + 
  geom_tile(aes(fill = value)) +
  geom_path(data=shp,aes(long, lat),colour="red") + 
  facet_wrap(~ variable, ncol = 7) +
  scale_fill_gradientn(colours = magma(30), na.value = "transparent") +
  theme_bw() +
  theme(axis.text = element_blank(), legend.position = "right", legend.direction = "vertical",
        legend.key.height = unit(5, "cm"), legend.key.width = unit(1, "cm"),
        legend.title = element_blank(), legend.text = element_text(size = 30), 
        strip.text = element_text(size = 28, face = "bold")) +
  coord_equal()