在地图上绘制经验正交函数（EOF）

Question

大家好，我正在对以下数据进行EOF分析：https://www.dropbox.com/s/b37gzandqv0kdu9/gpcc_spi6_1901_2013.rda?dl=0

这是我创建的一个文件，其中包含某个区域的SPI，我想对它进行EOF分析，这是我到目前为止所做的：

require(utils)
library(colorRamps)
library(sinkr)

Datos <- load("gpcc_spi6_1901_2013.rda")

dim(spi_c)#18,31,1356
#I cut some of the years off
spi1 <- spi_c[ , ,205:1356] #31, 18, 1152 lon ,lat, months

#Then I only take the months from May-Oct for the whole series
reform <- array(spi1, dim = c(31,18,12,96))
verano <- reform[ , ,5:10, ]

#Come back to lon, lat, month format
r_verano <- array(verano, dim=c(31,18,576))

#I know you can only apply EOF to a matrix so I convert the 3D array
#Into a 2D matrix

arreglo <- sapply(1:dim(r_verano)[3], function(x)t(r_verano[,,x])) #558, 576
arreglo <- t(arreglo)

#EOF
neofs <- 20
#I wanted to compare the recursive and non recursive cases

Et <- eof(arreglo, recursive=FALSE, metho="irlba", nu=neofs)

Er <- eof(arreglo, recursive=TRUE, method = "irlba", nu=neofs)

但是现在我不知道如何绘制一个地图中的第一个EOF，然后是另一个地图中的第二个EOF，该地区是墨西哥，所以我想把海岸线放在里面，这样可以看得更清楚

另外，就我对EOF分析的理解而言，时间应该在列上，而行上的天文台却看着帮助它在我看来是另一种方式，这就是为什么我使用{{ 1}}但你能证实吗？感谢。