# \donttest{
library(sf)
library(terra)
e <- ext(179407.8, 181087.9, 331134.4, 332332.1)
r1 <- rast(e, resolution=20)
r1[] <- sample(1:5, ncell(r1), replace=TRUE)
r2 <- rast(e, resolution=20)
r2[] <- sample(1:5, ncell(r2), replace=TRUE)
d = 5 # kernel
( r.kappa <- raster.change(r1, r2, s = d) )
( r.ttest <- raster.change(r1, r2, s = d, stat="t.test") )
( r.ent <- raster.change(r1, r2, s = d, stat="entropy") )
( r.cor <- raster.change(r1, r2, s = d, stat="cor") )
( r.ce <- raster.change(r1, r2, s = d, stat = "cross-entropy") )
( r.kl <- raster.change(r1, r2, s = d, stat = "divergence") )
opar <- par(no.readonly=TRUE)
par(mfrow=c(3,2))
plot(r.kappa, main="Kappa")
plot(r.ttest[[1]], main="Paired t-test")
plot(r.ent, main="Delta Entropy")
plot(r.cor, main="Rank Correlation")
plot(r.kl, main="Kullback-Leibler")
plot(r.ce, main="cross-entropy")
par(opar)
# }
Run the code above in your browser using DataLab