MOP Y EXDET

setwd("G:/CursoNichos2019/ExDet ejemplo/ejemplo")

#todos los ascci_bios deberan estar en la misma carpeta, en este caso
library(ntbox)

m_stack <- raster::stack("AusBio13.asc", "AusBio14.asc", "AusBio5.asc", "AusBio6.asc")

raster::plot(m_stack) #para verificar nuestras bios en "M"

g_stack <- raster::stack(c("SaBio13.asc", "SaBio14.asc", "SaBio5.asc", "SaBio6.asc"))

raster::plot(g_stack, ylim=c(-35,-20), xlim=c(15,35)) #para verificar nuestras bios en "proyected"

NT2 <- exdet_multvar(M_stack = m_stack,
                     G_stack = g_stack)
raster::plot(NT2, col=rainbow(100), ylim=c(-35,-20), xlim=c(15,35))

jpeg(filename="Mapa de ExDet (NT2).jpeg",   # Nombre del archivo y extension
     width = 16,    # Anchura
     height = 12,   # Altura
     res= 500,       # Resolucion 500 ppi es un estandar
     units = "cm")  # Unidades

raster::plot(NT2, col=rainbow(100), ylim=c(-35,-20), xlim=c(15,35))
dev.off() #Cierra la Figura, si no no abre en word

## png 
##   2

# Writing the raster from memory to disk (as a new file)
raster::writeRaster(NT2, filename="aqui esta el NT2 en ascii.asc", format="ascii",overwrite=TRUE)

#AQUI HACE EL MOP
mop_res <- mop(M_stack = m_stack,
               G_stack = g_stack, percent = 10,
               comp_each=2000)

## Computation progress:  25 % 
## Computation progress:  50 % 
## Computation progress:  75 %

raster::plot(mop_res, ylim=c(-35,-20), xlim=c(15,35))

   #exportar nuestro MOP          
jpeg(filename="Mapa de MOP.jpeg",   # Nombre del archivo y extension
     width = 16,    # Anchura
     height = 12,   # Altura
     res= 500,       # Resolucion 500 ppi es un estandar
     units = "cm")  # Unidades

raster::plot(mop_res, col=rainbow(100), ylim=c(-35,-20), xlim=c(15,35))
dev.off() #Cierra la Figura, si no no abre en wor

## png 
##   2

# Writing the raster from memory to disk (as a new file)
raster::writeRaster(mop_res, filename="aqui esta el MOP en ascii.asc", format="ascii",overwrite=TRUE)