Add risk surface to dynamic map

Source: R/plot.R

Add posterior probability surface to dynamic map

overlay_risk_map(
  myplot,
  project,
  K = NULL,
  source = NULL,
  threshold = 0.1,
  col = rev(col_hotcold()),
  opacity = 0.8,
  smoothing = 1,
  legend = FALSE,
  iterations = 50
)

Arguments

myplot

dynamic map produced by plot_map() function.
project

an RgeoProfile project, as produced by the function rgeoprofile_project().
K

which value of K to plot.
source

which source to plot. If NULL then plot combined surface.
threshold

what proportion of posterior probability surface to plot.
col

set of plotting colours.
opacity

opacity of posterior probability surface (that is not invisible due to being below threshold).
smoothing

what level of smoothing to apply to posterior probability surface. Smoothing is applied using the raster function disaggregate, with method = "bilinear".
legend

whether or not a legend is plotted
iterations

the number of random parameter sets to generate and combine risk maps

Examples

p <- rgeoprofile_file("tutorial1_project.rds")
plot1 <- plot_map()
plot1 <-overlay_surface(myplot = plot1, project = p, K = 2, source = NULL, threshold = 0.5)
plot1