Построение основных компонентов растрового стека в R - PullRequest
Новая
       31

Построение основных компонентов растрового стека в R

0 голосов
/ 28 мая 2018

Я пытаюсь воспроизвести рисунок 1 из здесь , который помечен как «Стандартизированные повреждения основных компонентов с 1 по 8, полученные из ежемесячного NDV; данные, январь 1986 г. по декабрь 1988 г.».

Я использую principal из psych пакета.Вот пример, где у меня есть растровый стек, значения и каждая координата.Предположим, что каждый растр представляет переменную в течение года. 

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, как я могу получить такой же график, как бумага с корреляциями?

...