在R中绘制rasterstack的主要组成部分

Question

我正在尝试从here复制图1，其被标记为“从每月NDV得出的标准化主成分图像1到8;数据，1986年1月到1988年12月”。

我正在使用principal包中的psych。这是一个例子，我有一个rasterstack，值和每个坐标。假设每个栅格代表一年内的变量。

library(psych) 
library(raster)

r1 <- r2 <- r3 <- r4<- raster(nrow=10, ncol=10)
values(r1) <- runif(ncell(r1))
values(r2) <- runif(ncell(r2))
values(r3) <- runif(ncell(r3))
values(r4) <- runif(ncell(r4))
x <- stack(r1, r2, r3, r4)

extract.value.from.raster.stack<-extract(x,1:ncell(x))
coord<-xyFromCell(x,1:ncell(x))
lon<-coord[1,]
lat<-coord[2,]

pcaan<-t(extract.value.from.raster.stack)
pca3<-principal(pcaan,nfactor=5,rotate="none")

在进行PCA分析之后，我如何获得类似于具有相关性的纸张的类似图？