terra (version 0.8-6)

plot: Plot a SpatRaster

Description

Plot (that is, make a map of) the values of a SpatRaster or SpatVector

Usage

# S4 method for SpatRaster,numeric
plot(x, y=1, col, 
  type, mar=c(5.1, 4.1, 4.1, 7.1), legend=TRUE, axes=TRUE,
  pal=list(), pax=list(), maxcell=50000, ...)
# S4 method for SpatRaster,missing
plot(x, y, maxcell=50000, nc, nr, main, 
   maxnl=16, ...)

Arguments

x

SpatRaster

y

missing or positive integer indicating the layer(s) to be plotted, or the name of the layer to be mapped

col

character. Colors

type

character. Type of map/legend. One of "continuous", "classes", or "internval"

mar

numeric vector of lenght 4 to set the margins of the plot

legend

logical. Draw a legend?

axes

logical. Draw axes?

pal

list. Parameters for the legend

pax

list. Parameters for drawing axes

maxcell

positive integer. Maximum number of cells to use for the plot

nc

positive integer. Optional. The number of columns to divide the plotting device in (when plotting multiple layers)

nr

positive integer. Optional. The number of rows to divide the plotting device in (when plotting multiple layers)

main

character. Main plot titles (one for each layer to be plotted)

maxnl

positive integer. Maximum number of layers to plot (for a multi-layer object)

...

additional graphical arguments

See Also

image, plotVector, scatterplot

Examples

# NOT RUN {
f <- system.file("ex/test.tif", package="terra") 
r <- rast(f)
plot(r)

plot(r, type="interval")

e <- c(178200,178400,331600,333600)
plot(r, pal=list(ext=e, title="Title\n", title.cex=1.5))

d <- (r > 400) + (r > 600)
plot(d)
plot(d, type="classes")

plot(d, type="interval", levels=0:3) 
plot(d, type="interval", levels=3, pal=list(legend=c("0-1", "1-2", "2-3"))) 
plot(d, type="classes", pal=list(legend=c("M", "X", "A")))

x <- trunc(r/600)
x <- as.factor(x)
levels(x) <- c("earth", "wind", "fire")
plot(x)


# two plots with the same legend
dev.new(width=6, height=4, noRStudioGD = TRUE)
par(mfrow=c(1,2))
plot(r, range=c(100,1800))
plot(r/2, range=c(100,1800))

# as you only need one legend:
par(mfrow=c(1,2))
plot(r, range=c(100,1800), mar=c(4, 3, 4, 3), pal=list(shrink=0.9, cex=.8), 
	pax=list(sides=1:2, cex.axis=.6))
text(182500, 335000, "Two maps, one plot", xpd=NA)
plot(r/2, range=c(100,1800), mar=c(4, 2, 4, 4), legend=FALSE, 
	pax=list(sides=c(1,4), cex.axis=.6))
# }