-
Notifications
You must be signed in to change notification settings - Fork 0
/
test.R
35 lines (27 loc) · 1.08 KB
/
test.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
setwd("C:/Users/Marlon/Documents/R/Heterogeneity/")
peru <- raster::stack(list.files("Peru/clima/", pattern = ".asc$", full.names = T))
raster::plot(peru[[1]])
ext <- raster::extent(-81, -78, -7, 0)
#var.stack <- raster::crop(peru, ext)
var.stack <- peru
cuba <- raster::stack(list.files("Bios/", pattern = ".asc$", full.names = T))
var.stack <- cuba
dist.window <- raster::xres(var.stack) * 2
#ncell.window <- 2
nlayers.mean <- dim(var.stack)[3]
var.scale <- TRUE
rescale.result <- TRUE
multi.parallel <- TRUE
comp.each <- 7000
n.cores <- 6
ex.fun <- 2
# simple parallel
system.time({hetero <- kuenm_heterogen(var.stack = var.stack, dist.window = dist.window)})
# multi parallel
system.time({hetero <- kuenm_heterogen(var.stack = var.stack, dist.window = dist.window, nlayers.mean = nlayers.mean,
var.scale = var.scale, multi.parallel = multi.parallel, comp.each = comp.each)})
x11()
raster::plot(hetero)
x11()
raster::plot(raster::raster("cuba_heterogeneity.tif"))
raster::writeRaster(hetero, filename = "cuba_heterogeneity.tif", format = "GTiff")