rasterComp: Comparative raster mapping

In voluModel: Modeling Species Distributions in Three Dimensions

Comparative raster mapping

Description

A convenient wrapper around terra::plot to generate formatted plots comparing two rasters. This is used in the context of voluModel to overlay semi-transparent distributions (coded as 1) in two different RasterLayers.

Usage

rasterComp(
  rast1 = NULL,
  rast2 = NULL,
  col1 = "#1b9e777F",
  col2 = "#7570b37F",
  rast1Name = "Set 1",
  rast2Name = "Set 2",
  land = NA,
  landCol = "black",
  title = "A Raster Comparison",
  graticule = TRUE,
  ...
)

Arguments

rast1	A single SpatRaster showing the distribution of the species corresponding to rast1Name. Should have values of 0 (absence) and 1 (presence). Can also be NULL.
rast2	A single SpatRaster showing the distribution of the species corresponding to rast2Name. Should have values of 0 (absence) and 1 (presence). Must match the extent and resolution of rast1. Can also be NULL.
col1	Color for rast1 presences
col2	Color for rast2 presences
rast1Name	An optional name for the first set of occurrences, which will be color-coded to occs1Col in the resulting plot.
rast2Name	An optional name for the first set of occurrences, which will be color-coded to occs2Col in the resulting plot.
land	An optional coastline polygon shapefile of types sf or SpatRaster to provide geographic context for the occurrence points.
landCol	Color for land on map.
title	A title for the plot.
graticule	Do you want a grid of lon/lat lines?
...	Additional optional arguments to pass to terra::plot().

Value

A plot of class recordedplot overlaying mapped, semitransparent extents of the input rasters

Note

The extents of rast1 and rast2 must match.

See Also

plot

Examples

library(terra)
rast1 <- rast(ncol=10, nrow=10)
values(rast1) <- rep(0:1, 50)

rast2 <- rast(ncol=10, nrow=10)
values(rast2) <- c(rep(0, 50), rep(1,50))

rasterComp(rast1 = rast1, rast2 = rast2)

voluModel documentation built on Aug. 8, 2025, 6:26 p.m.