if (FALSE) {
require(Rfssa)
load_github_data("https://github.com/haghbinh/Rfssa/blob/master/data/Jambi.RData")
## Raw image data
NDVI <- Jambi$NDVI
EVI <- Jambi$EVI
## Kernel density estimation of pixel intensity
D0_NDVI <- matrix(NA, nrow = 512, ncol = 448)
D0_EVI <- matrix(NA, nrow = 512, ncol = 448)
for (i in 1:448) {
D0_NDVI[, i] <- density(NDVI[, , i], from = 0, to = 1)$y
D0_EVI[, i] <- density(EVI[, , i], from = 0, to = 1)$y
}
d <- 11
u <- seq(0, 1, length.out = 512)
Y_1 <- Rfssa::fts(list(D0_NDVI), list(list(d, "bspline")), list(u))
Y_2 <- Rfssa::fts(list(D0_EVI), list(list(d, "bspline")), list(u))
out <- cor.fts(Y_1, Y_2)
print(out)
}
Run the code above in your browser using DataLab