Elevation raster aggregated using aggregate() function, elevation column added using dplyr package and mutate() function

Different sub-plots based on soil types using facet_wrap(~TYPE_)

Polygons mapped based on their properties

Shapefile's projection was transformed into UTM

Point, line, and polygons

The subtraction of a DTM from a DSM to result in a CHM (canopy height model).

Harvard Forest DTM Data