library(raster)
library(arc2r)
library(dplyr)
data("swissAlti3D")

plot(swissAlti3D)
histinfo <- hist(swissAlti3D)
# create a reclassification matrix
reclass_vec <- histinfo$breaks

reclass_m <- matrix(c(
  lag(reclass_vec,default = 0),
  reclass_vec,
  seq_along(reclass_vec)
),ncol = 3)


reclass_m
swissAlti3D_reclass <- reclassify(swissAlti3D, reclass_m)

plot(swissAlti3D_reclass)


arc2r/book documentation built on March 5, 2021, 2:10 p.m.