prioritizr (version 3.0.4)

add_locked_out_constraints: Add locked out constraints

Description

Add constraints to ensure that certain planning units are not prioritized in the solution. For example, it may be useful to lock out planning units that have been degraded and are not longer suitable for conserving species. If specific planning units should be locked in to the solution, use add_locked_in_constraints.

Usage

add_locked_out_constraints(x, locked_out)
# S4 method for ConservationProblem,numeric
add_locked_out_constraints(x, locked_out)
# S4 method for ConservationProblem,character
add_locked_out_constraints(x, locked_out)
# S4 method for ConservationProblem,Spatial
add_locked_out_constraints(x, locked_out)
# S4 method for ConservationProblem,Raster
add_locked_out_constraints(x, locked_out)

Arguments

x

ConservationProblem-class object.

locked_out

Object that determines which planning units that should be locked out. See details for more information.

Value

ConservationProblem-class object.

Details

The locked out planning units can be specified in several different ways:

integer

vector of indices pertaining to which planning units should be locked out.

character

column name in the attribute table values indicating if planning units should be locked out. This option is only available if the planning units in x are a Spatial-class object. The column in the attribute table should have logical (i.e. TRUE or FALSE) values indicating if the planning unit is to be locked out.

Raster-class object

planning units in x that intersect with cells in y are locked out. Specifically, only if the intersect with cells in y are that are not equal to zero or NA.

Spatial-class object.

planning units in x that spatially intersect with locked_in are locked out.

See Also

constraints.

Examples

# NOT RUN {
# create basic problem
p1 <- problem(sim_pu_polygons, sim_features, "cost") %>%
      add_min_set_objective() %>%
      add_relative_targets(0.2)

# create problem with added locked out constraints using integers
p2 <- p1 %>% add_locked_out_constraints(which(sim_pu_polygons$locked_in))

# create problem with added locked out constraints using a field name
p3 <- p1 %>% add_locked_out_constraints("locked_out")

# create problem with added locked out constraints using raster data
p4 <- p1 %>% add_locked_out_constraints(sim_locked_out_raster)

# create problem with added locked out constraints using spatial polygons
# data
locked_out <- sim_pu_polygons[sim_pu_polygons$locked_out == 1, ]
p5 <- p1 %>% add_locked_out_constraints(locked_out)
# }
# NOT RUN {
# solve problems
s1 <- solve(p1)
s2 <- solve(p2)
s3 <- solve(p3)
s4 <- solve(p4)
s5 <- solve(p5)

# plot solutions
par(mfrow = c(3, 2), mar = c(0, 0, 4.1, 0))
plot(s1, main = "none locked out")
plot(s1[s1$solution_1 == 1, ], col = "darkgreen", add = TRUE)

plot(s2, main="locked out (integer input)")
plot(s2[s2$solution_1 == 1, ], col = "darkgreen", add = TRUE)

plot(s3, main="locked out (character input)")
plot(s3[s3$solution_1 == 1, ], col = "darkgreen", add = TRUE)

plot(s4, main = "locked out (raster input)")
plot(s4[s4$solution_1 == 1, ], col = "darkgreen", add = TRUE)

plot(s5, main = "locked out (polygon input)")
plot(s5[s5$solution_1 == 1, ], col = "darkgreen", add = TRUE)
# }
# NOT RUN {
# }