inst/shiny-examples/respineapp/server.R
In ajpelu/respine: Evaluate RESilience of Mediterranean PINE plantations
library('shiny')
library('shinycssloaders')
library('raster')
library('landscapeR')
library('rasterVis')
library('RColorBrewer')
library('sp')
library('rgeos')
### -------------------------------
# Load functions
source('createLandscape.R', local=TRUE)
source('initRichness.R', local = TRUE)
source('dist2nf.R', local = TRUE)
source('disper.R', local=TRUE)
source('disper_time.R', local=TRUE)
### -------------------------------
# Initial configuration
## Create empty landscape
set.seed(123)
ancho <- 63 * 2
alto <- 53 * 2
m <- matrix(nrow=alto, ncol=ancho, byrow = T)
r <- raster(m)
extent(r) <- matrix(c(0, 0, ancho, alto), nrow=2)
r[] <- 0
## Some parameters
line_pol <- 2 ### Line width polygon
pp_value <- 1 ### Value for Pine plantation
nf_value <- 2 ### Value for Natural forest
# Richness range
ri_range <- as.data.frame(
cbind(value = c(0,1,2,3),
lowRich = c(0, 12.82, mean(13.72, 15.62), 1),
upRich = c(0, 13.34, mean(16.11, 19.66), 2)))
# Input year/m2
piBird = (3.7)/50
piMammal = (0.2)/50
# Themes for raster richness
richness_theme <- rasterTheme(region = brewer.pal(9, "YlGn"))
### -------------------------------
# SERVER
shinyServer(
function(input, output, session){
valores <- reactiveValues(
doPlotInitialMap = 0,
doRiqueza = 0,
doTime = 0)
observeEvent(input$doPaisaje, {
# 0 will be coerced to FALSE
# 1+ will be coerced to TRUE
valores$doPlotInitialMap <- input$doPaisaje
})
### ----------------------------------------------
# Density
den_pp <- reactive ({
list(
den = switch(input$density_pp, 'baja' = 100, 'media' = 1250, 'alta' = 3000),
col = switch(input$density_pp,'baja' = '#a1d99b', 'media' = '#238b45','alta' = '#00441b'))
})
# Create landscape
landscapeInit <- reactive({
createLandscape(r, size_pp = input$size_pp, size_nf = input$size_nf, n_nf = input$n_nf)
})
# Past Use
pastUse <- reactive({
switch(input$pp_pastUse, 'Bosque natural' = 'Oak', 'Matorral' = 'Shrubland',
'Pastizal' = 'Pasture','Cultivo' = 'Crop')
})
### ----------------------------------------------
# Dispersion table
## slider conditioned to small_bird slider (see ui.R)
output$mb <- renderUI({
sliderInput(inputId = "mb",
label = "Aves mediano tamaño",
min = 0, max = 100 - input$sb, value = 0)
})
disp <- reactive({
list(persb = input$sb,
permb = input$mb,
perma = (100-(input$sb + input$mb)))
})
output$disptable <- renderTable({
tabla <- cbind(SmallBirds = disp()$persb,
MediumBirds = disp()$permb,
Mammals = disp()$perma)
tabla},
hover = TRUE, spacing = 'xs', align = 'c', digits = 0)
### ----------------------------------------------
## Distance raster
dist_raster <- reactive({
dist2nf(landscapeInit(), nf_value = nf_value)
})
## Compute initial Richnness
rasterRich <- reactive({
initRichness(r = landscapeInit(), draster = dist_raster(),
r_range = ri_range, treedensity = den_pp()$den,
pastUse = pastUse(), rescale = FALSE)
})
## Get bouondary of pp
limit_pp <- reactive({
rasterToPolygons(landscapeInit(), fun=function(x){x==pp_value}, dissolve = TRUE)
})
## extension of Landscape Init
ext <- reactive({
list(
xmin = extent(landscapeInit())@xmin,
xmax = extent(landscapeInit())@xmax,
ymin = extent(landscapeInit())@ymin,
ymax = extent(landscapeInit())@ymax)
})
### ----------------------------------------------
## Compute dispersion rasters
rasterDisp <- reactive({
disper(x = landscapeInit(), xr = rasterRich(), nf_value = nf_value, pp_value = pp_value)
})
## Compute Richness pine plantations
rich_pp <- reactive({
calc(stack(landscapeInit(), rasterRich()), fun=function(x) ifelse(x[1] == pp_value, x[1]*x[2], NA))
})
## Input propagule
propagule <- reactive({
# Compute propagule input by cell
piBird * ((rasterDisp()[['msb']] * disp()$persb) + (rasterDisp()[['mmb']] * disp()$permb)) + (rasterDisp()[['mma']] * disp()$perma) * piMammal
})
## Richness statistics
rich_nf <- reactive({
rich_nf <- calc(stack(landscapeInit(), rasterRich()), fun=function(x) ifelse(x[1] == nf_value, (x[1]/nf_value)*x[2], NA))
})
output$rich_table_init <- renderTable({
tabla <- cbind(
Ecosistema = c("Repoblación de Pinar", "Bosques naturales"),
Media = c(round(cellStats(rich_pp(), mean),2),
round(cellStats(rich_nf(), mean),2)),
Min = c(round(cellStats(rich_pp(), min), 2),
round(cellStats(rich_nf(), min),2)),
Max = c(round(cellStats(rich_pp(), max),2),
round(cellStats(rich_nf(), max),2)))
tabla},
hover = TRUE, spacing = 'l', align = 'c',
digits = 2, striped = TRUE)
### ----------------------------------------------
# Endpoints
## Initial Map
output$initial_map <- renderPlot({
if (valores$doPlotInitialMap == 0) return()
isolate({
colores <- c('lightgoldenrod1', # Crops
'green', # Natural forests
'white', # Other
den_pp()$col) # Pine plantation
key_landuses <- list(text = list(lab = c("Cultivos", "Bosques Naturales","Matorrales", "Pinares")),
rectangles=list(col = colores), space='bottom', columns=4)
levelplot(landscapeInit(), att='landuse', scales=list(draw=FALSE),
col.regions = colores, colorkey=FALSE, key = key_landuses) +
spplot(limit_pp(), fill = "transparent", col = "black",
xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
colorkey = FALSE, lwd=line_pol)
})
})
## Richness Map (initial)
output$richness_map <- renderPlot({
if (valores$doPlotInitialMap == FALSE) return()
isolate({
mapa_riqueza <- rasterRich()
mapa_riqueza[mapa_riqueza == 0] <- NA
levelplot(mapa_riqueza, par.settings = richness_theme, margin = FALSE,
scales=list(draw=FALSE), pretty=TRUE,
colorkey = list(space = "bottom")) +
spplot(limit_pp(), fill = "transparent", col = "black",
xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
colorkey = FALSE, lwd=line_pol)
})
})
## Propagule Input
output$richness_disper <- renderPlot({
levelplot(propagule(),
margin=FALSE, par.settings = RdBuTheme)
})
## Richness End
rich_end <- reactive({
propagulo_time <- rich_pp() + propagule()*input$timeRange
rich_time <- calc(stack(landscapeInit(),
rasterRich(),
propagulo_time),
fun = function(x) ifelse(x[1] == pp_value, x[1]*x[3], x[2]))
rich_time[rich_time== 0] <- NA
list(
rich_pp_end = propagulo_time,
rich_time = rich_time)
})
## Evolution time dispersion
output$richness_disperTime <- renderPlot({
rend <- rich_end()$rich_time
levelplot(stack(rend),
par.settings = richness_theme, margin = FALSE, pretty=TRUE,
scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
spplot(limit_pp(), fill = "transparent", col = "black",
xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
colorkey = FALSE, lwd=line_pol)
})
output$rich_table_end <- renderTable({
tabla <- cbind(
Ecosistema = c("Repoblación de Pinar",
"Repoblación de Pinar final",
"Bosques naturales"),
Media = c(round(cellStats(rich_pp(), mean),2),
round(cellStats(rich_end()$rich_pp_end, mean),2),
round(cellStats(rich_nf(), mean),2)),
Min = c(round(cellStats(rich_pp(), min), 2),
round(cellStats(rich_end()$rich_pp_end, min),2),
round(cellStats(rich_nf(), min),2)),
Max = c(round(cellStats(rich_pp(), max),2),
round(cellStats(rich_end()$rich_pp_end, max),2),
round(cellStats(rich_nf(), max),2)))
tabla},
hover = TRUE, spacing = 'l', align = 'c',
digits = 2, striped = TRUE)
# invalidateLater(millis = 1000, session)
# valores$doTime = isolate(valores$doTime) + 1
#
#
# if(valores$doTime < input$timeRange){
# propagulo_time <- rich_pp() + propagule()*valores$doTime
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(x[1] == pp_value, x[1]*x[3], x[2]))
# rich_time[rich_time== 0] <- NA
#
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
# } else {
# propagulo_time <- rich_pp() + propagule()*input$timeRange
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(x[1] == pp_value, x[1]*x[3], x[2]))
# rich_time[rich_time== 0] <- NA
#
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
# }
# for (i in 1:input$timeRange){
#
# propagulo_time <- rich_pp() + propagule()*i
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(
# x[1] == pp_value, x[1]*x[3], x[2]))
#
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
#
# }
# # valores$doTime == vals$counter
# invalidateLater(millis = 500, session)
# valores$doTime = isolate(valores$doTime) + 1
#
#
#
#
#
#
# if(valores$doTime < input$timeRange) {
# propagulo_time <- rich_pp() + propagule()*valores$doTime
#
# # propagulo_time <- propagule()[['rich_pp']] + (propagule()[['seed_input']])*valores$doTime
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(
# x[1] == pp_value, x[1]*x[3], x[2]))
# names(rich_time) <- paste0('rich_y',valores$doTime)
# rich_time[rich_time == 0] <- NA
# limite <- rasterToPolygons(landscapeInit(), fun=function(x){x==1}, dissolve = TRUE)
# mytheme <- rasterTheme(region = brewer.pal(9, "YlGn"))
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
# } else {
#
# propagulo_time <- rich_pp() + propagule()*valores$doTime*input$timeRange
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(
# x[1] == pp_value, x[1]*x[3], x[2]))
#
# # names(rich_time) <- paste0('rich_y',valores$doTime)
# rich_time[rich_time == 0] <- NA
#
# limite <- rasterToPolygons(landscapeInit(), fun=function(x){x==1}, dissolve = TRUE)
# mytheme <- rasterTheme(region = brewer.pal(9, "YlGn"))
#
# levelplot(stack(rich_time),
# par.settings = mytheme, margin = FALSE,
# scales=list(draw=FALSE),
# colorkey = list(space = "bottom"),
# pretty=TRUE) +
# spplot(limite, fill = "transparent", col = "black",
# xlim = c(extent(landscapeInit())@xmin,
# extent(landscapeInit())@xmax),
# ylim = c(extent(landscapeInit())@ymin,
# extent(landscapeInit())@ymax),
# colorkey = FALSE, lwd=line_pol)
#
#
#
#
# }
}
)
ajpelu/respine documentation built on May 14, 2019, 8:19 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
library('shiny')
library('shinycssloaders')
library('raster')
library('landscapeR')
library('rasterVis')
library('RColorBrewer')
library('sp')
library('rgeos')
### -------------------------------
# Load functions
source('createLandscape.R', local=TRUE)
source('initRichness.R', local = TRUE)
source('dist2nf.R', local = TRUE)
source('disper.R', local=TRUE)
source('disper_time.R', local=TRUE)
### -------------------------------
# Initial configuration
## Create empty landscape
set.seed(123)
ancho <- 63 * 2
alto <- 53 * 2
m <- matrix(nrow=alto, ncol=ancho, byrow = T)
r <- raster(m)
extent(r) <- matrix(c(0, 0, ancho, alto), nrow=2)
r[] <- 0
## Some parameters
line_pol <- 2 ### Line width polygon
pp_value <- 1 ### Value for Pine plantation
nf_value <- 2 ### Value for Natural forest
# Richness range
ri_range <- as.data.frame(
cbind(value = c(0,1,2,3),
lowRich = c(0, 12.82, mean(13.72, 15.62), 1),
upRich = c(0, 13.34, mean(16.11, 19.66), 2)))
# Input year/m2
piBird = (3.7)/50
piMammal = (0.2)/50
# Themes for raster richness
richness_theme <- rasterTheme(region = brewer.pal(9, "YlGn"))
### -------------------------------
# SERVER
shinyServer(
function(input, output, session){
valores <- reactiveValues(
doPlotInitialMap = 0,
doRiqueza = 0,
doTime = 0)
observeEvent(input$doPaisaje, {
# 0 will be coerced to FALSE
# 1+ will be coerced to TRUE
valores$doPlotInitialMap <- input$doPaisaje
})
### ----------------------------------------------
# Density
den_pp <- reactive ({
list(
den = switch(input$density_pp, 'baja' = 100, 'media' = 1250, 'alta' = 3000),
col = switch(input$density_pp,'baja' = '#a1d99b', 'media' = '#238b45','alta' = '#00441b'))
})
# Create landscape
landscapeInit <- reactive({
createLandscape(r, size_pp = input$size_pp, size_nf = input$size_nf, n_nf = input$n_nf)
})
# Past Use
pastUse <- reactive({
switch(input$pp_pastUse, 'Bosque natural' = 'Oak', 'Matorral' = 'Shrubland',
'Pastizal' = 'Pasture','Cultivo' = 'Crop')
})
### ----------------------------------------------
# Dispersion table
## slider conditioned to small_bird slider (see ui.R)
output$mb <- renderUI({
sliderInput(inputId = "mb",
label = "Aves mediano tamaño",
min = 0, max = 100 - input$sb, value = 0)
})
disp <- reactive({
list(persb = input$sb,
permb = input$mb,
perma = (100-(input$sb + input$mb)))
})
output$disptable <- renderTable({
tabla <- cbind(SmallBirds = disp()$persb,
MediumBirds = disp()$permb,
Mammals = disp()$perma)
tabla},
hover = TRUE, spacing = 'xs', align = 'c', digits = 0)
### ----------------------------------------------
## Distance raster
dist_raster <- reactive({
dist2nf(landscapeInit(), nf_value = nf_value)
})
## Compute initial Richnness
rasterRich <- reactive({
initRichness(r = landscapeInit(), draster = dist_raster(),
r_range = ri_range, treedensity = den_pp()$den,
pastUse = pastUse(), rescale = FALSE)
})
## Get bouondary of pp
limit_pp <- reactive({
rasterToPolygons(landscapeInit(), fun=function(x){x==pp_value}, dissolve = TRUE)
})
## extension of Landscape Init
ext <- reactive({
list(
xmin = extent(landscapeInit())@xmin,
xmax = extent(landscapeInit())@xmax,
ymin = extent(landscapeInit())@ymin,
ymax = extent(landscapeInit())@ymax)
})
### ----------------------------------------------
## Compute dispersion rasters
rasterDisp <- reactive({
disper(x = landscapeInit(), xr = rasterRich(), nf_value = nf_value, pp_value = pp_value)
})
## Compute Richness pine plantations
rich_pp <- reactive({
calc(stack(landscapeInit(), rasterRich()), fun=function(x) ifelse(x[1] == pp_value, x[1]*x[2], NA))
})
## Input propagule
propagule <- reactive({
# Compute propagule input by cell
piBird * ((rasterDisp()[['msb']] * disp()$persb) + (rasterDisp()[['mmb']] * disp()$permb)) + (rasterDisp()[['mma']] * disp()$perma) * piMammal
})
## Richness statistics
rich_nf <- reactive({
rich_nf <- calc(stack(landscapeInit(), rasterRich()), fun=function(x) ifelse(x[1] == nf_value, (x[1]/nf_value)*x[2], NA))
})
output$rich_table_init <- renderTable({
tabla <- cbind(
Ecosistema = c("Repoblación de Pinar", "Bosques naturales"),
Media = c(round(cellStats(rich_pp(), mean),2),
round(cellStats(rich_nf(), mean),2)),
Min = c(round(cellStats(rich_pp(), min), 2),
round(cellStats(rich_nf(), min),2)),
Max = c(round(cellStats(rich_pp(), max),2),
round(cellStats(rich_nf(), max),2)))
tabla},
hover = TRUE, spacing = 'l', align = 'c',
digits = 2, striped = TRUE)
### ----------------------------------------------
# Endpoints
## Initial Map
output$initial_map <- renderPlot({
if (valores$doPlotInitialMap == 0) return()
isolate({
colores <- c('lightgoldenrod1', # Crops
'green', # Natural forests
'white', # Other
den_pp()$col) # Pine plantation
key_landuses <- list(text = list(lab = c("Cultivos", "Bosques Naturales","Matorrales", "Pinares")),
rectangles=list(col = colores), space='bottom', columns=4)
levelplot(landscapeInit(), att='landuse', scales=list(draw=FALSE),
col.regions = colores, colorkey=FALSE, key = key_landuses) +
spplot(limit_pp(), fill = "transparent", col = "black",
xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
colorkey = FALSE, lwd=line_pol)
})
})
## Richness Map (initial)
output$richness_map <- renderPlot({
if (valores$doPlotInitialMap == FALSE) return()
isolate({
mapa_riqueza <- rasterRich()
mapa_riqueza[mapa_riqueza == 0] <- NA
levelplot(mapa_riqueza, par.settings = richness_theme, margin = FALSE,
scales=list(draw=FALSE), pretty=TRUE,
colorkey = list(space = "bottom")) +
spplot(limit_pp(), fill = "transparent", col = "black",
xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
colorkey = FALSE, lwd=line_pol)
})
})
## Propagule Input
output$richness_disper <- renderPlot({
levelplot(propagule(),
margin=FALSE, par.settings = RdBuTheme)
})
## Richness End
rich_end <- reactive({
propagulo_time <- rich_pp() + propagule()*input$timeRange
rich_time <- calc(stack(landscapeInit(),
rasterRich(),
propagulo_time),
fun = function(x) ifelse(x[1] == pp_value, x[1]*x[3], x[2]))
rich_time[rich_time== 0] <- NA
list(
rich_pp_end = propagulo_time,
rich_time = rich_time)
})
## Evolution time dispersion
output$richness_disperTime <- renderPlot({
rend <- rich_end()$rich_time
levelplot(stack(rend),
par.settings = richness_theme, margin = FALSE, pretty=TRUE,
scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
spplot(limit_pp(), fill = "transparent", col = "black",
xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
colorkey = FALSE, lwd=line_pol)
})
output$rich_table_end <- renderTable({
tabla <- cbind(
Ecosistema = c("Repoblación de Pinar",
"Repoblación de Pinar final",
"Bosques naturales"),
Media = c(round(cellStats(rich_pp(), mean),2),
round(cellStats(rich_end()$rich_pp_end, mean),2),
round(cellStats(rich_nf(), mean),2)),
Min = c(round(cellStats(rich_pp(), min), 2),
round(cellStats(rich_end()$rich_pp_end, min),2),
round(cellStats(rich_nf(), min),2)),
Max = c(round(cellStats(rich_pp(), max),2),
round(cellStats(rich_end()$rich_pp_end, max),2),
round(cellStats(rich_nf(), max),2)))
tabla},
hover = TRUE, spacing = 'l', align = 'c',
digits = 2, striped = TRUE)
# invalidateLater(millis = 1000, session)
# valores$doTime = isolate(valores$doTime) + 1
#
#
# if(valores$doTime < input$timeRange){
# propagulo_time <- rich_pp() + propagule()*valores$doTime
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(x[1] == pp_value, x[1]*x[3], x[2]))
# rich_time[rich_time== 0] <- NA
#
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
# } else {
# propagulo_time <- rich_pp() + propagule()*input$timeRange
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(x[1] == pp_value, x[1]*x[3], x[2]))
# rich_time[rich_time== 0] <- NA
#
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
# }
# for (i in 1:input$timeRange){
#
# propagulo_time <- rich_pp() + propagule()*i
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(
# x[1] == pp_value, x[1]*x[3], x[2]))
#
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
#
# }
# # valores$doTime == vals$counter
# invalidateLater(millis = 500, session)
# valores$doTime = isolate(valores$doTime) + 1
#
#
#
#
#
#
# if(valores$doTime < input$timeRange) {
# propagulo_time <- rich_pp() + propagule()*valores$doTime
#
# # propagulo_time <- propagule()[['rich_pp']] + (propagule()[['seed_input']])*valores$doTime
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(
# x[1] == pp_value, x[1]*x[3], x[2]))
# names(rich_time) <- paste0('rich_y',valores$doTime)
# rich_time[rich_time == 0] <- NA
# limite <- rasterToPolygons(landscapeInit(), fun=function(x){x==1}, dissolve = TRUE)
# mytheme <- rasterTheme(region = brewer.pal(9, "YlGn"))
# levelplot(stack(rich_time),
# par.settings = richness_theme, margin = FALSE, pretty=TRUE,
# scales=list(draw=FALSE), colorkey = list(space = "bottom")) +
# spplot(limit_pp(), fill = "transparent", col = "black",
# xlim = c(ext()$xmin, ext()$xmax), ylim = c(ext()$ymin, ext()$ymax),
# colorkey = FALSE, lwd=line_pol)
# } else {
#
# propagulo_time <- rich_pp() + propagule()*valores$doTime*input$timeRange
#
# rich_time <- calc(stack(landscapeInit(), rasterRich(), propagulo_time),
# fun = function(x) ifelse(
# x[1] == pp_value, x[1]*x[3], x[2]))
#
# # names(rich_time) <- paste0('rich_y',valores$doTime)
# rich_time[rich_time == 0] <- NA
#
# limite <- rasterToPolygons(landscapeInit(), fun=function(x){x==1}, dissolve = TRUE)
# mytheme <- rasterTheme(region = brewer.pal(9, "YlGn"))
#
# levelplot(stack(rich_time),
# par.settings = mytheme, margin = FALSE,
# scales=list(draw=FALSE),
# colorkey = list(space = "bottom"),
# pretty=TRUE) +
# spplot(limite, fill = "transparent", col = "black",
# xlim = c(extent(landscapeInit())@xmin,
# extent(landscapeInit())@xmax),
# ylim = c(extent(landscapeInit())@ymin,
# extent(landscapeInit())@ymax),
# colorkey = FALSE, lwd=line_pol)
#
#
#
#
# }
}
)
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.