## ---- vect_1------------------------------------------------------------------
library(terra)
p <- vect(system.file("ex/lux.shp", package="terra"))


## ---- maps_1, fig.height=5----------------------------------------------------
plot(p)


## ---- maps_2, fig.height=5----------------------------------------------------
n <- nrow(p)
plot(p, col=rainbow(n))


## ---- maps_3------------------------------------------------------------------
plot(p, "NAME_2", col=rainbow(25))


## ---- maps_4a, fig.width=10---------------------------------------------------
plot(p, c("NAME_1", "NAME_2"), col=rainbow(25))


## ---- maps_4b, fig.width=10---------------------------------------------------
par(mfrow=c(1,2))
m <- c(3.1, 3.1, 2.1, 2.1)
plot(p, "NAME_1", col=rainbow(25), mar=m, plg=list(x="topright"), pax=list(las=1))
plot(p, "NAME_2", col=rainbow(25), mar=m, plg=list(x="topright", cex=.75), pax=list(las=1))


## ---- maps_5, fig.width=10----------------------------------------------------
par(mfrow=c(1,2))
m <- c(3.1, 3.1, 1.1, 1.1)
plot(p, "NAME_1", col=rainbow(25), mar=m, plg=list(x="topright", title="District", bty = "o"), main="", axes=FALSE)
axis(1, at=c(5,7)); axis(1)
axis(2, at=c(49,51)); axis(2, las=1)

plot(p, "NAME_2", col=rainbow(25), mar=m, plg=list(x="topright", cex=.75, title="Canton", bty = "o"), main="", axes=FALSE)
axis(1, at=c(5, 7)); axis(1)


## ---- maps_6------------------------------------------------------------------
d <- aggregate(p, "NAME_1")
plot(p, col="light blue", lty=2, border="red", lwd=2)
lines(d, lwd=5)
lines(d, col="white", lwd=1)
text(p, "NAME_2", cex=.8, halo=TRUE)


## ---- maps_10-----------------------------------------------------------------
f <- system.file("ex/elev.tif", package="terra")
r <- rast(f)


## ---- maps_11-----------------------------------------------------------------
plot(r)


## ---- maps_12-----------------------------------------------------------------
plot(r)

lines(p, lwd=2)
set.seed(12)
xy <- spatSample(r, 20, "random", na.rm=TRUE, xy=TRUE)
points(xy, pch=20, col="red", cex=2)


## ---- maps_13-----------------------------------------------------------------
m <- c(3.1, 3.1, 1.1, 1.1)
plot(r, type="interval", plg=list(x="topright"), mar=m)


## ---- maps_14-----------------------------------------------------------------
rr <- round(r/100)
plot(rr, plg=list(x="topright"), mar=m)


## ---- maps_15, fig.height=5, fig.width=9--------------------------------------
set.seed(0)
x <- rast(nrow=10, ncol=10)
values(x) <- sample(0:2, ncell(r), replace=TRUE)
cls <- c("forest", "water", "urban")
levels(x) <- cls
names(x) <- "land cover"
is.factor(x)
x
plot(x, col=c("green", "blue", "light gray"))
text(x, digits=3, cex=.75, halo=TRUE)


## ---- maps_16-----------------------------------------------------------------
library(terra)
b <- rast(system.file("ex/logo.tif", package="terra"))
plot(b)


## ---- maps_17-----------------------------------------------------------------
r <- rast(p, res=0.01 )
values(r) <- 1:ncell(r)
r <- mask(r, p)


## ---- maps_20-----------------------------------------------------------------
plotRGB(b, r=1, g=2, b=3)


## ---- maps_data---------------------------------------------------------------
library(terra)
f <- system.file("ex/lux.shp", package="terra")
p <- vect(f)


## ---- maps_21-----------------------------------------------------------------
library(maptiles)
bg <- get_tiles(ext(p))
plotRGB(bg)
lines(p, col="blue", lwd=3)


## ----leaflet, message=FALSE---------------------------------------------------
library(sf)
library(leaflet)
leaflet(sf::st_as_sf(p)) %>%
   addTiles() %>%
   addPolygons(fillOpacity=0.1)