Índice de Vegetação de Diferença Normalizada (NDVI) de uma célula de grade no Alasca central. NDVI é uma medida da verdura da vegetação e está relacionada à cobertura de vegetação verde, atividade fotossintética e biomassa verde. NDVI varia entre -1 e 1.

O solo descoberto e a neve geralmente tem valores de NDVI abaixo de 0.2, enquanto que as áreas vegetadas têm valores de NDVI acima de 0.2. Esta série temporal do NDVI é do conjunto de dados GIMMS (Global Inventory, Monitoring, and Modeling Studies), em Tucker et al. 2005, que fornece estimativas do NDVI a partir de observações de satélite AVHRR (Advanced Very High Resolution Radiometer).

Referência
Tucker, C.; Pinzon, J.; Brown, M.; Slayback, D.; Pak, E.; Mahoney, R.; Vermote, E.; El Saleous, N.,
An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT
VEGETATION NDVI data. International Journal of Remote Sensing 2005, 26, 4485-4498. 

ndvi.obs = ts(read.table("http://leg.ufpr.br/~lucambio/CE017/20212S/ndvi.txt", col.names = "ndvi"), 
              start = c(1982,1), freq = 12)
par(mar=c(4,4,1,1))
plot(ndvi.obs, pch=16, xlab="TEMPO", ylab="NDVI", type="b", col="black", main="")
abline(h=mean(ndvi.obs, na.rm = TRUE), col="red", lwd=2)
grid()

library(astsa)
monthplot(ndvi.obs, xlab="MESES", ylab="NDVI", type="l")
grid()

Observamos que, para cada mês, nos últimos anos sempre os valores de NDVI são menores do que nos primeiros anos observados. 

par(mar=c(4,4,1,1))
plot(ndvi.obs, pch=16, xlab="TEMPO", ylab="NDVI", type="b", col="black", main="")
lines(ksmooth(time(ndvi.obs), ndvi.obs, "normal", bandwidth=1), lwd=2, col=4)
grid()

Eliminando a tendência

ndvi.obs1 = ndvi.obs - ksmooth(time(ndvi.obs), ndvi.obs, "normal", bandwidth=1)$y
par(mar=c(4,4,1,1))
plot(ndvi.obs1, pch=16, xlab="TEMPO", ylab="NDVI", type="b", col="black", main="")
abline(h=mean(ndvi.obs1), col="red", lwd=2)
grid()

Observando a estrutura de autocorrelção

lag1.plot(ndvi.obs1, 12)

acf2(ndvi.obs1)

##      [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
## ACF  0.64  0.35 -0.01 -0.37 -0.67 -0.74 -0.65 -0.34 -0.01  0.36  0.62  0.69
## PACF 0.64 -0.11 -0.33 -0.38 -0.44 -0.29 -0.28 -0.07 -0.12 -0.03  0.03 -0.01
##      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
## ACF   0.61  0.33  0.00 -0.37 -0.63 -0.72 -0.60 -0.33  0.05  0.35  0.63  0.72
## PACF  0.09 -0.05  0.01 -0.13 -0.19 -0.17 -0.15 -0.11 -0.03 -0.17 -0.04 -0.01
##      [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36]
## ACF   0.60  0.32 -0.02 -0.35 -0.59 -0.69 -0.58 -0.33  0.01  0.35  0.57  0.65
## PACF -0.02 -0.04 -0.04  0.02 -0.01 -0.03 -0.01 -0.06 -0.02 -0.01 -0.03 -0.03
##      [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48]
## ACF   0.57  0.30  0.00 -0.36 -0.58 -0.65 -0.55 -0.31  0.03  0.33  0.57  0.65
## PACF  0.03 -0.05  0.09 -0.08  0.00 -0.01 -0.01  0.00  0.00 -0.03  0.01  0.04
monthplot(ndvi.obs1, xlab="MESES", ylab="NDVI", type="l")