Add posterior probability surface to dynamic map

Source: R/plot.R

Add posterior probability surface to dynamic map

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

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.
realised

if TRUE then plot surface for realised sources only.
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

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