README.md

In fgashakamba/hypsoLoop: Drawing hypsographic curves and calculate hypsometric integrals of the sub-catchments of a watershed

hypsoLoop

hypsoLoop has been developed as tool to automate the drawing of hypsographic curves and the calculation of hypsometric curves for a watershed that contains several sub-catchments.

Installation

You can install the released version of hypsoLoop from CRAN with:

install.packages("hypsoLoop")

And the development version from GitHub with:

# install.packages("devtools")
devtools::install_github("fgashakamba/hypsoLoop")

Example use of hypsoLoop

The workhorse function of hypsoLoop is called drawHypsocurves(). It accepts two arguments: a watersheds boundary layer (an sf or sp object) and an digital elevation model (DEM) covering the area of interest (a rasterLayer object). Additionally, the user may specify whether they want to get an output folder created in the root directory or not. If not, the hypsometric curves are drawn in the viewer and a summary table containing the hypsometric integers and other useful information on the sub-catchments is returned.

library(hypsoLoop)
drawHypsoCurves(watersheds, DEM)

#>   SUB_CATCHMENT_CODE MIN_ELEV MAX_ELEV TOTAL_AREA H_INTEGRAL
#> 1            Mulindi     1616     2225    2571.53      0.736
#> 2          Cyonyonyo     1614     2217    3666.87      0.791
#> 3              Yanze     1368     2071    3431.98      0.857

Other two functions are exposed from this hypsoLoop. These are the generateHypsoTables() which can be used to generate area-elevation tables for the sub-catchments of a watersheds and whose input requirements are the same as drawHypsoCurves(). The calc_areas() function is a handy function that be used to calculate the areas covered by each class of a categorical raster. It accepts one input argument which has to be a categorical raster of type rasterLayer and returns a table that summarizes the areas (in Hectares) covered by each class or category defined in the input raster.

library(hypsoLoop)
calc_areas(lulc_yanze)
#> # A tibble: 7 x 2
#>   CLASS   AREA
#>   <int>  <dbl>
#> 1     1 1379. 
#> 2     2  531. 
#> 3     3 6857. 
#> 4     5   94.3
#> 5     6    0.2
#> 6    10  424. 
#> 7    11  403.

fgashakamba/hypsoLoop documentation built on Nov. 13, 2022, 8:17 a.m.