Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.width = 6,
fig.height = 4)
## ----setup, include = FALSE---------------------------------------------------
library(topoDistance)
library(gdistance)
library(RColorBrewer)
library(scales)
## -----------------------------------------------------------------------------
plot(Yosemite$DEM, col = terrain.colors(99))
## -----------------------------------------------------------------------------
xy <- matrix(ncol = 2, byrow = TRUE,
c(-119.5566, 37.72474,
-119.5157, 37.76688,
-119.4718, 37.76078))
colnames(xy) <- c("longitude", "latitude")
xy
## -----------------------------------------------------------------------------
tdist <- topoDist(Yosemite$DEM, xy, paths = TRUE)
tdist
## -----------------------------------------------------------------------------
td.mat <- tdist[[1]]
td.mat[td.mat > 8000] <- NA
td.mat
## -----------------------------------------------------------------------------
topoPathMap(Yosemite$DEM, xy, topoPaths = tdist, type = "hillshade",
pathWidth = 4, cex = 2, bg = "blue")
## -----------------------------------------------------------------------------
sdmColors <- colorRampPalette(c("blue", "green", "yellow", "orange", "red"), space = "rgb", interpolate = "linear")
plot(Yosemite$SDM, col = sdmColors(99))
## -----------------------------------------------------------------------------
tLCP <- topoLCP(Yosemite$DEM, costSurface = Yosemite$SDM, pts = xy, paths = TRUE)
tLCP
## -----------------------------------------------------------------------------
topoPathMap(Yosemite$DEM, xy, topoPaths = tLCP, type = "hillshade",
costSurface = Yosemite$SDM, pathWidth = 4, pathColor = "purple")
## -----------------------------------------------------------------------------
topoPathMap(Yosemite$DEM, xy, topoPaths = tLCP, type = "hillshade",
costSurface = Yosemite$SDM, pathWidth = 4, pathColor = "purple")
lines(tdist[[2]], lty = 2, lwd = 2)
## ---- warning=FALSE-----------------------------------------------------------
topoProfile(Yosemite$DEM, topoPaths = tLCP, pts = 1000,
type = "base", singlePlot = TRUE)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.